Methods of mechanical contact constraint enforcement using the mortar finite element method
Explore the source record for details and available documents.
Engineering - Mechanical and civil engineering↗
SEARCH · Engineering Papers
Results for “Contact enforcement”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A simple levelset contact algorithm for large overlap removal and robust preloads
A simple approach to simulate contact between deformable objects is presented which relies on levelset descriptions of the Lagrangian geometry and an optimization-based solver. Modeling contact between objects remains a significant challenge for computational mechanics simulations. Common approaches are either plagued by lack of robustness or are exceedingly complex and require a significant number of heuristics. In contrast, the levelset contact approach presented herein is essentially heuristic free. Furthermore, the presented algorithm enables resolving and enforcing contact between objects with a significant amount of initial overlap. Examples demonstrating the feasibility of this approach are shown, including the standard Hertz contact problem, the robust removal of overlap between two overlapping blocks, and overlap-removal and pre-load for a bolted configuration.
42 ENGINEERING↗
Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method
Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.
58 GEOSCIENCES↗
On Practical Aspects of Variational Consistency in Contact Dynamics
Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.
42 ENGINEERING↗
Multiphasic droplet microfluidics platform for controlled bacteria and mammalian cell co-culture
Microfluidics has revolutionized high-throughput miniaturized biological assays. However, co-culture of mammalian cells and bacteria remains challenging in microfluidic systems due to incompatible growth requirements, limited spatial control, and the requirement for a mammalian cell adhesion matrix. Here, we present a microfluidic platform that generates multiphasic droplets which encapsulate mammalian and bacterial cells, enabling their direct and indirect co-culture. By combining photopolymerizable hydrogels with polymer phase separation, we generate core–shell droplets composed of a liquid and a hydrogel compartment. The hydrogel compartment supports mammalian cell adhesion and culture, while the liquid compartment sustains bacterial growth. We demonstrate two droplet architectures that allow physical bacteria–mammalian cell contacts or enforce complete physical separation, representing direct and indirect co-culture. Our multiphasic droplets are stable, customizable, able to sustain co-culture for over 24 hours, and compatible with fluorescence-based cell sorting technologies. Overall, our multiphasic droplet microfluidic platform provides a scalable and versatile tool for high-throughput co-culture and screening of host–microbe interactions.
59 BASIC BIOLOGICAL SCIENCES↗
BISON Robustness and Performance Improvements
BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Structure preserving and energy dissipative contact approaches for implicit dynamics
In this work, several structure preserving and energy dissipative contact approaches are proposed and evaluated. The time integration schemes considered are general with regard to the version of constraint type, but here the emphasis was on mortar contact. The proposed mortar contact approach conserves both linear and angular momentum for mortar contact in a novel way. The proposed time integration scheme can conserve energy or provide strict contact dissipation. In addition, the proposed scheme enforces both gap constraints and gap velocity constraints (i.e., persistency). Using a midstep time integrator often causes energy dissipation during initial impact, here this energy can be recovered in a novel way. Enforcing the gap velocity constraint mitigates the contact chatter of the contact pressure and nodes in many problems. Whereas some approaches enforce the persistency condition and gap constraints simultaneously during the solution of the equations of motion (EOM) requiring multipliers for both constraints included in the equation set, here the gap constraint is solved through the equations of motion and the persistency condition is satisfied in the time integration scheme by the velocity update after the equations of motion. It is shown that this approach is strictly dissipative in that a plastic contact condition can be achieved. Analogous to a coefficient of restitution for rigid bodies, any dissipated energy can then be returned upon release if energy conservation is desired. Structure preserving methods are good for long-time dynamics simulations and energy conserving and strictly dissipative methods can overcome stability issues associated with standard time integration algorithms such as the Newmark method.
42 ENGINEERING↗
Interface Problem Formulation Improvements with Application to Nuclear Fuel Performance Analysis
The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation Program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report documents new developments and robustness improvements in mechanical and thermal (gap heat transfer) contact formulations. The improvements range from the migration of industrial level (“assessment”) nuclear fuel model setups to the usage of mortar formulations, the addition of frictional contact to one-dimensional layered representations of fuel and cladding components, and the addition of the Petrov-Galerkin approach to dual mortar, which improves performance on curved, relatively coarse meshes. In addition, the Lagrange-multiplier enforcement of mechanical mortar contact constraints has been extended to two additional types of enforcement: penalty and augmented Lagrange-Uzawa. We show that the latter approach yields the same interface results as dual mortar in the Multiphysics Object-Oriented Simulation Environment, with the advantage of not worsening the condition number of the system matrix—thereby enabling the use of some general implementations of iterative preconditioners, at the expense of additional system evaluations (i.e., augmentations).
42 ENGINEERING↗
Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners
Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.
97 MATHEMATICS AND COMPUTING↗
Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro
Abstract In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like protease (3CLpro) can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.50 Å resolution crystal structure of 3CLpro C145S bound to NEMO 226–234 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for, in the pathology of COVID-19.
60 APPLIED LIFE SCIENCES↗
A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies
A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.
97 MATHEMATICS AND COMPUTING↗
The latent variable proximal point algorithm for variational problems with inequality constraints
The latent variable proximal point (LVPP) algorithm is a framework for solving infinite-dimensional variational problems with pointwise inequality constraints. The algorithm is a saddle point reformulation of the Bregman proximal point algorithm. At the continuous level, the two formulations are equivalent, but the saddle point formulation is more amenable to discretization because it introduces a structure-preserving transformation between a latent function space and the feasible set. Working in this latent space is much more convenient for enforcing inequality constraints than the feasible set, as discretizations can employ general linear combinations of suitable basis functions, and nonlinear solvers can involve general additive updates. LVPP yields numerical methods with observed mesh-independence for obstacle problems, contact, fracture, plasticity, and others besides; in many cases, for the first time. The framework also extends to more complex constraints, providing means to enforce convexity in the Monge–Ampère equation and handling quasi-variational inequalities, where the underlying constraint depends implicitly on the unknown solution. Here, in this paper, we describe the LVPP algorithm in a general form and apply it to ten problems from across mathematics.
Inequality constraints↗
Nonlinear convergence in contact mechanics: Immersed boundary finite volume
In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.
42 ENGINEERING↗
A mortar thermomechanical contact computational framework for nuclear fuel performance simulation
Nuclear fuel performance simulations involve the modeling of complex physical phenomena, ranging from fission gas release to fuel swelling and other temperature-induced effects. For light-water reactors (LWRs), swelling of the fuel and the pressure it imposes on the clad when they come into contact causes permanent clad deformation. Accurately characterizing the fuel-cladding interaction, which involves multiple physics, is essential to accurately simulate the fuel/cladding system. Thermomechanical modeling of this problem using a variationally consistent enforcement (e.g., a mortar approach) has been shown to improve the quality of results and facilitate convergence. Here, we present a general multiphysics computational framework for solving nuclear fuel problems using a mortar approach in BISON, a nuclear fuel performance code. In this study analyses show that using the mortar approach, which enables variationally consistent constraint enforcement, improves the quality of results as compared to the more commonly used node-on-face enforcement for representative LWR nuclear fuel simulations.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Non-conformal interface-cohesive modeling with the shifted boundary method
The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.
Yang, Cheng-Hau↗
Embedding space approach to Lorentzian CFT amplitudes and causal spherical functions
Conformal field theory in a Minkowski setting is discussed in an embedding space approach, paying special attention to causality constraints for four-point amplitudes. The physics of dilatation and Lorentz boost is emphasized in specifying the noncompact maximal Abelian subgroup of S O ( d , 2 ) . Reduction of a conformal field theory four-point amplitudes as functions of cross ratios is shown to be equivalent to enforcing H bi-invariance, i.e., F ( h g h ′ ) = F ( g ) , with g ∈ S O ( d , 2 ) and H an appropriate subgroup. Causality is imposed by introducing appropriate semigroups. Causal zonal spherical functions are constructed, making contact with Minkowski conformal blocks introduced previously. Published by the American Physical Society 2024
Agarwal, Pulkit (ORCID:0000000346581691)↗
Improved Fuel Particle and Compact Matrix Fission Product Release Prediction with Improved Mass Transfer Model
This study was motivated by the need to incorporate more physics-based mass transfer models in BISON in order to improve fuel particle and compact matrix fission product release prediction. This document details the integration of new modeling capabilities in BISON, including (1) development of sorption mass transfer models; improvement of thermal contact model by adding neon to the gas gap inventory, (3) implementation of PCEA and IG-110 graphite thermal models, (4) validation using Advanced Gas Reactor (AGR)-3/4 compact and capsule rings data. These new capabilities have been shown to enforce the desired physics with satisfactory accuracy. BISON’s predictions of fission product release of AGR-3/4 compacts compare favorably with PARFUME and the concentration profiles across capsule rings show reasonable trend to the experiment measurement.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Developing nanometer scale, atomically precise metallo-catalysts with Molecular Lego
We propose to create atomically precise, highly robust, nanometer scale macromolecules that organize biomimetic metal binding groups in designed three-dimensional pockets to accelerate specific polymerization reactions. These atomically precise catalysts will create atomically precise polymers with improved mechanical and environmental properties from renewable building blocks. We seek to create and study nanometer scale catalysts that mimic metallo-enzymes in their activity and selectivity because these larger catalysts will be able to make more extensive non-covalent contacts with the transition state of the growing polymer. These catalysts will have extended lifetimes because they will better protect the metals and higher activity because additional reactive groups can be brought in close contact with the active site to enhance reactivity. They will be far more tolerant to extremes of temperature and will not denature in non-aqueous solvents because they are held together by multiple, strong, covalent bonds. Our approach to catalysts could be used to develop catalysts to replace noble metal-based catalysts when targeting other reactions, because the geometry of the active site enforced by the Molecular Lego scaffolding will access new reactivity, mimicking how nature uses earth abundant metals for much of its catalysis. We will initially target Lewis acid-based catalysts that assemble aliphatic polyesters with improved stereoselectivity, tacticity and complex alternation of monomers. Aliphatic polyesters have received growing attention as attractive, environmentally benign, and sustainable alternatives to polymers developed from petroleum feedstocks. Aliphatic polyesters are environmentally friendly because they undergo facile hydrolytic degradation to benign products, and they have excellent properties and high biocompatibility. We will assemble these catalysts using unique “Spiroligomers” (aka Molecular Lego) chemistry developed in the Schafmeister group combined with the organometallic chemistry expertise of the Dobereiner group. Polymers are an attractive synthetic target of nanoscale atomically precise catalysts because their resulting structure, stereochemistry and tacticity has profound impact on their properties (crystalline vs amorphous) and provides a readout of the catalytic mechanism. Polymer synthesis also serves as an example of “atomically precise manufacturing” where atomically precise nanoscale catalysts construct polymers with atomic precision at increasing length scales.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗