12,708 research outputs found
Cut Finite Elements for Convection in Fractured Domains
We develop a cut finite element method (CutFEM) for the convection problem in
a so called fractured domain which is a union of manifolds of different
dimensions such that a dimensional component always resides on the boundary
of a dimensional component. This type of domain can for instance be used
to model porous media with embedded fractures that may intersect. The
convection problem can be formulated in a compact form suitable for analysis
using natural abstract directional derivative and divergence operators. The cut
finite element method is based on using a fixed background mesh that covers the
domain and the manifolds are allowed to cut through a fixed background mesh in
an arbitrary way. We consider a simple method based on continuous piecewise
linear elements together with weak enforcement of the coupling conditions and
stabilization. We prove a priori error estimates and present illustrating
numerical examples
On the role of confinement on solidification in pure materials and binary alloys
We use a phase-field model to study the effect of confinement on dendritic
growth, in a pure material solidifying in an undercooled melt, and in the
directional solidification of a dilute binary alloy. Specifically, we observe
the effect of varying the vertical domain extent () on tip selection,
by quantifying the dendrite tip velocity and curvature as a function of
, and other process parameters. As decreases, we find that the
operating state of the dendrite tips becomes significantly affected by the
presence of finite boundaries. For particular boundary conditions, we observe a
switching of the growth state from 3-D to 2-D at very small , in both
the pure material and alloy. We demonstrate that results from the alloy model
compare favorably with those from an experimental study investigating this
effect.Comment: 13 pages, 9 figures, 3 table
A Combined Preconditioning Strategy for Nonsymmetric Systems
We present and analyze a class of nonsymmetric preconditioners within a
normal (weighted least-squares) matrix form for use in GMRES to solve
nonsymmetric matrix problems that typically arise in finite element
discretizations. An example of the additive Schwarz method applied to
nonsymmetric but definite matrices is presented for which the abstract
assumptions are verified. A variable preconditioner, combining the original
nonsymmetric one and a weighted least-squares version of it, is shown to be
convergent and provides a viable strategy for using nonsymmetric
preconditioners in practice. Numerical results are included to assess the
theory and the performance of the proposed preconditioners.Comment: 26 pages, 3 figure
New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid
Seismic imaging of the mantle has revealed large and small scale
heterogeneities in the lower mantle; specifically structures known as large low
shear velocity provinces (LLSVP) below Africa and the South Pacific. Most
interpretations propose that the heterogeneities are compositional in nature,
differing in composition from the overlying mantle, an interpretation that
would be consistent with chemical geodynamic models. Numerical modeling of
persistent compositional interfaces presents challenges, even to
state-of-the-art numerical methodology. For example, some numerical algorithms
for advecting the compositional interface cannot maintain a sharp compositional
boundary as the fluid migrates and distorts with time dependent fingering due
to the numerical diffusion that has been added in order to maintain the upper
and lower bounds on the composition variable and the stability of the advection
method. In this work we present two new algorithms for maintaining a sharper
computational boundary than the advection methods that are currently openly
available to the computational mantle convection community; namely, a
Discontinuous Galerkin method with a Bound Preserving limiter and a
Volume-of-Fluid interface tracking algorithm. We compare these two new methods
with two approaches commonly used for modeling the advection of two distinct,
thermally driven, compositional fields in mantle convection problems; namely,
an approach based on a high-order accurate finite element method advection
algorithm that employs an artificial viscosity technique to maintain the upper
and lower bounds on the composition variable as well as the stability of the
advection algorithm and the advection of particles that carry a scalar quantity
representing the location of each compositional field. All four of these
algorithms are implemented in the open source FEM code ASPECT
A "poor man's" approach to topology optimization of natural convection problems
Topology optimization of natural convection problems is computationally
expensive, due to the large number of degrees of freedom (DOFs) in the model
and its two-way coupled nature. Herein, a method is presented to reduce the
computational effort by use of a reduced-order model governed by simplified
physics. The proposed method models the fluid flow using a potential flow
model, which introduces an additional fluid property. This material property
currently requires tuning of the model by comparison to numerical Navier-Stokes
based solutions. Topology optimization based on the reduced-order model is
shown to provide qualitatively similar designs, as those obtained using a full
Navier-Stokes based model. The number of DOFs is reduced by 50% in two
dimensions and the computational complexity is evaluated to be approximately
12.5% of the full model. We further compare to optimized designs obtained
utilizing Newton's convection law.Comment: Preprint version. Please refer to final version in Structural
Multidisciplinary Optimization https://doi.org/10.1007/s00158-019-02215-
- …