Function set 'adiab' - an Ideal Adiabatic
simulation of a specific Stirling engine configuration
From the flow diagram below we see that four different systems are invoked to
do an Ideal Adiabatic simulation. The main program 'stimadiab'
first defines the system to be simulated in terms of the set of global variables
set up by the 'define'
set of functions, as described previously. It then invokes function 'adiabatic'
which solves the set of differential equations (function 'dadiab')
over a number of cycles until convergence is attained (function 'adiab'),
and then fills in the solution matrix for a complete cycle (function 'filmat')
and finally displays the solution matrix (function 'prntad').
The differential equation set is solved by using the Classical Fourth Order
Runge-Kutta method (function 'rk4')
which is included in the system 'odes' (ordinary differential
equations).
The dynamics of the solution algorithm lies in
function 'adiab',
which initialises the variables, invokes the Runge-Kutta function over a number
of cycles, checks for cyclic convergence, then fills in the solution matrix. The
method of uniquely specifying the variables of solution is done by means of
symbolic constants in the header file 'adiabatic.h'.
Notice that the function 'volume'
includes only sinusoidal volume variations (function 'sinevl') and the Ross Yoke-drive volume variations (function
'yokevl'). As before, it is intended that the user
will modify and augment this system as required for specific systems.
Furthermore, to do the various energy and temperature plots in this chapter we
replaced the module 'prntad' with a function 'plotad' to create data files for off-line plotting.
__________________________________________________________________________
__________________________________________________________________________
On to the Simple
Analysis
Back to the Stirling Cycle
Machines home page