<I PROPERTY="SolverType"     VALUE="VelocityCorrectionScheme" />

            <I PROPERTY="AdvectionForm"  VALUE="Convective" />

            <I PROPERTY="SpectralhpDealiasing"  VALUE="True" />

            <I PROPERTY=“SpectralVanishingViscosity”  VALUE="DGKernel" />

I not sure wether Dealiasiing ON works (you should see if it is mentioned in the summary when you run)

I also think you may need Convective form advection (rather than Skey Symmetric for Dealiasing. Finally you may also require to have to increase the quadrature points to allow for geometric dealiasing. I am cc’ing Walid who has been running 3D most recently.

Finally you might want to have a lower polynomial order to start with.

Cheers,

Spencer

PS Are you also using the LowEnergyBlock preconditioners? 

On 23 Mar 2019, at 16:06, Moto Light <redamansyhh@gmail.com> wrote:

Hi All,

i am runing simple case at Re1000, i tried many meshs all of them give me the same error

//////////////////////

Steps: 91200    Time: 424.694      CPU Time: 4.48716s

Steps: 91210    Time: 424.695      CPU Time: 4.62339s

CG iterations made = 5001 using tolerance of 1e-09 (error = nan, rhs_mag = inf)

job ended at

/////////////////////////////

any suggestions?

my case file

///////////////////////////////////////////

<?xml version="1.0" encoding="UTF-8" ?>

<NEKTAR>

    <EXPANSIONS>

    <E COMPOSITE="C[1]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- 1 Domain inlet plate upper -->

    <E COMPOSITE="C[2]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- 2 Domain inlet plate lower -->

    <E COMPOSITE="C[3]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- ****3 Domain square faces up down front meshs**** -->

    <E COMPOSITE="C[4]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- 4 Domain upper domain lower in up stream -->

    <E COMPOSITE="C[5]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- 5 Domain upper lower downstream  -->

    <E COMPOSITE="C[6]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- ***6 Domain down stream back side square*** -->

    <E COMPOSITE="C[7]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- 7 Domain down stream after square big one -->

    <E COMPOSITE="C[8]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />        <!-- 7 Domain down stream after square big one -->

   </EXPANSIONS>

  <CONDITIONS>

    <SOLVERINFO>

      <I PROPERTY="SOLVERTYPE"              VALUE="VelocityCorrectionScheme"/>

      <I PROPERTY="EQTYPE"                  VALUE="UnsteadyNavierStokes"    />

      <I PROPERTY="Projection"              VALUE="Continuous"              />

      <I PROPERTY="EvolutionOperator"       VALUE="Nonlinear"               />

      <I PROPERTY="TimeIntegrationMethod"   VALUE="IMEXOrder2"              />

      <I PROPERTY="Driver"                  VALUE="Standard"                />

      <I PROPERTY="HOMOGENEOUS"             VALUE="1D"     />

      <I PROPERTY="AdvectionForm"           VALUE="SkewSymmetric"           /> 

      <I PROPERTY="Extrapolation"           VALUE="Standard"                />

      <I PROPERTY="USEFFT"                  VALUE="FFTW"                    />

 <I PROPERTY="DEALIASING"           VALUE="ON"                   />  

   <!--  <I PROPERTY="SpectralVanishingViscosity"   VALUE="True"                    /> -->     

  </SOLVERINFO> 

    <PARAMETERS>

      <P> TimeStep      = 0.0001     </P>

      <P> NumSteps      = 2500000    </P>  <!-- 500000  250sec -->

      <P> IO_CheckSteps = 5000       </P>  <!-- 120 check files-->

      <P> IO_InfoSteps  = 10         </P>

      <P> Re            = 200         </P>

      <P> Kinvis        = 1./Re      </P>

      <P> IO_CFLSteps   = 50         </P>

      <P> LZ    = 9        </P>  <!-- 3. * Fastest spanwise     -->

      <P> HomModesZ   = 32      </P>  <!-- how many cells for Z direction     -->

    </PARAMETERS>

    

    <VARIABLES>

      <V ID="0"> u </V> 

      <V ID="1"> v </V> 

      <V ID="2"> w </V>

      <V ID="3"> p </V> 

    </VARIABLES>

  <BOUNDARYREGIONS>

      <B ID="0"> C[9]   </B> <!-- Inlet1_upper  -->

      <B ID="1"> C[10]  </B> <!-- Inlet2_Lower  -->

      <B ID="2"> C[11]  </B> <!-- Splitter plate  -->

      <B ID="3"> C[12]  </B> <!-- Upper Boundary  -->

      <B ID="4"> C[13]  </B> <!-- Lower_Boundary  -->

      <B ID="5"> C[14]  </B> <!-- Outlet  -->

      <B ID="6"> C[15]  </B> <!-- Wall_square   -->

   </BOUNDARYREGIONS>

 <BOUNDARYCONDITIONS>

      <REGION REF="0">          <!-- Inlet_upper  -->

        <D VAR="u" VALUE="5" />

        <D VAR="v" VALUE="0" />

<D VAR="w" VALUE="0" />

        <N VAR="p" USERDEFINEDTYPE="H" VALUE="0"/>  

      </REGION>

      <REGION REF="1">          <!-- Inlet_lower  -->

        <D VAR="u" VALUE="1" />

        <D VAR="v" VALUE="0" />

<D VAR="w" VALUE="0" />

        <N VAR="p" USERDEFINEDTYPE="H" VALUE="0"/>  

      </REGION>

       <REGION REF="2">           <!--Splitter plate  -->

        <D VAR="u" VALUE="0" />

        <D VAR="v" VALUE="0" />

<D VAR="w" VALUE="0" />

        <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />  

      </REGION>

      <REGION REF="3">     <!-- Upper Boundary  -->

        <N VAR="u" USERDEFINEDTYPE="H" VALUE="0" />

        <D VAR="v" VALUE="0" />

<N VAR="w" USERDEFINEDTYPE="H" VALUE="0" />

        <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />  

      </REGION>

       <REGION REF="4">     <!-- Lower_Boundary  -->

        <N VAR="u" USERDEFINEDTYPE="H" VALUE="0" />

        <D VAR="v" VALUE="0" />

<N VAR="w" USERDEFINEDTYPE="H" VALUE="0" />

        <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />  

      </REGION>

       <REGION REF="5">     <!-- Outlet  BC -->

        <N VAR="u" USERDEFINEDTYPE="HOutflow" VALUE="0" />

        <N VAR="v" USERDEFINEDTYPE="HOutflow" VALUE="0" />

<N VAR="w" USERDEFINEDTYPE="HOutflow" VALUE="0" />

        <D VAR="p" USERDEFINEDTYPE="HOutflow" VALUE="0" />  

      </REGION>

       <REGION REF="6">      <!-- WallA_square  -->

        <D VAR="u" VALUE="0" />

        <D VAR="v" VALUE="0" />

<D VAR="w" VALUE="0" />

        <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />  

      </REGION>

    </BOUNDARYCONDITIONS>

    <!-- We initialise from pre-computed flow for speed -->

    <FUNCTION NAME="InitialConditions">

         <E VAR="w" VALUE="awgn(1e-3)*cos(2*PI*z/LZ)" />    

      <!--  <F VAR="u,v,p,w" FILE="SquareN3D_179.chk" /> -->  

       <!-- <E VAR="u,v,w,p" VALUE="0" />   -->

         <F VAR="u,v,p" FILE="SquareN2D_200.chk" />                 <!-- start from 2D results    -->

    </FUNCTION>

    <FUNCTION NAME="ExactSolution">

      <E VAR="u" VALUE="0" />

      <E VAR="v" VALUE="0" />

      <E VAR="w" VALUE="0" />

      <E VAR="p" VALUE="0" />

    </FUNCTION>

  </CONDITIONS>

    

  <FILTERS>

     <FILTER TYPE="AeroForces">

      <PARAM NAME="OutputFile">stats/stats1/CylinderForces3D</PARAM>

      <PARAM NAME="OutputFrequency">10</PARAM>

      <PARAM NAME="Boundary"> B[6] </PARAM>

    </FILTER>

</FILTERS>

</NEKTAR>

/////////////////////////////////

please help me in this

Regrads,

cheers,

Reda PhD student 

Aerospace department 

UPC Barcelona

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

Spencer Sherwin FREng, FRAeS
Head, Aerodynamics,
Director of Research Computing Service,
Professor of Computational Fluid Mechanics,
Department of Aeronautics,
s.sherwin@imperial.ac.uk                                   South Kensington Campus,
Phone: +44 (0)20 7594 5052                              Imperial College London,
Fax:   +44 (0)20 7594 1974                               London, SW7 2AZ,  UK
http://www.imperial.ac.uk/people/s.sherwin/