Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Force on Force”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Force field for calcium sulfate minerals to predict structural, hydration, and interfacial properties

Calcium sulfates such as anhydrite, hemihydrate, and gypsum find widespread use in building materials, implants, and tissue healing. We introduce a simple and compatible atomistic force field for all calcium sulfate phases that reproduces a wide range of experimental data including lattice parameters, surface, hydration, mechanical, and thermal properties in 1% to 5% accuracy relative to experiments. The performance is several times better than prior force fields and DFT methods, which lead to errors in structures and energies up to 100%. We explain (hkl) cleavage energies, the dynamics of (hkl) water interfaces, and new insights into molecular origins of crystal-facet specific hydration and solubility. Impressive agreement of computed and experimentally measured hydration energies is shown. The models add to the Interface force field (IFF) and are compatible with multiple force fields (CHARMM, AMBER, GROMOS, CVFF, PCFF, OPLS-AA) for property predictions of sulfate-containing materials from atoms to the large nanometer scale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular level study of uni/multi-axial deformation response of tobermorite 11 Å: A force field comparison study

Molecular deformation response of tobermorite 11 Å (a mineral analogue to C-S-H, the primary constituent of hydrated cement paste) under uni/tri-axial loading (both tension and compression) has been investigated in this manuscript using empirical and reactive force fields. A comparison has been chalked out amongst the responses using different force fields. Evolution of atomic environment during the deformation have studied using reactive force field. This study demonstrate the key observations on the application of different force fields to study the large deformation behaviour of the material.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗
DOE OSTI2020

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Refined Classical Force Field for Choline Chloride and Ethylene Glycol Mixtures over Wide Composition Range

In this work, a classical force field based on the General Amber Force Field (GAFF) was refined for the simulation of choline chloride (ChCl) and ethylene glycol (EG) mixtures over a wide composition range by scaling the partial charges and van der Waals parameters. The scaling factors were derived by fitting the simulation results to only eight experimental density and viscosity data points of pure EG, and ChCl/EG mixtures at 1:2, 1:4, and 1:6 molar ratios. Using the refined force field, properties essential for electrochemical applications such as density, viscosity, self-diffusion coefficient, and ionic conductivity were calculated, and excellent agreement to experimental results was found even for compositions and temperatures not used in the fitting procedure. In addition, new experimental data for density, viscosity, and ionic conductivity are reported as a function of temperature and composition for this mixture. To the best of our knowledge, this is the first classical force field developed for the study of ChCl/EG mixtures over a composition range that includes the eutectic point. Using the new model, the liquid dynamics was studied in terms of ionic conductivity. It was found that the dynamics in ChCl/EG mixtures with ChCl mole fraction higher than 20% is similar to that of ionic liquids, high temperature molten salts, and highly concentrated water-in-salt electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2021

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. Here, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

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

Reparameterizing a Lipid Force Field Using Small-Angle X-ray Scattering to Improve Predictions of Multicomponent Membranes under Organic Solvent Stress

Understanding and predicting properties of lipid-bilayer membranes are essential to elucidating organismal physiology and pathophysiology. Therefore, substantial efforts have been undertaken to derive accurate molecular-mechanics force fields (FF) to allow simulation of their properties; however, much of these past efforts focused on tuning force fields to accurately reproduce the properties of model single-component membranes and not multicomponent or stressed membranes. Here, we tuned the CHARMM36 FF for a simple 2-component model of a Gram-positive bacterium. This updated force field is parametrized against a single condensed-phase property, the small-angle X-ray scattering (SAXS) intensities of the 2-component lipid system, using a version of ForceBalance implemented previously (named ForceBalance-SAS) with SAXS and small-angle neutron scattering (SANS) intensities as optimization targets. After tuning, we observe improved agreement with the experimental SAXS of the membrane under 1-butanol and tetrahydrofuran solvent stresses, with a factor of 10.0 and 7.5 reduction, respectively, in χ 2 , the measure of the discrepancy between experimental and computed SAXS intensities. Furthermore, this reparametrized force field yielded improved agreement between experimental and simulated membrane thicknesses for the pure lipid systems. However, we note limitations in transferability and diagnose the source of such limitations, particularly the need for SANS in addition to SAXS.

Demerdash, Omar N. A. [Oak Ridge National Laborato↗
DOE OSTI2023

Solution Structure and Hydration Forces between Mica and Hydrophilic Versus Hydrophobic Surfaces

Solid-liquid interfaces are central to a range of interesting phenomena including catalysis, heterogeneous nucleation, water desalination, and biomolecular assembly. While three-dimensional Fast Force Mapping (3D FFM) has emerged as a technique for resolving interfacial solution structure at the molecular scale, key challenges for data interpretation persist, most notably regarding the influence of the probe on the measured structure. Using the mica-water system as a case study, we investigate the effect of hydrophilic and hydrophobic probes on interfacial solution structure measured by 3D FFM. Data from hydrophilic silicon-based probes are in good agreement with molecular dynamics simulations, wherein the innermost water molecules adsorb preferentially at the surface ditrigonal cavity sites followed by two subsequent hydration layers. In contrast, the hydrophobic carbon-based probes detect vertical oscillatory features, but do not show lateral patterning that matches the underlying mica lattice. At high ionic strength, up to six of these oscillatory features are observed extending 2 nm into the solution phase with an average spacing of 0.29 ± (0.04) nm. Further, we also determine that the repulsive hydration force between mica and the hydrophilic probe depends on the nature and concentration of ions in solution. Specifically, solutions with stronger ion-water and ion-ion interactions increase the interfacial viscosity, which results in a stronger repulsive hydration force as the probe approaches the surface. Based on these observations, we present a scheme for controlling the outcomes of particle attachment and aggregation by varying the solution conditions to tune the hydration force.

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

Indicator Patterns of Forced Change Learned by an Artificial Neural Network

Abstract Many problems in climate science require the identification of signals obscured by both the “noise” of internal climate variability and differences across models. Following previous work, we train an artificial neural network (ANN) to predict the year of a given map of annual‐mean temperature (or precipitation) from forced climate model simulations. This prediction task requires the ANN to learn forced patterns of change amidst a background of climate noise and model differences. We then apply a neural network visualization technique (layerwise relevance propagation) to visualize the spatial patterns that lead the ANN to successfully predict the year. These spatial patterns thus serve as “reliable indicators” of the forced change. The architecture of the ANN is chosen such that these indicators vary in time, thus capturing the evolving nature of regional signals of change. Results are compared to those of more standard approaches like signal‐to‐noise ratios and multilinear regression in order to gain intuition about the reliable indicators identified by the ANN. We then apply an additional visualization tool (backward optimization) to highlight where disagreements in simulated and observed patterns of change are most important for the prediction of the year. This work demonstrates that ANNs and their visualization tools make a powerful pair for extracting climate patterns of forced change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Revisiting Interior Water Mass Responses to Surface Forcing Changes and the Subsequent Effects on Overturning in the Southern Ocean

Abstract Two coupled climate models, differing primarily in horizontal resolution and treatment of mesoscale eddies, were used to assess the impact of perturbations in wind stress and Antarctic ice sheet (AIS) melting on the Southern Ocean meridional overturning circulation (SO MOC), which plays an important role in global climate regulation. The largest impact is found in the SO MOC lower limb, associated with the formation of Antarctic Bottom Water (AABW), which in both models is enhanced by wind and weakened by AIS meltwater perturbations. Even though both models under the AIS melting perturbation show similar AABW transport reductions of 4–5 Sv (50%–60%), the volume deflation of AABW south of 30°S is four times greater in the higher resolution simulation (−20 vs. −5 Sv). Water mass transformation (WMT) analysis reveals that surface‐forced dense water formation on the Antarctic shelf is absent in the higher resolution and reduced by half in the lower resolution model in response to the increased AIS melting. However, the decline of the AABW volume (and its inter‐model difference) far exceeds the surface‐forced WMT changes alone, which indicates that the divergent model responses arise from interactions between changes in surface forcing and interior mixing processes. This model divergence demonstrates an important source of uncertainty in climate modeling, and indicates that accurate shelf processes together with scenarios accounting for AIS melting are necessary for robust projections of the deep ocean's response to anthropogenic forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unsteady Land-Sea Breeze Circulations in the Presence of a Synoptic Pressure Forcing

Unsteady land-sea breezes (LSBs) that result from time-varying surface temperature contrasts Δθ(t) are explored in the presence of a constant synoptic pressure forcing, M g , oriented from sea to land (α = 0°) or land to sea (α = 180°). Large eddy simulations reveal the development of four distinctive regimes, depending on the joint interaction between M g , α, and Δθ(t) in modulating the fine-scale dynamics. Time lags, computed as the shifts that maximize correlation coefficients of the velocity between the unsteady and the corresponding steady scenarios at Δθ = Δθ max , are found to be significant and to extend 2 hr longer for α = 0° compared to α = 180°. These diurnal dynamics result in nonequilibrium conditions that are significantly affected by the flow history, and that behave differently over the two patches for the different α’s. Turbulence is found to be out of equilibrium with the mean flow, and the mean itself is found to be out of equilibrium with the thermal forcing. The sea surface heat flux is consistently more sensitive than its land counterpart to the time-varying external forcing Δθ(t), and more so for synoptic forcing from land to sea (α = 180°). Hence, although the land reaches equilibrium faster, the sea patch is found to exert a stronger control on the turbulence-mean flow equilibrium response. Finally, the vertical velocity profile at the shore and shore-normal velocity transects at the first grid level are shown to encode the multiscale regimes of the LSBs evolution and can thus be used to identify these regimes using k-means clustering.

58 GEOSCIENCES↗
DOE OSTI2020

A joint role for forced and internally-driven variability in the decadal modulation of global warming

Despite the observed monotonic increase in greenhouse-gas concentrations, global mean temperature displays important decadal fluctuations typically attributed to both external forcing and internal variability. Here, we provide a robust quantification of the relative contributions of anthropogenic, natural, and internally-driven decadal variability of global mean sea surface temperature (GMSST) by using a unique dataset consisting of 30-member large initial-condition ensembles with five Earth System Models (ESM-LE). We present evidence that a large fraction (~29–53%) of the simulated decadal-scale variance in individual timeseries of GMSST over 1950–2010 is externally forced and largely linked to the representation of volcanic aerosols. Comparison with the future (2010–2070) period suggests that external forcing provides a source of additional decadal-scale variability in the historical period. Given the unpredictable nature of future volcanic aerosol forcing, it is suggested that a large portion of decadal GMSST variability might not be predictable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Low cost and massively parallel force spectroscopy with fluid loading on a chip

Current approaches for single molecule force spectroscopy are typically constrained by low throughput and high instrumentation cost. Herein, a low-cost, high throughput technique is demonstrated using microfluidics for multiplexed mechanical manipulation of up to ~4000 individual molecules via molecular fluid loading on-a-chip (FLO-Chip). The FLO-Chip consists of serially connected microchannels with varying width, allowing for simultaneous testing at multiple loading rates. Molecular force measurements are demonstrated by dissociating Biotin-Streptavidin and Digoxigenin-AntiDigoxigenin interactions along with unzipping of double stranded DNA of varying sequence under different dynamic loading rates and solution conditions. Rupture force results under varying loading rates and solution conditions are in good agreement with prior studies, verifying a versatile approach for single molecule biophysics and molecular mechanobiology. FLO-Chip enables straightforward, rapid, low-cost, and portable mechanical testing of single molecules that can be implemented on a wide range of microscopes to broaden access and may enable new applications of molecular force spectroscopy.

42 ENGINEERING↗
DOE OSTI2023

Observation of electron orbital signatures of single atoms within metal-phthalocyanines using atomic force microscopy

Resolving the electronic structure of a single atom within a molecule is of fundamental importance for understanding and predicting chemical and physical properties of functional molecules such as molecular catalysts. However, the observation of the orbital signature of an individual atom is challenging. We report here the direct identification of two adjacent transition-metal atoms, Fe and Co, within phthalocyanine molecules using high-resolution noncontact atomic force microscopy (HR-AFM). HR-AFM imaging reveals that the Co atom is brighter and presents four distinct lobes on the horizontal plane whereas the Fe atom displays a “square” morphology. Pico-force spectroscopy measurements show a larger repulsion force of about 5 pN on the tip exerted by Co in comparison to Fe. Our combined experimental and theoretical results demonstrate that both the distinguishable features in AFM images and the variation in the measured forces arise from Co’s higher electron orbital occupation above the molecular plane. The ability to directly observe orbital signatures using HR-AFM should provide a promising approach to characterizing the electronic structure of an individual atom in a molecular species and to understand mechanisms of certain chemical reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
DOE OSTI2024

eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.

08 HYDROGEN↗