<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="RANS" solver_version="5.0" study="Mixing_Tee" version="2.0">
	<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 name="1" status="on">
				<probe_x>0.2802416</probe_x>
				<probe_y>-0.000360228</probe_y>
				<probe_z>0.06927766</probe_z>
			</probe>
			<probe name="2" status="on">
				<probe_x>0.279949</probe_x>
				<probe_y>0.06919297</probe_y>
				<probe_z>-0.000650659</probe_z>
			</probe>
			<probe name="3" status="on">
				<probe_x>0.2805034</probe_x>
				<probe_y>-0.000568508</probe_y>
				<probe_z>-0.06933761</probe_z>
			</probe>
			<probe name="4" status="on">
				<probe_x>0.2797282</probe_x>
				<probe_y>-0.690437</probe_y>
				<probe_z>-0.000258354</probe_z>
			</probe>
			<probe name="5" status="on">
				<probe_x>0.5613663</probe_x>
				<probe_y>-0.000248044</probe_y>
				<probe_z>0.06927146</probe_z>
			</probe>
			<probe name="6" status="on">
				<probe_x>0.5607841</probe_x>
				<probe_y>0.06920226</probe_y>
				<probe_z>-0.001397463</probe_z>
			</probe>
			<probe name="7" status="on">
				<probe_x>0.5583545</probe_x>
				<probe_y>-0.000566083</probe_y>
				<probe_z>-0.06933812</probe_z>
			</probe>
			<probe name="8" status="on">
				<probe_x>0.5605911</probe_x>
				<probe_y>-0.06904431</probe_y>
				<probe_z>-0.000937696</probe_z>
			</probe>
			<probe_format choice="DAT"/>
			<probe_recording_frequency>1</probe_recording_frequency>
			<writer id="-1" label="results">
				<directory name="postprocessing"/>
				<format name="med" options=""/>
				<frequency period="time_step">100</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>3000</iterations>
			<relaxation_coefficient>0.7</relaxation_coefficient>
			<zero_iteration status="off"/>
		</steady_management>
		<time_averages/>
		<time_parameters>
			<iterations>10</iterations>
			<time_passing>0</time_passing>
			<time_step_ref>0.1</time_step_ref>
		</time_parameters>
	</analysis_control>
	<boundary_conditions>
		<boundary label="BC_1" name="1" nature="inlet">inlet_1</boundary>
		<boundary label="BC_2" name="2" nature="inlet">inlet_2</boundary>
		<boundary label="BC_3" name="3" nature="outlet">outlet</boundary>
		<boundary label="BC_4" name="4" nature="wall">wall</boundary>
		<inlet field_id="none" label="BC_1">
			<scalar choice="dirichlet" name="temperature" type="thermal">
				<dirichlet>19</dirichlet>
			</scalar>
			<turbulence choice="hydraulic_diameter">
				<hydraulic_diameter>0.14</hydraulic_diameter>
			</turbulence>
			<velocity_pressure choice="norm" direction="normal">
				<norm>0.585</norm>
			</velocity_pressure>
		</inlet>
		<inlet field_id="none" label="BC_2">
			<scalar choice="dirichlet" name="temperature" type="thermal">
				<dirichlet>36</dirichlet>
			</scalar>
			<turbulence choice="hydraulic_diameter">
				<hydraulic_diameter>0.1</hydraulic_diameter>
			</turbulence>
			<velocity_pressure choice="norm" direction="normal">
				<norm>0.764</norm>
			</velocity_pressure>
		</inlet>
		<outlet field_id="none" label="BC_3">
			<scalar choice="neumann" name="temperature" type="thermal">
				<neumann>0</neumann>
			</scalar>
		</outlet>
		<wall field_id="none" label="BC_4">
			<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_rescue>0</restart_rescue>
			<restart_with_auxiliary 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="simple">
			<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="variable" label="Density" name="density">
				<formula>a = 1.4078e-5;
b = -5.5855e-3;
c = -2.8886e-3;
d = 1000.4;
density = temperature*(temperature*(a*temperature+b)+c)+d;</formula>
				<initial_value>996.0697</initial_value>
			</property>
			<property choice="variable" label="LamVisc" name="molecular_viscosity">
				<formula>a = -1.9296e-9;
b = 4.7256e-7;
c = -4.2088e-5;
d = 1.6947e-3;
molecular_viscosity = temperature*(temperature*(a*temperature+b)+c)+d;</formula>
				<initial_value>0.000877</initial_value>
			</property>
			<property choice="variable" label="SpecHeat" name="specific_heat">
				<formula>a = -1.0224e-4;
b = 2.9201e-2;
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;</formula>
				<initial_value>4181.252</initial_value>
			</property>
			<property choice="variable" label="ThermalCond" name="thermal_conductivity">
				<formula>a = -1.0224e-4;
b = 2.9201e-2;
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a2 = -3.0374e-13;
b2 = -2.1701e-9;
c2 = 4.7970e-7;
d2 = 1.3538e-4;
thermal_conductivity = (temperature*(temperature*(a*temperature+b)+c)+d)*specific_heat;</formula>
				<initial_value>0.613653</initial_value>
			</property>
			<reference choice="user_material"/>
		</fluid_properties>
		<gravity>
			<gravity_x>0</gravity_x>
			<gravity_y>0</gravity_y>
			<gravity_z>-9.81</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="Mesh_1.unv" path="C:/Users/stichko/Documents/Simulations/Code_Saturne/Mixing_Tee/MESH"/>
		</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>
		<fans/>
		<gas_combustion model="off" option="off"/>
		<groundwater_model model="off"/>
		<joule_effect model="off"/>
		<radiative_transfer model="off">
			<absorption_coefficient type="constant">0</absorption_coefficient>
			<restart status="off"/>
		</radiative_transfer>
		<reference_values>
			<pressure>101325</pressure>
			<temperature>19</temperature>
			<velocity>0.5</velocity>
		</reference_values>
		<solid_fuels model="off"/>
		<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">
				<blending_factor>1</blending_factor>
				<formula zone_id="1">temperature = 19.;</formula>
				<max_value>36</max_value>
				<min_value>19</min_value>
				<rhs_reconstruction>1</rhs_reconstruction>
				<solver_precision>1e-05</solver_precision>
				<turbulent_flux_model>SGDH</turbulent_flux_model>
			</variable>
		</thermal_scalar>
		<turbomachinery model="off">
			<joining/>
		</turbomachinery>
		<turbulence model="k-epsilon-PL">
			<gravity_terms status="on"/>
			<initialization choice="reference_value" zone_id="1"/>
			<property label="TurbVisc" name="turbulent_viscosity"/>
			<variable label="epsilon" name="epsilon">
				<blending_factor>0</blending_factor>
				<rhs_reconstruction>1</rhs_reconstruction>
				<solver_precision>1e-05</solver_precision>
			</variable>
			<variable label="k" name="k">
				<blending_factor>0</blending_factor>
				<rhs_reconstruction>1</rhs_reconstruction>
				<solver_precision>1e-05</solver_precision>
			</variable>
			<wall_function>3</wall_function>
		</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">
				<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>
				<solver_precision>1e-05</solver_precision>
			</variable>
			<variable dimension="3" label="Velocity" name="velocity">
				<blending_factor>1</blending_factor>
				<rhs_reconstruction>1</rhs_reconstruction>
				<solver_precision>1e-05</solver_precision>
			</variable>
		</velocity_pressure>
	</thermophysical_models>
</Code_Saturne_GUI>