No Arabic abstract
Recent advances in imaging techniques have enabled us to visualize lung tumors or nodules in early-stage cancer. However, the positions of nodules can change because of intraoperative lung deflation, and the modeling of pneumothorax-associated deformation remains a challenging issue for intraoperative tumor localization. In this study, we introduce spatial and geometric analysis methods for inflated/deflated lungs and discuss heterogeneity in pneumothorax-associated deformation. Contrast-enhanced CT images simulating intraoperative conditions were acquired from live Beagle dogs. Deformable mesh registration techniques were designed to map the surface and subsurface tissues of lung lobes. The developed framework addressed local mismatches of bronchial tree structures and achieved stable registration with a Hausdorff distance of less than 1 mm and a target registration error of less than 5 mm. Our results show that the strain of lung parenchyma was 35% higher than that of bronchi, and that subsurface deformation in the deflated lung is heterogeneous.
We propose a nonlinear registration-based model reduction procedure for rapid and reliable solution of parameterized two-dimensional steady conservation laws. This class of problems is challenging for model reduction techniques due to the presence of nonlinear terms in the equations and also due to the presence of parameter-dependent discontinuities that cannot be adequately represented through linear approximation spaces. Our approach builds on a general (i.e., independent of the underlying equation) registration procedure for the computation of a mapping $Phi$ that tracks moving features of the solution field and on an hyper-reduced least-squares Petrov-Galerkin reduced-order model for the rapid and reliable computation of the solution coefficients. The contributions of this work are twofold. First, we investigate the application of registration-based methods to two-dimensional hyperbolic systems. Second, we propose a multi-fidelity approach to reduce the offline costs associated with the construction of the parameterized mapping and the reduced-order model. We discuss the application to an inviscid supersonic flow past a parameterized bump, to illustrate the many features of our method and to demonstrate its effectiveness.
We present a general -- i.e., independent of the underlying equation -- registration procedure for parameterized model order reduction. Given the spatial domain $Omega subset mathbb{R}^2$ and the manifold $mathcal{M}= { u_{mu} : mu in mathcal{P} }$ associated with the parameter domain $mathcal{P} subset mathbb{R}^P$ and the parametric field $mu mapsto u_{mu} in L^2(Omega)$, our approach takes as input a set of snapshots ${ u^k }_{k=1}^{n_{rm train}} subset mathcal{M}$ and returns a parameter-dependent bijective mapping ${Phi}: Omega times mathcal{P} to mathbb{R}^2$: the mapping is designed to make the mapped manifold ${ u_{mu} circ {Phi}_{mu}: , mu in mathcal{P} }$ more amenable for linear compression methods. In this work, we extend and further analyze the registration approach proposed in [Taddei, SISC, 2020]. The contributions of the present work are twofold. First, we extend the approach to deal with annular domains by introducing a suitable transformation of the coordinate system. Second, we discuss the extension to general two-dimensional geometries: towards this end, we introduce a spectral element approximation, which relies on a partition ${ Omega_{q} }_{q=1} ^{N_{rm dd}}$ of the domain $Omega$ such that $Omega_1,ldots,Omega_{N_{rm dd}}$ are isomorphic to the unit square. We further show that our spectral element approximation can cope with parameterized geometries. We present rigorous mathematical analysis to justify our proposal; furthermore, we present numerical results for a heat-transfer problem in an annular domain, a potential flow past a rotating symmetric airfoil, and an inviscid transonic compressible flow past a non-symmetric airfoil, to demonstrate the effectiveness of our method.
We propose a model reduction procedure for rapid and reliable solution of parameterized hyperbolic partial differential equations. Due to the presence of parameter-dependent shock waves and contact discontinuities, these problems are extremely challenging for traditional model reduction approaches based on linear approximation spaces. The main ingredients of the proposed approach are (i) an adaptive space-time registration-based data compression procedure to align local features in a fixed reference domain, (ii) a space-time Petrov-Galerkin (minimum residual) formulation for the computation of the mapped solution, and (iii) a hyper-reduction procedure to speed up online computations. We present numerical results for a Burgers model problem and a shallow water model problem, to empirically demonstrate the potential of the method.
We propose and analyse the properties of a new class of models for the electromechanics of cardiac tissue. The set of governing equations consists of nonlinear elasticity using a viscoelastic and orthotropic exponential constitutive law (this is so for both active stress and active strain formulations of active mechanics) coupled with a four-variable phenomenological model for human cardiac cell electrophysiology, which produces an accurate description of the action potential. The conductivities in the model of electric propagation are modified according to stress, inducing an additional degree of nonlinearity and anisotropy in the coupling mechanisms; and the activation model assumes a simplified stretch-calcium interaction generating active tension or active strain. The influence of the new terms in the electromechanical model is evaluated through a sensitivity analysis, and we provide numerical validation through a set of computational tests using a novel mixed-primal finite element scheme.
We perform the linear stability analysis for a new model for poromechanical processes with inertia (formulated in mixed form using the solid deformation, fluid pressure, and total pressure) interacting with diffusing and reacting solutes convected in the medium. We find parameter regions that lead to spatio-temporal instabilities of the coupled system. The mutual dependences between deformation and diffusive patterns are of substantial relevance in the study of morphoelastic changes in biomaterials. We provide a set of computational examples in 2D and 3D (related to brain mechanobiology) that can be used to form a better understanding on how, and up to which extent, the deformations of the porous structure dictate the generation and suppression of spatial patterning dynamics, also related to the onset of mechano-chemical waves.