From vorticity to the velocity field

 <br /><br />\begin{equation*} \boldsymbol{\omega} = \sqrt{\frac{V_b k}{N \sigma^3}} \sum_{k=1}^{N} \boldsymbol{\alpha}^k g_{\sigma}(\mathbf{x}-\mathbf{x}^k)\end{equation*}<br /><br />(1)


This first equation is basicly the same we used in the SEM method, applied to the vorticity field instead of the velocity field. What we need now to calculate the velocity field from the generated vorticity.

 <br /><br />\begin{equation*} \nabla \times \boldsymbol{\omega} = \nabla (\nabla \cdot \mathbf{u}) - \nabla^2 \mathbf{u} \end{equation*}<br /><br />(2)


Using our divergence free hypotesis ($ \nabla \cdot \mathbf{u} = 0 $) we are allowed to delete the first term in the right part in equation (2), and we get:

 <br /><br />\begin{equation*} \nabla^2 \mathbf{u} = - \nabla \times \boldsymbol{\omega} \end{equation*}<br /><br />(3)


This is the equation we have to solve to get our velocity field. To solve it we can use the Green's function and the Biot-Savart law that relate circulation and velocity. All together we get as the final result:

 <br /><br />\begin{equation*} \mathbf{u}'(\mathbf{x}) = - \frac{1}{4 \pi} \int \frac{\mathbf{x - y}}{|\mathbf{x - y}|^3} \times \boldsymbol{\omega}(\mathbf{y})~d\mathbf{y} \end{equation*}<br /><br />(4)


Using (1) this formula, after some passages, can be simplified into:

 <br /><br />\begin{equation*} \mathbf{u}'(\mathbf{x}) = - \frac{1}{4 \pi} \int \frac{\mathbf{x - y}}{|\mathbf{x - y}|^3} \times \sqrt{\frac{V_b k}{N \sigma^3}} \sum_{k=1}^{N} \boldsymbol{\alpha}^k g_{\sigma}(\mathbf{y}-\mathbf{y}^k)~d\mathbf{y} \end{equation*}<br /><br />(5)

 <br /><br />\begin{equation*} \mathbf{u'}(\mathbf{x}) = \sqrt{\frac{V_b k}{N \sigma^3}} \sum_{k=1}^{N} \mathbf{K}_{\sigma}(\frac{\mathbf{x}-\mathbf{x}^k}{\sigma}) \times \boldsymbol{\alpha}^k \end{equation*}<br /><br />(6)

SEMDFC - SEM with Divergence Free Constraint

The main idea is to develop the old SEM method, forcing one component to satisfy the divergence free condition. To pursue such a scope, we decided to force the 3rd component to be Divergece Free, even if such a constraint implies loss of part of the Reynolds Stress information.

The first problem I dealt with is the stability one. In fact, since we are basically solving $ \nabla \cdot \mathbf{u} = 0 $ to calculate the third component, and since this is a differential equation, we have to pay a lot of care about it.

stability_problem.jpg The first algorithm used is the following one (note the pedix N,S,E,W,T,B have the meaning of the neighbours point):

   \begin{equation*} \frac{\Delta w}{\Delta z} = - (\frac{\Delta u}{\Delta x} + \frac{\Delta v}{\Delta y}) \end{equation*}  (7)

   \begin{equation*} \frac{\Delta w}{\Delta z} {|}_P = - (\frac{U_T - U_B}{2 \Delta x} + \frac{V_N - V_S}{2 \Delta y})~~~~~~ \end{equation*}  (8)

   \begin{equation*} W_P = W_E + \frac{\Delta w}{\Delta z} {|}_P ~ \Delta z \end{equation*}  (9)

As we can see from this plot, this method is totally divergent because the error can propagate with the calculation while solving equation (9). So I slightly modified this equation in a "back-solving" one:


   \begin{equation*} W_P = W_W - \frac{\Delta w}{\Delta z} {|}_P ~ \Delta z ~~~~~\end{equation*}  (10)

instabvel.jpgAs we can see from the following graph, this new algorithm did not solve the stability problem. In fact the new error propagation now is backward.

I tried then to run a simulation considering all the 27 sourroungind points (basicalli all the points forming a cube around the point P), but with mainly the same results.

Another problem that i saw is that the instantaneous velocity calculated has a shape totally different from the other two random generated components.

istant_vel.jpg


Current Tags:
create new tag
, view all tags
Topic revision: r7 - 2010-01-25 - 17:07:02 - RuggeroPoletto
Main Web
16 Dec 2017

Site

Manchester CfdTm
Code_Saturne

Ongoing Projects

ATAAC
KNOO

Previous Projects

DESider
FLOMANIA

Useful Links:

User Directory
Photo Wall
Upcoming Events
Add Event
 

Computational Fluid Dynamics and Turbulence Mechanics
@ the University of Manchester
Copyright © by the contributing authors. Unless noted otherwise, all material on this web site is the property of the contributing authors.