Checknig the DF condition in the DF-SEM





Using a general function %$ q(x,y,z) $%, we may calculate the divergence in this case, which of course in a DF method must be imposed to be 0.

%BEGINLATEX{label="eq:DIVERGENZA"}% \begin{equation*}


\nabla\cdot\mathbf{u'=} & \sigma^{2}\{ & [\alpha_{3}(x-x^{k})-\alpha_{1}(z-z^{k})][+\frac{1}{r^{3}}\frac{\partial q}{\partial y}]+...\\ & & [\alpha_{2}(x-x^{k})-\alpha_{1}(y-y^{k})][-\frac{1}{r^{3}}\frac{\partial q}{\partial z}]+...\\ & & [\alpha_{3}(y-y^{k})-\alpha_{2}(z-z^{k})][-\frac{1}{r^{3}}\frac{\partial q}{\partial x}]\}\end{array}

\end{equation*} %ENDLATEX%

where %$ r = \sqrt{(x - x^k)^2+(y - y^k)^2+(z-z^k)^2} $%



\frac{\partial q}{\partial x} - \frac{3 (x - x^k) q}{ r^2} = 0 \\ \frac{\partial q}{\partial y} - \frac{3 (y - y^k) q}{ r^2} = 0 \\ \frac{\partial q}{\partial z} - \frac{3 (z - z^k) q}{ r^2} = 0

\end{cases} \end{equation*}%ENDLATEX%

The final function is: %$ q = A_0 r^3 $%. This is basically a cubic function that, when used, simplifies with the denominator in equation %reflatex{eq:DFSEM1}%, leaving the i-th velocity component independent from the i-th axis (u independent of x, v independent of y and w independent of z).

Further investigation showed that the system %reflatex{eq:SISTDIFFERENZIALE}% is not the only solution of our problem. In fact basically in equation %reflatex{eq:DIVERGENZA}% may be grouped in a different way and lead us to a different system of solving equations. Anyway, in the Matlab provided at the bottom of the page, I tried (using a symbolic calculator) to manage a %$ q(\sqrt{(x - x^k)^2+(y - y^k)^2+(z-z^k)^2}) $% general function, and for each function depending only on the distance leads us to a divergence free velocity field [see matlab script].


After all these considerations I wanted to check everything in the DF-SEM. I calculate then the divergence in any simulations and I plotted the results in fig. 2.I calculate the derivative with a central difference 2nd order scheme. This picture shows that the method is divergence free everywhere but near the eddy edge.


After all these considerations I wanted to check everything in the DF-SEM. I calculate then the divergence in any simulations and I plotted the results in fig. 2.




u = \frac{q(\mathbf{r})}{\mathbf{r}^3} (\alpha_{3}(y-y^{k})-\alpha_{2}(z-z^{k}))\\ v=\frac{q(\mathbf{r})}{\mathbf{r}^{3}}(\alpha_{1}(z-z^{k})-\alpha_{3}(x-x^{k}))\\ w=\frac{q(\mathbf{r})}{\mathbf{r}^{3}}(\alpha_{2}(x-x^{k})-\alpha_{1}(y-y^{k}))







\frac{\partial u}{\partial x}=[\frac{\partial q(\mathbf{r})}{\partial\mathbf{r}}\frac{\partial\mathbf{r}}{\partial x}\mathbf{r}^{3}-q(\mathbf{r})3\mathbf{r}^{2}\frac{\partial\mathbf{r}}{x}](\alpha_{3}(y-y^{k})-\alpha_{2}(z-z^{k}))\ \frac{\partial v}{\partial y}=[\frac{\partial q(\mathbf{r})}{\partial\mathbf{r}}\frac{\partial\mathbf{r}}{\partial y}\mathbf{r}^{3}-q(\mathbf{r})3\mathbf{r}^{2}\frac{\partial\mathbf{r}}{y}](\alpha_{1}(z-z^{k})-\alpha_{3}(x-x^{k})) \ \frac{\partial w}{\partial z}=[\frac{\partial q(\mathbf{r})}{\partial\mathbf{r}}\frac{\partial\mathbf{r}}{\partial z}\mathbf{r}^{3}-q(\mathbf{r})3\mathbf{r}^{2}\frac{\partial\mathbf{r}}{z}](\alpha_{2}(x-x^{k})-\alpha_{1}(y-y^{k}))




Current Tags:
create new tag
, view all tags
Topic attachments
I Attachment Action Size Date Who Comment
elsem symcal2cart.m manage 1.2 K 2010-04-07 - 12:42 RuggeroPoletto  
Topic revision: r10 - 2010-04-12 - 15:55:39 - RuggeroPoletto
Main Web
16 Oct 2019


Manchester CfdTm

Ongoing Projects


Previous Projects


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.