Difference between revisions of "FlameModelling"

From SYSOS

 
(2 intermediate revisions by the same user not shown)
Line 5: Line 5:
 
[[Image:Figure1.JPG]]
 
[[Image:Figure1.JPG]]
  
Three files can be used to simulate the above response.
+
Three files can be used to simulate the response of this system.
 +
* [http://www.eng.ox.ac.uk/~labap/flame/simsimplegeometry.m simsimplegeometry.m] contains all parameters and executes the simulation.
 +
* [http://www.eng.ox.ac.uk/~labap/flame/sfunsimplegeometry.m sfunsimplegeometry.m] contains the S-function (for more details on how to write S-functions in MALTAB see MATLAB's user guide). There is a line in this code for switching between time-varying and fixed delays.
 +
* [http://www.eng.ox.ac.uk/~labap/flame/mdlsimplegeometry.mdl mdlsimplegeometry.mdl] is the SIMULINK file.
 +
 
 +
For analyzing the results, and to find the mean over time, [http://www.eng.ox.ac.uk/~labap/flame/meanovtime.m meanovtime.m] is used.
 +
 
 +
Using these files, the response for omega = 160rad/s for the cases of fixed and time-varying delays is:
  
 
[[Image:Omega150fixeddelay.jpg|frame|left|Fixed delay, omega = 160 rad/s.]][[Image:Omega150tvdelay.jpg|frame|right|Time-Varying delay, omega = 160 rad/s.]]
 
[[Image:Omega150fixeddelay.jpg|frame|left|Fixed delay, omega = 160 rad/s.]][[Image:Omega150tvdelay.jpg|frame|right|Time-Varying delay, omega = 160 rad/s.]]
  
 
==More Complicated Geometry==
 
==More Complicated Geometry==
 +
 +
Consider now the more complicated geometry shown below.
 +
 
[[Image:Figure2.JPG]]
 
[[Image:Figure2.JPG]]
 +
 +
For this case, we have the following MATLAB files.
 +
* [http://www.eng.ox.ac.uk/~labap/flame/simcomplexgeometry.m simcomplexgeometry.m] contains all parameters and executes the simulation. This saves everything in a .mat file which can be called by showsimresult.m (see below).
 +
* [http://www.eng.ox.ac.uk/~labap/flame/sfuncomplexgeometry.m sfuncomplexgeometry.m] contains the S-function (for more details on how to write S-functions in MALTAB see MATLAB's user guide). There are several changes that can be made, e.g., one can allow for different equivalence ratios at different radial positions at the entry point, etc. 
 +
* [http://www.eng.ox.ac.uk/~labap/flame/mdlcomplexgeometry.mdl mdlcomplexgeometry.mdl] is the SIMULINK file.
 +
* [http://www.eng.ox.ac.uk/~labap/flame/showsimresult.m showsimresult.m] takes the mat file that was produced in simcomplexgeometry and animates the flame.
 +
 +
Apart from the [http://www.eng.ox.ac.uk/~labap/flame/meanovtime.m meanovtime.m] file, the newer version of MATLAB does not have a file I used in the above simulations:
 +
* [http://www.eng.ox.ac.uk/~labap/flame/vec.m vec.m] reshapes a matrix into a vector by stacking its columns.

Latest revision as of 17:05, 30 August 2008

Simple Geometry

Consider the geometry shown in the figure below.

Figure1.JPG

Three files can be used to simulate the response of this system.

  • simsimplegeometry.m contains all parameters and executes the simulation.
  • sfunsimplegeometry.m contains the S-function (for more details on how to write S-functions in MALTAB see MATLAB's user guide). There is a line in this code for switching between time-varying and fixed delays.
  • mdlsimplegeometry.mdl is the SIMULINK file.

For analyzing the results, and to find the mean over time, meanovtime.m is used.

Using these files, the response for omega = 160rad/s for the cases of fixed and time-varying delays is:

Fixed delay, omega = 160 rad/s.
Time-Varying delay, omega = 160 rad/s.

More Complicated Geometry

Consider now the more complicated geometry shown below.

Figure2.JPG

For this case, we have the following MATLAB files.

  • simcomplexgeometry.m contains all parameters and executes the simulation. This saves everything in a .mat file which can be called by showsimresult.m (see below).
  • sfuncomplexgeometry.m contains the S-function (for more details on how to write S-functions in MALTAB see MATLAB's user guide). There are several changes that can be made, e.g., one can allow for different equivalence ratios at different radial positions at the entry point, etc.
  • mdlcomplexgeometry.mdl is the SIMULINK file.
  • showsimresult.m takes the mat file that was produced in simcomplexgeometry and animates the flame.

Apart from the meanovtime.m file, the newer version of MATLAB does not have a file I used in the above simulations:

  • vec.m reshapes a matrix into a vector by stacking its columns.