Manual DICOM
Manual DICOM
Manual DICOM
DICOM V5.3
Contents
1
INTRODUCTION ..................................................................................................................... 2
2.2
STEADY .......................................................................................................................... 3
2.2.1
General .................................................................................................................. 3
2.2.2
Numerics ............................................................................................................... 4
2.2.3
Morphology ........................................................................................................... 5
2.2.4
2.2.5
2.2.6
Example ............................................................................................................... 10
2.3
TRANSIENT .................................................................................................................. 11
2.3.1
General ................................................................................................................ 11
2.3.2
Numerics ............................................................................................................. 11
2.3.3
Discretization....................................................................................................... 12
2.3.4
Morphology ......................................................................................................... 14
2.3.5
2.3.6
2.3.7
Ejemplo................................................................................................................ 20
Relst.dat .............................................................................................................. 21
3.1.2
Integ.dat .............................................................................................................. 22
3.2
STEADY ........................................................................................................................ 23
3.2.1
3.3
Estacionario.dat................................................................................................... 23
TRANSIENT .................................................................................................................. 25
3.3.1
Temp1.dat ........................................................................................................... 26
3.3.2
Temp2.dat ........................................................................................................... 27
3.3.3
Dat.dat ................................................................................................................. 28
1 INTRODUCTION
The developed model possesses a main screen to enter input data. All the
results, with fuels characteristics are written in a file.dat. This type of file allows the
utilization of different software like Excel or Grapher to plot the data.
In the next chapter, we will explain more accurately the parameters needed to
create a test model and how to enter them in the software. And in the third and last
chapter, we will explain all the output data and the structure of output files.
2 INPUT DATA
In the next paragraph we will explain all the input interfaces screens and the
data and the corresponding data, which can be as well entered by the keyboard or with
a text file. But, we have to distinguish the two cases (steady and transient) in the data
introduction because its a different data processing routine.
Nevertheless, the input interface is similar in the both cases, sharing a lot of
data. In the both cases, there are three buttons on the right of the screen which are:
- Open: allows opening input data of a previous test which is kept on a text file.
- Save: allows saving test results in a text file.
- Start: starts the calculation.
DICOM_v53_C7H16.exe: n-heptane
DICOM_v53_C10H22.exe: decane
DICOM_v53_C12H26.exe: dodecane
DICOM_v53_C16H36.exe: hexadecane
DICOM_v53_DIESEL.exe: diesel
2.2 STEADY
The steady test studies one case in which, injection rate and chamber
properties are constants. Therefore, there will not be any saved files for injection rate,
density or pressure in function of the time. But all this information will interfere directly
from the interface of entry.
Another of the main differences will be explained hereinafter in the results of
exit, since this one does not possess temporary files.
The execution, as the user will be able to see, will be faster than the transient case, in
which the conditions change constantly.
2.2.1
General
The first tag allows selecting if the model is transient or steady. The only field on
this screen is the path where the data will be saved:
- Directory of output Files: In this field we have to indicate the address in which the
program will write the results. Moreover, there is a button, whose the name is
Browse, which allows to search in the folders.
Figure I.3. Interface of information of DICOM's entry, 'GENERAL' tab of steady model.
For the steady case it will not be necessary to provide additional information.
2.2.2
Numerics
In this tag we have to introduce parameters for discretization that are already in
exit tables or calculations. These parameters are:
- Maximum length for calculation [m]: Define the maximum length for calculation in
which it will come near during the study of the flame. We can choose a characteristic
length like the diameter of combustion chamber.
- Increase in x of the discretization dx [m]: provides the equi spread (dx) in the
spatial discretization of the spray. The smaller the value of dx, the longer it will take for
the whole calculation to finish.
- Mass fraction increase for State Relationship: Length between points in the table
of data for the steady model. We have to take into count that mass fraction data are
between 0 and 1.
4
- Mass fraction increase for integral tables: Like the previous parameter,
discretization of mass fraction but used in calculation and data storage of integrals
calculated by the program.
- Convergence boundary for iterative computation: Indicative value for calculation
end, its meaning, if the difference between results of 2 successive computations is less
than this limit value, the calculation stops.
Moreover, we can see a picture of the screen that is described with values that,
by default, the program introduces.
Figure I.4. Interface of information of DICOM's entry, 'NUMERICS' tab of steady model.
2.2.3
Morphology
In the next tag there are some parameters that define the spray, like geometry, speed
profile and distribution of mass fraction.
- Velocity fraction for spray boundary: Velocity limit at which it will be considered
there is no spray for being too low. It constitutes a way of delimiting the jet.
- Schmidt number: The Schmidt number is in the ratio of momentum diffusivity
(viscosity) and mass diffusivity. The value of the Schmidt number is normally 1 but we
can change it by a number between 0 and 1.
- Spray angle []: In this field we have to introduce the cone angle to define the steady
spray.
- Radial profile: This parameter allows choosing one of the four profiles available
(Exponential, Spalding, Hinze y Schlichting) to calculate speed and radial mass fraction
from its values in the central axis. The exponential-Gauss profile is the more used for
modeling.
Figure I.5. Interface of information of DICOM's entry, 'MORPHOLOGY' tab of steady model.
2.2.4
Injection rate
To introduce injection rate we can choose between direct input including the
momentum flow and mass flow; and derived input from the mass flow.
Figure I.6. Interface of information of DICOM's entry, 'INJECTION RATE' and 'DIRECT
INPUT' tab of steady model.
Figure I.7. Interface of information of DICOM's entry, 'INJECTION RATE' and 'DERIVED
INPUT' tab of steady model.
2.2.5
Mixing law
Finally, the last tag contains parameters which define the air/fuel mixture that depend
on both characteristics. There are three cases to be chosen:
2.2.5.1 Isothermal spray
In this case there is not combustion or heat Exchange. This is the simplest case which
requires only:
- Air density [kg/m3]: Air density which is in the chamber before de injection.
- Fuel density [kg/m3]: Density of injected fuel.
- Stoichiometric mass fraction: Relation of fuel/air mass necessary to realize the
combustion in stoichiometrics conditions.
Figure I.8. Interface of information of DICOM's entry, 'MIXING LAW' tab of steady model to an
Isothermal Spray.
Figure I.9. Interface of information of DICOM's entry, 'MIXING LAW' tab of steady model for a
Gas test.
Example
Finally, we can see an example of a file input.txt that saves all the input data of in
interface as we have explained earlier.
2.2.6
10
2.3 TRANSIENT
The transient models give more real tests because of the variation of conditions at
each moment, because it will have to include text files in which exists temporal relation
of the written variables. However, in the next paragraph we will see each tag and the
changes in comparison with steady tests.
2.3.1
General
The interface is similar to the steady case, the only difference is that it includes a field
to input the time step to save the temporary results:
- Time interval to saving results [s]: Define time interval where the software will
create a new result file .dat.
Figure I.11. Interface of information of DICOM's entry, 'GENERAL' tab of transient model.
2.3.2
Numerics
11
- Mass fraction increase for integral tables: Like the previous parameter,
discretization of mass fraction but used in calculation and data storage of integrals
calculated by the program.
- Convergence boundary for main equations: Indicative value for calculation end, its
meaning, if the difference between results of 2 successive computations is less than
this limit value, the calculation stops.
- Velocity value when the Spray finishes in the cell [m/s]: Boundary value of the
end of the spray, such cell takes the velocity value of the spray, it will be considered as
the last value of spray.
Figure I.12. Interface of information of DICOM's entry, 'NUMERICS' tab of transient model.
2.3.3
Discretization
In this tag, are defined temporal and spatial criteria of discretization as well as the
relation between both to avoid mistakes of calculation. Moreover its not necessary to
have a very precise spatial and temporal discretization which generates an excessive
number of calculations and take a much longer time. For that, the following parameters
are defined:
First, we will choose the mode selection which will be automatic, already optimized
inside the program or free, defined by user.
12
- Courant Number: This constant defines the relation between spray propagation,
velocity and temporal and spatial discretization.
- Spatial increase: From the previous equation we can find the Courant Number, and
after the only parameter needed to calculate the temporal increase is the spatial
increase. We can do the inverse, for example, give the temporal increase to have the
spatial increase but in both cases, in automatic or free discretization, we need one
more parameter:
- Spatial increase constant: This parameter allows modifications, independent of
mode selection chosen, the spatial increment, multiplying the calculation per the
mentioned constant.
13
2.3.4
Morphology
The parameters introduced in this part are the same as in the steady case. The only
difference between both cases is in the program of transient model which do not
contain the possibility to use different radial profiles to calculate the distribution field of
speed and mass fraction, and it use by default an exponential profile.
- Velocity fraction for spray boundary: Velocity limit at which it will be considered
there is no spray for being too low. It constitutes a way of delimiting the jet.
- Schmidt number: The Schmidt number is in the ratio of momentum diffusivity
(viscosity) and mass diffusivity. The value of the Schmidt number is normally 1 but we
can change it by a number between 0 and 1.
- Spray angle []: In this field we have to introduce the cone angle to define the steady
spray.
14
Figure I.13. Interface of information of DICOM's entry, 'MORPHOLOGY' tab of transient model.
2.3.5
Injection rate.
The only difference between introduction of injection rate in transient and steady test is
that, in transient, we can use variable injection rate along the injection. Therefore, we
will not introduce constant values in the injection rate parameters of mass flow and
momentum, but a table of these values in function of the time. To introduce these
values we will introduce the name and address of the file that has data.
15
Figure I.14. Interface of information of DICOM's entry, 'INJECTION RATE' and 'DIRECT
INPUT' tab of transient model.
16
Figure I.15. Interface of information of DICOM's entry, 'INJECTION RATE' and 'DERIVED
INPUT' tab of transient model.
After, we can see an example of a temporal file: mass flow rate, momentum rate,
pressure and density.
17
2.3.6
Mixing law
For the mixing law too, there is changes in the parameters, like the ambient density file
(in isothermal case), or pressure inside the camber (gases cases or spray not
isothermal).
2.3.6.1 Isothermal spray
In this case there is not combustion or heat Exchange. This is the simplest case which
precise only:
- Air density [kg/m3]: Air density which is in the chamber before de injection.
- Fuel density [kg/m3]: Density of injected fuel.
- Stoichiometric mass fraction: Relation of fuel/air mass necessary to realize the
combustion in stoichiometric conditions.
18
Figure I.17. Interface of information of DICOM's entry, 'INJECTION RATE' and 'MIXING
LAW' tab of transient model.
19
An instant previous to the indicated one, although the mixture fraction is inferior at f LOL
there will be not a reactive zone in the flame. We have to avoid incompatibility in the
choice of values. For example, if we have an oxygen mass fraction null, the spray will
be inert and the other values have to be coherent: fLOL null y tSOC bigger than the final
time of calculation.
Figure I.18. Interface of information of DICOM's entry, 'INJECTION RATE' and 'MIXING
LAW' tab of transient model.
2.3.7
Example
20
3 OUTPUT DATA
Output files created by both types of calculations, steady and transient, are very
different and we will do a chapter for each one of them. All the files are written in the
address introduced in the interface.
But there are two files which are the same for steady or transient. These two files are:
relst.dat and integ.dat.
3.1.1
Relst.dat
In this output file the relative results to states relations are written, temperatures,
densities and mass fraction. Nevertheless, there are many changes in function of the
modeled test:
- Isothermal: In this case its only the mixture fraction and density associated with a
given mixture that are calculated.
- Gaseous jet: If the test is inert the mixture fraction, the density, temperature and
mass fractions of each fuel and gas associated at each mixture fraction are tabulated.
If the test is reactive, all these columns double for the reactive case
- Spray: Like the previous file, the program includes two columns which contain mass
fractions of each fuel and their liquid and vapor phase.
21
3.1.2
Integ.dat
The last file of the steady calculation includes the tables of integrals calculated by the
model. In these, are tabulated with the column of mixture fraction, a group of columns
variable in function of the test type:
- Isothermal: Are calculated integrals of parameters 1, 2, Sc, 1+Sc y 0, besides the
integrals of mass under relative dosage.
- Gaseous jet: Are added, at the previous columns, the integrals of all of the species,
that's meaning, each unburned fuel and gases which take part
- Spray: Finally, to the previous integrals are added the results of the integrals of mass
evaporated of each fuel, liquid and vapor mass of the same one.
22
3.2 STEADY
The steady test creates three files that the program saves and names automatically:
estacionario.dat, relst.dat and integ.dat.
Now, we will describe the estacionario.dat, the other two have already been described
before. Its important to remind that the steady generates generic results, and not for
real instants like in transient conditions, since in an stationary test the jet is supposedly
stable and therefore the results would not change from a moment to other one in the
calculation. Moreover, the number of values depends on the type of test.
3.2.1
Estacionario.dat
This file has the main results of the model, in which there is the following data:
parameter
xeje
unit
m
ucl
fcl
tcx
m/s
s
definition
Position on the central axis, for all values its an axial
position.
Spray velocity on its axis.
Mixture fraction in the spray axis.
Time of the position on the spray axis.
23
Rini_u, Rini_f
R
x/deq
m
-
numiter
kg/s
umed
fmed
rho_cl
rho_med
rho_med_flujo
R_x
m/s
kg/m3
kg/m3
kg/m3
m
Revap_CXHy
mf
kg/m
mfmix_x
kg/m
mf_CXHymixevap
kg/m
ma
S_x
Sevap_CXHy
kg/m
m
m
miny
mfmix_x_ac
kg
kg
mfCXHy_mixevap
_ac
ma_ac
mfCXHy_sq
kg
mfCXHy_l
kg/m
mfCXHy_sq_ac
mf_q_ac
mfCXHy_l_ac,
mfCXHy_v_ac
kg/m
kg/m
kg/m
kg
kg/m
mO2_ac,
mCO2_ac,
mH2O_ac
kg/m
Tcl
Tmed
YfCXHy_cl
YfCXHy_v_cl,
YfCXHy_l_cl
YN2_cl,YO2_cl,Y
CO2_cl,YH2O_cl
K
K
-
3.3 TRANSIENT
In the transient test, we have already seen that the input data could be variable, in
consequence, a study is realized with a temporal follow up which is composed by the
following files: Temp1.dat, Temp2.dat, relest.dat, integ.dat and dat.dat. The first two
files are unique for each test which reference different variables in function of time in
different columns. The integrals file and state relations file are the same as in steady
calculations, with the only difference that the calculation for both is made two times,
one time before the combustion and a second at the beginning of the combustion. And
finally, the temporal files, with interval of time specified in the first screen of input data,
they gather concrete moments of the jet in studding a series of variables around its
position in the axis.
Now, we are going to see each text file:
25
3.3.1
Temp1.dat
Only one file of this type is created, which gathers variables that are used in the
calculation module at each moment. These are the variables of the calculation, with a
column of time variable:
parameter
u0
unit
m/s
definition
Velocity of injection at the exit of the nozzle.
I0
M0
kg/s
rho_a[]
do
Ta
Tst
fLOL
kg/m
3
kg/m
3
m
K
K
-
fevap_CxHy
rho_f
26
3.3.2
Temp2.dat
Variables tabulated in this archive represent some variables at the exit of the loop of
calculation. Are represented in function of the time:
parameter
S
S_x
S_evap_CxHy
S_evap_o_CxHy
unit
m
m
m
m
miny
ma
mfmix_x_ac
kg
kg
kg
mfCxHy_mixevap kg
_ac
definition
Penetration value of the jet.
Penetration value of the surface with relative dosage x.
Penetration value of the surface of evaporation of the fuel.
Penetration value of the surface of evaporation of the fuel,
calculated from the origin.
Fuel mass injected up to a given moment.
Air mass accumulated in the jet up to a given moment.
Fuel mass accumulated up to the moment under the relative
dosage x.
Fuel mass accumulated up to the moment under the dosage
of evaporation of each fuel.
Lift-off-length, distance since the exit of the nozzle until the lift
off of the flame
Mass of the unburned fuel accumulated up to the moment.
Global mass of burned fuel accumulated up to the moment.
Mass of gaseous species which take part, accumulated inside
the jet up to a given moment.
LOL
mfCxHy_sq_ac
mf_q_ac
mO2_ac,
mCO2_ac,
Mh2O_ac
kg
kg
kg
mfCxHy_l,
mfCxHy_v
uclmax
x_uclmax
fclmax
x_fclmax
27
3.3.3
Dat.dat
The program creates files at every time interval, determined by the user to save results.
This file contains variables that are represented in function of the axis position of the
jet, which define an immediate picture of jet.
These variables are:
parameter
uclmax
fcl
unit
m/s
-
definition
I
Mf
INT_u_XR_ant
INT_f_XR_ant
M
N
kg/s
kg*m/s
kg
kg/s
umed
fmed
rho_cl
rho_med
m/s
kg/m3
kg/m3
m
m
R
Rini
R_x
Revap_CXHy
mf
kg/m
mfmix_x
kg/m
mf_CXHymixeva kg/m
p
ma
mfCXHy_sq
kg/m
kg/m
mfCXHy_l,
mfCXHy_v
kg/m
mO2,
mCO2,
mH2O
Tcl
Tmed
YfCXHy_cl
YfCXHy_v_cl,
YfCXHy_l_cl
YN2_cl,
YO2_cl,
YCO2_cl,
YH2O_cl
kg/m
K
K
-
of x.
Radius of each section in which the spray passes from
totally liquid to totally gaseous. There are as many
columns as fuels and if the jet is not a spray the radius
gives itself for the relative dosage 3.
Total fuel quantity that there is by length in each section.
29
30