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

Instabilities related to the geometric source term

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 14, 2008, 22:24
Default Instabilities related to the geometric source term
  #1
jinwon park
Guest
 
Posts: n/a
For symmetric flows, I formulated symmetric 1-D Euler equations and axisymmetric 2-D Euler equations. To numerically solve the universal Euler equations with sources, I applied operator splitting method to the set of equations as dU/dt+dF/dr=S(r) in 1-D and dU/dt+dF/dr+dG/dz=S(r,z). For example, I first solve the standard Euler equation then the solution becomes the initial condition for source terms.

1. dUbar/dt+dF/dr=0 giving rise to Ubar

2. dU/dt=S

When I include the source terms for symmetric solutions, I suffered from stability troubles. However, when I do not include these, the code performed great. Could anyone advice me about the method I used?
  Reply With Quote

Old   February 17, 2008, 22:17
Default Re: Instabilities related to the geometric source
  #2
Patrick
Guest
 
Posts: n/a
Why dont you directly solve including the source term ? Is the source term an implicit function of the primitive variables? or just some terms of the form rho*v_r etc?

What time scheme do you use? A Runge-Kutta or a simple leap-frog? or an implicit method?

I am not sure I understand the reason behind performind such a splitting ?
  Reply With Quote

Old   February 17, 2008, 22:21
Default Re: Instabilities related to the geometric source
  #3
jinwon park
Guest
 
Posts: n/a
The source terms are jsut some terms such as rho*v as you told. The reason I solved it by splitting methods was based on the statement in the literature. In that literature(E.F. Toro, Riemann problem...), the operator splitting method is only one suggested. Can we solve it by direct coupling with advection terms? If so, it will be better to the simple implementation. Have you seen any literature recommending for coupled solution approach?
  Reply With Quote

Old   February 20, 2008, 10:44
Default Re: Instabilities related to the geometric source
  #4
Patrick
Guest
 
Posts: n/a
Hi, sorry for being late to respond...

Are solving the equations with an implicit time scheme or an explicit ?

If it is explicit then this is straightforward, if implicit then it involves something like inversing a matrix (or some iteration of some kind).

I have use spliting only to split space dimensions and to split "implicit" terms from "explicit" terms (again as reltated to time scheme).

I must confess that it is kind of difficult to follow the problem without seeing the original equations and the numerical equations.
  Reply With Quote

Old   February 20, 2008, 10:50
Default Re: Instabilities related to the geometric source
  #5
jinwon park
Guest
 
Posts: n/a
I am using an explicit scheme. Based on the previous post, I substitute the source term into the right hand side. I am testing many cases. In 1-D, it worked well. Thanks for your reply.

Jinwon
  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
momentum source term zwdi FLUENT 14 June 27, 2017 16:40
How to define a transient species source term? aleisia FLUENT 3 March 7, 2011 23:58
[Gmsh] Compiling gmshFoam with OpenFOAM-1.5 BlGene OpenFOAM Meshing & Mesh Conversion 10 August 6, 2009 05:26
UDF Source Term Units? Brian FLUENT 1 October 24, 2005 10:15
UDF's - linearizing a negative source term Steve Howell FLUENT 1 March 18, 2001 17:21


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