Add new test for variable P

solvers/IncNavierStokesSolver/CMakeLists.txt

@@ -57,6 +57,7 @@ IF( NEKTAR_SOLVER_INCNAVIERSTOKES )
     ADD_NEKTAR_TEST(KovaFlow_m3)
     ADD_NEKTAR_TEST(KovaFlow_m8)
     ADD_NEKTAR_TEST(KovaFlow_m8_short_HOBC)
+    ADD_NEKTAR_TEST(KovaFlow_varP)
     #ADD_NEKTAR_TEST(KovaFlow_Oseen_m8)
     ADD_NEKTAR_TEST_LENGTHY(KovaFlow_3DH1D_P5_20modes_MVM)
     ADD_NEKTAR_TEST_LENGTHY(KovaFlow_3DH1D_P5_20modes_MVM_Deal)

solvers/IncNavierStokesSolver/Tests/KovaFlow_varP.tst

+<?xml version="1.0" encoding="utf-8"?>
+<test>
+    <description>Kovasznay Flow variable P</description>
+    <executable>IncNavierStokesSolver</executable>
+    <parameters>KovaFlow_varP.xml</parameters>
+    <files>
+        <file description="Session File">KovaFlow_varP.xml</file>
+    </files>
+    <metrics>
+        <metric type="L2" id="1">
+            <value variable="u" tolerance="1e-12">8.95241e-05</value>
+            <value variable="v" tolerance="1e-12">6.41405e-06</value>
+	    <value variable="p" tolerance="1e-12">0.000932063</value>
+        </metric>
+        <metric type="Linf" id="2">
+            <value variable="u" tolerance="1e-12">0.000144453</value>
+            <value variable="v" tolerance="1e-12">9.04985e-06</value>
+	    <value variable="p" tolerance="1e-12">0.00283058</value>
+        </metric>
+    </metrics>
+</test>

solvers/IncNavierStokesSolver/Tests/KovaFlow_varP.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">
+
+    <EXPANSIONS>
+        <E COMPOSITE="C[7]" NUMMODES="8" FIELDS="u,v,p" TYPE="MODIFIED" />
+        <E COMPOSITE="C[8]" NUMMODES="9" FIELDS="u,v,p" TYPE="MODIFIED" />
+        <E COMPOSITE="C[9]" NUMMODES="7" FIELDS="u,v,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      = 1       </P>
+            <P> IO_CheckSteps = 1       </P>
+            <P> IO_InfoSteps  = 1       </P>
+            <P> Kinvis        = 0.025        </P>
+        </PARAMETERS>
+
+        <VARIABLES>
+            <V ID="0"> u </V>
+            <V ID="1"> v </V>
+            <V ID="2"> p </V>
+        </VARIABLES>
+
+        <BOUNDARYREGIONS>
+            <B ID="0"> C[1] </B>
+            <B ID="1"> C[2] </B>
+            <B ID="2"> C[3] </B>
+        </BOUNDARYREGIONS>
+
+        <BOUNDARYCONDITIONS>
+            <REGION REF="0">
+                <D VAR="u" VALUE="1-1.619099729265964*cos(2*PI*y)" />
+                <D VAR="v" VALUE="-0.248344108585656*sin(2*PI*y)" />
+                <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
+            </REGION>
+            <REGION REF="1">
+                <D VAR="u" VALUE="1-0.381463333531742*cos(2*PI*y)" />
+                <D VAR="v" VALUE="-0.058510399212408*sin(2*PI*y)" />
+                <D VAR="p" VALUE="0.427242862585425" />
+            </REGION>
+            <REGION REF="2">
+                <N VAR="u" VALUE="0" />
+                <D VAR="v" VALUE="0" />
+                <N VAR="p" VALUE="0" />
+            </REGION>
+        </BOUNDARYCONDITIONS>
+
+        <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="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
+        </FUNCTION>
+
+        <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="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
+        </FUNCTION>
+
+    </CONDITIONS>
+
+    <GEOMETRY DIM="2" SPACE="2">
+        <VERTEX COMPRESSED="B64Z-LittleEndian" BITSIZE="64">eJxjYEAGD/aj0gwMjFjlf9jDRJhQ5D/Yo8szY5VHmM+C1fwPcP2s2N0Hl2dDk18/MZnnf1aFLUyEHUpHO9a8Cq6u2YtuPweqfgz3caK5H2r+XpgIF5o8uvu4scoj/MeD1X5E+PGiut8WXZ4Pq36E+fxY5RHuE0CTdze//3OXWwLcf4JQ+vG+paKKtTm26OYLQemmsp+NV/7326CbL4zi/py9UPPh8QMANMBNrQAA</VERTEX>
+        <EDGE COMPRESSED="B64Z-LittleEndian" BITSIZE="64">eJx1kskSgjAQBUFckDUIigv6/5/pgelLVyWXV/3yMjWTpCj2dQjtQkvxoBzchlZickflx9CTuAk9i6l3UT6F1uJrRhupz7di5unE1Ou1Dw/iSX1PyiX5twyPYnTOKHUX+dzjXczcD/XPP1jFvOtTTN8v+bz3W8wcH/n8g03MHF/5rJ+Y+f7dGwRp</EDGE>
+        <ELEMENT>
+            <Q COMPRESSED="B64Z-LittleEndian" BITSIZE="64">eJx1zskOgyAABmFsRaq0bl3t4vs/ppfh4CTl8iVDAn8I+1PhAY9/eo0RG92XnvCErd4tvcOMZ/1zwR4H7YrqI056p1Gf8aq9Sf2Gd+0t+x/4xJf2tuoLvrW3U//gV3uz+g9X7d0A9rcC/gAA</Q>
+        </ELEMENT>
+        <COMPOSITE>
+            <C ID="7"> Q[0-3]             </C>     <!-- Domain -->
+            <C ID="8"> Q[4-7]             </C>     <!-- Domain -->
+            <C ID="9"> Q[8-11]             </C>     <!-- Domain -->
+            <C ID="1"> E[23,25,27,29] </C>
+            <C ID="2"> E[3,6,9,12] </C>
+            <C ID="3"> E[0,11,13,21-22,30] </C>
+        </COMPOSITE>
+        <DOMAIN> C[7-9] </DOMAIN>
+    </GEOMETRY>
+</NEKTAR>