<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="Case-2-MiddleStove" solver_version="5.0;6.0;5.0;6.0;5.2;6.0" study="Verify-VolumeHeatSource" version="2.0">
  <additional_scalars>
    <users/>
  </additional_scalars>
  <analysis_control>
    <error_estimator/>
    <output>
      <listing_printing_frequency>100</listing_printing_frequency>
      <mesh id="-1" label="Fluid domain" type="cells">
        <all_variables status="on"/>
        <location>all[]</location>
        <writer id="-1"/>
      </mesh>
      <mesh id="-2" label="Boundary" type="boundary_faces">
        <all_variables status="on"/>
        <location>all[]</location>
        <writer id="-1"/>
      </mesh>
      <probe_format choice="DAT"/>
      <probe_recording_frequency>1</probe_recording_frequency>
      <writer id="-1" label="results">
        <directory name="postprocessing"/>
        <format name="ensight" options="text"/>
        <frequency period="time_step">20</frequency>
        <output_at_end status="on"/>
        <output_at_start status="on"/>
        <time_dependency choice="fixed_mesh"/>
      </writer>
    </output>
    <profiles/>
    <scalar_balances/>
    <time_averages/>
    <time_parameters>
      <iterations>2000</iterations>
      <property label="CourantNb" name="courant_number">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <probes_recording status="off"/>
      </property>
      <property label="FourierNb" name="fourier_number">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <probes_recording status="off"/>
      </property>
      <thermal_time_step status="off"/>
      <time_passing>0</time_passing>
      <time_step_ref>0.001</time_step_ref>
    </time_parameters>
  </analysis_control>
  <boundary_conditions>
    <boundary label="WALL" name="1" nature="wall">WALL</boundary>
    <boundary label="FREE_INLET_OUTLET" name="2" nature="free_inlet_outlet">FREE_INLET_OUTLET</boundary>
    <boundary label="OUTFLOW" name="3" nature="imposed_p_outlet">OUTFLOW</boundary>
    <boundary label="SYMMETRY" name="4" nature="symmetry">SYMMETRY</boundary>
    <free_inlet_outlet field_id="none" label="FREE_INLET_OUTLET">
      <headLoss>
        <formula>K = 0;</formula>
      </headLoss>
    </free_inlet_outlet>
    <imposed_p_outlet field_id="none" label="OUTFLOW">
      <dirichlet name="pressure">101300</dirichlet>
    </imposed_p_outlet>
    <symmetry field_id="none" label="SYMMETRY"/>
    <wall field_id="none" label="WALL">
      <scalar choice="neumann" name="temperature" type="thermal">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off">
        <dirichlet component="0" name="velocity">0</dirichlet>
        <dirichlet component="1" name="velocity">0</dirichlet>
        <dirichlet component="2" name="velocity">0</dirichlet>
      </velocity_pressure>
    </wall>
  </boundary_conditions>
  <calculation_management>
    <block_io/>
    <partitioning/>
    <run_type>mesh preprocess</run_type>
    <start_restart>
      <frozen_field status="off"/>
    </start_restart>
  </calculation_management>
  <lagrangian model="off"/>
  <numerical_parameters>
    <gradient_reconstruction choice="0"/>
    <gradient_transposed status="on"/>
    <hydrostatic_pressure status="off"/>
    <pressure_relaxation>1</pressure_relaxation>
    <velocity_pressure_algo choice="simplec">
      <piso_sweep_number>1</piso_sweep_number>
    </velocity_pressure_algo>
    <velocity_pressure_coupling status="off"/>
    <wall_pressure_extrapolation>0</wall_pressure_extrapolation>
  </numerical_parameters>
  <physical_properties>
    <fluid_properties>
      <material choice="user_material"/>
      <method choice="user_properties"/>
      <property choice="user_law" label="Density" name="density">
        <formula>T0 = 298.83;	#at which rho0 is defined

#rho = rho0 [1 - beta0*(T-T0)]; beta0 = 1/T0
density = rho0*(2.0 - temperature/T0);</formula>
        <initial_value>1.17862</initial_value>
        <probes_recording status="off"/>
      </property>
      <property choice="constant" label="LamVisc" name="molecular_viscosity">
        <initial_value>1.83e-05</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <property choice="constant" label="SpecHeat" name="specific_heat">
        <initial_value>1017.24</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <property choice="constant" label="ThermalCond" name="thermal_conductivity">
        <initial_value>0.02495</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <reference choice="user_material"/>
      <reference_pressure>101325</reference_pressure>
      <reference_temperature>293.15</reference_temperature>
    </fluid_properties>
    <gravity>
      <gravity_x>0</gravity_x>
      <gravity_y>-9.8</gravity_y>
      <gravity_z>0</gravity_z>
    </gravity>
    <notebook/>
    <omega>
      <omega_x>0</omega_x>
      <omega_y>0</omega_y>
      <omega_z>0</omega_z>
    </omega>
  </physical_properties>
  <solution_domain>
    <extrusion/>
    <faces_cutting status="off"/>
    <joining>
      <face_joining name="1">
        <fraction>0.1</fraction>
        <plane>25</plane>
        <selector>INTERFACE</selector>
        <verbosity>1</verbosity>
        <visualization>1</visualization>
      </face_joining>
    </joining>
    <mesh_smoothing status="off"/>
    <meshes_list>
      <mesh name="MSH_Bot_Solid.med"/>
      <mesh name="MSH_Contrictor.med"/>
      <mesh name="MSH_Stove_Solid.med"/>
      <mesh name="MSH_Top_Solid.med"/>
      <meshdir name="../MESH/MiddleStove"/>
    </meshes_list>
    <periodicity/>
    <thin_walls/>
    <volumic_conditions>
      <zone groundwater_law="off" head_losses="off" id="1" initialization="on" label="all_cells" mass_source_term="off" momentum_source_term="off" porosity="off" scalar_source_term="off" thermal_source_term="off">all[]</zone>
      <zone groundwater_law="off" head_losses="off" id="2" initialization="off" label="VOL_STOVE" mass_source_term="off" momentum_source_term="off" porosity="off" scalar_source_term="off" thermal_source_term="on">VOL_STOVE</zone>
    </volumic_conditions>
  </solution_domain>
  <thermophysical_models>
    <ale_method status="off">
      <formula>mesh_viscosity_1 = 1;</formula>
      <mesh_viscosity type="isotrop"/>
    </ale_method>
    <atmospheric_flows model="off"/>
    <compressible_model model="off"/>
    <conjugate_heat_transfer>
      <external_coupling/>
    </conjugate_heat_transfer>
    <fans/>
    <gas_combustion model="off" option="off"/>
    <groundwater_model model="off"/>
    <hgn_model model="off"/>
    <joule_effect model="off"/>
    <porosities/>
    <radiative_transfer model="off">
      <absorption_coefficient type="constant">0</absorption_coefficient>
      <restart status="off"/>
    </radiative_transfer>
    <reference_values>
      <length/>
    </reference_values>
    <solid_fuels model="off"/>
    <source_terms>
      <thermal_formula name="temperature" zone_id="2">#S = q_dot = Watt
S = 10000.0;
dS = 0;</thermal_formula>
    </source_terms>
    <thermal_scalar model="temperature_kelvin">
      <property label="Dimensionless Thermal flux" name="boundary_layer_nusselt" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Boundary temperature" name="boundary_temperature" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Thermal flux" name="thermal_flux" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Tplus" name="tplus" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <variable label="TempK" name="temperature" type="thermal">
        <blending_factor>1</blending_factor>
        <formula zone_id="1">temperature = 300.;</formula>
        <listing_printing status="off"/>
        <max_value>1e+12</max_value>
        <min_value>0</min_value>
        <probes_recording status="off"/>
        <rhs_reconstruction>1</rhs_reconstruction>
        <turbulent_flux_model>SGDH</turbulent_flux_model>
      </variable>
    </thermal_scalar>
    <turbomachinery model="off">
      <joining/>
    </turbomachinery>
    <turbulence model="k-omega-SST">
      <gravity_terms status="on"/>
      <initialization choice="reference_value" zone_id="1"/>
      <property label="TurbVisc" name="turbulent_viscosity"/>
      <reference_velocity>1</reference_velocity>
      <variable label="k" name="k">
        <blending_factor>0</blending_factor>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
      <variable label="omega" name="omega">
        <blending_factor>0</blending_factor>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
      <wall_function>3</wall_function>
    </turbulence>
    <velocity_pressure>
      <initialization>
        <formula zone_id="1">velocity[0] = 0.;
velocity[1] = 1.0;
velocity[2] = 0.;</formula>
      </initialization>
      <property label="Stress" name="stress" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Stress, normal" name="stress_normal" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Stress, tangential" name="stress_tangential" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="total_pressure" name="total_pressure">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <probes_recording status="off"/>
      </property>
      <property label="Yplus" name="yplus" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <variable label="Pressure" name="pressure">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <probes_recording status="off"/>
        <rhs_reconstruction>2</rhs_reconstruction>
      </variable>
      <variable dimension="3" label="Velocity" name="velocity">
        <blending_factor>1</blending_factor>
        <listing_printing status="off"/>
        <probes_recording status="off"/>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
    </velocity_pressure>
  </thermophysical_models>
</Code_Saturne_GUI>