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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling Contact Angle vs. Temperature for the Quartz-Water-Decane System

Innovative approaches are needed to improve the efficiency of oil recovery technologies to meet the growing demands of fossil-fuel based energy consumption. Enhanced oil recovery (EOR) methods such as low-salinity waterflooding and chemically tuned waterflooding aim to optimize the reservoir’s wetting properties, detaching oil globules from rock surfaces and allowing easier oil flow through pore throats. This wetting behavior is commonly quantified by contact angle measurements of the rock-oil-brine interface, which have been thoroughly investigated and theorized for many systems at ambient temperatures and pressures. However, few studies exist for extending contact angle theories away from ambient conditions. In this paper, we model the contact angles of a quartz-water-decane system at elevated temperatures using the surface tension component (STC) approach. Temperature-dependent van der Waals [Lifshitz-van der Waals (LW)] interactions and hydrogen-bonding (acid-base) interactions were calculated and are incorporated into the model for the quartz-water-decane interface. Additionally, the Hough and White procedure was used to create temperature-dependent dielectric functions of quartz, water, and normal decane for calculations of Hamaker coefficients. Hamaker coefficients calculated this way are highly linear with temperature and agree well with Israelachvili’s approximation. The acid-base interactions likely contribute the most to system wettability changes. Resulting contact angles of the quartz-water-decane system shift from water-wet (16°) to slightly water-wet (57.4°) as temperature increases. The model was also successfully verified for the quartz-air-water system. Our results can be used in future studies to determine optimal injected water compositions for specific rock-oil-brine and other systems with consideration of reservoir temperature.

02 PETROLEUM↗
DOE OSTI2023

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

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

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

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

Population genomics and history of speciation reveal fishery management gaps in two related redfish species ( Sebastes mentella and Sebastes fasciatus )

Abstract Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus . We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus . Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter‐ and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

Benestan, Laura M.↗
DOE OSTI2023

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

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

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

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

Three-Dimensional Evaluation of Sand Particle Fracture Using Discrete-Element Method and Synchrotron Microcomputed Tomography Images

Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.

58 GEOSCIENCES↗
DOE OSTI2024

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

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↗