PERDIDAS POR FRICCION EN UNA TUBERÍA

FRICTION LOSSES IN A PIPE

Nicolás Rojas Arias

Camilo “El toche” Fernández Cristancho

Fecha realización: marzo 14 de 2016

fecha entrega: marzo 28 de2016

RESUMEN

Se observan y analizan los cambios y pérdidas por fricción de un fluido dentro de un conducto
vertical durante la fluencia de este. Determinando a mediada que el fluido atraviesa la tubería, las
pérdidas de energía ejercidas por la fricción, disminución de presión con base al gradiente
hidráulico y la velocidad, al igual que los números de Reynolds.

Se halló que ( … )

PALABRAS CLAVE

Perdidas por fricción, gradiente hidráulico, números de Reynolds.

ABSTRACT

They are observed and analyzed the changes and friction losses on a fluid into the vertical pipe
during the fluence. Searching as it traverses the pipe, the energy losses for the friction, presasure
decrease by the hydraulic gradient and velocity, as a Reynolds´ numbers.

key words: friction losses, hydraulic gradient, Reynolds numbers.

I INTRODUCCION

Las pérdidas por fricción se dan porque al estar el fluido en movimiento realiza resistencia que se
opone a dicho movimiento convirtiendo parte de la energía en energía térmica (calor), que se disipa a
través de las paredes de la tubería. Las válvulas se encargan de controlar la dirección o el flujo
volumétrico del fluido generando turbulencia local en el fluido, ocasionando pérdida de energía
transformándose en calor. Estas pérdidas son consideradas perdidas menores ya que en un sistema
grande las pérdidas por fricción en las tuberías son mayores en comparación a la delas válvulas
y accesorios.

Las pérdidas y ganancias de energía en un sistema se contabilizan en términos de energía por unidad de peso
del fluido que circula por él. Esto también se conoce como carga (h):

ʎa  Energía que se agrega al fluido con un dispositivo mecánico; es común que se le denomine carga
total sobre la bomba.

ʎr Energía que se remueve del fluido por medio de un dispositivo mecánico.

ʎl  Pérdidas de energía del sistema por la fricción en las tuberías o pérdidas menores por válvulas y otros
accesorios.

La magnitud de las pérdidas de energía que produce la fricción del fluido, las válvulas y accesorios, es
directamente proporcional a la carga de velocidad del fluido. Esto se expresa en forma matemática así:

v2
ʎl=k (1)
2g

El término K es el coeficiente de resistencia.

ECUACION GENERAL DE LA ENERGÍA

La ecuación general de la energía es la extensión de la ecuación de Bernoulli, lo que permite

resolver problemas en los que hay pérdidas y ganancias de energía.
Para un sistema, la expresión del principio de conservación de energía es:

E 1+ʎ a−ʎ r−ʎ L=E 2 (2)

E1 y E2 denotan la energía que posee el fluido por unidad de peso en las secciones 1 y 2.

La energía que posee el fluido por unidad de peso es:

p v2
E= + + z (3)
y 2g

Es esencial que la ecuación general de la energía se escriba en la dirección del flujo.

El comportamiento de un fluido en lo que se refiere a las pérdidas de energía depende de que el

fluido sea laminar o turbulento. Un medio para preceder este comportamiento en el flujo es con
el manejo de números de Reynolds. Demostrado por Os borne Reynolds, esta ecuación se define
como

VDP VD
ℜ= = (4)
n v

Donde V es velocidad, D diámetro de la tubería P densidad del fluido y n la viscosidad cinemática

del fluido.

Este número relaciona las fuerzas de inercia sobre un elemento de fluido a la fuerza viscosa .

En la práctica se tiene en cuenta que este número puede variar de la siguiente manera.
<2000 para fluidos laminares

2000-4000 para fluidos transitorios

>4000 para fluidos turbulentos.

Teniendo en cuanta la ecuación (3) es de observar que ʎl es la perdida de energía en el sistema. Y se

expresa mediante la ecuación de DARCY:

f∗L 2
∗v
D (5)
ʎ l=
2g

f= factor de fricción; depende del régimen

64
f= (6) Para régimen laminar
ℜ

F= f =
0.25 (7) para régimen turbulento
¿¿¿

D
ES LA RUGOCIDAD RELATIVA Y E ES RUGOSIDAD PROM.
E

L= longitud de la corriente

D= diámetro de la tubería

V= velocidad promedio del flujo

Para la toma de pérdidas de fricción se dan la siguiente relación

Kf ∗v 2
ʎ l= (8)
2g

Donde Kf es el factor de perdida

II PROCEDIMIENTO EXPERIMENTAL

Para la realización de la experiencia fue requerido el uso del banco hidráulico en el cual se
procedió a la toma de medidas de las mangueras presentes, inyectando aire dentro del sistema de
tuberías, se calibra el líquido a una presión de las mismas características dentro de cada una de
las mangueras. Vaciando cada 10ml de agua dentro de una probeta, se tomó en cuenta la
medición de cada una de las alturas del fluido perteneciente en cada una de las mangueras
(ténganse en cuenta que aunque el líquido desplazado del sistema era agua, se encuentra
presente en dos de las cuatro mangueras, mercurio, como fluido secundario).

Haciendo un comparativo con respecto al flujo del fluido, se procedió a realizar los anteriores
pasos en volúmenes de: 50ml para flujo laminar y 100 ml para flujos transitorios y turbulentos.
Véase los parámetros en tabla I

TABLA I parámetros del experimento

CARACTERÍSTICA
LONGITUD DE LAS MANGUERAS 524mm
DIAMETRO DE LA TUBERÍA 3mm
COEFICIENTE DE RUGOSIDAD 0
FLUIDOS UTILIZADOS H2O ; Hg
TOMA DE DATOS POR PRUEBA 5

La toma de datos para cada uno de los tubos y experiencias tuvo que realizarse de manera eficaz,
rápida, y directa teniendo en cuenta la escala presente dentro del sistema, ya que, la variación de
altura en cada uno de los casos variaba en lapsus de tiempo bastante cortos, haciendo
indispensable la destreza del observador en el momento de la toma y recolección de datos.

Se tuvo en cuenta que cada par de tubos o columnas capilares dentro del banco se encontraba
conectado directamente en su parte inferior haciendo que a media que disminuía el fluido en un
capilar A, el acenso del fluido en B era directamente proporcional con respecto a ser contraparte.
Concluida la experiencia y tomados los datos correspondientes a cada una de pruebas, se
procedió a la tabulación de estos mostrados en las tablas II a XI.

III RESULTADOS

Aunque los datos obtenidos son satisfactorios y óptimos para poder realizar un análisis, es
requerido y, preferible ver más a fondo observar diversas variables que pueden colaborar a un
mejor entendimiento de lo que sucede. Con base en esto se toma en cuenta:

A Determinación de la velocidad.

Teniendo en cuenta la diferencia de volumen el cual se ha despejado de manera considerable, es

posible apreciar que el aumento de volumen de un tubo es directamente influido por la
disminución de su contraparte, teniendo en cuenta esto:

Δl=l 2−l1 (9)

Teniendo una distancia de desplazamiento y el tiempo requerido para el desplazamiento de del

fluido en los capilares, asumimos la velocidad cómo:

Δl
V= (10)
t

Los datos de velocidad obtenidos para cada una de las condiciones se encuentran en las tablas VII a
XI.

B Determinación de los números de Reynolds.

Los números de Reynolds, son números adimensionales los cuales pueden ser hallados a partir de
la ecuación (4) en donde se asume que:

- Densidad del agua es igual a 1000 kg/m 3

- Densidad del mercurio es igual a 13600 kg/m 3

- Viscosidad dinámica del agua 0.001

- Viscosidad dinámica del mercurio 0.00155

Las demás variables se encuentran presentes en la tabla I

VDP VD
ℜ= = Donde,
n v
n = viscosidad dinámica
p= densidad del fluido
V= velocidad
D= diámetro del tubo o capilar

Los números de Reynolds fueron registrados en sus tablas correspondientes.

C Determinación del caudal.

Dentro de la lógica matemática y el estudio de los fenómenos de transporte se aprecian dos tipos de caudal,
en relación se ve observado dentro del fluido el volumen que circula o la masa del fluido desplazado dentro
del capilar. Con base en las dos variables, las ecuaciones correspondientes serán:

m
Qm=Qv∗p= (11)
t

¥
Q v=V ∗A= (12) donde,
t

Qm = Caudal másico
Qv = Caudal volumétrico
¥ = volumen
V = velocidad

Los resultados obtenidos, fueron tabulados dentro de las tablas VII a XI.

D Determinación de la velocidad crítica

E Determinación fórmula de Darcy

Es una fórmula básica para el cálculo de pérdidas de carga en tuberías en donde teniendo en cuenta la
ecuación (5) se toma en términos de caudal la siguiente transformación.

8 fL Q 2
hf = (13) Donde,
π 2 g D2
Hf= perdida de carga
f= coeficiente de rozamiento
L= longitud
V= velocidad
g= gravedad
Q= caudal volumétrico.

F Determinación de la pendiente de carga

Teniendo en cuenta la siguiente ecuación

 V =0.8492C∗R0.63 S0.54 (14)
V= velocidad
C=coeficiente de rugosidad
R= Radio d la tubería
S= pendiente de carga

Despejando la pendiente de carga, queda:

V
S= 0.54
(√ 0.8492C∗R ) 0.63 (15)

Teniendo en cuenta las anteriores ecuaciones y datos obtenidos a partir de estas, es posible tabular los
datos en las siguientes tablas.

TABLA II. Para flujo laminar

No. V t Columna de H2O T° Re

(ml) (s) h1(mm) h2(mm) (°C)
1 2.39 31.04 311 266 20 43
2 7.27 18.72 325 242 20 133
3 12.97 13.7 338 229 20 239
4 18.52 11.64 350 217 20 343
5 35.69 11.47 359 105 20 664

TABLA III. Para flujo laminar

No. V t Columna de Hg T° Re
(ml) (s) h1(mm) h2(mm) (°C)
1 0.04 31.04 179 175 20 34
2 0.06 18.72 180 174 20 84
3 0.09 13.7 183 173 20 192
4 0.11 11.64 184 172 20 382
5 0.13 11.47 185 171 20 321

TABLA IV. Para flujo transitorio

No. V t Columna de H2O T° Re

(ml) (s) h1(mm) h2(mm) (°C)
1 25.86 13.9 390 163 20 490
2 32.93 12.59 405 140 20 631
 3 34.32 13.51 422 122 20 666
4 43.21 12.12 443 100 20 849
5 44.92 12.59 460 85 20 894

TABLA V. Para flujo transitorio

No. V t Columna de Hg T° Re
(ml) (s) h1(mm) h2(mm) (°C)
1 0.16 13.9 185 168 20 322
2 0.20 12.59 187 166 20 439
3 0.25 13.51 180 164 20 507
4 0.28 12.12 192 162 20 652
5 0.31 12.59 193 160 20 690

TABLA VI. Para flujo turbulento

No. V t Columna de Hg T° Re
(ml) (s) h1(mm) h2(mm) (°C)
1 0.53 8.32 205 149 20 1772
2 0.58 8.19 210 148 20 1993
3 0.64 6.96 213 145 20 2572
4 0.70 5.44 216 142 20 3581
5 0.77 5.5 220 138 20 3924

TABLA VII. Para flujo laminar

No. V t Q V H2O hf f F hf Sf1 Re

(ml) (s) (m3/s) (m/s) h1(m) h2(m) Darcy Caudal
1 2.39 31.04 1.4 E-5 0.001 311 266 0.045 15 1.47 0.0003 1 43
2 7.27 18.72 4.2 E-5 0.004 325 242 0.083 5 0.48 0.0008 8 133
3 12.97 13.7 7.5 E-5 0.008 338 229 0.109 3 0.27 0.0015 25 239
4 18.52 11.64 1.1 E-4 0.011 350 217 0.133 2 0.19 0.0022 48 343
5 35.69 11.47 2.1 E-4 0.022 359 105 0.254 1 0.10 0.0042 16 664
3

TABLA VIII. Para flujo laminar

No. V t Q V Hg hf f F hf Sf1 Re
(ml) (s) (m3/s) (m/s) h1(m) h2(m) Darcy Caudal
1 0.04 31.04 1.2 E-6 0.0001 179 175 0.004 71 1.89 0.0000 0 34
2 0.06 18.72 3.0 E-6 0.0003 180 174 0.006 29 0.76 0.0000 0 84
3 0.09 13.7 6.9 E-6 0.0007 183 173 0.01 13 0.33 0.00001 0 198
4 0.11 11.64 9.7 E-6 0.0010 184 172 0.012 9 0.24 0.00001 1 271
5 0.13 11.47 1.2 E-5 0.0012 185 171 0.014 8 0.20 0.00001 1 321

TABLA IX. Para flujo transitorio

No. V t Q V H2O hf f F hf Sf1 Re

(ml) (s) (m3/s) (m/s) h1(m) h2(m) Darcy Caudal
1 25.86 13.9 1.5 E-4 0.016 390 163 0.23 0.09 0.10 0.0002 93 490
2 32.93 12.59 2.0 E-4 0.021 405 140 0.27 0.08 0.12 0.0003 14 631
9
3 34.32 13.51 2.1 E-4 0.22 422 122 0.3 0.08 0.12 0.0003 16 666
4
4 43.21 12.12 2.7 E-4 0.28 443 100 0.34 0.07 0.14 0.0005 25 849
7
5 44.92 12.59 2.8 E-4 0.30 460 85 0.38 0.07 0.15 0.0005 28 894
3

TABLA X. Para flujo transitorio

No. V t Q V Hg hf f F hf Sf Re
(ml) (s) (m3/s) (m/s) h1(m) h2(m) Darcy Caudal
1 0.16 13.9 1.2 E-5 0.001 185 168 0.02 0.11 0.03 0 1 322
2 0.20 12.59 1.6 E-5 0.002 187 166 0.02 0.10 0.03 0 1 439
3 0.25 13.51 1.8 E-5 0.002 190 164 0.03 0.09 0.04 0 2 507
4 0.28 12.12 2.3 E-5 0.002 192 162 0.03 0.08 0.04 0 3 652
5 0.31 12.59 2.5 E-5 0.003 193 160 0.03 0.08 0.04 0 3 690

TABLA XI. Para flujo turbulento

No. V t Q V Hg hf f F hf Sf Re
(ml) (s) (m3/s) (m/s) h1(m) h2(m) Darcy Caudal
1 0.53 8.32 6.3 E-5 0.007 205 149 0.06 0.05 0.08 0 18 1772
2 0.58 8.19 7.1 E-5 0.008 210 148 0.06 0.05 0.09 0 22 1993
3 0.64 6.96 9.2 E-5 0.10 213 145 0.07 0.05 0.11 0 16 572
4 0.70 5.44 1.3 E-4 0.014 216 142 0.07 0.04 0.15 0 66 3581
5 0.77 5.5 1.4 E-4 0.015 220 138 0.08 0.04 0.16 0 79 3924

IV ANÁLISIS DE RESULTADOS
 V CONCLUSIONES

REFERENCIAS.

NUMEROS DE REYNOLDS

REYNOLDS´ NUMBERS

Nicolás Rojas Arias

Camilo “ El toche” Fernández Cristancho

Fecha realización:
fecha entrega:

RESUMEN

Los fluidos pueden comportarse de distintas maneras teniendo en cuenta si el fluido es

newtoniano o no. Además, la influencia de esto, más la importancia del sitio y la velocidad por la
que se mueve el fluido, demuestra la particularidad de cómo se puede comportar el fluido en
dichas circunstancias, es decir, si este se mueve de forma laminar o torna a un flujo turbulento.
Para esto se analizan los cambios del fluido dentro de un capilar, las pérdidas por fricción,
disminución de presión y los números de Reynolds.
Se encontró que (…)

Palabras clave: números de Reynolds. Flujo, presión.

ABSTRACT

Fluids could play of different forms, it observing if the fluid is Newtonian or isn´t. so, the
influences, plus the importance of the site and the fluid´s velocity show how can play this in these
places, that is to say, if this is moving I laminar or turbulent flow. For this are analyzed the friction
losses and changes of the fluid into the capillary the pression decrease and the Reynolds´
numbers.

Key words: Reynolds numbers, flow, pression.

I INTRODUCCION

El comportamiento de un fluido en lo que se refiere a las pérdidas de energía depende de que el

fluido sea laminar o turbulento. Un medio para preceder este comportamiento en el flujo es con
el manejo de números de Reynolds. Demostrado por Os borne Reynolds, esta ecuación se define
como

VDP VD
ℜ= = (1)
n v

Donde V es velocidad, D diámetro de la tubería P densidad del fluido y n la viscosidad cinemática

del fluido.

Este número relaciona las fuerzas de inercia sobre un elemento de fluido a la fuerza viscosa.

En la práctica se tiene en cuenta que este número puede variar de la siguiente manera.

<2000 para fluidos laminares

2000-4000 para fluidos transitorios

>4000 para fluidos turbulentos.

Las bases para la comprensión del por qué se decide si un fluido puede ser laminar, transitorio o
turbulento según sea el caso puede verse como variable fundamental la velocidad a la que se
encuentra sometido el fluido, además esta no solo influye en el tipo de flujo presente en el ensayo
sino en el aumento o disminución del caudal presente al pasar a través de un tubo.
Para el análisis y control de esto se han creado dos ecuaciones las cuales pueden representar de
manera numérica el caudal de una manera ya sea volumétrica o másica.

m
Qm=Qv∗p= (2)
t

¥
Q v=V ∗A= (3) donde,
t

Qm = Caudal másico
Qv = Caudal volumétrico
¥ = volumen
V = velocidad

II PROCEDIMIENTO EXPERIMENTAL

Para la realización de la experiencia fue requerido el uso del banco hidráulico en el cual se
procedió a la toma de medidas de las mangueras presentes, inyectando aire dentro del sistema de
tuberías, se calibra el líquido a una presión de las mismas características dentro de cada una de
las mangueras. Vaciando cada 10ml de agua dentro de una probeta, se tomó en cuenta la
medición de cada una de las alturas del fluido perteneciente en cada una de las mangueras
(ténganse en cuenta que aunque el líquido desplazado del sistema era agua, se encuentra
presente en dos de las cuatro mangueras, mercurio, como fluido secundario).

Haciendo un comparativo con respecto al flujo del fluido, se procedió a realizar los anteriores
pasos en volúmenes de: 50ml para flujo laminar y 100 ml para flujos transitorios y turbulentos.
Véase los parámetros en tabla I

TABLA I parámetros del experimento

CARACTERÍSTICA
LONGITUD DE LAS MANGUERAS 524mm
DIAMETRO DE LA TUBERÍA 3mm
COEFICIENTE DE RUGOSIDAD 0
FLUIDOS UTILIZADOS H2O ; Hg
TOMA DE DATOS POR PRUEBA 5

La toma de datos para cada uno de los tubos y experiencias tuvo que realizarse de manera eficaz,
rápida, y directa teniendo en cuenta la escala presente dentro del sistema, ya que, la variación de
altura en cada uno de los casos variaba en lapsus de tiempo bastante cortos, haciendo
indispensable la destreza del observador en el momento de la toma y recolección de datos.

Se tuvo en cuenta que cada par de tubos o columnas capilares dentro del banco se encontraba
conectado directamente en su parte inferior haciendo que a media que disminuía el fluido en un
capilar A, el acenso del fluido en B era directamente proporcional con respecto a ser contraparte.
Concluida la experiencia y tomados los datos correspondientes a cada una de pruebas, se
procedió a la tabulación de estos mostrados en las tablas VII a XI.

Aunque los datos obtenidos son satisfactorios y óptimos para poder realizar un análisis, es
requerido y, preferible ver más a fondo observar diversas variables que pueden colaborar a un
mejor entendimiento de lo que sucede. Con base en esto se toma en cuenta:

A Determinación de la velocidad. Teniendo en cuenta la diferencia de volumen el cual se ha

despejado de manera considerable, es posible apreciar que el aumento de volumen de un tubo es
directamente influido por la disminución de su contraparte, teniendo en cuenta esto:

Δl=l 2−l1 (4)

Teniendo una distancia de desplazamiento y el tiempo requerido para el desplazamiento de del

fluido en los capilares, asumimos la velocidad cómo:

Δl
V= (5)
t

Los datos de velocidad obtenidos para cada una de las condiciones se encuentran en las tablas II a
VI

B Determinación de los números de Reynolds.

Los números de Reynolds, son números adimensionales los cuales pueden ser hallados a partir de
la ecuación (1) en donde se asume que:

- Densidad del agua es igual a 1000 kg/m 3

- Densidad del mercurio es igual a 13600 kg/m 3

- Viscosidad dinámica del agua 0.001

- Viscosidad dinámica del mercurio 0.00155

Las demás variables se encuentran presentes en la tabla I

VDP VD
ℜ= = Donde,
n v

n = viscosidad dinámica
p= densidad del fluido
V= velocidad
D= diámetro del tubo o capilar

Los números de Reynolds fueron registrados en sus tablas correspondientes.

C Determinación del caudal.

Dentro de la lógica matemática y el estudio de los fenómenos de transporte se aprecian dos tipos de caudal,
en relación se ve observado dentro del fluido el volumen que circula o la masa del fluido desplazado dentro
del capilar. Con base en las dos variables, las ecuaciones correspondientes serán (2) y (3) respectivamente.

m
Qm=Qv∗p=
t

¥
Q v=V ∗A=
t

Qm = Caudal másico
Qv = Caudal volumétrico
¥ = volumen
V = velocidad

Los resultados obtenidos, fueron tabulados dentro de las tablas II a VI

Concluidos los experimentos dentro del banco hidráulico, y el análisis matemático

correspondiente a cada una de las actividades, se continuó a registras los datos dentro de las
siguientes dos tablas.

TABLA II. Para H2O en fluido laminar.

V(ml) T(s) Re Q(m3/s) Qm(kg/s) V(m/s) Crítica

2.39 31.04 43 1.4 E-5 0.014 0.001
2.27 18.72 133 4.2 E-5 0.042 0.004
12.97 13.7 239 7.5 E-5 0.075 0.008
18.52 11.64 343 1.1 E-4 0.108 0.011
35.69 11.47 664 2.1 E-4 0.209 0.022

TABLA III. Para Hg en fluido laminar.

V(ml) T(s) Re Q(m3/s) Qm(kg/s) V(m/s) Crítica

0.04 31.04 34 1.2 E-6 0.017 0.0001
0.06 18.72 84 3.0 E-6 0.041 0.0003
0.09 13.7 198 6.9 E-6 0.094 0.0007
0.11 11.64 271 9.7 E-6 0.132 0.0010
0.13 11.47 321 1.2 E-5 0.156 0.0012

TABLA IV. Para H2O en fluido transitorio.

 V(ml) T(s) Re Q(m3/s) Qm(kg/s) V(m/s) Crítica
25.86 13.9 490 1.5 E-4 0.15 0.016
32.93 12.59 631 2.0 E-4 0.20 0.021
34.32 13.51 666 2.1 E-4 0.21 0.22
43.21 12.12 849 2.7 E-4 0.27 0.28
44.92 12.59 894 2.8 E-4 0.28 0.30

TABLA V. Para Hg en fluido transitorio.

V(ml) T(s) Re Q(m3/s) Qm(kg/s) V(m/s) Crítica

0.13 13.9 322 1.2 E-5 0.16 0.001
0.20 12.59 439 1.6 E-5 0.21 0.002
0.25 13.51 507 1.8 E-5 0.25 0.002
0.28 12.12 652 2.3 E-5 0.32 0.002
0.31 12.59 690 2.5 E-5 0.34 0.003

TABLA VI. Para Hg en fluido turbulento.

V(ml) T(s) Re Q(m3/s) Qm(kg/s) V(m/s) Crítica

0.53 8.32 1772 6.3 E-5 0.86 0.007
0.58 8.19 1993 7.1 E-5 0.97 0.008
0.64 6.96 572 9.2 E-5 1.25 0.10
0.70 5.44 3581 1.3 E-4 0.74 0.014
0.77 5.50 3924 1.4 E-4 0.91 0.015

III ANÁLISIS DE RESULTADOS

IV CONCLUSIONES

REFERENCIAS

 [1] Xilong Zhang, Yichun Wang, Peng Cang, Rui Wang, “Experimental investigation of thermal
hydraulic performance of heat exchangers with different Reynolds numbers on both air-side and
water-side”, Applied Thermal Engineering, January 2016

 [2] I. Rodríguez, O. Lehmkuhl, J. Chiva, R. Borrell, A. Oliva, “On the flow past a circular
cylinder from critical to super-critical Reynolds numbers: Wake topology and vortex shedding”,
International Journal of Heat and Fluid Flow, Volume 55,  Pages 91-103, October 2015.

 [3] Seong Mo Yeon, Jianming Yang, Frederick Stern,” Large-eddy simulation of the flow past
a circular cylinder at sub- to super-critical Reynolds numbers”, Applied Ocean Research,  Available
online 23 December 2015
 [4] ZhenZhen Tao, Tong Zhang,”Stability and convergence of two-level iterative methods for
the stationary incompressible magnetohydrodynamics with differentReynolds numbers”, Journal of
Mathematical Analysis and Applications,  Volume 428, Issue 1,  Pages 627-652, 1 August 2015.