Closure model blocks

Closure Model lists

The closure model block in the XML input file may contain multiple closure model factories. Each closure model factory typically corresponds to a specific physical process or model solved within Vertex-CFD. Every factory is defined as a sublist of the closure model block and is uniquely identified by a Model ID. The type of closure model factory is specified by the Closure Model Factory parameter, which defaults to Navier-Stokes. Consequently, this parameter may be omitted for closure model factories associated with the Navier–Stokes equations.

The closure model block enumerates all optional closure models required to fully define the selected physics. The closure models currently available in Vertex-CFD are described in the sections below, organized by physical module, and representative examples are provided to illustrate the required XML syntax.

Conduction equation

The sublist name model_id must match the value specified under Model ID in the Physics Blocks sublist. For example, if Model ID is defined as solid in the Physics Blocks, the corresponding sublist in Closure Models must also be named solid.

 1  <ParameterList name="Physics Blocks">
 2    <ParameterList name="ConductionPhysicsBlock">
 3      <ParameterList>
 4        <Parameter name="Type" type="string" value="Conduction"/>
 5        <Parameter name="Basis Order" type="int" value="1"/>
 6        <Parameter name="Integration Order" type="int" value="2"/>
 7        <Parameter name="Model ID" type="string" value="solid"/>
 8        <Parameter name="Build Source Residual" type="bool" value="true"/>
 9      </ParameterList> <!-- -->
10    </ParameterList> <!-- ConductionPhysicsBlock -->
11  </ParameterList> <!-- Physics Blocks -->

 1      </ParameterList> <!-- Temperature -->
 2    </ParameterList> <!-- eblock-0_0 -->
 3  </ParameterList> <!-- Initial Conditions -->
 4
 5  <ParameterList name="Closure Models">
 6    <ParameterList name="solid">
 7      <Parameter name="Closure Factory Type" type="string" value="Conduction"/>
 8      <ParameterList name="Error Norm">
 9        <Parameter name="Type" type="string" value="ConductionErrorNorms"/>
10      </ParameterList> <!-- Error Norm -->
11      <ParameterList name="Exact Solution">
12        <Parameter name="Type" type="string" value="ConductionExactSolution"/>
13        <Parameter name="Volumetric Heat Source Value" type="double" value="10.0"/>
14        <Parameter name="Thermal Conductivity Coefficient" type="double" value="1.0"/>
15        <Parameter name="Right Temperature Boundary Value" type="double" value="330.0"/>
16      </ParameterList> <!-- Exact Solution -->
17      <ParameterList name="Constant heat source">
18        <Parameter name="Type" type="string" value="ConductionVolumetricSource"/>
19        <Parameter name="Heat Source Type" type="string" value="xlinear"/>
20        <Parameter name="X-value of the left boundary" type="double" value="0.0"/>
21        <Parameter name="X-value of the right boundary" type="double" value="1.0"/>
22        <Parameter name="Volumetric Heat Source Value" type="double" value="10.0"/>
23      </ParameterList> <!-- Constant heat source -->
24      <ParameterList name="Material Properties">
25        <Parameter name="Type" type="string" value="ConstantMaterialProperties"/>
26        <Parameter name="Density Value" type="double" value="1.0"/>
27        <Parameter name="Specific Heat Capacity Value" type="double" value="1.0"/>
28        <Parameter name="Thermal Conductivity Type" type="string" value="inverse_proportional"/>

Conjugate heat transfer

The sublist names model_id must match the value specified under Model ID in the Physics Blocks sublist. For example, if Model ID is defined as solid_1 or fluid in the Physics Blocks, the corresponding sublist in Closure Models must also be named fluid or solid_1.

 1  <ParameterList name="Physics Blocks">
 2    <ParameterList name="FluidPhysicsBlock">
 3      <ParameterList>
 4        <Parameter name="Type" type="string" value="IncompressibleNavierStokes"/>
 5        <Parameter name="Basis Order" type="int" value="1"/>
 6        <Parameter name="Integration Order" type="int" value="2"/>
 7        <Parameter name="Model ID" type="string" value="fluid"/>
 8        <Parameter name="Build Temperature Equation" type="bool" value="true"/>
 9        <Parameter name="Continuity Model" type="string" value="EDAC"/>
10      </ParameterList> <!-- -->
11    </ParameterList> <!-- FluidPhysicsBlock -->
12    <ParameterList name="ConductionPhysicsBlock">
13      <ParameterList>
14        <Parameter name="Type" type="string" value="Conduction"/>
15        <Parameter name="Basis Order" type="int" value="1"/>
16        <Parameter name="Integration Order" type="int" value="2"/>
17        <Parameter name="Model ID" type="string" value="solid_1"/>
18      </ParameterList> <!-- -->
19    </ParameterList> <!-- ConductionPhysicsBlock -->
20  </ParameterList> <!-- Physics Blocks -->

 1  <ParameterList name="Closure Models">
 2    <ParameterList name="fluid">
 3      <ParameterList name="Fluid Properties">
 4        <Parameter name="Kinematic viscosity" type="double" value="1.0"/>
 5        <Parameter name="Density" type="double" value="1.0"/>
 6        <Parameter name="Thermal conductivity" type="double" value="0.2"/>
 7        <Parameter name="Specific heat capacity" type="double" value="1.0"/>
 8        <Parameter name="Artificial compressibility" type="double" value="10000000.0"/>
 9        <Parameter name="Build Temperature Equation" type="bool" value="true"/>
10        <Parameter name="Continuity Model" type="string" value="EDAC"/>
11        <Parameter name="Heat capacity ratio" type="double" value="1.0"/>
12      </ParameterList> <!-- Fluid Properties -->
13      <ParameterList name="DQDT">
14        <Parameter name="Type" type="string" value="IncompressibleTimeDerivative"/>
15      </ParameterList> <!-- DQDT -->
16      <ParameterList name="CONVECTIVE_FLUX">
17        <Parameter name="Type" type="string" value="IncompressibleConvectiveFlux"/>
18      </ParameterList> <!-- CONVECTIVE_FLUX -->
19      <ParameterList name="VISCOUS_FLUX">
20        <Parameter name="Type" type="string" value="IncompressibleViscousFlux"/>
21      </ParameterList> <!-- VISCOUS_FLUX -->
22      <ParameterList name="element_length">
23        <Parameter name="Type" type="string" value="ElementLength"/>
24      </ParameterList> <!-- element_length -->
25      <ParameterList name="local_time_step_size">
26        <Parameter name="Type" type="string" value="IncompressibleLocalTimeStepSize"/>
27      </ParameterList> <!-- local_time_step_size -->
28    </ParameterList> <!-- fluid -->
29    <ParameterList name="solid_1">
30      <Parameter name="Closure Factory Type" type="string" value="Conduction"/>
31      <ParameterList name="Material Properties">
32        <Parameter name="Type" type="string" value="ConstantMaterialProperties"/>
33        <Parameter name="Density Value" type="double" value="1.0"/>
34        <Parameter name="Specific Heat Capacity Value" type="double" value="1.0"/>
35        <Parameter name="Thermal Conductivity Value" type="double" value="0.2"/>
36      </ParameterList> <!-- Material Properties -->
37    </ParameterList> <!-- solid_1 -->
38  </ParameterList> <!-- Closure Models -->