[danhnam]

Originally Posted by eliotfoss

Thanks for your reply, I see, I did not realize there would be a problem with the pimple loop itself, I hope that you can share you findings soon!

I had another question, in your paper, you define the strain rate for the counterflow non-premixed flame with a multiplying factor of 2. In the textbooks I've read this 2 is not included in the definition. Is this due to a lack of standardization for the definition, or some other reason?

Thank you for your time!

Hi Eliot Foss!

We referred to this paper for the global strain rate definition. https://www.researchgate.net/publica...ffusion_flames

Danh Nam,

[danhnam]

Originally Posted by rmishra

In reference to the following file:

src/thermophysicalModels/specie/thermo/rfJanaf/rfJanafThermoI.H (Lines 286-295)

The dcpdT term is ideal. Since we are dealing with real gas combustion we need to solve for real gas Jacobian. So the temperature derivative will also include the derivative of departure function apart from the ideal gas part.

Please let me know if you have time to discuss this.

Hi Mishra!

You are right. When you retrieve the temperature from the enthalpy for real-gas cases, you need to derive dH/dT term (i.e., derivative of departure part w.t. the temperature). The existing code is only applicable for ideal gas problems.

Danh Nam,

[ShinyRiver]

Originally Posted by onofrio

Good morning,

I found that, in PengRobinson equationOfState model, the M term, present in both Cp and CpMCv functions, is the same as in the soaveRedlichKwong model when actually it is not so (as you can see in the original PengRobinsonGas formulation). This fact leads to negative Cp (I’m dealing with nitrogen at cryogenic conditions).
If you substitute back the original formulation for the M term everything works fine.

Nicola

Hi Nicola,

I have the same problem of negative CP in FLUENT. There I have the possibility to correct the Peng-Robinson equation. Could you explain this modification a bit more in detail?

Thanks in advance

[rmishra]

Hello,

I believe one of the ways the chemistry libraries are changed, the chemistry model has to be selected among the different options of real gas EOS. However, can you tell me how would one go about also implementing reduction/tabulation since in classic OF the option was there as one of the entries to chemistryModel. Thanks.

[araphi]

Hello,
When I install, the directory ($WM_PROJECT_DIR/wmake/scripts/AllwmakeParseArguments) is missing, the library cannot be installed, blockMesh and other operations cannot be performed.

[adorean]

Hello,

I saw in the tutorials that one is able to generate properties file for a temperature range but for a specific pressure.
How this library can be used for a case where one has 100 bar pressure difference between inlet and outlet?

Thank you!

Ervin

[danhnam]

Originally Posted by rmishra

Hello,

I believe one of the ways the chemistry libraries are changed, the chemistry model has to be selected among the different options of real gas EOS. However, can you tell me how would one go about also implementing reduction/tabulation since in classic OF the option was there as one of the entries to chemistryModel. Thanks.

Hi rmishra!

Sorry for my late reply. It's quite complicated. If you are still want to implement it for reduction/tabulation, please contact me via my email (danhnam11@gmail.com).

Danh Nam

[danhnam]

Originally Posted by araphi

Hello,
When I install, the directory ($WM_PROJECT_DIR/wmake/scripts/AllwmakeParseArguments) is missing, the library cannot be installed, blockMesh and other operations cannot be performed.

Hello araphi,

Did you solve your problem?

[danhnam]

Originally Posted by adorean

Hello,

I saw in the tutorials that one is able to generate properties file for a temperature range but for a specific pressure.
How this library can be used for a case where one has 100 bar pressure difference between inlet and outlet?

Thank you!

Ervin

Hello Ervin!
Since the test case in the tutorial is designated only for generating thermophysical properties. Theoretically, if you provide p, T, then you can get properties. Actually, I don't know why you need to generate properties at such conditions (i.e., having a pressure difference between the inlet and outlet). If you need properties at 100 bar, just change the pressure in the internal domain to 100 bar. If you want to get the properties for a range of pressure, you can run several times corresponding to different pressures.
For your case, I have no idea. It would be better if you provide more information about your case. So that I may help you.

Danh Nam,

[adorean]

Hello Danh Nam,

If I understand correctly, for a multiphase single component non reacting (CO2) simulation, for example if at time 0 the domain/outlet has 10 bar pressure, one needs to generate properties for 10 bar. And even if the inlet pressure BC is 100 bar, the model is able to get itself the correct properties, or I need to provide multiple .dat files (10,20, ... 100 bar)?

For example for interPhaseChangeFoam how should I change the thermophysicalProperties file and how to use the generated tabular properties file OF.p10.CO2.dat?

Is there somewhere such an example case, for whatever solver?

Thank you!

Best regards,

Ervin

[adorean]

Hello,

Anyone, please?

In Ansys, one can generate 1 real gas properties table for a range of temperatures and pressures. How can this be achieved here?

And how can one point the OpenFOAM solvers to use these properties files (OF.p10.CH4.dat, for example) in a case?

Thank you!

Best regards,

Ervin

[danhnam]

Hi Ervin!

If I understand correctly, for a multiphase single component non reacting (CO2) simulation, for example if at time 0 the domain/outlet has 10 bar pressure, one needs to generate properties for 10 bar. And even if the inlet pressure BC is 100 bar, the model is able to get itself the correct properties, or I need to provide multiple .dat files (10,20, ... 100 bar)?

If you used the model in my library (and solver), it will automatically update your properties regardless your pressure inside the domain is 10 bar or 100 bar. As I sad, if you provide value of p and T, you can get properties.

For your case, you cannot use my real-gas library in interPhaseChangeFoam if you don't modify the solver since the interPhaseChangeFoam solver does not use "psiThermo" or "rhoThermo" class. And to develop a new version of "interPhaseChangeFoam" solver to be able used real-gas library is a different story.


Best regards,
Danh Nam

[ricx]

Hi!
Thank you for the library.
I am using a in house solver based on twoPhaseEulerFoam that requires Thermophysical models but not any kind of reactions. I am planning to add an equation of state for CH4 and also I would need the data for CH4 in order to calibrate the thermo coefficient to be chosen in a context of two phase liquid-gas non newtonian flow.
Would that make sense if I don`t consider reactions (and so also temperature effects)?

Thank you.

[danhnam]

Originally Posted by ricx

Hi!
Thank you for the library.
I am using a in house solver based on twoPhaseEulerFoam that requires Thermophysical models but not any kind of reactions. I am planning to add an equation of state for CH4 and also I would need the data for CH4 in order to calibrate the thermo coefficient to be chosen in a context of two phase liquid-gas non newtonian flow.
Would that make sense if I don`t consider reactions (and so also temperature effects)?

Thank you.

Hi!
I've just quickly checked the tutorial test case in "twoPhaseEulerFoam" and I found that this solver uses pureMixture model.
Actually, my real-gas library is not applicable for pureMixture model. If you want to implement it, pls contact me. I will let you know. For pureMixture its quite easy, you just copy my source code and do some modification.

Danh Nam.

[ricx]

Originally Posted by danhnam

Hi!
I've just quickly checked the tutorial test case in "twoPhaseEulerFoam" and I found that this solver uses pureMixture model.
Actually, my real-gas library is not applicable for pureMixture model. If you want to implement it, pls contact me. I will let you know. For pureMixture its quite easy, you just copy my source code and do some modification.

Danh Nam.

Hi,

Thank you for your answer. Yes, I am interested, I saw you shared your e-mail, I will contact you there.

Riccardo

[Congren_Zheng]

Hi, the library is very nice. I got a problem during Allwmake.
The error code is like this:
In file included from specie/specie.H:155:0,
from specie/specie.C:26:
specie/specieI.H: In member function ‘void Foam::specie:perator+=(const Foam::specie&)’:
specie/specieI.H:108:21: error: ‘small’ was not declared in this scope
if (mag(sumY) > small)
^~~~~
specie/specieI.H:108:21: note: suggested alternative: ‘scalbl’
if (mag(sumY) > small)
^~~~~
scalbl
specie/specieI.H: In function ‘Foam::specie Foam:perator+(const Foam::specie&, const Foam::specie&)’:
specie/specieI.H:127:46: error: ‘small’ was not declared in this scope
const scalar sumY = max(st1.Y_ + st2.Y_, small);
^~~~~
specie/specieI.H:127:46: note: suggested alternative: ‘scalbl’
const scalar sumY = max(st1.Y_ + st2.Y_, small);
^~~~~
scalbl
specie/specieI.H: In function ‘Foam::specie Foam:perator==(const Foam::specie&, const Foam::specie&)’:
specie/specieI.H:157:22: error: ‘small’ was not declared in this scope
if (mag(diffY) < small)
^~~~~
specie/specieI.H:157:22: note: suggested alternative: ‘scalbl’
if (mag(diffY) < small)
^~~~~
scalbl
specie/specieI.H:165:24: error: ‘great’ was not declared in this scope
scalar molWeight = great;
^~~~~
specie/specieI.H:165:24: note: suggested alternative: ‘read’
scalar molWeight = great;
^~~~~
read
specie/specieI.H:166:23: error: ‘small’ was not declared in this scope
if (mag(diffRW) > small)
^~~~~
specie/specieI.H:166:23: note: suggested alternative: ‘scalbl’
if (mag(diffRW) > small)
^~~~~
scalbl
/usr/lib/openfoam/openfoam2306/wmake/rules/General/transform:37: recipe for target 'Make/linux64GccDPInt32Opt/specie/specie.o' failed
make: *** [Make/linux64GccDPInt32Opt/specie/specie.o] Error 1

The problem is that there is no definition of variable 'small'

[shrigames]

I got the same error for Allwake. And is the real Gas solver configured for openFoam 9 and above?

[al.csc]

Thank you for sharing such a great work. Does anyone plan to port it to ESI? If not, I can start and create a pull request

[Li Xuejiao]

When I wanted to replace methane with hydrogen in the case tutorials/laminarCounterFlowFlame-2D/GRI30/, I used H2_Li from inputThermoData_OF_format, but I encountered the following error. Does anyone know the reason for this? I would greatly appreciate your help.
Selecting chemistryReader foamChemistryReader
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam:imensionedField<double, Foam::volMesh>:perator/=(Foam:imensionedField<double, Foam::volMesh> const&) at ??:?
#4 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>:perator/=(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::SRKchungKineticMixture<Foam::chungKineticTra nsport<Foam::species::thermo<Foam::rfJanafThermo<F oam::soaveRedlichKwongFoam::rfSpecie >, Foam::sensibleEnthalpy> > >::correctMassFractions() at ??:?
#6 Foam::SRKchungKineticMixture<Foam::chungKineticTra nsport<Foam::species::thermo<Foam::rfJanafThermo<F oam::soaveRedlichKwongFoam::rfSpecie >, Foam::sensibleEnthalpy> > >::SRKchungKineticMixture(Foam::dictionary const&, Foam::ListFoam::word const&, Foam::HashPtrTable<Foam::chungKineticTransport<Foa m::species::thermo<Foam::rfJanafThermo<Foam::soave RedlichKwongFoam::rfSpecie >, Foam::sensibleEnthalpy> >, Foam::word, Foam::string::hash> const&, Foam::fvMesh const&, Foam::word const&) at ??:?
#7 Foam::SRKchungKineticReactingMixture<Foam::chungKi neticTransport<Foam::species::thermo<Foam::rfJanaf Thermo<Foam::soaveRedlichKwongFoam::rfSpecie >, Foam::sensibleEnthalpy> > >::SRKchungKineticReactingMixture(Foam::dictiona ry const&, Foam::fvMesh const&, Foam::word const&) at ??:?
#8 Foam::heThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::SRKchungKineticReactingM ixture<Foam::chungKineticTransport<Foam::species:: thermo<Foam::rfJanafThermo<Foam::soaveRedlichKwong Foam::rfSpecie >, Foam::sensibleEnthalpy> > > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#9 Foam::hePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::SRKchungKineticReactingM ixture<Foam::chungKineticTransport<Foam::species:: thermo<Foam::rfJanafThermo<Foam::soaveRedlichKwong Foam::rfSpecie >, Foam::sensibleEnthalpy> > > > >::hePsiThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#10 Foam:siReactionThermo::addfvMeshConstructorToTab le<Foam::hePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::SRKchungKineticReactingM ixture<Foam::chungKineticTransport<Foam::species:: thermo<Foam::rfJanafThermo<Foam::soaveRedlichKwong Foam::rfSpecie >, Foam::sensibleEnthalpy> > > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
Floating point exception (core dumped)

[wellzt]

Hi guys,

Did somebody ever tried BC "externalWallHeatFluxTemperature"?
I encounter some problems with this, when I use BC like
{
type externalWallHeatFluxTemperature;
mode flux;
q uniform 118000;

kappaMethod lookup;
kappa kappa;

value $internalField;
}
Solver will crush at EEqn.solve():
{
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::rfJanafThermo<Foam::PengRobinson<Foam::rfSpe cie> >::Ha(double, double) const at ??:?
#4 Foam::heThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::PRchungTakaReactingMixtu re<Foam::chungTakaTransport<Foam::species::thermo< Foam::rfJanafThermo<Foam::PengRobinson<Foam::rfSpe cie> >, Foam::sensibleEnthalpy> > > > >::he(Foam::Field<double> const&, Foam::Field<double> const&, int) const at ??:?
#5 Foam::mixedEnergyFvPatchScalarField::updateCoeffs( ) at ??:?
#6 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/home/wellzt/OpenFOAM/wellzt-6/platforms/linux64GccDPInt32Opt/bin/realFluidRhoReactingBuoyantFoam"
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fv:ptionList:perator()<double>(Foam::Geo metricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/wellzt/OpenFOAM/wellzt-6/platforms/linux64GccDPInt32Opt/bin/realFluidRhoReactingBuoyantFoam"
#8 ? in "/home/wellzt/OpenFOAM/wellzt-6/platforms/linux64GccDPInt32Opt/bin/realFluidRhoReactingBuoyantFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? in "/home/wellzt/OpenFOAM/wellzt-6/platforms/linux64GccDPInt32Opt/bin/realFluidRhoReactingBuoyantFoam"
Floating point exception
}
Why should it happend? externalWallHeatFluxTemperatureFvPatchScalarField. C located at src\TurbulenceModels\compressible\turbulentFluidTh ermoModels\derivedFvPatchFields\externalWallHeatFl uxTemperature. As I see, the turbulentFluidThermoModels were recompiled in realFluidThermophysicalModels-6

best,
jeff