MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM

19 al 21 DE SEPTIEMBRE DE 2018CAMPECHE, CAMPECHE, MÉXICO

Tema A4 Termofluidos: Mecánica respiratoria.

"Modelo matemático de parámetros concentrados para describir la histéresis

hidrodinámica de las curvas rinomanométricas"

Pérez Mota J. E., Solorio Ordaz F. J., Cervates de Gortari J. G.

Departamento de Termofluidos, DIMEI, Facultad de Ingeniería
UNAM, Av. Universidad no. 3000, Ciudad Universitaria, C.P. 04510, CDMX, México
[email protected], [email protected], [email protected]

RESUMEN

En la otorrinolaringología es importante medir la capacidad de respirar, para este propósito se emplea la

rinomanometría, prueba en la que, como función del tiempo, se miden el flujo volumétrico y la diferencia de presión.
Cuando se grafican estas dos variables aparecen bucles que denotan la existencia de un fenómeno histerético asociado
con el flujo. El propósito de esta investigación es brindar una explicación sencilla acerca de la aparición del fenómeno
de histéresis hidrodinámica en las curvas rinomanométricas. Para ello se resolvió la ecuación de balance de momentum
y se tomó en cuenta la compresibilidad del fluido en un volumen de control con dos fronteras y un flujo oscilatorio. Se
obtuvo una ecuación que deja ver que la histéresis hidrodinámica es producto de la aceleración de un flujo pulsátil y la
compresibilidad del aire. Estos resultados podrían ser útiles para rediseñar los rinomanómetros y anticipar los
resultados de una turbinectomía.
Palabras Clave: Histéresis hidrodinámica, Rinomanometría, Modelo de parámetros concentrados, Balance de momentum.

ABSTRACT

To measure the ability to breathe is important in the otorhinolaryngology, for this purpose, the rhinomanometry is used.
In this test, as a function of time, the flow rate and the pressure difference are measured; Loops appear when both
variables are plotted, which means that exist a hysteretic phenomenon associated with the flow. The purpose of this
research was to provide a simple explanation about the appearance of the hydrodynamic hysteresis phenomenon in the
rhinomanometric plots. For this, in a control volume with two boundaries and an oscillatory flow, the compressibility of
the air was taken into account and the momentum balance equation was solved. As a result, it was obtained an equation
which shows that the hydrodynamic hysteresis is a consequence of the air compressibility and the acceleration in a
pulsatile flow. These findings could be useful to redesign the rhinomanometers and to predict the results of a
turbinectomy.
Keywords: Hydrodynamic hysteresis, Rhinomanometry, Concentrated parameters model, Momentum balance.

Nomenclatura t Tiempo (s)

p A , p B Presión manométrica estática en las fronteras A y B (Pa)
A Área transversal promedio en el volumen de control (m2)
AA , AB Área de la sección transversal en las fronteras A y B (m2) po Presión de referencia (Pa)

DA , DB Diámetro hidráulico en A y B (m)

p Diferencia de presión (Pa)

HA Histéresis hidrodinámica por acumulación de masa (Pa)

pM Diferencia de presión máxima (Pa)

HI Histéresis hidrodinámica por efectos inerciales (Pa)

, Constantes empíricas

Tr Periodo de tiempo de un ciclo respiratorio (s)

A, B Coeficientes de flujo en A y en B
 Densidad media en el volumen de control (kgm-3)
U Velocidad media en el volumen de control (ms-1)
U A ,UB Velocidad media en las fronteras A y B (ms-1)
 A , B Densidad media en las fronteras A y B (kgm-3)

V Volumen del sistema abierto (m-3) o Densidad de referencia (kgm-3)

VB VB Flujo volumétrico en las fronteras A y B (m3s-1)  Tiempo normalizado

,
A , B Coeficientes en A y en B
VBM Flujo volumétrico máximo en B (m3s-1)
Fx Vector de fuerza en la dirección x (N)
ao Velocidad del sonido (ms-1)
n Vector normal unitario
k Relación de calores específicos u Vector de velocidad (ms-1)
n Número entero y par ux Vector de velocidad en la dirección x (ms-1)

ISSN 2448-5551 TF 137 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

. .
1. Introducción
.
V2
.. .
.
V2

.
Los flujos pulsátiles son relativamente comunes y, para su
análisis o estudio, parte de su complejidad proviene de su ΔP3
naturaleza transitoria. De manera general, se pueden
encontrar en la dinámica sanguínea y en la respiración.
ΔP4

.. .
ΔP2

Para resolver problemas que involucran flujos pulsátiles se

pueden emplear los métodos tradicionales en el campo de
.
.
.
V1
ΔP1
. .
V1
la mecánica de fluidos, estos son: la dinámica de fluidos
(a)
computacional (CFD), soluciones analíticas para casos
particulares de las ecuaciones que gobiernan el flujo de
fluidos, o bien, la experimentación.

..
Las simulaciones numéricas por medio de CFD de
cualquiera de los casos señalados involucran grandes .
tiempos de cómputo. Por otra parte, resolver las ecuaciones V2, ΔP3

..
de Navier-Stokes para geometrías complejas y en régimen .
. V2, ΔP4
transitorio es un verdadero reto. Sin embargo, los métodos V1, ΔP2
de parámetros concentrados son una alternativa para
.
rápidamente obtener estimaciones de las magnitudes más V1, ΔP1
importantes involucradas en el flujo de fluidos, es decir,
valores promedio la presión y la velocidad.
(b)
De manera específica, en el campo de otorrinolaringología . .
existe la rinomanometría, una prueba que sirve para
diagnosticar la obstrucción nasal. Dicha prueba permite . .
medir la diferencia de presión y la velocidad del flujo de .
V2.. .
.

.
aire en la vía nasal durante la respiración, Dibildox [1]. V2
Dado que ambas mediciones, la diferencia de presión y la

..
razón de flujo, se toman al mismo tiempo existe una ΔP3 ΔP4

..
correspondencia biunívoca entre ellas. Para ciertos casos,
cuando se grafica la diferencia de presión vs. flujo ΔP1 ΔP2
volumétrico, se observa que para un valor del flujo . . .
volumétrico existen dos valores distintos de diferencia de V1 V1
presión. A este fenómeno se le ha llamado histéresis y en .
este trabajo se le denominará histéresis hidrodinámica. (c)

Para ilustrar dicho fenómeno, se puede dividir el ciclo

respiratorio en 4 lapsos, los dos primeros durante la fase de
inhalación (un lapso de aceleración y otro de
desaceleración) y los otros dos lapsos en la fase de
exhalación (también con un lapso de aceleración y otro de
desaceleración);
. entonces, para cierto valor del flujo
.
V2, ΔP3 ..
volumétrico, V1, se presenta durante el lapso de inhalación-
aceleración un valor de diferencia de presión . ΔP1, luego
para el mismo valor del flujo volumétrico, V1, se presenta
.
V1, ΔP2 ..
V1, ΔP1
V2, ΔP4

durante el lapso de inhalación-desaceleración un valor de

diferencia de presión ΔP2; durante la fase de exhalación
ocurre algo similar,
. para un mismo valor del flujo (d)
volumétrico, V2, se presentarían las diferencias de presión
ΔP3 y ΔP4, ver Fig. 1. De acuerdo con lo descrito, la . .
diferencia de presión, ΔP, no sigue la misma trayectoria
durante el lapso inhalación-aceleración que durante el Figura 1 - (a) y (b) Histéresis debida a efectos inerciales [3], las
lapso inhalación-desaceleración; ΔP, tampoco sigue la gráficas (a) y (b) representan la misma información. (c) y (d)
misma trayectoria durante el lapso exhalación-aceleración Histéresis debida a efectos de acumulación de masa [7], las gráficas
que durante el lapso exhalación-desaceleración. (c) y (d) contienen los mismos datos.

ISSN 2448-5551 TF 138 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

En la literatura se encontró que la histéresis hidrodinámica, conservación de masa a un flujómetro de orificio y

se atribuye principalmente a tres motivos: el primero radica obtuvieron la Ec. 2, graficada en las Figs. 1 (c) y (d) junto
. .
en los efectos inerciales (aceleración), el segundo se refiere con VB = VBM sin(2πτ). Groß y Peters [7] hacen notar que
a los efectos de acumulación de masa (elasticidad de los incluso una variación pequeña de densidad puede producir
tejidos vivos, sellos inadecuados y compresibilidad del una derivada temporal que, multiplicada por un volumen,
aire) y el tercero reside en los errores de medición tiene un fuerte impacto en la diferencia de presión.
(suciedad en el equipo, ruido electrónico, desfasamiento de
la señal eléctrica en los sensores, entre otras). Actualmente,
VBM sin2  
2 2
p k  2 BVVBM 
existen sensores de presión y otros elementos usados para Δp = o A   1  cos2  (2)
la instrumentación que tienen una alta resolución temporal 2  ao AA   ao AB  Tr
2

pero los bucles en las curvas rinomanométricas aún se
manifiestan [2], de este modo los efectos de acumulación Debido a que existen dos puntos de vista sobre el origen de
de masa y los efectos inerciales pueden ser la causa de la la histéresis hidrodinámica, en el año 2005 Clement y
histéresis hidrodinámica, Fig. 1. Gordts [8] en una adenda titulada "Reporte de consenso
sobre rinometría acústica y rinomanometría (Consensus
Sosteniendo que la histéresis hidrodinámica se debe a los report on acoustic rhinometry and rhinomanometry)" dejan
efectos inerciales se encuentra Hoffrichter et al. [3] quienes claro que la histéresis se puede deber tanto a efectos de
haciendo uso de la ecuación de balance de cantidad de acumulación de masa como a efectos inerciales.
movimiento llegaron a la Ec. 1a; si para esta ecuación se
considera un . flujo
. volumétrico dado por un patrón A continuación se presenta un estudio de parámetros
sinusoidal, VB = VBM sin(2πτ), se llega a la Ec. 1b, misma concentrados del flujo compresible y pulsátil para un
que se grafica en las Figs. 1 (a) y (b). En la Ec. 1a se puede volumen de control con dos fronteras. De esto se deriva
ver que el segundo término del lado derecho está una expresión analítica con la que se muestra que las
multiplicado por (V/A2), este factor ronda un valor de histéresis hidrodinámica se puede manifestar por una
0.67 x 106 m-1 y provoca que la derivada del flujo combinación de efectos inerciales y de acumulación de
volumétrico con respecto al tiempo adquiera relevancia. De masa, tal y como señalaron Clement y Gordts [8].
acuerdo con Hoffrichter et al. [3], se puede mostrar que
algunos rinomanómetros que, en lugar de un bucle, siempre
2. Método
muestran la curva presión-flujo pasando por el origen, han
hecho algún manejo de los verdaderos datos medidos de
Para la siguiente deducción se emplea un sistema abierto
presión y flujo. En su estudio numérico Hörschler et al. [4]
similar al que plantean Groß y Peters [7], como se muestra
señalan que la histéresis (debida a efectos inerciales) se
en la Fig. 2. El sistema abierto está conformado por un
atenúa a bajas frecuencias respiratorias.
volumen V que representa el volumen de aire de la
máscara, adaptadores y mangueras usadas en los
o  2 o V dVB
Δp = VB  (1a) rinomanómetros, cabe señalar que este volumen también
2 A2 Tr A2 d podría incluir el volumen de la cavidad nasal. El orificio en
A produce una resistencia hidrodinámica similar al de una
2
o  VBM  nariz real. El orificio en B sirve para medir el flujo
Δp =   sin2  sin2   volumétrico. Además, la. bomba
2  A  . puede producir un flujo
(1b) pulsátil en B dado por: VB = VBM sin(2πτ). Las propiedades
 V
2 o BM 2 cos2 
V ρo y po son la densidad y la presión de de referencia. A
Tr A través de la ecuación de balance de cantidad de
movimiento en su forma integral, Ec. 3, se mostrará que la
Defendiendo la postura de que la histéresis hidrodinámica histéresis hidrodinámica puede darse por una combinación
se debe a los efectos de acumulación de masa se de efectos inerciales y efectos de acumulación de masa.
encuentran Fodil et al. [5] quienes reportaron la variación
del área transversal de tres regiones fisiológicas (válvula
.
nasal, cornete inferior y meato medio). Sus mediciones Δp(t) VB(t)
indican que el área de la sección transversal de las regiones p(t), ρ(t), V, A, U(t) po, ρo
señaladas decrece extremadamente rápido cuando se aplica pA pB
una presión negativa. Este comportamiento introduce un
fenómeno de histéresis hidrodinámica dentro de la curva A B
Figura 2 - Volumen de control con dos fronteras.
presión-flujo. Schumacher et al. [6] señalan que la
histéresis hidrodinámica es consecuencia de la elasticidad y

 F   u
que se acentúa con el incremento de la obstrucción nasal.
Por su parte, Groß y Peters [7] atribuyen la histéresis x x u  ndA   u dV x (3)
t
hidrodinámica a variaciones de densidad, ellos aplicaron la sc vc

ISSN 2448-5551 TF 139 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

Se llega a la Ec. 4 aplicando el método de parámetros Los cambios de la presión, p, son pequeños al compararlos
concentrados a la Ec. 3 en donde se tomaron en cuenta las con la presión de referencia, po, por lo que la relación de
variables de la Fig. 2 y se consideró que: el flujo es presiones, (p/po), es cercana a la unidad, entonces, los
unidireccional, las velocidades en A y en B son uniformes factores (p/po)1/k y (p/po)(1-k)/k también lo son. Sin embargo,
y para el volumen de control se tomó una velocidad media, no es el cambio de presión lo que cuenta, es más bien la
U, y una densidad media, ρ. derivada de la presión respecto al tiempo. Considerando lo
anterior, la Ec. (8) se simplifica a la Ec. (9)
p A AA  pB AB 
 VA VA VB VB 
  A AAU A U A   B ABU B U B  V dtd U 
(4)
1 o  p A AA  pB AB     
 AA AB 
(9)
Al  dU 1 d  p 
. escribir
. la Ec. 4 en términos de los flujos volumétricos, V U  
VA y VB, y al realizar la derivada del producto ρU se obtiene  dt k dt  po 
la Ec. 5.
Si se define el tiempo adimensional, τ, considerando Tr
p A AA  pB AB  como el tiempo de duración de un ciclo, se logra la Ec.
 VA VA V V   dU d  (5) (11).
  A  B B B   V  U
 A A AB   dt dt 
  t Tr (10)

Dado que se presentan pequeñas variaciones de densidad,

es posible realizar las aproximaciones de ρA ≈ ρo y ρB ≈ ρo.  VA VA VB VB 
La Ec. (5) se convierte en la (6).
1 o  p A AA  pB AB     
 AA AB 
(11)
 dU 1 d  p 
1 o  p A AA  pB AB   V
 U  
Tr  d k d  po 
 VA VA VB VB    dU d    (6)
   V   U  
 A A A B  
 o dt dt   o  Por otra parte, de la de conservación de energía, se tiene la
ecuación para calcular el flujo volumétrico en un
Además, considerando que los cambios de densidad son flujómetro de orificio. Los caudales en A y B están dados
pequeños, se puede hacer la consideración de flujo por las Ecs. (12) y (13), donde ζA = γA-2 y ζB = γB-2, siendo
isentropico, (ρ/ρo) = (p/po)1/k. Con esto en cuenta, se tiene la γA y γB, los coeficientes de flujo.
Ec. (7)
2p
VA  AA (12)
1 o  p A AA  pB AB    A o
 1 k 1k 
 VA VA VB VB  d   p   (7)
2 p  po 
 p dU
    V    U  
 AA AB   po  dt dt   po   VB  AB (13)
    B o
Con la derivada
Conocido el flujo volumétrico en B se sigue que
 1k 
d   p   1  p 
1k  
d p 2 p  po 
 
k 
VBM sin2   AB
      B o
dt   po   k  po  dt  po 
 
Y resolviendo para p/po se consigue la Ec. (14)
La Ec. (7) puede reescribirse como sigue
 V 2 
 o BM 2 B sin 2 2   1
p
(14)
 V V V V  po po 2 AB
1 o  p A AA  pB AB    A A  B B  
 AA AB 
En términos de las variables de referencia, ρo y po, para el
 1k
 1 k 
   (8)
flujo isentropico se tiene de White [10]
 p  dU 1  p  k  d  p  
V    U    
dt  po 
 o 
p dt k  po 
 o po  k ao2 (15)
 

ISSN 2448-5551 TF 140 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

Y la Ec. (14) se convierte en la Ec. (16), cuya derivada con 3. Resultados

respecto al tiempo adimensional es la Ec. (17).
Se indicó en la introducción que el fenómeno de histéresis
kV 2  B
sin 2 2   1
p hidrodinámica puede manifestarse por efectos inerciales o
 BM (16)
po 2 2
2ao AB por efectos de acumulación de masa. Con la Ec. (23) queda
claro que la histéresis hidrodinámica puede manifestarse
por una combinación de ambos efectos. Esto es más
 p  2kVBM 2
B
sin2 cos2 
d evidentes desde la Ec. (11); en dicha ecuación hay dos
   (17)
d  po 
2
ao AB2 derivadas. La primera de ellas, dU/dτ, es una aceleración y
representa la contribución de los efectos inerciales. La
segunda, d(p/po)/dτ, proviene de d(ρ/ρo)/dt, por lo que
Con las Ecs. (16) y (17), la Ec. (11) se transforma en la Ec.
representa la variación de la compresibilidad del fluido
(18).
respecto al tiempo, y es la contribución de los efectos de
acumulación de masa.
 VA VA VB VB 
1 o  p A AA  pB AB     
 AA AB  Una vez ubicados el término inercial y el termino de
(18) acumulación de masa, se hace uso del tiempo normalizado,
 dU 2VBM
2
B 
sin2 cos2 
V τ, para ubicar
 U . algunos puntos relevantes en las gráficas de
 d
2 2
Tr ao AB  Δp vs. VB, normalizadas. estás variables empleando sus
valores máximos, ΔpM y VBM. En la Tabla 1 se presentan
. los
Sustituyendo valores de τ para los que las Ecs. (22) y (17) y VB alcanzan
. . en esta ultima ecuación los caudales dados un máximo, un mínimo o adquieren un valor de cero.
por VB = VBM sin(2πτ) y la Ec. (12) se consigue la Ec. (19).
Tabla 1 – Para n par: valores de τ donde
. se alcanzan los valores
 p A AA  pB AB   2 AAp VBM
2
sin2  sin2   máximo, mínimo y cero de las funciones VB, dU/dτ, y d(p/po)/dτ.
      Máximo Mínimo Cero
 o    A  o AB  .
(19) VB (2n+1)/4 (2n-1)/4 n/4
V  dU 2VBM
2
B  dU/dτ n/2 (n+1)/2 (n+1)/4
  U sin2 cos2  d(p/po)/dτ (2n+1)/8 (2n-1)/8 n/8
Tr  d 2 2
ao AB  . .
Tomando en cuenta VB y dU/dτ, el flujo volumétrico, VB,
El término (pAAA - pBAB) se puede estimar como la adquiere el valor de cero, por ejemplo, en τ = 1; para ese
diferencia de presión en el flujómetro en A multiplicada mismo tiempo normalizado, dU/dτ alcanza el valor
por un factor α y el área promedio del volumen de control. máximo. En otras palabras, la aceleración es máxima
cuando se cambia de la fase de inhalación a la fase. de
 p A AA  pB AB   Ap exhalación y por este motivo en la gráfica de Δp vs. VB de
(20)
la Fig. 1(d) el bucle no pasa por el origen (0,0).
.
Se asume que la velocidad media en el volumen de control Considerando ahora VB y d(p/po)/dτ, se nota que el término
sigue el patrón sinusoidal dado por la Ec. (21) , en la que β de acumulación de masa alcanza su máximo por ejemplo
es una constante. La derivada con respecto al tiempo en τ = 1/8 y .un mínimo .en τ = 3/8, estos tiempos
adimensional de la Ec. (21) es la Ec. (22). corresponden a VB/21/2 ≈ 0.71VB en la fase de exhalación en
los lapsos de aceleración y desaceleración respectivamente.
.
VB VBM sin2  Con esto, el bucle de la gráfica de Δp vs. VB de de la Fig. 1
U   (21)
A A (b) se .ensancha al máximo para dichos tiempos. Además,
tanto VB como d(p/po)/dτ se hacen cero para τ = 0, 1/8, y 1,
dU 2 VBM cos2  y por esto el bucle sí pasa por el origen (0,0).
 (22)
d A .
En las gráficas de Δp vs. VB se habla de histéresis cundo
Con las Ecs. (21) y (22), el balance de momentum, Ec. hay un ensanchamiento en el bucle, por ello, se define la
(19), finalmente se convierte en la Ec. (23). histéresis hidrodinámica debida a efectos inerciales con la
Ec. (24). Esta ecuación representa la diferencia de presión
 A 2 AA  VBM 2 en el punto en el que ocurre el cambio de la fase de
p    sin2 sin2  exhalación a la fase de inhalación menos la diferencia de
  o  A  o  AB presión en el punto en el que ocurre el cambio de la fase de
(23)
2 VVBM   V  
2 inhalación a la fase de exhalación.
 cos2 1   B  BM sin2  
Tr A 
  ao AB   H I  p 2n1 4  p 2n2 4 (24)

ISSN 2448-5551 TF 141 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

La histéresis hidrodinámica debida a efectos de 4. Discusión

acumulación de masa queda definida con la Ec. (25)

H A  p 2n1 8  p 2n3 8

En la Sección 2 se dedujo la Ec. (23), en esta queda claro
(25)
que la histéresis hidrodinámica puede ocurrir por la
. aceleración del flujo dentro del volumen de control y
A continuación, en la Fig. (3), se grafican VB y Δp dada por también por la compresibilidad del fluido. Para encontrar
la Ec. (23). Para generar la gráfica se usaron los parámetros dicha ecuación se partió del método de parámetros
de la Tabla 2, es importante decir que estos datos también concentrados que no considera lo siguiente: detalles de los
sirvieron para realizar las gráficas de la Fig. 1. vectores de velocidad, características del flujo, existencia
de remolinos, cambios de dirección del flujo y el régimen
En la Fig. 3(b) se puede ver que el bucle no atraviesa el del flujo (laminar, transición o turbulento). No obstante, el
origen, este hecho denota la presencia de efectos inerciales, método de parámetros concentrados permite estimar las
además,
. el bucle
. presenta un marcado ensanchamiento en magnitudes más importantes del flujo en un volumen de
VB/21/2 ≈ 0.71VB que indica la existencia de efectos de control. Para la presente investigación estos datos fueron el
acumulación de masa. También, en la Fig. 3(a), se nota. un flujo volumétrico y la diferencia de presión.
desfasamiento entre las gráficas de Δp vs. VB, un
comportamiento característico de los flujos pulsátiles. Es importante señalar que la contribución de la presente
investigación radica en que de manera conjunta sí se
Tabla 2 – Valores de distintos parámetros empleados en las gráficas tomaron en cuenta: la compresibilidad del aire y la inercia.
de las Figs. (1), (3) y (4). A diferencia de Hoffrichter et al. [3] que únicamente
Variable Valor Variable Valor consideraron la inercia, y Groß y Peters [7] que solo
ao α
consideraron la compresibilidad del aire.
343 m/s 1
ρo 0.92 kg/m 3 β 0.04 Como se observa en la Fig. 4 (a), para un valor de
ζA 2.4 A A=( AA + AB)/2 V = 0.002 m3 se producen una bucle que es más ancho
ζB
cerca del origen. De acuerdo con la Ec. (24) se tendría una
2.4 AA AA=π DA2/4
histéresis hidrodinámica por efectos inerciales mayor a la
V 0.001 m 3 AB AB=π DB2/4 reportada en la Fig. 3 (b). En la Fig. 4 (b), para valores de
.
VBM 0.001 m3/s DA y DB 0.007 m ζA = 0.75 y ζB = 0.45 x 103 se. obtiene una curva en la que
Tr 5s po 101325 Pa el bucle es más estrecho en VB/21/2. De acuerdo con la Ec.
(25) se tendría una histéresis hidrodinámica por
. . acumulación de masa menor a la reportada en la Fig. 3 (b).

.. .
De este resultado se deriva que es posible atenuar la
. . histéresis hidrodinámica por efectos de acumulación de
V2 V2

.
masa al cambiar las propiedades geométricas del
ΔP3

.
flujómetro señalado en la Fig. 2.
ΔP4

. .
..
V1
ΔP1
ΔP2

..
V1
En la mayoría de las investigaciones, en las que se utiliza el
CFD como instrumento de análisis, se asume que la
naturaleza periódica del flujo de aire en el sistema
. respiratorio se puede despreciar. La suposición de que el
(a) régimen de flujo es permanente, se basa en que el número
de Womersley ronda el valor de 1, esto significa que la
frecuencia del ciclo respiratorio no es muy alta y permite
que el flujo volumétrico se encuentre en fase con la

..
diferencia de presión. Como se muestra en la Fig. 3 sí
. puede existir desfasamiento aún a bajas frecuencias
V2, ΔP3 respiratorias. Por lo tanto, para el futuro estudio de la
.
V1, ΔP2 .. .
.
V2, ΔP4
mecánica respiratoria es importante considerar flujo
transitorio y compresible. En este sentido, autores como
Doorly et al. [9] señalan que es importante tener en cuenta
V1, ΔP1 estado no permanente en el ciclo respiratorio, cuando se
presenta el fenómeno de histéresis, y en especial cuando se
está modelando algún proceso acumulativo (que guarda
(b) memoria del proceso), por indicar algunos: la deposición
. . de partículas y la transferencia de calor.

Figura 3 - (a) y (b) Histéresis debida a una combinación de efectos de

acumulación de masa y efectos inerciales. ζB = 1.5 x 103.

ISSN 2448-5551 TF 142 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

de 0.67 x 106 de acuerdo con los datos de la Tabla 2,

provoca que la aceleración del fluido en el volumen de
control adquiera relevancia

Los efectos de acumulación de masa, se originan por la

compresibilidad del aire, la elasticidad de los tejidos vivos,
la deformación de diversos componentes flexibles del
rinomanómetro, entre otras. Por otra parte, para tener
diagnósticos médicos veraces, es deseable contar con
rinomanómetros que minimicen los efectos de la histéresis
(a) hidrodinámica debida a efectos de acumulación de masa.
. . Se pueden atenuar dichos efectos al realizar cambios en las
propiedades geométricas de los flujómetros, como se vio
en la Sección 4. Por lo tanto, las ecuaciones expuestas en
este trabajo podrían ser útiles para diseñar rinomanómetros
en los que los efectos de la histéresis hidrodinámica por
acumulación de masa sean mínimos.

Por último, debido a la influencia que tiene la velocidad en

los coeficientes de: transferencia de calor por convección y
transporte de especies químicas; se tiene la sospecha de
que también existen fenómenos de histéresis para la
temperatura e histéresis para la concentración molar en el
(b) sistema respiratorio.
. .
REFERENCIAS

Figura 4 - Histéresis por combinación de efectos de acumulación de

masa y efectos inerciales. Se emplearon los valores de la Tabla 2 con [1] Dibildox Martínez Javier, Temas Selectos en
ligeros cambios (a) V = 0.002 m3 y (b) ζA = 0.75 y ζB = 0.45 x 103. Otorrinolaringología y Cirugía de Cabeza y Cuello,
(Cap. 20). Bloomington, IN: Palibrio LLC (2017). ISBN:
5. Conclusiones 978-1-5065-2167-1.
[2] Klaus Vogt, Klaus-Dieter Wernecke, Hans Behrbohm,
En el presente trabajo se realizó un estudio de parámetros Wolfgang Gubisch, Mara Argale, Four-phase
concentrados del flujo compresible y pulsátil para un rhinomanometry: a multicentric retrospective analysis
volumen de control con dos fronteras. Dicho volumen de of 36,563 clinical measurements. European Archives of
control representa el volumen de aire contenido en los Oto-Rhino-Laryngology 5 (2015).
diversos componentes de un rinomanómetro y en la [3] Hellmut Hoffrichter, Klaus Vogt, Wolfgang Althaus,
cavidad nasal. Se derivó un modelo matemático que mostró Wolfgang Hasse, 4-Phase-Rhinomanometry Basics and
claramente que la histéresis hidrodinámica que causa los Practice 2010 (The mathematical-physical concept of 4-
bucles en las curvas rinomanométricas es un fenómeno que phase-rhinomanometry). Rhinology 21 (2010).
puede aparecer como resultado de la combinación de [4] I. Hörschler, W. Schröder, M.Meinke, On the
efectos de acumulación de masa y de efectos inerciales. assumption of steadiness of nasal cavity flow. Journal of
Biomechanics 43 (2010).
Durante la respiración, es un hecho que el aire manifieste [5] Redouane Fodil, Lydia Brugel-Ribere, Céline Croce,
velocidades subsónicas. Además, los efectos inerciales Gabriela Sbirlea-Apiou, Christian Larger, Jean-François
están asociados a la periodicidad del ciclo respiratorio Papon, Christophe Delclaux, André Coste, Daniel
(aceleración local) y también a los cambios de dirección Isabey, and Bruno Louis, Inspiratory flow in the nose: a
del flujo inherentes al diseño del rinomanómetro y a la model coupling flow and vasoerectile tissue
morfología del sistema respiratorio (aceleración distensibility. Journal of Applied Physiology 98 (2005).
convectiva). Por los motivos mencionados, hay cierto [6] Michael J. Schumacher, Nasal Congestion and Airway
debate sobre la presencia de histéresis hidrodinámica Obstruction: the Validity of Available Objective and
debida a efectos inerciales. Sin embargo, en esta Subjective Measures. Current Allergy and Asthma
investigación se observó que la histéresis hidrodinámica Reports 2 (2002).
debida a efectos inerciales, HI, presenta una fuerte [7] T. F. Groß and F. Peters, A FluidMechanical
dependencia con la geometría, es decir, una dependencia Interpretation of Hysteresis in Rhinomanometry.
directa con la relación V/A2, misma que, al tener un valor International Scholarly Research Network 2011 (2011).

ISSN 2448-5551 TF 143 Derechos Reservados © 2018, SOMIM

 MEMORIAS DEL XXIV CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
19 al 21 DE SEPTIEMBRE DE 2018 CAMPECHE, CAMPECHE, MÉXICO

[8] P.A.R. Clement and F. Gordts, Consensus report on

acoustic rhinometry and rhinomanometry. Rhinology,
43 (2005).
[9] D.J. Doorly, D.J. Taylor, R.C. Schroter, Mechanics of
airflow in the human nasal airways. Respiratory
Physiology & Neurobiology 163 (2008).
[10] Frank M. White Fluid Mechanics, 7th ed. McGraw-
Hill (2009).

ISSN 2448-5551 TF 144 Derechos Reservados © 2018, SOMIM