Search NASA⌕ Search

SEARCH · Search NASA

Results for “force method”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Application of the Most Likely Extreme Response Method for Wave Energy Converters: Preprint

Extreme loads are often a key cost driver for wave energy converters (WECs). As an alternative to exhaustive Monte Carlo or long-term simulations, the most likely extreme response (MLER) method allows mid- and high-fidelity simulations to be used more efficiently in evaluating WEC response to events at the edges of the design envelope, and is therefore applicable to system design analysis. The study discussed in this paper applies the MLER method to investigate the maximum heave, pitch, and surge force of a point absorber WEC. Most likely extreme waves were obtained from a set of wave statistics data based on spectral analysis and the response amplitude operators (RAOs) of the floating body; the RAOs were computed from a simple radiation-and-diffraction-theory-based numerical model. A weakly nonlinear numerical method and a computational fluid dynamics (CFD) method were then applied to compute the short-term response to the MLER wave. Effects of nonlinear wave and floating body interaction on the WEC under the anticipated 100-year waves were examined by comparing the results from the linearly superimposed RAOs, the weakly nonlinear model, and CFD simulations. Overall, the MLER method was successfully applied. In particular, when coupled to a high-fidelity CFD analysis, the nonlinear fluid dynamics can be readily captured.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Residually Stressed Bimaterial Beam Specimen for Measuring Environmentally Assisted Crack Growth

Background: Subcritical crack growth can occur in a brittle material when the stress intensity factor is smaller than the fracture toughness if an oxidizing agent (such as water) is present at the crack tip. Objective: Here we present a novel bi-material beam specimen which can measure environmentally assisted crack growth rates. The specimen is “self-loaded” by residual stress and requires no external loading. Methods: Two materials with different coefficient of thermal expansion are diffusion bonded at high temperature. After cooling to room temperature a subcritical crack is driven by thermal residual stresses. A finite element model is used to design the specimen geometry in terms of material properties in order to achieve the desired crack tip driving force. Results: The specimen is designed so that the crack driving force decreases as the crack extends, thus enabling the measurement of the crack velocity versus driving force relationship with a single test. The method is demonstrated by measuring slow crack growth data in soda lime silicate glass and validated by comparison to previously published data. Conclusions: The self-loaded nature of the specimen makes it ideal for measuring the very low crack velocities needed to predict brittle failure at long lifetimes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shape-constrained Input Estimation for Efficient Multi-shaker Vibration Testing

Multi-shaker vibration testing is gaining interest from structural dynamics test engineers as it can provide a much more accurate match to complicated field vibration responses than traditional single-axis shaker tests. However, the force capabilities of the small modal shakers typically used in multi-shaker vibration tests has limited the achievable response levels. To date, most multi-shaker vibration tests have been performed using a variety of standard, commercially-available control systems. While these control systems are adequate for a wide range of multiple-input/multiple-output tests, their control algorithms have not been tailored for the specific problem of multi-shaker vibration tests: efficiently coordinating the various shakers to work together to achieve a desired response. Here, a new input estimation algorithm is developed and demonstrated using simulations and actual test data. This algorithm, dubbed shape-constrained input estimation, is shown to effectively coordinate multiple shakers using a set of constraint vectors based on the deflection shapes of the test structure. This is accomplished by using the singular vector shapes of the system frequency response matrix, which allows the constraint vectors to automatically change as a function of frequency. Simulation and test results indicate a significant reduction in the input forces required to achieve a desired response. Finally, the results indicate that shape-constrained input estimation is an effective method to achieve higher response levels from limited shaker forces which will enable higher level multi-shaker vibration tests to be performed.

42 ENGINEERING↗
DOE OSTI2020

Surprising similarities in model and observational aerosol radiative forcing estimates

Abstract. The radiative forcing from aerosols (particularly through their interaction with clouds) remains one of the most uncertain components of the human forcing of the climate. Observation-based studies have typically found a smaller aerosol effective radiative forcing than in model simulations and were given preferential weighting in the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). With their own sources of uncertainty, it is not clear that observation-based estimates are more reliable. Understanding the source of the model and observational differences is thus vital to reduce uncertainty in the impact of aerosols on the climate. These reported discrepancies arise from the different methods of separating the components of aerosol forcing used in model and observational studies. Applying the observational decomposition to global climate model (GCM) output, the two different lines of evidence are surprisingly similar, with a much better agreement on the magnitude of aerosol impacts on cloud properties. Cloud adjustments remain a significant source of uncertainty, particularly for ice clouds. However, they are consistent with the uncertainty from observation-based methods, with the liquid water path adjustment usually enhancing the Twomey effect by less than 50 %. Depending on different sets of assumptions, this work suggests that model and observation-based estimates could be more equally weighted in future synthesis studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Development of a New Method to Investigate the Dynamic Friction Behavior of Interfaces Using a Kolsky Tension Bar

In order to understand interfacial interaction of a bi-material during an impact loading event, the dynamic friction coefficient is one of the key parameters that must be characterized and quantified. In this study, a new experimental method to determine the dynamic friction coefficient between two metals was developed by using a Kolsky tension bar and a custom-designed friction fixture. Polyvinylidene fluoride (PVDF) force sensors were used to measure the normal force applied to the friction tribo pairs and the friction force was measured with conventional Kolsky tension bar method. To evaluate the technique, the dynamic friction coefficient between 4340 steel and 7075-T6 aluminum was investigated at an impact speed of approximately 8 m/s. Additionally, the dynamic friction coefficient of the tribo pairs with varied surface roughness was also investigated. The data suggest that higher surface roughness leads to higher friction coefficients at the same speed of 8 m/s.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A practical perspective on the implementation of hyperdynamics for accelerated simulation

Consideration is given to several practical issues arising during the implementation of hyperdynamics, a methodology that extends the time scale of the conventional molecular dynamics simulation potentially by orders of magnitude. First, the methodology is reformulated in terms of the transition rate based on the buffer region approach (buffer rate), which can describe transitions in more general contexts than the transition state theory (TST). It will be shown that hyperdynamics can exactly preserve the buffer rate as well as the TST rate, which broadens the scope of the method. Next, the originally proposed scheme to compute the boost factor on-the-fly is reviewed and some alternative methods, one of which uses the umbrella sampling method, are presented. Finally, the methodology is validated in the context of a 1-dimensional example potential and a 3-dimensional simulation of the motion of an atomic force microscope tip moving along a surface.

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

Unrealized Global Temperature Increase: Implications of Current Uncertainties

Unrealized increase in global mean surface air temperature (GMST) may result from the climate system not being in steady state with forcings and/or from cessation of negative aerosol forcing that would result from decreases in emissions. An observation-constrained method is applied to infer the dependence of Earth's climate sensitivity on forcing by anthropogenic aerosols within the uncertainty on that forcing given by the Fifth (2013) Assessment Report of the Intergovernmental Panel on Climate Change. Within these uncertainty ranges the increase in GMST due to temperature lag for future forcings held constant is slight (0.09–0.19 K over 20 years; 0.12–0.26 K over 100 years). However the incremental increase in GMST that would result from a hypothetical abrupt cessation of sources of aerosols could be quite large, but is highly uncertain, 0.1–1.3 K over 20 years. Decrease in CO2 abundance and forcing following abrupt cessation of emissions would offset these increases in GMST over 100 years by as little as 0.09 K to as much as 0.8 K. The uncertainties quantified here greatly limit confidence in projections of change in GMST that would result from any strategy for future reduction of emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Materials perspective on Casimir and van der Waals interactions

Interactions induced by electromagnetic fluctuations, such as van der Waals and Casimir forces, are of universal nature present at any length scale between any types of systems. In such interactions these are important not only for the fundamental science of materials behavior, but also for the design and improvement of micro- and nanostructured devices. In the past decade, many new materials have become available, which has stimulated the need for understanding their dispersive interactions. The field of van der Waals and Casimir forces has experienced an impetus in terms of developing novel theoretical and computational methods to provide new insights into related phenomena. The understanding of such forces has far reaching consequences as it bridges concepts in materials, atomic and molecular physics, condensed-matter physics, high-energy physics, chemistry, and biology. Our review summarizes major breakthroughs and emphasizes the common origin of van der Waals and Casimir interactions. Progress related to novel ab initio modeling approaches and their application in various systems, interactions in materials with Dirac-like spectra, force manipulations through nontrivial boundary conditions, and applications of van der Waals forces in organic and biological matter are examined. Finally, the outlook of the review is to give the scientific community a materials perspective of van der Waals and Casimir phenomena and stimulate the development of experimental techniques and applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Robust and scalable uncertainty estimation with conformal prediction for machine-learned interatomic potentials

Abstract Uncertainty quantification (UQ) is important to machine learning (ML) force fields to assess the level of confidence during prediction, as ML models are not inherently physical and can therefore yield catastrophically incorrect predictions. Established a-posteriori UQ methods, including ensemble methods, the dropout method, the delta method, and various heuristic distance metrics, have limitations such as being computationally challenging for large models due to model re-training. In addition, the uncertainty estimates are often not rigorously calibrated. In this work, we propose combining the distribution-free UQ method, known as conformal prediction (CP), with the distances in the neural network’s latent space to estimate the uncertainty of energies predicted by neural network force fields. We evaluate this method (CP+latent) along with other UQ methods on two essential aspects, calibration, and sharpness, and find this method to be both calibrated and sharp under the assumption of independent and identically-distributed (i.i.d.) data. We show that the method is relatively insensitive to hyperparameters selected, and test the limitations of the method when the i.i.d. assumption is violated. Finally, we demonstrate that this method can be readily applied to trained neural network force fields with traditional and graph neural network architectures to obtain estimates of uncertainty with low computational costs on a training dataset of 1 million images to showcase its scalability and portability. Incorporating the CP method with latent distances offers a calibrated, sharp and efficient strategy to estimate the uncertainty of neural network force fields. In addition, the CP approach can also function as a promising strategy for calibrating uncertainty estimated by other approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Force Reconstruction from Ejection Tests of Stores from Aircraft Used for Model Predictions and Missing/Bad Gages

One of the more severe environments for a store on an aircraft is during the ejection of the store. During this environment it is not possible to instrument all component responses, and it is also likely that some instruments may fail during the environment testing. This work provides a method for developing these responses from failed gages and uninstrumented locations. First, the forces observed by the store during the environment are reconstructed. A simple sampling method is used to reconstruct these forces given various parameters. Then, these forces are applied to a model to generate the component responses. Validation is performed on this methodology.

42 ENGINEERING↗
DOE OSTI2022

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Accuracy of tip-sample interaction measurements using dynamic atomic force microscopy techniques: Dependence on oscillation amplitude, interaction strength, and tip-sample distance

Atomic force microscopy (AFM) is a versatile surface characterization method that can map a sample’s topography with high spatial resolution while simultaneously interrogating its surface chemistry through the site-specific high-resolution quantification of the forces acting between the sample and the probe tip. Thanks to considerable advances in AFM measurement technology, such local measurements of chemical properties have gained much popularity in recent years. To this end, dynamic AFM methodologies are implemented where either the oscillation frequency or the oscillation amplitude and phase of the vibrating cantilever are recorded as a function of tip-sample distance and subsequently converted to reflect tip-sample forces or interaction potentials. Such conversion has, however, been shown to produce non-negligible errors when applying the most commonly used mathematical conversion procedures if oscillation amplitudes are of the order of the decay length of the interaction. Extending on these earlier findings, the computational study presented in this paper reveals that the degree of divergence from actual values may also critically depend on both the overall strength of tip-sample interaction and the distance at which the interaction is obtained. These systematic errors can, however, be effectively eliminated by using oscillation amplitudes that are sufficiently larger than the decay length of the interaction potential.

Electrical properties and parameters↗
DOE OSTI2017

Theoretical Calculation and Simulation Studies of Balancing Forces and Axisymmetric Poloidal Halo Current in EAST

A new technology for developing fusion energy is to use hydrogen isotopes i.e., deuterium (D) and tritium (T) which is a combine effort for building up of International Thermonuclear Reactor (ITER) named as Tokamak, which will come into operation in 2020. Handsome amount of work has already been performed by many scientists contingent with plasma shape, halo current and plasma equilibrium properties by numerical techniques. Halo current calculation is the most important work in fusion reactors. In this paper, halo current and plasma cross-section energy have been calculated. The research covers calculation of plasma energy and determines the force on the vessel by magnetic field transformation method, which is the essential part of calculating halo current. Further, vertical displacement event (VDE) up/down halo current models have been calculated. These models give total magnetic field, asymmetric forces, conducting points and poloidal halo current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

Metal organic frameworks (MOFs) are a large class of porous materials and have garnered significant interest due to their large surface areas and their tunable physical and chemical properties. Numerous prior studies have been performed to screen large databases of this material class for promising DAC sorbent materials. These studies have often relied on classical model potentials. While density functional theory (DFT) calculations have been shown to be very accurate for modeling the interaction of CO2 with MOFs, such calculations are too computationally demanding for statistically significant adsorption predictions. To overcome this barrier, we developed methods for training models to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using machine learned force fields (MLFFs). These methods were parametrized based on DFT calculations of CO2 in a flexible MOF and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
DOE OSTI2016

Population-based 3D genome structure analysis reveals driving forces in spatial genome organization

Significance We provide a method for population-based structure modeling of whole diploid genomes using Hi-C data. The method considers the stochastic nature of chromosome structures, which allows a detailed analysis of the dynamic landscape of genome organizations. We predict and experimentally validate the presence of chromosome-specific higher-order centromere clusters, which can play a key role in the spatial organization of the human genome, specifically influencing the overall chromosome positioning, as well as the preference of specific chromosome conformations. Our approach generate predictive structural models of diploid genomes from Hi-C data, which can provide insights into the guiding principles of 3D genome organizations.

3D genome organization↗