In this paper we consider the numerical solution of the three-dimensional fluid-structure interaction problem in haemodynamics, in the case of real geometries, physiological data and finite elasticity vessel deformati...
详细信息
In this paper we consider the numerical solution of the three-dimensional fluid-structure interaction problem in haemodynamics, in the case of real geometries, physiological data and finite elasticity vessel deformations. We study some new inexact schemes, obtained from semi-implicit approximations, which treat exactly the physical interface conditions while performing just one or few iterations for the management of the interface position and of the fluid and structure non-linearities. We show that such schemes allow to improve the efficiency while preserving the accuracy of the related exact (implicit) scheme. To do this we consider both a simple analytical test case and two real cases of clinical interest in haemodynamics. We also provide an error analysis for a simple differential model problem when a BDF method is considered for the time discretization and only few Newton iterations are performed at each temporal instant. (C) 2014 Elsevier Inc. All rights reserved.
Osteocytes are believed to be the primary sensor of mechanical stimuli in bone, which orchestrate osteoblasts and osteoclasts to adapt bone structure and composition to meet physiological loading demands. Experimental...
详细信息
Osteocytes are believed to be the primary sensor of mechanical stimuli in bone, which orchestrate osteoblasts and osteoclasts to adapt bone structure and composition to meet physiological loading demands. Experimental studies to quantify the mechanical environment surrounding bone cells are challenging, and as such, computational and theoretical approaches have modelled either the solid or fluid environment of osteocytes to predict how these cells are stimulated in vivo. Osteocytes are an elastic cellular structure that deforms in response to the external fluid flow imposed by mechanical loading. This represents a most challenging multi-physics problem in which fluid and solid domains interact, and as such, no previous study has accounted for this complex behaviour. The objective of this study is to employ fluid-structure interaction (FSI) modelling to investigate the complex mechanical environment of osteocytes in vivo. Fluorescent staining of osteocytes was performed in order to visualise their native environment and develop geometrically accurate models of the osteocyte in vivo. By simulating loading levels representative of vigorous physiological activity ( compression and 300 Pa pressure gradient), we predict average interstitial fluid velocities and average maximum shear stresses surrounding osteocytes in vivo. Interestingly, these values occur in the canaliculi around the osteocyte cell processes and are within the range of stimuli known to stimulate osteogenic responses by osteoblastic cells in vitro. Significantly our results suggest that the greatest mechanical stimulation of the osteocyte occurs in the cell processes, which, cell culture studies have indicated, is the most mechanosensitive area of the cell. These are the first computational FSI models to simulate the complex multi-physics mechanical environment of osteocyte in vivo and provide a deeper understanding of bone mechanobiology.
This article presents an interface-energy-conserving coupling strategy for transient fluid-structure interaction. The solid sub-domain is discretized by finite element method with Newmark time integrator, whereas for ...
详细信息
This article presents an interface-energy-conserving coupling strategy for transient fluid-structure interaction. The solid sub-domain is discretized by finite element method with Newmark time integrator, whereas for the fluid sub-domain we use the mesh-less method SPH-ALE with 2nd order Runge-Kutta scheme. This paper proposes a method to impose a mean interface normal velocity continuity in such a way that the algorithmic interface energy is zero during the whole period of numerical simulation. This coupling method thus ensures that the coupled problem shall be stable in time. Secondly it will converge in time with the rate of convergence of the worst time integrator chosen for each problem. The proposed method is first applied to a mono-dimensional problem by which we investigate the phenomena of propagation of shock waves across the fluid-structure interface. A good agreement is observed between the numerical result and the analytical solution in the 1D shock wave propagation test cases. Finally, a multidimensional example is presented and compared to another coupling approach. (C) 2014 Elsevier Ltd. All rights reserved.
We consider the extension of the Nitsche method to the case of fluidstructureinteraction problems on unfitted meshes. We give a stability analysis for the space semi-discretized problem and show how this estimate ma...
详细信息
We consider the extension of the Nitsche method to the case of fluidstructureinteraction problems on unfitted meshes. We give a stability analysis for the space semi-discretized problem and show how this estimate may be used to derive optimal error estimates for smooth solutions, irrespectively of the mesh/interface intersection. We also discuss different strategies for the time discretization, using either fully implicit or explicit coupling (loosely coupled) schemes. Some numerical examples illustrate the theoretical discussion. (C) 2014 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://***/licenses/by-nc-nd/3.0/).
A computational human mitral valve (MV) model under physiological pressure loading is developed using a hybrid finite element immersed boundary method, which incorporates experimentally-based constitutive laws in a th...
详细信息
A computational human mitral valve (MV) model under physiological pressure loading is developed using a hybrid finite element immersed boundary method, which incorporates experimentally-based constitutive laws in a three-dimensional fluid-structure interaction framework. A transversely isotropic material constitutive model is used to characterize the mechanical behaviour of the MV tissue based on recent mechanical tests of healthy human mitral leaflets. Our results show good agreement, in terms of the flow rate and the closing and opening configurations, with measurements from in vivo magnetic resonance images. The stresses in the anterior leaflet are found to be higher than those in the posterior leaflet and are concentrated around the annulus trigons and the belly of the leaflet. The results also show that the chordae play an important role in providing a secondary orifice for the flow when the valve opens. Although there are some discrepancies to be overcome in future work, our simulations show that the developed computational model is promising in mimicking the in vivo MV dynamics and providing important information that are not obtainable by in vivo measurements. (c) 2014 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd.
For the large-scale offshore wind turbine blades, the governing equations in fluid domain and the motion equations in structural domain with geometric nonlinearity were developed based on ALE description, and the corr...
详细信息
For the large-scale offshore wind turbine blades, the governing equations in fluid domain and the motion equations in structural domain with geometric nonlinearity were developed based on ALE description, and the corresponding discrete equations were obtained. Blade entity model was built up using Pro/E, and the blade vibration characteristics under fluid-structure interaction (FSI) were simulated using ANSYS. Numerical results show that wind shear effect greatly increases the peak values of response curves for displacement and stress, makes the effect of bi-directional fluid-structure interaction (BFSI) more obvious, and also accelerates the attenuation of vibration curves. The displacement of blade airfoil increases nonlinearly along the span direction, and reaches the maximum at the blade tip. The maximum Mises stress appears in the middle of the blade, and reduces gradually towards each end of the blade. Furthermore, the contribution of wind shear effect (WSE) to displacement and Mises stress is much greater than that of FSI. (C) 2014 Elsevier Ltd. All rights reserved.
The aim of this paper is to analyze an elastoacoustic vibration problem employing a dual-mixed formulation in the solid domain. The Cauchy stress tensor and the rotation are the primary variables in the elastic struct...
详细信息
The aim of this paper is to analyze an elastoacoustic vibration problem employing a dual-mixed formulation in the solid domain. The Cauchy stress tensor and the rotation are the primary variables in the elastic structure while the standard pressure formulation is considered in the acoustic fluid. The resulting mixed eigenvalue problem is approximated by a conforming Galerkin scheme based on the lowest order Lagrange and Arnold-Falk-Winther finite element subspaces in the fluid and solid domains, respectively. We show that the scheme provides a correct approximation of the spectrum and prove quasi-optimal error estimates. Finally, we report some numerical experiments. (C) 2014 Elsevier Ltd. All rights reserved.
In this paper, we present a stable second-order time accurate scheme for solving fluid-structure interaction problems. The scheme uses so-called Combined Field with Explicit Interface (CFEI) advancing formulation base...
详细信息
In this paper, we present a stable second-order time accurate scheme for solving fluid-structure interaction problems. The scheme uses so-called Combined Field with Explicit Interface (CFEI) advancing formulation based on the Arbitrary Lagrangian-Eulerian approach with finite element procedure. Although loosely-coupled partitioned schemes are often popular choices for simulating FSI problems, these schemes may suffer from inherent instability at low structure to fluid density ratios. We show that our second-order scheme is stable for any mass density ratio and hence is able to handle strong added-mass effects. Energy-based stability proof relies heavily on the connections among extrapolation formula, trapezoidal scheme for second-order equation, and backward difference method for first-order equation. Numerical accuracy and stability of the scheme is assessed with the aid of two-dimensional fluid-structure interaction problems of increasing complexity. We confirm second-order temporal accuracy by numerical experiments on an elastic semi-circular cylinder problem. We verify the accuracy of coupled solutions with respect to the benchmark solutions of a cylinder-elastic bar and the Navier-Stokes flow system. To study the stability of the proposed scheme for strong added-mass effects, we present new results using the combined field formulation for flexible flapping motion of a thin-membrane structure with low mass ratio and strong added-mass effects in a uniform axial flow. Using a systematic series of fluid-structure simulations, a detailed analysis of the coupled response as a function of mass ratio for the case of very low bending rigidity has been presented. (C) 2014 Elsevier Inc. All rights reserved.
Background: Atherosclerotic plaque progression and rupture are believed to be associated with mechanical stress conditions. In this paper, patient-specific in vivo intravascular ultrasound (IVUS) coronary plaque image...
详细信息
Background: Atherosclerotic plaque progression and rupture are believed to be associated with mechanical stress conditions. In this paper, patient-specific in vivo intravascular ultrasound (IVUS) coronary plaque image data were used to construct computational models with fluid-structure interaction (FSI) and cyclic bending to investigate correlations between plaque wall thickness and both flow shear stress and plaque wall stress conditions. Methods: IVUS data were acquired from 10 patients after voluntary informed consent. The X-ray angiogram was obtained prior to the pullback of the IVUS catheter to determine the location of the coronary artery stenosis, vessel curvature and cardiac motion. Cyclic bending was specified in the model representing the effect by heart contraction. 3D anisotropic FSI models were constructed and solved to obtain flow shear stress (FSS) and plaque wall stress (PWS) values. FSS and PWS values were obtained for statistical analysis. Correlations with p < 0.05 were deemed significant. Results: Nine out of the 10 patients showed positive correlation between wall thickness and flow shear stress. The mean Pearson correlation r-value was 0.278 +/- 0.181. Similarly, 9 out of the 10 patients showed negative correlation between wall thickness and plaque wall stress. The mean Pearson correlation r-value was -0.530 +/- 0.210. Conclusion: Our results showed that plaque vessel wall thickness correlated positively with FSS and negatively with PWS. The patient-specific IVUS-based modeling approach has the potential to be used to investigate and identify possible mechanisms governing plaque progression and rupture and assist in diagnosis and intervention procedures. This represents a new direction of research. Further investigations using more patient follow-up data are warranted.
"Experimental and numerical studies of fluid–structureinteraction phenomena inside the head when subjected to a dynamical loading."Computer Methods in Biomechanics and Biomedical Engineering, 17(sup1), pp....
详细信息
"Experimental and numerical studies of fluid–structureinteraction phenomena inside the head when subjected to a dynamical loading."Computer Methods in Biomechanics and Biomedical Engineering, 17(sup1), pp. 46–47
暂无评论