pinecone/core-2020-05-10-deduplication huggingface dataset:
A decomposition-based approach to uncertainty analysis of feed-forward multicomponent systems

To support effective decision making, engineers should comprehend and manage various uncertainties throughout the design process. Unfortunately, in today's modern systems, uncertainty analysis can become cumbersome and computationally intractable for one individual or group to manage. This is particularly true for systems comprised of a large number of components. In many cases, these components may be developed by different groups and even run on different computational platforms. This paper proposes an approach for decomposing the uncertainty analysis task among the various components comprising a feed-forward system and synthesizing the local uncertainty analyses into a system uncertainty analysis. Our proposed decomposition-based multicomponent uncertainty analysis approach is shown to be provably convergent in distribution under certain conditions. The proposed method is illustrated on quantification of uncertainty for a multidisciplinary gas turbine system and is compared to a traditional system-level Monte Carlo uncertainty analysis approach.SUTD-MIT International Design CentreUnited States. Defense Advanced Research Projects Agency. META Program (United States. Air Force Research Laboratory Contract FA8650-10-C-7083)Vanderbilt University (Contract VU-DSR#21807-S7)United States. Federal Aviation Administration. Office of Environment and Energy (FAA Award 09-C-NE-MIT, Amendments 028, 033, and 038
Real-time monitoring of thermal processes by reduced order modeling

This work presents a simple technique for real-time monitoring of thermal processes. Real-time simulationbased control of thermal processes is a big challenge because high-fidelity numerical simulations are costly and cannot be used, in general, for real-time decision making. Very often, processes are monitored or controlled with a few measurements at some specific points. Thus, the strategy presented here is centered on fast evaluation of the response only where it is needed. To accomplish this, classical harmonic analysis is combined with recent model reduction techniques. This leads to an advanced harmonic methodology, which solves in real-time the transient heat equation at the monitored point. In order to apply the reciprocity principle, harmonic analysis is used in the space-frequency domain. Then, Proper Generalized Decomposition, a reduced order approach, pre-computes a transfer function able to produce the output response for a given excitation. This transfer function is computed offline and only once. The response at the monitoring point can be recovered performing a computationally inexpensive post-processing step. This last step can be performed online for real-time monitoring of the thermal process. Examples show the applicability of this approach for a wide range of problems ranging from fast temperature evaluation to inverse problems.Peer ReviewedPostprint (author’s final draft
Generalized Robin-Neumann explicit coupling schemes for incompressible fluid-structure interaction: stability analysis and numerics

International audienceWe introduce a new class of explicit coupling schemes for the numerical solution of fluid-structure interaction problems involving a viscous incompressible fluid and an elastic structure. These methods generalize the arguments reported in [13, 10] to the case of the coupling with thick-walled structures. The basic idea lies in the derivation of an intrinsic interface Robin consistency at the space semi-discrete level, using a lumped-mass approximation in the structure. The fluid-solid splitting is then performed through appropriate extrapolations of the solid velocity and stress on the interface. Based on these methods, a new, parameter-free, Robin-Neumann iterative procedure is also proposed for the partitioned solution of implicit coupling. A priori energy estimates, guaranteeing the (added-mass free) stability of the schemes and the convergence of the iterative procedure, are established. The accuracy and robustness of the methods are illustrated in several numerical examples
Mixed enrichment for the finite element method in heterogeneous media.

Problems of multiple scales of interest or of locally nonsmooth solutions may often involve heterogeneous media. These problems are usually very demanding in terms of computations with the conventional finite element method. On the other hand, different enriched finite element methods such as the partition of unity, which proved to be very successful in treating similar problems, are developed and studied for homogeneous media. In this work, we present a new idea to extend the partition of unity finite element method to treat heterogeneous materials. The idea is studied in applications to wave scattering and heat transfer problems where significant advantages are noted over the standard finite element method. Although presented within the partition of unity context, the same enrichment idea can also be extended to other enriched methods to deal with heterogeneous materials
Adaptive $h$-refinement for reduced-order models

This work presents a method to adaptively refine reduced-order models \emph{a
posteriori} without requiring additional full-order-model solves. The technique
is analogous to mesh-adaptive $h$-refinement: it enriches the reduced-basis
space online by `splitting' a given basis vector into several vectors with
disjoint support. The splitting scheme is defined by a tree structure
constructed offline via recursive $k$-means clustering of the state variables
using snapshot data. The method identifies the vectors to split online using a
dual-weighted-residual approach that aims to reduce error in an output quantity
of interest. The resulting method generates a hierarchy of subspaces online
without requiring large-scale operations or full-order-model solves. Further,
it enables the reduced-order model to satisfy \emph{any prescribed error
tolerance} regardless of its original fidelity, as a completely refined
reduced-order model is mathematically equivalent to the original full-order
model. Experiments on a parameterized inviscid Burgers equation highlight the
ability of the method to capture phenomena (e.g., moving shocks) not contained
in the span of the original reduced basis.Comment: submitted to the International Journal for Numerical Methods in
  Engineering, Special Issue on Model Reductio
Weak imposition of essential boundary conditions in the finite element approximation of elliptic problems with non-matching meshes

This is the accepted version of the following article: [Codina, R., and Baiges, J. (2015) Weak imposition of essential boundary conditions in the finite element approximation of elliptic problems with non-matching meshes. Int. J. Numer. Meth. Engng, 104: 624–654. doi: 10.1002/nme.4815.], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.4815/abstractIn this work, we propose a method to prescribe essential boundary conditions in the finite element approximation of elliptic problems when the boundary of the computational domain does not match with the element boundaries. The problems considered are the Poisson problem, the Stokes problem, and the Darcy problem, the latter both in the primal and in the dual formulation. The formulation proposed is of variational type. The key idea is to start with the variational form that defines the problem and treat the boundary condition as a constraint. The particular feature is that the Lagrange multiplier is not defined on the boundary where the essential condition needs to be prescribed but is taken as a certain trace of a field defined in the computational domain, either in all of it or just in a region surrounding the boundary. When approximated numerically, this may allow one to condense the DOFs of the new field and end up with a problem posed only in terms of the original unknowns. The nature of the field used to weakly impose boundary conditions depends on the problem being treated. For the Poisson problem, it is a flux; for the Stokes problem, a stress; for the Darcy problem in primal form, a velocity field; and for the Darcy problem in dual form, it is a potential.In this work, we propose a method to prescribe essential boundary conditions in the finite element approximation of elliptic problems when the boundary of the computational domain does not match with the element boundaries. The problems considered are the Poisson problem, the Stokes problem, and the Darcy problem, the latter both in the primal and in the dual formulation. The formulation proposed is of variational type. The key idea is to start with the variational form that defines the problem and treat the boundary condition as a constraint. The particular feature is that the Lagrange multiplier is not defined on the boundary where the essential condition needs to be prescribed but is taken as a certain trace of a field defined in the computational domain, either in all of it or just in a region surrounding the boundary. When approximated numerically, this may allow one to condense the DOFs of the new field and end up with a problem posed only in terms of the original unknowns. The nature of the field used to weakly impose boundary conditions depends on the problem being treated. For the Poisson problem, it is a flux; for the Stokes problem, a stress; for the Darcy problem in primal form, a velocity field; and for the Darcy problem in dual form, it is a potential.Peer ReviewedPostprint (author's final draft
The effects of element shape on the critical time step in explicit time integrators for elasto-dynamics

In this paper, the effects of element shape on the critical time step are investigated. The common rule-of-thumb, used in practice, is that the critical time step is set by the shortest distance within an element divided by the dilatational (compressive) wave speed, with a modest safety factor. For regularly shaped elements, many analytical solutions for the critical time step are available, but this paper focusses on distorted element shapes. The main purpose is to verify whether element distortion adversely affects the critical time step or not. Two types of element distortion will be considered, namely aspect ratio distortion and angular distortion, and two particular elements will be studied: four-noded bilinear quadrilaterals and three-noded linear triangles. The maximum eigenfrequencies of the distorted elements are determined and compared to those of the corresponding undistorted elements. The critical time steps obtained from single element calculations are also compared to those from calculations based on finite element patches with multiple elements. Copyright (c) 2014 John Wiley & Sons, Ltd.Peer ReviewedPostprint (author’s final draft
Selective mass scaling for distorted solid-shell elements in explicit dynamics: optimal scaling factor and stable time step estimate

The use of solid-shell elements in explicit dynamics has been so far limited by the small critical time step resulting from the small thickness of these elements in comparison with the in-plane dimensions. To reduce the element highest eigenfrequency in inertia dominated problems, the selective mass scaling approach previously proposed in [G. Cocchetti, M. Pagani and U. Perego, Comp. \& Struct. 2013; 127:39-52.] for parallelepiped elements is here reformulated for distorted solid-shell elements. The two following objectives are achieved: the critical time step is governed by the smallest element in-plane dimension and not anymore by the thickness; the mass matrix remains diagonal after the selective mass scaling. The proposed approach makes reference to one Gauss point, trilinear brick element, for which the maximum eigenfrequency can be computed analytically. For this element, it is shown that the proposed mass scaling can be interpreted as a geometric thickness scaling, obtaining in this way a simple criterion for the definition of the optimal mass scaling factor. A strategy for the effective computation of the element maximum eigenfrequency is also proposed. The considered mass scaling preserves the element translational inertia, while it modifies the rotational one, leading to errors in the kinetic energy when the motion rotational component is dominant. The error has been rigorously assessed for an individual element, and a simple formula for its estimate has been derived. Numerical tests, both in small and large displacements and rotations, using a state-of-the-art solid-shell element taken from the literature, confirm the effectiveness and accuracy of the proposed approach. Copyright {\copyright} 2014 John Wiley \& Sons, Ltd
Dynamic finite element formulation for Cosserat elastic plates

International audienceA displacement and rotation-based dynamic finite element formulation for Cosserat plates is discussed in detail in this paper. Special attention is given to the validation of the element through adequate benchmarks and patch tests. The choice of the interpolation functions and the order of integration of the stiffness and the mass matrices are extensively argued. The possibility of local system deficiencies is explored, and a shear locking investigation specifically elaborated for Cosserat plates is carried out. It is shown how the present formulation has interesting computational properties as compared to a classical continuum-based formulation and how it can provide suitable results despite the use of reduced integration. An example of application of the finite element is given, in which the natural frequencies of a masonry panel modelled by means of discrete elements are computed and compared with the finite elements solution
On the effect of the bulk tangent matrix in partitioned solution schemes for nearly incompressible fluids

The purpose of this paper is to study the effect of the bulk modulus on the iterative solution of free surface quasi-incompressible fluids using a mixed partitioned scheme. A practical rule to set up the value of a pseudo-bulk modulus a priori in the tangent bulk stiffness matrix for improving the conditioning of the linear system of algebraic equations is also given. The efficiency of the proposed strategy is tested in several problems analyzing the advantage of the modified bulk tangent matrix with regard to the stability of the pressure field, the convergence rate and the computational speed of the analyses. The technique has been tested on a finite calculus/particle finite element method Lagrangian formulation, but it can be easily extended to other quasi-incompressible stabilized finite element formulations. Copyright (C) 2014 John Wiley & Sons, Ltd.Peer ReviewedPostprint (author's final draft
Hourglass stabilization and the virtual element method

In this paper, we establish the connections between the virtual element method (VEM) and the hourglass control techniques that have been developed since the early 1980s to stabilize underintegrated C0 Lagrange finite element methods. In the VEM, the bilinear form is decomposed into two parts: a consistent term that reproduces a given polynomial space and a correction term that provides stability. The essential ingredients of inline image-continuous VEMs on polygonal and polyhedral meshes are described, which reveals that the variational approach adopted in the VEM affords a generalized and robust means to stabilize underintegrated finite elements. We focus on the heat conduction (Poisson) equation and present a virtual element approach for the isoparametric four-node quadrilateral and eight-node hexahedral elements. In addition, we show quantitative comparisons of the consistency and stabilization matrices in the VEM with those in the hourglass control method of Belytschko and coworkers. Numerical examples in two and three dimensions are presented for different stabilization parameters, which reveals that the method satisfies the patch test and delivers optimal rates of convergence in the L2 norm and the H1 seminorm for Poisson problems on quadrilateral, hexahedral, and arbitrary polygonal meshes.A.C. gratefully acknowledges support from the College of Science and Engineering of the University of Leicester and the support of the EPSRC (grant EP/L022745/1). G.M. gratefully acknowledges the support of the LDRD-ER project #20140270 ‘From the finite element method to the virtual element method’ at the Los Alamos National Laboratory and of the DOE Office of Science Advanced Scientific Computing Research (ASCR) Program in Applied Mathematics Research. A.R. gratefully acknowledges the research support from the University of Milano–Bicocca. N.S. gratefully acknowledges the research support of the National Science Foundation through contract grant CMMI-1334783 to the University of California at Davis.Peer-reviewedPost-prin
A dual domain decomposition method for finite element digital image correlation

The computational burden associated to finite element based digital image correlation methods is mostly due to the inversion of finite element systems and to image interpolations. A non-overlapping dual domain decomposition method is here proposed to rationalise the computational cost of high resolution finite element digital image correlation measurements when dealing with large images. It consists in splitting the global mesh into submeshes and the reference and deformed states images into subset images. Classic finite element digital image correlation formulations are first written in each subdomain independently. The displacement continuity at the interfaces is enforced by introducing a set of Lagrange multipliers. The problem is then condensed on the interface and solved by a conjugate gradient algorithm. Three different preconditioners are proposed to accelerate its convergence. The proposed domain decomposition method is here exemplified with real high resolution images. It is shown to combine the metrological performances of finite element based digital image correlation and the parallelisation ability of subset based methods
The extended finite element method with new crack-tip enrichment functions for an interface crack between two dissimilar piezoelectric materials

This paper studies the static fracture problems of an interface crack in linear piezoelectric bimaterial by means of the extended finite element method (X-FEM) with new crack-tip enrichment functions. In the X-FEM, crack modeling is facilitated by adding a discontinuous function and crack-tip asymptotic functions to the classical finite element approximation within the framework of the partition of unity. In this work, the coupled effects of an elastic field and an electric field in piezoelectricity are considered. Corresponding to the two classes of singularities of the aforementioned interface crack problem, namely, E class and class, two classes of crack-tip enrichment functions are newly derived, and the former that exhibits oscillating feature at the crack tip is numerically investigated. Computation of the fracture parameter, i.e., the J-integral, using the domain form of the contour integral, is presented. Excellent accuracy of the proposed formulation is demonstrated on benchmark interface crack problems through comparisons with analytical solutions and numerical results obtained by the classical FEM. Moreover, it is shown that the geometrical enrichment combining the mesh with local refinement is substantially better in terms of accuracy and efficiency.postprin
Optimization of a regularized distortion measure to generate curved high-order unstructured tetrahedral meshes

We present a robust method for generating high-order nodal tetrahedral curved meshes. The approach consists of modifying an initial linear mesh by first, introducing high-order nodes, second, displacing the boundary nodes to ensure that they are on the CAD surface, and third, smoothing and untangling the mesh obtained after the displacement of the boundary nodes to produce a valid curved high-order mesh. The smoothing algorithm is based on the optimization of a regularized measure of the mesh distortion relative to the original linear mesh. This means that whenever possible, the resulting mesh preserves the geometrical features of the initial linear mesh such as shape, stretching and size. We present several examples to illustrate the performance of the proposed algorithm. Furthermore, the examples show that the implementation of the optimization problem is robust and capable of handling situations in which the mesh before optimization contains a large number of invalid elements. We consider cases with polynomial approximations up to degree ten, large deformations of the curved boundaries, concave boundaries, and highly stretched boundary layer elements. The meshes obtained are suitable for high-order finite element analyses.Peer ReviewedPostprint (author’s final draft
Proper generalized decomposition of a geometrically parametrized heat problem with geophysical applications

The solution of a steady thermal multiphase problem is assumed to be dependent on a set of parameters describing the geometry of the domain, the internal interfaces and the material properties. These parameters are considered as new independent variables. The problem is therefore stated in a multidimensional setup. The proper generalized decomposition (PGD) provides an approximation scheme especially well suited to preclude dramatically increasing the computational complexity with the number of dimensions. The PGD strategy is reviewed for the standard case dealing only with material parameters. Then, the ideas presented in [Ammar et al., Parametric solutions involving geometry: A step towards efficient shape optimization. Comput. Methods Appl. Mech. Eng., 2014; 268:178-193] to deal with parameters describing the domain geometry are adapted to a more general case including parametrization of the location of internal interfaces. Finally, the formulation is extended to combine the two types of parameters. The proposed strategy is used to solve a problem in applied geophysics studying the temperature field in a cross section of the Earth crust subsurface. The resulting problem is in a 10-dimensional space, but the PGD solution provides a fairly accurate approximation (error 1%) using less that 150 terms in the PGD expansion. Copyright (c) 2015John Wiley & Sons, Ltd.Peer ReviewedPostprint (author’s final draft
Locking-free discontinuous finite elements for the upper bound yield design of thick plates

International audienceThis work investigates the formulation of finite elements dedicated to the upper bound kinematic approach of yield design or limit analysis of Reissner-Mindlin thick plates in shear-bending interaction. The main novelty of this paper is to take full advantage of the fundamental difference between limit analysis and elasticity problems as regards the class of admissible virtual velocity fields. In particular, it has been demonstrated for 2D plane stress, plane strain or 3D limit analysis problems that the use of discontinuous velocity fields presents a lot of advantages when seeking for accurate upper bound estimates. For this reason, discontinuous interpolations of the transverse velocity and the rotation fields for Reissner-Mindlin plates are proposed. The subsequent discrete minimization problem is formulated as a second-order cone programming (SOCP) problem and is solved using the industrial software package Mosek. A comprehensive study of the shear-locking phenomenon is also conducted and it is shown that discontinuous elements avoid such a phenomenon quite naturally, whereas continuous elements cannot without any specific treatment. This particular aspect is confirmed through numerical examples on classical benchmark problems and the so-obtained upper bound estimates are confronted to recently developed lower bound equilibrium finite elements for thick plates