******************* This email originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list https://spam.ic.ac.uk/SpamConsole/Senders.aspx to disable email stamping for this address. ******************* Hi, I've been trying to work with this solver in a high sub-sonic compressible regime (Ma=0.87), however this appears: ======================================================================= EquationType: NavierStokesCFE Session Name: Cilindro Spatial Dim.: 2 Max SEM Exp. Order: 3 Num. Processes: 3 Expansion Dim.: 2 Riemann Solver: HLLC Advection Type: Projection Type: Discontinuous Galerkin Diffusion Type: Advection: explicit AdvectionType: WeakDG Diffusion: explicit Time Step: 1e-05 No. of Steps: 60000 Checkpoints (steps): 400 Integration Type: ClassicalRungeKutta4 ======================================================================= Initial Conditions: - Field rho: rhoInf - Field rhou: rhoInf * uInf - Field rhov: rhoInf * vInf - Field E: pInf / (Gamma - 1) + 0.5 * rhoInf * (uInf * uInf + vInf * vInf) Writing: "Cilindro_0.chk" (0.00527309s, XML) Fatal : Level 0 assertion violation NaN found during time integration. Fatal : Level 0 assertion violation NaN found during time integration. Fatal : Level 0 assertion violation NaN found during time integration. -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun detected that one or more processes exited with non-zero status, thus causing the job to be terminated. The first process to do so was: Process name: [[28957,1],1] Exit code: 1 -------------------------------------------------------------------------- ************************************************************************************************************************************************************************ This is the input: ************************************************************************************************************************************************************************ <EXPANSIONS> <E COMPOSITE="C[0]" NUMMODES="3" FIELDS="rho,rhou,rhov,E" TYPE="MODIFIED" /> </EXPANSIONS> <CONDITIONS> <!-- M = 0.2, Re = 200, Pr = 0.72 --> <PARAMETERS> <P> TimeStep = 0.00001 </P> <P> FinTime = 0.6 </P> <P> NumSteps = FinTime/TimeStep </P> <P> IO_CheckSteps = 400 </P> <P> IO_InfoSteps = 400 </P> <P> GasConstant = 287.058 </P> <P> Gamma = 1.4 </P> <P> pInf = 22600 </P> <P> rhoInf = 0.364 </P> <P> Mach = 0.87 </P> <P> cInf = sqrt(Gamma * pInf / rhoInf) </P> <P> uInf = Mach*cInf </P> <P> vInf = 0.0 </P> <P> Twall = 258.5 </P> <P> Re = 400000 </P> <P> mu = 0.00001167 </P> <P> Pr = 0.72 </P> </PARAMETERS> <SOLVERINFO> <I PROPERTY="EQType" VALUE="NavierStokesCFE" /> <I PROPERTY="Projection" VALUE="DisContinuous" /> <I PROPERTY="AdvectionType" VALUE="WeakDG" /> <I PROPERTY="DiffusionType" VALUE="LDGNS" /> <I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4"/> <I PROPERTY="UpwindType" VALUE="HLLC" /> <I PROPERTY="ProblemType" VALUE="General" /> <I PROPERTY="ViscosityType" VALUE="Constant" /> </SOLVERINFO> ************************************************************************************************************************************************************************** Is the problem related to the quality of the mesh or it is related to the Solver info module? Any help is extremely appreciated. Thanks for your time. Morales González Octavio Augusto. U-MOFABI Facultad de Ingeniería, Universidad Nacional Autónoma de México.