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

Fractional step in ALE form - possible?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 29, 2006, 23:07
Default Fractional step in ALE form - possible?
  #1
Wee
Guest
 
Posts: n/a
hello,

i'm using the fractional step on a structured grid to solve flow past a cylinder. discretization is based on FVM, CN2 for viscous, AB2 for convective term. the non-dimensionalized eqns in integral form are:

f[ (u(*)-u(n))/dt ]dV + f[ 1.5*u(n)*U(n)-0.5*u(n-1))*U(n-1)) ]dS = 0.5*f[ L(u(n)+u(*)).n ]dS /Re

f[ dp/dn ]dS = f[ U(*).n ]dS /dt

u(n+1)=u(*)- dt*(G(p)) (cc)

U(n+1)=U(*)- dt*(dp/dn) (fc)

f[ ] is the integration summation

U is the face normal velocity, u is the cell center

L is the discrete laplace operator, G is the discrete gradient.

cc mean G(p) is evaluated at the cell center

fc mean dp/dn is evaluated at the face center to prevent pressure oscillation.

i intent to modified it to the ALE (Arbitrary Lagrangian-Eulerian) form to simulate oscillating cylinder.

searching the net shows that a lot of ppl use FEM or FVM with artificial compressibility. there are some who use SIMPLE/FVM with ALE. however i can't seem to find ppl using ALE in fractional step. i wonder why? is it just preference or...?

from my understanding, the main modification in the ALE form is the convective velocity, which in this case is the changing of U (face normal velocity) to U - Ub, where Ub is the face normal mesh velocity.

hence, the overall change is only in the momentum and poisson eqn, where all the U will become U - Ub.

also, the change in BC will be wall velocity from 0 to cylinder velocity, that's all.

with these modification, i'll be able to simulate oscillating cylinder.

  Reply With Quote

Old   September 30, 2006, 18:47
Default Re: Fractional step in ALE form - possible?
  #2
Mani
Guest
 
Posts: n/a
I suppose you'll have to define some grid velocity associated with the intermediate position (*). Also, in the case of a deforming grid, don't forget that dV belongs inside the time derivative term (you may not need this if your cylinder grid is rigid).
  Reply With Quote

Old   October 3, 2006, 10:20
Default Re: Fractional step in ALE form - possible?
  #3
Wee
Guest
 
Posts: n/a
hmm... i just tried moving the whole grid to simulate the oscillating cylinder but the answer's not correct.

BC for (*) is supposed to be at t=n+1, so the grid velocity is the one at t=n+1. other than that, i don't see any other changes required.

Btw, i found that the lift/drag coefficient computed using tecplot and my own extrapolation subroutine are different, both in value and frequency, although both are wrong anyway. strangely, for sll other cases (stationary,rotating,unsteady), the cl/cd are quite similar.

i believe there's no need to modify my cl/cd formula for the . that's what i was told in another post.

  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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
OpenFOAM-1.6 install cookbook MadsR OpenFOAM Installation 372 November 20, 2010 12:57
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


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