Commit 63734a5e authored by Dave Moxey's avatar Dave Moxey

Add IncNavierStokesSolver pyramid regression test

parent 9a04ce33
......@@ -42,6 +42,7 @@ IF( NEKTAR_SOLVER_INCNAVIERSTOKES )
ADD_NEKTAR_TEST(ChanFlow_m8_singular)
ADD_NEKTAR_TEST(Channel_Flow_3modes_rad)
ADD_NEKTAR_TEST(Hex_channel_m3)
ADD_NEKTAR_TEST(Pyr_channel_m3)
ADD_NEKTAR_TEST_LENGTHY(Hex_channel_m8)
ADD_NEKTAR_TEST_LENGTHY(Hex_channel_m8_srhs)
ADD_NEKTAR_TEST_LENGTHY(Hex_Kovasnay_SubStep)
......@@ -89,6 +90,7 @@ IF( NEKTAR_SOLVER_INCNAVIERSTOKES )
ADD_NEKTAR_TEST(ChanFlow_m3_par)
ADD_NEKTAR_TEST_LENGTHY(ChanFlow_m8_BodyForce_par)
ADD_NEKTAR_TEST_LENGTHY(Hex_channel_m8_par)
ADD_NEKTAR_TEST_LENGTHY(Pyr_channel_m6_par)
ADD_NEKTAR_TEST_LENGTHY(Hex_channel_m8_srhs_par)
ADD_NEKTAR_TEST(Tet_channel_m8_par)
ADD_NEKTAR_TEST_LENGTHY(Tet_channel_m8_iter_ml_par)
......
<?xml version="1.0" encoding="utf-8" ?>
<test>
<description>3D channel flow, Pyramidic elements, P=3</description>
<executable>IncNavierStokesSolver</executable>
<parameters>Pyr_channel_m3.xml</parameters>
<files>
<file description="Session File">Pyr_channel_m3.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12">0</value>
<value variable="v" tolerance="1e-12">0</value>
<value variable="w" tolerance="1e-12">0</value>
<value variable="p" tolerance="1e-8">1.16849e-14</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12">1.75668e-16</value>
<value variable="v" tolerance="1e-12">1.07627e-16</value>
<value variable="w" tolerance="1e-12">5.55112e-16</value>
<value variable="p" tolerance="1e-8">4.04121e-14</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="3" FIELDS="u,v,w,p" TYPE="MODIFIED" />
</EXPANSIONS>
<CONDITIONS>
<SOLVERINFO>
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="AdvectionForm" VALUE="Convective" />
<I PROPERTY="Projection" VALUE="Galerkin" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder1" />
</SOLVERINFO>
<PARAMETERS>
<P> TimeStep = 0.001 </P>
<P> NumSteps = 100 </P>
<P> IO_InfoSteps = 100 </P>
<P> Kinvis = 1 </P>
</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[1] </B> <!-- Inlet -->
<B ID="1"> C[6] </B> <!-- Outlet -->
<B ID="2"> C[2] </B> <!-- Wall -->
<B ID="3"> C[3] </B> <!-- Wall left -->
<B ID="4"> C[4] </B> <!-- Wall -->
<B ID="5"> C[5] </B> <!-- Wall right -->
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<D VAR="w" VALUE="y*(1-y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="1">
<N VAR="u" VALUE="0" />
<N VAR="v" VALUE="0" />
<N VAR="w" VALUE="0" />
<D VAR="p" VALUE="0" />
</REGION>
<REGION REF="2">
<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">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<D VAR="w" VALUE="y*(1-y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="4">
<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="5">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<D VAR="w" VALUE="y*(1-y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="w" VALUE="y*(1-y)" />
<E VAR="p" VALUE="-2*Kinvis*(z-1)" />
</FUNCTION>
<FUNCTION NAME="InitialConditions">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="w" VALUE="y*(1-y)" />
<E VAR="p" VALUE="-2*Kinvis*(z-1)" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="3" SPACE="3">
<VERTEX>
<V ID="0">0.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="1">1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="2">1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="3">0.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="4">5.00000000e-01 5.00000000e-01 5.00000000e-01</V>
<V ID="5">0.00000000e+00 0.00000000e+00 1.00000000e+00</V>
<V ID="6">1.00000000e+00 0.00000000e+00 1.00000000e+00</V>
<V ID="7">1.00000000e+00 1.00000000e+00 1.00000000e+00</V>
<V ID="8">0.00000000e+00 1.00000000e+00 1.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 3 2 </E>
<E ID="3"> 0 3 </E>
<E ID="4"> 0 4 </E>
<E ID="5"> 1 4 </E>
<E ID="6"> 2 4 </E>
<E ID="7"> 3 4 </E>
<E ID="8"> 0 5 </E>
<E ID="9"> 5 6 </E>
<E ID="10"> 1 6 </E>
<E ID="11"> 5 4 </E>
<E ID="12"> 6 4 </E>
<E ID="13"> 6 7 </E>
<E ID="14"> 2 7 </E>
<E ID="15"> 7 4 </E>
<E ID="16"> 7 8 </E>
<E ID="17"> 3 8 </E>
<E ID="18"> 8 4 </E>
<E ID="19"> 8 5 </E>
</EDGE>
<FACE>
<Q ID="0"> 0 1 2 3</Q>
<T ID="1"> 0 5 4</T>
<T ID="2"> 1 6 5</T>
<T ID="3"> 2 6 7</T>
<T ID="4"> 3 7 4</T>
<Q ID="5"> 8 9 10 0</Q>
<T ID="6"> 8 11 4</T>
<T ID="7"> 9 12 11</T>
<T ID="8"> 10 12 5</T>
<Q ID="9"> 10 13 14 1</Q>
<T ID="10"> 13 15 12</T>
<T ID="11"> 14 15 6</T>
<Q ID="12"> 14 16 17 2</Q>
<T ID="13"> 16 18 15</T>
<T ID="14"> 17 18 7</T>
<Q ID="15"> 17 19 8 3</Q>
<T ID="16"> 19 11 18</T>
<Q ID="17"> 19 16 13 9</Q>
</FACE>
<ELEMENT>
<P ID="0"> 0 1 2 3 4 </P>
<P ID="1"> 5 6 7 8 1 </P>
<P ID="2"> 9 8 10 11 2 </P>
<P ID="3"> 12 11 13 14 3 </P>
<P ID="4"> 15 14 16 6 4 </P>
<P ID="5"> 17 16 13 10 7 </P>
</ELEMENT>
<COMPOSITE>
<C ID="0"> P[0-5] </C>
<C ID="1"> F[0] </C>
<C ID="2"> F[5] </C>
<C ID="3"> F[9] </C>
<C ID="4"> F[12] </C>
<C ID="5"> F[15] </C>
<C ID="6"> F[17] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="6" FIELDS="u,v,w,p" TYPE="MODIFIED" />
</EXPANSIONS>
<CONDITIONS>
<SOLVERINFO>
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="AdvectionForm" VALUE="Convective" />
<I PROPERTY="Projection" VALUE="Galerkin" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder1" />
</SOLVERINFO>
<PARAMETERS>
<P> TimeStep = 0.001 </P>
<P> NumSteps = 10 </P>
<P> IO_InfoSteps = 100 </P>
<P> Kinvis = 1 </P>
</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[1] </B> <!-- Inlet -->
<B ID="1"> C[6] </B> <!-- Outlet -->
<B ID="2"> C[2] </B> <!-- Wall -->
<B ID="3"> C[3] </B> <!-- Wall left -->
<B ID="4"> C[4] </B> <!-- Wall -->
<B ID="5"> C[5] </B> <!-- Wall right -->
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<D VAR="w" VALUE="y*(1-y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="1">
<N VAR="u" VALUE="0" />
<N VAR="v" VALUE="0" />
<N VAR="w" VALUE="0" />
<D VAR="p" VALUE="0" />
</REGION>
<REGION REF="2">
<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">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<D VAR="w" VALUE="y*(1-y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="4">
<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="5">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<D VAR="w" VALUE="y*(1-y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="w" VALUE="y*(1-y)" />
<E VAR="p" VALUE="-2*Kinvis*(z-1)" />
</FUNCTION>
<FUNCTION NAME="InitialConditions">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="w" VALUE="y*(1-y)" />
<E VAR="p" VALUE="-2*Kinvis*(z-1)" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="3" SPACE="3">
<VERTEX>
<V ID="0">0.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="1">1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="2">1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="3">0.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="4">5.00000000e-01 5.00000000e-01 5.00000000e-01</V>
<V ID="5">0.00000000e+00 0.00000000e+00 1.00000000e+00</V>
<V ID="6">1.00000000e+00 0.00000000e+00 1.00000000e+00</V>
<V ID="7">1.00000000e+00 1.00000000e+00 1.00000000e+00</V>
<V ID="8">0.00000000e+00 1.00000000e+00 1.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 3 2 </E>
<E ID="3"> 0 3 </E>
<E ID="4"> 0 4 </E>
<E ID="5"> 1 4 </E>
<E ID="6"> 2 4 </E>
<E ID="7"> 3 4 </E>
<E ID="8"> 0 5 </E>
<E ID="9"> 5 6 </E>
<E ID="10"> 1 6 </E>
<E ID="11"> 5 4 </E>
<E ID="12"> 6 4 </E>
<E ID="13"> 6 7 </E>
<E ID="14"> 2 7 </E>
<E ID="15"> 7 4 </E>
<E ID="16"> 7 8 </E>
<E ID="17"> 3 8 </E>
<E ID="18"> 8 4 </E>
<E ID="19"> 8 5 </E>
</EDGE>
<FACE>
<Q ID="0"> 0 1 2 3</Q>
<T ID="1"> 0 5 4</T>
<T ID="2"> 1 6 5</T>
<T ID="3"> 2 6 7</T>
<T ID="4"> 3 7 4</T>
<Q ID="5"> 8 9 10 0</Q>
<T ID="6"> 8 11 4</T>
<T ID="7"> 9 12 11</T>
<T ID="8"> 10 12 5</T>
<Q ID="9"> 10 13 14 1</Q>
<T ID="10"> 13 15 12</T>
<T ID="11"> 14 15 6</T>
<Q ID="12"> 14 16 17 2</Q>
<T ID="13"> 16 18 15</T>
<T ID="14"> 17 18 7</T>
<Q ID="15"> 17 19 8 3</Q>
<T ID="16"> 19 11 18</T>
<Q ID="17"> 19 16 13 9</Q>
</FACE>
<ELEMENT>
<P ID="0"> 0 1 2 3 4 </P>
<P ID="1"> 5 6 7 8 1 </P>
<P ID="2"> 9 8 10 11 2 </P>
<P ID="3"> 12 11 13 14 3 </P>
<P ID="4"> 15 14 16 6 4 </P>
<P ID="5"> 17 16 13 10 7 </P>
</ELEMENT>
<COMPOSITE>
<C ID="0"> P[0-5] </C>
<C ID="1"> F[0] </C>
<C ID="2"> F[5] </C>
<C ID="3"> F[9] </C>
<C ID="4"> F[12] </C>
<C ID="5"> F[15] </C>
<C ID="6"> F[17] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<test>
<description>3D channel flow, Pyramidic elements, P=6</description>
<executable>IncNavierStokesSolver</executable>
<parameters>Pyr_channel_m6.xml</parameters>
<processes>3</processes>
<files>
<file description="Session File">Pyr_channel_m6.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12">3.77879e-14</value>
<value variable="v" tolerance="1e-12">4.29463e-14</value>
<value variable="w" tolerance="1e-12">3.96346e-13</value>
<value variable="p" tolerance="1e-8">5.15486e-12</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12">2.38835e-13</value>
<value variable="v" tolerance="1e-12">2.2169e-13</value>
<value variable="w" tolerance="1e-12">1.92338e-12</value>
<value variable="p" tolerance="1e-8">2.83573e-11</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