求指导compressibleInterfoam初始压力场设置

在边界处输入负压,能够在场内传播,但直接对场内进行初始化赋值一个负压,就报错。报错信息如下:

[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;
        }
    }