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

FOAM warning: Step size underflow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 3, 2020, 11:11
Default FOAM warning: Step size underflow
  #1
New Member
 
Jiaui Zhang
Join Date: Aug 2018
Location: Changsha, China
Posts: 17
Rep Power: 8
charryzzz is on a distinguished road
hi Foamers:


I am running a self-developed solver with combustion using seulex as OED solver.


The case is running with a time step of [deltaT 5e-07];


Then, I will get such warning:


FOAM Warning :
From function seulex::solve(scalar& x, scalarField& y, stepState&
in file ODESolvers/seulex/seulex.C at line 255
step size underflow :8.470329473e-22



So does anyone know why this warning appears? Will this warning affect the calculation result?


Here is the link for seulex.c

https://github.com/OpenFOAM/OpenFOAM...eulex/seulex.C
charryzzz is offline   Reply With Quote

Old   May 3, 2020, 13:50
Default
  #2
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 13
HPE is on a distinguished road
It looks like the step size is very small to the extent that it cannot be represented in floating-point computation. Is varying time-step size setting on, so that your input for the time-step (i.e. deltaT 5e-07) has been changed during the simulation?
HPE is offline   Reply With Quote

Old   May 3, 2020, 14:31
Default
  #3
New Member
 
Jiaui Zhang
Join Date: Aug 2018
Location: Changsha, China
Posts: 17
Rep Power: 8
charryzzz is on a distinguished road
Quote:
Originally Posted by HPE View Post
It looks like the step size is very small to the extent that it cannot be represented in floating-point computation. Is varying time-step size setting on, so that your input for the time-step (i.e. deltaT 5e-07) has been changed during the simulation?
Yes, it seems to be the case, and there is no such warning when I reduce the time step from 5E-6 to 5E-7, but the case will cost a much longer time.
charryzzz is offline   Reply With Quote

Old   May 3, 2020, 14:49
Default
  #4
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 13
HPE is on a distinguished road
Use fixed time step, and fix the time-step 1e-7 or lower instead of using variable time step.
HPE is offline   Reply With Quote

Old   June 5, 2023, 10:35
Default step size underflow :0
  #5
Member
 
bany
Join Date: Nov 2019
Posts: 50
Rep Power: 7
bany is on a distinguished road
Quote:
Originally Posted by HPE View Post
Use fixed time step, and fix the time-step 1e-7 or lower instead of using variable time step.

Sorry for recalling this Thread.

When i use chemFoam to calculate the ignition delay time for PRF80 with seulex_KLU, the infinite warnings "step size underflow :0" occur, as shown in the attachment. The mechanism creates a Jacobian with a sparsity of 98%. The calculation is performed by chemFoam with fixed time step 1E-7. In fact, small time steps like 1E-8, 1E-9 are also attempted but with the same warnings.

seulex_KLU is the solver seulex with KLU sparse matrix library, provided in the article https://doi.org/10.1016/j.combustflame.2016.09.018.
Attached Images
File Type: png seulex_KLU-warnings.png (46.0 KB, 11 views)
bany is offline   Reply With Quote

Old   May 9, 2024, 02:23
Default
  #6
Member
 
bany
Join Date: Nov 2019
Posts: 50
Rep Power: 7
bany is on a distinguished road
Quote:
Originally Posted by bany View Post
Sorry for recalling this Thread.

When i use chemFoam to calculate the ignition delay time for PRF80 with seulex_KLU, the infinite warnings "step size underflow :0" occur, as shown in the attachment. The mechanism creates a Jacobian with a sparsity of 98%. The calculation is performed by chemFoam with fixed time step 1E-7. In fact, small time steps like 1E-8, 1E-9 are also attempted but with the same warnings.

seulex_KLU is the solver seulex with KLU sparse matrix library, provided in the article https://doi.org/10.1016/j.combustflame.2016.09.018.

The problem has been solved. We changed another computer and cluster and the warning disappeared.
bany is offline   Reply With Quote

Reply

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
[Other] refineWallLayer Error Yuby OpenFOAM Meshing & Mesh Conversion 2 November 11, 2021 11:04
CFD by anderson, chp 10.... supersonic flow over flat plate varunjain89 Main CFD Forum 18 May 11, 2018 07:31
is internalField(U) equivalent to zeroGradient? immortality OpenFOAM Running, Solving & CFD 7 March 29, 2013 01:27
[Commercial meshers] ST_Malloc: out of memory.malloc_storage: unable to malloc Velocity SA, cfdproject OpenFOAM Meshing & Mesh Conversion 0 April 14, 2009 15:45
[Gmsh] Import gmsh msh to Foam adorean OpenFOAM Meshing & Mesh Conversion 24 April 27, 2005 08:19


All times are GMT -4. The time now is 22:23.