Commit 2da0e078 authored by Dave Moxey's avatar Dave Moxey

Remove problematic Couette flow cases and add an exact case that gives more...

Remove problematic Couette flow cases and add an exact case that gives more consistent results in parallel. Update tutorial subrepo to latest master.
parent 0556c661
Pipeline #1211 failed with stages
in 13 minutes and 12 seconds
Subproject commit 851f7583b288fb264e1d95db647eb541c054962a
Subproject commit 2b3a87f9ad78bf4457ba6906725901a054fa66bc
......@@ -119,10 +119,7 @@ IF( NEKTAR_SOLVER_COMPRESSIBLE_FLOW )
#ADD_NEKTAR_TEST(Perturbation_NS_M05_square_CBC_par LENGTHY)
#ADD_NEKTAR_TEST(Perturbation_NS_M15_square_CBC_par LENGTHY)
ADD_NEKTAR_TEST(RAE5240_BSF_LE_bcsfromfile_par)
ADD_NEKTAR_TEST(Couette_FRDG_LDG_GAUSS_par)
ADD_NEKTAR_TEST(Couette_FRHU_LDG_SEM_par)
ADD_NEKTAR_TEST(Couette_WeakDG_LDG_SEM_VariableMu_par)
ADD_NEKTAR_TEST(CylinderSubsonic_NS_WeakDG_LDG_SEM_VariableMu_par)
ADD_NEKTAR_TEST(Couette_WeakDG_LDG_Exact_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>--use-scotch 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-7">0.0885179</value>
<value variable="rhou" tolerance="1e-4">60.3076</value>
<value variable="rhov" tolerance="1e-6">0.211622</value>
<value variable="E" tolerance="1e-2">4866.11</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-7">0.0738853</value>
<value variable="rhou" tolerance="1e-4">60.9894</value>
<value variable="rhov" tolerance="2e-6">0.226922</value>
<value variable="E" tolerance="1e-2">4357.43</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>
<?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>--use-scotch 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-7">0.0889264</value>
<value variable="rhou" tolerance="1e-4">62.0963</value>
<value variable="rhov" tolerance="1e-6">0.175956</value>
<value variable="E" tolerance="1e-2">4903.95</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-7">0.0760421</value>
<value variable="rhou" tolerance="1e-4">56.0835</value>
<value variable="rhov" tolerance="1e-6">0.265899</value>
<value variable="E" tolerance="1e-2">4333.92</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>
<?xml version="1.0" encoding="utf-8" ?>
<test>
<description>NS, Couette flow with periodic BCs, exact solution, parallel</description>
<executable>CompressibleFlowSolver</executable>
<parameters>-I parts=6,7:2,4,5:0,1,3 Couette_WeakDG_LDG_Exact_par.xml</parameters>
<processes>3</processes>
<files>
<file description="Session File">Couette_WeakDG_LDG_Exact_par.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">5.26e-11</value>
<value variable="rhou" tolerance="1e-12">1.08199e-09</value>
<value variable="rhov" tolerance="1e-12">1.13107e-09</value>
<value variable="E" tolerance="1e-10">1.06807e-06</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">1.70021e-10</value>
<value variable="rhou" tolerance="1e-12">9.76513e-10</value>
<value variable="rhov" tolerance="1e-12">9.81855e-10</value>
<value variable="E" tolerance="1e-10">8.31031e-07</value>
</metric>
</metrics>
</test>
......@@ -61,96 +61,93 @@
</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>
<EXPANSIONS>
<E COMPOSITE="C[0]" FIELDS="rho,rhou,rhov,E" BASISTYPE="GLL_Lagrange,GLL_Lagrange" NUMMODES="4,4"
POINTSTYPE="GaussLobattoLegendre,GaussLobattoLegendre" NUMPOINTS="5,5" />
</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>
<CONDITIONS>
<PARAMETERS>
<P> TimeStep = 0.00005 </P>
<P> FinTime = 0.01 </P>
<P> NumSteps = FinTime/TimeStep </P>
<P> IO_CheckSteps = 500000 </P>
<P> IO_InfoSteps = 100 </P>
<P> GasConstant = 287.1428571 </P>
<P> Gamma = 1.4 </P>
<P> Cp = 1005 </P>
<P> Uwall = 69.445 </P>
<P> Twall = 300 </P>
<P> pc = 100000 </P>
<P> Pr = 0.72 </P>
<P> mu = 0.417 </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>
<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="Constant" />
</SOLVERINFO>
<VARIABLES>
<V ID="0"> rho </V>
<V ID="1"> rhou </V>
<V ID="2"> rhov </V>
<V ID="3"> E </V>
</VARIABLES>
<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>
<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>
<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="Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall)" />
<D VAR="rhou" VALUE="(Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall)) * Uwall" />
<D VAR="rhov" VALUE="0" />
<D VAR="E" VALUE="pc / (Gamma - 1) + 0.5 * (Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall)) * (Uwall*Uwall)" />
</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>
<FUNCTION NAME="InitialConditions">
<E VAR="rho" VALUE="Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall + Pr * Uwall * Uwall * y * (1-y))" />
<E VAR="rhou" VALUE="(Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall + Pr * Uwall * Uwall * y * (1-y))) * Uwall * y" />
<E VAR="rhov" VALUE="0" />
<E VAR="E" VALUE="pc / (Gamma - 1) + 0.5 * (Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall + Pr * Uwall * Uwall * y * (1-y))) * (Uwall * Uwall * y * y)" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="rho" VALUE="Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall + Pr * Uwall * Uwall * y * (1-y))" />
<E VAR="rhou" VALUE="Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall + Pr * Uwall * Uwall * y * (1-y)) * Uwall * y" />
<E VAR="rhov" VALUE="0" />
<E VAR="E" VALUE="pc / (Gamma - 1) + 0.5 * (Gamma / (Gamma-1) * (2 * pc) / (2 * Cp * Twall + Pr * Uwall * Uwall * y * (1-y))) * (Uwall * Uwall * y * y)" />
</FUNCTION>
</CONDITIONS>
</NEKTAR>
<?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>--use-scotch 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-7">0.0805614</value>
<value variable="rhou" tolerance="1e-4">51.8867</value>
<value variable="rhov" tolerance="1e-6">0.222213</value>
<value variable="E" tolerance="1e-2">4415.64</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-7">0.0716895</value>
<value variable="rhou" tolerance="1e-4">55.3017</value>
<value variable="rhov" tolerance="1e-6 ">0.424398</value>
<value variable="E" tolerance="1e-2">4227.23</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>NS, Subsonic Cylinder, mixed bcs, WeakDG advection and LDG diffusion, variable viscosity, parallel</description>
<executable>CompressibleFlowSolver</executable>
<parameters>--use-scotch CylinderSubsonic_NS_WeakDG_LDG_SEM_VariableMu_par.xml</parameters>
<processes>3</processes>
<files>
<file description="Session File">CylinderSubsonic_NS_WeakDG_LDG_SEM_VariableMu_par.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="2e-6">5.50797</value>
<value variable="rhou" tolerance="2e-6">2016.24</value>
<value variable="rhov" tolerance="2e-4">25.8675</value>
<value variable="E" tolerance="2e-6">6.27025e+06</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="2e-6">0.530811</value>
<value variable="rhou" tolerance="2e-4">91.9676</value>
<value variable="rhov" tolerance="2e-4">58.1354</value>
<value variable="E" tolerance="2e-1">346847</value>
</metric>
</metrics>
</test>
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment