|
This
email from oct_amg1289@hotmail.com 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 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.