CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Gas dispersion in air error message

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 28, 2013, 08:53
Default
  #21
New Member
 
Vivek
Join Date: Apr 2013
Posts: 9
Rep Power: 5
grk34 is on a distinguished road
Thanks Bruno!!

will do it! and i am using Openfoam version 2.2.0.
grk34 is offline   Reply With Quote

Old   April 28, 2013, 10:25
Default
  #22
New Member
 
Vivek
Join Date: Apr 2013
Posts: 9
Rep Power: 5
grk34 is on a distinguished road
Hi Brunos,

When i try to modify the reactingFoam model with my mesh i get the following error.

[COLOR="Red"]
Create time

Create mesh for time = 0


Reading g
Creating reaction model

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


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

Selecting chemistryReader foamChemistryReader
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 Uninterpreted:
#3 Foam::fvPatchField<double>:perator/=(Foam::fvPatchField<double> const&) at ??:?
...
.
.
.
.
.
.
at ??:?
#26 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#27
at ??:?
Floating point exception (core dumped)

Where is the problem??? What should i change?
grk34 is offline   Reply With Quote

Old   May 1, 2013, 23:48
Default
  #23
New Member
 
Vivek
Join Date: Apr 2013
Posts: 9
Rep Power: 5
grk34 is on a distinguished road
Quote:
Originally Posted by redhoax View Post
Hi. I'm running a reactingFoam case. kOmegaSST for turbulence, CH4 dispersion in air. no chemistry, no combustion. I get an error message that I'm not capable of interpreting, and was hoping maybe one of you guys were able to, or even if you could tip me on how to learn it. Anyways, here's the error:

Code:
#0  Foam::error::printStack(Foam::Ostream&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#6  Foam::compressible::RASModels::kOmegaSST::F2() const in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#7  Foam::compressible::RASModels::kOmegaSST::kOmegaSST(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&, Foam::word const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#8  Foam::compressible::RASModel::adddictionaryConstructorToTable<Foam::compressible::RASModels::kOmegaSST>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#9  Foam::compressible::RASModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#10  Foam::compressible::turbulenceModel::addturbulenceModelConstructorToTable<Foam::compressible::RASModel>::NewturbulenceModel(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#11  Foam::compressible::turbulenceModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libcompressibleTurbulenceModel.so"
#12  main in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/reactingFoam"
#13  __libc_start_main in "/lib/libc.so.6"
#14  Foam::UOPstream::write(char) in "/home/user/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/reactingFoam"
Thank you in advance.
Hi redhoax,
to understand the problem I tried my geometry in rhoSimpleFOAM. when i encountered this i changed my underrelaxation factors of rho and p. the simulation was running fine.

But when i run the case in reactingFOAM the problem is there. I am wondering is there any way i can change parameters ?(as i did for rhoSimpleFOAM) so i can run my simulations.
grk34 is offline   Reply With Quote

Old   May 5, 2013, 12:11
Default
  #24
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,655
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Lasse,
Quote:
Originally Posted by redhoax View Post
I updated the case and shared it with you on Dropbox. At the moment I'm experimenting with different variables, such as wind inlet velocity and fuel rate, and trying to decide how to make the platform geometry.

Lasse
I have finally managed to look into the folder you shared with me! And I noticed that your geometry has the "outlet" patch embedded into the "frontAndBack" patch! Therefore, the case is not properly defined, since the "outlet" patch is not present in the mesh

By the way, you mentioned in the share message, that when you activate the gravity, that the solver/solution would diverge. There are two possibilities:
  1. The missing "outlet" patch is critical. Without it, it's a wonder that the solver didn't crash sooner!
  2. The initial internal fields for "epsilon" and "k" and very low. Try starting with higher values, because (if I'm not mistaken) you're dealing with a compressible fluid solver and the pressure is absolute, not relative. Therefore, k-Epsilon values can be rather high and a good initial field is critical for the solution to not diverge!
Best regards,
Bruno
__________________
wyldckat is online now   Reply With Quote

Old   May 5, 2013, 18:32
Default
  #25
New Member
 
Join Date: Feb 2013
Posts: 17
Rep Power: 5
redhoax is on a distinguished road
wow, i'm so sorry. i discovered that myself the same day(just after) i updated the dropbox case and corrected it. mustve forgot to upload the corrected one. I will however test higher k-Epsilon values. i'll update that blockmeshdict asap
redhoax is offline   Reply With Quote

Old   May 6, 2013, 20:54
Default
  #26
New Member
 
Join Date: Feb 2013
Posts: 17
Rep Power: 5
redhoax is on a distinguished road
Hello!:-)

i updated the case. the solution is not diverging now, however...

i'm struggling with what seems to me as reversed flow from the outlet patch. at both top and bottom corner along the outlet(although one is located on the top patch) patch are strong velocity fields, oppositely directed. at the top is inflow, bottom outflow. see attached pictures for visual on how ridiculous my situation is. All taken at t=1s:
This is full mesh 3D: Ufields5.jpg
This is top and outlet only. The blue field on the top patch is -23m/s in negative X-direction...topandoutlet_UX_may7.png
Finally a picture of my outlet patch. How did it come to this? outlet_UX_May7.jpg

sine the issue arised after extracting the outlet patch from the top patch, and the problem shows on top and outlet patch, these boundary conditions are what i've looked into and played with the most. i made a table with an overview of my current conditions in the 0-folder. bcandparameters_table_may7.jpg (This image is in the dropbox folder in full size. Also the full case is updated.)

with this new issue i feel like i'm stumbling around in the dark looking for the solution in tutorials while not able to rationalize where it comes from, but i continue to look into it of course! And thanks a lot for all help you have provided me already. It's really generous


Regards,
Lasse
redhoax is offline   Reply With Quote

Old   May 7, 2013, 04:08
Default
  #27
New Member
 
Join Date: Feb 2013
Posts: 17
Rep Power: 5
redhoax is on a distinguished road
Update: i've tried to eliminate possible sources of error, but in failing to locate the root cause i have now made the case a laminar one with a viscosity a 1000 times higher than normal. still, after just t=1s the entire domain is showing unrealistic velocities..

Full domain from side: laminar_full.jpg

Air inlet only: airpatch_laminar.jpg

So the question is, what is causing the domain to obtain these velocity fields.. I'm in pursuit. Dropbox updated accordingly.

Lasse
redhoax is offline   Reply With Quote

Old   October 25, 2013, 05:31
Default Gas dispersion
  #28
New Member
 
Alex Lee
Join Date: Sep 2012
Posts: 16
Rep Power: 6
Alex Lee is on a distinguished road
Quote:
Originally Posted by redhoax View Post
I updated the case and shared it with you on Dropbox. At the moment I'm experimenting with different variables, such as wind inlet velocity and fuel rate, and trying to decide how to make the platform geometry.

Lasse

Can you share the case files with me too..... my email is alexleekh@singnet.com.sg Thanks
Alex Lee is offline   Reply With Quote

Reply

Tags
ch4, error, gas dispersion, komegasst, reactingfoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
error message cuteapathy CFX 14 March 20, 2012 07:45
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 07:59
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02


All times are GMT -4. The time now is 16:28.