11,597 research outputs found
A "poor man's" approach for high-resolution three-dimensional topology optimization of natural convection problems
This paper treats topology optimization of natural convection problems. A
simplified model is suggested to describe the flow of an incompressible fluid
in steady state conditions, similar to Darcy's law for fluid flow in porous
media. The equations for the fluid flow are coupled to the thermal
convection-diffusion equation through the Boussinesq approximation. The coupled
non-linear system of equations is discretized with stabilized finite elements
and solved in a parallel framework that allows for the optimization of high
resolution three-dimensional problems. A density-based topology optimization
approach is used, where a two-material interpolation scheme is applied to both
the permeability and conductivity of the distributed material. Due to the
simplified model, the proposed methodology allows for a significant reduction
of the computational effort required in the optimization. At the same time, it
is significantly more accurate than even simpler models that rely on convection
boundary conditions based on Newton's law of cooling. The methodology discussed
herein is applied to the optimization-based design of three-dimensional heat
sinks. The final designs are formally compared with results of previous work
obtained from solving the full set of Navier-Stokes equations. The results are
compared in terms of performance of the optimized designs and computational
cost. The computational time is shown to be decreased to around 5-20% in terms
of core-hours, allowing for the possibility of generating an optimized design
during the workday on a small computational cluster and overnight on a high-end
desktop
A CutFEM method for two-phase flow problems
In this article, we present a cut finite element method for two-phase
Navier-Stokes flows. The main feature of the method is the formulation of a
unified continuous interior penalty stabilisation approach for, on the one
hand, stabilising advection and the pressure-velocity coupling and, on the
other hand, stabilising the cut region. The accuracy of the algorithm is
enhanced by the development of extended fictitious domains to guarantee a well
defined velocity from previous time steps in the current geometry. Finally, the
robustness of the moving-interface algorithm is further improved by the
introduction of a curvature smoothing technique that reduces spurious
velocities. The algorithm is shown to perform remarkably well for low capillary
number flows, and is a first step towards flexible and robust CutFEM algorithms
for the simulation of microfluidic devices
A partition of unity approach to fluid mechanics and fluid-structure interaction
For problems involving large deformations of thin structures, simulating
fluid-structure interaction (FSI) remains challenging largely due to the need
to balance computational feasibility, efficiency, and solution accuracy.
Overlapping domain techniques have been introduced as a way to combine the
fluid-solid mesh conformity, seen in moving-mesh methods, without the need for
mesh smoothing or re-meshing, which is a core characteristic of fixed mesh
approaches. In this work, we introduce a novel overlapping domain method based
on a partition of unity approach. Unified function spaces are defined as a
weighted sum of fields given on two overlapping meshes. The method is shown to
achieve optimal convergence rates and to be stable for steady-state Stokes,
Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for
FSI in the case of a 2D mock aortic valve simulation. These initial results
point to the potential applicability of the method to a wide range of FSI
applications, enabling boundary layer refinement and large deformations without
the need for re-meshing or user-defined stabilization.Comment: 34 pages, 15 figur
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
Introduction to discrete functional analysis techniques for the numerical study of diffusion equations with irregular data
We give an introduction to discrete functional analysis techniques for
stationary and transient diffusion equations. We show how these techniques are
used to establish the convergence of various numerical schemes without assuming
non-physical regularity on the data. For simplicity of exposure, we mostly
consider linear elliptic equations, and we briefly explain how these techniques
can be adapted and extended to non-linear time-dependent meaningful models
(Navier--Stokes equations, flows in porous media, etc.). These convergence
techniques rely on discrete Sobolev norms and the translation to the discrete
setting of functional analysis results
DOLFIN: Automated Finite Element Computing
We describe here a library aimed at automating the solution of partial differential equations using the finite element method. By employing novel techniques for automated code generation, the library combines a high level of expressiveness with efficient computation. Finite element variational forms may be expressed in near mathematical notation, from which low-level code is automatically generated, compiled and seamlessly integrated with efficient implementations of
computational meshes and high-performance linear algebra. Easy-to-use object-oriented interfaces to the library are provided in the form of a C++ library and a Python module. This paper discusses the mathematical abstractions and methods used in the design of the library and its implementation. A number of examples are presented to demonstrate the use of the library in application code
- ā¦