CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   upwinding of "curl"-type convection (http://www.cfd-online.com/Forums/main/261-upwinding-curl-type-convection.html)

 Maxim Olshanskii October 19, 1998 10:19

upwinding of "curl"-type convection

Hi,

Thank you for any advice and/or reference on upwinding of convection written in "curl" form: \curl u \times u, where \times stands for vector product.

Generaly, what is the essential reading and references on the topic of "upwinding"?

Thanks

 John C. Chien October 23, 1998 13:16

Re: upwinding of "curl"-type convection

In the finite-difference approximation of the first order derivative term, such as u * ( df/dx ), there are several ways to approximate it. The central difference form is: u(i,j,k) * ( (f(i+1,j,k)-f(i-1,j,k))/ (x(i+1,j,k)-x(i-1,j,k)) ). The other form is: IF ( u(i,j,k) is Positive ) , u(i,j,k) * ( (f(i,j,k)-f(i-1,j,k))/ (x(i,j,k)-x(i-1,j,k)) ) ,AND IF ( u(i,j,k) is NEGATIVE ) , u(i,j,k) * ( (f(i+1,j,k)-f(i,j,k))/ (x(i+1,j,k)-x(i,j,k)) ). This is the so-called 1st order up-wind one-sided difference form. It is less accurate than the central difference form, but it is more stable when used with the second-order diffusion term in the Navier-Stokes equations. In general, the variable f can be anything, and u does not have to be the physical velocity at all.

 Adrin Gharakhani October 23, 1998 16:05

Re: upwinding of "curl"-type convection

>Thank you for any advice and/or reference on upwinding of >convection written in "curl" form: \curl u
>\times u, where \times stands for vector product.

Basically, this term is equal to vorticity^velocity, where ^ implies vector product as you mention above. However, this term is not just convection even though it does show up subsequent to taking the curl of the "velocity convection" term! If you expand it out, you'll see that it is equal to the convection of vorticity+stretch of vorticity+ (2 other terms having to do with the divergence of vorticity and velocity). Then the question becomes whether you wish to obtain the convection of vorticity (u.dot.del(vorticity)) in the upwind form, which also would require you to obtain the stretch term (vorticity.dot.del(velocity)) in a consitent manner.

Or alternatively, you can just difference vorticity^velocity but you cannot call this a convection term anymore (it includes at least 2 physically different processes in it!)

For a consistent finite differencing of the above term, check a paper by Weinan E (the last name is just that one letter E) in Journal of Computational Physics in the 90's (I don't remember the exact date or title). As keyword for your search you can use velicity (as in velocity-vorticity), or vorticity (and perhaps impulse and magnetization)

Hope it helps

 Adrin Gharakhani October 23, 1998 20:34

Correction: upwinding of "curl"...

Correction!

I was not clear (and was thus incorrect) in my terminology in my previous message.

I meant to say curl(vorticity^velocity) in my explanation and not just vorticity^velocity! The former results from taking the curl of the convection term. The latter.... what does it mean?