6,297 research outputs found
An unconditionally stable algorithm for generalized thermoelasticity based on operator-splitting and time-discontinuous Galerkin finite element methods
An efficient time-stepping algorithm is proposed based on operator-splitting and the space–time discontinuous Galerkin finite element method for problems in the non-classical theory of thermoelasticity. The non-classical theory incorporates three models: the classical theory based on Fourier’s law of heat conduction resulting in a hyperbolic–parabolic coupled system, a non-classical theory of a fully-hyperbolic extension, and a combination of the two. The general problem is split into two contractive sub-problems, namely the mechanical phase and the thermal phase. Each sub-problem is discretized using the space–time discontinuous Galerkin finite element method. The sub-problems are stable which then leads to unconditional stability of the global product algorithm. A number of numerical examples are presented to demonstrate the performance and capability of the method
von Neumann Stability Analysis of Globally Constraint-Preserving DGTD and PNPM Schemes for the Maxwell Equations using Multidimensional Riemann Solvers
The time-dependent equations of computational electrodynamics (CED) are
evolved consistent with the divergence constraints. As a result, there has been
a recent effort to design finite volume time domain (FVTD) and discontinuous
Galerkin time domain (DGTD) schemes that satisfy the same constraints and,
nevertheless, draw on recent advances in higher order Godunov methods. This
paper catalogues the first step in the design of globally constraint-preserving
DGTD schemes. The algorithms presented here are based on a novel DG-like method
that is applied to a Yee-type staggering of the electromagnetic field variables
in the faces of the mesh. The other two novel building blocks of the method
include constraint-preserving reconstruction of the electromagnetic fields and
multidimensional Riemann solvers; both of which have been developed in recent
years by the first author. We carry out a von Neumann stability analysis of the
entire suite of DGTD schemes for CED at orders of accuracy ranging from second
to fourth. A von Neumann stability analysis gives us the maximal CFL numbers
that can be sustained by the DGTD schemes presented here at all orders. It also
enables us to understand the wave propagation characteristics of the schemes in
various directions on a Cartesian mesh. We find that the CFL of DGTD schemes
decreases with increasing order. To counteract that, we also present
constraint-preserving PNPM schemes for CED. We find that the third and fourth
order constraint-preserving DGTD and P1PM schemes have some extremely
attractive properties when it comes to low-dispersion, low-dissipation
propagation of electromagnetic waves in multidimensions. Numerical accuracy
tests are also provided to support the von Neumann stability analysis
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
Numerical simulation of electrophoresis separation processes
A new Petrov-Galerkin finite element formulation has been proposed for transient convection-diffusion problems. Most Petrov-Galerkin formulations take into account the spatial discretization, and the weighting functions so developed give satisfactory solutions for steady state problems. Though these schemes can be used for transient problems, there is scope for improvement. The schemes proposed here, which consider temporal as well as spatial discretization, provide improved solutions. Electrophoresis, which involves the motion of charged entities under the influence of an applied electric field, is governed by equations similiar to those encountered in fluid flow problems, i.e., transient convection-diffusion equations. Test problems are solved in electrophoresis and fluid flow. The results obtained are satisfactory. It is also expected that these schemes, suitably adapted, will improve the numerical solutions of the compressible Euler and the Navier-Stokes equations
Symplectic Model Reduction of Hamiltonian Systems
In this paper, a symplectic model reduction technique, proper symplectic
decomposition (PSD) with symplectic Galerkin projection, is proposed to save
the computational cost for the simplification of large-scale Hamiltonian
systems while preserving the symplectic structure. As an analogy to the
classical proper orthogonal decomposition (POD)-Galerkin approach, PSD is
designed to build a symplectic subspace to fit empirical data, while the
symplectic Galerkin projection constructs a reduced Hamiltonian system on the
symplectic subspace. For practical use, we introduce three algorithms for PSD,
which are based upon: the cotangent lift, complex singular value decomposition,
and nonlinear programming. The proposed technique has been proven to preserve
system energy and stability. Moreover, PSD can be combined with the discrete
empirical interpolation method to reduce the computational cost for nonlinear
Hamiltonian systems. Owing to these properties, the proposed technique is
better suited than the classical POD-Galerkin approach for model reduction of
Hamiltonian systems, especially when long-time integration is required. The
stability, accuracy, and efficiency of the proposed technique are illustrated
through numerical simulations of linear and nonlinear wave equations.Comment: 25 pages, 13 figure
- …