Hex_channel_m3_Flowrate.xml 10.7 KB
Newer Older
Dave Moxey's avatar
Dave Moxey committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
    <GEOMETRY DIM="3" SPACE="3">
        <VERTEX>
            <V ID="0">0.00000000e+00 0.00000000e+00 0.00000000e+00</V>
            <V ID="1">5.00000000e-01 0.00000000e+00 0.00000000e+00</V>
            <V ID="2">5.00000000e-01 5.00000000e-01 0.00000000e+00</V>
            <V ID="3">0.00000000e+00 5.00000000e-01 0.00000000e+00</V>
            <V ID="4">0.00000000e+00 0.00000000e+00 5.00000000e-01</V>
            <V ID="5">5.00000000e-01 0.00000000e+00 5.00000000e-01</V>
            <V ID="6">5.00000000e-01 5.00000000e-01 5.00000000e-01</V>
            <V ID="7">0.00000000e+00 5.00000000e-01 5.00000000e-01</V>
            <V ID="8">0.00000000e+00 0.00000000e+00 1.00000000e+00</V>
            <V ID="9">5.00000000e-01 0.00000000e+00 1.00000000e+00</V>
            <V ID="10">5.00000000e-01 5.00000000e-01 1.00000000e+00</V>
            <V ID="11">0.00000000e+00 5.00000000e-01 1.00000000e+00</V>
            <V ID="12">5.00000000e-01 1.00000000e+00 0.00000000e+00</V>
            <V ID="13">0.00000000e+00 1.00000000e+00 0.00000000e+00</V>
            <V ID="14">5.00000000e-01 1.00000000e+00 5.00000000e-01</V>
            <V ID="15">0.00000000e+00 1.00000000e+00 5.00000000e-01</V>
            <V ID="16">5.00000000e-01 1.00000000e+00 1.00000000e+00</V>
            <V ID="17">0.00000000e+00 1.00000000e+00 1.00000000e+00</V>
            <V ID="18">1.00000000e+00 0.00000000e+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.00000000e+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.00000000e+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.00000000e+00 0.00000000e+00</V>
            <V ID="25">1.00000000e+00 1.00000000e+00 5.00000000e-01</V>
            <V ID="26">1.00000000e+00 1.00000000e+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="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" />
Dave Moxey's avatar
Dave Moxey committed
157
            <I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder3" />
Dave Moxey's avatar
Dave Moxey committed
158 159 160 161
        </SOLVERINFO>

        <PARAMETERS>
            <P> TimeStep      = 0.01 </P>
Dave Moxey's avatar
Dave Moxey committed
162
            <P> NumSteps      = 500 </P>
Dave Moxey's avatar
Dave Moxey committed
163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242
            <P> IO_CheckSteps = 0    </P>
            <P> IO_InfoSteps  = 100    </P>
            <P> Kinvis        = 1   </P>
            <P> Flowrate      = 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">
                <P VAR="u" VALUE="[1]" />
                <P VAR="v" VALUE="[1]" />
                <P VAR="w" VALUE="[1]" />
                <P VAR="p" VALUE="[1]" />
            </REGION>
            <REGION REF="1">
                <P VAR="u" VALUE="[0]" USERDEFINEDTYPE="Flowrate" />
                <P VAR="v" VALUE="[0]" />
                <P VAR="w" VALUE="[0]" />
                <P 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">
                <P VAR="u" VALUE="[5]" />
                <P VAR="v" VALUE="[5]" />
                <P VAR="w" VALUE="[5]" />
                <P VAR="p" VALUE="[5]" />
            </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">
                <P VAR="u" VALUE="[3]" />
                <P VAR="v" VALUE="[3]" />
                <P VAR="w" VALUE="[3]" />
                <P VAR="p" VALUE="[3]" />
            </REGION>
        </BOUNDARYCONDITIONS>

        <FUNCTION NAME="FlowrateForce">
            <E VAR="ForceX" VALUE="0" />
            <E VAR="ForceY" VALUE="0" />
            <E VAR="ForceZ" VALUE="1" />
        </FUNCTION>
        <FUNCTION NAME="ExactSolution">
            <E VAR="u" VALUE="0" />
            <E VAR="v" VALUE="0" />
            <E VAR="w" VALUE="6*y*(1-y)" />
            <E VAR="p" VALUE="0" />
        </FUNCTION>
        <FUNCTION NAME="InitialConditions">
            <E VAR="u" VALUE="0" />
            <E VAR="v" VALUE="0" />
            <E VAR="w" VALUE="0" />
            <E VAR="p" VALUE="0" />
        </FUNCTION>
    </CONDITIONS>
</NEKTAR>