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

Solving the Psystem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 2, 2005, 09:15
Default I wish to test the numerical m
  #1
Member
 
Marco Kupiainen
Join Date: Mar 2009
Posts: 31
Rep Power: 7
kupiainen is on a distinguished road
I wish to test the numerical methods in OpenFOAM on a well known hyperbolic test-problem, namely the p-system:
[v u]_t + [-u A*v^-gamma]_x = 0
I have written a piece of code that solves this in 1D, using the OpenFOAM space-discretization-operator div. The time-stepping is made with a 2:nd order TVD Runge-Kutta. My question is the following:
The convection-specific schemes calculate the interpolation based on the flux of the flow velocity. What should I consider as the flow velocity here?
kupiainen is offline   Reply With Quote

Old   May 2, 2005, 11:27
Default I have solved the "shallow-wat
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 12
henry is on a distinguished road
I have solved the "shallow-water" equations in various ways but never the "p-system". In the former the choice of flux variable is obvious but in the later it is definitely not obvious. What physical problem can the "p-system" be applied to?
henry is offline   Reply With Quote

Old   May 2, 2005, 12:31
Default The p-system is a model for is
  #3
Member
 
Marco Kupiainen
Join Date: Mar 2009
Posts: 31
Rep Power: 7
kupiainen is on a distinguished road
The p-system is a model for isentropic (=constant entropy) or polytropic gas given in a Lagrangian coordinate system.
v is the specific volume i.e. 1/rho and u is the velocity, and p=A*v^-gamma is the pressure (A>0 and gamma=1.4 typically for air).
I have defined both v and u as volScalarFields, and by this choice I can use the linear interpolation scheme. However since I want to solve a Riemann problem I want to test all the different TVD methods already implemented in OpenFOAM. What is the easiest solution to this?
kupiainen is offline   Reply With Quote

Old   May 2, 2005, 12:34
Default If u is the velocity then you
  #4
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 12
henry is on a distinguished road
If u is the velocity then you should construct a flux from it somehow, ideally by constructing a pressure equation and getting the flux from that but if that is not appropriate then by linear interpolation.
henry 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
Cannot see residuals while solving Dre CFX 3 January 29, 2009 11:00
Why solving kar OpenFOAM Running, Solving & CFD 2 April 30, 2008 15:15
FvVectorMatrixresidual before solving cosimobianchini OpenFOAM Running, Solving & CFD 0 January 28, 2007 11:02
solving a Matrix Behafarid Main CFD Forum 3 September 9, 2006 23:14
MHD solving? Rajil Saraswat Main CFD Forum 1 March 21, 2003 11:26


All times are GMT -4. The time now is 22:27.