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

coding a steady multiphase solver twophaseeulerfoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 13, 2011, 13:16
Default coding a steady multiphase solver twophaseeulerfoam
  #1
Member
 
Miles
Join Date: Sep 2011
Posts: 48
Rep Power: 5
miles_davis is on a distinguished road
Hi,

I am using the twophaseeeulerfoam solver to simulate air bubbles in water.
I want to have a steady state version of the solver.

Solution1
According to what I have read, I should first try to under-relax the resolution and raise the time step.
I do not understand why? In fact I do not understand the role of under-relaxtion factor in an unsteady solver.
Does anyone have a point on the subject or some documentation of the relaxation process in Openfoam (and in general, reading is never bad!)? This could help me to understand.

Solution2
On the other hand if I want to change the solver to be steady witch library I was thinking to do so
- Remove time dependence in the solver
- Change the vel-press scheme (I was thinking using simple scheme)
- Modify the fvsol and fvscheme dict.
Does this seem correct?
Does this seem an affordable task for an OF newbee?

Thanks for your opinions and help.

Regards,
miles_davis is offline   Reply With Quote

Old   October 13, 2011, 13:16
Default
  #2
Member
 
Miles
Join Date: Sep 2011
Posts: 48
Rep Power: 5
miles_davis is on a distinguished road
Hi

I have tried the first solution but it seems that putting a subdict for relaxation at the end of fvSolution has no effect.
I am not sure if it is taken into account.

So I switched to the coding of a steady twophaseeluerfoam solver.
I am tryning to understand how the solver is coded.


I have a few questions and my sorry in advance if it seems to easy but my C++ is very bad.
1.fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime)
Is it possible to modify the definition of fvVectorMatrix like that:
fvVectorMatrix UbEqn(Ub, Ub.dimensions())


  1. volScalarField rUaA(1.0/UaEqn.A())
I do not understand what is done here.
What is the definition and the purpose of rUaA and UaEqn.A() ?
Is this related to the first guess on momentum equation in the PIMPLE algorithm?


  1. Ua = rUaA*UaEqn.H()
Same question for UaEqn.H()


4. DDTU.h
What is the purpose of having a special file for the equation on UA and UB?
Can it be included in my SIMPLE loop -since I will move the scheme from pimple to simple.


Thanks for any help.


Regards,
miles_davis is offline   Reply With Quote

Old   October 14, 2011, 20:52
Default
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Just a suggestion: do not post the same message multiple times, or you won't receive any answer.

Best,

Alberto
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   October 15, 2011, 19:45
Default
  #4
Member
 
Miles
Join Date: Sep 2011
Posts: 48
Rep Power: 5
miles_davis is on a distinguished road
Sorry for that

That was not my intention.
I wanted to move the thread to the proper section and delete the older. But I did'nt manage to find how to delete a thread...

orry again.

Miles
miles_davis is offline   Reply With Quote

Old   October 15, 2011, 19:49
Default
  #5
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by miles_davis View Post
Sorry for that

That was not my intention.
I wanted to move the thread to the proper section and delete the older. But I did'nt manage to find how to delete a thread...

orry again.

Miles
Ooops, I did not pay attention to the sections Sorry. I answered in the other thread.

I don't know if it's possible to delete a thread without administrative rights. :-?
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
convergence issues on steady state solver icemaniac178 CFX 1 March 30, 2011 19:11
time step and iterations in steady state problem using transient solver jing113cn FLUENT 2 January 15, 2010 04:18
Steady state solver for chemical reaction hamburgFoam OpenFOAM Running, Solving & CFD 2 January 8, 2010 09:34
Looking for a steady state supersonic solver achilleslaststand OpenFOAM Running, Solving & CFD 0 December 20, 2009 10:23
steady state and multiphase... ben Phoenics 5 April 25, 2003 05:48


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