rotationVelocity in rotatingTotalPressure and SRF

deepblue17 · July 18, 2009, 06:20

Hello,

I think the computation of the rotationVelocity in rotatingTotalPressure and SRF could be simplified. The piece of code is:

Code:

vector axisHat = omega_/mag(omega_);
vectorField rotationVelocity =
omega_ ^ (patch().Cf() - axisHat*(axisHat & patch().Cf()));

In mathematical formulation this yields to:

$\vec{axis} = \frac{\vec{\omega} }{ |\vec{\omega}|}$
$\vec{V} = \vec{\omega} \times ( \vec{r} - \vec{axis} \bullet (\vec{axis} \bullet \vec{r} ) ) = \vec{\omega} \times \vec{r} - \vec{\omega} \times \vec{axis} \bullet (\vec{axis} \bullet \vec{r} )$

Reformulation of the right term:
$\vec{\omega} \times \vec{axis} \bullet (\vec{axis} \bullet \vec{r} ) = ( \vec{\omega} \times \frac{ \vec{\omega} }{ |\vec{\omega}|} ) \bullet (\vec{axis} \bullet \vec{r} ) = ( \vec{\omega} \times \vec{\omega} ) \bullet \frac{ (\vec{axis} \bullet \vec{r} )}{ |\vec{\omega}|} = \vec{0}$

This reduces the computation of the rotationVelocity to:

$\vec{V} = \vec{\omega} \times \vec{r}$

or in C++

Code:

vectorField rotationVelocity = omega_ ^ patch().Cf();

immortality · December 7, 2012, 08:11

can we use rotatingTotalPressure in a 2-D case?

July 18, 2009, 06:20	rotationVelocity in rotatingTotalPressure and SRF	#1
deepblue17 Member Oliver Borm Join Date: Mar 2009 Posts: 60 Rep Power: 17	Hello, I think the computation of the rotationVelocity in rotatingTotalPressure and SRF could be simplified. The piece of code is: Code: vector axisHat = omega_/mag(omega_); vectorField rotationVelocity = omega_ ^ (patch().Cf() - axisHat*(axisHat & patch().Cf())); In mathematical formulation this yields to: $\vec{axis} = \frac{\vec{\omega} }{ \|\vec{\omega}\|}$ $\vec{V} = \vec{\omega} \times ( \vec{r} - \vec{axis} \bullet (\vec{axis} \bullet \vec{r} ) ) = \vec{\omega} \times \vec{r} - \vec{\omega} \times \vec{axis} \bullet (\vec{axis} \bullet \vec{r} )$ Reformulation of the right term: $\vec{\omega} \times \vec{axis} \bullet (\vec{axis} \bullet \vec{r} ) = ( \vec{\omega} \times \frac{ \vec{\omega} }{ \|\vec{\omega}\|} ) \bullet (\vec{axis} \bullet \vec{r} ) = ( \vec{\omega} \times \vec{\omega} ) \bullet \frac{ (\vec{axis} \bullet \vec{r} )}{ \|\vec{\omega}\|} = \vec{0}$ This reduces the computation of the rotationVelocity to: $\vec{V} = \vec{\omega} \times \vec{r}$ or in C++ Code: vectorField rotationVelocity = omega_ ^ patch().Cf();

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
SRF and interFoam	rcastilla	OpenFOAM Running, Solving & CFD	3	November 9, 2018 09:56
Rotating Reference Frame	paullees	OpenFOAM Running, Solving & CFD	14	May 3, 2018 00:24
Rotating objects time dependend SRF	naish	OpenFOAM Running, Solving & CFD	0	January 24, 2009 03:55
Periodic boundary and SRF	Alex H	FLUENT	3	September 18, 2006 12:25
Strange things in SRF (urgent)	Cem	FLUENT	0	December 19, 2005 10:37

December 7, 2012, 08:11		#2
immortality Senior Member Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,208 Rep Power: 26	can we use rotatingTotalPressure in a 2-D case?