search

duncanam / rhoReactingCentralFoam

An OpenFOAM solver with AMR that combines rhoCentralFoam and rhoReactingFoam for high-speed reactive flows.
GNU General Public License v3.0
47 stars 22 forks source link

Something goes wrong with the dynamicRefineFvMesh #3

Closed EarlFan closed 3 years ago

EarlFan commented 3 years ago

Dear all,

Hi! I find this compressible AMR reacting solver this week but cannot successfully run the tutorial case with the AMR setup. To remove any possibly parallel cause, I try to run it in serial but still encounter the same problem. The outputs are

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 7
Exec   : rhoReactingCentralFoam
Date   : Jul 03 2021
Time   : 14:36:28
Host   : "FanE-Dell-P7910"
PID    : 32468
I/O    : uncollated
Case   : /mnt/d/01Code/OpenFOAM-v7-deton/rhoReactingCentralFoam-master/examples/simple_detonation
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicRefineFvMesh
Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleInternalEnergy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader foamChemistryReader
    elements not defined in "/mnt/d/01Code/OpenFOAM-v7-deton/rhoReactingCentralFoam-master/examples/simple_detonation/constant/reactions"
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating reaction model

Selecting combustion model laminar
Selecting chemistry solver 
{
    solver          ode;
    method          standard;
}

StandardChemistryModel: Number of species = 4 and reactions = 1
Selecting ODE solver RKF45
    using integrated reaction rate
No finite volume options present
fluxScheme: Kurganov

Starting time loop

Central Courant Number mean: 0.00030816535 max: 0.00048901458 acoustic max: 0.00082955963
Time = 1.19904e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for Ux, Initial residual = 6.298955e-09, Final residual = 2.0589634e-24, No Iterations 1
DICPBiCGStab:  Solving for Uy, Initial residual = 3.7303124e-08, Final residual = 1.1957133e-23, No Iterations 1
DICPBiCGStab:  Solving for Uz, Initial residual = 1.8931358e-08, Final residual = 9.2659172e-24, No Iterations 1

min / max mag(U) ; 5.852886e-20 / 0.077089524

DILUPBiCGStab:  Solving for O2, Initial residual = 1, Final residual = 2.8512151e-21, No Iterations 1
DILUPBiCGStab:  Solving for H2, Initial residual = 0.99999998, Final residual = 2.8512082e-21, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 1, Final residual = 2.8536087e-21, No Iterations 1
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for e, Initial residual = 3.6335901e-08, Final residual = 3.2187083e-23, No Iterations 1
min/max(T) = 300.31668, 1340.4837
ExecutionTime = 0.3 s  ClockTime = 1 s

Selected 120 cells for refinement out of 2500.
Refined from 2500 to 3340 cells.
Selected 0 split points out of a possible 120.
Central Courant Number mean: 0.00039263623 max: 0.001235783 acoustic max: 0.0021152427
Time = 2.63616e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for Ux, Initial residual = 2.9268134e-09, Final residual = 2.110105e-23, No Iterations 1
DICPBiCGStab:  Solving for Uy, Initial residual = 5.0411427e-08, Final residual = 5.9207122e-23, No Iterations 1
DICPBiCGStab:  Solving for Uz, Initial residual = 0.99992818, Final residual = 2.8884616e-28, No Iterations 1

min / max mag(U) ; 1.6635645e-08 / 0.13952852

DILUPBiCGStab:  Solving for O2, Initial residual = 0.093971111, Final residual = 6.4021242e-23, No Iterations 1
DILUPBiCGStab:  Solving for H2, Initial residual = 0.093971683, Final residual = 6.4015162e-23, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0.093488949, Final residual = 6.3969465e-23, No Iterations 1
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for e, Initial residual = 3.0438376e-08, Final residual = 2.6195759e-22, No Iterations 1
min/max(T) = 300.3853, 1340.7741
ExecutionTime = 0.39 s  ClockTime = 1 s

Selected 980 cells for refinement out of 3340.
Refined from 3340 to 10200 cells.
Selected 0 split points out of a possible 980.
Central Courant Number mean: 0.00052878705 max: 0.0029699838 acoustic max: 0.0051565472
Time = 4.35873e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for Ux, Initial residual = 1.0976751e-08, Final residual = 1.937347e-21, No Iterations 1
DICPBiCGStab:  Solving for Uy, Initial residual = 6.0840305e-08, Final residual = 2.0745518e-21, No Iterations 1
DICPBiCGStab:  Solving for Uz, Initial residual = 0.56786968, Final residual = 6.2046291e-26, No Iterations 1

min / max mag(U) ; 6.032677e-08 / 0.33428277

DILUPBiCGStab:  Solving for O2, Initial residual = 0.023432807, Final residual = 1.0843779e-22, No Iterations 1
DILUPBiCGStab:  Solving for H2, Initial residual = 0.02343298, Final residual = 1.0844077e-22, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0.02328678, Final residual = 1.0844476e-22, No Iterations 1
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for e, Initial residual = 2.0772664e-08, Final residual = 7.8286085e-21, No Iterations 1
min/max(T) = 300.4675, 1341.1923
ExecutionTime = 0.66 s  ClockTime = 1 s

Selected 6400 cells for refinement out of 10200.
Refined from 10200 to 55000 cells.
Selected 0 split points out of a possible 6580.
Central Courant Number mean: 0.00074354821 max: 0.0071236897 acoustic max: 0.012368057
Time = 6.42437e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for Ux, Initial residual = 2.9213247e-08, Final residual = 1.6807225e-19, No Iterations 1
DICPBiCGStab:  Solving for Uy, Initial residual = 7.6714466e-08, Final residual = 7.2321903e-20, No Iterations 1
DICPBiCGStab:  Solving for Uz, Initial residual = 5.912485e-05, Final residual = 2.9147644e-16, No Iterations 1

min / max mag(U) ; 1.4660595e-07 / 0.8049728

DILUPBiCGStab:  Solving for O2, Initial residual = 0.017900724, Final residual = 1.323338e-21, No Iterations 1
DILUPBiCGStab:  Solving for H2, Initial residual = 0.017900876, Final residual = 1.3233439e-21, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0.017772355, Final residual = 1.3238858e-21, No Iterations 1
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
DICPBiCGStab:  Solving for e, Initial residual = 2.6775532e-08, Final residual = 2.9563792e-19, No Iterations 1
min/max(T) = 300, 1341.6952
ExecutionTime = 1.47 s  ClockTime = 2 s

Selected 0 cells for refinement out of 55000.
Selected 3380 split points out of a possible 6580.
Unrefined from 55000 to 31340 cells.
#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::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/fane/OpenFOAM/fane-7/platforms/linux64GccDPInt32Opt/bin/rhoReactingCentralFoam"
#5  ? in "/home/fane/OpenFOAM/fane-7/platforms/linux64GccDPInt32Opt/bin/rhoReactingCentralFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? in "/home/fane/OpenFOAM/fane-7/platforms/linux64GccDPInt32Opt/bin/rhoReactingCentralFoam"
Floating point exception

The problem arises from the first unrefined process.

By the way, when using the static grid by switch the dynamicFvMesh in dynamicMeshDict to staticFvMesh, the code works OK. In addition, I use the OpenFOAM 7 with these information in the README.org

#+TITLE:     README for [[http://openfoam.org/download/source][OpenFOAM-7]]
#+AUTHOR:               The OpenFOAM Foundation
#+DATE:                     8th July 2019
#+LINK:                  http://openfoam.org
#+OPTIONS: author:nil ^:{}
# Copyright (c) 2015-2019 OpenFOAM Foundation.

Any advice will be appreciate!

Regards, Fan E

EarlFan commented 3 years ago

Dear all,

I solved this problem by installing a latest OpenFOAM-v7, so I will close this issue.

Have a good day.

Regards, Fan E