|
[Sponsors] |
January 26, 2012, 04:42 |
Problem with Surface Interpolation
|
#1 |
New Member
Giuseppe La Spina
Join Date: Jan 2012
Posts: 1
Rep Power: 0 |
Hello,
I'm a new user of OpenFOAM (2.0.1) and I started writing a new solver for two phase compressible model. The numerical scheme that I used is the one proposed by Kurganov and Tadmor, that is implemented in rhoCentralFoam but for a single phase. The problem that I found is related to the surface interpolation of the pressure. When I run the solver and it reaches the instruction Code:
surfaceScalarField p1_pos ( fvc::interpolate(p1, pos, "reconstruct(p)") ); Code:
#0 Foam::error::printStack(Foam::Ostream&) in "/home/giuseppe/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/giuseppe/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::LimitedScheme<double, Foam::MinmodLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/giuseppe/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #4 Foam::limitedSurfaceInterpolationScheme<double>::weights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/giuseppe/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #5 Foam::surfaceInterpolationScheme<double>::interpolate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/giuseppe/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #6 at rhoCentralFoam2.C:0 #7 in "/home/giuseppe/OpenFOAM/giuseppe-2.0.1/platforms/linux64GccDPOpt/bin/rhoCentralFoam2" #8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #9 in "/home/giuseppe/OpenFOAM/giuseppe-2.0.1/platforms/linux64GccDPOpt/bin/rhoCentralFoam2" Floating point exception So I tried to use different values of pressure, like 10^8 in the left section and 10^6 in the right section. With this choice of pressures I don't get any error... Thus I tried another way to compute the interpolation, using the following code Code:
surfaceScalarField p1_pos ( fvc::interpolate(0.001*p1, pos, "reconstruct(p)") ); surfaceScalarField p1_neg ( fvc::interpolate(0.001*p1, neg, "reconstruct(p)") ); p1_pos = p1_pos / 0.001; p1_neg = p1_neg / 0.001; Is it normal that I get this kind of errors when I use high values of pressure? Does exist a different way to avoid them? Thank you very much. Giuseppe Last edited by glaspina; January 31, 2012 at 03:48. |
|
Tags |
interpolate, pos, rhocentralfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem generating a surface mesh | danieldiver | STAR-CCM+ | 5 | August 8, 2011 08:34 |
urgent help needed (rhie-chow interpolation problem) | Ardalan | Main CFD Forum | 2 | March 18, 2011 15:22 |
Problem about a 2D surface rotation | skyblue_mech | Main CFD Forum | 0 | June 8, 2010 04:49 |
Problem with capturing water-spreading for free surface flow | devesh.baghel | OpenFOAM | 2 | December 10, 2009 01:21 |
Zero Volume Fraction in Free Surface Problem | marega | CFX | 1 | September 10, 2009 06:31 |