|
[Sponsors] |
Error solving Energy Eqs - Compressible Solver |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 11, 2017, 20:21 |
Error solving Energy Eqs - Compressible Solver
|
#1 |
Member
|
Dear FOAMers,
I am working on a compressible solver (density based) that I created using the pre-existing chtMultiRegionFoam solver. As stated in the title, I get (see hereunder) an error when solving for the following energy equation: fvScalarMatrix EEqn ( fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K) == fvm::laplacian(kEff, T) + HS ); HS is just a volumetric source term, and kEff is the effective conductivity. he is the internal energy (I compute it using an external Fortran code) while K is set as being: const volScalarField& K = 0.5 * (U & U); From the log I can tell that the error most probably arises from the laplacian term, and a division by 0 occurs somewhere, but after investigating I cannot find out the source of the error. Thanks a lot for your help ! #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 at sigaction.c:? #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam:perator/(Foam::UList<double> const&, Foam::tmp<Foam::Field<double> > const&) at ??:? #5 Foam::compressible::customturbulentTemperatureCoup ledBaffleMixedFvPatchScalarField::updateCoeffs() at ??:? #6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() at ??:? #7 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:? #8 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #9 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #10 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #11 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) at ~/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvmLaplacian.C:219 (discriminator 1) #12 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ~/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvmLaplacian.C:251 (discriminator 5) #13 at ~/OpenFOAM/OpenFOAM-2.3.0/applications/solvers/heatTransfer/CompressibleHeSolver/./fluid/EEqn.H:24 #14 __libc_start_main at ??:? #15 at ??:? Floating point exception (core dumped) |
|
Tags |
compressible flow, energy equation, multi region |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Density-based, fully-coupled, compressible solver in OpenFOAM? | francesco_capuano | OpenFOAM Running, Solving & CFD | 7 | November 20, 2018 13:12 |
[ANSYS Meshing] Help with element size | sandri_92 | ANSYS Meshing & Geometry | 14 | November 14, 2018 07:54 |
Star cd es-ice solver error | ernarasimman | STAR-CD | 2 | September 12, 2014 00:01 |
Near compressible fluid - denstiy-based solver | Hooman | Main CFD Forum | 8 | December 19, 2011 03:03 |
Using pressure based solver for compressible flows | aamer | Main CFD Forum | 2 | June 15, 2011 12:45 |