Commit 373d0ec0 authored by Chris Cantwell's avatar Chris Cantwell

Ported more IncNavierStokesSolver tests.

Minor updates to Tester.
parent 601ed819
......@@ -17,6 +17,10 @@ ADD_SOLVER_EXECUTABLE(IncNavierStokesSolver solvers
${IncNavierStokesSolverSource})
ADD_NEKTAR_TEST(Channel_Flow_3modes)
ADD_NEKTAR_TEST(Channel_Flow_8modes_Body_Force)
ADD_NEKTAR_TEST(Channel_Flow_8modes)
ADD_NEKTAR_TEST(Channel_Flow_2D_BC_from_file)
ADD_NEKTAR_TEST(Kovasznay_Flow_3modes)
IF (NEKTAR_USE_MPI)
ADD_NEKTAR_TEST(Channel_Flow_3modes_Parallel)
......
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>Channel Flow boundary conditions from file</description>
<executable>IncNavierStokesSolver</executable>
<parameters>Channel_Flow_2D_BC_from_file.xml</parameters>
<files>
<file description="Session File">Channel_Flow_2D_BC_from_file.xml</file>
<file description="Session File">Channel_Flow_2D_BC_from_file_u_1.bc</file>
<file description="Session File">Channel_Flow_2D_BC_from_file_u_3.bc</file>
<file description="Session File">Channel_Flow_2D_BC_from_file_v_3.bc</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12">1.85224e-13</value>
<value variable="v" tolerance="1e-12">2.65241e-13</value>
<value variable="p" tolerance="1e-8">4.8272e-11</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12">8.31335e-13</value>
<value variable="v" tolerance="1e-12">2.54448e-13</value>
<value variable="p" tolerance="1e-8">1.03323e-10</value>
</metric>
</metrics>
</test>
<?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[0]" NUMMODES="6" FIELDS="u,v,p" TYPE="MODIFIED" />
</EXPANSIONS>
<CONDITIONS>
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="AdvectionForm" VALUE="Convective" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder1" />
<I PROPERTY="solvertype" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="Projection" VALUE="Galerkin" />
</SOLVERINFO>
<PARAMETERS>
<P> TimeStep = 0.001 </P>
<P> NumSteps = 1000 </P>
<P> IO_CheckSteps = 10000 </P>
<P> IO_InfoSteps = 10 </P>
<P> Kinvis = 1 </P>
<P> gamma =2 </P>
<P> C =20 </P>
<P> p =3.14 </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>
<B ID="3"> C[4] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="u" VALUE="0" />
<D VAR="v" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="1">
<D VAR="u" FILE="Channel_Flow_2D_BC_from_file_u_1.bc" />
<!-- y*(10-y) -->
<D VAR="v" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="2">
<N VAR="u" VALUE="0" />
<N VAR="v" VALUE="0" />
<D VAR="p" VALUE="0" />
</REGION>
<REGION REF="3">
<N VAR="u" FILE="Channel_Flow_2D_BC_from_file_u_3.bc"/>
<!-- 10-2*y (y=5.0) -->
<D VAR="v" FILE="Channel_Flow_2D_BC_from_file_v_3.bc" />
<!-- zero field -->
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="BodyForce">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
</FUNCTION>
<FUNCTION NAME="InitialConditions">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="p" VALUE="0" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="y*(10-y)" />
<E VAR="v" VALUE="0" />
<E VAR="p" VALUE="-2*Kinvis*(x-1)" />
</FUNCTION>
a
</CONDITIONS>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0"> 0.0 0.0 0.0 </V>
<V ID="1"> 0.5 0.0 0.0 </V>
<V ID="2"> 1.0 0.0 0.0 </V>
<V ID="3"> 0.0 5.0 0.0 </V>
<V ID="4"> 0.5 5.0 0.0 </V>
<V ID="5"> 1.0 5.0 0.0 </V>
<V ID="6"> 0.0 10.0 0.0 </V>
<V ID="7"> 0.5 10.0 0.0 </V>
<V ID="8"> 1.0 10.0 0.0 </V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 0 3 </E>
<E ID="3"> 1 4 </E>
<E ID="4"> 2 5 </E>
<E ID="5"> 3 4 </E>
<E ID="6"> 4 5 </E>
<E ID="7"> 3 6 </E>
<E ID="8"> 4 7 </E>
<E ID="9"> 5 8 </E>
<E ID="10"> 6 7 </E>
<E ID="11"> 7 8 </E>
</EDGE>
<ELEMENT>
<Q ID="0"> 0 3 5 2 </Q>
<Q ID="1"> 1 4 6 3 </Q>
<Q ID="2"> 5 8 10 7 </Q>
<Q ID="3"> 6 9 11 8 </Q>
</ELEMENT>
<COMPOSITE>
<C ID="0"> Q[0-3] </C>
<C ID="1"> E[0,1,10,11] </C> <!-- Walls -->
<C ID="2"> E[2,7] </C> <!-- Inflow -->
<C ID="3"> E[4,9] </C> <!-- Outflow -->
<C ID="4"> E[5] </C> <!-- Horizontal -->
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<ELEMENTS FIELDS="u" SHAPE="Segment" BASIS="Modified_A" NUMMODESPERDIR="UNIORDER:6" ID="2,7">eJxjYEAGlg4w+sa23NtqXly2//6DwPs9Bbe35e6+wLAXVR0DAyOU36/iOWkKK4stRJTHdn3u7W29idy2AHMwGGEA</ELEMENTS>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<ELEMENTS FIELDS="u" SHAPE="Segment" BASIS="Modified_A" NUMMODESPERDIR="UNIORDER:6" ID="5">eJwBMADP/zKKOwltWo08StLQHhxLpDy+zPYYx46FPKshYM7xNoQ8AMVtoqmiSTy9LREC+Il9vAD2FcwA</ELEMENTS>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<ELEMENTS FIELDS="v" SHAPE="Segment" BASIS="Modified_A" NUMMODESPERDIR="UNIORDER:6" ID="5">eJxjYCANAAAAMAAB</ELEMENTS>
</NEKTAR>
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>Channel Flow P=7</description>
<executable>IncNavierStokesSolver</executable>
<parameters>Channel_Flow_8modes.xml</parameters>
<files>
<file description="Session File">Channel_Flow_8modes.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12">4.15938e-16</value>
<value variable="v" tolerance="1e-12">1.73924e-16</value>
<value variable="p" tolerance="1e-12">9.13284e-15</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12">1.50435e-14</value>
<value variable="v" tolerance="1e-12">1.0547e-15</value>
<value variable="p" tolerance="1e-12">8.199e-14</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8"?>
<NEKTAR xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
xsi:noNamespaceSchemaLocation="http://www.nektar.info/nektar.xsd">
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="8" 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 = 1000 </P>
<P> IO_CheckSteps = 1000 </P>
<P> IO_InfoSteps = 1000 </P>
<P> Kinvis = 1 </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="0" />
<D VAR="v" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="1">
<D VAR="u" VALUE="y*(1-y)" />
<D VAR="v" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="2">
<N VAR="u" VALUE="0" />
<N VAR="v" VALUE="0" />
<D VAR="p" VALUE="0" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="InitialConditions">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="p" VALUE="0" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="y*(1-y)" />
<E VAR="v" VALUE="0" />
<E VAR="p" VALUE="-2*Kinvis*(x-1)" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0"> 0.0 0.0 0.0 </V>
<V ID="1"> 0.5 0.0 0.0 </V>
<V ID="2"> 1.0 0.0 0.0 </V>
<V ID="3"> 0.0 0.5 0.0 </V>
<V ID="4"> 0.5 0.5 0.0 </V>
<V ID="5"> 1.0 0.5 0.0 </V>
<V ID="6"> 0.0 1.0 0.0 </V>
<V ID="7"> 0.5 1.0 0.0 </V>
<V ID="8"> 1.0 1.0 0.0 </V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 0 3 </E>
<E ID="3"> 1 4 </E>
<E ID="4"> 2 5 </E>
<E ID="5"> 3 4 </E>
<E ID="6"> 4 5 </E>
<E ID="7"> 3 6 </E>
<E ID="8"> 4 7 </E>
<E ID="9"> 5 8 </E>
<E ID="10"> 6 7 </E>
<E ID="11"> 7 8 </E>
</EDGE>
<ELEMENT>
<Q ID="0"> 0 3 5 2 </Q>
<Q ID="1"> 1 4 6 3 </Q>
<Q ID="2"> 5 8 10 7 </Q>
<Q ID="3"> 6 9 11 8 </Q>
</ELEMENT>
<COMPOSITE>
<C ID="0"> Q[0-3] </C>
<C ID="1"> E[0,1,10,11] </C> <!-- Walls -->
<C ID="2"> E[2,7] </C> <!-- Inflow -->
<C ID="3"> E[4,9] </C> <!-- Outflow -->
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>Channel Flow P=7 Body Force</description>
<executable>IncNavierStokesSolver</executable>
<parameters>Channel_Flow_8modes_Body_Force.xml</parameters>
<files>
<file description="Session File">Channel_Flow_8modes_Body_Force.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12">9.90432e-06</value>
<value variable="v" tolerance="1e-12">0</value>
<value variable="p" tolerance="1e-12">5.171e-14</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12">1.40068e-05</value>
<value variable="v" tolerance="1e-12">4.82146e-16</value>
<value variable="p" tolerance="1e-12">8.575e-14</value>
</metric>
</metrics>
</test>
<?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[0]" NUMMODES="8" 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 = 1000 </P>
<P> IO_CheckSteps = 1000 </P>
<P> IO_InfoSteps = 1000 </P>
<P> Kinvis = 1 </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="0" />
<D VAR="v" VALUE="0" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="1">
<P VAR="u" VALUE="[2]" />
<P VAR="v" VALUE="[2]" />
<P VAR="p" VALUE="[2]" />
</REGION>
<REGION REF="2">
<P VAR="u" VALUE="[1]" />
<P VAR="v" VALUE="[1]" />
<P VAR="p" VALUE="[1]" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="InitialConditions">
<E VAR="u" VALUE="0" />
<E VAR="v" VALUE="0" />
<E VAR="p" VALUE="0" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="y*(1-y)" />
<E VAR="v" VALUE="0" />
<E VAR="p" VALUE="0" />
</FUNCTION>
<FUNCTION NAME="BodyForce">
<E VAR="u" VALUE="2" />
<E VAR="v" VALUE="0" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0"> 0.0 0.0 0.0 </V>
<V ID="1"> 0.5 0.0 0.0 </V>
<V ID="2"> 1.0 0.0 0.0 </V>
<V ID="3"> 0.0 0.5 0.0 </V>
<V ID="4"> 0.5 0.5 0.0 </V>
<V ID="5"> 1.0 0.5 0.0 </V>
<V ID="6"> 0.0 1.0 0.0 </V>
<V ID="7"> 0.5 1.0 0.0 </V>
<V ID="8"> 1.0 1.0 0.0 </V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 0 3 </E>
<E ID="3"> 1 4 </E>
<E ID="4"> 2 5 </E>
<E ID="5"> 3 4 </E>
<E ID="6"> 4 5 </E>
<E ID="7"> 3 6 </E>
<E ID="8"> 4 7 </E>
<E ID="9"> 5 8 </E>
<E ID="10"> 6 7 </E>
<E ID="11"> 7 8 </E>
</EDGE>
<ELEMENT>
<Q ID="0"> 0 3 5 2 </Q>
<Q ID="1"> 1 4 6 3 </Q>
<Q ID="2"> 5 8 10 7 </Q>
<Q ID="3"> 6 9 11 8 </Q>
</ELEMENT>
<COMPOSITE>
<C ID="0"> Q[0-3] </C>
<C ID="1"> E[0,1,10,11] </C> <!-- Walls -->
<C ID="2"> E[2,7] </C> <!-- Inflow -->
<C ID="3"> E[4,9] </C> <!-- Outflow -->
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<ELEMENTS FIELDS="u,v,p" SHAPE="Quadrilateral" BASIS="Modified_A,Modified_A" NUMMODESPERDIR="UNIORDER:3,3" ID="0,1,2,3,4,5,6,7,8,9,10,11">eJx1kX0w1Hkcx3+IkXDu0GWPUm1hCLk6YfUhpLNyictDhZLa6oxOVNxxHqO5iUSRHpCikZ60F5X13ZsrXJ42D2nHriublUhsKmdXbrP33Zpmvr8/v/Oa13zer1+aRx8zjvYWhtfc6aibmAAb8eOKrdbd6F2ME0unSQS2y7upN0ZCiOzxZUVE9MD76Q9fS21ScaCe1kAZGjVfHdKr20jk0/73z6e1WRyR+3XcncYuOXcjhUfIoF5bgGy6DCWMjru5pzQiEo/9lb6R9v/ShLAvOF0/naW457eh0xzG4aBLo/wypHrfI79Z3EDksb+2ZhPTX+7fO82vafJW3CP3OCeNSN8OScrQfMqq2jS4EeEuT6L8HXwWTUPhabMSljcP4Z0sy6A8d002NBiNmV2YbgfchTrXnaY/pwrleeiwaZFcIo93elCSF7Zyv138hoKBQB7CXXju+TdHllUhA8cFx20nuYjE451P/Hovz9ZiA5VVaryC6gDcxSU9/FbU1A2Uv/nGhJUbl8jjLqM0SaGB3F+pXqAHETyEu1ClK2PDTauQ73j0G18eF+EuKz2/zH7xepZLZ2X5ghaLJoR3dhp7MCraROjc93FzkzV4gLvM/P6iWZxEyQMr7UscIo93Yn/jHIuvZSualF1mPEjCsE98+juHhxCJxzux/6yEsWyLGg9wlw///3IRxzm9dfqvCH8OkcddsN/k+hpvLjQpu3zwLC7mcLZb5O5akogQRfg8B7yYnPR9DIoyvr28UAttq3w5RdfXgFmDPcE38pqRwx6JQV8xH40FS3kFJmIg8Yr3sNph7VdJRQJV0Lqg6ufTrg5qCg+Ip1bKeKV8+OHq0s1Wi8SIxFumiAZMx7c5U0luXYsfaKFYn+ETFiMacIUlmLknK7CtKKiUj4aC6NKERWIiT9obvFtxz7PmRzbMMj4YrLvwVM9MTOyDdxa/HDLw1lyCKP+lOWPaDMBdssyrauKiRSg/YayuPVkEJB7vDC3Id/zn1hdQsfDahGs9HeEuxq5Cgf0BEQhTvBuuHxUhEo93nrN9x+w1kvtrnASqlgzAXcLGulUzY0SobF+jnV+m/B4CT1qLuyQyh7W8DokAmZ2xZeSJiH3wTq76I32HzBxGVE6J/4EgN8BdKs90sdv9RtCzwlTray5CIPF4p+I9uTbDSWbsYGOPcJeLCg+8N6NNndgoRCQe77RNzi6hdFI5fdWnXAPC3AB3WWaouEdwsiRgrZcQSDxpL+5iRH80c8+ded4l67cIUfHqY9GSgB7wCelaM5bAh9y2WFF2TDvgd0ejI45pht1QLHzXbfjj049+zStHO+r8oTe80fPK2klE4vG7bktObrncX6DpcvtPuR9rkuwsYs8MeUKTdetjIXz0fM7jdyttdm+k3L+6lv6zyaZP7tFomKQ6AyAr6mYQWvfR8zmP39kPYrzD5P6F65Od33xyD2UqO9+m4QPhyHLLoOskwl1WzbbL6x+8B/ST7YfiLnYodyacf140kMsFboXNlJ5KH+AuYb0urNr4ThAyE9X6F0gRicc7520tf734+T3Ys+PqCP9yB+Auwu39+ro7O4Fe3777yBwpIvF45zfsvpG7cr9N733jnWrye/7vordwR6RpSidYl0z85GUmRSQed6k0Kj+2X3wPbtOuHTSp6gDcZVXE1ruhkZ1QOr4XCb6SItxlfnOsWuigCIVID/9dJutR7lySEd9q3y9G5gUH1a8mDSq7ZAayOc1HZeC2IdSNb6fCJfF4p9qhV/vLn4tQRH2qzv3pHmWX0N2udqa/yuDkrVMbf5mrwiXxeGe97q6uh3I/41iIeUvqoLKLoft326KzZVC3/P1DD3sVLonHXaqP6wn6BkSIVfLH2QxVgbJLkF98eHWCDDK+Pctvpalw/wPKM9E+</ELEMENTS>
</NEKTAR>
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>Kovasznay Flow P=2</description>
<executable>IncNavierStokesSolver</executable>
<parameters>Kovasznay_Flow_3modes.xml</parameters>
<files>
<file description="Session File">Kovasznay_Flow_3modes.xml</file>
<file description="Session File">Kovasznay_Flow_3modes.rst</file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-3">0.343606</value>
<value variable="v" tolerance="1e-3">0.0255658</value>
<value variable="p" tolerance="1e-3">0.347747</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-3">0.322039</value>
<value variable="v" tolerance="1e-3">0.0243043</value>
<value variable="p" tolerance="1e-3">0.519425</value>
</metric>
</metrics>
</test>
<?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[0]" NUMMODES="3" 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 = 100 </P>
<P> IO_CheckSteps = 100 </P>
<P> IO_InfoSteps = 100 </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*3.141592653589793*y)" />
<D VAR="v" VALUE="-0.248344108585656*sin(2*3.141592653589793*y)" />
<N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
</REGION>
<REGION REF="1">
<D VAR="u" VALUE="1-0.381463333531742*cos(2*3.141592653589793*y)" />
<D VAR="v" VALUE="-0.058510399212408*sin(2*3.141592653589793*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">
<F VAR="u,v,p" FILE="Kovasznay_Flow_3modes.rst" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="(1-exp(-0.963740544195769*x)*cos(2*3.141592653589793*y))" />
<E VAR="v"
VALUE="(-0.963740544195769/(2*3.141592653589793))*exp(-0.963740544195769*x)*sin(2*3.141592653589793*y)" />
<E VAR="p" VALUE="0.5*(1-exp(-2*0.963740544195769*x))" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="2" SPACE="2">
<VERTEX>
<V ID="0">-5.000e-01 -5.000e-01 0.000e+00</V>
<V ID="1">-5.000e-01 1.500e+00 0.000e+00</V>
<V ID="2">1.000e+00 1.500e+00 0.000e+00</V>
<V ID="3">1.000e+00 -5.000e-01 0.000e+00</V>
<V ID="4">-5.000e-01 1.000e+00 0.000e+00</V>
<V ID="5">-5.000e-01 5.000e-01 0.000e+00</V>
<V ID="6">-5.000e-01 1.388e-12 0.000e+00</V>
<V ID="7">-1.619e-12 -5.000e-01 0.000e+00</V>
<V ID="8">5.000e-01 -5.000e-01 0.000e+00</V>
<V ID="9">1.000e+00 -1.388e-12 0.000e+00</V>
<V ID="10">1.000e+00 5.000e-01 0.000e+00</V>
<V ID="11">1.000e+00 1.000e+00 0.000e+00</V>
<V ID="12">5.000e-01 1.500e+00 0.000e+00</V>
<V ID="13">1.619e-12 1.500e+00 0.000e+00</V>
<V ID="14">5.000e-01 1.000e+00 0.000e+00</V>
<V ID="15">5.000e-01 5.000e-01 0.000e+00</V>
<V ID="16">5.000e-01 -4.626e-13 0.000e+00</V>
<V ID="17">8.097e-13 1.000e+00 0.000e+00</V>
<V ID="18">5.551e-17 5.000e-01 0.000e+00</V>
<V ID="19">-8.095e-13 4.626e-13 0.000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 2 12 </E>
<E ID="1"> 12 14 </E>
<E ID="2"> 14 11 </E>
<E ID="3"> 11 2 </E>
<E ID="4"> 14 15 </E>
<E ID="5"> 15 10 </E>
<E ID="6"> 10 11 </E>
<E ID="7"> 15 16 </E>
<E ID="8"> 16 9 </E>
<E ID="9"> 9 10 </E>
<E ID="10"> 16 8 </E>
<E ID="11"> 8 3 </E>
<E ID="12"> 3 9 </E>
<E ID="13"> 12 13 </E>
<E ID="14"> 13 17 </E>
<E ID="15"> 17 14 </E>
<E ID="16"> 17 18 </E>
<E ID="17"> 18 15 </E>
<E ID="18"> 18 19 </E>
<E ID="19"> 19 16 </E>
<E ID="20"> 19 7 </E>
<E ID="21"> 7 8 </E>
<E ID="22"> 13 1 </E>
<E ID="23"> 1 4 </E>
<E ID="24"> 4 17 </E>
<E ID="25"> 4 5 </E>
<E ID="26"> 5 18 </E>
<E ID="27"> 5 6 </E>
<E ID="28"> 6 19 </E>
<E ID="29"> 6 0 </E>
<E ID="30"> 0 7 </E>
</EDGE>