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

Custom BC (Temperature, cyclic) in compressible solver (rhoSimpleFOAM) crashing

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 14, 2022, 06:10
Post Custom BC (Temperature, cyclic) in compressible solver (rhoSimpleFOAM) crashing
  #1
New Member
 
D. Isik
Join Date: Jan 2019
Location: Karlsruhe, Germany
Posts: 1
Rep Power: 0
memento is on a distinguished road
Hi Foamers,

I implemented my own custom Temperature boundary condition in OpenFOAM. While the incompressible case is running fine and the solution is indeed plausible, my problem emerges when I port the bc to a compressible solver.


My BC is for cyclic cases, I'm using AMI.
I use the fixedJumpAMI condition and developed my BC based on that.

As I said above when runnning in simpleFOAM everything works fine and as intended.
But during the thermo.correct() in EEqn.H the solver crashes with SigFpe.


So as I debugged, multiple things seemed odd:
  • In one case I decomposed the domain and usually the solver crashes as mentioned above during the thermo.correct().
    Here I started writing the gMin/gMax values of the used fields (T, p, psi, rho, mu, he and alpha) on all boundaries. Now the solvers crashes with SigFpe during evaluation of gMin/gMax on the processor patches.
  • I let the solver write those fields before the thermo.correct() and the values of he seemed odd (the he transport eq. is solved at this point). All values of e are in the range of -1e8 but on the inlet patch there are 2 patch face values with +4e6 and +2e6. The list reads as follows:
Code:
inlet
 {
  type         energyJumpAMI;
  patchType cyclicAMI;
  jump         uniform -17957379;
  value        nonuniform List<scalar>
80
(
4771773.9
2337384.7
-1849468.6
-4666132.8
...
)
}
Does anybody already have experience implementing a custom BC for Temperature in compressible cases?
Is there something I haven't considered?
I'm really stuck and moving in circles.
memento is offline   Reply With Quote

Reply

Tags
boundary condition, compressible

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
Divergence detected in AMG solver. VOF. Mr.Mister Fluent Multiphase 4 September 25, 2023 07:10
Fail to converge when solving with a fabricated solution zizhou FLUENT 0 March 22, 2021 06:33
[ANSYS Meshing] Help with element size sandri_92 ANSYS Meshing & Geometry 14 November 14, 2018 07:54
Ansys CFX problem: unexpected very high temperatures in premix laminar combustion faizan_habib7 CFX 4 February 1, 2016 17:00
compressible unsteady high Re and Temperature Flow: which solver to use? loreasr OpenFOAM Running, Solving & CFD 0 November 25, 2015 09:15


All times are GMT -4. The time now is 19:03.