******************* 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.
Hi Octavio, It could be that you need to decrease the time step size. Sometimes the required time step is very small when using explicit time integration. I would try to reduce it to 1e-6, 1e-7 etc. and see if it works. To compute an approximate of the time step see eq. (32) in the attached paper (CFL can be smaller than 0.1 sometimes). Best wishes, Margarida On Fri, Apr 16, 2021 at 1:05 AM Octavio Morales <oct_amg1289@hotmail.com> wrote:
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 <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.
_______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/nektar-users
Thank you for the information and the help Margarida, you were right! Best Wishes, Octavio A. Morales ________________________________ De: Margarida Moragues Ginard <margarida.moragues@gmail.com> Enviado: viernes, 16 de abril de 2021 01:13 a. m. Para: Octavio Morales <oct_amg1289@hotmail.com> CC: nektar-users@imperial.ac.uk <nektar-users@imperial.ac.uk> Asunto: Re: [Nektar-users] Compressible flow solver Hi Octavio, It could be that you need to decrease the time step size. Sometimes the required time step is very small when using explicit time integration. I would try to reduce it to 1e-6, 1e-7 etc. and see if it works. To compute an approximate of the time step see eq. (32) in the attached paper (CFL can be smaller than 0.1 sometimes). Best wishes, Margarida On Fri, Apr 16, 2021 at 1:05 AM Octavio Morales <oct_amg1289@hotmail.com<mailto:oct_amg1289@hotmail.com>> wrote: This email from oct_amg1289@hotmail.com<mailto: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<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. _______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk<mailto:Nektar-users@imperial.ac.uk> https://mailman.ic.ac.uk/mailman/listinfo/nektar-users
participants (2)
-
Margarida Moragues Ginard
-
Octavio Morales