Filters

2011 - 2019202020 - 20227

More filters

2019 - 2019162020 - 202211
Reset filters

Settings

Author: Caiazzo, Alfonso ×

Search Results

Now showing 1 - 10 of 27
  • Thumbnail Image
    Item
    Inversion-recovery MR elastography of the human brain for improved stiffness quantification near fluid-solid boundaries
    (New York, NY [u.a.] : Wiley-Liss, 2021) Lilaj, Ledia; Herthum, Helge; Meyer, Tom; Shahryari, Mehrgan; Bertalan, Gergely; Caiazzo, Alfonso; Braun, Jürgen; Fischer, Thomas; Hirsch, Sebastian; Sack, Ingolf
    Purpose: In vivo MR elastography (MRE) holds promise as a neuroimaging marker. In cerebral MRE, shear waves are introduced into the brain, which also stimulate vibrations in adjacent CSF, resulting in blurring and biased stiffness values near brain surfaces. We here propose inversion-recovery MRE (IR-MRE) to suppress CSF signal and improve stiffness quantification in brain surface areas. Methods: Inversion-recovery MRE was demonstrated in agar-based phantoms with solid-fluid interfaces and 11 healthy volunteers using 31.25-Hz harmonic vibrations. It was performed by standard single-shot, spin-echo EPI MRE following 2800-ms IR preparation. Wave fields were acquired in 10 axial slices and analyzed for shear wave speed (SWS) as a surrogate marker of tissue stiffness by wavenumber-based multicomponent inversion. Results: Phantom SWS values near fluid interfaces were 7.5 ± 3.0% higher in IR-MRE than MRE (P =.01). In the brain, IR-MRE SNR was 17% lower than in MRE, without influencing parenchymal SWS (MRE: 1.38 ± 0.02 m/s; IR-MRE: 1.39 ± 0.03 m/s; P =.18). The IR-MRE tissue–CSF interfaces appeared sharper, showing 10% higher SWS near brain surfaces (MRE: 1.01 ± 0.03 m/s; IR-MRE: 1.11 ± 0.01 m/s; P <.001) and 39% smaller ventricle sizes than MRE (P <.001). Conclusions: Our results show that brain MRE is affected by fluid oscillations that can be suppressed by IR-MRE, which improves the depiction of anatomy in stiffness maps and the quantification of stiffness values in brain surface areas. Moreover, we measured similar stiffness values in brain parenchyma with and without fluid suppression, which indicates that shear wavelengths in solid and fluid compartments are identical, consistent with the theory of biphasic poroelastic media. © 2021 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine
  • Thumbnail Image
    Item
    Multiscale Coupling of One-dimensional Vascular Models and Elastic Tissues
    (Dordrecht [u.a.] : Springer Science + Business Media B.V, 2021) Heltai, Luca; Caiazzo, Alfonso; Müller, Lucas O.
    We present a computational multiscale model for the efficient simulation of vascularized tissues, composed of an elastic three-dimensional matrix and a vascular network. The effect of blood vessel pressure on the elastic tissue is surrogated via hyper-singular forcing terms in the elasticity equations, which depend on the fluid pressure. In turn, the blood flow in vessels is treated as a one-dimensional network. Intravascular pressure and velocity are simulated using a high-order finite volume scheme, while the elasticity equations for the tissue are solved using a finite element method. This work addresses the feasibility and the potential of the proposed coupled multiscale model. In particular, we assess whether the multiscale model is able to reproduce the tissue response at the effective scale (of the order of millimeters) while modeling the vasculature at the microscale. We validate the multiscale method against a full scale (three-dimensional) model, where the fluid/tissue interface is fully discretized and treated as a Neumann boundary for the elasticity equation. Next, we present simulation results obtained with the proposed approach in a realistic scenario, demonstrating that the method can robustly and efficiently handle the one-way coupling between complex fluid microstructures and the elastic matrix. © 2021, The Author(s).
  • Thumbnail Image
    Item
    Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers
    (New York, NY [u.a.] : Springer Science + Business Media B.V., 2020) Blank, Laura; Rioseco, Ernesto Meneses; Caiazzo, Alfonso; Wilbrandt, Ulrich
    Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and numerically investigate the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and a variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geological configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems. © 2020, The Author(s).
  • Thumbnail Image
    Item
    On the feasibility of using open source solvers for the simulation of a turbulent air flow in a dairy barn
    (Amsterdam [u.a.] : Elsevier, 2020) Janke, David; Caiazzo, Alfonso; Ahmed, Naveed; Alia, Najib; Knoth, Oswald; Moreau, Baptiste; Wilbrandt, Ulrich; Willink, Dilya; Amon, Thomas; John, Volker
    Two transient open source solvers, OpenFOAM and ParMooN, and the commercial solver Ansys Fluent are assessed with respect to the simulation of the turbulent air flow inside and around a dairy barn. For this purpose, data were obtained in an experimental campaign at a 1:100 scaled wind tunnel model. All solvers used different meshes, discretization schemes, and turbulence models. The experimental data and numerical results agree well for time-averaged stream-wise and vertical-wise velocities. In particular, the air exchange was predicted with high accuracy by both open source solvers with relative differences less than 4% and by the commercial solver with a relative difference of 9% compared to the experimental results. With respect to the turbulent quantities, good agreements at the second (downwind) half of the barn inside and especially outside the barn could be achieved, where all codes accurately predicted the flow separation and, in many cases, the root-mean-square velocities. Deviations between simulations and experimental results regarding turbulent quantities could be observed in the first part of the barn. These deviations can be attributed to the utilization of roughness elements between inlet and barn in the experiment that were not modeled in the numerical simulations. Both open source solvers proved to be promising tools for the accurate prediction of time-dependent phenomena in an agricultural context, e.g., like the transport of particulate matter or pathogen-laden aerosols in and around agricultural buildings. © 2020 The Authors
  • Thumbnail Image
    Item
    Multiscale modeling of vascularized tissues via non-matching immersed methods
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Heltai, Luca; Caiazzo, Alfonso
    We consider a multiscale approach based on immersed methods for the efficient computational modeling of tissues composed of an elastic matrix (in two or three-dimensions) and a thin vascular structure (treated as a co-dimension two manifold) at a given pressure. We derive different variational formulations of the coupled problem, in which the effect of the vasculature can be surrogated in the elasticity equations via singular or hyper-singular forcing terms. These terms only depends on information defined on co-dimension two manifolds (such as vessel center line, cross sectional area, and mean pressure over cross section), thus drastically reducing the complexity of the computational model. We perform several numerical tests, ranging from simple cases with known exact solutions to the modeling of materials with random distributions of vessels. In the latter case, we use our immersed method to perform an in silico characterization of the mechanical properties of the effective biphasic material tissue via statistical simulations.
  • Thumbnail Image
    Item
    ParMooN - a modernized program package based on mapped finite elements
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Wilbrandt, Ulrich; Bartsch, Clemens; Ahmed, Naveed; Alia, Najib; Anker, Felix; Blank, Laura; Caiazzo, Alfonso; Ganesa, Sashikumaar; Giere, Swetlana; Matthies, Gunar; Meesala, Raviteja; Shamim, Abdus; Venkatesan, Jagannath; John, Volker
    PARMOON is a program package for the numerical solution of elliptic and parabolic partial differential equations. It inherits the distinct features of its predecessor MOONMD [28]: strict decoupling of geometry and finite element spaces, implementation of mapped finite elements as their definition can be found in textbooks, and a geometric multigrid preconditioner with the option to use different finite element spaces on different levels of the multigrid hierarchy. After having presented some thoughts about in-house research codes, this paper focuses on aspects of the parallelization, which is the main novelty of PARMOON. Numerical studies, performed on compute servers, assess the efficiency of the parallelized geometric multigrid preconditioner in comparison with parallel solvers that are available in external libraries. The results of these studies give a first indication whether the cumbersome implementation of the parallelized geometric multigrid method was worthwhile or not.
  • Thumbnail Image
    Item
    Reconstruction of flow domain boundaries from velocity data via multi-step optimization of distributed resistance
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2022) Pártl, Ondřej; Wilbrandt, Ulrich; Mura, Joaquín; Caiazzo, Alfonso
    We reconstruct the unknown shape of a flow domain using partially available internal velocity measurements. This inverse problem is motivated by applications in cardiovascular imaging where motion-sensitive protocols, such as phase-contrast MRI, can be used to recover three-dimensional velocity fields inside blood vessels. In this context, the information about the domain shape serves to quantify the severity of pathological conditions, such as vessel obstructions. We consider a flow modeled by a linear Brinkman problem with a fictitious resistance accounting for the presence of additional boundaries. To reconstruct these boundaries, we employ a multi-step gradient-based variational method to compute a resistance that minimizes the difference between the computed flow velocity and the available data. Afterward, we apply different post-processing steps to reconstruct the shape of the internal boundaries. To limit the overall computational cost, we use a stabilized equal-order finite element method. We prove the stability and the well-posedness of the considered optimization problem. We validate our method on three-dimensional examples based on synthetic velocity data and using realistic geometries obtained from cardiovascular imaging.
  • Thumbnail Image
    Item
    Analysis of a stabilized penalty-free Nitsche method for the Brinkman, Stokes, and Darcy problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Blank, Laura; Caiazzo, Alfonso; Chouly, Franz; Lozinski, Alexei; Mura, Joaquin
    In this paper we study the Brinkman model as a unified framework to allow the transition between the Darcy and the Stokes problems. We propose an unconditionally stable low-order finite element approach, which is robust with respect to the whole range of physical parameters, and is based on the combination of stabilized equal-order finite elements with a non-symmetric penalty free Nitsche method for the weak imposition of essential boundary conditions. In particular, we study the properties of the penalty-free Nitsche formulation for the Brinkman setting, extending a recently reported analysis for the case of incompressible elasticity (T. Boiveau & E. Burman. IMA J. Numer. Anal. 36 (2016), no.2, 770-795). Focusing on the two-dimensional case, we obtain optimal a priori error estimates in a mesh-dependent norm, which, converging to natural norms in the cases of Stokes or Darcy flows, allows to extend the results also to these limits. Moreover, we show that, in order to obtain robust estimates also in the Darcy limit, the formulation shall be equipped with a Grad-Div stabilization and an additional stabilization to control the discontinuities of the normal velocity along the boundary. The conclusions of the analysis are supported by numerical simulations.
  • Thumbnail Image
    Item
    Computational modelling and simulation of cancer growth and migration within a 3D heterogeneous tissue: The effects of fibre and vascular structure
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Macnamara, Cicely K.; Caiazzo, Alfonso; Ramis-Conde, Ignacio; Chaplain, Mark A.J.
    The term cancer covers a multitude of bodily diseases, broadly categorised by having cells which do not behave normally. Since cancer cells can arise from any type of cell in the body, cancers can grow in or around any tissue or organ making the disease highly complex. Our research is focused on understanding the specific mechanisms that occur in the tumour microenvironment via mathematical and computational modeling. We present a 3D individual-based model which allows one to simulate the behaviour of, and spatio-temporal interactions between, cells, extracellular matrix fibres and blood vessels. Each agent (a single cell, for example) is fully realised within the model and interactions are primarily governed by mechanical forces between elements. However, as well as the mechanical interactions we also consider chemical interactions, for example, by coupling the code to a finite element solver to model the diffusion of oxygen from blood vessels to cells. The current state of the art of the model allows us to simulate tumour growth around an arbitrary blood-vessel network or along the striations of fibrous tissue.
  • Thumbnail Image
    Item
    Computational haemodynamics in stenotic internal jugular veins
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Caiazzo, Alfonso; Montecinos, Gino; Müller, Lucas O.; Toro, Eleuterio F.; Haacke, E. Mark
    Stenosis in internal jugular veins (IJVs) are frequently associated to pathological venous circulation and insufficient cerebral blood drainage. In this work, we set up a computational framework to assess the relevance of IJV stenoses through numerical simulation, combining medical imaging, patient-specific data and a mathematical model for venous occlusions. Coupling a three-dimensional (3D) description of blood flow in IJVs with a reduced one-dimesional model (1D) for major intracranial veins, we are able to model different anatomical configurations, an aspect of importance to understand the impact of IJV stenosis in intracranial venous haemodynamics. We investigate several stenotic configurations in a physiologic patient-specific regime, quantifying the effect of the stenosis in terms of venous pressure increase and wall shear stress patterns. Simulation results are in qualitative agreement with reported pressure anomalies in pathological cases. Moreover, they demonstrate the potential of the proposed multiscale framework for individualbased studies and computer-aided diagnosis.