Discrete Element Modeling of Sea Ice Dynamics and Failure with a Damage-Based Cohesive Contact Law
Explore the source record for details and available documents.
O'Connor, Devin Thomas↗
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.
Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets
Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.
97 MATHEMATICS AND COMPUTING↗
Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations
Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.
97 MATHEMATICS AND COMPUTING↗
Thermoelectric coolers for high-power-density 3D electronics heat management
Future advancements in three-dimensional (3D) electronics require robust thermal management methodology. Thermoelectric coolers (TECs) are reliable and solid-state heat pumping devices with high cooling capacity that can meet the requirements of emerging 3D microelectronic devices. In this work, we first provide the design of TECs for electronics cooling using a computational model and then experimentally validate the main predictions. Key device parameters such as device thickness, leg density, and contact resistance were studied to understand their influence on the performance of TECs. Our results show that it is possible to achieve high cooling power density through optimization of TE leg height and packing density. Scaling of TECs is shown to provide ultra-high cooling power density.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
A particle finite element-based model for droplet spreading analysis
Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.
42 ENGINEERING↗
Cladding Profilometry Analysis of Experimental Breeder Reactor-II Metallic Fuel Pins with HT9, D9, and SS316 Cladding
BISON finite element method fuel performance simulations were conducted using an existing automated process that couples the Fuels Irradiation & Physics Database (FIPD) and the Integral Fast Reactor Materials Information System database by writing input files and comparing the BISON output to post-irradiation fuel pin profilometry measurements contained within the databases. The importance of this work is to demonstrate the ability to benchmark fuel performance metallic fuel models within BISON using Experimental Breeder Reactor-II fuel pin data for a number of similar pins, while building off previous modeling efforts. Changes to the generic BISON input file include implementing pin specific axial power and flux profiles, pin specific fluences, frictional contact, and irradiation-induced volumetric swelling models for cladding. A statistical analysis of irradiation-induced volumetric swelling models for HT9, D9, and SS316 was performed for experiments X421/X421A, X441/X441A, and X486. Between these three experiments, there were 174 post-irradiation examination (PIE) profilometries used for validating the swelling models presented using a standard error of the estimate (SEE) method. Implementation of the volumetric swelling models for D9 and SS316 claddings was found to have a significant impact on the BISON profilometry simulated, where HT9 clad pins had an insignificant change due to low fluence values. BISON profilometry simulated for HT9, D9, and SS316 fuel pins agreed with PIE profilometry measurements, with assembly SEE values being 4.4 × 10−3 for X421A, 2.0 × 10−3 for X441A, and 2.8 × 10−3 for X486. D9 clad pins in X421/X421A had the highest SEE values, which is due to the BISON simulated profilometry being shifted axially. While this work accomplished its purpose to demonstrate the modeling of multiple fuel pins from the databases to help validate models, the results suggest that the continued development of metallic fuel models is necessary for qualifying new metallic fuel systems to better capture some physical performance phenomena, such as the hot pressing of U-Pu-Zr and the fuel cladding chemical interaction.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Spherical and sessile droplet dynamics by fluctuating hydrodynamics
We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.
Bell, John B↗
An Analytical Solution for the Initiation and Early Progression of Fretting Wear in Spherical Contacts
Abstract This article derives analytical solutions to calculate the wear volume at the initiation of fretting motion and its early progression over the first few oscillation cycles. The Archard-based model considers a deformable hemisphere that is contact with a deformable flat block. The material pairs investigated are special alloys, the Inconel 617/Incoloy 800H, and Inconel 617/Inconel 617. The analytical study begins with a unidirectional frictional sliding contact, where the local interfacial sliding distance and the nominal sliding distance at the initiation of gross slip are derived. The obtained analytical expressions for unidirectional sliding are then used to derive the corresponding wear volume for the initiation and early progression of gross slip and the wear volume for a general fretting cycle under elastic conditions. These analytical derivations are all verified by the finite element analysis (FEA). The FEA method and the analytical solutions render virtually identical results for both similar and dissimilar material pairs. The effects of plasticity on the wear volume under elastic–plastic conditions are also investigated. It is found that the fretting wear volumes obtained from the FEA simulations, which include plasticity, are close to those obtained from the analytical expressions for purely elastic regimes. All the results are presented in normalized forms, which can be easily generalized and applied to three-dimensional fretting wear of other material pairs.
Engineering↗
Phase transitions in Schloegl's second model for autocatalysis on a Bethe lattice
Schloegl's second model (also known as the quadratic contact process) on a lattice involves spontaneous particle annihilation at rate p and autocatalytic particle creation at empty sites with n ≥ 2 occupied neighbors. The particle creation rate for exactly n occupied neighbors is selected here as n(n - 1)/[z(z - 1)] for lattice coordination number z. We analyze this model on a Bethe lattice. Precise behavior for stochastic models on regular periodic infinite lattices is usually surmised from kinetic Monte Carlo simulation on a finite lattice with periodic boundary conditions. However, the persistence of boundary effects for a Bethe lattice complicates this process, e.g., by inducing spatially heterogenous states. This motivates the exploration of various boundary conditions and unconventional simulation ensembles on the Bethe lattice to predict behavior for infinite size. Here, we focus on z = 3, and predict a discontinuous transition to the vacuum state on the infinite lattice when p exceeds a threshold value of around 0.053.
97 MATHEMATICS AND COMPUTING↗
Spontaneous Imbibition Velocities in Porous Metal Layers on Ceramic Substrates Calculated Using Microstructural Analyses
In this paper, the imbibition velocity of molten silver into a wettable porous nickel interlayer on a non-wetting sapphire substrate is investigated. Several imbibition models for porous media that use bulk characteristics and material properties are employed to predict the imbibition front location as a function of time. The model predictions are compared to experimental observations. It is found that pore size distribution is a better predictor of imbibition velocity than correlations based on permeability and porosity. The tortuosity of the microstructure is also found to have a significant effect and should be considered. A new model accounting for the different liquid contact angles on the underlying substrate and the porous interlayer material is proposed and achieves better agreement with the experimental observations.
36 MATERIALS SCIENCE↗
Lithium vanadium oxide (Li 1.1 V 3 O 8 ) thick porous electrodes with high rate capacity: utilization and evolution upon extended cycling elucidated via operando energy dispersive X-ray diffraction and continuum simulation
The phase distribution of lithiated LVO in thick (~500 μm) porous electrodes (TPEs) designed to facilitate both ion and electron transport was determined using synchrotron-based operando energy dispersive X-ray diffraction (EDXRD). Probing 3 positions in the TPE while cycling at a 1C rate revealed a homogeneous phase transition across the thickness of the electrode at the 1st and 95th cycles. Additionally, continuum modelling indicated uniform lithiation across the TPE in agreement with the EDXRD results and ascribed decreasing accessible active material to be the cause of loss in delivered capacity between the 1st and 95th cycles. The model was supported by the observation of significant particle fracture by SEM consistent with loss of electrical contact. Overall, the combination of operando EDXRD, continuum modeling, and ex situ measurements enabled a deeper understanding of lithium vanadium oxide transport properties under high rate extended cycling within a thick highly porous electrode architecture.
25 ENERGY STORAGE↗
Differential multi-probe thermal transport measurements of multi-walled carbon nanotubes grown by chemical vapor deposition
Carbon nanotubes (CNTs) are quasi-one dimensional nanostructures that display both high thermal conductivity for potential thermal management applications and intriguing low-dimensional phonon transport phenomena. In comparison to the advances made in the theoretical calculation of the lattice thermal conductivity of CNTs, thermal transport measurements of CNTs have been limited by either the poor temperature sensitivity of Raman thermometry technique or the presence of contact thermal resistance errors in sensitive two-probe resistance thermometry measurements. Here we report advances in a multi-probe measurement of the intrinsic thermal conductivity of individual multi-walled CNT samples that are transferred from the growth substrate onto the measurement device. The sample-thermometer thermal interface resistance is directly measured by this multi-probe method and used to model the temperature distribution along the contacted sample segment. The detailed temperature profile helps to eliminate the contact thermal resistance error in the obtained thermal conductivity of the suspended sample segment. A differential electro-thermal bridge measurement method is established to enhance the signal-to-noise ratio and reduce the measurement uncertainty by over 40%. The obtained thermal resistances of multiple suspended segments of the same MWCNT samples increase nearly linearly with increasing length, revealing diffusive phonon transport as a result of phonon-defect scattering in these MWCNT samples. The measured thermal conductivity increases with temperature and reaches up to 390 ± 20 W m –1 K –1 at room temperature for a 9-walled MWCNT. In conclusion, theoretical analysis of the measurement results suggests submicron phonon mean free paths due to extrinsic phonon scattering by extended defects such as grain boundaries. The obtained thermal conductivity is decreased by a factor of 3 upon electron beam damage and surface contamination of the CNT sample.
42 ENGINEERING↗
Experimental and theoretical description of the process of contact laser surgery with a titanium-doped optothermal fibre converter
In an in vitro experiment simulating a surgeon’s actions in the process of contact laser surgery of soft biological tissue, the dependences of the temperature of a titanium-doped optothermal fibre converter (TOTFC) and the depths of coagulation and ablation of biological tissue on the average radiation power of a diode laser with a wavelength of 980 nm and on the speed of the converter movement along biological tissue are obtained. The structural, optical, and thermophysical models of TOTFC are discussed, as well as the thermophysical model of the interaction of a laser-heated converter with biological tissue, which takes into account the temperature dependences of the basic thermophysical parameters of the converter and biological tissue, as well as the contribution of the thickness h{sub int} of the water vapour layer between the converter and biotissue. It is shown that the proposed model allows describing the result of contact laser surgery of soft biotissue with TOTFC adequately to the experiment. (laser biophotonics)
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase
Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.
Chemistry↗
Mechanisms influencing and prediction of tool influence function spots during hemispherical sub-aperture tool polishing on fused silica
Sub-aperture tool polishing of precision optics requires a detailed understanding of the local material removal [tool influence function (TIF)] at the contact spot between the workpiece and tool to achieve high removal determinism and hence precision of the optic relative to the desired/design surface figure. In this study, the mechanisms influencing and the quantitative prediction of the removal rate and shape of TIF spots during polishing of fused silica glass with cerium oxide slurry using a rotating hemispherical pad-foam tool for a wide variety of process conditions (including tool properties, kinematics, and applied displacements) are investigated. The TIF volumetric removal rate can be estimated utilizing the average relative velocity and contact area using a simple analytical model. In addition, stability of the volumetric removal rate for fixed process conditions is shown to be greatly dependent on the pad preparation and amount of tool use (affecting both pad topography and slurry buildup), whose general behavior shows an increase in removal rate followed by stabilization with polishing time. The determination of the TIF removal shape is more complex. An extended version of the Preston removal model is developed to explain a comprehensive set of measured TIF removal shapes to within ∼ <#comment/> 22 % <#comment/> . This model incorporates a number of phenomena impacting the TIF removal shape including: (a) temporal and spatial dependent relative velocity between the workpiece and tool; (b) an elastic mechanics based, as well as hydrodynamic, pressure distribution; (c) a spatially dependent friction coefficient possibly caused by both reduced slurry replenishment in low velocity regions and pad slurry islands (100 µm scale) and porosity (millimeter scale); and (d) a shear-based removal mechanism on the periphery of the contact spot.
Suratwala, T. (ORCID:0000000190861039)↗
Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning
While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.
36 MATERIALS SCIENCE↗
Role of water model on ion dissociation at ambient conditions
In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗