|
[Sponsors] |
[swak4Foam] Aeroacoustic modelling using groovyBC |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 28, 2013, 08:27 |
Aeroacoustic modelling using groovyBC
|
#1 |
New Member
Duncan Weatherhead
Join Date: Feb 2012
Location: University of Exeter
Posts: 11
Rep Power: 14 |
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 |
|
January 28, 2013, 08:41 |
|
#2 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37 |
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 |
|
January 29, 2013, 07:28 |
|
#3 |
New Member
Duncan Weatherhead
Join Date: Feb 2012
Location: University of Exeter
Posts: 11
Rep Power: 14 |
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 |
|
January 29, 2013, 08:08 |
|
#4 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37 |
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 |
|
Tags |
aeroacoustics, groovybc, printstack, rhopimplefoam |
|
|
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 |