<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="Thermal_coupling" solver_version="6.0;;7.0" study="Forum_cases" version="2.0">
  <additional_scalars>
    <users/>
  </additional_scalars>
  <analysis_control>
    <output>
      <listing_printing_frequency>1</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="CSV"/>
      <probe_recording_frequency>1</probe_recording_frequency>
      <probes_interpolation choice=""/>
      <probes_snap choice=""/>
      <writer id="-1" label="results">
        <directory name="postprocessing"/>
        <format name="med" options="separate_meshes"/>
        <frequency period="time_step">10</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>100</iterations>
      <property label="CourantNb" name="courant_number"/>
      <property label="FourierNb" name="fourier_number"/>
      <time_passing>0</time_passing>
      <time_step_ref>0.00125</time_step_ref>
    </time_parameters>
  </analysis_control>
  <boundary_conditions>
    <boundary label="Face1" name="1" nature="wall">Face1</boundary>
    <boundary label="Face2" name="2" nature="wall">Face2</boundary>
    <boundary label="FluxH" name="3" nature="wall">FluxH</boundary>
    <boundary label="Inlet" name="4" nature="inlet">Inlet</boundary>
    <boundary label="Outlet" name="5" nature="outlet">Outlet</boundary>
    <boundary label="SymFluid" name="6" nature="symmetry">SymFluid</boundary>
    <boundary label="SymSolid" name="7" nature="symmetry">SymSolid</boundary>
    <boundary label="Wall" name="8" nature="wall">Wall</boundary>
    <inlet field_id="none" label="None">
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <norm>1</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="Inlet">
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>35</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <norm>0.1</norm>
      </velocity_pressure>
    </inlet>
    <outlet field_id="none" label="Outlet">
      <scalar choice="neumann" name="temperature" type="thermal">
        <neumann>0</neumann>
      </scalar>
    </outlet>
    <symmetry field_id="none" label="SymFluid"/>
    <symmetry field_id="none" label="SymSolid"/>
    <wall field_id="none" label="Face2">
      <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="FluxH">
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>100</dirichlet>
      </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="Face1">
      <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="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>
        <roughness>0.001</roughness>
      </velocity_pressure>
    </wall>
  </boundary_conditions>
  <calculation_management>
    <block_io/>
    <partitioning/>
    <start_restart>
      <frozen_field status="off"/>
    </start_restart>
  </calculation_management>
  <lagrangian model="off"/>
  <numerical_parameters>
    <algo_density_variation choice=""/>
    <extended_neighborhood choice=""/>
    <gradient_reconstruction choice="green_iter"/>
    <gradient_transposed status="on"/>
    <pressure_relaxation>0.9</pressure_relaxation>
    <velocity_pressure_algo choice="simplec">
      <piso_sweep_number>1</piso_sweep_number>
    </velocity_pressure_algo>
    <velocity_pressure_coupling status="off"/>
  </numerical_parameters>
  <physical_properties>
    <fluid_properties>
      <material choice="user_material"/>
      <method choice="user_properties"/>
      <property choice="constant" label="Density" name="density">
        <initial_value>998.2</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <zone zone_id="3">
          <formula>density = 7850;</formula>
        </zone>
      </property>
      <property choice="constant" label="DiffDyn" name="dynamic_diffusion">
        <initial_value>0.01</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <property choice="constant" label="LamVisc" name="molecular_viscosity">
        <initial_value>0.0010016</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <zone zone_id="3">
          <formula>molecular_viscosity = mu0;</formula>
        </zone>
      </property>
      <property choice="constant" label="SpecHeat" name="specific_heat">
        <initial_value>4180</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <zone zone_id="3">
          <formula>specific_heat = 420;</formula>
        </zone>
      </property>
      <property choice="constant" label="ThermalCond" name="thermal_conductivity">
        <initial_value>0.598</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
        <zone zone_id="3">
          <formula>thermal_conductivity = 45;</formula>
        </zone>
      </property>
      <reference_pressure>101325</reference_pressure>
      <reference_temperature>35</reference_temperature>
    </fluid_properties>
    <gravity>
      <gravity_x>0</gravity_x>
      <gravity_y>0</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>
      <extrude_mesh name="1">
        <layers_number>5</layers_number>
        <reason>1.1</reason>
        <selector>Outlet or Face1</selector>
        <thickness>0.0025</thickness>
      </extrude_mesh>
    </extrusion>
    <faces_cutting status="off"/>
    <joining>
      <face_joining name="1">
        <fraction>0.5</fraction>
        <plane>45.0</plane>
        <selector>all[]</selector>
        <verbosity>1</verbosity>
        <visualization>1</visualization>
      </face_joining>
    </joining>
    <mesh_cartesian>
      <x_direction law="constant" max="1.0" min="0.0" ncells="1" prog="1.0"/>
      <y_direction law="constant" max="1.0" min="0.0" ncells="1" prog="1.0"/>
      <z_direction law="constant" max="1.0" min="0.0" ncells="1" prog="1.0"/>
    </mesh_cartesian>
    <mesh_smoothing status="off"/>
    <meshes_list>
      <mesh name="TuboM2.med" path="/home/lgarelli/RUNS_local/Code_Saturne/Forum_cases/Thermal_coupling"/>
    </meshes_list>
    <periodicity/>
    <thin_walls/>
    <volumic_conditions>
      <zone groundwater_law="off" head_losses="off" id="1" initialization="on" label="all_cells" momentum_source_term="off" physical_properties="on" porosity="off" scalar_source_term="off" solid="off" thermal_source_term="off">all[]</zone>
      <zone groundwater_law="off" head_losses="off" id="2" initial_value="off" initialization="on" label="Zone_fluid" momentum_source_term="off" physical_properties="on" porosity="off" scalar_source_term="off" solid="off" thermal_source_term="off">Fluid</zone>
      <zone groundwater_law="off" head_losses="off" id="3" initial_value="off" initialization="on" label="Zone_solid" momentum_source_term="off" physical_properties="on" porosity="off" scalar_source_term="off" solid="on" thermal_source_term="off">Solid</zone>
    </volumic_conditions>
  </solution_domain>
  <thermophysical_models>
    <ale_method/>
    <atmospheric_flows model="off">
      <large_scale_meteo status="off"/>
    </atmospheric_flows>
    <compressible_model model="off"/>
    <conjugate_heat_transfer>
      <external_coupling>
        <syrthes_instances/>
      </external_coupling>
    </conjugate_heat_transfer>
    <gas_combustion model="off">
      <soot_model model="off"/>
      <thermodynamical_pressure status="off"/>
    </gas_combustion>
    <groundwater_model model="off"/>
    <hgn_model model="off"/>
    <internal_coupling>
      <coupled_scalars>
        <scalar name="temperature"/>
      </coupled_scalars>
      <solid_zones>
        <zone name="Zone_2"/>
        <zone name="Zone_solid"/>
      </solid_zones>
    </internal_coupling>
    <joule_effect model="off"/>
    <porosities/>
    <radiative_transfer model="off">
      <absorption_coefficient type="constant">0</absorption_coefficient>
      <restart status="on"/>
    </radiative_transfer>
    <reference_values>
      <length/>
    </reference_values>
    <solid_fuels model="off"/>
    <source_terms/>
    <thermal_scalar model="temperature_celsius">
      <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"/>
      <property label="Thermal flux" name="thermal_flux" support="boundary"/>
      <property label="Tplus" name="tplus" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <variable label="TempC" name="temperature" type="thermal">
        <formula zone_id="1">temperature = 35.;</formula>
        <formula zone_id="3">temperature = 35.;</formula>
        <formula zone_id="2">temperature = 35.;</formula>
        <max_value>100</max_value>
        <min_value>20</min_value>
        <rhs_reconstruction>1</rhs_reconstruction>
        <time_step_factor>10</time_step_factor>
        <turbulent_flux_model>SGDH</turbulent_flux_model>
      </variable>
    </thermal_scalar>
    <turbomachinery model="off">
      <joining/>
    </turbomachinery>
    <turbulence model="k-omega-SST">
      <gravity_terms status="off"/>
      <initialization choice="reference_value" zone_id="1"/>
      <initialization choice="reference_value" zone_id="2"/>
      <initialization choice="reference_value" zone_id="3"/>
      <property label="TurbVisc" name="turbulent_viscosity"/>
      <reference_velocity>0.1</reference_velocity>
      <variable label="k" name="k">
        <blending_factor>0</blending_factor>
        <order_scheme choice="upwind"/>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
      <variable label="omega" name="omega"/>
      <wall_function>7</wall_function>
    </turbulence>
    <velocity_pressure>
      <initialization>
        <formula zone_id="1">velocity[0] = 0.1;
velocity[1] = 0.;
velocity[2] = 0.;</formula>
        <formula zone_id="2">velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;</formula>
        <formula zone_id="3">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>2</rhs_reconstruction>
      </variable>
      <variable dimension="3" label="Velocity" name="velocity">
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
    </velocity_pressure>
  </thermophysical_models>
</Code_Saturne_GUI>