<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="SpalAl" solver_version="5.0;6.0" study="Stab" version="2.0">
  <additional_scalars>
    <users/>
    <variable label="BFG" name="BFG" type="user">
      <formula zone_id="1">BFG = 0;</formula>
      <max_value>1</max_value>
      <min_value>0</min_value>
      <property choice="constant" label="Dscal1" name="BFG_diffusivity">
        <formula>BFG_diffusivity = 1.e-5;</formula>
        <initial_value>1e-05</initial_value>
      </property>
      <rhs_reconstruction>1</rhs_reconstruction>
      <solver_precision>1e-05</solver_precision>
      <turbulent_flux_model>SGDH</turbulent_flux_model>
    </variable>
    <variable label="O2" name="O2" type="user">
      <formula zone_id="1">O2 = 0;</formula>
      <max_value>1</max_value>
      <min_value>0</min_value>
      <property choice="constant" label="Dscal2" name="O2_diffusivity">
        <initial_value>1e-05</initial_value>
      </property>
      <rhs_reconstruction>1</rhs_reconstruction>
      <solver_precision>1e-05</solver_precision>
      <turbulent_flux_model>SGDH</turbulent_flux_model>
    </variable>
    <variable label="N2" name="N2" type="user">
      <formula zone_id="1">N2 = 0;</formula>
      <max_value>1</max_value>
      <min_value>0</min_value>
      <property choice="constant" label="Dscal3" name="N2_diffusivity">
        <initial_value>1e-05</initial_value>
      </property>
      <rhs_reconstruction>1</rhs_reconstruction>
      <solver_precision>1e-05</solver_precision>
      <turbulent_flux_model>SGDH</turbulent_flux_model>
    </variable>
  </additional_scalars>
  <analysis_control>
    <error_estimator/>
    <output>
      <listing_printing_frequency>200</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>10</probe_recording_frequency>
      <writer id="-1" label="results">
        <directory name="postprocessing"/>
        <format name="ensight" options="separate_meshes"/>
        <frequency period="time_step">1000</frequency>
        <output_at_end status="on"/>
        <output_at_start status="on"/>
        <time_dependency choice="fixed_mesh"/>
      </writer>
    </output>
    <profiles/>
    <scalar_balances/>
    <steady_management status="on">
      <iterations>10</iterations>
      <relaxation_coefficient>0.7</relaxation_coefficient>
      <zero_iteration status="off"/>
    </steady_management>
    <time_averages/>
    <time_parameters>
      <iterations>100000</iterations>
      <max_courant_num>1</max_courant_num>
      <max_fourier_num>10</max_fourier_num>
      <property label="CourantNb" name="courant_number"/>
      <property label="FourierNb" name="fourier_number"/>
      <property label="LocalTime" name="local_time_step"/>
      <relaxation_coefficient>0.7</relaxation_coefficient>
      <time_passing>2</time_passing>
      <time_step_max_factor>1000</time_step_max_factor>
      <time_step_min_factor>0.1</time_step_min_factor>
      <time_step_ref>5e-06</time_step_ref>
      <time_step_var>0.1</time_step_var>
    </time_parameters>
  </analysis_control>
  <boundary_conditions>
    <boundary label="BC_1" name="1" nature="wall">AXIAL_1_3_2D</boundary>
    <boundary label="BC_2" name="2" nature="wall">AXIAL_1_4_2D</boundary>
    <boundary label="BC_3" name="3" nature="wall">AXIAL_1_7_2D</boundary>
    <boundary label="BC_4" name="4" nature="wall">AXIAL_1_8_2D</boundary>
    <boundary label="BC_5" name="5" nature="wall">AXIAL_1_9_2D</boundary>
    <boundary label="BC_6" name="6" nature="inlet">INLET_GAZ_2D</boundary>
    <boundary label="BC_7" name="7" nature="inlet">INLET_O2_2D</boundary>
    <boundary label="BC_8" name="8" nature="inlet">INLET_TRANSPORT_2D</boundary>
    <boundary label="BC_9" name="9" nature="imposed_p_outlet">OUTLET_VENT_2D</boundary>
    <boundary label="BC_10" name="10" nature="imposed_p_outlet">OUTLET_VENT_RESISTANT_2D</boundary>
    <boundary label="BC_11" name="11" nature="wall">WALL_INT_TUYERE_2D</boundary>
    <boundary label="BC_12" name="12" nature="wall">WALL_REFROIDIES_2D</boundary>
    <imposed_p_outlet field_id="none" label="BC_9">
      <dirichlet name="pressure">340000</dirichlet>
    </imposed_p_outlet>
    <imposed_p_outlet field_id="none" label="BC_10">
      <dirichlet name="pressure">340000</dirichlet>
    </imposed_p_outlet>
    <inlet field_id="none" label="None">
      <scalar choice="dirichlet" label="BFG" name="BFG" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" label="O2" name="O2" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" label="N2" name="N2" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <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="BC_6">
      <scalar choice="dirichlet" name="BFG" type="user">
        <dirichlet>1</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="O2" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="N2" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>0.18</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <fraction>0</fraction>
        <gas_type choice="fuel"/>
        <norm>32.9262</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="BC_7">
      <scalar choice="dirichlet" name="BFG" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="O2" type="user">
        <dirichlet>1</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="N2" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>0.0316</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <fraction>1</fraction>
        <gas_type choice="oxydant"/>
        <norm>89.3665</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="BC_8">
      <scalar choice="dirichlet" name="BFG" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="O2" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="N2" type="user">
        <dirichlet>1</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>0.0156</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <flow1>26.9205</flow1>
        <fraction>1</fraction>
        <gas_type choice="oxydant"/>
        <norm>26.92</norm>
      </velocity_pressure>
    </inlet>
    <wall field_id="none" label="BC_1">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="BC_2">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="BC_3">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="BC_4">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="BC_5">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="BC_11">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
    <wall field_id="none" label="BC_12">
      <scalar choice="neumann" name="BFG" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="O2" type="user">
        <neumann>0</neumann>
      </scalar>
      <scalar choice="neumann" name="N2" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off"/>
    </wall>
  </boundary_conditions>
  <calculation_management>
    <block_io/>
    <partitioning/>
    <start_restart>
      <frozen_field status="off"/>
      <restart path="/gpfs1/home/2006003/cbarna01/IGAR_project/7045908/FastNew/NewG/RESU/20200428-0933/checkpoint"/>
      <restart_rescue>0</restart_rescue>
      <restart_with_auxiliary status="off"/>
    </start_restart>
  </calculation_management>
  <lagrangian model="off"/>
  <numerical_parameters>
    <density_relaxation>0.95</density_relaxation>
    <gradient_reconstruction choice="4"/>
    <gradient_transposed status="on"/>
    <hydrostatic_pressure status="on"/>
    <pressure_relaxation>1</pressure_relaxation>
    <velocity_pressure_algo choice="simplec">
      <piso_sweep_number>1</piso_sweep_number>
    </velocity_pressure_algo>
    <velocity_pressure_coupling status="on">
      <property label="VPsolve1" name="weight_matrix_X">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <property label="VPsolve2" name="weight_matrix_Y">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <property label="VPsolve3" name="weight_matrix_Z">
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
    </velocity_pressure_coupling>
    <wall_pressure_extrapolation>0</wall_pressure_extrapolation>
  </numerical_parameters>
  <physical_properties>
    <fluid_properties>
      <material choice="user_material"/>
      <method choice="user_properties"/>
      <property choice="predefined_law" label="Density" name="density">
        <formula>density = 1.0;</formula>
        <initial_value>1</initial_value>
      </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">
        <formula>CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
temperature = enthalpy / 1040.;

if ( enthalpy &gt; 0) {
molecular_viscosity = mu_ref * (T0+CST / temperature+CST) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}</formula>
        <initial_value>1e-05</initial_value>
      </property>
      <property choice="constant" label="SpecHeat" name="specific_heat">
        <formula>specific_heat = -1.;</formula>
        <initial_value>1017.24</initial_value>
      </property>
      <property choice="constant" label="ThermalCond" name="thermal_conductivity">
        <formula>temperature = enthalpy / 1040.;
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;</formula>
        <initial_value>0.01</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <reference choice="user_material"/>
      <reference_fuel_temperature>1223.13</reference_fuel_temperature>
      <reference_molar_mass>0.028966</reference_molar_mass>
      <reference_oxydant_temperature>298.13</reference_oxydant_temperature>
      <reference_pressure>350000</reference_pressure>
      <reference_temperature>300</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/>
    <faces_cutting status="off"/>
    <joining/>
    <mesh_smoothing status="off"/>
    <meshes_list>
      <mesh name="TuyereAxiale.med"/>
    </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>
    </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>
    <gas_combustion model="d3p" option="adiabatic">
      <data_file>dp_C3P_RIGHT</data_file>
      <property label="temperature" name="temperature"/>
      <property label="ym_fuel" name="ym_fuel"/>
      <property label="ym_oxyd" name="ym_oxyd"/>
      <property label="ym_prod" name="ym_prod"/>
      <variable label="mixture_fraction" name="mixture_fraction" type="model">
        <blending_factor>0</blending_factor>
        <order_scheme choice="upwind"/>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
      <variable label="mixture_fraction_variance" name="mixture_fraction_variance" type="model">
        <blending_factor>0</blending_factor>
        <order_scheme choice="upwind"/>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
    </gas_combustion>
    <groundwater_model model="off"/>
    <head_losses/>
    <hgn_model model="off"/>
    <joule_effect model="off"/>
    <porosities/>
    <radiative_transfer model="off">
      <absorption_coefficient type="modak">0</absorption_coefficient>
      <directions_number>3</directions_number>
      <frequency>1</frequency>
      <intensity_resolution_listing_printing>0</intensity_resolution_listing_printing>
      <property label="Luminance" name="luminance">
        <listing_printing status="off"/>
      </property>
      <property label="Absorption" name="rad_absorption">
        <listing_printing status="off"/>
      </property>
      <property label="Absorption_coefficient" name="rad_absorption_coeff">
        <listing_printing status="off"/>
      </property>
      <property label="Coeff_ech_convectif" name="rad_exchange_coefficient" support="boundary"/>
      <property label="Emission" name="rad_emission">
        <listing_printing status="off"/>
      </property>
      <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">
        <listing_printing status="off"/>
      </property>
      <property label="Radiative_source_term" name="radiative_source_term">
        <listing_printing status="off"/>
      </property>
      <property label="Srad" name="radiative_st">
        <listing_printing status="off"/>
      </property>
      <property label="Thickness" name="thickness" support="boundary"/>
      <property label="Thermal_conductivity" name="wall_thermal_conductivity" support="boundary"/>
      <property label="Thickness" name="wall_thickness" 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="off">
      <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"/>
    </thermal_scalar>
    <turbomachinery model="off">
      <joining/>
    </turbomachinery>
    <turbulence model="Spalart-Allmaras">
      <gravity_terms status="off"/>
      <initialization choice="formula" 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.2</reference_length>
      <reference_velocity>30</reference_velocity>
      <variable label="nu_tilda" name="nu_tilda"/>
    </turbulence>
    <velocity_pressure>
      <initialization>
        <formula zone_id="1">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"/>
      <property label="Stress, tangential" name="stress_tangential" support="boundary"/>
      <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>