Commit 8c7a7c6a authored by Spencer Sherwin's avatar Spencer Sherwin
Browse files

Added regression test for 3D substepping


git-svn-id: https://gforge.sci.utah.edu/svn/nektar/trunk@4033 305cdda6-5ce1-45b3-a98d-dfc68c8b3305
parent 0f1169d1
......@@ -110,6 +110,8 @@ int main(int argc, char* argv[])
//Test Substepping scheme on Kovasnay solution
Execute("IncNavierStokesSolver","Test_KovaFlow_SubStep_2order.xml","Kovasnay solution using sub-stepping");
Execute("IncNavierStokesSolver","Test_Hex_Kovasnay_SubStep.xml","3D Hex Kovasnay solution using sub-stepping");
//Test Modified Arnoldi adjoint stability (VelCorrectionScheme)
Execute("IncNavierStokesSolver","ChanStability_adj.xml","Adjoint stability (Mod. Arnoldi): Channel");
......
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<GEOMETRY DIM="3" SPACE="3">
<VERTEX>
<V ID="0">-0.5000000e+00 -0.5000000e+00 0.00000000e+00</V>
<V ID="1">2.50000000e-01 -0.5000000e+00 0.00000000e+00</V>
<V ID="2">2.50000000e-01 5.00000000e-01 0.00000000e+00</V>
<V ID="3">-0.5000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="4">-0.5000000e+00 -0.5000000e+00 5.00000000e-01</V>
<V ID="5">2.50000000e-01 -0.5000000e+00 5.00000000e-01</V>
<V ID="6">2.50000000e-01 5.00000000e-01 5.00000000e-01</V>
<V ID="7">-0.5000000e+00 5.00000000e-01 5.00000000e-01</V>
<V ID="8">-0.5000000e+00 -0.5000000e+00 1.00000000e+00</V>
<V ID="9">2.50000000e-01 -0.5000000e+00 1.00000000e+00</V>
<V ID="10">2.50000000e-01 5.00000000e-01 1.00000000e+00</V>
<V ID="11">-0.5000000e+00 5.00000000e-01 1.00000000e+00</V>
<V ID="12">2.50000000e-01 1.50000000e+00 0.00000000e+00</V>
<V ID="13">-0.5000000e+00 1.50000000e+00 0.00000000e+00</V>
<V ID="14">2.50000000e-01 1.50000000e+00 5.00000000e-01</V>
<V ID="15">-0.5000000e+00 1.50000000e+00 5.00000000e-01</V>
<V ID="16">2.50000000e-01 1.50000000e+00 1.00000000e+00</V>
<V ID="17">-0.5000000e+00 1.50000000e+00 1.00000000e+00</V>
<V ID="18">1.00000000e+00 -0.5000000e+00 0.00000000e+00</V>
<V ID="19">1.00000000e+00 5.00000000e-01 0.00000000e+00</V>
<V ID="20">1.00000000e+00 -0.5000000e+00 5.00000000e-01</V>
<V ID="21">1.00000000e+00 5.00000000e-01 5.00000000e-01</V>
<V ID="22">1.00000000e+00 -0.5000000e+00 1.00000000e+00</V>
<V ID="23">1.00000000e+00 5.00000000e-01 1.00000000e+00</V>
<V ID="24">1.00000000e+00 1.50000000e+00 0.00000000e+00</V>
<V ID="25">1.00000000e+00 1.50000000e+00 5.00000000e-01</V>
<V ID="26">1.00000000e+00 1.50000000e+00 1.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 0 3 </E>
<E ID="2"> 0 4 </E>
<E ID="3"> 1 2 </E>
<E ID="4"> 1 5 </E>
<E ID="5"> 2 3 </E>
<E ID="6"> 2 6 </E>
<E ID="7"> 3 7 </E>
<E ID="8"> 4 5 </E>
<E ID="9"> 4 7 </E>
<E ID="10"> 5 6 </E>
<E ID="11"> 6 7 </E>
<E ID="12"> 4 8 </E>
<E ID="13"> 5 9 </E>
<E ID="14"> 6 10 </E>
<E ID="15"> 7 11 </E>
<E ID="16"> 8 9 </E>
<E ID="17"> 8 11 </E>
<E ID="18"> 9 10 </E>
<E ID="19"> 10 11 </E>
<E ID="20"> 3 13 </E>
<E ID="21"> 2 12 </E>
<E ID="22"> 12 13 </E>
<E ID="23"> 12 14 </E>
<E ID="24"> 13 15 </E>
<E ID="25"> 7 15 </E>
<E ID="26"> 6 14 </E>
<E ID="27"> 14 15 </E>
<E ID="28"> 14 16 </E>
<E ID="29"> 15 17 </E>
<E ID="30"> 11 17 </E>
<E ID="31"> 10 16 </E>
<E ID="32"> 16 17 </E>
<E ID="33"> 1 18 </E>
<E ID="34"> 18 19 </E>
<E ID="35"> 18 20 </E>
<E ID="36"> 19 2 </E>
<E ID="37"> 19 21 </E>
<E ID="38"> 5 20 </E>
<E ID="39"> 20 21 </E>
<E ID="40"> 21 6 </E>
<E ID="41"> 20 22 </E>
<E ID="42"> 21 23 </E>
<E ID="43"> 9 22 </E>
<E ID="44"> 22 23 </E>
<E ID="45"> 23 10 </E>
<E ID="46"> 19 24 </E>
<E ID="47"> 24 12 </E>
<E ID="48"> 24 25 </E>
<E ID="49"> 21 25 </E>
<E ID="50"> 25 14 </E>
<E ID="51"> 25 26 </E>
<E ID="52"> 23 26 </E>
<E ID="53"> 26 16 </E>
</EDGE>
<FACE>
<Q ID="0"> 0 3 5 1</Q>
<Q ID="1"> 0 4 8 2</Q>
<Q ID="2"> 3 6 10 4</Q>
<Q ID="3"> 5 7 11 6</Q>
<Q ID="4"> 1 2 9 7</Q>
<Q ID="5"> 8 10 11 9</Q>
<Q ID="6"> 8 13 16 12</Q>
<Q ID="7"> 10 14 18 13</Q>
<Q ID="8"> 11 15 19 14</Q>
<Q ID="9"> 9 12 17 15</Q>
<Q ID="10"> 16 18 19 17</Q>
<Q ID="11"> 5 21 22 20</Q>
<Q ID="12"> 21 23 26 6</Q>
<Q ID="13"> 22 24 27 23</Q>
<Q ID="14"> 20 7 25 24</Q>
<Q ID="15"> 11 26 27 25</Q>
<Q ID="16"> 26 28 31 14</Q>
<Q ID="17"> 27 29 32 28</Q>
<Q ID="18"> 25 15 30 29</Q>
<Q ID="19"> 19 31 32 30</Q>
<Q ID="20"> 33 34 36 3</Q>
<Q ID="21"> 33 35 38 4</Q>
<Q ID="22"> 34 37 39 35</Q>
<Q ID="23"> 36 6 40 37</Q>
<Q ID="24"> 38 39 40 10</Q>
<Q ID="25"> 38 41 43 13</Q>
<Q ID="26"> 39 42 44 41</Q>
<Q ID="27"> 40 14 45 42</Q>
<Q ID="28"> 43 44 45 18</Q>
<Q ID="29"> 36 46 47 21</Q>
<Q ID="30"> 46 48 49 37</Q>
<Q ID="31"> 47 23 50 48</Q>
<Q ID="32"> 40 49 50 26</Q>
<Q ID="33"> 49 51 52 42</Q>
<Q ID="34"> 50 28 53 51</Q>
<Q ID="35"> 45 52 53 31</Q>
</FACE>
<ELEMENT>
<H ID="0"> 0 1 2 3 4 5 </H>
<H ID="1"> 5 6 7 8 9 10 </H>
<H ID="2"> 11 3 12 13 14 15 </H>
<H ID="3"> 15 8 16 17 18 19 </H>
<H ID="4"> 20 21 22 23 2 24 </H>
<H ID="5"> 24 25 26 27 7 28 </H>
<H ID="6"> 29 23 30 31 12 32 </H>
<H ID="7"> 32 27 33 34 16 35 </H>
</ELEMENT>
<COMPOSITE>
<C ID="0"> H[0-7] </C>
<C ID="1"> F[0,11,20,29] </C>
<C ID="2"> F[1,6,21,25] </C>
<C ID="3"> F[22,26,30,33] </C>
<C ID="4"> F[13,17,31,34] </C>
<C ID="5"> F[4,9,14,18] </C>
<C ID="6"> F[10,19,28,35] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
<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="Mixed_CG_Discontinuous" />
<I PROPERTY="SubSteppingScheme" VALUE="True" />
<I PROPERTY="TimeIntegrationMethod" VALUE="BackwardEuler" />
</SOLVERINFO>
<PARAMETERS>
<P> TimeStep = 0.005 </P>
<P> NumSteps = 10 </P>
<P> IO_CheckSteps = 10 </P>
<P> IO_InfoSteps = 1 </P>
<P> Kinvis = 0.025 </P>
<P> MinSubSteps = 4 </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>
<B ID="1"> C[6] </B>
<B ID="2"> C[2] </B>
<B ID="3"> C[3] </B>
<B ID="4"> C[4] </B>
<B ID="5"> C[5] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<D VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<D VAR="w" VALUE="0" />
<D VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</REGION>
<REGION REF="1">
<D VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<D VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<D VAR="w" VALUE="0" />
<D VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</REGION>
<REGION REF="2">
<D VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<D VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<D VAR="w" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="3">
<D VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<D VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<D VAR="w" VALUE="0" />
<D VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</REGION>
<REGION REF="4">
<D VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<D VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<D VAR="w" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="5">
<D VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<D VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<D VAR="w" VALUE="0" />
<D VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<E VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<E VAR="w" VALUE="0" />
<E VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</FUNCTION>
<FUNCTION NAME="InitialConditions">
<E VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*PI*y))" />
<E VAR="v" VALUE="(-0.963740544195769/(2*PI))*exp(-0.963740544195769*x)*sin(2*PI*y)" />
<E VAR="w" VALUE="0" />
<E VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</FUNCTION>
</CONDITIONS>
</NEKTAR>
......@@ -85,14 +85,6 @@ L 2 error (variable v) : 2.91568e-14
L inf error (variable v) : 2.14091e-13
----------------------------------------
ChanStability.xml
L 2 error (variable u) : 3.63474e-11
L inf error (variable u) : 6.53318e-10
L 2 error (variable v) : 1.03201e-10
L inf error (variable v) : 1.55026e-09
L 2 error (variable w) : 7.74576e-09
L inf error (variable w) : 2.59194e-08
L 2 error (variable p) : 8.18866e-08
L inf error (variable p) : 3.42898e-07
L 2 error (variable u) : 0.0372818
L inf error (variable u) : 0
L 2 error (variable v) : 0.0228428
......@@ -169,14 +161,14 @@ L 2 error (variable v) : 2.96601e-14
L inf error (variable v) : 1.25209e-13
----------------------------------------
Test_ChanFlow_3DH1D_MVM.xml
L 2 error (variable u) : 4.67619e-16
L inf error (variable u) : 2.498e-15
L 2 error (variable v) : 9.18537e-17
L inf error (variable v) : 4.53618e-16
L 2 error (variable u) : 3.34833e-16
L inf error (variable u) : 1.94289e-15
L 2 error (variable v) : 1.21337e-16
L inf error (variable v) : 4.75925e-16
L 2 error (variable w) : 0
L inf error (variable w) : 1.53665e-16
L 2 error (variable p) : 4.76821e-14
L inf error (variable p) : 2.11386e-13
L inf error (variable w) : 1.02571e-17
L 2 error (variable p) : 2.64454e-14
L inf error (variable p) : 1.4011e-13
----------------------------------------
Test_ChanFlow_3DH1D_FFT.xml
L 2 error (variable u) : 3.61998e-16
......@@ -382,3 +374,13 @@ L inf error (variable v) : 0.00586633
L 2 error (variable p) : 0.0123206
L inf error (variable p) : 0.0501207
----------------------------------------
Test_Hex_Kovasnay_SubStep.xml
L 2 error (variable u) : 0.048372
L inf error (variable u) : 0.0616737
L 2 error (variable v) : 0.00253099
L inf error (variable v) : 0.00726816
L 2 error (variable w) : 0.000796309
L inf error (variable w) : 0.00562612
L 2 error (variable p) : 0.00900028
L inf error (variable p) : 0.0453188
----------------------------------------
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