-
Recovery by discretization corrected particle strength exchange (DC PSE) operators
Authors:
Benjamin F. Zwick,
George C. Bourantas,
Farah Alkhatib,
Adam Wittek,
Karol Miller
Abstract:
A new recovery technique based on discretization corrected particle strength exchange (DC PSE) operators is developed in this paper. DC PSE is a collocation method that can be used to compute derivatives directly at nodal points, instead of by projection from Gauss points as is done in many finite element-based recovery techniques. The proposed method is truly meshless and does not require patches…
▽ More
A new recovery technique based on discretization corrected particle strength exchange (DC PSE) operators is developed in this paper. DC PSE is a collocation method that can be used to compute derivatives directly at nodal points, instead of by projection from Gauss points as is done in many finite element-based recovery techniques. The proposed method is truly meshless and does not require patches of elements to be defined, which makes it generally applicable to point clouds and arbitrary element topologies. Numerical examples show that the proposed method is accurate and robust.
△ Less
Submitted 15 April, 2022;
originally announced April 2022.
-
Intracranial hemodynamics simulations: An efficient and accurate immersed boundary scheme
Authors:
D. S. Lampropoulos,
G. C. Bourantas,
B. F. Zwick,
G. C. Kagadis,
A. Wittek,
K. Miller,
V. C. Loukopoulos
Abstract:
Computational fluid dynamics (CFD) studies have been increasingly used for blood flow simulations in intracranial aneurysms (ICAs). However, despite the continuous progress of body-fitted CFD solvers, generating a high quality mesh is still the bottleneck of the CFD simulation, and strongly affects the accuracy of the numerical solution. To overcome this challenge, which will allow preforming CFD…
▽ More
Computational fluid dynamics (CFD) studies have been increasingly used for blood flow simulations in intracranial aneurysms (ICAs). However, despite the continuous progress of body-fitted CFD solvers, generating a high quality mesh is still the bottleneck of the CFD simulation, and strongly affects the accuracy of the numerical solution. To overcome this challenge, which will allow preforming CFD simulations efficiently for a large number of aneurysm cases we use an Immersed Boundary (IB) method. The proposed scheme relies on Cartesian grids to solve the incompressible Navier-Stokes (N-S) equations, using a finite element solver, and Lagrangian points to discretize the immersed object. All grid generations are conducted through automated algorithms which require no user input. Consequently, we verify the proposed method by comparing our numerical findings (velocity values) with published experimental results. Finally, we test the ability of the scheme to efficiently handle hemodynamic simulations on complex geometries on a sample of four patient-specific intracranial aneurysms.
△ Less
Submitted 27 July, 2020;
originally announced July 2020.
-
Immersed boundary finite element method for blood flow simulation
Authors:
G. C. Bourantas,
D. L. Lampropoulos,
B. F. Zwick,
V. C. Loukopoulos,
A. Wittek,
K. Miller
Abstract:
We present an efficient and accurate immersed boundary (IB) finite element (FE) solver for numerically solving incompressible Navier--Stokes equations. Particular emphasis is given to internal flows with complex geometries (blood flow in the vasculature system). IB methods are computationally costly for internal flows, mainly due to the large percentage of grid points that lie outside the flow dom…
▽ More
We present an efficient and accurate immersed boundary (IB) finite element (FE) solver for numerically solving incompressible Navier--Stokes equations. Particular emphasis is given to internal flows with complex geometries (blood flow in the vasculature system). IB methods are computationally costly for internal flows, mainly due to the large percentage of grid points that lie outside the flow domain. In this study, we apply a local refinement strategy, along with a domain reduction approach in order to reduce the grid that covers the flow domain and increase the percentage of the grid nodes that fall inside the flow domain. The proposed method utilizes an efficient and accurate FE solver with the incremental pressure correction scheme (IPCS), along with the boundary condition enforced IB method to numerically solve the transient, incompressible Navier--Stokes flow equations. We verify the accuracy of the numerical method using the analytical solution for Poiseuille flow in a cylinder. We further examine the accuracy and applicability of the proposed method by considering flow within complex geometries, such as blood flow in aneurysmal vessels and the aorta, flow configurations which would otherwise be extremely difficult to solve by most IB methods. Our method offers high accuracy, as demonstrated by the verification examples, and high efficiency, as demonstrated through the solution of blood flow within complex geometry on an off-the-shelf laptop computer.
△ Less
Submitted 4 July, 2020;
originally announced July 2020.
-
Simple and robust element-free Galerkin method with interpolating shape functions for finite deformation elasticity
Authors:
George Bourantas,
Benjamin F. Zwick,
Grand Joldes,
Adam Wittek,
Karol Miller
Abstract:
In this paper, we present a meshless method belonging to the family of element-free Galerkin (EFG) methods. The distinguishing feature of the presented meshless method is that it allows accurate enforcement of essential boundary conditions. The method uses total Lagrangian formulation with explicit time integration to facilitate code simplicity and robust computations in applications that involve…
▽ More
In this paper, we present a meshless method belonging to the family of element-free Galerkin (EFG) methods. The distinguishing feature of the presented meshless method is that it allows accurate enforcement of essential boundary conditions. The method uses total Lagrangian formulation with explicit time integration to facilitate code simplicity and robust computations in applications that involve large deformations and non-linear materials. We use a regularized weight function, which closely approximates the Kronecker delta, to generate interpolating shape functions. The imposition of the prescribed displacements on the boundary becomes as straightforward as in the finite element (FE) method. The effectiveness and accuracy of the proposed method is demonstrated using 3D numerical examples that include cylinder indentation by 70% of its initial height, and indentation of the brain.
△ Less
Submitted 14 May, 2020;
originally announced May 2020.