在边界处输入负压,能够在场内传播,但直接对场内进行初始化赋值一个负压,就报错。报错信息如下:
[0] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
[0] #3 ? in /lib/x86_64-linux-gnu/libm.so.6
[0] #4 pow in /lib/x86_64-linux-gnu/libm.so.6
[0] #5 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::adiabaticPerfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::init(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) at ??:?
[0] #6 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::adiabaticPerfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
[0] #7 Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::adiabaticPerfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
[0] #8 Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
[0] #9 Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
[0] #10 Foam::phaseModel::phaseModel(Foam::word const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[0] #11 void Foam::LPtrList<Foam::DLListBase, Foam::phaseModel>::readIstream<Foam::phaseModel::iNew>(Foam::Istream&, Foam::phaseModel::iNew const&) at ??:?
[0] #12 Foam::DictionaryBase<Foam::LPtrList<Foam::DLListBase, Foam::phaseModel>, Foam::phaseModel>::DictionaryBase<Foam::phaseModel::iNew>(Foam::Istream&, Foam::phaseModel::iNew const&) at ??:?
[0] #13 Foam::multiphaseMixtureThermo::multiphaseMixtureThermo(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
[0] #14 ? at ??:?
[0] #15 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[0] #16 ? at ??:?
[hengzhu-PC:101729] *** Process received signal ***
[hengzhu-PC:101729] Signal: Floating point exception (8)
[hengzhu-PC:101729] Signal code: (-6)
[hengzhu-PC:101729] Failing at address: 0x3e800018d61
[hengzhu-PC:101729] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x12730)[0x7fc1e4967730]
[hengzhu-PC:101729] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0x10b)[0x7fc1e49675cb]
[hengzhu-PC:101729] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x12730)[0x7fc1e4967730]
[hengzhu-PC:101729] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x78856)[0x7fc1e4a0a856]
[hengzhu-PC:101729] [ 4] /lib/x86_64-linux-gnu/libm.so.6(pow+0x14)[0x7fc1e49a2004]
[hengzhu-PC:101729] [ 5] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam8heThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_21adiabaticPerfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE4initERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEESL_RSJ_+0xea)[0x7fc1e763a57a]
[hengzhu-PC:101729] [ 6] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam8heThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_21adiabaticPerfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEEC2ERKNS_6fvMeshERKNS_4wordE+0x2b2)[0x7fc1e766c342]
[hengzhu-PC:101729] [ 7] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo27addfvMeshConstructorToTableINS_11heRhoThermoIS0_NS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_21adiabaticPerfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEEEE3NewERKNS_6fvMeshERKNS_4wordE+0x2f)[0x7fc1e766c4bf]
[hengzhu-PC:101729] [ 8] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11basicThermo3NewINS_9rhoThermoEEENS_7autoPtrIT_EERKNS_6fvMeshERKNS_4wordE+0x11d)[0x7fc1e755c39d]
[hengzhu-PC:101729] [ 9] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo3NewERKNS_6fvMeshERKNS_4wordE+0x9)[0x7fc1e7556b69]
[hengzhu-PC:101729] [10] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libmultiphaseMixtureThermo.so(_ZN4Foam10phaseModelC1ERKNS_4wordERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES9_+0x402)[0x7fc1e7853132]
[hengzhu-PC:101729] [11] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libmultiphaseMixtureThermo.so(_ZN4Foam8LPtrListINS_10DLListBaseENS_10phaseModelEE11readIstreamINS2_4iNewEEEvRNS_7IstreamERKT_+0x13b)[0x7fc1e7884ffb]
[hengzhu-PC:101729] [12] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libmultiphaseMixtureThermo.so(_ZN4Foam14DictionaryBaseINS_8LPtrListINS_10DLListBaseENS_10phaseModelEEES3_EC2INS3_4iNewEEERNS_7IstreamERKT_+0x24)[0x7fc1e788f934]
[hengzhu-PC:101729] [13] /home/hengzhu/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/lib/libmultiphaseMixtureThermo.so(_ZN4Foam23multiphaseMixtureThermoC1ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEE+0xd3)[0x7fc1e7872ac3]
[hengzhu-PC:101729] [14] MultiphaseCavBubbleFoam[0x44ac1b]
[hengzhu-PC:101729] [15] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xeb)[0x7fc1e47b809b]
[hengzhu-PC:101729] [16] MultiphaseCavBubbleFoam[0x45423a]
[hengzhu-PC:101729] *** End of error message ***
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node hengzhu-PC exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
场内压力对应一个sine函数,速度由$U=p/(\rho_0*c_0)$获得。
边界条件如下:outside指两个自由面,wall 指压力输入面
outside
{
type waveTransmissive;
field p_rgh;
phi phi;
rho rho;
fieldInf 101325;
gamma 7.15;
lInf 10e-3;
value uniform 101325;
correctSupercritical false;
}
wall
{
type uniformFixedValue;
uniformValue sine;
uniformValueCoeffs
{
frequency 3333;
amplitude 8; //A
scale 100000; //s
level 101325; //l
t0 150e-6;
}
}