Hello, everyone,
Could anyone tell me how can the variance of mixture fraction be normalized by this: Code:
Zvar_normalized = Zvar[celli] / (Z[celli]*(1.-Z[celli])); Thank you in advance! Best regards, Wen |
Its Not the mixture fraction, it's its Varianz ... Zvar is defined between 0 and 0.25 (parabolic profile). No area weight ed interpolation possible.
See also http://www.holzmann-cfd.de/index.php...ameletModel2.2 |
Thank you for your reply. And sorry for my mistake. I have corrected it.
Regards, Wen |
Dear Tobi,
May I kindly ask two questions about the libOpenSMOKE: (1) About the enthalpy defect: I do not know how you generate the flamelet library with different enthalpy defects, that means do you modify the boundary condition of the temperature flamelet equation or add a source term to the that equation? (2) About the actual enthalpy (H): How to calculate the actual enthalpy at the inlet as that in the H initial file? Quote:
Best, Wen |
Dear Wen,
(1) the enthalpy defect is a source term added in the flamelet equations (not in the openFOAM solver). The flamelet equations are not solved in openFOAM. These eqn are solved within the binary files and you can add a defect, that means a delta enthalpy to the adiabatic condition. (2) The actual enthalpy at some cell is calculated using a passive scalar transport equation (not what you mentioned). The equation is (in SIMPLE): Code:
fvScalarMatrix hEqn The adiabatic values of the inlet are the adiabatic enthalpies of fuel and oxidizer and is shown at the beginning of the libOpenSMOKE solver. Run it once, abort, get the values of the adiabatic enthalpy, insert them into the H file and restart. |
I do not know whether the enthalpy defect is set as a source term or as the boundary condition in the fuel and oxidizer sides (the common choice), because it is a binary file. But I know that the total enthalpy on the fuel and oxidizer side is a calculated value if we use FlameMaster solver. Of course, the method you implemented is right.
Thank you for you hints, I already know that. Best regards, Wen |
Dear Wenxu,
you ment the enthalpy defect within the flamelet equations (so flame-master, binaries of libOpenSMOKE, cantera)? I think it is a source term and not a boundary condition because if you compare flamelets of the same scalar dissipation rates and different enthalpy defects you can see that the temperature profile is moved (somehow with a defined value --> delta distribiution). Therefore I think its not a boundary condition |
Hello everyone,
i have compiled flamelet model in openfoam-4.x but while running the tutorial case i am getting the following error : Code:
--> FOAM FATAL ERROR: |
Hi Tobi, Hope your holidays were awesome,
I have compiled flamelet model in openfoam-4.x but while running the tutorial case i am getting the following error : --> FOAM FATAL ERROR: Not implemented From function flameletThermo::Z() in file flameletThermo/flameletThermo.C at line 160. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::flameletThermo::Z() at ??:? #3 ? at ??:? #4 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #5 ? at ??:? Aborted (core dumped) |
You do not have to post it everywhere - I already got the notification ;)
|
Hi Tobi,
I am using ur flameletPimplefoam in openfoam-2.3.x for supersonic combustion but after some iterations it is diverging. Z and Zvar values go in thousands. Can u help me out ? |
Maybe it is based on the fact that the solver has no supersonic implementation or it is based on your mesh, boundary conditions or whatever. In addition I do not support old versions anymore.
If it is not possible, try the flameletFoam from Hagen Müller. |
Density calculation
Hi Tobi,
First of all: thank you for all the work you have done on the flamelet solver! I have a question about the rho calculation in the thermophysical model. Am I correct in assuming that rho is looked up in the flamelet table? I am a bit confused due to the "flameletThermo.C" file in which rho is returned as p*psi() and psi is calculated as 1/RT in the perfect gas equation of state. But in pdfFlameletThermo.C "density_reynolds" is looked up and used in calculations for psi. However, I do not seem to be able to determine the place where the rho-field is connected to the looked up value. A related question: in the validation syngas flame you use the following values for the pure mixture: mixture { specie { nMoles 1; molWeight 28.9; } thermodynamics { Cp 1005; Hf 0; } transport { As 1.67212e-06; Ts 170.672; } } How did you determine the molWeight and Cp value? Or are they not used in the simulation at all? |
Hi,
Code:
|
Thank you for the clarification!
So that means that also none of the following parameters of thermotype matter, except "type pdfFlameletThermo;" thermoType { type pdfFlameletThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } |
Correct.
It is just needed for the object we construct but in the new thermodynamic we made in that model, everything is based on the Look-Up-Table:
|
Wall boundary condition
Hi Tobi,
I have a question about the implemented boundary conditions in the pdfFlameletThermo.C file. I don't get why the boundary condition for fixed temperature should coïncide with a fixed enthalpy BC with value 0 as the following error message suggests: FatalErrorIn ( "pdfFlameletThermo<BasicFlameletThermo, MixtureType>::update()" ) << "Boundary conditions are wrong: " << "fixed temperature BC must be fixed enthaplie BC with value 0;" << abort(FatalError); I think that a wall would have a zero gradient BC for H and a fixedValue one for T, because the enthalpy H is not necessarily 0 at a wall. It is an easy alteration (for me/you/...) to change/add this condition, but maybe I'm overseeing something. What do you think about this requirement? Thank you! |
Hi,
it was my implementation. The reason for that is very easy. It is based on the fact that libOpenSMOKE has some special treatment for fixedValue BC for the temperature. Based on the temperature the application calculates an enthalpydefect and thus choose the appropriate LUT. It was not working in the official one so I debugged a few things here. The enthalpy fixed value = 0 does not influence the solution. it is just a dummy for further treatment in the code. I have everything in my thesis. Unfortunately it is in german. Even the usage of the log-normal in the properties will change a lot of things :). To make it short: Fixed T: With T, Z, Z'', chi at face find enthalpy defect at face. Based on the enthalpy defect, calculate the enthalpy at the face. |
1 Attachment(s)
|
All times are GMT -4. The time now is 21:39. |