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

Momentum Source in an accelerating frame

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 24, 2006, 03:22
Default Momentum Source in an accelerating frame
  #1
fakir
Guest
 
Posts: n/a
Hi

i am trying to add a momentum term to account for an accelerating (sinusoidal translation) reference frame

I have added a sub domain as follows

FLOW: DOMAIN: Ship

SUBDOMAIN: Body Forces Translation

Coord Frame = local

Location = Assembly

BULK SOURCES:

Option = Use Volume Fraction

MOMENTUM SOURCE:

GENERAL MOMENTUM SOURCE:

Momentum Source X Component = - x_accel * density

Momentum Source Y Component = - y_accel * density

Momentum Source Z Component = - z_accel * density

Option = Cartesian Components

END

END

END

END END END

but this doesn't seem to change the simulation at all, even when the amplitude of the accelerations is very large. Two questions

1) Is my CFX set-up correct? 2) Is my physics correct?

Any advice on the advice would be appreciated
  Reply With Quote

Old   November 24, 2006, 10:09
Default Re: Momentum Source in an accelerating frame
  #2
Robin
Guest
 
Posts: n/a
Hi Fakir,

If there is no density variation, the acceleration is purely linear (i.e. no rotation) and you have a closed control volume without pressure boundaries, these will have no effect because all the forces balance out. If your accelerations account for rotation, then you might see something.

If it is purely acceleration you are trying to add, I would do so by using adding gravitional acceleration in the domain object instead.

Regards, Robin
  Reply With Quote

Old   November 24, 2006, 19:04
Default Re: Momentum Source in an accelerating frame
  #3
fakir
Guest
 
Posts: n/a
Hi,

I am adding both translation and rotation.

The rotation is calculated by a similar momentum source addition, as well as a rotation of the gravity vectors.

The density is non-constant, as it is a multi-phase flow

I may just have to try setting one of the momentum source terms to a constant, and see if that has an effect

Thanks for the advice
  Reply With Quote

Old   November 27, 2006, 11:43
Default Re: Momentum Source in an accelerating frame
  #4
Robin
Guest
 
Posts: n/a
Hi Fakir,

Are you adding momentum sources AND gravity? This would not be advisable and you may be cancelling your source terms out. You should be able to do the entire thing by adding the appropriate acceleration vectors in the gravity section. It is perhaps unfortunate that CFX calls this gravity, as it is just acceleration.

-Robin
  Reply With Quote

Old   November 29, 2006, 03:17
Default Re: Momentum Source in an accelerating frame
  #5
fakir
Guest
 
Posts: n/a
Hi Robin,

That is an excellent suggestion, and one I had thought of, but possibly not applicable just here. Let me digress a bit.

There are basically two interacting mechanisms here

1) Accelerations set up due to the translational acceleration. This can easily be included in the gravity terms

2) Acceleration set up due to rotational acceleration, which I can't so easily incorporate.

(There is also gravity, but I guess you could include that in the above)

The basically formula for (2) is well known and is given by

a = d^2r/dt^2-2*(omega cross v)-omega cross (omega cross r)- (d(omega)dt cross r)

Basically I have evaluated this component wise in order to get a matrix [a_x,a_y_a_z]

For example a_x = pitch_accel * z - yaw_accel * y + pitch_rate * (roll_rate * y - pitch_rate * x) - yaw_rate+(-roll_rate*z+yaw_rate*x)+2*(pitch_rate*Water.w - yaw_rate*Water.v)

assuming the physics and maths is correct!

There are similar a_y and a_z

Then for the rotation for we add a momentum source in the x component "-density * (Xx * ax + Yx * ay + Zx * az)" where Xx, Yx, Zx are rotation components (eg Zx = sin(Pitch) * cos(Roll) )

If the expression Xx * ax + Yx * ay + Zx * az is attempted to be added to gravity this is flagged as an error, as ax contains x,y,z,u,w, which cannot be used in gravity expressions.

I hope there is an easier way to calculate these!!
  Reply With Quote

Old   November 30, 2006, 02:52
Default Re: Momentum Source in an accelerating frame
  #6
fakir
Guest
 
Posts: n/a
Here is the (possible) source of the problem

If I define the momentum source in the x direction as

-density * (a * omega^2 * sin(omega*t))

I get (possibly) correct results.

If I define an expression 'Surge Accel = a * omega^2 * sin(omega*t)' then set the momentum source to -density * (Surge Accel)

I get nothing. This doesn't seem to make sense! In the end the accel's will be read in from a file, so I may just try that.
  Reply With Quote

Old   December 1, 2006, 20:13
Default Re: Momentum Source in an accelerating frame
  #7
Robin
Guest
 
Posts: n/a
Both can be included in the gravity box. Just decompose your centrifugal acceleration into the x, y and z components and add them to your translational acceleration.

An alternative to all of this is to use the moving mesh and specify the actual motion of your domain in a subdomain over your entire mesh. The mesh will end up moving as a solid body and the fluid will react accordingly. Nice thing about this method is you will see the whole thing moving in Post. Just remember to specify the mesh motion at your boundaries with either the same equations as the subdomain, or as "unspecified".

Regards, Robin
  Reply With Quote

Old   December 2, 2006, 05:00
Default Re: Momentum Source in an accelerating frame
  #8
fakir
Guest
 
Posts: n/a
Hi Robin,

Thanks again for the help. I will give it another go, but CFX wouldn't let me add the rotational terms to gravity as they are functions of x,y,z (the 'r' vector'), and the liquid velocity (the 'v' vector).

As an aside I tried reading the whole thing in from a text file, and got the type of motion I expected (albeit larger amplitudes then I thought, but that may be a whole other story)
  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
swak4foam building problem GGerber OpenFOAM Installation 54 April 24, 2015 16:02
Simulation of Axial Fan Flow using A Momentum Source Subdomain Liam CFX 28 July 16, 2013 08:24
pisoFoam compiling error with OF 1.7.1 on MAC OSX Greg Givogue OpenFOAM Programming & Development 3 March 4, 2011 18:18
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51


All times are GMT -4. The time now is 13:29.