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

Forcing update on boundary fields after each iteration

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 10, 2018, 20:43
Default Forcing update on boundary fields after each iteration
  #1
Member
 
Pavan
Join Date: Jan 2016
Posts: 53
Rep Power: 9
pvpnrao is on a distinguished road
I am using pimplefoam in Openfoam 5.x.
The boundary condition seems to update only once per timestep.



I need this updating to take place after each pimple loop as well as at the end of the time-step.



Can some one suggest how to edit pimpleFoam to achieve this?
pvpnrao is offline   Reply With Quote

Old   September 11, 2018, 11:18
Default
  #2
Member
 
Pavan
Join Date: Jan 2016
Posts: 53
Rep Power: 9
pvpnrao is on a distinguished road
Any clues?
pvpnrao is offline   Reply With Quote

Old   September 11, 2018, 11:42
Default
  #3
Senior Member
 
anonymous
Join Date: Jan 2016
Posts: 416
Rep Power: 13
simrego is on a distinguished road
Hi!


As i know it is updated in every equation solving. Check the code or just simply write a pretty dummy coded BC, for example:

https://cpp.openfoam.org/v6/classFoa...d.html#details
just set operator==(0);
and use as a pressure BC at outlet.

and also write out that the BC is evaluated before operator==;

ie:
Info<< "Here is my BC again..." << endl;


If you see this message, your BC is evaluated. You will see this message a lot in the logs!
simrego is offline   Reply With Quote

Old   September 11, 2018, 12:27
Default
  #4
Member
 
Pavan
Join Date: Jan 2016
Posts: 53
Rep Power: 9
pvpnrao is on a distinguished road
simrego


Thank you very much for the quick response.



I am printing the area averaged pressure at inlet boundary after each pimple loop. The inlet pressure fluctualtes in the first few pimple loops then settles for a fixed negative value for rest of the PIMPLE iterations.



I am intializing the case with converveged solution from a steady state run . The areaAverageInletPressure at the begining of the first time step shows the initial pressure. But the solution does not reach the same pressure at the end of the time step. At least in the first timestep I expect the inlet pressure to be equal to the initialized value.


My cluster system is under maintenence today, i will post the iteration log ASAP.
Meanwhile, any comments or suggestions will be helpful.
pvpnrao is offline   Reply With Quote

Old   September 11, 2018, 12:53
Default
  #5
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 351
Rep Power: 14
Santiago is on a distinguished road
Quote:
Originally Posted by pvpnrao View Post
simrego


Thank you very much for the quick response.



I am printing the area averaged pressure at inlet boundary after each pimple loop. The inlet pressure fluctualtes in the first few pimple loops then settles for a fixed negative value for rest of the PIMPLE iterations.



I am intializing the case with converveged solution from a steady state run . The areaAverageInletPressure at the begining of the first time step shows the initial pressure. But the solution does not reach the same pressure at the end of the time step. At least in the first timestep I expect the inlet pressure to be equal to the initialized value.


My cluster system is under maintenence today, i will post the iteration log ASAP.
Meanwhile, any comments or suggestions will be helpful.
There is a class member in geometrciFields called correctBoundaryConditions()

So, if you write

U.correctBoundaryConditions();

Youll get what you intend.
Santiago is offline   Reply With Quote

Old   September 11, 2018, 13:10
Default
  #6
Member
 
Pavan
Join Date: Jan 2016
Posts: 53
Rep Power: 9
pvpnrao is on a distinguished road
Santiago


I tried this option, please see the code below. I dont know whether I placed the code within the proper loop or not. Nonetheless, it did not work. Please suggest if needs to be placed somewhere else:


// --- Pressure-velocity PIMPLE corrector loop
while (pimple.loop())
{
#include "UEqn.H"

// --- Pressure corrector loop
while (pimple.correct())
{
#include "pEqn.H"
}

if (pimple.turbCorr())
{
laminarTransport.correct();
turbulence->correct();
}
U.correctBoundaryConditions();
}


Furthermore this systax is present in the list one of pEqn.H. So, I guess it is already being executed, but I could not see any difference.
pvpnrao 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
3D Windturbine simulation in SU2 k.vimalakanthan SU2 14 February 8, 2019 15:43
[snappyHexMesh] sHM layer process keeps getting killed MBttR OpenFOAM Meshing & Mesh Conversion 4 August 15, 2016 04:21
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 05:39
Waterwheel shaped turbine inside a pipe simulation problem mshahed91 CFX 3 January 10, 2015 12:19
an odd(at least for me!) reconstructPar error on a field immortality OpenFOAM Running, Solving & CFD 3 June 3, 2013 23:36


All times are GMT -4. The time now is 18:25.