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

Static equation solving

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By HamzauTT

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 28, 2023, 05:58
Default Static equation solving
  #1
New Member
 
Hamza El Fathi
Join Date: Aug 2023
Posts: 6
Rep Power: 3
HamzauTT is on a distinguished road
Hello,
I'm trying to solve a static equation in OpenFOAM (without time derivatives). Initially, I removed the time derivative from the scalarTransportFoam solver and the simple time loop with no results. I also tried retaining the time loop and running the simulation for one step, yet didnt work either. I'm seeking assistance here to know whether anyone is aware of a static solver (no time derivative term in the equation) within OpenFOAM that I could potentially adapt for my specific case.
HamzauTT is offline   Reply With Quote

Old   August 28, 2023, 11:42
Default
  #2
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: https://olesenm.github.io/
Posts: 1,714
Rep Power: 40
olesen has a spectacular aura aboutolesen has a spectacular aura about
Quote:
Originally Posted by HamzauTT View Post
Hello,
I'm trying to solve a static equation in OpenFOAM (without time derivatives). Initially, I removed the time derivative from the scalarTransportFoam solver and the simple time loop with no results. I also tried retaining the time loop and running the simulation for one step, yet didnt work either. I'm seeking assistance here to know whether anyone is aware of a static solver (no time derivative term in the equation) within OpenFOAM that I could potentially adapt for my specific case.

laplacianFoam? It doesn't get much simpler than that.
olesen is offline   Reply With Quote

Old   August 29, 2023, 10:45
Default
  #3
New Member
 
Hamza El Fathi
Join Date: Aug 2023
Posts: 6
Rep Power: 3
HamzauTT is on a distinguished road
I checked laplacian foam and it seems to have a time loop.
HamzauTT is offline   Reply With Quote

Old   August 29, 2023, 12:38
Default
  #4
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 798
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
laplacianFoam does indeed have a while(simple.loop()).

See e.g. line 85 of https://develop.openfoam.com/Develop...aplacianFoam.C

Even if ddtSchemes in systems/fvSchemes is set steadySteady, the loop remains (most obviously).

Two possible scenarios exist.

1/ The source term is linear in the state variable solved for, and the residual of the SIMPLE iteration drops to a very small number in only one iteration. Only one iteration is required.

2/ The source term is non-linear in the state variable solved for, and various iterations are required to bring the residual of the SIMPLE iteration to a small number.
dlahaye is offline   Reply With Quote

Old   August 30, 2023, 06:14
Default
  #5
New Member
 
Hamza El Fathi
Join Date: Aug 2023
Posts: 6
Rep Power: 3
HamzauTT is on a distinguished road
The simulation does indeed stop after the first iteration.But it does not provide any notification of the solving process for that iteration, and I do not receive any output.
HamzauTT is offline   Reply With Quote

Old   August 30, 2023, 08:22
Default
  #6
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 798
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
Copy and paste here output from terminal that you receive.
dlahaye is offline   Reply With Quote

Old   August 30, 2023, 09:41
Default
  #7
New Member
 
Hamza El Fathi
Join Date: Aug 2023
Posts: 6
Rep Power: 3
HamzauTT is on a distinguished road
Here is it and thanks for helping:


Create time

Create mesh for time = 0


SIMPLE: no convergence criteria found. Calculations will run for 1 steps.

Reading field T

Reading field U

Reading transportProperties

Reading diffusivity DT

Reading/calculating face flux field phi

No finite volume options present

Calculating scalar transport

Courant Number mean: 0.0495 max: 0.0495
Time = 0.001

hamza469@DESKTOP:~/1D$
HamzauTT is offline   Reply With Quote

Old   August 30, 2023, 10:41
Default
  #8
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 798
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
Thanks.

Pls share controlDict, fvSolution and fvSchemes.
dlahaye is offline   Reply With Quote

Old   September 1, 2023, 11:33
Default
  #9
New Member
 
Hamza El Fathi
Join Date: Aug 2023
Posts: 6
Rep Power: 3
HamzauTT is on a distinguished road
Thank you for your help. The steadyState option is indeed more practical than modifying the solver. I figured out that my mistake was that the boundary conditions were not defined correctly.
dlahaye likes this.
HamzauTT is offline   Reply With Quote

Reply


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
chtMultiRegionSimpleFoam turbulent case Aditya Patil OpenFOAM Running, Solving & CFD 6 April 24, 2017 23:13
HeatSource BC to the whole region in chtMultiRegionHeater xsa OpenFOAM Running, Solving & CFD 3 November 7, 2016 06:07
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 05:03
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 04:49.