[niklas]

> Can you please tell me how and where can I turn on the turbulent dispersion?

have you looked in sprayProperties?

dispersionModel off;
or
dispersionModel stochasticDispersionRAS;
or
dispersionModel gradientDispersionRAS;

although I dont recommend using turbulent dispersion in 2D with non evaporation sprays.
It looks awful.

[markusrehm]

Hello,

if I use another Schmidt number say Sc_t=0.7 and put it like Niklas proposed

fvm::laplacian(turbulence->muEff()/Sc, Yi)

my temperature drops quickly to 0K.

Do I have to adjust something similiar in hEqn?

Regards Markus.

[haghajani]

Hi,

I am using "reactingFoam", to model high pressure carbon dioxide release in gas phase (near saturation line).

I added new Equation of state (peng Robinson). To make reactingfoam to work with Peng Robinson EOS, I have created a new chemistryReader with sutherlandTransport<specieThermo<janafThermo<pengR obinson>>> in typedef.

In consequence, I created new chemistry Model, new reacting mixture and new chemistrySolver to make openfoam read from pengRobinson Eos.

after changing the relevant libraries in reactingfoam and compiling, I tried to run a case with updated reactingfoam. I got this error;
...
Constructing chemical mechanism
terminate called after throwing an instance of 'std::bad_cast'
what(): std::bad_cast
Aborted

Would you please let me know where should i correct?
Is there any better approach to work with a new EoS in reactingFoam

Best regards,
Hamed

[gschaider]

Originally Posted by haghajani

Hi,

I am using "reactingFoam", to model high pressure carbon dioxide release in gas phase (near saturation line).

I added new Equation of state (peng Robinson). To make reactingfoam to work with Peng Robinson EOS, I have created a new chemistryReader with sutherlandTransport<specieThermo<janafThermo<pengR obinson>>> in typedef.

In consequence, I created new chemistry Model, new reacting mixture and new chemistrySolver to make openfoam read from pengRobinson Eos.

after changing the relevant libraries in reactingfoam and compiling, I tried to run a case with updated reactingfoam. I got this error;
...
Constructing chemical mechanism
terminate called after throwing an instance of 'std::bad_cast'
what(): std::bad_cast
Aborted

Would you please let me know where should i correct?
Is there any better approach to work with a new EoS in reactingFoam

Best regards,
Hamed

The source of your problem is the line in the chemistryModel-constructor where it says:


reactions_(dynamic_cast<const reactingMixture&>(thermo)),

To cut a long story short: reactingFoam only works with thermoPhysical-models that are subclasses of reactingMixture (and I think, but am not sure, that some stuff also assumes that your EoS is perfect gas).

So there is no trivial way to do what you want, I'm afraid

Bernhard

[haghajani]

Dear Bernhard,
You are right,

As far as I try new ways, i get few success!

As you mentioned, the source of problem is the line "reactions_(dynamic_cast<const nreactingMixture&>(thermo))", in my case.

Is it possible to develope thermophysicalModels/combustion/hCombustionThermo/hCombustionThermos.C ; where currently is,
defineTemplateTypeNameAndDebug(hMixtureThermo<reactingMixture>, 0)

typedef hMixtureThermo<reactingMixture> hMixtureThermoReactingMixture;

addToRunTimeSelectionTable
(
hCombustionThermo,
hMixtureThermoReactingMixture,
fvMesh
);
addToRunTimeSelectionTable
(
basicThermo,
hMixtureThermoReactingMixture,
fvMesh
);
for reactingFoam.

I have found the link http://www.cfd-online.com/Forums/ope...eal-gases.html, where prof.Jasak post a comment for already all solvers except reactingFoam.

Thanks a lot,
Hamed

[gschaider]

Originally Posted by haghajani

Dear Bernhard,
You are right,

As far as I try new ways, i get few success!

As you mentioned, the source of problem is the line "reactions_(dynamic_cast<const nreactingMixture&>(thermo))", in my case.

Is it possible to develope thermophysicalModels/combustion/hCombustionThermo/hCombustionThermos.C ; where currently is,
defineTemplateTypeNameAndDebug(hMixtureThermo<reactingMixture>, 0)

typedef hMixtureThermo<reactingMixture> hMixtureThermoReactingMixture;

addToRunTimeSelectionTable
(
hCombustionThermo,
hMixtureThermoReactingMixture,
fvMesh
);
addToRunTimeSelectionTable
(
basicThermo,
hMixtureThermoReactingMixture,
fvMesh
);
for reactingFoam.

I have found the link http://www.cfd-online.com/Forums/ope...eal-gases.html, where prof.Jasak post a comment for already all solvers except reactingFoam.

Thanks a lot,
Hamed

Hi Hamed!

Something happened between my last posting and your post: OpenfOAM 1.6 was released!

Have a look at section 4.3 of http://www.opencfd.co.uk/openfoam/do...Notes-1.6.html
This might be of interest for you

Bernhard

[haghajani]

Dear Bernhard,
Thanks a lot for the clue,

I have changed the Thread, because it is not related to std_cast error anymore, Would you please do me a favor and reply my question in the this Thread:
http://www.cfd-online.com/Forums/ope...eal-gases.html


Best regards,
Hamed

[vijay480]

Hello all,

I tried simulating sprayEngineFoam and I have the following errror. I am working in 2.2.0 version ..

Build : 2.2.0-5be49240882f
Exec : sprayEngineFoam
Date : May 28 2013
Time : 22:00:01
Host : "ubuntu"
PID : 2971
Case : /home/sriksm/OpenFOAM/sriksm-2.2.0/run/forte
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create engine time

Create mesh for time = -180

Selecting engineMesh layered
deckHeight: 0.11
piston position: 0.0179

Reading g
Creating combustion model

Selecting combustion model PaSR<psiChemistryCombustion>
Selecting chemistry type
{
chemistrySolver ode;
chemistryThermo psi;
}

Selecting thermodynamics package
{
type hePsiThermo;
mixture inhomogeneousMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

terminate called after throwing an instance of 'std::bad_cast'
what(): std::bad_cast
Aborted (core dumped)



Can anyone please help me ?