<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="CASE-201" study="21_PISTON" version="2.0">
	<solution_domain>
		<volumic_conditions>
			<zone head_losses="off" initialization="on" label="all_cells" mass_source_term="off" momentum_source_term="off" name="1" scalar_source_term="off" thermal_source_term="off">all[]</zone>
		</volumic_conditions>
		<meshes_list>
			<mesh format="ideas" name="pistonfinal.unv"/>
		</meshes_list>
		<join_meshes status="off"/>
		<faces_cutting status="off"/>
		<reorientation status="off"/>
		<periodic_boundary/>
		<standalone/>
	</solution_domain>
	<thermophysical_models>
		<velocity_pressure>
			<variable label="Pression" name="pressure">
				<reference_pressure>101325</reference_pressure>
				<solveur_precision>1e-05</solveur_precision>
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
			</variable>
			<variable label="VitesseX" name="velocity_U">
				<blending_factor>1</blending_factor>
				<solveur_precision>1e-05</solveur_precision>
			</variable>
			<variable label="VitesseY" name="velocity_V">
				<blending_factor>1</blending_factor>
				<solveur_precision>1e-05</solveur_precision>
			</variable>
			<variable label="VitesseZ" name="velocity_W">
				<blending_factor>1</blending_factor>
				<solveur_precision>1e-05</solveur_precision>
			</variable>
			<property label="total_pressure" name="total_pressure"/>
			<property label="Yplus" name="yplus" support="boundary"/>
			<property label="Efforts" name="effort" support="boundary"/>
			<property label="all_variables" name="all_variables" support="boundary"/>
		</velocity_pressure>
		<ale_method status="on">
			<variable label="mesh_u" name="mesh_velocity_U">
				<solveur_precision>1e-05</solveur_precision>
			</variable>
			<variable label="mesh_v" name="mesh_velocity_V">
				<solveur_precision>1e-05</solveur_precision>
			</variable>
			<variable label="mesh_w" name="mesh_velocity_W">
				<solveur_precision>1e-05</solveur_precision>
			</variable>
			<property choice="user_function" label="mesh_vi1" name="mesh_viscosity_1">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
			</property>
			<mesh_viscosity type="isotrop"/>
			<fluid_initialization_sub_iterations>0</fluid_initialization_sub_iterations>
			<formula>mesh_vi1 = 1;</formula>
			<max_iterations_implicitation>1</max_iterations_implicitation>
			<implicitation_precision>1e-05</implicitation_precision>
			<external_coupling_post_synchronization status="off"/>
		</ale_method>
		<turbulence model="off">
			<initialization choice="reference_velocity">
				<reference_velocity>1</reference_velocity>
			</initialization>
		</turbulence>
		<thermal_scalar model="off"/>
		<radiative_transfer model="off"/>
		<conjugate_heat_transfer>
			<external_coupling status="off"/>
			<external_coupling status="off"/>
			<external_coupling status="on"/>
		</conjugate_heat_transfer>
		<gas_combustion model="off"/>
		<pulverized_coal model="off"/>
		<joule_effect model="off"/>
		<atmospheric_flows model="off"/>
	</thermophysical_models>
	<numerical_parameters>
		<multigrid status="on"/>
		<gradient_transposed status="off"/>
		<velocity_pressure_coupling status="off"/>
		<pressure_relaxation>1</pressure_relaxation>
		<wall_pressure_extrapolation>0</wall_pressure_extrapolation>
		<gradient_reconstruction choice="0"/>
	</numerical_parameters>
	<physical_properties>
		<fluid_properties>
			<property choice="constant" label="MasseVol" name="density">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
				<initial_value>1000</initial_value>
			</property>
			<property choice="constant" label="ViscLam" name="molecular_viscosity">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
				<initial_value>0.004</initial_value>
			</property>
			<property choice="constant" label="ChSpec" name="specific_heat">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
				<initial_value>4177</initial_value>
			</property>
			<property choice="constant" label="CondTherm" name="thermal_conductivity">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
				<initial_value>0.02495</initial_value>
			</property>
		</fluid_properties>
		<gravity>
			<gravity_x>0</gravity_x>
			<gravity_y>0</gravity_y>
			<gravity_z>0</gravity_z>
		</gravity>
		<hydrostatic_pressure status="off"/>
	</physical_properties>
	<additional_scalars/>
	<boundary_conditions>
		<boundary label="inlet" name="1" nature="wall">7</boundary>
		<boundary label="outlet" name="2" nature="outlet">5</boundary>
		<boundary label="wall6" name="3" nature="wall">6</boundary>
		<boundary label="sym1" name="4" nature="symmetry">1 or 4</boundary>
		<boundary label="wall3" name="5" nature="wall">3</boundary>
		<outlet label="outlet"/>
		<wall label="wall6">
			<velocity_pressure choice="off"/>
			<ale choice="fixed_boundary"/>
		</wall>
		<wall label="wall3">
			<velocity_pressure choice="off"/>
			<ale choice="fixed_boundary"/>
		</wall>
		<wall label="inlet">
			<ale choice="fixed_velocity">
				<formula>TT    = iter * dt;
TTDC=3.265434084;
TBDC=4.200151065;
PERIO = TBDC;
TIME2 = mod(TT+PERIO,PERIO);

if (TIME2 &lt;= 0.05)
{
    VPISTON=(124.4*TIME2)*1.0e-03;
}
else if ((TIME2 &gt; 0.05) &amp;&amp; (TIME2 &lt;= (TTDC -0.05)))
{
    VPISTON=6.22e-03;
}
else if ((TIME2 &gt; (TTDC -0.05)) &amp;&amp; (TIME2 &lt;= TTDC))
{
    VPISTON=(-124.4*(TIME2 - TTDC))*1.0e-03;
}
else if ((TIME2 &gt; TTDC) &amp;&amp; (TIME2 &lt;= (TTDC+0.18)))
{
    VPISTON=(-147.222*(TIME2 - TTDC))*1.0e-03;
}
else if ((TIME2 &gt; (TTDC + 0.18)) &amp;&amp; (TIME2 &lt;= (TBDC -0.18)))
{
    VPISTON=-26.5e-03;
}
else if (TIME2 &gt; (TBDC -0.18))
{
    VPISTON=(147.222*(TIME2 - TBDC))*1.0e-03;
}

mesh_w = -VPISTON;
mesh_u=0; mesh_v=0;</formula>
			</ale>
			<velocity_pressure choice="off"/>
		</wall>
		<wall label="outlet">
			<ale choice="fixed_boundary"/>
		</wall>
		<boundary label="wall2" name="6" nature="wall">2</boundary>
		<wall label="wall2">
			<velocity_pressure choice="off"/>
			<ale choice="sliding_boundary"/>
		</wall>
		<boundary label="sym2" name="7" nature="symmetry">11 or 41</boundary>
		<symmetry label="sym1"/>
		<wall label="sym1">
			<ale choice="fixed_boundary"/>
		</wall>
		<symmetry label="sym2"/>
		<wall label="sym2">
			<ale choice="sliding_boundary"/>
		</wall>
		<variable/>
		<scalar/>
	</boundary_conditions>
	<analysis_control>
		<output>
			<postprocessing_mesh_options choice="10"/>
			<syrthes_boundary status="off"/>
			<listing_printing_frequency>1</listing_printing_frequency>
			<postprocessing_frequency>10</postprocessing_frequency>
			<probe_recording_frequency>1</probe_recording_frequency>
			<fluid_domain status="on"/>
			<domain_boundary status="off"/>
			<postprocessing_format choice="EnSight"/>
			<postprocessing_options choice="binary"/>
		</output>
		<time_parameters>
			<time_step_ref>0.01</time_step_ref>
			<iterations>1000</iterations>
			<time_passing>0</time_passing>
			<property label="NbCourant" name="courant_number">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
			</property>
			<property label="NbFourier" name="fourier_number">
				<listing_printing status="off"/>
				<postprocessing_recording status="off"/>
			</property>
			<zero_time_step status="off"/>
		</time_parameters>
		<steady_management status="off"/>
		<time_averages/>
		<profiles/>
	</analysis_control>
	<calcul_management>
		<start_restart>
			<restart status="off"/>
		</start_restart>
	</calcul_management>
	<lagrangian model="off"/>
</Code_Saturne_GUI>