
[Sponsors] 
August 17, 2012, 15:29 
simultaneous conservation of angular and linear momenta

#1 
New Member
Join Date: Aug 2012
Posts: 1
Rep Power: 0 
I presently have written a code that evolves Euler's equations with Newtonian selfgravity on an Eulerian Cartesian grid. I use the KurganovTadmor central difference scheme for the advection.
I would like to be able to numerically conserve both linear and angular momenta. Does anyone know of any way to do this? I have figured out the fundamental reason the KT method does not conserve angular momentum when advecting linear momenta  it does not preserve the symmetry of the stress tensor. But I can't figure out how to alter the scheme to account for this effect. Any suggestions? 

August 23, 2012, 10:23 
conserving momenta

#2 
New Member
Patrick Godon
Join Date: Apr 2010
Posts: 19
Rep Power: 9 
Hello Peter Parker,
one way to help conserve angular momentum is to actually solve the equations for the momenta, so instead of solving for the "primitive" variables, you solve for the momenta: U = rho*vx V = rho*vy W = rho*vz where rho=density and the (vx,vy,vz) is the velocity in cartesian coordinates. A way to help conserve angular momentum is (here too) to solve for the angular momentum, so assuming for example that you have some kind of rotating flow, you used cylindrical (or spherical) coordinates and write U = rho*vr U = rho*vf W = rho*vz where the cylindrical coordinates are (r, phi, z) and the velocity is then given by (vr, vf, vz) respectively. If you have some rotating flow with "rotating boundaries" (or rotating boundary conditions) you can solve in the rotating frame of reference, so that you actually solve only for the fluctuations of the flow about the speed of rotation. For that you put in the equation vf ==> dvf + omega x r where omega is the angular velocity of the rotating system "x" is the cross product and r is the radial coordinates in the cylindrical coordinates. So you solve only for dvf, the difference between vf and omega x r . You can look in a few books for rotating references frames, Landau Lifshitz Hydrodynamics is one of them, there are more oriented one for selfgravitating rotating (celestial) objects such as the book of JeanLouis Tassoul (I think it is something like rotating stars or something like that). I hope this helps. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Multiple floating objects  CKH  OpenFOAM Running, Solving & CFD  12  March 21, 2016 14:05 
Angular momentum conservation?  Barry  Main CFD Forum  19  April 11, 2002 10:37 