CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Bugs

chtMultiRegionFoam: crash with fvDOM and unstructured mesh

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

Like Tree1Likes
  • 1 Post By Gidra

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 23, 2017, 13:56
Unhappy chtMultiRegionFoam: crash with fvDOM and unstructured mesh
  #1
New Member
 
Andrey P.
Join Date: Sep 2016
Location: Russian Federation
Posts: 3
Rep Power: 8
Gidra is on a distinguished road
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:rintStack(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)
Attached Files
File Type: gz test_structured.tar.gz (46.3 KB, 14 views)
File Type: gz test_unstructured.tar.gz (100.8 KB, 9 views)
Gidra is offline   Reply With Quote

Old   April 25, 2017, 11:53
Default
  #2
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,969
Blog Entries: 45
Rep Power: 126
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Quick questions:
  1. Which OpenFOAM version and/or fork are you using?
  2. Are you running in serial mode or parallel mode?
__________________
wyldckat is offline   Reply With Quote

Old   April 26, 2017, 12:58
Default
  #3
New Member
 
Andrey P.
Join Date: Sep 2016
Location: Russian Federation
Posts: 3
Rep Power: 8
Gidra is on a distinguished road
1. OpenFOAM-4.1 (www.openfoam.org) and OpenFOAM-v1612+ (www.openfoam.com)
2. serial mode
Gidra is offline   Reply With Quote

Old   May 6, 2017, 10:15
Smile
  #4
New Member
 
Andrey P.
Join Date: Sep 2016
Location: Russian Federation
Posts: 3
Rep Power: 8
Gidra is on a distinguished road
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;
}
wyldckat likes this.
Gidra is offline   Reply With Quote

Reply

Tags
fvdom, radiation, unstructured mesh

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On



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