Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic contacts”

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 OSTI2020

Pressure Drop with Moving Contact Lines and Dynamic Contact Angles in a Hydrophobic Round Minichannel: Visualization via Synchrotron X-ray Imaging and Verification of Experimental Correlations

In hydrophobic mini- and microchannels, slug flow with moving contact lines is typically generated under various two-phase flow conditions. There is a significant pressure drop in this flow pattern with moving contact lines, which is closely related to the dynamic contact angles. Researchers have investigated dynamic contact angles experimentally for decades, but due to the limitations of visualization techniques, these experiments have typically been conducted in low Weber number regions (We < 10 –3 ). In this study, we clearly visualized the dynamic contact angles of a liquid slug in high Weber number regions (10 –3 < We <1) via synchrotron X-ray imaging with high temporal (~1000 fps) and spatial (~ 2 μm/pixel) resolutions. We precisely measured the pressure drop with moving contact lines in a hydrophobic minichannel (inner diameter = 1.018 mm). On the basis of our experimental data, we verified previous correlations for dynamic contact angles and explored the relationship between pressure drop with moving contact lines and dynamic contact angles.

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

Hydrodynamics of countercurrent flows in a structured packed column: effects of initial wetting and dynamic contact angle

The structured packed column has gained considerable attention as a prospective technology to mitigate carbon dioxide emissions from thermal power plants. The computational flow investigation of such systems is a multiscale problem, and local hydrodynamics plays a key role in overall column efficiency. Accordingly, multiphase flow studies in the representative elementary unit (REU) of a packed column can offer insights into the local flow behavior such as flow pattern, interfacial area, etc. The interfacial area plays a key role in the overall interphase mass transfer and thereby in separation efficiency. The predicted effective area matches well with corresponding ones in experiments for Mellapak 250.Y packings. The contact angle (i.e., solid substrate behavior) is one of the critical factors that dictates wettings, and thereby the interfacial area. The dynamic contact angle (DCA) was considered in order to explore the effects of contact angle hysteresis on the interfacial area. DCA has more pronounced impact on the interfacial area for a solvent possessing a higher value of surface tension and equilibrium contact angle. The interfacial area shows temporal undulation and does not achieve a pseudo steady state due to the capillary force. In contrary, the interfacial area gets a net value for a solvent having low surface tension value. Further, specification of the static contact angle at the solid substrate shows the least value of the interfacial area. We also extensively studied the effect of the initial sheet condition (dry vs wet) on the interfacial area. The initially wetted sheets show slightly higher value of the interfacial area as compared to the initially dry sheets at a fixed flow rate for a solvent having high surface tension value. Wetting hysteresis decreases with increasing liquid loads. On the other hand, wetting hysteresis does not occur for a solvent having value of low surface tension and equilibrium contact angle.

structured packings, dynamic contact angle, interf↗
DOE OSTI2023

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

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

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2026

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗
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 OSTI2021

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
DOE OSTI2021

Drop impact tester development for spent nuclear fuel vibration integrity study

The accumulated damage of SNF rod is defined as the integrated damage induced by the SNF assembly harmonic vibrations and SNF assembly & cask system contact interactions induced transient shocks, experienced by SNF rods during normal conditions of transportation (NCT). The CIRFT technology has provided means to resolve/understand the SNF rod harmonic vibration reliability issues associated with maximum axial clad tubing stress induced by SNF inertia related dynamic bending loading. However, the SNF rod integrity in association with the transient shock loads during NCT through dynamic contact interactions among SNF assembly cannot be investigated systematically through CIRFT technology alone. Thus, the proposed drop impact tester development become a viable tool for SNF system transient shock effect investigation. Furthermore, the transient shocks induced loading mode is normally perpendicular to the SNF rod axial orientation, i.e., in a transverse orientation to the SNF rod. To proceed with such loading mode investigation, a new device, “SNF Drop Impact Tester” that provides the transverse impact shock load onto a SNF rod through a free-drop projectile was developed in FY20 for SNF System Vibration Reliability Investigation; where the SNF system contact impact intensity can be calibrated accordingly with projectile weight and projectile head geometry. This newly developed research tool can further provide a detailed understanding about the effect of dynamic contact-interaction loadings; and combined harmonic vibration and transient shock loading modes on the fatigue damage evolution of the HBU SNF under NCT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2020

Current measurements in the intermittent-contact mode of atomic force microscopy using the Fourier method: a feasibility analysis

Atomic force microscopy (AFM) is an important tool for measuring a variety of nanoscale surface properties, such as topography, viscoelasticity, electrical potential and conductivity. Some of these properties are measured using contact methods (static contact or intermittent contact), while others are measured using noncontact methods. Some properties can be measured using different approaches. Conductivity, in particular, is mapped using the contact-mode method. However, this modality can be destructive to delicate samples, since it involves continuously dragging the cantilever tip on the surface during the raster scan, while a constant tip–sample force is applied. In this paper we discuss a possible approach to develop an intermittent-contact conductive AFM mode based on Fourier analysis, whereby the measured current response consists of higher harmonics of the cantilever oscillation frequency. Such an approach may enable the characterization of soft samples with less damage than contact-mode imaging. To explore its feasibility, we derive the analytical form of the tip–sample current that would be obtained for attractive (noncontact) and repulsive (intermittent-contact) dynamic AFM characterization, and compare it with results obtained from numerical simulations. Although significant instrumentation challenges are anticipated, the modelling results are promising and suggest that Fourier-based higher-harmonics current measurement may enable the development of a reliable intermittent-contact conductive AFM method.

14 SOLAR ENERGY↗
DOE OSTI2022

Interlaminar improvement of carbon fiber/epoxy composites via fluorine‐containing high‐epoxy‐value sizing agent

Abstract To improve the interfacial bonding effect of carbon fibers (CFs) and epoxy resins, we prepared a fluorine‐containing high‐epoxy‐value sizing agent (named TS) in this article. It is an emulsion with a more uniform particle size distribution and higher stability better than a commercial Takemoto sizing agent. In addition, sulfonic acid groups can improve the hydrophilicity of the TS, and the presence of epoxy groups enhances the wettability between CFs and the matrix. The properties of the sizing agents were examined by Fourier transform infrared spectroscopy, dynamic light scattering, and thermogravimetric analysis. The sized CFs were evaluated by scanning electron microscopy, x‐ray photoelectron spectroscopy, monofilament tensile test, and dynamic contact angle. Afterward, CF/epoxy composites were prepared, and their fracture morphology and interlaminar shear strength (ILSS) were examined. Results indicated that the TS sizing significantly increased the content of activated carbon on the surface of the CFs and improved the bonding effect between the fiber and the matrix. After TS was sized, the ILSS increased from 49.64 to 58.74 MPa, 18.34% higher than unsized CF/epoxy composites, showing the prominent interfacial reinforcement better than commercial Takemoto.

Qiu, Zhijie↗
DOE OSTI2024

Synthesizing realistic sand assemblies with denoising diffusion in latent space

Abstract The shapes and morphological features of grains in sand assemblies have far‐reaching implications in many engineering applications, such as geotechnical engineering, computer animations, petroleum engineering, and concentrated solar power. Yet, our understanding of the influence of grain geometries on macroscopic response is often only qualitative, due to the limited availability of high‐quality 3D grain geometry data. In this paper, we introduce a denoising diffusion algorithm that uses a set of point clouds collected from the surface of individual sand grains to generate grains in the latent space. By employing a point cloud autoencoder, the three‐dimensional point cloud structures of sand grains are first encoded into a lower‐dimensional latent space. A generative denoising diffusion probabilistic model is trained to produce synthetic sand that maximizes the log‐likelihood of the generated samples belonging to the original data distribution measured by a Kullback‐Leibler divergence. Numerical experiments suggest that the proposed method is capable of generating realistic grains with morphology, shapes and sizes consistent with the training data inferred from an F50 sand database. We then use a rigid contact dynamic simulator to pour the synthetic sand in a confined volume to form granular assemblies in a static equilibrium state with targeted distribution properties. To ensure third‐party validation, 50,000 synthetic sand grains and the 1542 real synchrotron microcomputed tomography (SMT) scans of the F50 sand, as well as the granular assemblies composed of synthetic sand grains are made available in an open‐source repository.

Vlassis, Nikolaos N.↗
DOE OSTI2020

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗