File mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_input.txt artifact 6836ac416b part of check-in 8ae9b04b7f

# Numerical parameter file (DieselCycle_input.txt)
# Generated by MTT at Thu Dec  4 11:17:09 GMT 1997

# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% Version control history
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% $Id$
# %% $Log$
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


# Set the inputs

if ((t>=0.0)&&(t<1.0))		#Compression
  u(1) = 0.0;			# Entropy flow
  r = 1e10;			# Big R to make a flow source 
  u(2) = 0.8*r;			#- Volume rate-of-change
  TopPressure = (gamma_g-1)*(x(1)/x(2))
endif;

if ((t>=1.0)&&(t<2.0))		#Heating
  u(1) = 1000;			# Entropy flow
  r    = 1e4;			# small r constant pressure
  u(2) = TopPressure;		# Pressure source
  Volume = x(4);
endif;

if ((t>=2.0)&&(t<3.0))		#Expansion
  u(1) = 0.0;			# Entropy flow
  r = 1e10;			# Big R to make a flow source 
  u(2) = -(1-Volume)*r;		# Volume rate-of-change
endif;

if (t>=3.0)			#Cooling
  Pressure = (gamma_g-1)*x(1)/x(2);
  u(1) = -1000*(Pressure>1e5);	# Entropy flow
  r = 1e10;			# Big R to make a flow source 
  u(2) = 0.0*r;			# Volume rate-of-change
endif;