Skip to content

Example for ibm + cylinder with the current periodic poisson solver#260

Merged
mathrack merged 3 commits intomainfrom
257-example-for-ibm-cylinder-with-the-current-periodic-poisson-solver
Feb 16, 2026
Merged

Example for ibm + cylinder with the current periodic poisson solver#260
mathrack merged 3 commits intomainfrom
257-example-for-ibm-cylinder-with-the-current-periodic-poisson-solver

Conversation

@mathrack
Copy link
Contributor

No description provided.

@mathrack mathrack self-assigned this Dec 11, 2025
@mathrack mathrack linked an issue Dec 11, 2025 that may be closed by this pull request
@mathrack mathrack marked this pull request as ready for review December 11, 2025 13:45
@ia267 ia267 added this to the Immersed Boundary Method milestone Dec 11, 2025
@ia267 ia267 added the examples Example simulation case label Dec 11, 2025
@ia267 ia267 self-requested a review December 11, 2025 14:28
Copy link
Collaborator

@ia267 ia267 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Approved - I was running this case and it diverged around iteration 3170. I think for now we could set the n_iters = 3000

@mathrack
Copy link
Contributor Author

mathrack commented Feb 16, 2026

Approved - I was running this case and it diverged around iteration 3170. I think for now we could set the n_iters = 3000

I have reproduced the issue on CPU. The number of time steps is now limited to 3000 to avoid divergence. If one wants to perform a periodic simulation with a longer simulation time, I think there are two solutions :

  1. Extend the domain. A domain twice longer would allow a simulation twice longer.
  2. Add a buffer layer to damp fluctuations at the outflow
  3. Reduce the time step. The IBM might absorb the incoming velocity fluctuations with a smaller time step

Of course, the best solution would be to use inflow-outflow instead of periodicity. Hopefully, this will be ready shortly.

@mathrack mathrack requested a review from ia267 February 16, 2026 14:57
Copy link
Collaborator

@ia267 ia267 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Make sense to limit it to 3000

@mathrack mathrack merged commit 930cf52 into main Feb 16, 2026
2 checks passed
@mathrack mathrack deleted the 257-example-for-ibm-cylinder-with-the-current-periodic-poisson-solver branch February 16, 2026 15:41
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

examples Example simulation case

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Example for IBM + cylinder with the current periodic poisson solver

2 participants