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.

_______________________________________________
Nektar-users mailing list
Nektar-users@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/nektar-users