Jagged-time-step technique improving convergence order of Fernandezs Explicit Robin-Neumann scheme for the coupling of incompressible fluid with thin-walled structure
Inspired by Rybaks multiple-time-step technique, jagged-time-step technique is proposed and applied to Fernandezs Explicit Robin-Neumann scheme. For some instances, numerical experiments demonstrate higher convergence orders and accuracy with lower computation cost as time and space get refined.
We consider a fully discrete loosely coupled scheme for incompressible fluid-structure interaction based on the time semi-discrete splitting method introduced in {emph{[Burman, Durst & Guzman, arXiv:1911.06760]}}. The splittling method uses a Robin-Robin type coupling that allows for a segregated solution of the solid and the fluid systems, without inner iterations. For the discretisation in space we consider piecewise affine continuous finite elements for all the fields and ensure the inf-sup condition by using a Brezzi-Pitkaranta type pressure stabilization. The interfacial fluid-stresses are evaluated in a variationally consistent fashion, that is shown to admit an equivalent Lagrange multiplier formulation. We prove that the method is unconditionally stable and robust with respect to the amount of added-mass in the system. Furthermore, we provide an error estimate that shows the error in the natural energy norm for the system is $mathcal Obig(sqrt{T}(sqrt{Delta t} + h)big)$ where $T$ is the final time, $Delta t$ the time-step length and $h$ the space discretization parameter.
This work aims at providing some novel and practical ideas to improve accuracy of some partitioned algorithms, precisely Fernandezs Explicit Robin-Neumann and fully decoupled schemes, for the coupling of incompressible fluid with thin-walled structure. Inspired by viscosity of fluid and justified by boundary layer theory, the force between fluid and structure corresponding to viscosity is increased. Numerical experiments demonstrate improvement of accuracy under such modification. To improve accuracy of fully decoupled schemes further, the underlying projection method is replaced.
We present a novel preconditioning technique for Krylov subspace algorithms to solve fluid-structure interaction (FSI) linearized systems arising from finite element discretizations. An outer Krylov subspace solver preconditioned with a geometric multigrid (GMG) algorithm is used, where for the multigrid level sub-solvers, a field-split (FS) preconditioner is proposed. The block structure of the FS preconditioner is derived using the physical variables as splitting strategy. To solve the subsystems originated by the FS preconditioning, an additive Schwarz (AS) block strategy is employed. The proposed field-split preconditioner is tested on biomedical FSI applications. Both 2D and 3D simulations are carried out considering aneurysm and venous valve geometries. The performance of the FS preconditioner is compared with that of a second preconditioner of pure domain decomposition type.
This paper continues to study the explicit two-stage fourth-order accurate time discretiza- tions [5, 7]. By introducing variable weights, we propose a class of more general explicit one-step two-stage time discretizations, which are different from the existing methods, such as the Euler methods, Runge-Kutta methods, and multistage multiderivative methods etc. We study the absolute stability, the stability interval, and the intersection between the imaginary axis and the absolute stability region. Our results show that our two-stage time discretizations can be fourth-order accurate conditionally, the absolute stability region of the proposed methods with some special choices of the variable weights can be larger than that of the classical explicit fourth- or fifth-order Runge-Kutta method, and the interval of absolute stability can be almost twice as much as the latter. Several numerical experiments are carried out to demonstrate the performance and accuracy as well as the stability of our proposed methods
We develop a novel iterative solution method for the incompressible Navier-Stokes equations with boundary conditions coupled with reduced models. The iterative algorithm is designed based on the variational multiscale formulation and the generalized-$alpha$ scheme. The spatiotemporal discretization leads to a block structure of the resulting consistent tangent matrix in the Newton-Raphson procedure. As a generalization of the conventional block preconditioners, a three-level nested block preconditioner is introduced to attain a better representation of the Schur complement, which plays a key role in the overall algorithm robustness and efficiency. This approach provides a flexible, algorithmic way to handle the Schur complement for problems involving multiscale and multiphysics coupling. The solution method is implemented and benchmarked against experimental data from the nozzle challenge problem issued by the US Food and Drug Administration. The robustness, efficiency, and parallel scalability of the proposed technique are then examined in several settings, including moderately high Reynolds number flows and physiological flows with strong resistance effect due to coupled downstream vasculature models. Two patient-specific hemodynamic simulations, covering systemic and pulmonary flows, are performed to further corroborate the efficacy of the proposed methodology.
Yiyi Huang
.
(2020)
.
"Jagged-time-step technique improving convergence order of Fernandezs Explicit Robin-Neumann scheme for the coupling of incompressible fluid with thin-walled structure"
.
Yiyi Huang
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا