Modeling Sea Ice Fracture in a Discrete Element Framework with a Cohesive Contact Model.
Abstract not provided.
O'Connor, Devin↗
SEARCH · Search NASA
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.
Flow Simulations of Size Disperse Granular Particles with Realistic Contact Models.
Abstract not provided.
Santos, Andrew↗
Transient pattern formation in an active matter contact poisoning model
Abstract One of the most notable features in repulsive particle based active matter systems is motility-induced-phase separation (MIPS) where a dense, often crystalline phase and low density fluid coexist. Most active matter studies involve time-dependent activity; however, there are many active systems where individual particles transition from living or moving to dead or nonmotile due to lack of fuel, infection, or poisoning. Here we consider an active matter particle system at densities where MIPS does not occur. When we add a small number of infected particles that can poison other particles, rendering them nonmotile, we find a rich variety of time dependent pattern formation, including MIPS, a wetting phase, and a fragmented state formed when mobile particles plow through a nonmotile packing. We map the patterns as a function of time scaled by epidemic duration, and show that the pattern formation is robust for a wide range of poisoning rates and activity levels. We also show that pattern formation does not occur in a random death model, but requires the promotion of nucleation by contact poisoning. Our results should be relevant to biological and active matter systems where there is some form of poisoning, death, or transition to nonmotility.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Sensitivity analysis of thermal contact conductance modeling to inform MiniFuel irradiation capsule designs
The MiniFuel irradiation platform has been developed by Oak Ridge National Laboratory as a flexible, high-throughput separate effects testing capability within the High Flux Isotope Reactor (HFIR). Finite element thermal models are relied upon to design MiniFuel experiments to achieve a specific time-averaged irradiation temperature for experimental objectives. A previous study identified that uncertainty in the component heat generation rates and thermal contact conductance (TCC) model are the most significant contributors to predicted fuel temperature variance. To address both sources of uncertainty, this work performs sensitivity analysis on the TCC model to identify high-impact, high-uncertainty parameters that contribute to fuel temperature variance. The TCC model is analyzed in increasing detail, first using a standalone Python code, then again after coupling Python to the BISON fuel performance code. Furthermore, the parameters with the largest contributions to fuel temperature variance which can be reduced through design changes are identified as the initial subcapsule gas pressure, contact pressure between the fuel and dish, and the effective surface roughness of the interface. A set of design recommendations for future capsule designs has been established and applied to reduce the previously quantified average fuel temperature uncertainty ranges of ± 40 °C in the HFIR vertical experiment facilities (VXF) and ± 80 °C in the removable beryllium (RB) reflector to approximately ± 32 °C and ± 53 °C, respectively. This equates to a 21 % and 33 % reduction in the uncertainty range of the average fuel temperature for VXF and RB, respectively.
BISON↗
Implementation of low-burnup swelling focused assessment case and enhancement of FIPD-BISON integration framework
This report summarizes the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification (V&V) project in FY2023, with a focus on the capability development and assessment case development for steady state irradiation experiments. The X423 IFR experiment assessment case, which is focused on the evaluation of low burnup swelling behavior of both binary and ternary metallic fuel, has been developed. It covers 58 fuel pins with various fuel compositions and 158 sets of post-irradiation examination (PIE) data with burnup ranging from ~0.5 at.% to ~5.0 at.%, providing a powerful tool within BISON for the calibration, evaluation, verification and validation of both existing and future metallic fuel swelling models at low burnup. Meanwhile, the X447 IFR experiment assessment case established in FY2022 was enhanced to utilize the latest advanced features in BISON, such as automatic differentiation and Mortar contact models. With such enhancements, especially the Mortar mechanical contact model with friction, BISON is now capable of predicting axial fuel growth that is consistent with experimental measurements. Additionally, during the development, gaps and potential scopes were identified and explorative efforts were made, providing insight for future works.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Dynamical responses of constrained pipe conveying fluids and its dependence on the modeling of the contact force
Accurately modeling the impact force used in the analysis of loosely constrained cantilevered pipes conveying fluid is imperative. If little information is known of the motion-limiting constraints used in experiments, the analysis of the system may yield inaccurate predictions. Here in this work, multiple forcing representations of the impact force are defined and analyzed for a cantilevered pipe that conveys fluid. Depending on the representation of the impact force, the dynamics of the pipe can vary greatly when only the stiffness of the constraints is known from experiments. Three gap sizes of the constraints are analyzed, and the representation of the impact force used to analyze the system is found to significantly affect the response of the pipe at each gap size. An investigation on the effects of the vibro-impact force representation is performed through using basin of attraction analysis and nonlinear characterization of the system’s response.
42 ENGINEERING↗
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↗
A high-resolution pseudo-polygon discrete element model for regional sea ice
Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.
54 ENVIRONMENTAL SCIENCES↗
Statistical Analysis of Intertube Tunneling Contacts in the Macroscopic Electrical Conductivity of Carbon Nanotube Fibers
Here, this study investigates the influence of tunneling contact resistances between carbon nanotubes (CNTs) on electron transport and electrical conductivity of macroscopic carbon nanofibers (CNFs), which profoundly impacts the performance of CNT thin film electronics, CNF electron emitters and cathodes, and energy conversion and storage devices. Utilizing a self-consistent electrical contact model coupling a transmission line model with tunneling current, we calculate the contact resistances of a plethora of CNT-CNT contacts within a CNF fiber, which consists of aligned, densely packed CNTs. A statistical analysis is conducted, using Gaussian distributions to account for variations in contact lengths, tunneling gap distances, and single CNT aspect ratios, to calculate the CNT-CNT contact resistance and the overall resistance of CNT fiber. By scaling our model to a macroscopic level, our results are in good agreement with experimental measurements. Our calculation suggests that while increasing the contact overlap length diminishes individual CNT-CNT contact resistance, it could paradoxically increase macroscopic CNT fiber resistance for a given constant CNF mass density, which is due to that fact that a larger overlap length allows more CNTs to pack along an electrical conduction path per unit length, leading to more tunneling contact junctions connected in series and thus less number of parallel conduction paths within the fiber cross section. Increasing tunneling gap distance increases both individual contact and overall fiber resistance. This research provides a simple design tool for tailoring CNT fiber electrical properties to promote real-world applications using CNTs or similar low-dimensional materials.
36 MATERIALS SCIENCE↗
Computational design of metamaterials with self contact
Inverse homogenization in combination with contact modeling, topology optimization and shape optimization is used to design metamaterials with optimized macroscopic response. The homogenization assumes length scale separation which allows the non-linear macroscopic behavior to be obtained by analyzing a single unit cell in a lattice structure. Self contact in the unit cell, which is modeled using a third medium contact method, is leveraged to obtain a complex homogenized response. The inverse homogenization problem is initially formulated as a topology optimization problem, where the macroscopic stress–strain behavior is tuned to our liking. However, it is well known that boundary phenomena are difficult to model in topology optimization and that interface modeling is crucial to accurately analyze contact. For that reason, the boundary representation of the topology optimized design is extracted and used as initial design in a subsequent shape optimization. The behaviors of our designs are verified by performing rigorous post-processing analyzes using conforming meshes and conventional contact formulations.
42 ENGINEERING↗
Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes
Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.
36 MATERIALS SCIENCE↗
A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks
We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.
bioenergy↗
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↗
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↗
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↗
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↗
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↗