CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Wrong figure in Hrv. thesis? (https://www.cfd-online.com/Forums/openfoam-programming-development/115764-wrong-figure-hrv-thesis.html)

santiagomarquezd April 6, 2013 03:05

Wrong figure in Hrv. thesis?
 
1 Attachment(s)
Hi there, a few days ago Dr. Axel Larreteguy told me about some concerns about Figure 3.7 in Hrv. thesis (see attachment). Here we have the vectors \vec{\Delta} and \vec{S}, which have units of area and then vector \vec{k}, the correction, which seems to take in account the displacement of \vec{\Delta} from the face centroid b, so that it has units of distance. The concerns are:

1. Why the differences in units?
2. Do the vector decomposition make sense as is drawn?

Below the figure it is stated that vector \vec{k} is not used

3. The correction vector is not used because simplification or because it is zero (I guess so)?

Following the work of Mathur and Murthy

http://www.tandfonline.com/doi/abs/1...915105#preview

it leads to a zero correction. If it is the wrong interpretation

4. How do we interpret the figure?

Thanks in advance,

bigphil April 6, 2013 14:53

Quote:

Originally Posted by santiagomarquezd (Post 418701)
Hi there, a few days ago Dr. Axel Larreteguy told me about some concerns about Figure 3.7 in Hrv. thesis (see attachment). Here we have the vectors \vec{\Delta} and \vec{S}, which have units of area and then vector \vec{k}, the correction, which seems to take in account the displacement of \vec{\Delta} from the face centroid b, so that it has units of distance. The concerns are:

1. Why the differences in units?
2. Do the vector decomposition make sense as is drawn?

Below the figure it is stated that vector \vec{k} is not used

3. The correction vector is not used because simplification or because it is zero (I guess so)?

Following the work of Mathur and Murthy

http://www.tandfonline.com/doi/abs/1...915105#preview

it leads to a zero correction. If it is the wrong interpretation

4. How do we interpret the figure?

Thanks in advance,

Hi,

Not sure if I can answer your questions, but

In Hrv's thesis (and in OpenFOAM), non-orthogonal correction is not employed on the boundary. This assumption is fine if we assume that phi does not vary along the boundary.
Although, if we assume that phi varies linearly along the boundary then non-orthgonal correction is required to maintain accuracy, as outlined in
E. De Villiers thesis (Fig. 4.7).

I believe that this correction is not performed in OpenFOAM because it may lead to mass flux across walls.

But boundary non-orthogonal correction can be employed through use of custom boundary conditions, such as fixedDisplacement (corrected version of fixedValue) in the solidMechanics branch of OpenFOAM-1.6-ext.

Hope it helped,
Philip

santiagomarquezd April 7, 2013 01:37

Thanks Philip, it helped! I'm going to elaborate some ideas to post.

Regards.


All times are GMT -4. The time now is 04:37.