317,706 research outputs found
A multiscale-multiphysics strategy for numerical modeling of thin piezoelectric sheets
Flexible piezoelectric devices made of polymeric materials are widely used
for micro- and nano-electro-mechanical systems. In particular, numerous recent
applications concern energy harvesting. Due to the importance of computational
modeling to understand the influence that microscale geometry and constitutive
variables exert on the macroscopic behavior, a numerical approach is developed
here for multiscale and multiphysics modeling of piezoelectric materials made
of aligned arrays of polymeric nanofibers. At the microscale, the
representative volume element consists in piezoelectric polymeric nanofibers,
assumed to feature a linear piezoelastic constitutive behavior and subjected to
electromechanical contact constraints using the penalty method. To avoid the
drawbacks associated with the non-smooth discretization of the master surface,
a contact smoothing approach based on B\'ezier patches is extended to the
multiphysics framework providing an improved continuity of the
parameterization. The contact element contributions to the virtual work
equations are included through suitable electric, mechanical and coupling
potentials. From the solution of the micro-scale boundary value problem, a
suitable scale transition procedure leads to the formulation of a macroscopic
thin piezoelectric shell element.Comment: 11 pages, 6 pages, 21 reference
Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
We present a new family of high order accurate fully discrete one-step
Discontinuous Galerkin (DG) finite element schemes on moving unstructured
meshes for the solution of nonlinear hyperbolic PDE in multiple space
dimensions, which may also include parabolic terms in order to model
dissipative transport processes. High order piecewise polynomials are adopted
to represent the discrete solution at each time level and within each spatial
control volume of the computational grid, while high order of accuracy in time
is achieved by the ADER approach. In our algorithm the spatial mesh
configuration can be defined in two different ways: either by an isoparametric
approach that generates curved control volumes, or by a piecewise linear
decomposition of each spatial control volume into simplex sub-elements. Our
numerical method belongs to the category of direct
Arbitrary-Lagrangian-Eulerian (ALE) schemes, where a space-time conservation
formulation of the governing PDE system is considered and which already takes
into account the new grid geometry directly during the computation of the
numerical fluxes. Our new Lagrangian-type DG scheme adopts the novel a
posteriori sub-cell finite volume limiter method, in which the validity of the
candidate solution produced in each cell by an unlimited ADER-DG scheme is
verified against a set of physical and numerical detection criteria. Those
cells which do not satisfy all of the above criteria are flagged as troubled
cells and are recomputed with a second order TVD finite volume scheme. The
numerical convergence rates of the new ALE ADER-DG schemes are studied up to
fourth order in space and time and several test problems are simulated.
Finally, an application inspired by Inertial Confinement Fusion (ICF) type
flows is considered by solving the Euler equations and the PDE of viscous and
resistive magnetohydrodynamics (VRMHD).Comment: 39 pages, 21 figure
An adaptive Cartesian embedded boundary approach for fluid simulations of two- and three-dimensional low temperature plasma filaments in complex geometries
We review a scalable two- and three-dimensional computer code for
low-temperature plasma simulations in multi-material complex geometries. Our
approach is based on embedded boundary (EB) finite volume discretizations of
the minimal fluid-plasma model on adaptive Cartesian grids, extended to also
account for charging of insulating surfaces. We discuss the spatial and
temporal discretization methods, and show that the resulting overall method is
second order convergent, monotone, and conservative (for smooth solutions).
Weak scalability with parallel efficiencies over 70\% are demonstrated up to
8192 cores and more than one billion cells. We then demonstrate the use of
adaptive mesh refinement in multiple two- and three-dimensional simulation
examples at modest cores counts. The examples include two-dimensional
simulations of surface streamers along insulators with surface roughness; fully
three-dimensional simulations of filaments in experimentally realizable
pin-plane geometries, and three-dimensional simulations of positive plasma
discharges in multi-material complex geometries. The largest computational
example uses up to million mesh cells with billions of unknowns on
computing cores. Our use of computer-aided design (CAD) and constructive solid
geometry (CSG) combined with capabilities for parallel computing offers
possibilities for performing three-dimensional transient plasma-fluid
simulations, also in multi-material complex geometries at moderate pressures
and comparatively large scale.Comment: 40 pages, 21 figure
A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting
An explicit moving boundary method for the numerical solution of
time-dependent hyperbolic conservation laws on grids produced by the
intersection of complex geometries with a regular Cartesian grid is presented.
As it employs directional operator splitting, implementation of the scheme is
rather straightforward. Extending the method for static walls from Klein et
al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme
calculates fluxes needed for a conservative update of the near-wall cut-cells
as linear combinations of standard fluxes from a one-dimensional extended
stencil. Here the standard fluxes are those obtained without regard to the
small sub-cell problem, and the linear combination weights involve detailed
information regarding the cut-cell geometry. This linear combination of
standard fluxes stabilizes the updates such that the time-step yielding
marginal stability for arbitrarily small cut-cells is of the same order as that
for regular cells. Moreover, it renders the approach compatible with a wide
range of existing numerical flux-approximation methods. The scheme is extended
here to time dependent rigid boundaries by reformulating the linear combination
weights of the stabilizing flux stencil to account for the time dependence of
cut-cell volume and interface area fractions. The two-dimensional tests
discussed include advection in a channel oriented at an oblique angle to the
Cartesian computational mesh, cylinders with circular and triangular
cross-section passing through a stationary shock wave, a piston moving through
an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil
profile.Comment: 30 pages, 27 figures, 3 table
Mechanistic and pathological study of the genesis, growth, and rupture of abdominal aortic aneurysms
Postprint (published version
A computational framework for the morpho-elastic development of molluskan shells by surface and volume growth
Mollusk shells are an ideal model system for understanding the morpho-elastic
basis of morphological evolution of invertebrates' exoskeletons. During the
formation of the shell, the mantle tissue secretes proteins and minerals that
calcify to form a new incremental layer of the exoskeleton. Most of the
existing literature on the morphology of mollusks is descriptive. The
mathematical understanding of the underlying coupling between pre-existing
shell morphology, de novo surface deposition and morpho-elastic volume growth
is at a nascent stage, primarily limited to reduced geometric representations.
Here, we propose a general, three-dimensional computational framework coupling
pre-existing morphology, incremental surface growth by accretion, and
morpho-elastic volume growth. We exercise this framework by applying it to
explain the stepwise morphogenesis of seashells during growth: new material
surfaces are laid down by accretive growth on the mantle whose form is
determined by its morpho-elastic growth. Calcification of the newest surfaces
extends the shell as well as creates a new scaffold that constrains the next
growth step. We study the effects of surface and volumetric growth rates, and
of previously deposited shell geometries on the resulting modes of mantle
deformation, and therefore of the developing shell's morphology. Connections
are made to a range of complex shells ornamentations.Comment: Main article is 20 pages long with 15 figures. Supplementary material
is 4 pages long with 6 figures and 6 attached movies. To be published in PLOS
Computational Biolog
- …