Entradas

An empiric profile for a turbulent flow in a circular pipe

Imagen
  Abstract Here I propose an empiric profile for the turbulent flow in a circular pipe. The profile is a simple continuous function, similar but different to the power law equation. It's only parameter is fixed based on the stress on the wall and a flat profile in the center. Equation development We all know that there is not an analytical solution for the velocity profile in a circular pipe in turbulent flow. There is, of course, for developed laminar flow:  \[v(r) = v_\text{c} \left[ 1 - \left(r \over R\right)^2 \right] \tag{1} \] where \(v_\text{c}\) is the velocity at the center, \(r\) is the radial coordinate, \(R\) the radius of the tube. On the other hand, there's a well-know empirical average profile for turbulent flow, which is know as the power-law profile: \[\bar{v}(r) = \bar{v}_\text{c} \left( 1 - {r \over R} \right)^n \tag{2} \] where \(n\) is an empirical parameter, usually around 1/7. This model represents acceptably well the average profile for the averag...

An approximation for the error function

Imagen
 The error function, erf( x ), is defined as the following function: \[ \operatorname{erf}( x) = {2 \over \pi} \int _0 ^x { e^{- t^2} \operatorname{d}t } \tag{1}\] This integral has no analytical solution through any of the conventional functions. Most modern software has a function to evaluate it. However, sometimes it is convenient to have an approximation at hand. It is known that it slowly converges as a series of powers. For this reason, the way to approximate it proposed here is through a geometrically similar function. The proposed base function is: \[y(x) = \tanh \left({{a x}\over 2}\right) = { {2 e^{ a x }} \over { 1 + e^{ a x } }} - 1 \tag{2}\] where \(a\) is a constant to be determined. This particular function is chosen because it satisfies the same function limits; that is to say, \[ \lim_{x \to \infty }{\operatorname{erf}( x)} = \lim_{x \to \infty } { y (x) } = 1 \tag{3.a}\] \[ \lim_{x \to -\infty } {\operatorname{erf}( x)} = \lim_{x \to -\infty }{ y (x) } = -1 \tag{...