Hi Fabian,


I don't know what is causing these oscillations (maybe you need a lower CFL), but regarding your other question: EulerADCFE allows using artificial diffusion for shock capture, while EulerCFE is just the standard Euler solver. For EulerADCFE you need to set the "ShockCaptureType" property (possible values are "Smooth" and "NonSmooth") and the parameters described in section 9.4.1 of the user guide.


In the most recent version of the master branch (this was changed just a few days ago), it is also possible to use shock capture directly with EulerCFE. In this newer version, EulerCFE and EulerADCFE are the same thing, with the second name kept for backwards compatibility.


Cheers,

Douglas 


From: nektar-users-bounces@imperial.ac.uk <nektar-users-bounces@imperial.ac.uk> on behalf of Selbach, Fabian <fabian.selbach@student.uni-siegen.de>
Sent: 30 November 2016 13:21:42
To: nektar-users
Subject: [Nektar-users] Nonphysical terminated simulation - CompressibleFlowSolver
 
Dear users,

I hope you can help one more time. I got nonphysical oscillations during my airfoil-simulation. At that time my simulations terminate with the message: "NaN found during time integration". For more information of the error and parameters see pictures below. I tried to solve this by a higher polynomial degree up to p=7 (this generates 1.5 Mio grid-points here for a 2D simulation), but the simulation terminates still. The time step itself is computed by a fixed CFL=0.9 condition. Maybe this is the wrong way, because this condition is not sufficient?

So please can somebody give me a hint about that? Do I have to use another kind of grid, because the resolution of the boundary at the peak of the airfoil is not high enough or the grid lines have no be orthogonal to the surface? I am already using the spherigons-tool and tried to get a smooth boundary of my airfoil.

I am using the CompressibleFlowSolver on a parallel run up to 96 CPUs by solving the Euler-Equations (EulerCFE) and the Navier-Stokes-Equations (NavierStokesCFE). The simulation terminates for both. 

By the way what is the difference between EulerCFE and EulerADCFE (there are no information about it anywhere)? I tried to run the example NACA0012 of the paper "Nektar++: An open-source spectral/hp element framework". The conditions use the EulerADCFE, but by running this example the simulation does not terminate, but also does not change (values of every variable in the whole domain are constant over the time - no changes?).

Thanks for any help!

Best regards

Fabian


Visualization of the grid at a polynomial degree of p=2 just before terminating:


With auto-scale:


Same solution with scaling by user:


Please take into account the scaling above.


Error and conditions for example p=2:




Visualization of the grid at a polynomial degree of p=7 just before terminating:


(a few time steps after that the simulation terminates and the values go against inf...)