We study a recent formulation for fluid-structure interaction problems based on the use of a distributed Lagrange multiplier in the spirit of the fictitious domain approach. In this paper, we focus our attention on a ...
详细信息
We study a recent formulation for fluid-structure interaction problems based on the use of a distributed Lagrange multiplier in the spirit of the fictitious domain approach. In this paper, we focus our attention on a crucial computational aspect regarding the interface matrix for the finite element discretization: it involves integration of functions supported on two different meshes. Several numerical tests show that accurate computation of the interface matrix has to be performed in order to ensure the optimal convergence of the method. (c) 2022 Elsevier B.V. All rights reserved.
In this article, we formulate a monolithic optimal control method for general timedependent fluid-structure interaction (FSI) systems with large solid deformation: we consider a displacement-tracking type of objective...
详细信息
In this article, we formulate a monolithic optimal control method for general timedependent fluid-structure interaction (FSI) systems with large solid deformation: we consider a displacement-tracking type of objective with a constraint of the solid velocity, and tackle the time-dependent control problems by a piecewise-in-time control method;we cope with the large solid displacement using a one-velocity fictitious domain method, and solve the fully-coupled FSI and the corresponding adjoint equations in a monolithic manner. The proposed method is implemented in open-source software package FreeFEM++ and assessed by three numerical experiments, in the aspects of stability of the numerical scheme for different regularisation parameters, and efficiency of reducing the objective function with control of the solid velocity. (c) 2022 Elsevier Ltd. All rights reserved.
Blood flow, dam or ship construction and numerous other problems in biomedical and general engineering involve incompressible flows interacting with elastic structures. Such interactions heavily influence the deformat...
详细信息
Blood flow, dam or ship construction and numerous other problems in biomedical and general engineering involve incompressible flows interacting with elastic structures. Such interactions heavily influence the deformation and stress states which, in turn, affect the engineering design process. Therefore, any reliable model of such physical processes must consider the coupling of fluids and solids. However, complexity increases for non-Newtonian fluid models, as used, e.g., for blood or polymer flows. In these fluids, subtle differences in the local shear rate can have a drastic impact on the flow and hence on the coupled problem. There, existing (semi-) implicit solution strategies based on split-step or projection schemes for Newtonian fluids are not applicable, while extensions to non-Newtonian fluids can lead to substantial numerical overhead depending on the chosen fluid solver. To address these shortcomings, we present here a higher-order accurate, added-mass-stable fluid-structure interaction scheme centered around a split-step fluid solver. We compare several implicit and semi-implicit variants of the algorithm and verify convergence in space and time. Numerical examples show good performance in both benchmarks and an idealised setting of blood flow through an abdominal aortic aneurysm considering physiological parameters. (C) 2021 The Author(s). Published by Elsevier Ltd.
A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calcu...
详细信息
A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calculation of fluid and structure deformation by adding Lagrangian particles to Eulerian grids. The overall computational domain is solved by the Eulerian method, and the Lagrangian particles with specified volume and influence domains are used to track structural deformations. The bidirectional mapping of physical quantities is achieved using the weighted average of the influence domain, which are based on the topological relationship between Eulerian grids and Lagrangian particles. Then, the data dependence solution and parallel algorithm realization are presented for the large-scale numerical calculations of explosion and impact problems. Additionally, the corresponding parallel program is developed based on the message passing interface(MPI) standard, and the parallel efficiency of parallel hydrocode are tested. The numerical results of typical explosion and impact problems are compared with corresponding experimental data to verify the effectiveness of the method. These comparisons show that the embedded EulerianLagrangian method successfully combine the advantages of both the Eulerian and Lagrangian methods to efficiently calculate the processes of large deformation and dynamic damage to the materials. The results presented in this work provide a useful reference point for further research on explosion and impact problems.
Objective: To explore the role of hemodynamic factors in the occurrence of superior mesenteric artery (SMA) dissection (SMAD) using a fluid-structure interaction (FSI) simulation method, and to identify histopathologi...
详细信息
Objective: To explore the role of hemodynamic factors in the occurrence of superior mesenteric artery (SMA) dissection (SMAD) using a fluid-structure interaction (FSI) simulation method, and to identify histopathologic changes occurring in the wall of the SMA. Methods: A total of 122 consecutive patients diagnosed with SMAD and 122 controls were included in this study. Hemodynamic factors were calculated using a FSI simulation method. Additionally, SMA spec-imens obtained from 12 cadavers were stained for histological quantitative ***: The mean aortomesenteric angle (59.7 degrees +/- 21.4 degrees vs 48.2 degrees +/- 16.8 degrees;p < .001) and SMA maximum curvature (0.084 +/- 0.078 mm-1 vs 0.032 +/- 0.023 mm -1;p < .001) were higher in SMAD patients than the controls. Larger aortomesenteric angles and SMA curvatures were associated with higher and more concentrated wall shear stress at anterior wall of the SMA curve segment, co-located with the dissec-tion origins. The mean thickness of media (325.18 +/- 44.87 mu m vs 556.92 +/- 138.32 mu m;p = .003) was thinner in the anterior wall of the SMA curve than in the posterior wall. The area fractions of elastin (17.96% +/- 3.36% vs 27.06% +/- 4.18%;p = .002) and collagen (45.43% +/- 6.89% vs 55.57% +/- 7.57%;p = .036) were lower in anterior wall of the SMA curve than in posterior ***: Increased aortomesenteric angle and SMA curvature are risk factors for SMAD. Both of these factors can cause local hemodynamic abnormalities, which can lead to histopathologic changes in anterior wall of SMA.(c) 2022 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/4.0/ )
The paper presents the multidisciplinary design optimization (MDO) of a deep-V planing-hull grillage panel subject to slamming loads in regular waves. Namely, fluidstructureinteraction (FSI) experiments, computation...
详细信息
The paper presents the multidisciplinary design optimization (MDO) of a deep-V planing-hull grillage panel subject to slamming loads in regular waves. Namely, fluidstructureinteraction (FSI) experiments, computations, and MDO are presented and discussed for a bottom-panel grillage of a high-speed Generic Prismatic Planing Hull in regular waves. Computations are per-formed via one-and tightly coupled two-way computational fluid and structural dynamics (CFD/ CSD) using unsteady Reynolds-averaged Navier-Stokes equation solvers to compute the hydro-dynamic loads. The structural assessment of the original/traditional grillage is performed using a fully parametric finite element (FE) model, showing the significant effects of the FE boundary conditions on the structural response. Firstly, an equivalent static and uniform load is identified via CFD and applied during optimization using two design spaces. The selected optimized design provides a grillage-weight reduction of 35% and an overall factor of safety equal to 1.72. The optimized design presents variations of stiffeners dimensions across the grillage with the largest stiffener at the middle, distributing the stress more uniformly among the stiffeners. The effects of one-versus two-way coupling are negligible for both the original/traditional and optimized grillages (as per the hydroelasticity factor R), whereas the effects of FE boundary conditions on the analysis and optimization outcomes are significant, confirming the need for proper calibration of the FE model in FSI and MDO studies. Secondly, MDO is performed with a dynamic load applied via one-way coupling FSI. An additional 5% weight reduction is identified, achieving a 40% weight reduction compared to the traditional grillage. The optimal design presents the largest stiffener close to the keel, which is significantly different than the design obtained for uniform/static load. Comparison of computational and experimental data is very good, indicating that th
We present a stability analysis of the Discontinuous Galerkin method on polygonal and polyhedral meshes (PolyDG) for the Stokes problem. In particular, we analyze the discrete inf-sup condition for different choices o...
详细信息
We present a stability analysis of the Discontinuous Galerkin method on polygonal and polyhedral meshes (PolyDG) for the Stokes problem. In particular, we analyze the discrete inf-sup condition for different choices of the polynomial approximation order of the velocity and pressure approximation spaces. To this aim, we employ a generalized inf-sup condition with a pressure stabilization term. We also prove a priori hp-version error estimates in suitable norms. We numerically check the behaviour of the inf-sup constant and the order of convergence with respect to the mesh configuration, the mesh-size, and the polynomial degree. Finally, as a relevant application of our analysis, we consider the PolyDG approximation for a 2D fluid-structure interaction problem and we numerically explore the stability properties of the method.
The multi-resolution method, e.g., the Adaptive Particle Refinement (APR) method, has been developed to increase the local particle resolution and therefore the solution quality within a pre-defined refinement zone in...
详细信息
The multi-resolution method, e.g., the Adaptive Particle Refinement (APR) method, has been developed to increase the local particle resolution and therefore the solution quality within a pre-defined refinement zone instead of using a globally uniform resolution for Smoothed Particle Hydrodynamics (SPH). However, sometimes, the targeted zone of interest can be varying, and the corresponding topology is very complex, thus the conventional APR method is not able to track these characteristics adaptively. In this study, a novel Block-based Adaptive Particle Refinement (BAPR) method is developed, which is able to provide the necessary local refinement flexibly for any targeted characteristic, and track it adaptively. In BAPR, the so-called activation status of the block array defines the refinement regions, where the transition and activated zones are determined accordingly. A regularization method for the generated particles in the newly activated blocks is developed to render an isotropic distribution of these new particles. The proposed method has been deployed for simulating fluid-structure interaction
Technical mitigation devices have been demonstrated to be efficient for decreasing the peak and mean wind pressures on the roof systems of low-rise buildings through disrupting vortex for-mation. Previous investigatio...
详细信息
Technical mitigation devices have been demonstrated to be efficient for decreasing the peak and mean wind pressures on the roof systems of low-rise buildings through disrupting vortex for-mation. Previous investigations have regarded mitigation devices as rigid bodies and neglected aeroelastic behavior. However, the relatively large uplift force induces small-stiffness mitigation devices vibration and introduces a considerable aeroelastic effect. These effects might result in an increase in negative wind pressure on the building roof and the aeroelastic instability problems of mitigation devices. Therefore, the current study conducted two-way fluid-structure interaction (FSI) numerical simulation for a low-rise flat-roof building with mitigation devices (metal and polymer spoilers) to investigate the influence of the aeroelastic effect of spoiler on the wind pressure on the building roof. The large eddy simulation (LES) method and the finite element analysis (FEA) model are used to solve the flow field and the structural responses of the spoiler, respectively. The computational model and method are verified by comparing the calculation results with previous experimental results and benchmark computational tests. The dynamic response of spoiler, the flow instabilities and the wind pressure on the building roof are analyzed. The influence mechanism of the spoiler vibration on the flow field over the building, the wind pressure on the building roof and the uplift force on the spoiler were revealed considering the influencing factors, including the vibrational amplitude of spoiler and the wind direction. The results show that the aeroelastic effect of spoiler increases the negative wind pressure on the roof, and the larger the vibrational amplitude of spoiler is, the more the negative wind pressure in-creases. The negative wind pressure on the eave also increases with increasing wind direction angle. The main reason for this phenomenon is that the spoiler vibration incr
This paper presents a partitioned framework for the numerical simulation of fluid-structure interactions by coupling the lattice Boltzmann method (LBM) and the finite element method (FEM). The two numerical methods LB...
详细信息
This paper presents a partitioned framework for the numerical simulation of fluid-structure interactions by coupling the lattice Boltzmann method (LBM) and the finite element method (FEM). The two numerical methods LBM and FEM are coupled with an implicit immersed boundary method (IBM) in a strong way, which ensures exactly the no-slip condition and the continuities of velocity and stress at the fluid-solid interface and each instant in time. In the proposed partitioned coupling procedure, the coupling system of equations are first established and then condensed to the interface. By solving the condensed coupling system of equations, the interface force field is obtained and sent to both solvers to accomplish time integrations in each sub-domain. In addition, two strategies based on linear interpolation in time are proposed to handle the cases with non-conforming time-steps in the fluid and solid sub-domains. Through several 2D and 3D numerical test-cases on the mechanical heart valve, the fluid-induced vibration of a deformable solid beam, the flapping flag, the proposed coupling framework is validated with good agreements with references. Finally, a test-case on the interaction between the blood flow and the aortic valve is carried out, showing the applicability of the present framework in realistic biomechanical applications. (C) 2021 Elsevier Inc. All rights reserved.
暂无评论