Add new tests for 3DH2D

solvers/IncNavierStokesSolver/CMakeLists.txt

@@ -44,6 +44,7 @@ IF( NEKTAR_SOLVER_INCNAVIERSTOKES )
     ADD_NEKTAR_TEST(ChanFlow_m8_BodyForce)
     ADD_NEKTAR_TEST(ChanFlow_m8_singular)
     ADD_NEKTAR_TEST(Channel_Flow_3modes_rad)    
+    ADD_NEKTAR_TEST(Couette_3DH2D_MVM)
     ADD_NEKTAR_TEST(Hex_channel_m3)
     ADD_NEKTAR_TEST(Pyr_channel_m3)
     ADD_NEKTAR_TEST_LENGTHY(Hex_channel_m8)
@@ -98,6 +99,7 @@ IF( NEKTAR_SOLVER_INCNAVIERSTOKES )
         ADD_NEKTAR_TEST(CylFlow_MovBody)
         ADD_NEKTAR_TEST(KovaFlow_3DH1D_MVM_FFTW_Consistency)
         ADD_NEKTAR_TEST(KovaFlow_m8_short_HOBC_3D1H)
+        ADD_NEKTAR_TEST(Couette_3DH2D_FFT)
     ENDIF (NEKTAR_USE_FFTW)
 
     IF (NEKTAR_USE_MPI)

solvers/IncNavierStokesSolver/Tests/Couette_3DH2D.xml

+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
+    xsi:noNamespaceSchemaLocation="http://www.nektar.info/schema/nektar.xsd">
+    <GEOMETRY DIM="1" SPACE="1">
+
+        <VERTEX>
+            <V ID="0">  0.00    0.0    0.0 </V>
+            <V ID="1">  0.20    0.0    0.0 </V>
+            <V ID="2">  0.50    0.0    0.0 </V>
+            <V ID="3">  0.80    0.0    0.0 </V>
+            <V ID="4">  1.00    0.0    0.0 </V>
+        </VERTEX>
+
+        <ELEMENT>
+            <S ID="0"> 0 1 </S>
+            <S ID="1"> 1 2 </S>
+            <S ID="2"> 2 3 </S>
+            <S ID="3"> 3 4 </S>
+        </ELEMENT>
+        <COMPOSITE>
+            <C ID="0"> S[0-3]</C>
+            <C ID="1"> V[0]</C>
+            <C ID="2"> V[4]</C>
+        </COMPOSITE>
+
+        <DOMAIN> C[0] </DOMAIN>
+
+    </GEOMETRY>
+
+    <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"/>
+          <I PROPERTY="HOMOGENEOUS" VALUE="2D"/>
+        </SOLVERINFO>
+
+        <PARAMETERS>
+            <P> TimeStep      = 0.1       </P>
+            <P> NumSteps      = 100         </P>
+            <P> IO_CheckSteps = NumSteps         </P>
+            <P> IO_InfoSteps  = NumSteps          </P>
+            <P> Kinvis        = 1         </P>
+            <P> HomModesY = 4               </P>
+            <P> LY        = 0.5             </P>
+            <P> HomModesZ = 4               </P>
+            <P> LZ        = 0.5             </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[2] </B>
+        </BOUNDARYREGIONS>
+
+        <BOUNDARYCONDITIONS>
+            <REGION REF="0">
+              <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="1">
+               <D VAR="u" VALUE="0" />
+               <D VAR="v" VALUE="1" />
+               <D VAR="w" VALUE="0.5" />
+               <D VAR="p" VALUE="0"/>
+            </REGION>
+        </BOUNDARYCONDITIONS>
+
+        <FUNCTION NAME="InitialConditions">
+            <E VAR="u" VALUE="0" />
+            <E VAR="v" VALUE="0" />
+            <E VAR="w" VALUE="0" />
+            <E VAR="p" VALUE="0" />
+        </FUNCTION>
+
+        <FUNCTION NAME="ExactSolution">
+            <E VAR="u" VALUE="0" />
+            <E VAR="v" VALUE="x" />
+            <E VAR="w" VALUE="0.5*x" />
+            <E VAR="p" VALUE="0" />
+        </FUNCTION>
+
+    </CONDITIONS>
+
+</NEKTAR>

solvers/IncNavierStokesSolver/Tests/Couette_3DH2D_FFT.tst

+<?xml version="1.0" encoding="utf-8"?>
+<test>
+    <description>Plane Couette Flow 3D homogeneous 2D, P=3, 4x4 Fourier modes (FFT)</description>
+    <executable>IncNavierStokesSolver</executable>
+    <parameters>-I USEFFT=FFTW Couette_3DH2D.xml</parameters>
+    <files>
+        <file description="Session File">Couette_3DH2D.xml</file>
+    </files>
+    <metrics>
+        <metric type="L2" id="1">
+            <value variable="u" tolerance="1e-6">0</value>
+            <value variable="v" tolerance="1e-6">0</value>
+            <value variable="w" tolerance="1e-6">0</value>
+            <value variable="p" tolerance="1e-6">0</value>
+        </metric>
+        <metric type="Linf" id="2">
+            <value variable="u" tolerance="1e-6">0</value>
+            <value variable="v" tolerance="1e-6">0</value>
+            <value variable="w" tolerance="1e-6">0</value>
+            <value variable="p" tolerance="1e-6">0</value>
+        </metric>
+    </metrics>
+</test>
+
+

solvers/IncNavierStokesSolver/Tests/Couette_3DH2D_MVM.tst

+<?xml version="1.0" encoding="utf-8"?>
+<test>
+    <description>Plane Couette Flow 3D homogeneous 2D, P=3, 4x4 Fourier modes (MVM)</description>
+    <executable>IncNavierStokesSolver</executable>
+    <parameters>Couette_3DH2D.xml</parameters>
+    <files>
+        <file description="Session File">Couette_3DH2D.xml</file>
+    </files>
+    <metrics>
+        <metric type="L2" id="1">
+            <value variable="u" tolerance="1e-6">0</value>
+            <value variable="v" tolerance="1e-6">0</value>
+            <value variable="w" tolerance="1e-6">0</value>
+            <value variable="p" tolerance="1e-6">0</value>
+        </metric>
+        <metric type="Linf" id="2">
+            <value variable="u" tolerance="1e-6">0</value>
+            <value variable="v" tolerance="1e-6">0</value>
+            <value variable="w" tolerance="1e-6">0</value>
+            <value variable="p" tolerance="1e-6">0</value>
+        </metric>
+    </metrics>
+</test>
+
+