Search Results
Multiscale Coupling of One-dimensional Vascular Models and Elastic Tissues
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).
Model hierarchies and higher-order discretisation of time-dependent thin-film free boundary problems with dynamic contact angle
We present a mathematical and numerical framework for the physical problem of thin-film fluid flows over planar surfaces including dynamic contact angles. In particular, we provide algorithmic details and an implementation of higher-order spatial and temporal discretisation of the underlying free boundary problem using the finite element method. The corresponding partial differential equation is based on a thermodynamic consistent energetic variational formulation of the problem using the free energy and viscous dissipation in the bulk, on the surface, and at the moving contact line. Model hierarchies for limits of strong and weak contact line dissipation are established, implemented and studied. We analyze the performance of the numerical algorithm and investigate the impact of the dynamic contact angle on the evolution of two benchmark problems: gravity-driven sliding droplets and the instability of a ridge.
Multiscale coupling of one-dimensional vascular models and elastic tissues
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. The pressure and velocity of the blood in the vessels 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.
Multiscale modeling of vascularized tissues via non-matching immersed methods
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.
Variational approach to fluid-structure interaction via GENERIC
We present a framework to systematically derive variational formulations for fluid-structure interaction problems based on thermodynamical driving functionals and geometric structures in different coordinate systems by suitable transformations within this formulation. Our approach provides a promising basis to construct structure-preserving discretization strategies.
Error estimates in weighted Sobolev norms for finite element immersed interface methods
When solving elliptic partial differential equations in a region containing immersed interfaces (possibly evolving in time), it is often desirable to approximate the problem using a uniform background discretisation, not aligned with the interface itself. Optimal convergence rates are possible if the discretisation scheme is enriched by allowing the discrete solution to have jumps aligned with the surface, at the cost of a higher complexity in the implementation. A much simpler way to reformulate immersed interface problems consists in replacing the interface by a singular force field that produces the desired interface conditions, as done in immersed boundary methods. These methods are known to have inferior convergence properties, depending on the global regularity of the solution across the interface, when compared to enriched methods. In this work we prove that this detrimental effect on the convergence properties of the approximate solution is only a local phenomenon, restricted to a small neighbourhood of the interface. In particular we show that optimal approximations can be constructed in a natural and inexpensive way, simply by reformulating the problem in a distributionally consistent way, and by resorting to weighted norms when computing the global error of the approximation.
Error estimates of B-spline based finite-element method for the wind-driven ocean circulation
We present the error analysis of a B-spline based finite-element approximation of the stream-function formulation of the large scale winddriven ocean circulation. In particular, we derive optimal error estimates for h-refinement using a Nitsche-type variational formulations of the two simplied linear models of the stationary quasigeostrophic equations, namely the Stommel and StommelMunk models. Numerical results on rectangular and embedded geometries confirm the error analysis.