16,917 research outputs found
Efficient computation of high index Sturm-Liouville eigenvalues for problems in physics
Finding the eigenvalues of a Sturm-Liouville problem can be a computationally
challenging task, especially when a large set of eigenvalues is computed, or
just when particularly large eigenvalues are sought. This is a consequence of
the highly oscillatory behaviour of the solutions corresponding to high
eigenvalues, which forces a naive integrator to take increasingly smaller
steps. We will discuss some techniques that yield uniform approximation over
the whole eigenvalue spectrum and can take large steps even for high
eigenvalues. In particular, we will focus on methods based on coefficient
approximation which replace the coefficient functions of the Sturm-Liouville
problem by simpler approximations and then solve the approximating problem. The
use of (modified) Magnus or Neumann integrators allows to extend the
coefficient approximation idea to higher order methods
Randomized Local Model Order Reduction
In this paper we propose local approximation spaces for localized model order
reduction procedures such as domain decomposition and multiscale methods. Those
spaces are constructed from local solutions of the partial differential
equation (PDE) with random boundary conditions, yield an approximation that
converges provably at a nearly optimal rate, and can be generated at close to
optimal computational complexity. In many localized model order reduction
approaches like the generalized finite element method, static condensation
procedures, and the multiscale finite element method local approximation spaces
can be constructed by approximating the range of a suitably defined transfer
operator that acts on the space of local solutions of the PDE. Optimal local
approximation spaces that yield in general an exponentially convergent
approximation are given by the left singular vectors of this transfer operator
[I. Babu\v{s}ka and R. Lipton 2011, K. Smetana and A. T. Patera 2016]. However,
the direct calculation of these singular vectors is computationally very
expensive. In this paper, we propose an adaptive randomized algorithm based on
methods from randomized linear algebra [N. Halko et al. 2011], which constructs
a local reduced space approximating the range of the transfer operator and thus
the optimal local approximation spaces. The adaptive algorithm relies on a
probabilistic a posteriori error estimator for which we prove that it is both
efficient and reliable with high probability. Several numerical experiments
confirm the theoretical findings.Comment: 31 pages, 14 figures, 1 table, 1 algorith
Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches
A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed
High performance reduced order modeling techniques based on optimal energy quadrature: application to geometrically non-linear multiscale inelastic material modeling
A High-Performance Reduced-Order Model (HPROM) technique, previously presented by the authors in the context of hierarchical multiscale models for non linear-materials undergoing infinitesimal strains, is generalized to deal with large deformation elasto-plastic problems. The proposed HPROM technique uses a Proper Orthogonal Decomposition procedure to build a reduced basis of the primary kinematical variable of the micro-scale problem, defined in terms of the micro-deformation gradient fluctuations. Then a Galerkin-projection, onto this reduced basis, is utilized to reduce the dimensionality of the micro-force balance equation, the stress homogenization equation and the effective macro-constitutive tangent tensor equation. Finally, a reduced goal-oriented quadrature rule is introduced to compute the non-affine terms of these equations. Main importance in this paper is given to the numerical assessment of the developed HPROM technique. The numerical experiments are performed on a micro-cell simulating a randomly distributed set of elastic inclusions embedded into an elasto-plastic matrix. This micro-structure is representative of a typical ductile metallic alloy. The HPROM technique applied to this type of problem displays high computational speed-ups, increasing with the complexity of the finite element model. From these results, we conclude that the proposed HPROM technique is an effective computational tool for modeling, with very large speed-ups and acceptable accuracy levels with respect to the high-fidelity case, the multiscale behavior of heterogeneous materials subjected to large deformations involving two well-separated scales of length.Peer ReviewedPostprint (author's final draft
A 4-node assumed-stress hybrid shell element with rotational degrees of freedom
An assumed-stress hybrid/mixed 4-node quadrilateral shell element is introduced that alleviates most of the deficiencies associated with such elements. The formulation of the element is based on the assumed-stress hybrid/mixed method using the Hellinger-Reissner variational principle. The membrane part of the element has 12 degrees of freedom including rotational or drilling degrees of freedom at the nodes. The bending part of the element also has 12 degrees of freedom. The bending part of the element uses the Reissner-Mindlin plate theory which takes into account the transverse shear contributions. The element formulation is derived from an 8-node isoparametric element. This process is accomplished by assuming quadratic variations for both in-plane and out-of-plane displacement fields and linear variations for both in-plane and out-of-plane rotation fields along the edges of the element. In addition, the degrees of freedom at midside nodes are approximated in terms of the degrees of freedom at corner nodes. During this process the rotational degrees of freedom at the corner nodes enter into the formulation of the element. The stress field are expressed in the element natural-coordinate system such that the element remains invariant with respect to node numbering
Adaptive meshless centres and RBF stencils for Poisson equation
We consider adaptive meshless discretisation of the Dirichlet problem for Poisson equation based on numerical differentiation stencils obtained with the help of radial basis functions. New meshless stencil selection and adaptive refinement algorithms are proposed in 2D. Numerical experiments show that the accuracy of the solution is comparable with, and often better than that achieved by the mesh-based adaptive finite element method
An algebraic multigrid method for mixed discretizations of the Navier-Stokes equations
Algebraic multigrid (AMG) preconditioners are considered for discretized
systems of partial differential equations (PDEs) where unknowns associated with
different physical quantities are not necessarily co-located at mesh points.
Specifically, we investigate a mixed finite element discretization of
the incompressible Navier-Stokes equations where the number of velocity nodes
is much greater than the number of pressure nodes. Consequently, some velocity
degrees-of-freedom (dofs) are defined at spatial locations where there are no
corresponding pressure dofs. Thus, AMG approaches leveraging this co-located
structure are not applicable. This paper instead proposes an automatic AMG
coarsening that mimics certain pressure/velocity dof relationships of the
discretization. The main idea is to first automatically define coarse
pressures in a somewhat standard AMG fashion and then to carefully (but
automatically) choose coarse velocity unknowns so that the spatial location
relationship between pressure and velocity dofs resembles that on the finest
grid. To define coefficients within the inter-grid transfers, an energy
minimization AMG (EMIN-AMG) is utilized. EMIN-AMG is not tied to specific
coarsening schemes and grid transfer sparsity patterns, and so it is applicable
to the proposed coarsening. Numerical results highlighting solver performance
are given on Stokes and incompressible Navier-Stokes problems.Comment: Submitted to a journa
- ā¦