CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[swak4Foam] Aeroacoustic modelling using groovyBC

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By ngj

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 28, 2013, 08:27
Default Aeroacoustic modelling using groovyBC
  #1
New Member
 
Duncan Weatherhead
Join Date: Feb 2012
Location: University of Exeter
Posts: 11
Rep Power: 14
duncan21187 is on a distinguished road
Hi All
I am attempting to simulate pressure wave propagation in a pipe, using an oscillating pressure boundary condition defined using the groovyBC utility. When attempting to run I get the following error message:

Create time

Create mesh for time = 0


PIMPLE: max iterations = 50
field "(U|k|epsilon)" : relTol 0, tolerance 0.0001

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#6 Foam::compressible::RASModels::kOmegaSST::F2() const in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#7 Foam::compressible::RASModels::kOmegaSST::kOmegaSS T(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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#8 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kOmegaS ST>::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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#9 Foam::compressible::RASModel::New(Foam::GeometricF ield<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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#10 Foam::compressible::turbulenceModel::addturbulence ModelConstructorToTable<Foam::compressible::RASMod el>::NewturbulenceModel(Foam::GeometricField<doubl e, 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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#11 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleTurbulenceModel.so"
#12
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/rhoPimpleFoam"
#13 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#14
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/rhoPimpleFoam"
Floating point exception (core dumped)

I am aware that printStack errors often hint to a lack of sufficient memory to run the calculation (I am running this on a laptop) but I was wondering if there was anything else it could be?

Many thanks

Duncan
duncan21187 is offline   Reply With Quote

Old   January 28, 2013, 08:41
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Duncan,

If you read the error message from the bottom and upward, i.e. starting with "#14", it first tells you that somewhere in the turbulence model, a bad operation occurs (#11). Later on, it specifically tells you that it is in the compressible implementation of kOmegaSST (#7) and even more specifically in the F2() method of the turbulence model (#6).

The important thing then is told, namely that the operation, which goes wrong is a division (#4 and #3), which suggests division by 0 (read: zero). Thus conclusion:

Check that neither the internal field nor the boundary values in k and omega are not 0. I am unsure how the printStack would look like, but the zero could originate from a correct evaluation of groovyBC, which returns zero, which then goes wrong internally in the turbulence model.

Good luck,

Niels
immortality likes this.
ngj is offline   Reply With Quote

Old   January 29, 2013, 07:28
Default
  #3
New Member
 
Duncan Weatherhead
Join Date: Feb 2012
Location: University of Exeter
Posts: 11
Rep Power: 14
duncan21187 is on a distinguished road
Thanks Niels.
I have checked both of them and have improved the situation somewhat. It now appears to be a case of choosing the pressure boundary conditions. The simulation will crash after the first iteration with the following printStack error:
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::correct() in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/rhoPimpleFoam"
#6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#7
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/rhoPimpleFoam"
Floating point exception (core dumped)

if I use anything other than 'zeroGradient' for both inlet and outlet 'p' fields. I have tried using a 'groovyBCFixedValue' field and a simple 'fixedValue' field but both yield the same result. Can anyone suggest why this might be please?

Many thanks

Duncan
duncan21187 is offline   Reply With Quote

Old   January 29, 2013, 08:08
Default
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Duncan,

I have never run anything with the termodynamic models, so I am unfortunately not the right person to ask.

Though, you could go to the method

[CODE]
Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate()
[CODE]

and see, whether some of the operations could produce a crash (reason: latest reported point in the print stack). This could e.g. be sqrt of a negative number, negative numbers to a scalar power, tanh to a large number, etc.

Good luck,

Niels
ngj is offline   Reply With Quote

Reply

Tags
aeroacoustics, groovybc, printstack, rhopimplefoam


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[swak4Foam] groovyBC for oscillatory flow liybzd OpenFOAM Community Contributions 5 November 12, 2018 08:53
CFD in Naval Hydrodynamics, Off-Shore and Wave Modelling with OpenFOAM hjasak OpenFOAM Announcements from Other Sources 2 February 13, 2017 05:59
[swak4Foam] reactingMultiPhaseEulerFoam problems with groovyBC zanilu70 OpenFOAM Community Contributions 4 December 13, 2016 07:46
[swak4Foam] Change in alpha and U with groovyBC in twoPhaseEulerFoam dani2702 OpenFOAM Community Contributions 0 November 17, 2016 04:30
error message cuteapathy CFX 14 March 20, 2012 07:45


All times are GMT -4. The time now is 01:48.