<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="Desktop" solver_version="4.0;5.0;5.3;5.0;6.0" study="xhuang" version="2.0">
  <additional_scalars/>
  <analysis_control>
    <error_estimator/>
    <output>
      <listing_printing_frequency>2000</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 name="1" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0</probe_z>
      </probe>
      <probe name="2" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.1</probe_z>
      </probe>
      <probe name="3" status="on">
        <probe_x>0.15</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.12</probe_z>
      </probe>
      <probe name="10D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.17</probe_z>
      </probe>
      <probe name="5" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>-0.15</probe_z>
      </probe>
      <probe name="6" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.4</probe_z>
      </probe>
      <probe name="7" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.14</probe_z>
      </probe>
      <probe name="8" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.14</probe_z>
      </probe>
      <probe name="1D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.017</probe_z>
      </probe>
      <probe name="2D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.034</probe_z>
      </probe>
      <probe name="3D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.051</probe_z>
      </probe>
      <probe name="4D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.068</probe_z>
      </probe>
      <probe name="5D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.085</probe_z>
      </probe>
      <probe name="6D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.102</probe_z>
      </probe>
      <probe name="7D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.119</probe_z>
      </probe>
      <probe name="8D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.136</probe_z>
      </probe>
      <probe name="9D" status="on">
        <probe_x>0</probe_x>
        <probe_y>0</probe_y>
        <probe_z>0.153</probe_z>
      </probe>
      <probe_format choice="DAT"/>
      <probe_recording_frequency>1</probe_recording_frequency>
      <writer id="-1" label="results">
        <directory name="postprocessing"/>
        <format name="ensight" options="binary"/>
        <frequency period="time_step">2000</frequency>
        <output_at_end status="on"/>
        <time_dependency choice="fixed_mesh"/>
      </writer>
    </output>
    <profiles/>
    <scalar_balances/>
    <time_averages/>
    <time_parameters>
      <iterations>1500000</iterations>
      <property label="CourantNb" name="courant_number"/>
      <property label="FourierNb" name="fourier_number"/>
      <thermal_time_step status="off"/>
      <time_passing>0</time_passing>
      <time_step_ref>1e-05</time_step_ref>
    </time_parameters>
  </analysis_control>
  <boundary_conditions>
    <boundary label="top" name="1" nature="wall">top</boundary>
    <boundary label="bottom" name="2" nature="wall">bottom</boundary>
    <boundary label="heater" name="3" nature="wall">heater</boundary>
    <boundary label="inlet" name="4" nature="inlet">inlet</boundary>
    <boundary label="inlet_side" name="5" nature="wall">inlet_side</boundary>
    <boundary label="v5" name="6" nature="wall">v5</boundary>
    <boundary label="plate_top" name="7" nature="wall">plate_top</boundary>
    <boundary label="plate_side" name="8" nature="wall">plate_side</boundary>
    <boundary label="plate_bottom" name="9" nature="wall">plate_bottom</boundary>
    <boundary label="outlet" name="10" nature="outlet">outlet</boundary>
    <boundary label="outlet_side" name="11" nature="wall">outlet_side</boundary>
    <inlet field_id="none" label="inlet">
      <mapped_inlet status="on">
        <translation_x>0</translation_x>
        <translation_y>0</translation_y>
        <translation_z>0.085</translation_z>
      </mapped_inlet>
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>508.15</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>0.017</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="flow1" direction="normal">
        <flow1>1.3</flow1>
      </velocity_pressure>
    </inlet>
    <outlet field_id="none" label="outlet">
      <scalar choice="neumann" name="temperature" type="thermal">
        <neumann>0</neumann>
      </scalar>
    </outlet>
    <variable/>
    <wall field_id="none" label="top">
      <radiative_data choice="itpimp">
        <emissivity>0.65</emissivity>
        <external_temperature_profile>300</external_temperature_profile>
        <flux>0</flux>
        <internal_temperature_profile>433</internal_temperature_profile>
        <thickness>0.1</thickness>
        <wall_thermal_conductivity>3</wall_thermal_conductivity>
      </radiative_data>
      <scalar choice="neumann" label="temperature" name="temperature" type="thermal">
        <neumann>933.04</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>
    <wall field_id="none" label="bottom">
      <radiative_data choice="itpimp">
        <emissivity>1</emissivity>
        <external_temperature_profile>300</external_temperature_profile>
        <flux>0</flux>
        <internal_temperature_profile>823</internal_temperature_profile>
        <thickness>0.1</thickness>
        <wall_thermal_conductivity>3</wall_thermal_conductivity>
      </radiative_data>
      <scalar choice="neumann" label="temperature" name="temperature" type="thermal">
        <neumann>933.04</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>
    <wall field_id="none" label="heater">
      <radiative_data choice="itpimp">
        <emissivity>0.45</emissivity>
        <external_temperature_profile>300</external_temperature_profile>
        <flux>0</flux>
        <internal_temperature_profile>595</internal_temperature_profile>
        <thickness>0.1</thickness>
        <wall_thermal_conductivity>3</wall_thermal_conductivity>
      </radiative_data>
      <scalar choice="neumann" label="temperature" name="temperature" type="thermal">
        <neumann>-35654.31</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>
    <wall field_id="none" label="inlet_side">
      <radiative_data choice="itpimp">
        <emissivity>1</emissivity>
        <external_temperature_profile>300</external_temperature_profile>
        <flux>0</flux>
        <internal_temperature_profile>823</internal_temperature_profile>
        <thickness>0.1</thickness>
        <wall_thermal_conductivity>3</wall_thermal_conductivity>
      </radiative_data>
      <scalar choice="neumann" label="temperature" 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>
    <wall field_id="none" label="v5">
      <scalar choice="neumann" name="temperature" type="thermal">
        <neumann>933.04</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="plate_top">
      <scalar choice="neumann" name="temperature" type="thermal">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="plate_side">
      <scalar choice="neumann" name="temperature" type="thermal">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="plate_bottom">
      <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>
    <wall field_id="none" label="outlet_side">
      <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>standard</run_type>
    <start_restart>
      <frozen_field status="off"/>
      <restart path="RESU/ncouple52/checkpoint"/>
      <restart_rescue>0</restart_rescue>
      <restart_with_auxiliary status="on"/>
    </start_restart>
  </calculation_management>
  <lagrangian model="off">
    <output/>
  </lagrangian>
  <numerical_parameters>
    <gradient_reconstruction choice="0"/>
    <gradient_transposed status="on"/>
    <hydrostatic_pressure status="off"/>
    <pressure_relaxation>0.9</pressure_relaxation>
    <time_scheme_order>2</time_scheme_order>
    <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>density = 11065.-1.293*temperature;</formula>
        <initial_value>1.17862</initial_value>
      </property>
      <property choice="user_law" label="LamVisc" name="molecular_viscosity">
        <formula>molecular_viscosity = 4.94/10^4*exp(754.1/temperature);</formula>
        <initial_value>1.83e-05</initial_value>
      </property>
      <property choice="user_law" label="SpecHeat" name="specific_heat">
        <formula>specific_heat = 164.8-3.94/10^2*temperature+1.25/10^5*temperature^2-4.56*10^5/temperature^2;</formula>
        <initial_value>520</initial_value>
      </property>
      <property choice="user_law" label="ThermalCond" name="thermal_conductivity">
        <formula>thermal_conductivity = 3.284 + 1.67/10^2*temperature-2.305/10^6.*temperature^2.;</formula>
        <initial_value>0.02495</initial_value>
      </property>
      <reference choice="user_material"/>
      <reference_pressure>101325</reference_pressure>
      <reference_temperature>535.15</reference_temperature>
    </fluid_properties>
    <gravity>
      <gravity_x>0</gravity_x>
      <gravity_y>0</gravity_y>
      <gravity_z>-9.8</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/>
    <mesh_input path="../MESH/mesh_output1"/>
    <mesh_smoothing status="off"/>
    <meshes_list>
      <mesh name="cylinde_4_64m.med" path="C:/Users/xiaoxue/Desktop/cylinder-natural-5.0/MESH"/>
    </meshes_list>
    <periodicity/>
    <thin_walls/>
    <volumic_conditions>
      <zone darcy_law="off" 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="on" label="solid" mass_source_term="off" momentum_source_term="off" porosity="off" scalar_source_term="off" thermal_source_term="off">all[]</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>
      <directions_number>3</directions_number>
      <frequency>1</frequency>
      <intensity_resolution_listing_printing>0</intensity_resolution_listing_printing>
      <property label="Absorption" name="rad_absorption"/>
      <property label="Absorption_coefficient" name="rad_absorption_coeff"/>
      <property label="Coeff_ech_convectif" name="rad_exchange_coefficient" support="boundary"/>
      <property label="Emission" name="rad_emission"/>
      <property label="Emissivity" name="emissivity" support="boundary"/>
      <property label="Flux_convectif" name="rad_convective_flux" support="boundary"/>
      <property label="Flux_incident" name="rad_incident_flux" support="boundary"/>
      <property label="Flux_net" name="rad_net_flux" support="boundary"/>
      <property label="Qrad" name="radiative_flux"/>
      <property label="Radiative_source_term" name="radiative_source_term"/>
      <property label="Thickness" name="thickness" support="boundary"/>
      <property label="Thermal_conductivity" name="wall_thermal_conductivity" support="boundary"/>
      <quadrature>1</quadrature>
      <restart status="off"/>
      <temperature_listing_printing>1</temperature_listing_printing>
      <thermal_radiative_source_term>2</thermal_radiative_source_term>
    </radiative_transfer>
    <reference_values>
      <length/>
    </reference_values>
    <solid_fuels model="off"/>
    <source_terms/>
    <thermal_scalar model="temperature_kelvin">
      <property label="Dimensionless Thermal flux" name="boundary_layer_nusselt" support="boundary"/>
      <property label="Boundary temperature" name="boundary_temperature" support="boundary"/>
      <property label="Thermal flux" name="thermal_flux" support="boundary"/>
      <property label="Tplus" name="tplus" support="boundary"/>
      <variable label="TempK" name="temperature" type="thermal">
        <blending_factor>1</blending_factor>
        <formula zone_id="1">temperature = 523;</formula>
        <formula zone_id="2">temperature = 523;</formula>
        <max_value>1e+12</max_value>
        <min_value>0</min_value>
        <rhs_reconstruction>10</rhs_reconstruction>
        <time_step_factor>50</time_step_factor>
        <turbulent_flux_model>SGDH</turbulent_flux_model>
      </variable>
    </thermal_scalar>
    <turbomachinery model="off">
      <joining/>
    </turbomachinery>
    <turbulence model="LES_dynamique">
      <gravity_terms status="off"/>
      <initialization choice="reference_value" zone_id="1">
        <formula>cmu = 0.09;
k = 1.5*(0.02*uref)^2;
epsilon = k^1.5*cmu/almax;</formula>
      </initialization>
      <property label="TurbVisc" name="turbulent_viscosity"/>
      <reference_length choice="prescribed">0.017</reference_length>
      <reference_velocity>0.62</reference_velocity>
      <wall_function>3</wall_function>
    </turbulence>
    <velocity_pressure>
      <initialization>
        <formula zone_id="1">velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;</formula>
        <formula zone_id="2">velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;</formula>
      </initialization>
      <property label="Stress" name="stress" support="boundary"/>
      <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"/>
      <property label="Yplus" name="yplus" support="boundary"/>
      <variable label="Pressure" name="pressure">
        <rhs_reconstruction>10</rhs_reconstruction>
      </variable>
      <variable dimension="3" label="Velocity" name="velocity">
        <blending_factor>1</blending_factor>
        <rhs_reconstruction>10</rhs_reconstruction>
      </variable>
    </velocity_pressure>
  </thermophysical_models>
</Code_Saturne_GUI>