CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

About theory:how to get the pressure equation when solving u-p simultaneously?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 28, 2012, 22:33
Default About theory:how to get the pressure equation when solving u-p simultaneously?
  #1
New Member
 
zhiwei liu
Join Date: Mar 2009
Posts: 22
Rep Power: 8
lzw2003 is on a distinguished road
Hi,

Now i want to solve u,v,w,p simultaneously using coupled solver. But there is no equation for pressure. So i think the important is to get a pressure euqation.A method i know is to use mass conservation equation and equation of state to get a pressure equation.
I think this method is not suitable for my situation.

So,
Need the mass conservation equation special handle?Or
is there any other way to get the pressure equation ?
lzw2003 is offline   Reply With Quote

Old   March 28, 2012, 22:38
Default
  #2
New Member
 
zhiwei liu
Join Date: Mar 2009
Posts: 22
Rep Power: 8
lzw2003 is on a distinguished road
Or how to solve u,v,w,p simultaneously compared to segregated solver?
lzw2003 is offline   Reply With Quote

Old   March 29, 2012, 06:14
Default
  #3
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 372
Rep Power: 10
arjun is on a distinguished road
I am afraid not many people can help you on it here.

Anyway to answer your question, there are two main class of treating this whole thing, but both of them will converge to similar approaches in the end.

First approach is not to construct pressure equation but rather create a saddle system.

| Au G | | u | = src_u
| D 0 | | p | = src_p

And solve this. Solving this could be done with the help of GMRES or BiCGStab. To run Kyrlov method you will need preconditioner and this preconditioner could be constructed in many ways.

Search for Preconditioner for saddle system or for Stokes problem.

Have a look at Ales Janka's work too. He has some ppts that Summarizes it very well.


Second approach is to create a velocity and pressure correction matrix based on finite volume approach. Here velocity correction could be written in terms of pressure correction gradients (SIMPLE Method). And this system will come out to be of form:


| Au G | | u' | = src_u
| D Ap | | p' | = src_p

Again you can apply same methods as above to solve this system or you could construct coupled AMG on it. Look at Ales Janka's work.


If you read this whole business in deep you will find that ultimately after some re-arragements you are solving the same thing in both approaches.

PS: If you are using AMG, traditional Gauss-Seidel does not work well as smoother. Braess Sarazin is good smoother but is costly. (Try it first because it is surefire method, if you get it working experiment with other things).


I hope this helps.
arjun is offline   Reply With Quote

Old   March 29, 2012, 21:36
Default
  #4
New Member
 
zhiwei liu
Join Date: Mar 2009
Posts: 22
Rep Power: 8
lzw2003 is on a distinguished road
Hi, arjun:
thank you for this. I think i should study the first approach,maybe it will be helpful.

And i also want to know that is it identity when i transfer the mass conservation equation base on density (Eq.9.1)to that base on pressure(Eq. 2.13). Will the mass conserve if i used FVM to handle (2.13)?




Attached Images
File Type: jpg massmomentum.jpg (74.9 KB, 60 views)
File Type: jpg mass.jpg (4.6 KB, 59 views)
lzw2003 is offline   Reply With Quote

Old   March 29, 2012, 22:57
Default
  #5
New Member
 
zhiwei liu
Join Date: Mar 2009
Posts: 22
Rep Power: 8
lzw2003 is on a distinguished road
Hi,arjun
Now I reading some papers, there is a method for coupled solver, using the Rhie-Chow interpolation
(momentum interpolation method, MIM),
proposed by Rhie-Chow
(Paper name: Numerical study of the turbulent flow past an airfloil with trailing edge separation) .
it is a method for coupled solver in co_located grid. I think this is suitable for me.
I dont know if you are familar in this field. But I want to know is there any other method for the coupled solver in co-located grid?

Thanks!
lzw2003 is offline   Reply With Quote

Old   March 31, 2012, 15:29
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,651
Rep Power: 23
FMDenaro will become famous soon enough
but is your goal to fulfill div V=0 ? Uzawa method was an old unsplit procedure...
FMDenaro is offline   Reply With Quote

Old   April 2, 2012, 00:17
Default
  #7
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 372
Rep Power: 10
arjun is on a distinguished road
Quote:
Originally Posted by lzw2003 View Post
Hi, arjun:
thank you for this. I think i should study the first approach,maybe it will be helpful.

And i also want to know that is it identity when i transfer the mass conservation equation base on density (Eq.9.1)to that base on pressure(Eq. 2.13). Will the mass conserve if i used FVM to handle (2.13)?





I have not implemented coupled solver for compressible case so I can not say much. But my understanding is that density based solvers are already coupled in nature. (the way they are implemented).
And for pressure based coupled I am sure with using Rhie Chow momentum interpolation and other terms coupled solver could be created. (Fluent it seems has one). I have no time to dwell onto it so did not do it.


Quote:
Originally Posted by lzw2003 View Post
Hi,arjun
Now I reading some papers, there is a method for coupled solver, using the Rhie-Chow interpolation
(momentum interpolation method, MIM),
proposed by Rhie-Chow
(Paper name: Numerical study of the turbulent flow past an airfloil with trailing edge separation) .
it is a method for coupled solver in co_located grid. I think this is suitable for me.
I dont know if you are familar in this field. But I want to know is there any other method for the coupled solver in co-located grid?

Thanks!
I implemented coupled solver based on this approach and it is not difficult to do.

This is my break of time for implementing coupled version in iNavier.
3 days to create coupled multigrid for coupled system of equation.
1 day to convert segregated algorithm into coupled one.
2 hours to create 2 D version of coupled solver.

So my advise first concentrate of matrix solver for this system. Once you have it working things are very easy after that.


Edited to add: Look for work done by M. Darwish also on coupled solver. His papers give very good explanation on implementing it.
arjun is offline   Reply With Quote

Old   April 8, 2012, 21:14
Default
  #8
New Member
 
zhiwei liu
Join Date: Mar 2009
Posts: 22
Rep Power: 8
lzw2003 is on a distinguished road
Yes, my goal is to fulfill the div V=0.
lzw2003 is offline   Reply With Quote

Old   April 8, 2012, 22:13
Default
  #9
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 372
Rep Power: 10
arjun is on a distinguished road
Quote:
Originally Posted by lzw2003 View Post
Yes, my goal is to fulfill the div V=0.

All the methods are trying to do exactly this. Uzawa's method is slow and can not compete with AMG based methods.
arjun 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
lift and drag on ship superstructures vaina74 OpenFOAM Running, Solving & CFD 3 June 8, 2010 12:30
MRFSimpleFOAM goes divergenced! renyun0511 OpenFOAM Running, Solving & CFD 0 November 19, 2009 03:11
Parallel rasInterFoam openfoam_user OpenFOAM Running, Solving & CFD 4 November 1, 2008 05:14
MRFSimpleFoam amp cyclic patches david OpenFOAM Running, Solving & CFD 36 October 21, 2008 21:55
Negative value of k causing simulation to stop velan OpenFOAM Running, Solving & CFD 1 October 17, 2008 05:36


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