[gkhokhar]

I'm trying to run a case where I'm modeling combustion in a RAMJET.

There are options for wall heat flux, but I'd also like to use heat generation per unit volume.

This term is often included in the differential form of the energy equation, does SU2 allow a method for specifying this term?

Any help is appreciated.

[bigfootedrockmidget]

We do not have exactly what you are asking for, although addition of general source terms has many applications.

Something like this has been implemented for the radiation model in incompressible flow. There is the config option HEAT_SOURCE, as part of the radiation model:

% Apply a volumetric heat source as a source term (NO, YES) in the form of an ellipsoid (YES, NO)
HEAT_SOURCE = NO
% Value of the volumetric heat source
HEAT_SOURCE_VAL = 1.0E6
% Rotation of the volumetric heat source respect to Z axis (degrees)
HEAT_SOURCE_ROTATION_Z = 0.0
% Position of heat source center (Heat_Source_Center_X, Heat_Source_Center_Y, Heat_Source_Center_Z)
HEAT_SOURCE_CENTER = ( 0.0, 0.0, 0.0 )
% Vector of heat source radii (Heat_Source_Radius_A, Heat_Source_Radius_B, Heat_Source_Radius_C)
HEAT_SOURCE_AXES = ( 1.0, 1.0, 1.0 )

But if you would like to use something like this in general to create a heat source for the energy equation, you need to add it to the compressible solver as well, and make its implementation independent of the radiation solver.


This method is also available through the python wrapper, if you would like to have many different heat sources available, it is probably best to implement it through the python wrapper.

[bigfootedrockmidget]

You could also use the NEMO solver for it, then you have detailed chemistry available in the compressible solver, but you can not separate the energy source from the chemistry. So it will be much more expensive.

[gkhokhar]

Thanks. I think I'll try the NEMO solver. It would be great to add a radiative heat source, but since that's not possible for the compressible solver without implementing things I'm not sure how to do, I'll go with NEMO despite the added computational cost.

[bigfootedrockmidget]

Did you manage to get NEMO running for your testcase? Which fuel are you using?

[AeroHank]

Originally Posted by bigfootedrockmidget

Did you manage to get NEMO running for your testcase? Which fuel are you using?

We have not tried it yet. Need to get the implicit solver running a bit better first. But I believe our first go will be a two-step global mechanism with HTTB: https://doi.org/10.2514/6.2023-0161

[bigfootedrockmidget]

Is your case supersonic? We have FGM available for the incompressible solver. When you have the mechanism, it is relatively easy to add it, it's mainly these ionisation properties that are difficult to obtain, but you can just put some rough estimates in the files.