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

New solver: ddtFoam

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

Like Tree20Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 21, 2023, 03:56
Default convergence pb with ddtFoam
  #61
New Member
 
HÚlo´se Magliano
Join Date: Feb 2023
Posts: 4
Rep Power: 3
Helo´seM is on a distinguished road
Hi everyone, I have a problem with the following case :

I try to simulate a rotational detonation combustor with hydrogen combustion. I'm using the ddtFoam solver from OpenFoam.
The computational domain is an annular cylinder surface which i have created with blockMesh and the checkMesh indicates that everything is okay.
I want to inject uniformly stoichiometric H2ľair at 293 K at the inlet. I have part of my domains also that have specific settings of chemicals, pressure and temperture. The result of my setFieldsDict is the image below:

setfields.png

The ignition zone is set with the mixture at high temperature and pressure in order to auto-ignate.
Deflagration is modeled through a reaction progress variable c for each cell (with c 0 representing an unburnt particle and c 1 representing a burnt particle).

At the inlet the pressure should be 1bar and temperature 293K.
My courant number is set to be 0.3 max.

My simulations are crashing but i don't know why: my courant number is increasing and so DeltaT is decreasing to compensate.
Here is the error message i get :

Acoustic Courant Number mean: 9.22991e-05 max: 0.300978
deltaT = 3.65653e-10
Time = 0.00000962374, time step = 275
Quenching 9.57776e-07 % of the turbulent reaction
Calculating Riemann fluxes, second order = true
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for Ux, Initial residual = 0.00726394, Final residual = 3.66489e-12, No Iterations 1
DICPCG: Solving for Uy, Initial residual = 0.00455063, Final residual = 2.40089e-12, No Iterations 1
DICPCG: Solving for Uz, Initial residual = 0.00314578, Final residual = 1.56327e-12, No Iterations 1
DICPCG: Solving for fH, Initial residual = 4.19317e-05, Final residual = 1.14321e-11, No Iterations 1
diagonal: Solving for tauHI, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for tauLO, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for c, Initial residual = 3.79914e-05, Final residual = 9.35827e-14, No Iterations 1
Combustion progress = 99.13%
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoEu, Initial residual = 0, Final residual = 0, No Iterations 0
[7]
[7]
[7] --> FOAM FATAL ERROR:
[7] Maximum number of iterations exceeded
Test = 11198.4 K,
eu_ = 1.20943e+07 J/kg, e(Test,fH) = 1.06602e+07 J/kg
[7]
[7] From function eUnburnedThermo::correct(const volScalarField& fH)
[7] in file eUnburnedThermo/eUnburnedThermo.C at line 413.
[7]
FOAM parallel run aborting
[7]
[7] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[7] #1 Foam::error::abort() at ??:?
[7] #2 Foam::eUnburnedThermo::correct() at ??:?
[7] #3
[7] at ??:?
[7] #4 __libc_start_main in "/lib64/libc.so.6"
[7] #5
[7] at ??:?
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 7 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.


My files are attached to this post.
If someone know what my problem could be please ? (maybe with my pressure BC ?)



system.zip

0.zip

constant.zip
Helo´seM is offline   Reply With Quote

Reply

Tags
combustion, compressible, density based, detonation, transonic

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Fluent Adjoint Solver? ex10148 FLUENT 16 September 28, 2018 08:11
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
Divergence problem Smaras FLUENT 13 February 21, 2013 05:03
3d vof Smaras FLUENT 2 February 19, 2013 06:58
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


All times are GMT -4. The time now is 17:59.