Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

velocity becomes nan after few steps when initial velocity set to 0.2 #8

Open
Comma0103 opened this issue Dec 21, 2023 · 1 comment

Comments

@Comma0103
Copy link

I set the initial inlet velocity of Karman Vortex Street to 0.2, and kept other configs as below:

lbm = lbm_solver(
            "Karman Vortex Street",
            801,
            201,
            0.01,
            [0, 0, 1, 0],
            [[0.2, 0.0], [0.0, 0.0], [0.0, 0.0], [0.0, 0.0]],
            1,
            [160.0, 100.0, 20.0])

Then I ran it by python lbm_solver.py, but got velocity and vorticity nan after first few steps:
image
Is this an internal problem of the solver? How can I make it more robust? Thanks.

@hietwll
Copy link
Owner

hietwll commented Dec 27, 2023

increasing the inlet velocity leads to a larger Re number. to perform a stable simulation, you need to increase the grid resolution and decrease the time step.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants