search

HiFiLES / HiFiLES-solver

High Fidelity Large Eddy Simulation Solver
Other
172 stars 131 forks source link

About 2D Taylor-Green Vortex #59

Closed arctanx closed 9 years ago

arctanx commented 9 years ago

Hi all, I tried to run a 2D Taylor-Green Vortex case with the initial conditions originally given at line 555 (Simple 2D div-free vortex) in eles.cpp. The input file was the one provided for the 3D TGV test case. The 2D triangular mesh was generated using Gmsh. The problem is that I kept running into the "Fatal error 'NaN in residual, exiting.' at ../src/eles.cpp:1681" after writing Tecplot file number 0. Could anyone shed a light on what I should do to get it running? Thanks in advance! Best Regards arctanx

mlopez14 commented 9 years ago

Perhaps the time-step size is too large. What time-step and polynomial order are you using?

arctanx commented 9 years ago

I've tried time-step size as small as 1E-9, and the polynomial order down to 1, still not working.

mlopez14 commented 9 years ago

You've found a bug! Pressure is being initialized to zero. See line 556 in eles.cpp. I'll make the quick change and let you know when it's ready.

mlopez14 commented 9 years ago

Fixed! Let us know if it works now. Remember to pull the latest change.

arctanx commented 9 years ago

It is working now! Thanks a lot. I should have noticed the bug.

mlopez14 commented 9 years ago

Thanks for reporting it! On Mon, May 18, 2015 at 10:45 PM arctanx notifications@github.com wrote:

It is working now! Thanks a lot. I should have noticed the bug.

— Reply to this email directly or view it on GitHub https://github.com/HiFiLES/HiFiLES-solver/issues/59#issuecomment-103346471 .