CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

upwinding of "curl"-type convection

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

LinkBack Thread Tools Display Modes
Old   October 19, 1998, 10:19
Default upwinding of "curl"-type convection
Maxim Olshanskii
Posts: n/a

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"?

  Reply With Quote

Old   October 23, 1998, 13:16
Default Re: upwinding of "curl"-type convection
John C. Chien
Posts: n/a
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.
  Reply With Quote

Old   October 23, 1998, 16:05
Default Re: upwinding of "curl"-type convection
Adrin Gharakhani
Posts: n/a
>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 ( in the upwind form, which also would require you to obtain the stretch term ( 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
  Reply With Quote

Old   October 23, 1998, 20:34
Default Correction: upwinding of "curl"...
Adrin Gharakhani
Posts: n/a

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?

Adrin Gharakhani
  Reply With Quote

Old   October 24, 1998, 03:57
Default Re: Correction: upwinding of "curl"...
Maxim Olshanskii
Posts: n/a
Thank you very much for the comments. I feel that I was not qiute paticular in my question. You talk about the terms after application of \curl to the momentum equations. Indeed these results in the term like \curl(vorticity * velocity). However, I am thinking about another trick. That is consider convection u\dot\nabla u as a sum: (\curl u) * u + \nabla (u^2)/2.

The seccond term goes to pressure and the first one is what remains from convection. From some reasons this another form of convection can be more convinient. But my first attempts of simply discretizing \curl by second oder finite difference in upwind derection do not give satisfactory results...

Thank you for the reference, I'll try to find.
  Reply With Quote


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
boundary conditions for simpleFoam calculation foam_noob OpenFOAM Running, Solving & CFD 8 July 1, 2015 08:07
second order schemes marine OpenFOAM 66 October 26, 2014 12:10
Instable natural convection case Peter88 OpenFOAM 5 August 18, 2011 01:23
Pressure instability with rhoSimpleFoam daniel_mills OpenFOAM Running, Solving & CFD 44 February 17, 2011 18:08
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50

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