CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Error solving Energy Eqs - Compressible Solver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 11, 2017, 21:21
Default Error solving Energy Eqs - Compressible Solver
  #1
Member
 
fouad abi
Join Date: Nov 2015
Location: Geneva, Switzerland
Posts: 38
Rep Power: 3
faab is on a distinguished road
Send a message via Skype™ to faab
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)
faab is offline   Reply With Quote

Old   February 12, 2017, 09:58
Default
  #2
Member
 
fouad abi
Join Date: Nov 2015
Location: Geneva, Switzerland
Posts: 38
Rep Power: 3
faab is on a distinguished road
Send a message via Skype™ to faab
Hi,

After I initialized "kEff" to a default value I was able to get past the first error I described in the previous post. However, now I get the following error:

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2
at sigaction.c:?
#3 memcpy at interp.c:?
#4 std::string::append(std::string const&) in "/home/faabid/OpenFOAM/ThirdParty-2.3.0/platforms/linux64/gcc-4.8.1/lib64/libstdc++.so.6"
#5 std::basic_string<char, std::char_traits<char>, std::allocator<char> > std:perator+<char, std::char_traits<char>, std::allocator<char> >(std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) at ~/OpenFOAM/ThirdParty-2.3.0/platforms/linux64/gcc-4.8.1/include/c++/4.8.1/bits/basic_string.h:2370
#6
at ~/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvcDiv.C:240 (discriminator 2)
#7
at ~/OpenFOAM/OpenFOAM-2.3.0/applications/solvers/heatTransfer/CompressibleHeSolver/./fluid/EEqn.H:24 (discriminator 1)
#8 __libc_start_main at ??:?
#9
at ??:?
Segmentation fault (core dumped)

Where fluid/EEqn.H:24 corresponds to:
+ fvc::ddt(rho, K) + fvc::div(phi, K)
I also remind that K is set to: const volScalarField& K = 0.5*magSqr(U);

Could anyone please help me ??
faab is offline   Reply With Quote

Reply

Tags
compressible flow, energy equation, multi region

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Density-based, fully-coupled, compressible solver in OpenFOAM? francesco_capuano OpenFOAM Running, Solving & CFD 5 July 11, 2017 18:35
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 00:01
[ANSYS Meshing] Help with element size sandri_92 ANSYS Meshing & Geometry 13 May 26, 2014 09:05
Near compressible fluid - denstiy-based solver Hooman Main CFD Forum 8 December 19, 2011 04:03
Using pressure based solver for compressible flows aamer Main CFD Forum 2 June 15, 2011 12:45


All times are GMT -4. The time now is 21:14.