72,339 research outputs found
Recommended from our members
A hybrid stabilization technique for simulating water wave - Structure interaction by incompressible Smoothed Particle Hydrodynamics (ISPH) method
The Smoothed Particle Hydrodynamics (SPH) method is emerging as a potential tool for studying water wave related problems, especially for violent free surface flow and large deformation problems. The incompressible SPH (ISPH) computations have been found not to be able to maintain the stability in certain situations and there exist some spurious oscillations in the pressure time history, which is similar to the weakly compressible SPH (WCSPH). One main cause of this problem is related to the non-uniform and clustered distribution of the moving particles. In order to improve the model performance, the paper proposed an efficient hybrid numerical technique aiming to correct the ill particle distributions. The correction approach is realized through the combination of particle shifting and pressure gradient improvement. The advantages of the proposed hybrid technique in improving ISPH calculations are demonstrated through several applications that include solitary wave impact on a slope or overtopping a seawall, and regular wave slamming on the subface of open-piled structure
Intimal and medial contributions to the hydraulic resistance of the arterial wall at different pressures: a combined computational and experimental study
The hydraulic resistances of the intima and media determine water flux and the advection of macromolecules into and across the arterial wall. Despite several experimental and computational studies, these transport processes and their dependence on transmural pressure remain incompletely understood. Here, we use a combination of experimental and computational methods to ascertain how the hydraulic permeability of the rat abdominal aorta depends on these two layers and how it is affected by structural rearrangement of the media under pressure. Ex vivo experiments determined the conductance of the whole wall, the thickness of the media and the geometry of medial smooth muscle cells (SMCs) and extracellular matrix (ECM). Numerical methods were used to compute water flux through the media. Intimal values were obtained by subtraction. A mechanism was identified that modulates pressure-induced changes in medial transport properties: compaction of the ECM leading to spatial reorganization of SMCs. This is summarized in an empirical constitutive law for permeability and volumetric strain. It led to the physiologically interesting observation that, as a consequence of the changes in medial microstructure, the relative contributions of the intima and media to the hydraulic resistance of the wall depend on the applied pressure; medial resistance dominated at pressures above approximately 93 mmHg in this vessel
Diffusion-Based Coarse Graining in Hybrid Continuum-Discrete Solvers: Applications in CFD-DEM
In this work, a coarse-graining method previously proposed by the authors in
a companion paper based on solving diffusion equations is applied to CFD-DEM
simulations, where coarse graining is used to obtain solid volume fraction,
particle phase velocity, and fluid-particle interaction forces. By examining
the conservation requirements, the variables to solve diffusion equations for
in CFD-DEM simulations are identified. The algorithm is then implemented into a
CFD-DEM solver based on OpenFOAM and LAMMPS, the former being a
general-purpose, three-dimensional CFD solver based on unstructured meshes.
Numerical simulations are performed for a fluidized bed by using the CFD-DEM
solver with the diffusion-based coarse-graining algorithm. Converged results
are obtained on successively refined meshes, even for meshes with cell sizes
comparable to or smaller than the particle diameter. This is a critical
advantage of the proposed method over many existing coarse-graining methods,
and would be particularly valuable when small cells are required in part of the
CFD mesh to resolve certain flow features such as boundary layers in wall
bounded flows and shear layers in jets and wakes. Moreover, we demonstrate that
the overhead computational costs incurred by the proposed coarse-graining
procedure are a small portion of the total costs in typical CFD-DEM simulations
as long as the number of particles per cell is reasonably large, although
admittedly the computational overhead of the coarse graining often exceeds that
of the CFD solver. Other advantages of the present algorithm include more
robust and physically realistic results, flexibility and easy implementation in
almost any CFD solvers, and clear physical interpretation of the computational
parameter needed in the algorithm. In summary, the diffusion-based method is a
theoretically elegant and practically viable option for CFD-DEM simulations
High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
We present a new family of very high order accurate direct
Arbitrary-Lagrangian-Eulerian (ALE) Finite Volume (FV) and Discontinuous
Galerkin (DG) schemes for the solution of nonlinear hyperbolic PDE systems on
moving 2D Voronoi meshes that are regenerated at each time step and which
explicitly allow topology changes in time.
The Voronoi tessellations are obtained from a set of generator points that
move with the local fluid velocity. We employ an AREPO-type approach, which
rapidly rebuilds a new high quality mesh rearranging the element shapes and
neighbors in order to guarantee a robust mesh evolution even for vortex flows
and very long simulation times. The old and new Voronoi elements associated to
the same generator are connected to construct closed space--time control
volumes, whose bottom and top faces may be polygons with a different number of
sides. We also incorporate degenerate space--time sliver elements, needed to
fill the space--time holes that arise because of topology changes. The final
ALE FV-DG scheme is obtained by a redesign of the fully discrete direct ALE
schemes of Boscheri and Dumbser, extended here to moving Voronoi meshes and
space--time sliver elements. Our new numerical scheme is based on the
integration over arbitrary shaped closed space--time control volumes combined
with a fully-discrete space--time conservation formulation of the governing PDE
system. In this way the discrete solution is conservative and satisfies the GCL
by construction.
Numerical convergence studies as well as a large set of benchmarks for
hydrodynamics and magnetohydrodynamics (MHD) demonstrate the accuracy and
robustness of the proposed method. Our numerical results clearly show that the
new combination of very high order schemes with regenerated meshes with
topology changes lead to substantial improvements compared to direct ALE
methods on conforming meshes
Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures
The voxel finite cell method employs unfitted finite element meshes and voxel quadrature rules to seamlessly
transfer CT data into patient-specific bone discretizations. The method, however, still requires the explicit
parametrization of boundary surfaces to impose traction and displacement boundary conditions, which
constitutes a potential roadblock to automation. We explore a phase-field based formulation for imposing
traction and displacement constraints in a diffuse sense. Its essential component is a diffuse geometry model
generated from metastable phase-field solutions of the Allen-Cahn problem that assumes the imaging data as
initial condition. Phase-field approximations of the boundary and its gradient are then employed to transfer
all boundary terms in the variational formulation into volumetric terms. We show that in the context of the
voxel finite cell method, diffuse boundary conditions achieve the same accuracy as boundary conditions
defined over explicit sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field,
the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method
by analyzing stresses in a human femur and a vertebral body
- …