Commit 50bf09e7 authored by Gianmarco Mengaldo's avatar Gianmarco Mengaldo

Regression test results modified accordingly to LDGNS fix. 4 parallel RT added.

parent 527f1587
......@@ -24,10 +24,6 @@ IF( NEKTAR_SOLVER_COMPRESSIBLE_FLOW )
ADD_SOLVER_EXECUTABLE(CompressibleFlowSolver solvers
${CompressibleFlowSolverSource})
ADD_NEKTAR_TEST (CylinderSubsonic_FRDG_SEM)
ADD_NEKTAR_TEST (CylinderSubsonic_FRDG_GAUSS)
......@@ -76,6 +72,10 @@ IF( NEKTAR_SOLVER_COMPRESSIBLE_FLOW )
#ADD_NEKTAR_TEST_LENGTHY(Perturbation_NS_M05_square_CBC_par)
#ADD_NEKTAR_TEST_LENGTHY(Perturbation_NS_M15_square_CBC_par)
ADD_NEKTAR_TEST(RAE5240_BSF_LE_bcsfromfile_par)
ADD_NEKTAR_TEST(CylinderSubsonic_NS_WeakDG_LDG_SEM_VariableMu_par)
ADD_NEKTAR_TEST(Couette_FRDG_LDG_GAUSS_par)
ADD_NEKTAR_TEST(Couette_WeakDG_LDG_SEM_VariableMu_par)
ADD_NEKTAR_TEST(Couette_FRHU_LDG_SEM_par)
ENDIF (NEKTAR_USE_MPI)
IF (NEKTAR_USE_FFTW)
......
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>NS, Couette flow, mixed bcs, FRHU advection and LDG diffusion, SEM, parallel</description>
<executable>CompressibleFlowSolver</executable>
<parameters>Couette_FRDG_LDG_GAUSS_par.xml</parameters>
<processes>3</processes>
<files>
<file description="Session File">Couette_FRDG_LDG_GAUSS_par.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0885112</value>
<value variable="rhou" tolerance="1e-12">60.3078</value>
<value variable="rhov" tolerance="1e-8">0.212138</value>
<value variable="E" tolerance="1e-12">4868.44</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0738533</value>
<value variable="rhou" tolerance="1e-12">60.9698</value>
<value variable="rhov" tolerance="2e-6">0.226637</value>
<value variable="E" tolerance="1e-12">4360.52</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0">-1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="1">-5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
<V ID="2">-5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="3">-1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="4">-5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
<V ID="5">-1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="6">-2.08166817e-12 0.00000000e+00 0.00000000e+00</V>
<V ID="7">-2.77555797e-17 5.00000000e-01 0.00000000e+00</V>
<V ID="8">2.08166817e-12 1.00000000e+00 0.00000000e+00</V>
<V ID="9">5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
<V ID="10">5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="11">5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
<V ID="12">1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="13">1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="14">1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 2 3 </E>
<E ID="3"> 3 0 </E>
<E ID="4"> 2 4 </E>
<E ID="5"> 4 5 </E>
<E ID="6"> 5 3 </E>
<E ID="7"> 1 6 </E>
<E ID="8"> 6 7 </E>
<E ID="9"> 7 2 </E>
<E ID="10"> 7 8 </E>
<E ID="11"> 8 4 </E>
<E ID="12"> 6 9 </E>
<E ID="13"> 9 10 </E>
<E ID="14"> 10 7 </E>
<E ID="15"> 10 11 </E>
<E ID="16"> 11 8 </E>
<E ID="17"> 9 12 </E>
<E ID="18"> 12 13 </E>
<E ID="19"> 13 10 </E>
<E ID="20"> 13 14 </E>
<E ID="21"> 14 11 </E>
</EDGE>
<ELEMENT>
<Q ID="0"> 0 1 2 3 </Q>
<Q ID="1"> 2 4 5 6 </Q>
<Q ID="2"> 7 8 9 1 </Q>
<Q ID="3"> 9 10 11 4 </Q>
<Q ID="4"> 12 13 14 8 </Q>
<Q ID="5"> 14 15 16 10 </Q>
<Q ID="6"> 17 18 19 13 </Q>
<Q ID="7"> 19 20 21 15 </Q>
</ELEMENT>
<COMPOSITE>
<C ID="0"> Q[0-7] </C>
<C ID="100"> E[0,7,12,17] </C>
<C ID="200"> E[18,20] </C>
<C ID="300"> E[5,11,16,21] </C>
<C ID="400"> E[3,6] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="3" FIELDS="rho,rhou,rhov,E" TYPE="GAUSS_LAGRANGE" />
</EXPANSIONS>
<CONDITIONS>
<!-- Castonguay Test-Case -->
<!-- M = 0.2, Re = 200, Pr = 0.72 -->
<PARAMETERS>
<P> TimeStep = 0.0001 </P>
<P> NumSteps = 200 </P>
<P> FinTime = 0 </P>
<P> IO_CheckSteps = 500000 </P>
<P> IO_InfoSteps = 500000 </P>
<P> GasConstant = 287.058 </P>
<P> Gamma = 1.4 </P>
<P> pInf = 101325 </P>
<P> rhoInf = 1.225 </P>
<P> uInf = 68.0588 </P>
<P> vInf = 0.0 </P>
<P> Twall = 300.15 </P>
<P> mu = 0.4169 </P>
<P> thermalConductivity = 581.6936 </P>
</PARAMETERS>
<SOLVERINFO>
<I PROPERTY="EQType" VALUE="NavierStokesCFE" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="AdvectionType" VALUE="FRDG" />
<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>
<VARIABLES>
<V ID="0"> rho </V>
<V ID="1"> rhou </V>
<V ID="2"> rhov </V>
<V ID="3"> E </V>
</VARIABLES>
<BOUNDARYREGIONS>
<B ID="0"> C[100] </B>
<B ID="1"> C[200] </B>
<B ID="2"> C[300] </B>
<B ID="3"> C[400] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="rho" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="rhou" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="rhov" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="E" USERDEFINEDTYPE="WallViscous" VALUE="0" />
</REGION>
<REGION REF="1">
<P VAR="rho" VALUE="[3]" />
<P VAR="rhou" VALUE="[3]" />
<P VAR="rhov" VALUE="[3]" />
<P VAR="E" VALUE="[3]" />
</REGION>
<REGION REF="2">
<D VAR="rho" VALUE="rhoInf" />
<D VAR="rhou" VALUE="rhoInf * uInf" />
<D VAR="rhov" VALUE="rhoInf * vInf" />
<D VAR="E" VALUE="1.05 * Twall * rhoInf * GasConstant / (Gamma - 1)" />
</REGION>
<REGION REF="3">
<P VAR="rho" VALUE="[1]" />
<P VAR="rhou" VALUE="[1]" />
<P VAR="rhov" VALUE="[1]" />
<P VAR="E" VALUE="[1]" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="InitialConditions">
<E VAR="rho" VALUE="rhoInf"/>
<E VAR="rhou" VALUE="rhoInf * uInf" />
<E VAR="rhov" VALUE="rhoInf * vInf" />
<E VAR="E" VALUE="pInf / (Gamma - 1) + 0.5 * rhoInf * (uInf * uInf + vInf * vInf)"/>
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="rho" VALUE="rhoInf" />
<E VAR="rhou" VALUE="rhoInf * uInf * y" />
<E VAR="rhov" VALUE="rhoInf * vInf" />
<E VAR="E" VALUE="109100 * ((1 / (Gamma - 1)) + (uInf * uInf * y * y / (2 * GasConstant)) / (Twall + 0.05 * y * Twall + 0.72 * uInf * uInf * y * (1 - y) / (2 * 1004.7)))" />
</FUNCTION>
</CONDITIONS>
</NEKTAR>
......@@ -9,15 +9,15 @@
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0889271</value>
<value variable="rhou" tolerance="1e-12">62.1128</value>
<value variable="rhov" tolerance="1e-8">0.175956</value>
<value variable="E" tolerance="1e-12">4905.04</value>
<value variable="rhou" tolerance="1e-12">62.1008</value>
<value variable="rhov" tolerance="1e-8">0.17592</value>
<value variable="E" tolerance="1e-12">4904.01</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0760154</value>
<value variable="rhou" tolerance="1e-12">56.0464</value>
<value variable="rhov" tolerance="2e-6">0.265763</value>
<value variable="E" tolerance="1e-12">4381.12</value>
<value variable="rho" tolerance="1e-12">0.076032</value>
<value variable="rhou" tolerance="1e-12">56.1036</value>
<value variable="rhov" tolerance="2e-6">0.265578</value>
<value variable="E" tolerance="1e-12">4361.42</value>
</metric>
</metrics>
</test>
......
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>NS, Couette flow, mixed bcs, FRHU advection and LDG diffusion, SEM, parallel</description>
<executable>CompressibleFlowSolver</executable>
<parameters>Couette_FRHU_LDG_SEM_par.xml</parameters>
<processes>2</processes>
<files>
<file description="Session File">Couette_FRHU_LDG_SEM_par.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0889265</value>
<value variable="rhou" tolerance="1e-12">62.1008</value>
<value variable="rhov" tolerance="1e-8">0.17592</value>
<value variable="E" tolerance="1e-12">4904.01</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.076032</value>
<value variable="rhou" tolerance="1e-12">56.1036</value>
<value variable="rhov" tolerance="2e-6">0.265578</value>
<value variable="E" tolerance="1e-12">4361.42</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0">-1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="1">-5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
<V ID="2">-5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="3">-1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="4">-5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
<V ID="5">-1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="6">-2.08166817e-12 0.00000000e+00 0.00000000e+00</V>
<V ID="7">-2.77555797e-17 5.00000000e-01 0.00000000e+00</V>
<V ID="8">2.08166817e-12 1.00000000e+00 0.00000000e+00</V>
<V ID="9">5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
<V ID="10">5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="11">5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
<V ID="12">1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="13">1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="14">1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 2 3 </E>
<E ID="3"> 3 0 </E>
<E ID="4"> 2 4 </E>
<E ID="5"> 4 5 </E>
<E ID="6"> 5 3 </E>
<E ID="7"> 1 6 </E>
<E ID="8"> 6 7 </E>
<E ID="9"> 7 2 </E>
<E ID="10"> 7 8 </E>
<E ID="11"> 8 4 </E>
<E ID="12"> 6 9 </E>
<E ID="13"> 9 10 </E>
<E ID="14"> 10 7 </E>
<E ID="15"> 10 11 </E>
<E ID="16"> 11 8 </E>
<E ID="17"> 9 12 </E>
<E ID="18"> 12 13 </E>
<E ID="19"> 13 10 </E>
<E ID="20"> 13 14 </E>
<E ID="21"> 14 11 </E>
</EDGE>
<ELEMENT>
<Q ID="0"> 0 1 2 3 </Q>
<Q ID="1"> 2 4 5 6 </Q>
<Q ID="2"> 7 8 9 1 </Q>
<Q ID="3"> 9 10 11 4 </Q>
<Q ID="4"> 12 13 14 8 </Q>
<Q ID="5"> 14 15 16 10 </Q>
<Q ID="6"> 17 18 19 13 </Q>
<Q ID="7"> 19 20 21 15 </Q>
</ELEMENT>
<COMPOSITE>
<C ID="0"> Q[0-7] </C>
<C ID="100"> E[0,7,12,17] </C>
<C ID="200"> E[18,20] </C>
<C ID="300"> E[5,11,16,21] </C>
<C ID="400"> E[3,6] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="3" FIELDS="rho,rhou,rhov,E" TYPE="GLL_LAGRANGE_SEM" />
</EXPANSIONS>
<CONDITIONS>
<!-- Castonguay Test-Case -->
<!-- M = 0.2, Re = 200, Pr = 0.72 -->
<PARAMETERS>
<P> TimeStep = 0.0001 </P>
<P> NumSteps = 200 </P>
<P> FinTime = 0 </P>
<P> IO_CheckSteps = 500000 </P>
<P> IO_InfoSteps = 500000 </P>
<P> GasConstant = 287.058 </P>
<P> Gamma = 1.4 </P>
<P> pInf = 101325 </P>
<P> rhoInf = 1.225 </P>
<P> uInf = 68.0588 </P>
<P> vInf = 0.0 </P>
<P> Twall = 300.15 </P>
<P> mu = 0.4169 </P>
<P> thermalConductivity = 581.6936 </P>
</PARAMETERS>
<SOLVERINFO>
<I PROPERTY="EQType" VALUE="NavierStokesCFE" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="AdvectionType" VALUE="FRHU" />
<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>
<VARIABLES>
<V ID="0"> rho </V>
<V ID="1"> rhou </V>
<V ID="2"> rhov </V>
<V ID="3"> E </V>
</VARIABLES>
<BOUNDARYREGIONS>
<B ID="0"> C[100] </B>
<B ID="1"> C[200] </B>
<B ID="2"> C[300] </B>
<B ID="3"> C[400] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="rho" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="rhou" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="rhov" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="E" USERDEFINEDTYPE="WallViscous" VALUE="0" />
</REGION>
<REGION REF="1">
<P VAR="rho" VALUE="[3]" />
<P VAR="rhou" VALUE="[3]" />
<P VAR="rhov" VALUE="[3]" />
<P VAR="E" VALUE="[3]" />
</REGION>
<REGION REF="2">
<D VAR="rho" VALUE="rhoInf" />
<D VAR="rhou" VALUE="rhoInf * uInf" />
<D VAR="rhov" VALUE="rhoInf * vInf" />
<D VAR="E" VALUE="1.05 * Twall * rhoInf * GasConstant / (Gamma - 1)" />
</REGION>
<REGION REF="3">
<P VAR="rho" VALUE="[1]" />
<P VAR="rhou" VALUE="[1]" />
<P VAR="rhov" VALUE="[1]" />
<P VAR="E" VALUE="[1]" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="InitialConditions">
<E VAR="rho" VALUE="rhoInf"/>
<E VAR="rhou" VALUE="rhoInf * uInf" />
<E VAR="rhov" VALUE="rhoInf * vInf" />
<E VAR="E" VALUE="pInf / (Gamma - 1) + 0.5 * rhoInf * (uInf * uInf + vInf * vInf)"/>
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="rho" VALUE="rhoInf" />
<E VAR="rhou" VALUE="rhoInf * uInf * y" />
<E VAR="rhov" VALUE="rhoInf * vInf" />
<E VAR="E" VALUE="109100 * ((1 / (Gamma - 1)) + (uInf * uInf * y * y / (2 * GasConstant)) / (Twall + 0.05 * y * Twall + 0.72 * uInf * uInf * y * (1 - y) / (2 * 1004.7)))" />
</FUNCTION>
</CONDITIONS>
</NEKTAR>
......@@ -8,16 +8,16 @@
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0889271</value>
<value variable="rhou" tolerance="1e-12">62.1128</value>
<value variable="rhov" tolerance="1e-8">0.175956</value>
<value variable="E" tolerance="1e-12">4905.04</value>
<value variable="rho" tolerance="1e-12">0.0889265</value>
<value variable="rhou" tolerance="1e-12">62.1008</value>
<value variable="rhov" tolerance="1e-8">0.17592</value>
<value variable="E" tolerance="1e-12">4904.01</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0760154</value>
<value variable="rhou" tolerance="1e-12">56.0464</value>
<value variable="rhov" tolerance="2e-6">0.265763</value>
<value variable="E" tolerance="1e-12">4381.12</value>
<value variable="rho" tolerance="1e-12">0.076032</value>
<value variable="rhou" tolerance="1e-12">56.1036</value>
<value variable="rhov" tolerance="2e-6">0.265578</value>
<value variable="E" tolerance="1e-12">4361.42</value>
</metric>
</metrics>
</test>
......
......@@ -11,14 +11,14 @@
<value variable="rho" tolerance="1e-12">0.000397706</value>
<value variable="rhou" tolerance="1e-12">48.1295</value>
<value variable="rhov" tolerance="1e-12">0.145836</value>
<value variable="rhow" tolerance="1e-12">8.79652e-06</value>
<value variable="rhow" tolerance="1e-8">8.79652e-06</value>
<value variable="E" tolerance="1e-12">17519.9</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.00139611</value>
<value variable="rhou" tolerance="1e-12">83.3516</value>
<value variable="rhov" tolerance="1e-12">0.505196</value>
<value variable="rhow" tolerance="1e-12">3.177e-05</value>
<value variable="rhow" tolerance="1e-8">3.177e-05</value>
<value variable="E" tolerance="1e-12">18953</value>
</metric>
</metrics>
......
......@@ -11,14 +11,14 @@
<value variable="rho" tolerance="1e-12">0.000397706</value>
<value variable="rhou" tolerance="1e-12">48.1295</value>
<value variable="rhov" tolerance="1e-12">0.145836</value>
<value variable="rhow" tolerance="1e-12">8.79652e-06</value>
<value variable="rhow" tolerance="1e-12">8.79677e-06</value>
<value variable="E" tolerance="1e-12">17519.9</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.00139611</value>
<value variable="rhou" tolerance="1e-12">83.3516</value>
<value variable="rhov" tolerance="1e-12">0.505196</value>
<value variable="rhow" tolerance="1e-12">3.177e-05</value>
<value variable="rhow" tolerance="1e-12">3.17716e-05</value>
<value variable="E" tolerance="1e-12">18953</value>
</metric>
</metrics>
......
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>NS, Couette flow, mixed bcs, WeakDG advection and LDG diffusion, Variable Viscosity, parallel</description>
<executable>CompressibleFlowSolver</executable>
<parameters>Couette_WeakDG_LDG_SEM_VariableMu_par.xml</parameters>
<processes>3</processes>
<files>
<file description="Session File">Couette_WeakDG_LDG_SEM_VariableMu_par.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0809477</value>
<value variable="rhou" tolerance="1e-12">51.9087</value>
<value variable="rhov" tolerance="1e-8">0.221151</value>
<value variable="E" tolerance="1e-12">4414.9</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0724519</value>
<value variable="rhou" tolerance="1e-12">55.3223</value>
<value variable="rhov" tolerance="1e-8">0.425157</value>
<value variable="E" tolerance="1e-12">4226.24</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0">-1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="1">-5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
<V ID="2">-5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="3">-1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="4">-5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
<V ID="5">-1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="6">-2.08166817e-12 0.00000000e+00 0.00000000e+00</V>
<V ID="7">-2.77555797e-17 5.00000000e-01 0.00000000e+00</V>
<V ID="8">2.08166817e-12 1.00000000e+00 0.00000000e+00</V>
<V ID="9">5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
<V ID="10">5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="11">5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
<V ID="12">1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="13">1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="14">1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 2 3 </E>
<E ID="3"> 3 0 </E>
<E ID="4"> 2 4 </E>
<E ID="5"> 4 5 </E>
<E ID="6"> 5 3 </E>
<E ID="7"> 1 6 </E>
<E ID="8"> 6 7 </E>
<E ID="9"> 7 2 </E>
<E ID="10"> 7 8 </E>
<E ID="11"> 8 4 </E>
<E ID="12"> 6 9 </E>
<E ID="13"> 9 10 </E>
<E ID="14"> 10 7 </E>
<E ID="15"> 10 11 </E>
<E ID="16"> 11 8 </E>
<E ID="17"> 9 12 </E>
<E ID="18"> 12 13 </E>
<E ID="19"> 13 10 </E>
<E ID="20"> 13 14 </E>
<E ID="21"> 14 11 </E>
</EDGE>
<ELEMENT>
<Q ID="0"> 0 1 2 3 </Q>
<Q ID="1"> 2 4 5 6 </Q>
<Q ID="2"> 7 8 9 1 </Q>
<Q ID="3"> 9 10 11 4 </Q>
<Q ID="4"> 12 13 14 8 </Q>
<Q ID="5"> 14 15 16 10 </Q>
<Q ID="6"> 17 18 19 13 </Q>
<Q ID="7"> 19 20 21 15 </Q>
</ELEMENT>
<COMPOSITE>
<C ID="0"> Q[0-7] </C>
<C ID="100"> E[0,7,12,17] </C>
<C ID="200"> E[18,20] </C>
<C ID="300"> E[5,11,16,21] </C>
<C ID="400"> E[3,6] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
<EXPANSIONS>
<E COMPOSITE="C[0]" FIELDS="rho,rhou,rhov,E" BASISTYPE="GLL_Lagrange,GLL_Lagrange" NUMMODES="3,3"
POINTSTYPE="GaussLobattoLegendre,GaussLobattoLegendre" NUMPOINTS="4,4" />
</EXPANSIONS>
<CONDITIONS>
<!-- Castonguay Test-Case -->
<!-- M = 0.2, Re = 200, Pr = 0.72 -->
<PARAMETERS>
<P> TimeStep = 0.0001 </P>
<P> FinTime = 0.05 </P>
<P> NumSteps = FinTime/TimeStep </P>
<P> IO_CheckSteps = 500000 </P>
<P> IO_InfoSteps = 1000 </P>
<P> GasConstant = 287.058 </P>
<P> Gamma = 1.4 </P>
<P> pInf = 101325 </P>
<P> rhoInf = 1.225 </P>
<P> uInf = 68.0588 </P>
<P> vInf = 0.0 </P>
<P> Twall = 300.15 </P>
<P> mu = 0.4169 </P>
<P> thermalConductivity = 581.6936 </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="RungeKutta2_ImprovedEuler"/>
<I PROPERTY="UpwindType" VALUE="ExactToro" />
<I PROPERTY="ProblemType" VALUE="General" />
<I PROPERTY="ViscosityType" VALUE="Variable" />
</SOLVERINFO>
<VARIABLES>
<V ID="0"> rho </V>
<V ID="1"> rhou </V>
<V ID="2"> rhov </V>
<V ID="3"> E </V>
</VARIABLES>
<BOUNDARYREGIONS>
<B ID="0"> C[100] </B>
<B ID="1"> C[200] </B>
<B ID="2"> C[300] </B>
<B ID="3"> C[400] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="rho" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="rhou" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="rhov" USERDEFINEDTYPE="WallViscous" VALUE="0" />
<D VAR="E" USERDEFINEDTYPE="WallViscous" VALUE="0" />
</REGION>
<REGION REF="1">
<P VAR="rho" VALUE="[3]" />
<P VAR="rhou" VALUE="[3]" />
<P VAR="rhov" VALUE="[3]" />
<P VAR="E" VALUE="[3]" />
</REGION>
<REGION REF="2">
<D VAR="rho" VALUE="rhoInf" />
<D VAR="rhou" VALUE="rhoInf * uInf" />
<D VAR="rhov" VALUE="rhoInf * vInf" />
<D VAR="E" VALUE="1.05 * Twall * rhoInf * GasConstant / (Gamma - 1)" />
</REGION>
<REGION REF="3">
<P VAR="rho" VALUE="[1]" />
<P VAR="rhou" VALUE="[1]" />
<P VAR="rhov" VALUE="[1]" />
<P VAR="E" VALUE="[1]" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="InitialConditions">
<E VAR="rho" VALUE="rhoInf"/>
<E VAR="rhou" VALUE="rhoInf * uInf" />
<E VAR="rhov" VALUE="rhoInf * vInf" />
<E VAR="E" VALUE="pInf / (Gamma - 1) + 0.5 * rhoInf * (uInf * uInf + vInf * vInf)"/>
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="rho" VALUE="rhoInf" />
<E VAR="rhou" VALUE="rhoInf * uInf * y" />
<E VAR="rhov" VALUE="rhoInf * vInf" />
<E VAR="E" VALUE="109100 * ((1 / (Gamma - 1)) + (uInf * uInf * y * y / (2 * GasConstant)) / (Twall + 0.05 * y * Twall + 0.72 * uInf * uInf * y * (1 - y) / (2 * 1004.7)))" />
</FUNCTION>
</CONDITIONS>
</NEKTAR>
......@@ -16,7 +16,7 @@
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.332458</value>
<value variable="rhou" tolerance="1e-12">87.3753</value>