We develop a mixed formulation for incompressible hyperelastodynamics based on a continuum modeling framework recently developed in the work of Liu and Marsden and smooth generalizations of the Taylor‐Hood element based on nonuniform rational B‐splines (NURBSs). This continuum formulation draws a link between computational fluid dynamics and computational solid dynamics. This link inspires an energy stability estimate for the spatial discretization, which favorably distinguishes the formulation from the conventional mixed formulations for finite elasticity. The inf‐sup condition is utilized to provide a bound for the pressure field. The generalized‐
We present a general framework to solve elastodynamic problems by means of the virtual element method (VEM) with explicit time integration. In particular, the VEM is extended to analyze nearly incompressible solids using the B‐bar method. We show that, to establish a B‐bar formulation in the VEM setting, one simply needs to modify the stability term to stabilize only the deviatoric part of the stiffness matrix, which requires no additional computational effort. Convergence of the numerical solution is addressed in relation to stability, mass lumping scheme, element size, and distortion of arbitrary elements, either convex or nonconvex. For the estimation of the critical time step, two approaches are presented, ie, the maximum eigenvalue of a system of mass and stiffness matrices and an effective element length. Computational results demonstrate that small edges on convex polygonal elements do not significantly affect the critical time step, whereas convergence of the VEM solution is observed regardless of the stability term and the element shape in both two and three dimensions. This extensive investigation provides numerical recipes for elastodynamic VEMs with explicit time integration and related problems.
more » « less- NSF-PAR ID:
- 10460569
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Engineering
- Volume:
- 121
- Issue:
- 1
- ISSN:
- 0029-5981
- Page Range / eLocation ID:
- p. 1-31
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Summary α method is applied for temporal discretization, and a nested block preconditioner is invoked for the solution procedure. The inf‐sup stability for different pairs of NURBS elements is elucidated through numerical assessment. The convergence rate of the proposed formulation with various combinations of mixed elements is examined by the manufactured solution method. The numerical scheme is also examined under compressive and tensile loads for isotropic and anisotropic hyperelastic materials. Finally, a suite of dynamic problems is numerically studied to corroborate the stability and conservation properties. -
Although the full form of the Rayleigh–Plesset (RP) equation more accurately depicts the bubble behavior in a cavitating flow than its reduced form, it finds much less application than the latter in the computational fluid dynamic (CFD) simulation due to its high stiffness. The traditional variable time-step scheme for the full form RP equation is difficult to be integrated with the CFD program since it requires a tiny time step at the singularity point for convergence and this step size may be incompatible with time marching of conservation equations. This paper presents two stable and efficient numerical solution schemes based on the finite difference method and Euler method so that the full-form RP equation can be better accepted by the CFD program. By employing a truncation bubble radius to approximate the minimum bubble size in the collapse stage, the proposed schemes solve for the bubble radius and wall velocity in an explicit way. The proposed solution schemes are more robust for a wide range of ambient pressure profiles than the traditional schemes and avoid excessive refinement on the time step at the singularity point. Since the proposed solution scheme can calculate the effects of the second-order term, liquid viscosity, and surface tension on the bubble evolution, it provides a more accurate estimation of the wall velocity for the vaporization or condensation rate, which is widely used in the cavitation model in the CFD simulation. The legitimacy of the solution schemes is manifested by the agreement between the results from these schemes and established ones from the literature. The proposed solution schemes are more robust in face of a wide range of ambient pressure profiles.more » « less
-
more » « less
Abstract As the use of spectral/
hp element methods, and high-order finite element methods in general, continues to spread, community efforts to create efficient, optimized algorithms associated with fundamental high-order operations have grown. Core tasks such as solution expansion evaluation at quadrature points, stiffness and mass matrix generation, and matrix assembly have received tremendous attention. With the expansion of the types of problems to which high-order methods are applied, and correspondingly the growth in types of numerical tasks accomplished through high-order methods, the number and types of these core operations broaden. This work focuses on solution expansion evaluation at arbitrary points within an element. This operation is core to many postprocessing applications such as evaluation of streamlines and pathlines, as well as to field projection techniques such as mortaring. We expand barycentric interpolation techniques developed on an interval to 2D (triangles and quadrilaterals) and 3D (tetrahedra, prisms, pyramids, and hexahedra) spectral/hp element methods. We provide efficient algorithms for their implementations, and demonstrate their effectiveness using the spectral/hp element libraryNektar++ by running a series of baseline evaluations against the ‘standard’ Lagrangian method, where an interpolation matrix is generated and matrix-multiplication applied to evaluate a point at a given location. We present results from a rigorous series of benchmarking tests for a variety of element shapes, polynomial orders and dimensions. We show that when the point of interest is to be repeatedly evaluated, the barycentric method performs at worst$$50\%$$ $$7\times $$ $$30\%$$ $$6\times $$ $${\mathcal {O}}(k)$$ $${\mathcal {O}}(k^2)$$ -
We present a virtual element method (VEM)-based topology optimization framework using polyhedral elements, which allows for convenient handling of non-Cartesian design domains in three dimensions. We take full advantage of the VEM properties by creating a unified approach in which the VEM is employed in both the structural and the optimization phases. In the structural problem, the VEM is adopted to solve the three-dimensional elasticity equation. Compared to the finite element method, the VEM does not require numerical integration (when linear elements are used) and is less sensitive to degenerated elements (e.g., ones with skinny faces or small edges). In the optimization problem, we introduce a continuous approximation of material densities using the VEM basis functions. When compared to the standard element-wise constant approximation, the continuous approximation enriches the geometrical representation of structural topologies. Through two numerical examples with exact solutions, we verify the convergence and accuracy of both the VEM approximations of the displacement and material density fields. We also present several design examples involving non-Cartesian domains, demonstrating the main features of the proposed VEM-based topology optimization framework. The source code for a MATLAB implementation of the proposed work, named PolyTop3D, is available in the (electronic) Supplementary Material accompanying this publication.more » « less
-
null (Ed.)Abstract Combining the classical theory of optimal transport with modern operator splitting techniques, we develop a new numerical method for nonlinear, nonlocal partial differential equations, arising in models of porous media, materials science, and biological swarming. Our method proceeds as follows: first, we discretize in time, either via the classical JKO scheme or via a novel Crank–Nicolson-type method we introduce. Next, we use the Benamou–Brenier dynamical characterization of the Wasserstein distance to reduce computing the solution of the discrete time equations to solving fully discrete minimization problems, with strictly convex objective functions and linear constraints. Third, we compute the minimizers by applying a recently introduced, provably convergent primal dual splitting scheme for three operators (Yan in J Sci Comput 1–20, 2018). By leveraging the PDEs’ underlying variational structure, our method overcomes stability issues present in previous numerical work built on explicit time discretizations, which suffer due to the equations’ strong nonlinearities and degeneracies. Our method is also naturally positivity and mass preserving and, in the case of the JKO scheme, energy decreasing. We prove that minimizers of the fully discrete problem converge to minimizers of the spatially continuous, discrete time problem as the spatial discretization is refined. We conclude with simulations of nonlinear PDEs and Wasserstein geodesics in one and two dimensions that illustrate the key properties of our approach, including higher-order convergence our novel Crank–Nicolson-type method, when compared to the classical JKO method.more » « less
