[virothi]

The jet is pressed through a cylindrical nozzle of a length of 70 mm and an
inner diameter of 1 mm. Afterward, it expands in initially quiescent air of atmo-
spheric pressure. The jet consists of a Helium/Nitrogen mixture of a density of 0.31 kg/m3. The velocity at the nozzle outlet is approximately 200 m/s. Both
the air and the jet are of a temperature of T = 293 K. How we can determine the velocity and mass fraction fields?

[LuckyTran]

You just need to use the species transport model. Other than that, nothing fancy here and just about any solver settings will work.

[virothi]

I appreciate your reply LuckyTran, but I'm sorry for not mentioning OpenFoam in my question. Could you please suggest a solver for that solving that situation and how can determine the velocity and mass fractions?

[LuckyTran]

Use a solver with a YEqn.H. Try reactingFOAM


The velocity is the U field and mass fraction is Y field. You don't have to do any advanced post-processing to get them. If you can setup the case properly in OpenFOAM, they are right there.

[virothi]

I'm not sure what went wrong, can some one help....
Thanks in advance

Selecting combustion model laminar<psiChemistryCombustion>
Selecting chemistry type
{
    chemistrySolver EulerImplicit;
    chemistryThermo psi;
}

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

Selecting chemistryReader foamChemistryReader
--> FOAM Warning :
    From function entry::getKeyword(keyType&, Istream&)
    in file db/dictionary/entry/entryIO.C at line 80
    Reading C:/blueCFD-Core-2016/msys64/home/ofuser/jetspray/constant/reactions
    found on line 11 the punctuation token '{'
    expected either } or EOF


--> FOAM FATAL ERROR:
HE not found in table.  Valid entries:
5
(
N2
CO2
O2
CH4
H2O
)


    From function const T& Foam::HashTable<T, Key, Hash>::operator[](const Key&) const [with T = Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >*; Key = Foam::word; Hash = Foam::string::hash]
    in file X:/OpenFOAM-4.x/src/OpenFOAM/containers/HashTables/HashTable/HashTableI.T.H at line 134.

FOAM exiting

[wyldckat]

Quick answer: Without being able to see the file "jetspray/constant/reactions", it's not possible for us to properly diagnose why this happened:

found on line 11 the punctuation token '{'
    expected either } or EOF

The best we can tell you is exactly what the error message is telling you: You opened a bracket "{" but never closed it with another "}" somewhere else later.

Please follow the instructions given here, in order to make it easier for people to help you: How to give enough info to get help

[virothi]

Hi Bruno, I've corrected it but now the simulation crashed. I don't know why. Could you please suggest what changes do I need to make.

The jet is pressed through a cylindrical nozzle of a length of 70 mm and an
inner diameter of 1 mm. Afterward, it expands in initially quiescent air of atmospheric pressure. The jet consists of a Helium/Nitrogen mixture of a density of 0.31 kg/m3. The velocity at the nozzle outlet is approximately 200 m/s. Both the air and the jet are of a temperature of T = 293 K. Discuss the resulting velocity and mass fraction fields.

I'm trying to solve it with reactingFoam by making the chemistry off and combustion as well. but the solution is getting crashed.
thermophysicalProperties.cchemistryProperties.txt

combustionProperties.txt

thermo.txt

reactions.txt

thermophysicalProperties.c

Originally Posted by wyldckat

Quick answer: Without being able to see the file "jetspray/constant/reactions", it's not possible for us to properly diagnose why this happened:

The best we can tell you is exactly what the error message is telling you: You opened a bracket "{" but never closed it with another "}" somewhere else later.

Please follow the instructions given here, in order to make it easier for people to help you: How to give enough info to get help

[raumpolizei]

Hey,


since you don't mention the error message, I assume you still have the same problem:

in your thermoPhysicalProperties.txt

foamChemistryThermoFile "$FOAM_CASE/constant/thermo.compressibleGas";

Looks like you are using the wrong file here. Could you change it to

foamChemistryThermoFile "$FOAM_CASE/constant/thermo";

(according to the name of the files you've posted)?

If the solver is crashing for another reason and you would like to get some help, you should consider posting at least the error message.

Cheers

RP

[virothi]

>thanks,, I did it, still getting error as shown

virothi@COOLTEK09 MINGW64 OpenFOAM-4.x ~/react
$ reactingFoam.exe
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*   Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt   *\
|  Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com   |
\*---------------------------------------------------------------------------*/
Build  : 4.x-ed69f631ce88
Exec   : C:/blueCFD-Core-2016/OpenFOAM-4.x/platforms/mingw_w64GccDPInt32Opt/bin/reactingFoam.exe
Date   : Dec 28 2018
Time   : 20:33:29
Host   : "COOLTEK09"
PID    : 6912
Case   : C:/blueCFD-Core-2016/msys64/home/ofuser/react
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Creating reaction model

Selecting combustion model noCombustion<psiThermoCombustion>
Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader foamChemistryReader


--> FOAM FATAL IO ERROR:
file "C:/blueCFD-Core-2016/msys64/home/ofuser/react/constant/Thermo" does not exist

file: C:/blueCFD-Core-2016/msys64/home/ofuser/react/constant/Thermo at line 1.

    From function Foam::IFstream& Foam::IFstream::operator()() const
    in file db/IOstreams/Fstreams/IFstream.C at line 188.

FOAM exiting

Originally Posted by raumpolizei

Hey,


since you don't mention the error message, I assume you still have the same problem:

in your thermoPhysicalProperties.txt

foamChemistryThermoFile "$FOAM_CASE/constant/thermo.compressibleGas";

Looks like you are using the wrong file here. Could you change it to

foamChemistryThermoFile "$FOAM_CASE/constant/thermo";

(according to the name of the files you've posted)?

If the solver is crashing for another reason and you would like to get some help, you should consider posting at least the error message.

Cheers

RP

[wyldckat]

Quick answers:

1. raumpolizei's diagnose is unfortunately incorrect.
2. @virothi: Please undo the change that raumpolizei suggested.
3. Please re-adapt the tutorial case "combustion/reactingFoam/laminar/counterFlowFlame2D" to use the same chemical reactions you want to use.
4. Then compress the whole case, for example if the case is named "counterFlowFlame2D.modified", with the following command:
   

   tar -czf counterFlowFlame2D.modified.tar.gz counterFlowFlame2D.modified

5. Then attach the file "counterFlowFlame2D.modified.tar.gz" to your next post.

My request here is because it's not clear what is the exact problem, because I'm not familiar with this solver.
The original error message indicated in post #5 is not clear on the origin of the problem, given that it's still listing an outdated list of chemical reactions. So it feels like there is more than one file in the "constant" folder... and this separate custom tutorial case should help divide the problems apart.