Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “contact model”

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Towards a predictive, physics-based friction model for the dynamics of jointed structures

Bolted connections are ubiquitous in mechanical designs and pose a significant challenge to understanding and predicting the vibration response of assembled structures. The present paper develops a physics-based rough contact model of the frictional interactions within a joint. This model sums over the probable interactions of asperities – defined as locally maximum surface features – to determine the contact forces. Here, the tangential contact forces vary smoothly between sticking and slipping and allow the model to better capture the qualitative trends of experimental amplitude dependent frequency and damping than previous studies. Furthermore, the novel model is generalized to allow for arbitrarily varying normal pressure including potential separation to better represent the interfacial dynamics. This includes developing a new, computationally tractable approximation to the analytical Mindlin partial slip solution for tangential loading of contacting spheres. Further, the results highlight the importance of accurately characterizing the as-built topology of the interface, the plastic behavior of the contacting asperities, the relevant length scale of asperities, and the eccentricity of asperities. A predictive friction coefficient based on plasticity provides a poor match to experiments, so fitting the friction coefficient is also considered. Numerical results are compared to experiments on the Brake-Reuß Beam to assess the predictive potential of the models. While blind predictions over-predict the slip limit, the current model presents a significant improvement in physics-based modeling and highlights areas for ongoing research.

42 ENGINEERING↗
DOE OSTI2024

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗
DOE OSTI2021

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↗
DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗
DOE OSTI2024

Machine Learning Aided Modeling of Granular Materials: A Review

Artificial intelligence (AI) has become a buzzy word since Google’s AlphaGo beat a world champion in 2017. In the past five years, machine learning as a subset of the broader category of AI has obtained considerable attention in the research community of granular materials. This work offers a detailed review of the recent advances in machine learning-aided studies of granular materials from the particle-particle interaction at the grain level to the macroscopic simulations of granular flow. This work will start with the application of machine learning in the microscopic particle-particle interaction and associated contact models. Then, different neural networks for learning the constitutive behaviour of granular materials will be reviewed and compared. Finally, the macroscopic simulations of practical engineering or boundary value problems based on the combination of neural networks and numerical methods are discussed. We hope readers will have a clear idea of the development of machine learning-aided modelling of granular materials via this comprehensive review work.

42 ENGINEERING↗
DOE OSTI2021

Low‐Density 2D Superlattices Assembled via Directional DNA Bonding

Abstract It is critical to assemble nanoparticles (NPs) into superlattices with controlled symmetries and spacings on substrates for metamaterials applications, where such structural parameters dictate their properties. Here, we use DNA to assemble anisotropic NPs of three shapes—cubes, octahedra, and rhombic dodecahedra—on substrates and investigate their thermally induced reorganization into two‐dimensional (2D) crystalline films. We report two new low‐density 2D structures, including a honeycomb lattice based on octahedral NPs. The low‐density lattices favored here are not usually seen when particles are crystallized via other bottom‐up assembly techniques. Furthermore, we show that, consistent with the complementary contact model, a primary driving force for crystallization is the formation of directional, face‐to‐face DNA bonds between neighboring NPs and between NPs and the substrate. Our results can be used to deliberately prepare crystalline NP films with novel morphologies.

Miao, Ziyi↗
DOE OSTI2021

Low-Density 2D Superlattices Assembled via Directional DNA Bonding

It is critical to assemble nanoparticles (NPs) into superlattices with controlled symmetries and spacings on substrates for metamaterials applications, where such structural parameters dictate their properties. Here, we use DNA to assemble anisotropic NPs of three shapes-cubes, octahedra, and rhombic dodecahedra-on substrates and investigate their thermally induced reorganization into two-dimensional (2D) crystalline films. We report two new low-density 2D structures, including a honeycomb lattice based on octahedral NPs. The low-density lattices favored here are not usually seen when particles are crystallized via other bottom-up assembly techniques. Furthermore, we show that, consistent with the complementary contact model, a primary driving force for crystallization is the formation of directional, face-to-face DNA bonds between neighboring NPs and between NPs and the substrate. Finally, our results can be used to deliberately prepare crystalline NP films with novel morphologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Variational quantum and neural quantum states algorithms for the linear complementarity problem

Variational quantum algorithms (VQAs) are promising hybrid quantum-classical methods designed to leverage the computational advantages of quantum computing while mitigating the limitations of current noisy intermediate-scale quantum (NISQ) hardware. Although VQAs have been demonstrated as proofs of concept, their practical utility in solving real-world problems—and whether quantum-inspired classical algorithms can match their performance—remains an open question. We present a novel application of the variational quantum linear solver (VQLS) and its classical neural quantum states-based counterpart, the variational neural linear solver (VNLS), as key components within a minimum map Newton solver for a complementarity-based rigid-body contact model. We demonstrate using the VNLS that our solver accurately simulates the dynamics of rigid spherical bodies during collision events. These results suggest that quantum and quantum-inspired linear algebra algorithms can serve as viable alternatives to standard linear algebra solvers for modelling certain physical systems.

neural quantum states↗
DOE OSTI2024

DNA-mediated assembly of Au bipyramids into anisotropic light emitting kagome superlattices

Colloidal crystal engineering with DNA allows one to design diverse superlattices with tunable lattice symmetry, composition, and spacing. Most of these structures follow the complementary contact model, maximizing DNA hybridization on building blocks and producing relatively close-packed lattices. Here, low-symmetry kagome superlattices are assembled from DNA-modified gold bipyramids that can engage only in partial DNA surface matching. The bipyramid dimensions and DNA length can be engineered for two different superlattices with rhombohedral unit cells, including one composed of a periodic stacking of kagome lattices. Enabled by the partial facet alignment, the kagome lattices exhibit lattice distortion, bipyramid twisting, and planar chirality. When conjugated with Cy-5 dyes, the kagome lattices serve as cavities with high-density optical states and large Purcell factors along lateral directions, leading to strong dipole radiation along thezaxis and facet-dependent light emission. Such complex optical properties make these materials attractive for lasers, displays, and quantum sensing constructs.

Science & Technology - Other Topics↗
DOE OSTI2020

PLEXOS Inputs for Voisin et al. 2020 [Data set]

These data files provide MOSART-WM monthly flows at hydropower dam and HUC4 outlet locations. The flows are transformed into monthly generation constraints scenarios for hydro (generation potential) and thermal power plants (derated capacity) for input to the PLEXOS model. Contact: Sean Turner , Nathalie Voisin Authors: Nathalie Voisin, Ana Dyreson, Tao Fu, Matt O'Connell, Sean Turner, Tian Zhou, Jordan Macknick

13 HYDRO ENERGY↗
DOE OSTI2022

Analytical modeling and simulation of electrical contact resistance for elastic rough electrode surface contact including frictional temperature rise

An improved electrical contact resistance (ECR) model for elastic rough electrode contact is proposed, incorporating the effects of asperity interactions and temperature rise by frictional and joule heating. The analytical simulation results show that the ECR decreases steeply at the beginning of the contact between Al and Cu. However, it becomes stabilized after reaching a specific contact force. It is also found that the longer elapsed sliding contact time, the higher ECR due to the increase in electrical resistivity of electrode materials by the frictional temperature rise at the interface. The effects of surface roughness parameters on ECR are studied through the 3 2 full-factorial design-of-experiment analysis. Based on the two representative roughness parameters, i.e., root-mean-square (rms) roughness and asperity radius, their individual and coupled effects on the saturated ECR are examined. The saturated ECR increases with the rms roughness for a rough machined surface condition, but it is hardly affected by the asperity radius. On the other hand, the saturated ECR increases with both the rms roughness and the asperity radius under a smooth thin film surface condition.

42 ENGINEERING↗
DOE OSTI2023

An approach to grid-to-rod fretting wear modeling using dynamic mortar contact

This article presents the development of a grid-to-rod fretting (GTRF) modeling methodology that utilizes dynamic mortar contact. Here, we leverage a recently developed computational framework for modeling thermomechanical contact in dynamic simulations from the nuclear fuel performance code Bison. Usage of mortar contact ensures a smooth distribution of normal and frictional contact forces, displacements, and velocities on the contact interface, thus facilitating an accurate computation of wear. The integration within an advanced nuclear fuel performance code enables analysis of complex interactions between the fuel, cladding, and spacer grid. Such interactions include swelling, creep, fuel fragmentation, burnup, plenum pressure, and fission gas release. Our methodology is demonstrated via the modeling of two spacer grid geometries and the influence of fuel-cladding mechanical contact on the generation of fretting wear.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Design optimization of semiconductor piezoresistors with Schottky diode contacts

A modeling theory is developed to predict the performance of piezoresistors which incorporate Schottky diode electrical contacts. This new theory allows the design of high performance gauges which can be fabricated using NonLithographically-Based Microfabrication (NLBM) techniques. These semiconductor piezoresistors can be designed in customizable sizes and fabricated in parallel in order to integrate position sensing into MEMS flexural positioners. Customizable nanopositioning platforms will enable advances in a range of nano-scale fabrication and metrology applications. A semiconductor piezoresistor with Schottky diode contacts was fabricated and attached to a titanium flexure. This device is shown to match predicted electrical performance within about 8% and to show a gauge factor of 116, within 2% of the predicted value. Performance limits for Schottky diode semiconductor piezoresistors are identified to be about 127dB full noise dynamic range for a quarter bridge over a 10kHz sensor bandwidth on a 600µm width titanium flexure, making them ideal for sensing on meso-/micro-scale flexural positioners. Methods of reaching these performance bound are suggested and their impact on the sensor dynamic range are studied.

42 ENGINEERING↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

An Anisotropic Yield and Damage Material Model to Improve the Contact Pressure Analysis in a Biomass Shredding System

Size reduction systems used in biomass processing break biomass into smaller pieces by utilizing the kinetic energy from the sharp rotating blades. Abrasive and/or erosive wear caused by biomass comminution results in blade wear of the sharp edged cutters, deteriorating the process efficiency. Here, this study aims to optimize the blade design and improve the system efficiency by attempting to understand the interactions between the blades and biomass particles. Since real-time monitoring of these interactions is impractical during operation, mechanical simulations offer a viable alternative for investigating the shredding process. Yet, the irregular geometry and complex mechanical properties of biomass—such as the anisotropic nature of woodchips and their nonlinear fracture behavior—pose significant challenges for accurately simulating contact pressure. In this work an anisotropic yield material model, along with a damage initiation and evolution function, is applied to the woodchip particle to study the contact pressure on shredder blade, offering a scientific basis for improved blade design and process efficiency. This approach can be extended to other biomass processing systems with similar anisotropic feedstocks, making it a valuable tool for advancing sustainable biomass utilization.

09 - BIOMASS FUELS↗
DOE OSTI2024

Efficient Scalable Contact Network Generation from Population Data

Modeling the contacts among a population is critical to understanding the dynamics of a disease outbreak. Contact networks, where nodes are individuals and edges are contacts among them, are used to represent these complex individual-level interactions. In this work, we are given the daily activity schedules of an urban population that represent the activity location and time of individuals in a population during a single twenty four hour period over multiple days. Using collocation to determine contact between individuals, our goal is to extract hourly contact networks from large-scale activity data. We improve upon the existing adjacency matrix-based method by implementing our custom sparse matrix multiplication algorithm. Starting with a Python implementation, we achieve a 1600x speed up in the computation with a fast custom designed sparse matrix multiplier algorithm implemented in the C++ language. This work is central to future parallel designs of the problem.

97 MATHEMATICS AND COMPUTING↗