Latest revision as of 14:53, 18 March 2020

Resumen

Mediante la interconexión de elementos eléctricos tales como una resistencia, un inductor y un condensador se ha desarrollado un modelo de orden reducido del sistema arterial capaz de describir el comportamiento del flujo sanguíneo. En dicho modelado están representados implícitamente las diferentes ecuaciones diferenciales que describen el comportamiento de este sistema. Esto permite obtener y regular rápidamente, la actuación del modelo conjunto, tanto en condiciones fisiológicas normales como ante situaciones patológicas, lo que lo convierte en una herramienta útil para la docencia o investigación en fisiología cardiovascular.

Se implementaron dos modelos arteriales diferentes, un modelo elástico y un modelo visco elástico, con el propósito de analizar el flujo y la presión sanguínea que se podrían producir ante diferentes condiciones fisiológicas. Con el objetivo de analizar ambos modelos, se contrastó usando diferentes órdenes de elasticidad, se estudió la patología de disección aórtica y se comparó los resultados de manera cualitativa con soluciones similares observadas en la literatura. Para poder realizar dichas simulaciones se implementó un código de resolución de ecuaciones diferenciales ordinarias (DAE) en Matlab.

Abstract

Through the interconnection of electrical elements such as a resistor, inductor and capacitor, has been developed a reduced-order model of the arterial system capable of describing the behaviour of blood flow. In this modelling, the different differential equations that describe the behaviour of this system are represented implicitly. This allows us to obtain quickly and regulate the performance of the joint model, both in normal physiological conditions and in pathological situations, which makes it a useful tool for teaching or research in cardiovascular physiology.

It was implemented two different arterial models, an elastic model and a viscoelastic model, with the aim of analysing the flow and blood pressure that could occur under different physiological conditions. In order to analyse both models, it was contrasted using different orders of elasticity, the pathology of aortic dissection was studied and the results were compared qualitatively with similar solutions observed in the literature. In order to perform these simulations, it was implemented in Matlab a code for solving ordinary differential equations (DAE).

Agradecimientos

Quería agradecer a mi tutor Eduardo Soudah, por darme toda la ayuda posible a su abasto cuando lo he precisado y por guiarme para coger las riendas del trabajo.

En especial a mis padres, por su infinita paciencia y su apoyo incondicional en cualquier situación imaginable.

Glosario

1. Prefacio

1.1. Iniciativa

El gran interés mostrado por el funcionamiento y el comportamiento del sistema cardiovascular y, además, el poder recrear el comportamiento fisiológico de éste mediante la simulación, fueron los principales motivos para el desarrollo de este trabajo. También, sumado al hecho de poder poner en práctica los conocimientos adquiridos en Matlab durante el grado de ingeniería biomédica.

1.2. Requisitos previos

El conocimiento y uso del programa Matlab era de vital importancia, puesto que el trabajo se realiza con él en su totalidad. Por ello, era necesario contar con un sólido conocimiento del lenguaje de programación de éste para poder hacer un buen uso de esta herramienta. A parte, también era importante tener un conocimiento previo de conceptos eléctricos y de hemodinámica sanguínea para recrear los modelos en RLC y sus respectivas ecuaciones de presión y flujo. Por último, y no menos importante, tener una base matemática sobre ecuaciones diferenciales ordinarias y de cómo resolverlas, ayuda en este caso, a la determinación y resolución de éstas mediante métodos numéricos.

2. Introducción

El presente trabajo está dentro del campo de la biomecánica¹. El objetivo de esta rama de la ciencia es utilizar el conocimiento de las leyes mecánicas, herramientas y formulación matemática para elaborar modelos lo más realistas posibles del comportamiento de estructuras biológicas. Por ejemplo, el estudio del flujo sanguíneo, análisis de órganos, etc. La primera característica de la biomecánica es que es una ciencia multidisciplinar que involucran a la biología y mecánica. En este caso, la colaboración entre ambas disciplinas consiste en que la biología aporte el estudio de los fenómenos vitales, y la mecánica proporcione la teoría, los modelos matemáticos y los métodos de solución. Estos modelos matemáticos permiten hacer la simulación y a partir de éstas, se pueden hacer predicciones sobre el estado futuro de los elementos de estudio.

A su vez, dentro de la biomecánica, se encuentra la hemodinámica, la rama que tiene por objetivo estudiar el flujo sanguíneo. Este estudio no se focaliza sólo en el movimiento de la sangre, como obtener el caudal, flujo y presión, sino que también en las fuerzas que genera dicho flujo sobre los elementos sobre los que interactúa, como los vasos sanguíneos o el corazón. Sin embargo, este trabajo se centrará en el análisis y evaluación de los primeros.

(¹) Según la Real Academia Española (RAE), es la ciencia que estudia la aplicación de la mecánica en las estructuras y órganos de los seres vivos.

2.1. Motivación

Hoy en día las enfermedades cardiovasculares representan la principal causa de muerte entre los países desarrollados y contribuye a casi el 30% de todas las enfermedades mundiales mortales según la Organización Mundial de la Salud (OMS). Esas patologías están estrictamente relacionadas con las condiciones de la circulación sanguínea, por lo tanto, una compresión de la hemodinámica local puede tener aplicaciones útiles para la investigación médica.

No obstante, la extrapolación directa de estos parámetros hemodinámicos es complicada, pues la circulación de la sangre está dentro del cuerpo vivo, por lo que debe haber un compromiso entre la precisión de las mediciones y su nivel de invasividad al determinarlas [30].

La necesidad de seguir investigando acerca de estos sistemas debido a la trascendencia e impacto en la salud mundial que tienen, hace que se realicen estudios y descubrimientos de manera continuada con el fin de mejorar los tratamientos ya existentes, descubrir factores de relevancia y prevenir el efecto de diferentes hábitos o costumbres en la fisiología de este sistema tan importante. Es en este marco donde se sitúa este trabajo, en servir de utilidad recreando mediante una simulación numérica el flujo sanguíneo en una situación característica, pero aplicable a otras circunstancias [16,22,23].

Este hecho permite una interpretación sencilla y rápida, mediante la analogía a un sistema eléctrico menos complejo –al tratarse de un sistema de orden reducido-, de lo que está sucediendo y anticipándose al comportamiento de una red arterial en concreto en una situación específica.

2.2. Objetivo del trabajo

El propósito de este trabajo es desarrollar una herramienta numérica basada en modelos de orden reducido (0-D) para resolver modelos matemáticos que describan la presión y el flujo sanguíneos. Estos modelos deben son capaces de predecir las características hemodinámicas dentro del sistema arterial ante diferentes situaciones fisiopatológicas. La sangre se simula como un fluido homogéneo, incompresible y newtoniano², mientras que las arterias serán como tubos cilíndricos, rectos y deformables. El comportamiento arterial se representará mediante pequeños circuitos eléctricos compuestos por una resistencia, un inductor y un condensador (RLC).

Estos modelos se aplicarán, primeramente, sobre una sección de la aorta, simulando el primer flujo de entrada, y después se ampliará hacia el tronco aórtico, modelando una patología denominada disección aórtica. Se representará esta patología sobre un modelo elástico y visco elástico, donde se colocará una resistencia extra en serie con el condensador para simular la propiedad de viscoelasticidad del vaso.

Para poder realizar la simulación, se utilizará el programa Matlab, con el que se creará un código computacional en el que se introducirán los datos fisiológicos de la sección arterial con la que se va a trabajar, y mediante un sistema matricial de ecuaciones diferenciales, que se resolverán, se obtendrán los parámetros requeridos que son las ondas de flujo y presión sanguíneas de cada uno de estos modelos y se compararan. También se someterá al modelo visco elástico a una evaluación cuantitativa, aumentando y disminuyendo los valores de la elasticidad para ver el comportamiento de la presión sanguínea. Finalmente se contrastarán cualitativamente ambos modelos y con soluciones similares observadas en la literatura.

(²) Un fluido newtoniano es aquel cuya viscosidad puede considerarse constante. Este hecho hace que sea más sencillo describir su comportamiento.

2.3. Estado del arte

La irrupción de la aplicación de las matemáticas para representar fenómenos físicos supuso un cambio de perspectiva, pasando de descripciones y explicaciones cualitativas del sistema cardiovascular a una cuantificación del comportamiento de éste mediante ecuaciones matemáticas. A continuación, se da una visión histórica breve de esta evolución.

Pinceladas históricas

Partiendo del siglo IV a.C. Aristóteles describió la comunicación del corazón con los vasos sanguíneos, aunque no identificó todavía la circulación la circulación sanguínea. Siguiendo al siglo III a.C. Praxágoras, distinguió las funciones de las arterias y las venas, observando que las arterias impulsaban sangre y las venas no [28,29]. A finales del siglo II, el griego Galeno propuso que se propagaban los pulsos de presión originados en el corazón a través de la red arterial [24].

Durante un período más avanzado, alrededor del siglo XIII, Ibn Al-Nafi fue el primer médico que describió correctamente la circulación pulmonar [24]. Antes de eso, la teoría común y popular del flujo sanguíneo era que la comida se convertía en sangre en el hígado y luego funcionaba como un combustible.

Más tarde, tocando el siglo XVII, el inglés William Harvey, mediante unos estudios donde midió el caudal de la sangre que pasa por las venas, concluyó que la sangre bombea desde el corazón y por el cuerpo en sentido unidireccional [27].

A mediados del siglo XVIII, Malpighi y Van Leeuwenhoek describieron la existencia de capilares que unían venas y arterias, demostrando definitivamente que el sistema circulatorio era cerrado y el flujo unidireccional.

Durante el mismo siglo, Daniel Bernoulli investigó las leyes que rigen la presión arterial y relacionó la presión sanguínea con la velocidad de la sangre en la siguiente ecuación:

donde V es la velocidad del fluido en la sección considerada, ${\textstyle \rho }$ es la densidad del fluido –tomando como un flujo incomprensible, por lo que este valor es constante-, P es presión del fluido, tomando z como la diferencia de altura desde una cota de referencia y g como la gravedad. Además, esta ecuación era aplicable en flujos laminares con un caudal constante y viscosidad nula.

En el siglo XIX, el francés Jean Lonard Marie Poiseuille estableció la ley que lleva su nombre, dónde relaciona la ley física que explicaba la relación de la caída de presión y el flujo sanguíneo en condiciones de flujo estacionario. La ecuación rige del siguiente modo:

donde ${\textstyle Q}$ es el caudal, ${\textstyle \Delta p}$ la diferencia de presiones, ${\textstyle r}$ el radio del tubo, ${\textstyle L}$ la longitud del tubo y ${\textstyle \mathrm {\mu } }$ la viscosidad del fluido.

A lo largo del siglo XX, Otto Frank estableció el modelo “Efecto Windkessel”, donde la aorta se concibe como un tubo elástico con capacidad de almacenar fluido [25]. Así, en el extremo del corazón el fluido es introducido de forma intermitente, mientras que en otro extremo el fluido sale de forma aproximadamente constante. Fue entonces cuando, el sistema circulatorio fue concebido como un embalse elástico dónde el corazón bombea sangre y desde el cual sale una red de conductos no elásticos para abastecer el cuerpo.

Se elaboraron modelos matemáticos que simularon esta teoría y más tarde, Noordagraaf y Westerhof [25], utilizaron el concepto de Windkessel y construyeron un modelo analógico eléctrico de las principales arterias de la circulación sistémica y lo aplicaron para paredes de tubos elásticos. En su estudio, las principales arterias de la circulación sistémica se dividieron en un numero finito de segmentos arteriales. Cada segmento se representa por un circuito eléctrico, compuesto por una resistencia, un inductor y un condensador.

Hoy en día, el flujo sanguíneo puede ser modelado con diferentes niveles de precisión, que incluyen: modelos de orden reducido (0-D), modelos unidimensionales (1-D) [22], y modelos tridimensionales (3-D) [23]. La principal diferencia entre los modelos de orden reducido y los unidimensionales y tridimensionales es que los primeros, no tienen dependencia espacial, y solo dependen del factor tiempo.

Los modelos 0-D vienen definidos por una ecuación diferencial ordinaria, ya que es una función de una sola variable –el tiempo-Por el contrario, los otros modelos incluyen dimensiones espaciales, en el caso del 1-D una dimensión espacial, y en el caso del 3-D, tres dimensiones espaciales, x, y y/o z, dentro de un sistema de coordenadas elegido.

No obstante, para el presente trabajo, se ha escogido desarrollar el modelo en un sistema de orden reducido, 0-D.

Con el fin de entender mejor el estudio de la mecánica vascular, se brinda a continuación, una breve descripción del sistema y del funcionamiento del mismo y algunos aspectos más generales de la mecánica del corazón, el sistema vascular e insistiendo en el sistema arterial, pues será el sistema en el que se trabaja.

2.4. El sistema circulatorio

El sistema circulatorio se trata fundamentalmente de un elemento de transporte que facilita el desplazamiento de diferentes substancias a través del flujo sanguíneo mediante una bomba muscular –el corazón- en un circuito cerrado de tubos elásticos – vasos sanguíneos-.

De esta manera, el flujo sanguíneo presenta un papel crucial dentro del sistema circulatorio (SC) tanto en el transporte y distribución de nutrientes, oxígeno y agua al organismo, ajuste de la temperatura corporal, como en la eliminación de subproductos metabólicos (CO₂ y agua) de los tejidos corporales.

Existen dos tipos de circulación que sigue el flujo, que se corresponden con dos trayectorias distintas a través de dos sistemas de vasos sanguíneos independientes:

a. La circulación sistémica, que es la que parte del corazón (concretamente del ventrículo izquierdo) con sangre oxigenada y limpia, y a través de las arterias, arteriolas y capilares, llega a todos los tejidos para realizar el suministro. Después regresa por las venas hasta el corazón, entrando en éste por la aurícula derecha.

b. La circulación pulmonar, que se encarga de eliminar los gases de desecho, como el dióxido de carbono (CO₂) y tomar nuevo oxígeno. Parte del ventrículo derecho del corazón y lleva la sangre empobrecida en oxigeno hacia los pulmones, donde se realiza el cambio de gases. Posteriormente la sangre ya enriquecida con dicho oxígeno, vuelve al corazón entrando por la aurícula izquierda.

La trayectoria de la circulación sistémica es más larga y presenta más ramificaciones respecto a la pulmonar. Por ello, las presiones necesarias para cada una de ella son diferentes, siendo mayores en la sistémica.

Figura 1. Circuito de circulación arterial (rojo) y circulación venosa (azul) en el cuerpo humano Fuente: https://accessmedicina.mhmedical.com/content.aspx?bookid=1858&sectionid=134366226

2.4.1. El corazón

El corazón es el órgano muscular del sistema circulatorio que actúa como doble bomba que impulsa la sangre a lo largo de las dos circulaciones pulmonar y sistémica. Para mantener estos dos circuitos separados, el corazón se divide en dos partes -derecha e izquierda-, cada una de ellas con dos cámaras conocidas como aurículas y ventrículos: ventrículo izquierdo (LV), ventrículo derecho (RV), aurícula izquierda (LA) y aurícula derecha (RA) (véase Figura 2). El bombeo de la sangre se proporciona a través de dos fases que se alternan en cada latido, una fase de contracción –dónde se produce la eyección de la sangre-, llamada sístole, con una fase de relajación, llamada diástole. La repetición continua de estas dos fases representa el ciclo cardíaco (véase Figura 3). Durante la sístole, la sangre se bombea hacia la aorta y durante la diástole, la válvula mitral se abre, y la sangre fluye hacia la cavidad ventricular desde la aurícula, para ser impulsada nuevamente.

Figura 2. Recorrido del flujo sanguíneo a través del corazón Fuente: http://miaulaenlared.blogspot.com/2018/01/6-primaria-el-corazon-y-el-recorrido-de.html

Figura 3. Fase sístole y diástole durante un ciclo cardíaco Fuente: D.Corridini. A numerical tool for the simulation of the blood flow, 2014

La duración media común del ciclo cardíaco es de aproximadamente 0.80 segundos, que equivale a unos 70 latidos por minuto. El caudal en condiciones normales puede ser de 4,9 litros/minuto, aunque en situaciones de máximo esfuerzo puede llegar a 35 litros/minuto. El valor del caudal en reposo de 4900 ml/minuto se ha calculado suponiendo 70 latidos/minuto y un volumen de bombeado de 70 ml. Entonces se obtiene 70·70=4900 ml/min.

2.4.2. El sistema vascular

Los vasos sanguíneos que componen el sistema vascular (SV) son las arterias, las venas y los capilares, los cuales se diferencian por su tamaño, el tejido que los compone y el diámetro del conducto. Las arterias-que forman el sistema arterial (SA)-, son las de mayor diámetro, son vasos elásticos fuertes, capaces de transportar la sangre a grandes presiones. Inician el circuito llevando la sangre oxigenada desde el corazón al resto de los órganos.

Las funciones principales del SA dentro del SV son, primeramente, preservar una presión sanguínea adecuada. Como el corazón es una bomba intermitente, el avance de la sangre que se produciría sería únicamente en la fase de eyección o fase sistólica. Sin embargo, la estructura elástica de la pared arterial, permite almacenar parte de la energía cinética proporcionada durante la contracción del ventrículo, como energía potencial que actúa expandiendo la pared. Así, posteriormente, durante la diástole, la pared vuelve a su estado de reposo, transmitiendo energía cinética de la sangre acumulada en su interior. Debido a estos cambios de presión, el SA actúa como un amortiguador pulsátil, pues las oscilaciones generadas por estas variaciones de presión y flujo son atenuadas según la sangre avanza por la red arterial hacia el corazón. Esta disminución de las fluctuaciones es realizada principalmente por arterias menores y arteriolas. Esto asegura continuidad tanto de la presión periférica arterial como del flujo sanguíneo a nivel capilar. Del mismo modo, es capaz de regular la distribución del flujo sanguíneo. Esta función recae básicamente en los las arteriolas o vasos de resistencia³. La pared arteriolar presenta un mayor contenido de fibras musculares y, por lo tanto, es el vaso principal que puede modificar su diámetro ampliamente y controlar así el flujo que llega a un determinado segmento vascular. Este cambio diametral de las arteriolas –vasoconstricción⁴ y vasodilatación⁵- está regulado a través de factores nerviosos y hormonales, atendiendo siempre, a las necesidades del organismo.

Figura 4. Diagrama simplificado del sistema vascular. Las secciones en rojo pertenecen al sistema arterial (sangre oxigenada) y las secciones en azul al sistema venoso (sangre desoxigenada) Fuente: https://es.wikibooks.org/wiki/Anatom%C3%ADa/Aparato_Circulatorio

No obstante, la pared arterial no es homogénea. Consta de tres capas concéntricas, que de dentro a fuera se conocen como: la túnica íntima o interna –capa interna-, constituida por endotelio y está presente en arterias y venas. La túnica media –capa media-, compuesta por fibras musculares lisas y fibras de colágeno, y la túnica externa o adventicia –capa externa-, formada por tejido conectivo laxo, compuesta por fibroblastos y colágeno.

Además, está el lumen, que es el espacio interior de la arteria, por donde fluye la sangre.

También se encuentran diferencias entre la composición de cada capa en los distintos vasos, encontrando arterias elásticas que conforman las grandes arterías como la aorta, la arteria pulmonar o la carótida y suelen presentar un gran diámetro que puede variar de 1 cm a 2.5 cm. Presentan una gran área del lumen, y son conocidos como vasos de baja resistencia.

Figura 5. Partes de la arteria en comparación con la vena Fuente: https://periodicosalud.com/diferencias-arterias-venas/

Las arterias elásticas tienen una gran cantidad de filamentos de colágeno y elastina en la túnica media, que ayudan a mover y almacenar el flujo de sangre durante la sístole y a mantener una presión relativamente constante a su paso.

Las arterias musculares cuya función básica es distribuir la sangre a los diferentes órganos del cuerpo. Su diámetro oscila entre 1 cm y 0.2 cm. Tienen una gran cantidad de fibras musculares lisas en la túnica media que ayudan a regular el flujo sanguíneo.

Y finalmente las arteriolas. Las arteriolas tienen un diámetro del lumen pequeño respecto a las mencionadas con anterioridad que varía entre 10 μm a 0.2 mm y juegan un papel fundamental en el mantenimiento y la regulación de la presión arterial mediante la contracción variable del músculo liso de sus paredes y la regulación del aporte sanguíneo a los capilares.

Después, el flujo fluye por las arterias y llega a las arteriolas y es dirigido a los capilares, una red de vasos de diámetro muy pequeño, donde se lleva a cabo un intercambio de nutrientes hacia los tejidos. Los capilares unen las arteriolas y las pequeñas ramificaciones procedentes de las venas, las vénulas.

Finalmente, la sangre se recoge en el sistema venoso, el que los vasos se fusionan rápidamente en vasos de mayor tamaño que recogen la sangre de los tejidos y órganos y la llevan de vuelta al corazón. El volumen total del sistema venoso es mucho mayor que el del sistema arterial, aunque las amplitudes de presión y velocidades de flujo son menores respecto a las arterias.

Cada segmento del sistema vascular tiene una combinación ideal de tamaño, composición de pared, espesor y área para trabajar adecuadamente en cada circunstancia. Por ejemplo, la aorta y las arterias de mayor tamaño tienen una pared con un mayor contenido de colágeno y elastina debido que su principal función es la de enviar sangre a alta presión a todas las regiones del cuerpo que abastece. Gracias a ello, pueden deformarse y recuperar su tamaño original después de la eliminación de la carga. Véase en la Tabla 1, como la elasticidad de las principales arterias y venas en función de su tamaño y sección.

Cabe añadir que, como se ha mencionado en párrafos anteriores, en términos fisiológicos, el comportamiento de las paredes de los vasos sanguíneos presenta una propiedad elástica llamada viscoelasticidad. Dicha propiedad, proporciona el denominado “mecanismo de almacenamiento de energía”, como un resorte comprimido o un condensador eléctrico para almacenar la energía de la sangre durante la sístole (Figura 5).

Figura 6. Comportamiento de presión (p) y flujo (q) en la pared arterial durante el latido Fuente: R.Gul. Mathematical Modeling and Sensitivity Analysis of Lumped-Parameter Model of the Human Cardiovascular System. PhD thesis , 2016.

La propiedad viscosa disipa la energía transmitida a las paredes del vaso durante el latido. En la diástole, la sangre almacenada se empuja hacia la parte inferior de la red arterial y como resultado, se produce un flujo sanguíneo suave y pulsátil a través el sistema circulatorio.

Estas características de comportamiento se tendrán en cuenta para la realización de este trabajo, pues se determinarán flujos y presiones, por una parte, con variables elásticas y por otra, visco elásticas.

(³) Se conoce como resistencia vascular a la dificultad que opone un vaso sanguíneo al paso de la sangre a su través. Si la resistencia vascular disminuye, el flujo aumenta, por el contrario, si la resistencia vascular aumenta, el flujo disminuye

(⁴) Disminución del diámetro de un vaso por contracción de las fibras musculares.

(⁵) Aumento del diámetro de un vaso por la relajación de las fibras musculares.

2.4.3. La sangre

La sangre es un tejido líquido encargado de suministrar sustancias necesarias como nutrientes y oxígeno a las células que forman los tejidos y de recoger los desechos metabólicos de cada uno. Está formada por partículas en suspensión conocidas los glóbulos rojos, los glóbulos blancos y las plaquetas (45% del volumen total, y una parte líquida llamada plasma que representa el otro 55% restante y está compuesto mayoritariamente de agua).

Los glóbulos rojos que componen la sangre, también conocidos como eritrocitos, son los responsables del intercambio de oxígeno y dióxido de carbono en las células. Por otro lado, los glóbulos o leucocitos, juegan un papel importante en el sistema inmunológico humano, pues activan la respuesta inmunitaria del organismo, interviniendo en su defensa contra sustancias extrañas o agentes infecciosos. Las plaquetas o trombocitos, principales responsables de la coagulación sanguínea.

El plasma a su vez, es una solución electrolítica diluida que contiene varias proteínas y diferentes sustancias disueltas. También transporta nutrientes, azucares, hormonas y substancias de desecho contenidos en la sangre. Al contener una gran cantidad de agua, y al ser ésta incomprensible, da las propiedades de incompresibilidad a la sangre.

Figura 7. Composición de la sangre Fuente:https://www.google.es/search?q=composicion+de+la+sangre&source=lnms&tbm=isch&sa=X&ved=0ahUKEwi8_7_ViePdAhWSK1AKHeXCC98Q_AUICigB&biw=1280&bih=617&dpr=1.5

Desde un punto de vista mecánico, la sangre está compuesta por glóbulos rojos flexibles suspendidos en el plasma, que hacen que se comporte como un fluido no newtoniano debido a la heterogeneidad de su composición. Por tanto, cualquier modelo constitutivo que se proponga ha de reflejar como mínimo estas dos características, viscosidad e incomprensibilidad. No obstante, para el desarrollo de este trabajo, se toma la sangre como fluido newtoniano con una viscosidad propuesta de valor constante.

2.4.4. Ondas de presión y flujo

La presión sanguínea es el elemento que se ejecutara durante este trabajo. El parámetro de presión es la fuerza por unidad de superficie que ejerce la sangre contra las paredes vasculares, que viene determinada por el volumen de sangre que contiene el SA y por las propiedades de las paredes, midiéndose normalmente en mmHg, dina/cm² o Pascales. Por otro lado, el flujo o caudal es el volumen de sangre por unidad de tiempo que atraviesa el vaso, que se puede expresar en ml/s o l/min.

El corazón envía sangre de un modo pulsátil a las arterias. En cada sístole entra un volumen de sangre en la aorta provocando un incremento de presión, dando lugar a su punto máximo, conocido como presión sistólica. Del mismo modo, durante cada diástole, la presión disminuye representando su punto mínimo, presión diastólica. Estos cambios cíclicos generan una onda de presión también llamada pulso arterial.

En esta onda, hay una sección ascendente que se desarrolla durante la sístole ventricular, y una fase descendente, que se inicia con una breve disminución de presión, volviendo a aumentar y generando una muesca que presenta el nombre de incisura dicrota, producto del cierre de la válvula aortica. –la válvula aórtica se cierra cuando la presión en la aorta supera a la ventricular-. La presión arterial media es un parámetro cardiovascular importante ya que proporciona el valor de presión con que la sangre llega a los tejidos.

La diferencia entre la presión arterial sistólica y la diastólica se denomina presión diferencial o presión de pulso.

Figura 8. Curva de presión durante la sístole y la diástole Fuente: https://ocw.unican.es/mod/page/view.php?id=538

Las formas de la onda de presión y su velocidad dependen tanto del rendimiento del corazón como de las características de la red vascular. La Figura 9 muestra cómo la presión disminuye a medida que se aumenta la distancia respecto al corazón.

Figura 9. Presión pulsátil en el árbol vascular Fuente: D.Corridini. A numerical tool for the simulation of the blood flow, 2014

Durante el período de diástole, la presión en el ventrículo izquierdo es cercana a cero y la válvula aórtica se mantiene cerrada a causa de la diferencia de presión de aproximadamente unos 80 mmHg. En la aorta, la presión arterial media es aproximadamente de unos 100 mmHg con una presión sistólica de 120 mmHg y una presión diastólica de 80 mmHg. En las arteriolas la presión disminuye considerablemente y tanto la presión diastólica como la sistólica convergen gradualmente. Cuando llega el flujo sanguíneo al sistema de capilares, el flujo es aproximadamente constante y más lento, lo que facilidad el intercambio mediante difusión, de nutrientes y productos de desecho entre la sangre y el tejido enervado.

Observación: Para un adulto joven y sano los valores aceptables de presión sistólica y diastólica oscilan entre 120 mmHg y 80 mmHg respectivamente. No obstante, factores constitucionales (sexo, raza, peso) y del estilo de vida (dieta, hábitos como el consumo de tabaco o alcohol, etc.) influyen de forma muy importante en la presión arterial haciendo variar los valores y la forma de la onda.

3. Modelo matemático

En el capítulo anterior se proporcionó un conocimiento básico preliminar sobre el sistema cardiovascular. A lo largo de este capítulo veremos los modelos y las ecuaciones matemáticas utilizadas que describen el flujo y la presión sanguínea mediante un modelo elástico lineal y otro visco elástico [12]. Los modelos matemáticos permiten, no solo estudiar la propagación de las ondas en un segmento arterial aislado, sino también el sistema de circulación global, o al menos, aproximar su comportamiento.

3.1. Modelos de orden reducido 0-D

El modelado en 3D del flujo sanguíneo en el sistema cardiovascular puede proporcionar simulaciones muy detalladas, pero las derivaciones e implementaciones y la complejidad de dicho sistema lo hacen computacionalmente muy costoso y complejo.

Sin embargo, es posible derivarlo a modelos más simples, sencillos y menos costosos utilizando modelos de orden reducido, 0-D.

Estos modelos describen la física del sistema muy simple y el cambio de parámetros da una idea del comportamiento del sistema implementado. Se pueden encontrar diferentes modelos según los parámetros utilizados y si trabajan en un solo bloque –modelo unicompartimentado- o que se trabajan en varios bloques – modelo policompartimentado-.

3.1.1. Modelo unicompartimentado

En un modelo unicompartimentado, la red arterial o venosa de todo el sistema se representa con un solo circuito, compuesto de una resistencia, una inductancia y/o un condensador. Este tipo de modelo se empezó a utilizar para modelar la circulación del cuerpo entero, –un bloque para la circulación arterial, otro bloque para la circulación pulmonar, etc.- pero estaba limitado puesto que no puede describir los cambios de presión y de caudal en ubicaciones específicas de la red arterial (Figura 10).

Figura 10. Modelo esquemático tipo unicompartimentado simulando la circulación en cada sistema. Cada mini-modelo equivale a un sistema específico

3.1.2. Modelo policompartimentado

En los modelos policompartimentados, por el contrario, la red arterial se divide en múltiples segmentos; cada arteria, como en el caso de este trabajo, está representada por un circuito eléctrico RLC (Figura 11).

Figura 11. Modelo esquemático tipo policompartimentado. Cada mini-modelo equivale a un vaso

Los valores de estos parámetros se calculan a partir de las propiedades de los vasos sanguíneos –módulo de Young, diámetro del vaso, longitud del vaso, grosor, viscosidad del fluido, entre otros.

Con el fin de construir una red de vasos llenos, los circuitos eléctricos de cada segmento se conectan entre sí utilizando las condiciones adecuadas de salida y de bifurcación. Las condiciones de bifurcación se derivan de leyes de conservación de masa y del momento. Por eso, antes de construir la red, es fundamental conocer y considerar las características de cada segmento, es decir, su comportamiento mecánico, condiciones de bifurcación de flujo, etc.

Por lo tanto, los modelos que contienen parámetros agrupados en bloques poli-compartimentados, se consideran más adecuados para simular, estudiar y analizar el comportamiento general del sistema cardiovascular. Este tipo de modelo será el que se empleará en este trabajo.

3.2. Modelo matemático en 0-D

El modelado matemático es una herramienta útil y eficiente para la predicción e investigación de enfermedades cardiovasculares. El objetivo básico de estos modelos es una mejor compresión del comportamiento sistema cardiovascular de una manera rápida y no invasiva.

Las variables de estado hemodinámico están representadas por un numero finito de variables, definidas por segmentos vasculares específicos. Aunque dicho modelo tiene baja resolución espacial, es útil para evaluar el rendimiento general del sistema o subsistemas cardiovasculares, como las arterias del brazo, la bifurcación de la carótida, entre otros.

Por lo tanto, un segmente del vaso sanguíneo se puede representar por un circuito eléctrico, combinando una resistencia, un inductor y un condensador. Existe una analogía entre el flujo en un tubo y el flujo de corriente en un circuito eléctrico. El potencial eléctrico o voltaje (V) corresponde a la diferencia de presión de conducción (Δp) en el segmento determinado del vaso, la corriente (A) es análoga al flujo sanguíneo (q), la resistencia (R) es la representación de la resistencia del flujo viscoso, el capacitor es análogo a la complianza del tubo (C), el inductor es la inercia de la sangre (L) y la impedancia corresponde a la resistencia terminal de los vasos terminales (R_b). (Véase tabla 2)

Dichos parámetros se pueden calcular usando ecuaciones matemáticas-sistema de ecuaciones diferenciales ordinarias acopladas que describen la variación en el tiempo- que dependen de las propiedades físicas y dimensiones de dichos vasos arteriales para la determinación de parámetros como presión y flujo.

En este apartado, se presenta la formulación matemática de los modelos con sus parámetros correspondientes, pues el trabajo se centrará en la representación de varios segmentos específicos del sistema arterial del cuerpo humano, diferenciándolos en su comportamiento visco elástico y elástico y comparándolos entre ellos (m³·Pa^-1).

3.2.1. Modelo equivalente 0-D al segmento vascular

3.2.1.1. Modelo elástico

Los modelos elásticos en estas mismas dimensiones se basan en la idea de que el flujo atraviesa un tubo cilíndrico rígido y recto de longitud establecida l, espesor h y diámetro Ø. Teniendo en cuenta que los parámetros que lo componen simulan el comportamiento elástico de éste, su representación se adhiere a la Figura 12.

Figura 12. Dimensiones y características de un vaso sanguíneo representado

El circuito que modula un solo vaso arterial se ve representado por una resistencia, un condensador y una bobina, en el que se podrá calcular la presión que equivale al voltaje que lo atraviesa y el flujo, que equivale al amperaje. El modelo elástico se puede representar en un circuito L-invertido tal que así,

Figura 13. Circuito equivalente (RLC) de un vaso. Qin y pin representa la corriente y el voltaje del fluido en la entrada del vaso respectivamente. Qout y pout es la corriente y el voltaje de salida del vaso

Los parámetros a determinar en el circuito se describirán a continuación.

3.2.1.1.1 Parámetros de flujo y presión

a) Ecuaciones de presión y flujo en segmentos no terminales

Como se ha precisado con anterioridad, mediante la analogía electrica se determinaran los parámetros fisiologicos de presión y flujo. Para ello, se definiran, dentro de cada segmento, dos únicas incognitas que, tras aplicar las leyes de corriente y voltage de Kirchhoff se obtendrá un sistema de dos ecuaciones diferenciales para presión ${\textstyle p}$, y flujo ${\textstyle q}$.

De esto modo, serán las incógnitas en cada vaso a determinar, ${\textstyle {\overset {\cdot }{{q}_{in}}}}$ y ${\textstyle {\overset {\cdot }{{p}_{out}}}}$, que se definirán como:

No obstante, en el presente trabajo, la presión de entrada se calcula a partir de un flujo de entrada teórico definido por una senoide, entonces se hace un pequeño cambio en las ecuaciones para que puedan aplicarse en todas las situaciones.

Por un lado, sabiendo que, la ecuación diferencial del cálculo del flujo en la arteria i, donde ${\textstyle {\overset {\cdot }{{q}_{i}}}=}$${\displaystyle \,{\frac {d{q}_{i}}{dt}}}$, se tiene que:

siendo ${\textstyle {p}_{(i-1)}}$ la presión de salida del segmento anterior, ${\textstyle {p}_{i}}$ la presión, ${\textstyle {R}_{i}}$ la resistencia, ${\textstyle {L}_{i}}$ la inductancia y ${\textstyle {q}_{i}}$ el flujo que se quiere determinar en un segmento indicado.

Figura 14. Analogía eléctrica entre dos vasos arteriales y el circuito eléctrico en segmentos no terminales

Por otro, la ecuación diferencial del cálculo de la presión en la arteria i, donde ${\textstyle {\overset {\cdot }{{p}_{i}}}=}$${\displaystyle \,{\frac {d{p}_{i}}{dt}}}$ , se saca:

siendo ${\textstyle {q}_{i}}$ el flujo calculado en el segmento indicado, ${\textstyle {q}_{(i+1)}}$ el flujo calculado en el segmento siguiente y ${\textstyle {C}_{i}}$ la complianza en ese segmento.

b) Ecuaciones de presión y flujo en segmentos terminales

El segmento de red arterial que se utiliza en este trabajo se divide en un número finito de vasos que incluyen arteriolas y capilares.

Debido al gran gasto informático y a la complejidad que esto presenta, se cierra la red arterial mediante una conexión en cascada de dos circuitos elásticos L y L-invertidos.

Como en el caso anterior, aplicando las leyes de Kirchhoff, se obtiene un sistema de primer orden de 3 ecuaciones diferenciales:

donde ${\textstyle {\overset {\cdot }{{q}_{in}}}}$, ${\textstyle {\overset {\cdot }{{q}_{out}}}}$y ${\textstyle {\overset {\cdot }{p\,}}}$son las variables a determinar. En cada segmento terminal, ${\textstyle {p}_{in\,}}$es la presión obtenida en el segmento anterior y ${\textstyle {p}_{out}\,}$ mantiene un valor constante de 15 mmHg que representa la presión venosa media que se usa para calcular el flujo de salida en el límite, ${\textstyle {R}_{b}}$es la resistencia terminal fija obtenida de la literatura [4].

Figura 15. Representación esquemática de la analogía eléctrica en un vaso terminal que fluye hacia capilares

3.2.1.2. Modelo viscoelástico

3.2.1.2.1 Parámetros de flujo y presión

a) Ecuaciones de presión y flujo en segmentos no terminales

El sistema visco elástico presenta las mismas características que la anterior en segmentos no terminales, es decir, las incógnitas son igualmente, ${\textstyle {\overset {\cdot }{q}}}$ y ${\textstyle {\overset {\cdot }{p}}}$ y se podrán obtener partiendo de las mismas características anatómicas de las arterias, como su longitud, diámetro, espesor, etc. Así pues, el flujo que atraviesa el vaso significativo se podrá calcular con:

siendo igualmente, ${\textstyle {p}_{(i-1)}}$ la presión de salida del segmento anterior, ${\textstyle {p}_{i}}$ la presión, ${\textstyle {R}_{i}}$ la resistencia, ${\textstyle {L}_{i}}$ la inductancia y ${\textstyle {q}_{i}}$ el flujo requerido.

No obstante, antes de definir la ecuación de la presión, es importante conocer la mecánica del comportamiento del vaso, puesto que al contrario que antes, no es solamente elástico, sino que presenta un comportamiento visco elástico en su naturaleza dentro del S.C. Esta característica de la pared puede ser representada mediante el modelo de Voigt [12].

La representación mecánica de este modelo se interpreta mediante un amortiguador (η) y un muelle ( ${\textstyle \epsilon }$ ) que sigue la ley de Hooke⁶ conectado en paralelo al amortiguador, que representa un sólido que sufre una deformación viscoplástica reversible. Bajo la aplicación de una tensión constante, el material se deforma a un ritmo progresivamente lento. Una vez se le retira esa tensión exterior, el material se relaja hasta su estado no deformado original.

Se define matemáticamente por una ecuación que relaciona el comportamiento elástico y viscoso del material:

donde ${\textstyle \sigma }$ es la tensión a la que está sometido el material, E es el módulo de Young, ${\textstyle \epsilon }$ la deformación que se le aplica, ${\textstyle \eta }$ la viscosidad del material, y ${\textstyle {\frac {d\epsilon (t)}{dt}}}$ es la derivada en función del tiempo de la tensión.

Con esto, la analogía eléctrica del modelo de Voigt se puede obtener considerando la resistencia R y la capacitancia C en serie [12]. Aplicando la ley de Kirchhoff sobre voltaje y corriente, se tiene:

donde ${\textstyle V}$ es voltaje, C la capacitancia, ${\textstyle Q\,\,}$ la carga, R la resistencia y ${\textstyle {\frac {dQ}{dt}}}$ la corriente eléctrica I.

Figura 16. Representación esquemática mecánica y eléctrica equivalente del modelo de Voigt

En el caso de los vasos arteriales, la relación tensión-deformación o carga-voltaje se puede simplificar en términos de presión-volumen. La ecuación correspondiente explica la naturaleza visco elástica de la pared de los vasos siendo:

donde ${\textstyle p}$ es la presión, ${\textstyle C}$ es la complianza del vaso, ${\textstyle {Q}_{c}}$ indica el volumen de sangre, ${\textstyle {R}_{d}}$ es la resistencia que mide viscosidad del vaso y ${\textstyle {\frac {dQc}{dt}}}$= ${\textstyle {q}_{c}}$ es el flujo de sangre que atraviesa el vaso. Sabiendo que la incógnita ${\textstyle p}$ equivale a ${\textstyle p}$_out , tomando la derivada de (ec.13) se obtiene,

donde se debe cumplir que,

y substituyendo en (ec.10), (ec.15) y (ec.16) en la (ec.14) para el cálculo de la presión en segmentos no terminales, finalmente obtenemos que,

que queda para la formula general,

Hay que tener en cuenta que ${\textstyle {R}_{di}}$ es la resistencia que determina la viscosidad en particular que cada vaso específico, no es fija. Por lo tanto, cada segmento tendrá la suya propia, establecida por [12].

Así pues, la analogía eléctrica que representa el vaso no terminal y el sistema viscoelástico presenta la siguiente forma (Figura 17).

Figura 17. Analogía eléctrica simplificada entre un vaso sanguíneo y su sistema viscoelástico

b) Ecuaciones de presión y flujo en segmentos terminales

En el caso visco elástico, en los segmentos terminales también se cierra la red arterial mediante una conexión en cascada de dos circuitos L y L invertidos.

Aplicando pues, las leyes de Kirchhoff, se obtienen un total de 3 ecuaciones diferenciales de primer orden,

como ecuación de presión.

Figura 18. Representación esquemática de la analogía eléctrica en un vaso terminal viscoelástico que fluye hacia capilares

La analogía eléctrica que representa el vaso no terminal y el sistema visco elástico presenta la siguiente forma.

Una vez habiendo tenido claro la estructura del circuito, se definen los parámetros eléctricos que lo componen, tales como el resistor R, la inductancia L, y el condensador C.

3.2.1.3. Parámetros eléctricos

3.2.1.3.1 Resistencia del flujo sanguíneo (R)

El parámetro R representa la fuerza de resistencia a la que el fluido está sometido a medida que avanza por el vaso de determinadas características, tales como el diámetro interno del vaso (d), su longitud (l), entre otros. Un parámetro importante a tener en cuenta en la R es la viscosidad del fluido, puesto que está presente cuando el fluido es estable u oscilatorio, disipando la energía de éste, aunque el flujo este acelerando o desacelerando. La viscosidad sanguínea que en este caso toma un valor constante de µ=0.004 Pa·s.

La resistencia total del flujo en un segmento concreto se modela mediante resistores y puede calcularse a partir de la ecuación de Poiseuille (ec.2), que como vimos en apartados anteriores, describe la relación entre la caída de presión, Δp y el flujo sanguíneo constante y q a través de un vaso sanguíneo uniforme y rígido.

Presenta la siguiente ecuación:

donde aislamos R

donde ${\textstyle l}$ es la longitud del vaso, ${\textstyle \mathrm {\mu } }$ la viscosidad sanguínea y ${\textstyle r}$, el radio interno del vaso. La ecuación presente también es conocida como la ley hidrodinámica de Ohm. Generalmente, R se obtiene de la ley de Poiseuille, pero en la práctica esta expresión se denomina ley hidrodinámica de Ohm. Las unidades que se usará con el resistor será en Pa·s·m^-3.

3.2.1.3.2 Complianza del vaso sanguíneo (C)

Cuando el flujo atraviesa un vaso elástico, el volumen del vaso puede variar. La clave del efecto de capacitancia sobre el flujo en un vaso elástico es que afecta al volumen total del vaso, por lo tanto, la velocidad de flujo en la entrada puede no ser la misma que en la salida debido a la variación de presión. Así pues, la capacidad de la pared de un vaso sanguíneo para expandirse y contraerse con los cambios de presión se conoce como capacitancia o complianza del vaso, C. La relación entre presión-volumen se puede cuantificar de dos modos. Por una parte, con la capacitancia C y por otra, con la elasticidad, es decir, el módulo de la elasticidad, ${\textstyle {E}_{v}=}$${\displaystyle \,{\frac {1}{C}}}$.

La capacitancia se calcula usando la siguiente ecuación, donde ${\textstyle \Delta V}$ es el cambio de volumen de sangre, y ${\textstyle \Delta p}$ la variación de presión en el segmento, quedando como:

donde añadimos las propiedades de los vasos sanguíneos siendo ${\textstyle E}$ el módulo de Young y ${\textstyle h}$ el grosor de la pared del vaso, tal que:

Las unidades que se usarán con el condensador serán en m³·Pa^-1.

3.2.1.3.3 Inercia de la sangre (L)

En cada latido del corazón, el fluido se acelera y se desacelera en el tiempo de un modo no inmediato debido a su inercia, L.

Dicho parámetro relaciona la caída de presión : ${\textstyle \Delta p}$ con el cambio de caudal ${\textstyle q}$ en un segmento, que se expresa mediante: ${\textstyle \Delta p=}$${\displaystyle L{\frac {dq}{dt}}}$.

Del mismo modo, el valor numérico de L, se puede obtener utilizando las propiedades de los vasos sanguíneos, aplicando la ley de Newton relativa al movimiento que implica, fuerza (F), masa (m) y aceleración donde ${\textstyle a}$= ${\textstyle {\frac {dv}{dt}}}$ y siendo ${\textstyle F}$ =m ${\textstyle {\frac {dv}{dt}}}$.

Partiendo de ahí, se aplica a un vaso de longitud, ${\textstyle l}$:

donde ${\textstyle A}$ es el área de sección transversal del vaso cilíndrico, ${\textstyle \rho =1050}$ kg/m³ es la densidad de la sangre y V es el volumen del fluido que atraviesa ese segmento. Aplicando la definición de flujo volumétrico ${\textstyle q=}$${\displaystyle v\cdot A}$, donde ${\textstyle v}$ es la velocidad en ese tramo, se deriva y se obtiene:

aislamos el diferencial

Así, substituyendo ${\textstyle F}$ y ${\textstyle m}$ en la ecuación tal que,

que como se ha mencionado, L es la relación entre ${\textstyle \Delta p}$ y el caudal, quedando finalmente como:

Las unidades utilizadas con el inductor serán en Kg·m^-4.

3.2.1.4. Condiciones del modelo

3.2.1.4.1 Condiciones de bifurcación

La descomposición de segmentos vasculares requiere tener en cuenta la relación entre varios de ellos dentro de una misma red arterial, es decir, la bifurcación de vasos. Por ello, se definen también condiciones de bifurcación tanto para las presiones como para los flujos, quedando la estructura y las ecuaciones algebraicas tal que:

en presiones tanto de fusión como divergencia del vaso.

Figura 19. Representación esquemática de bifurcación de la red arterial

Estas ecuaciones se obtienen de la conservación de la masa y del momento, que se traduce en que, la presión siempre será constante y el caudal del flujo debe conservarse en una bifurcación.

Los flujos divergentes se dan a cabo durante la bifurcación de los vasos sanguíneos. De acuerdo a la conservación de masa mencionada, el flujo del nodo 1 sería igual a la suma del flujo del nodo 2 y el nodo 3 (${\textstyle {q}_{1}=}$${\displaystyle {q}_{3}+\,{q}_{2}}$), siendo la presión constante en todo el segmento, es decir, la presión de salida en el nodo 1 ha de ser igual a la de entrada del nodo 2 y 3. Del mismo modo ocurre con la conservación de masa en zonas donde el flujo se fusiona, el flujo del nodo 3 ha de ser la suma del flujo del nodo 1 y 2 (${\textstyle {q}_{3}=}$${\displaystyle {q}_{1}+\,{q}_{2}}$) y, de acuerdo con la ley de conservación del momento, para este caso la presión total también permanece continua en todos los nodos.

3.2.1.4.2 Condiciones de salida

La longitud total de los vasos sanguíneos en el sistema cardiovascular humano de una persona adulta es aproximadamente de unos 95.000 kilómetros. Con estas cifras, es prácticamente imposible modelar decenas de miles de capilares de arteriolas y capilares del sistema circulatorio.

Por eso, una solución viable existente es la de agrupar todos los nodos posteriores que componen el sistema arterial, en segmentos terminales. Los capilares son considerados vasos de resistencia, de modo que se usa una resistencia terminal R_b para presentar el efecto acumulativo de la microcirculación más allá de los segmentos terminales. Este resistor será constante y presentará un valor de 1.89E+09 Pa·s·m^-3[4].

Además, se aplica una condición de contorno de presión constante, p_out, conocida como condición de frontera de Dirichlet⁷, para la resolución de las ecuaciones diferenciales. Entonces, el valor de p_out –que se trata de la presión venosa media para calcular el flujo de salida en el límite p_out -será de 15 mmHg (2000 Pa) [12]. Este valor, junto con el circuito terminal, se usó igualmente para los cálculos de las ecuaciones diferenciales, aunque el segmento no fuese realmente terminal en el cuerpo humano.

3.2.1.5. Resolución numérica

Para que se pueda resolver el modelo matemático hay que tener en cuenta que el número total de ecuaciones ha de corresponder al número total de variables desconocidas puesto que ha de ser un sistema determinado. Tanto el modelo elástico como el visco elástico, hay dos ecuaciones diferenciales por cada segmento arterial representado no terminal: aorta torácica, lágrima proximal, lágrima distal, el FL y el TL, haciendo un total de 8 ecuaciones diferenciales. Juntamente con el segmento terminal, que añade otra ecuación diferencial más referente a la aorta torácica.

Además, como para cada bifurcación se deben cumplir las leyes de cantidad y momento y conservación de masa, se añaden al sistema 4 ecuaciones algebraicas más, haciendo un total de 13 ecuaciones.

Debido a que este sistema contiene tanto ecuaciones diferenciales como ordinarias, se obtiene un sistema de ecuaciones algebraicas diferenciales (DAE). Hay que remarcar que, tanto el modelo elástico como el visco elástico, tienen el mismo número de ecuaciones diferenciales y algebraicas, pero con la diferencia, que como se mencionó en el Apartado 3.2.1.2, la presión arterial en la salida de cada segmento para el modelo elástico, solo depende de la resistencia terminal, pero para el modelo viscoso, depende de la complianza del vaso (C) y de la resistencia de la arteria.

El número total de incógnitas será de 13 (presión y caudal de entrada y de salida en cada segmento no terminal, y caudal de entrada, de salida y presión para el segmento terminal), pero debido a que se conoce el caudal en la entrada del sistema en la aorta ascendente –se explicará en apartados siguientes-se obtiene un sistema de 12 ecuaciones con 12 incógnitas.

Para poder resolver este sistema de ecuaciones algebraicas diferenciales se necesita definir una matriz de masa. Esta matriz de masa tiene que ser diagonal constante con 1 para las variables diferenciales y 0 para las variables algebraicas, representándose: ${\displaystyle M(t,y)y'=f(t,y)}$, con ${\displaystyle M(t,y)}$ definida.

Por ejemplo, un sistema DAE con 2 ecuaciones diferenciales y una variable algebraica se expresaría como:

Por lo tanto, para el modelo desarrollado, la matriz de masa será de [12x8], de los cuales 1 son los primeros 8 elementos de la diagonal que serán igual a 1 y el resto igual a 0 correspondiendo a las ecuaciones algebraicas. Una vez definido la matriz de masa y las ecuaciones diferenciales y algebraicas, se utiliza el método de diferenciación numérica – ode15s – implementado en Matlab [21].

Destacar que este método es muy eficiente para sistemas en los que existen variación en la solución, como el presente. Es decir, el caudal de entrada posee zonas en las que el flujo varía lentamente(diástole) y zonas donde el flujo varía rápidamente (sístole). El modelo también se implementó mediante el lenguaje de programación Matlab.

(⁶) La ley de Hooke establece que el alargamiento unitario que experimenta un material elástico es directamente proporcional a la fuerza que se le aplica. La forma más común de representarla es mediante un muelle, donde se relaciona la fuerza F ejercida por el muelle con la elongación δ provocada por la fuerza aplicada, quedando F=-k δ, donde k es la constante elástica del muelle.

(⁷) La condición de frontera de Dirichlet se emplea cuando en una ecuación diferencial ordinaria o en una derivada parcial, se le especifican los valores de la solución que necesita la frontera del dominio

4. Resolución de los modelos

En esta sección se aplica la herramienta numérica explicada en apartados anteriores, a un modelo experimental. Primeramente, se modificará levemente la elasticidad del modelo visco elástico para ver el comportamiento de la presión, se relacionará cualitativamente este modelo con una versión experimental, y finalmente se compararán cualitativamente y cuantitativamente los modelos visco elástico y elástico.

Los datos que se utilizan en este caso son clínicamente relevantes, pues las resistencias periféricas, los condensadores y las inductancias y las propiedades geométricas son calculadas a partir de datos reales [9] y, por lo tanto, tienen valores fisiológicos realistas. Se han considerado los siguientes casos significativos:

• Sección aorta en sistema elástico

• Disección aórtica en sistema elástico

• Disección aórtica en sistema visco elástico

4.1. Modelo elástico de la arteria aorta

4.1.1. Parámetros de entrada

Se exponen a continuación, las condiciones de entrada al sistema para poder llevar a cabo la simulación, como será la recreación del primer flujo y presión de entrada.

Como primera validación del modelo utilizado se ha de asumir que:

• El ciclo cardiaco empieza en sístole.

• El período de sístole es 2/5 respecto al periodo total del ciclo cardíaco.

• La complianza arterial, la resistencia periférica y la inercia son representadas por un condensador, una resistencia y una bobina respectivamente.

• El flujo inicial que se obtendrá será el calculado en la aorta ascendente.

Con ello, se establecen los datos fisiológicos generales obtenidos en que la densidad de la sangre, ρ, será ρ=1050 ${\textstyle {\frac {Kg}{{m}^{3}}}}$ y la viscosidad de la sangre µ, será µ=4 mPa·s.

Para la entrada de flujo del primer latido (qi), se impone una función de caudal periódico. Cada período consiste en una onda sinusoidal, para reproducir la fase de la sístole y la fase de la diástole [20].

El flujo de sangre que emite el ventrículo hacia la aorta será representado durante el ciclo cardiaco completo-siendo este de 0.8 s- como I(t) con amplitud ${\textstyle {I}_{0}}$ durante la sístole, pues durante la diástole, los ventrículos están relajados y no hay flujo emergente, por lo tanto, I(t)=0. Sin embargo, la contracción ventricular durante el primer periodo, la sangre sí que fluye hacia la aorta definiendo la siguiente ecuación:

donde t es tiempo en segundos, ${\textstyle {T}_{c}}$ es el periodo de ciclo completo en segundos, ${\textstyle {T}_{s}}$ es el período de sístole en segundos, y ${\textstyle mod\left(t,{T}_{c}\right)}$ representa el vector ${\textstyle t}$ divido entre ${\textstyle {T}_{c}}$. El parámetro ${\textstyle {T}_{s}}$ se ha determina como 2/5 ${\textstyle {T}_{c}}$ de acuerdo a los parámetros dinámicos del ciclo cardiaco. Se establece también, un flujo entrante en el primer latido de 70 cm³.

Usando ésta información obtenemos la constante ${\textstyle {I}_{0}}$ en:

Gráfica 1. Primer latido simulado del flujo en la entrada de la arteria aorta

Usando el software MATLAB y el código desarrollado, se obtiene que ${\textstyle {I}_{0}=}$${\displaystyle 297,77\,ml}$ en el bombeo de inyección durante la sístole. Ese será el vector inicial utilizado, q_i.

Una vez se obtiene el vector de entrada, se debe obtener la presión primera para poner iniciar. Se representa a continuación una pequeña simulación de la presión inicial con la que se emprende el circuito obtenida a partir del flujo de entrada. Se obtienen un total de 20 ciclos cardiacos, de 0,8 segundos cada uno.

Gráfica 2. Primer parámetro de presión simulado en la entrada de la arteria aorta

Los valores obtenidos en p_in se consideran aceptables pues como se aprecia en la Figura 12, se obtienen valores máximos de ± 120 mmHg como presión sistólica, y valores mínimos de ± 80 mmHg como presión diastólica. Los valores fisiológicos reales utilizados de la aorta ascendente obtenidos de para esta iniciación son:

A partir de estos datos obtenemos los valores eléctricos con las formulas mencionadas en el capítulo tres:

4.1.2. Inicio de la aorta

Partiendo del primer latido del apartado anterior, se procederá a la simulación del comportamiento de presión y flujo en el primer segmento aórtico mediante el modelo elástico. Dicha sección, está compuesta por la aorta ascendente, el arco aórtico y la parte superior de la aorta descendente, la aorta torácica, respectivamente. Esto significa que se empleará un mini modelo para cada parte descrita, es decir, uno para la aorta ascendente con sus datos fisiológicos, otro para el arco aórtico, y otro para la aorta torácica.

Aunque anatómicamente sea un segmento curvado, al hacerlo sobre un modelo 0-D, se tomarán como varios cilindros rectos unidos entre sí, sin tener en cuenta los factores geométricos.

Los datos fisiológicos del modelo serán los siguientes:

La analogía eléctrica del primer segmento queda tal como muestra la Figura 20.

Figura 20. Analogía eléctrica simplificada de la aorta Fuente: https://es.wikipedia.org/wiki/Arco_a%C3%B3rtico

Se determina en cada segmento, la presión y el flujo del vaso con los datos fisiológicos expuestos en la página anterior. Cabe destacar que, q_in y p_in, son los vectores calculados al inicio de la sección y que la presión y el flujo que se calculará en el tercer segmento de la aorta torácica, es el que se empleará para la siguiente sección sobre la disección aórtica.

No obstante, hay que tener en cuenta que, en una persona, esta arteria, solo en los inicios presenta varias bifurcaciones, en la aorta ascendente hay dos, la coronaria derecha y coronaria izquierda, y en arco hay 4, subclavia derecha e izquierda y carótida derecha e izquierda.

Aún y así, para representar este segmento, se ha tomado como un único tubo a representar, el conductor principal, despreciando las bifurcaciones, teniendo en cuenta que la presión sanguínea será la misma.

Lo que se pretende establecer con este modelo, es que, a parte de la funcionalidad del primer latido, la simulación de sección anatómica arterial siguiente sea igual de funcional y recree con la mayor veracidad posible–dentro de lo que los modelos 0-D permiten hacer y representar- las características hemodinámicas de este segmento arterial tan fundamental, como es la aorta.

A continuación, se muestra una gráfica esquemática (Figura 21) [23] de cómo se propaga el flujo (m³/s) a través de la aorta durante la primera sístole, en un período de 0.3 segundos aproximadamente–período de eyección-. Siendo un flujo ínfimo al inicio –unos 4.753E-05 m³/s-, y va aumentando a medida que avanza la primera fase, alcanzando casi los 0.00025⁸ m³/s.

Figura 21. Representación del flujo aórtico durante la sístole [30]

Tras graficar, se observa como el flujo a través de esta primera sección no varía prácticamente, puesto que el caudal de entrada ha de igual siempre al de salida, es decir, la cantidad de flujo que sale y entra y atraviesa las arterias ha de mantenerse constante, para asegurar un buen funcionamiento de las funciones corporales y asegurar que el riego sanguíneo –en situaciones sin patologías- llegue a todos los extremos sin excepciones.

El pequeño desfase se podría deber a la propagación de onda a través de los tubos. El flujo máximo que se aprecia en la Figura 21, se ve representado en su punto máximo –presión sistólica- en la Gráfica 3.

Gráfica 3. Curvas de flujo del primer segmento de la aorta, indicando el punto máximo

Con la presión sanguínea ocurre el mismo fenómeno, debe mantenerse aproximadamente constante a lo largo de todos los vasos arteriales, desde que sale del ventrículo izquierdo, hacia la aorta ascendente, pasando por todo el tronco aórtico y sus ramificaciones hasta los órganos, hasta llegar a las extremidades y los capilares. La propagación de presión sobre todo el sistema arterial durante un latido entero (0.8 segundos aproximadamente) incluyendo sístole y diástole, asumiendo que el latido empieza en sístole, se puede observar simplificado en la Figura 22.

Figura 22. Propagación de la presión sanguínea durante la sístole [30]

Se puede apreciar como en ambas fases la presión de la sangre se muestra prácticamente constante a través de todo el sistema arterial. En la fase de sístole, el flujo es mayor, con comportamiento creciente, llegando a casi 16500 Pascales (120 mmHg) en los puntos más altos, puesto que es el período de eyección sanguínea y de contracción de los ventrículos (Figura 22). Durante la fase de diástole, donde se produce la dilatación ventricular, la presión sanguínea baja hasta estar entre los 12000 y 10000 Pascales (90-75 mmHg, respectivamente) en los puntos más bajos como las extremidades (Figura 23).

Figura 23. Propagación sanguínea durante la diástole [30]

En la Gráfica 4, que se anexa a continuación perteneciente a la primera sección de la aorta, se comprueba que la presión de sangre no varía prácticamente en ninguno de los tres puntos. El desfase presente puede deber a la propagación de la onda de presión a lo largo del tubo.

Gráfica 4. Presión en los primeros segmentos de la aorta

(⁸) Hay que añadir que los valores representados en esta figura, son representativos. Son valores medios obtenidos con datos fisiológicos específicos. La cantidad de caudal de una persona puede variar en función de sus características y su actividad física en algún momento determinado. Lo mismo se aplica a los valores de presión.

4.2. Disección aórtica

4.2.1. Patología

Hoy en día, las patologías aorticas representa un subgrupo muy importante dentro de las enfermedades cardiovasculares y, aunque su prevalencia es limitada, se asocian con una mortalidad muy alta (>50%) [12]. A pesar de las técnicas de diagnóstico actuales y las terapias mejoradas, el manejo y el seguimiento de las disecciones aorticas continúan siendo un desafío en la práctica clínica. La disección aórtica clásica empieza con la formación de una pequeña incisión llamada lágrima –tear- en la capa íntima de la arteria, que expone a la capa media siguiente a la presión pulsátil del flujo sanguíneo (Figura 24). Esto conlleva a una escisión en dirección longitudinal de ésta capa media a lo largo de la pared aórtica causando la disección. Este proceso de disección se extiende impulsado por la fuerza hacia delante del flujo sanguíneo aórtico.

Por ello, el lumen quedará dividido en dos partes, el canal principal con el verdadero lumen, TL, y el canal falso, fruto de esta patología, FL. En la mayoría de los casos, el TL y el FL se comunican a través del sitio de entrada y de salida. –véase anatomía de una artería en el Apartado 2.4.2.

Figura 24. Disección de la aorta ascendente (izquierda). Imagen de la aorta descendente en una resonancia magnética (centro) y en tomografía (derecha). TL: verdadero lumen. FL: Falso lumen Fuente: Validation of numerical flow simulations against in vitro phantom measurements in different type B aortic dissection scenarios, 2013

Los pacientes que lo padecen presentan una alta mortalidad a medio/largo plazo, debido a la dilatación progresiva del tejido de la aorta y su posterior ruptura. En la práctica actual, la predicción del resultado se basa principalmente en la determinación del diámetro aórtico máximo.

Sin embargo, se ha visto que este parámetro no es un determinante confiable en la progresión y en la ruptura de la aorta, pues variables hemodinámicas como la presión, las condiciones de flujo, la elasticidad de la pared y factores geométricos como la forma de curvatura de la aorta, la sección de comunicación entre falso lumen y verdadero lumen y el tamaño de las lágrimas juegan un papel importante e influyen directamente en el progreso de dilatación del tejido [17].

Dentro de esta patología podemos encontrar dos clasificaciones en función de la extensión en la aorta ascendente:

• Tipo A: Se origina en la aorta ascendente y se puede extender hacia todo el tronco de la aorta o se detiene en el arco aórtico.

• Tipo B: Se origina en la aorta descendente y rara vez se extiende. No involucra la sección ascendente.

Figura 25. Tipos de disección aórtica Fuente:https://es.wikipedia.org/wiki/Disecci%C3%B3n_a%C3%B3rtica

Así pues, el objetivo será la elaboración de un circuito eléctrico en dimensión 0-D análogo a la estructura arterial para simular el comportamiento del flujo sanguíneo dentro de una disección aórtica. Los datos utilizados para recrear las dimensiones de la arteria y sus extremos son propuestos por Alastruey [4].

Se recreará un escenario de tres secciones principalmente, donde se hará hincapié en la sección central perteneciente a la disección aórtica. La simulación se iniciará en la aorta torácica, seguida de una parte aorta descendente donde se hará la disección aórtica, y se finalizará con la aorta abdominal.

Dentro de la disección aórtica, se distinguirá el verdadero lumen (TL) y el falso lumen (FL) alternando diferentes tamaños de las lágrimas simulando hipotéticas situaciones. En cada escenario se determinará la presión en el TL y el FL en la sección proximal y distal tanto en el modelo elástico como en el visco elástico para ver la diferencia del comportamiento de la onda entre ambas y se contrastará desde un punto de vista fisiológico con otros modelos consultados en la literatura [9].

Cada lágrima se representará mediante un resistor con un valor particular. La sección de la aorta abdominal, al ser nuestro segmento terminal se recreará con el modelo de ecuaciones de flujo y presión para segmentos terminales explicados en el capítulo tres.

Los datos fisiológicos de la arteria abdominal y torácica serán los mismos para los tres escenarios y los dos modelos. Presentan los siguientes valores:

El modelo a representar se localiza y se presenta de la siguiente manera (Figura 26).

Figura 26. Analogía virtual en 3D de la disección aórtica del modelo representado

Así pues, los escenarios trabajados serán,

1. Escenario primero

Figura 27. Representación de la arteria. Primer escenario

Las condiciones de la segunda situación serán las siguientes: se determinará la presión en el falso y verdadero lumen en la sección proximal y lo mismo para la sección distal. Los diámetros de las lágrimas en este caso serán del mismo tamaño, siendo de 10 mm la primera lágrima y 10 mm la segunda.

Los datos fisiológicos con los que se simula el modelo tanto elástico como visco elástico del primer escenario y de los demás escenarios de [9] presentan los siguientes valores (Tabla 9).

2. Escenario segundo

Figura 28. Representación de la arteria. Segundo escenario

Las condiciones de la segunda situación serán las siguientes; se determinará la presión en el falso y verdadero lumen en la sección proximal y lo mismo para la sección distal. Los diámetros de las lágrimas en este caso serán de 10 mm la primera lágrima y 4 mm la segunda.

Los datos fisiológicos con los que se simula el modelo tanto elástico como visco elástico del escenario segundo presentan los siguientes valores:

3. Escenario tercero

Figura 29. Esquema sobre la arteria en el tercer escenario

Las condiciones de la tercera situación serán las siguientes; se determinará la presión en el falso y verdadero lumen en la sección proximal y lo mismo para la sección distal. Los diámetros de las lágrimas en este caso serán de 4 mm la primera lágrima y 0 mm la segunda, es decir, no habrá segunda lágrima.

Los datos fisiológicos con los que se simula el modelo tanto elástico como visco elástico del escenario primero presentan los siguientes valores:

4.2.2. Disección aórtica con modelo elástico

El modelo eléctrico que representará el segmento aórtico patológico se presenta en la Figura 30.

Figura 30. Modelo elástico de la disección aórtica con los parámetros principales a determinar

Donde las incógnitas principales serán las presiones en las sección proximal y distal tales como P_pTL que equivale a la presión proximal del verdadero lumen (TL), P_pFL que es la presión proximal a través del falso lumen (FL), P_dTL como presión del TL en la sección distal y finalmente, P_dFL que representa la presión distal en el FL.

El modelo elástico recrea el comportamiento de los vasos como cilindros, rectos y rígidos. Se divide en cuatro segmentos compuestos por circuitos RLC que pretenden simular una disección aórtica en una sección de la aorta descendente. El primer segmento viene con los valores de q_in y p_in preestablecidos que se explica en el Apartado 4.1.1.

El modelo que pertenece a esta patología se le añaden dos resistencias en paralelo para imitar el comportamiento de las lágrimas. El cambio de diámetro de las lágrimas se verá reflejado en un cambio en el valor de las resistencias que lo definen –véase apartado anterior-.

La resolución numérica empleada es la mencionada en el Apartado 3.2.1.5, que trata de definir un sistema de ecuaciones diferenciales de primer orden para obtener las variables requeridas. Las ecuaciones que las determinan se presentan de la siguiente manera,

La presión restante se calculará mediante la ecuación algebraica,

El parámetro ${\textstyle {R}_{DT}}$ es la resistencia específica para la lágrima en la sección distal. En estas ecuaciones, cada presión viene determinada indirectamente por los valores de flujo calculados en cada segmento (Figura 29).

4.2.3. Disección aórtica con modelo visco elástico

A diferencia del modelo elástico, el modelo visco elástico añade una resistencia extra (R_d) en serie con el condensador para cuantificar la viscoelasticidad del vaso. El modelo eléctrico se representa en el siguiente esquema.

Figura 31. Modelo viscoelástico de la disección aórtica con los parámetros principales a determinar

De igual modo que en la Figura 30, las incógnitas principales serán las presiones en las sección proximal y distal, P_pTL, P_pFL que son las presiones proximales a través del verdadero (TL) y falso lumen (FL) respectivamente, P_dTL como presión del TL en la sección distal y, P_dFL que representa la presión distal en el FL.

El primer segmento también viene con los valores de q_in y p_in preestablecidos obtenidos del Apartado 4.1.1.

A parte de las dos resistencias que simulan las lágrimas, presenta la resistencia R_d, añadida en serie al condensador, para recrear un modelo de orden reducido con propiedades visco elásticas. La importancia de añadir este resistor extra se encuentra en el hecho de que el comportamiento de la onda de presión que se obtiene, sea más real y que capture algunos comportamientos fisiológicos, como, por ejemplo, la muesca dicrota.

También, al suponer que la pared será más distensible, el pico máximo de la onda será más tardío respecto al elástico. El valor de esta resistencia se obtiene de [12] y varía según la sección, es decir, la aorta torácica tiene una, la abdominal otra, y la parte de aorta descendente en la que se aplica el modelo de la disección aórtica, tanto para el FL como el TL tendrán otra –en este caso se llama R_d as la perteneciente a la de la disección aórtica, y R_d ab, la que se emplea para la aorta abdominal -.

Las ecuaciones, también resultas mediante el método mencionado en el Apartado 3.2.1.5, que rigen la presión especificada en este modelo, se expresan como,

y como ecuación algebraica restante aplicando las leyes de Kirchhoff,

El parámetro ${\textstyle {R}_{DT}}$ es la resistencia específica para la lágrima en la sección distal. En estas ecuaciones, cada presión viene determinada indirectamente por los valores de flujo calculados en cada segmento (Figura 30).

4.3. Resultados

Con los modelos mencionados anteriormente, se exponen los resultados obtenidos tras su evaluación en distintas circunstancias. Primeramente, se representa el modelo visco elástico con varias modificaciones del módulo de Young, para ver cómo evoluciona la onda de presión en función del valor de su rigidez. Por otro lado, se someterá a este mismo modelo a una comparación cualitativa con el modelo experimental [9] y finalmente se mostrarán los modelos elástico y visco elástico conjuntamente, para resaltar las diferentes y similitudes, y remarcar la importancia de la variación de la forma de onda que presenta cada modelo. A esta sección se le añadirá una pequeña comparación cualitativa, mediante el cálculo del coeficiente de variación entre ambas ondas.

4.3.1. Variación del módulo de Young del modelo visco elástico

I. Escenario primero ${\displaystyle \phi p=10}$ mm, ${\displaystyle \phi d=10}$ mm)

El valor E de referenia es de 1070 Pascales en el TL. Se ha expuesto el sistema al doble, con 2140 Pascales y a la mitad, 530 Pascales. Para el falso lumen, la E de referencia es de 2490 Pascales, el máximo de 4980 Pascales, y el mínimo de 1250 Pascales, tanto para la sección proximal y distal.

a. Sección proximal

Gráfica 5. Presión en el FL variando el módulo de Young. Sección proximal

Gráfica 6. Presión en el TL variando el módulo de Young. Sección proximal

b. Sección distal

Gráfica 7. Presión en el FL variando el módulo de Young. Sección distal

Gráfica 8. Presión en el TL variando el módulo de Young. Sección distal

II. Escenario segundo ${\displaystyle \phi p=10}$ mm, ${\displaystyle \phi d=4}$ mm)

El valor E de referencia es de 1070 Pascales. Se ha expuesto el sistema al doble, con 2140 Pascales y a la mitad, 530 Pascales, en el caso del verdadero lumen. Para el falso lumen, la E de referencia es de 1510 Pascales, el máximo de 3020 Pascales, y el mínimo de 750 Pascales, para la sección proximal y distal.

c. Sección proximal

Gráfica 9. Presión en el TL variando el módulo de Young. Seccion distal

Gráfica 10. Presión en el FL variando el módulo de Young. Seccion distal

d. Sección distal

Gráfica 11. Presión en el FL variando el módulo de Young. Sección distal

Gráfica 12. Presión en el TL variando el módulo de Young. Sección distal

III. Escenario tercero ${\displaystyle \phi p=4}$ mm, ${\displaystyle \phi d=0}$ mm)

El valor E de referenia es de 1070 Pascales. Se ha expuesto el sistema al doble, con 2140 Pascales y a la mitad, 530 Pascales, en el caso del verdadero lumen. Para el falso lumen, la E de referencia es de 3910 Pascales, el máximo de 7640 Pascales, y el mínimo de 1910 Pascales, para la sección proximal y distal.

e. Sección proximal

Gráfica 13. Presión en el FL variando el módulo de Young. Sección proximal

Gráfica 14. Presión en el TL variando el módulo de Young. Sección distal

f. Sección distal

Gráfica 15. Presión en el FL variando el módulo de Young. Sección distal

Gráfica 16. Presión en el TL variando el módulo de Young. Sección distal

El objetivo de alterar el módulo de Young (E) predefinido por otros valores, es la observación del efecto de la elasticidad en la pared arterial en el verdadero y falso lumen y ver las variaciones que se producen en las curvas. Tras observar las gráficas, se puede definir que, un aumento del valor del módulo elástico en el vaso, se traduce a un vaso con mayor rigidez. Este hecho se plasma en un pico de onda más elevado y una presión ligeramente mayor. A la inversa, si E disminuye, el vaso es más deformable y la onda tiene una presión menor.

Esto sucede porque, el comportamiento visco elástico de los vasos proporciona un mecanismo de almacenamiento de energía, que se almacena durante la eyección sanguínea. La energía se disipa por las pareces del vaso. Durante la relajación, la energía almacenada empuja la sangre hacia los vasos siguientes, hecho que provoca un flujo sanguíneo pulsátil, una onda que lo representa más extensa y una ligera disminución en el valor de la presión. La onda que se obtiene en este caso, es más progresiva y con un punto máximo inferior.

Cabe destacar que, la diferencia de la forma de onda entre el TL y el FL se debe al caudal, es decir, el TL presenta un caudal mayor respecto al segundo. Esto puede verse reflejado puesto que en las ondas del TL se aprecian mejor características hemodinámicas destacables, como la muesca dicrota. La onda del FL es más suave y redondeada porque el flujo que lo atraviesa es mucho menor y cuesta percibir alguna característica destacable.

4.3.2. Comparación del modelo visco elástico versus modelo experimental

En esta sección se hará una comparación cualitativa del modelo visco elástico cogiendo como referencia un modelo experimental [9]. El objetivo de este apartado es comparar los cambios relativos de presión que se producen entre el TL y el FL del modelo viscoso desarrollado contra dicho modelo. No se han podido comprobar directamente ambos modelos, puesto que se desconocían ciertos datos del modelo experimental.

El modelo experimental se realizó in vitro –los datos con los que trabajan también son experimentales- sobre un sistema elástico, recreando una disección aórtica con un phantom de látex. Se espera observar similitudes significativas como la propagación de onda, valores parecidos de presiones en el FL y TL y cierta similitud con el comportamiento de la onda, teniendo en cuenta que la comparación se realiza desde un modelo visco elástico y la forma de onda que éste presenta.

Algunos de los datos con los se ha recreado el sistema visco elástico, -módulo de Young para el TL y el FL en cada escenario y las resistencias que representan las lágrimas-fueron recogidas de este modelo experimental. Sin embargo, los parámetros de entrada que se usaron en cada caso eran diferentes, puesto que en el caso del modelo visco elástico del trabajo, q_in y p_in eran calculadas mediante una senoide, y en el modelo in vitro fueron obtenidas experimentalmente.

I. Escenario primero (${\displaystyle \phi p=10}$ mm, ${\displaystyle \phi d=10}$ mm)

a. Sección proximal

Gráfica 17. Curva de presión en modelo viscoelástico. sección proximal. Primer escenario

Gráfica 18. Curva de presión en el modelo elástico experimental. sección proximal. Primer escenario

La evaluación de la primera onda en la parte proximal, muestran claras similitudes. Aunque las ondas son ciertamente diferentes, se puede apreciar como ambas curvas empiezan bajo el mismo comportamiento. Empiezan alrededor de unos 80 mmHg, yendo la presión del TL en cabeza respecto al FL. En el período de crecimiento de onda, ambos casos se separan notablemente, volviéndose a unir en el punto máximo, llegando a la par a unos 120 mmHg.

En la onda visco elástica se muestra el cierre de la válvula aórtica –muesca dicrota-, mientras que, en el modelo elástico experimental, esa característica fisiológica no se representa. Al decrecer la onda, en ambos casos los valores de presión del TL y el FL son muy semejantes y se unifican en una sola línea.

b. Sección distal

Gráfica 19. Curva de presión en modelo viscoelástico. sección distal. Primer escenario

Gráfica 20. Curva de presión en modelo experimental. sección distal. Primer escenario

En esta sección distal, se aprecia primeramente como la onda de presión del FL de salida es mayor respecto al TL. Durante el crecimiento de la onda, las dos presiones se cruzan, llegando al punto máximo de presión en ambas curvas al mismo instante. –aproximadamente en el segundo 12,3-Durante el período de decrecimiento, la curva del FL se adelanta hasta llegar a su punto más bajo y muestra un área bajo la curva ligeramente más extensa que TL. Cabe decir que respecto el comportamiento general, no se observa apenas desfase de onda entre el FL y TL en ninguno de los dos modelos.

II. Escenario segundo (${\displaystyle \phi p=10}$ mm, ${\displaystyle \phi d=4}$ mm)

c. Sección proximal

Gráfica 21. Curva de presión en modelo viscoelástico. sección proximal. Segundo escenario

Gráfica 22. Curva de presión en modelo experimental. sección proximal. Segundo escenario

En el segundo escenario la curva del falso lumen inicia en un valor mayor respecto al verdadero lumen. Inmediatamente se cruzan, siendo el verdadero lumen el que aumenta más rápidamente. Al llegar casi al punto máximo de presión, se vuelven a cruzar intercambiando papeles. El falso lumen alcanza una presión mayor en ambos casos. Se aprecia suavemente un leve salto durante el decrecimiento de las ondas, que pertenece a la incisura dicrota. Se atisba un pequeño desfase de onda en ambos casos, con el falso lumen por delante respecto al verdadero lumen. Durante el decrecimiento perteneciente a la fase de sístole, el TL disminuye más rápidamente.

La variación entre los valores de presión entre el FL y el TL puede deberse simplemente a una cuestión numérica durante la resolución de las ecuaciones diferenciales que componen el sistema. Como se ha especificado en otros apartados, se representan las lágrimas con resistencias. Según el tamaño de la lágrima el valor de dicha resistencia varia. En este caso, como ambas lágrimas poseen diámetros diferentes, la resistencia con la que se trabajará también será diferente, y en este caso, el resistor que representa la lágrima de 10 mm es menor que la de 4 mm, quedando R₁₀ < R₄. Por ese motivo, como el cálculo de la presión del falso lumen depende de estas variables indirectamente, puede haber sido un factor que haya hecho que aumente la presión en esa sección.

d. Sección distal