CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (https://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   chtMultiRegionFoam: crash with fvDOM and unstructured mesh (https://www.cfd-online.com/Forums/openfoam-bugs/186603-chtmultiregionfoam-crash-fvdom-unstructured-mesh.html)

Gidra April 23, 2017 12:56

chtMultiRegionFoam: crash with fvDOM and unstructured mesh
 
2 Attachment(s)
Hello everyone,

It seems like a bug in fvDOM in chtMultiRegionFoam solver. The Simulation is running without problems with a structured mesh. However I got the following errors when I tried to use a unstructured mesh:

Quote:

Solving for fluid region air
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.99990091, Final residual = 8.102764e-10, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.99989569, Final residual = 6.1266637e-10, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.99999999, Final residual = 3.4002234e-13, No Iterations 1
Radiation solver iter: 0
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#7 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#9 Foam::radiation::radiativeIntensityRay::correct() at ??:?
#10 Foam::radiation::fvDOM::calculate() at ??:?
#11 Foam::radiation::radiationModel::correct() at ??:?
#12 ? at ??:?
#13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14 ? at ??:?
Floating point exception (core dumped)

wyldckat April 25, 2017 10:53

Quick questions:
  1. Which OpenFOAM version and/or fork are you using?
  2. Are you running in serial mode or parallel mode?

Gidra April 26, 2017 11:58

1. OpenFOAM-4.1 (www.openfoam.org) and OpenFOAM-v1612+ (www.openfoam.com)
2. serial mode

Gidra May 6, 2017 09:15

I think that I have solved the problem

1. change air fvScheme as: div(Ji,Ii_h) Gauss upwind;

2. change constant/air/radiationPropertie as:

absorptionEmissionModel constantAbsorptionEmission;

constantAbsorptionEmissionCoeffs
{
absorptivity absorptivity [ m^-1 ] 0.01;
emissivity emissivity [ m^-1 ] 0.01;
E E [ kg m^-1 s^-3 ] 0;
}


All times are GMT -4. The time now is 05:06.