Nuestros investigadores

Jorge Aramburu Montenegro

Departamento de Ingeniería Mecánica y Materiales
Escuela de Ingeniería (TECNUN). Universidad de Navarra
Líneas de investigación
Hemodinámica computacional
Índice H
4, (Scopus, 05/01/2021)
6, (Google Scholar, 05/01/2021)
4, (WoS, 05/01/2021)

Publicaciones científicas más recientes (desde 2010)

Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl (Autor de correspondencia); Rivas Nieto, Alejandro; et al.
ISSN 0929-8673  Vol. 27  Nº 10  2020  págs. 1600 - 1615
Radioembolization (RE) is a valuable treatment for liver cancer. It consists of administering radioactive microspheres via an intra-arterially placed catheter with the aim of lodging these microspheres, which are driven by the bloodstream, in the tumoral bed. Even though it is a safe treatment, some radiation-induced complications may arise. In trying to detect or solve the possible incidences that cause nontarget irradiation, simulating the particle-hemodynamics in hepatic arteries during RE via computational fluid dynamics (CFD) tools has become a valuable approach. This paper reviews the parameters that influence the outcome of RE that have been studied via numerical simulations. In this numerical approach, the outcome of RE is regarded as successful if particles reach the artery branches that feed tumor-bearing liver segments. Up to 10 parameters have been reviewed. The variation of each parameter actually alters the hemodynamics pattern in the vicinities of the catheter tip, and locally alters the incorporation of the particles into the bloodstream. Therefore, in general, the local influences of these parameters should result in global differences in terms of particle distribution in the hepatic artery branches. However, it has been observed that under some (qualitatively described) appropriate conditions where particles align with blood streamlines, the local influence resulting from a variation of a given parameter vanishes and no global differences are observed. Furth
Autores: Ramos González, Juan Carlos; Antón Remírez, Raúl; Aramburu Montenegro, Jorge; et al.
Revista: DYNA
ISSN 0012-7361  Vol. 95  2020  págs. 640-645
The air conditioning systems used in transport vehicles (railway passenger cars) are constrained by two issues: the consumed energy and the occupied space. In order to assess how both restrictions can be fulfilled a parametric analysis and an optimization have been carried out by means of Design of Experiments (DoE) techniques applied to a mathematical model of a real AC system. To evaluate the energy efficiency of the system four parameters have been used: the COP (Coefficient of Performance), the refrigeration power of the system and its dimensionless expression and the effectiveness of the evaporator. The occupied space of the AC system has been characterized through the dimensionless volume of the evaporator. The same compressor has been used in the analysis and the input parameters varied have been the dimensionless evaporator volume and four operating conditions: the evaporator inlet air temperature and mass flow rate, the condenser inlet air temperature and the air temperature increment in the condenser. Through a Central Composite Surface Response design, the regression equations of the four output energy efficiency parameters have been obtained in a wide range of operating conditions. Results show that the dimensionless refrigeration power is the best parameter because is independent of the operating conditions and only depends on the dimensionless volume of the evaporator: the lower the latter, the higher the former. An optimal dimensionless volume has been identified for the simultaneous optimization of the four output variables of the system that provides an increase of 26% in the dimensionless refrigeration power with respect to the baseline system.
Autores: Ortega Calderón, Julio; Antón Remírez, Raúl (Autor de correspondencia); Ramos González, Juan Carlos; et al.
ISSN 2040-7939  Vol. 36  Nº 6  2020  págs. e3337
In the last decades, the numerical studies on hemodynamics have become a valuable explorative scientific tool. The very first studies were done over idealized geometries, but as numerical methods and the power of computers have become more affordable, the studies tend to be patient specific. We apply the study to the numerical analysis of tumor-targeting during liver radioembolization (RE). RE is a treatment for liver cancer, and is performed by injecting radiolabeled microspheres via a catheter placed in the hepatic artery. The objective of the procedure is to maximize the release of radiolabeled microspheres into the tumor and avoid a healthy tissue damage. Idealized virtual arteries can serve as a generalist approach that permits to separately analyze the effect of a variable in the microsphere distribution with respect to others. However, it is important to use proper physiological boundary conditions (BCs). It is not obvious, the need to account for the effect of tortuosity when using an idealized virtual artery. We study the use of idealized geometry of a hepatic artery as a valid research tool, exploring the importance of using realistic spiral-flow inflow BC. By using a literature-based cancer scenario, we vary two parameters to analyze the microsphere distribution through the outlets of the geometry. The parameters varied are the type of microspheres injected and the microsphere injection velocity. The results with realistic inlet velocity profile showed that the par
Autores: Mariscal-Harana, Jorge; Charlton, Peter; Vennin, Samuel; et al.
ISSN 0363-6135  2020 
Central blood pressure (cBP) is a highly prognostic cardiovascular (CV) risk factor whose accurate, invasive assessment is costly and carries risks to patients. We developed and assessed novel algorithms for estimating cBP from non-invasive aortic haemodynamic data and a peripheral BP measurement. These algorithms were created using three blood flow models: the 2-and-3-element Windkessel (0-D) models and a one-dimensional (1-D) model of the thoracic aorta. We tested new and existing methods for estimating CV parameters (ejection time, outflow BP, arterial resistance, compliance, pulse wave velocity, characteristic impedance) required for the cBP algorithms, using 'virtual' subjects (n=19,646) for which reference CV parameters were known exactly. We then tested the cBP algorithms using 'virtual' subjects (n=4,064), for which reference cBP were available free-of-measurement error, and clinical datasets containing invasive (n=10) and non-invasive (n=171) reference cBP waves across a wide-range of CV conditions. The 1-D algorithm outperformed the 0-D algorithms when the aortic vascular geometry was available, achieving central systolic blood pressure (cSBP) errors ¿2.1±9.7mmHg and root-mean-square-errors (RMSEs) ¿6.4±2.8mmHg against invasive reference cBP waves (n=10). When the aortic geometry was unavailable, the 3-element 0-D algorithm achieved cSBP errors ¿6.0±4.7mmHg and RMSEs ¿5.9±2.4mmHg against non-invasive reference cBP waves (n=171), outperforming the 2-element 0-D algor
Autores: Aramburu Montenegro, Jorge (Autor de correspondencia); Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 1025-5842  Vol. 22  Nº 5  2019  págs. 518 - 532
Balloon-occluded transarterial chemoembolisation (B-TACE) is an intraarterial transcatheter treatment for liver cancer. In B-TACE, an artery-occluding microballoon catheter occludes an artery and promotes collateral circulation for drug delivery to tumours. This paper presents a methodology for analysing the haemodynamics during B-TACE, by combining zero-dimensional and three-dimensional modelling tools. As a proof of concept, we apply the methodology to a patient-specific hepatic artery geometry and analyse two catheter locations. Results show that the blood flow redistribution can be predicted in this proof-of-concept study, suggesting that this approach could potentially be used to optimise catheter location.
Autores: Aramburu Montenegro, Jorge (Autor de correspondencia); Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 2040-7939  Vol. 34  Nº e2983  2018 
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 2040-7939  2017  págs. e02895
Liver radioembolization is a promising treatment option for combating liver tumors. It is performed by placing a microcatheter in the hepatic artery and administering radiation-emitting microspheres through the arterial bloodstream so that they get lodged in the tumoral bed. In avoiding nontarget radiation, the standard practice is to conduct a pretreatment, in which the microcatheter location and injection velocity are decided. However, between pretreatment and actual treatment some of the parameters that influence the particle distribution in the liver can vary, resulting in radiation-induced complications. The present study aims to analyze the influence of a commercially available microcatheter with an angled tip and particle injection velocity in terms of segment-to-segment particle distribution. Specifically, four tip orientations and two injection velocities are combined to yield a set of eight numerical simulations of the particle-hemodynamics in a patient-specific truncated hepatic artery. For each simulation, four cardiac pulses are simulated. Particles are injected during the first cycle, and the remaining pulses enable the majority of the injected particles to exit the computational domain. Results indicate that, in terms of injection velocity, particles are more spread out in the cross-sectional lumen areas as the injection velocity increases. The tip's orientation also plays a role because it influences the near-tip hemodynamics, therefore altering the particle travel through the hepatic artery. However, results suggest that particle distribution tries to match the blood flow split, therefore particle injection velocity and microcatheter tip orientation playing a minor role in segment-to-segment particle distribution.
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 2040-7939  Vol. 33  Nº 2  2017 
Liver radioembolization (RE) is a treatment option for patients with unresectable and chemorefractory primary and metastatic liver tumours. RE consists of intra-arterially administering via catheter radioactive microspheres that locally attack the tumours, sparing healthy tissue. Prior to RE, the standard practice is to conduct a treatment-mimicking pretreatment assessment via the infusion of Tc-99m-labelled macroaggregated albumin microparticles. The usefulness of this pretreatment has been debated in the literature, and thus, the aim of the present study is to shed light on this issue by numerically simulating the liver RE pretreatment and actual treatment particle-haemodynamics in a patient-specific hepatic artery under two different literature-based cancer scenarios and two different placements of a realistic end-hole microcatheter in the proper hepatic artery. The parameters that are analysed are the following: microagent quantity and size (accounting for RE pretreatment and treatment), catheter-tip position (near the proper hepatic artery bifurcation and away from it), and cancer burden (10% and 30% liver involvement). The conclusion that can be reached from the simulations is that when it comes to mimicking RE in terms of delivering particles to tumour-bearing segments, the catheter-tip position is much more important (because of the importance of local haemodynamic pattern alteration) than the infused microagents (i.e. quantity and size). Cancer burden is another important feature because the increase in blood flow rate to tumour-bearing segments increases the power to drag particles. These numerical simulation-based conclusions are in agreement with clinically observed events reported in the literature. Copyright (c) 2016 John Wiley & Sons, Ltd.
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Borro Yagüez, Diego; et al.
ISSN 2057-1976  Vol. 2  Nº 1  2016  págs. 015001
The analysis of the progression of cardiovascular diseases is an active area of ongoing research. This paper develops an image registration-based methodology to quantify the patient-specific local blood vessel shape variations that occur in the radial direction (i.e. expansion or shrinkage) over an imaging follow-up period, and an example is presented as proof of principle. The methodology can be used for complex vessels with bifurcations, and it is able to identify and address vessel deformations if changes in tortuosity or longitudinal direction are small. The methodology consists of (a) overlapping the baseline and follow-up vessel surfaces by matching the lumen centerline, (b) dividing the region of interest into slices perpendicular to the centerline and centering each slice, and (c) dividing each centered slice into sectors. The local approach consists of analyzing a representative point in each sector of each slice (i.e. each patch). In this paper the algorithm is applied to a patient-specific abdominal aortic aneurysm (AAA) as a proof of principle of the method. Six patient-specific image reconstructions from a single subject followed for 28 months are analyzed in pairs, yielding five time spans to which the algorithm was applied. The algorithm was able to quantify the AAA radial growth. The average AAA radial growths for the five case studies are ¿2.13 mm, 3.43 mm, ¿0.25 mm, 1.41 mm, and 0.84 mm, whereas the maximum local growths are 4.76 ± 0.15 mm, 9.30 ± 1.13 mm, 2.08 ± 0.05 mm, 4.10 ± 0.14 mm, and 4.16 ± 0.45 mm. The tolerance of the geometric local measurements is related to the matching processes (i.e. overlapping the geometries and centering each slice) because of the vessel deformation that took place over time. Thus, this methodology has been used to quantify the average AAA growth and the maximum local AAA growth (± the tolerance) as metrics of the vessel's radial growth.
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 0021-9290  Vol. 49   Nº 15   2016  págs. 3705 - 3713
Radioembolization, which consist of the implantation of radioactive microspheres via intra-arterially placed microcatheter, is a safe and effective treatment for liver cancer. Nevertheless, radioembolization-related complications and side effects may arise, which are an active area of ongoing research. The catheter design has been claimed as an option in reducing these complications. In this paper, the influence of catheter type and location are investigated. The study was undertaken by numerically simulating the particle¿hemodynamics in a patient-specific hepatic artery during liver radioembolization. The parameters modified were cancer scenario (30% liver involvement in the right lobe, `scenario A¿, and in both lobes, `scenario B¿), catheter type (standard end-hole microcatheter, SMC, and antireflux catheter, ARC), and the location of the tip in the proper hepatic artery (in the straight part, `inlet¿, and near the bifurcation, `bifurcation¿). Comparing ARC with SMC, the maximum and average (over segments) absolute difference in the percentage of particles that reached each segment were 19.62% and 9.06% when injecting near the inlet for scenario A; 3.54% and 1.07% injecting near the bifurcation for scenario A; and 18.31% and 11.85% injecting near the inlet for scenario B. It seems, therefore, that the location of the catheter tip in the artery is crucial in terms of particle distribution. Moreover, even though the near-tip blood flow was altered due to the presence of a catheter, the particle distribution matched the flow split if the distance between the injection point and the first bifurcation encountered enabled the alignment of particles with blood flow.
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 2040-7939  Vol. 32  Nº 11  2016  págs. e02764
Some of the latest treatments for unresectable liver malignancies (primary or metastatic tumours), which include bland embolisation, chemoembolisation, and radioembolisation, among others, take advantage of the increased arterial blood supply to the tumours to locally attack them. A better understanding of the factors that influence this transport may help improve the therapeutic procedures by taking advantage of flow patterns or by designing catheters and infusion systems that result in the injected beads having increased access to the tumour vasculature. Computational analyses may help understand the haemodynamic patterns and embolic-microsphere transport through the hepatic arteries. In addition, physiological inflow and outflow boundary conditions are essential in order to reliably represent the blood flow through arteries. This study presents a liver cancer arterial perfusion model based on a literature review and derives boundary conditions for tumour-bearing liver-feeding hepatic arteries based on the arterial perfusion characteristics of normal and tumorous liver segment tissue masses and the hepatic artery branching configuration. Literature-based healthy and tumour-bearing realistic scenarios are created and haemodynamically analysed for the same patient-specific hepatic artery. As a result, this study provides boundary conditions for computational fluid dynamics simulations that will allow researchers to numerically study, for example, various intravascular devices
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Rivas Nieto, Alejandro; et al.
ISSN 0021-9290  Vol. 49  Nº 15  2016  págs. 3714 - 3721
Liver radioembolization is a treatment option for patients with primary and secondary liver cancer. The procedure consists of injecting radiation-emitting microspheres via an intra-arterially placed microcatheter, enabling the deposition of the microspheres in the tumoral bed. The microcatheter location and the particle injection rate are determined during a pretreatment work-up. The purpose of this study was to numerically study the effects of the injection characteristics during the first stage of microsphere travel through the bloodstream in a patient-specific hepatic artery (i.e., the near-tip particle¿hemodynamics and the segment-to-segment particle distribution). Specifically, the influence of the distal direction of an end-hole microcatheter and particle injection point and velocity were analyzed. Results showed that the procedure targeted the right lobe when injecting from two of the three injection points under study and the remaining injection point primarily targeted the left lobe. Changes in microcatheter direction and injection velocity resulted in an absolute difference in exiting particle percentage for a given liver segment of up to 20% and 30%, respectively. It can be concluded that even though microcatheter placement is presumably reproduced in the treatment session relative to the pretreatment angiography, the treatment may result in undesired segment-to-segment particle distribution and therefore undesired treatment outcomes due to modifications of any of the parameters studied, i.e., microcatheter direction and particle injection point and velocity.
Autores: Aramburu Montenegro, Jorge; Antón Remírez, Raúl; Bernal, N.; et al.
ISSN 0954-4119  Vol. 229  Nº 4  2015  págs. 291 - 306
Physiological outflow boundary conditions are necessary to carry out computational fluid dynamics simulations that reliably represent the blood flow through arteries. When dealing with complex three-dimensional trees of small arteries, and therefore with multiple outlets, the robustness and speed of convergence are also important. This study derives physiological outflow boundary conditions for cases in which the physiological values at those outlets are not known (neither in vivo measurements nor literature-based values are available) and in which the tree exhibits symmetry to some extent. The inputs of the methodology are the three-dimensional domain and the flow rate waveform and the systolic and diastolic pressures at the inlet. The derived physiological outflow boundary conditions, which are a physiological pressure waveform for each outlet, are based on the results of a zero-dimensional model simulation. The methodology assumes symmetrical branching and is able to tackle the flow distribution problem when the domain outlets are at branches with a different number of upstream bifurcations. The methodology is applied to a group of patient-specific arteries in the liver. The methodology is considered to be valid because the pulsatile computational fluid dynamics simulation with the inflow flow rate waveform (input of the methodology) and the derived outflow boundary conditions lead to physiological results, that is, the resulting systolic and diastolic pressures at the inlet match the inputs of the methodology, and the flow split is also physiological.