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 1,351 records · Page 75

Simulated vs. empirical weather responsiveness of crop yields: US evidence and implications for the agricultural impacts of climate change

Global gridded crop models (GGCMs) are the workhorse of assessments of the agricultural impacts of climate change. Yet the changes in crop yields projected by different models in response to the same meteorological forcing can differ substantially. Through an inter-method comparison, we provide a first glimpse into the origins and implications of this divergence—both among GGCMs and between GGCMs and historical observations. We examine yields of rainfed maize, wheat, and soybeans simulated by six GGCMs as part of the Inter-Sectoral Impact Model Intercomparison Project-Fast Track (ISIMIP-FT) exercise, comparing 1981–2004 hindcast yields over the coterminous United States (US) against US Department of Agriculture (USDA) time series for about 1000 counties. Leveraging the empirical climate change impacts literature, we estimate reduced-form econometric models of crop yield responses to temperature and precipitation exposures for both GGCMs and observations. We find that up to 60% of the variance in both simulated and observed yields is attributable to weather variation. A majority of the GGCMs have difficulty reproducing the observed distribution of percentage yield anomalies, and exhibit aggregate responses that show yields to be more weather-sensitive than in the observational record over the predominant range of temperature and precipitation conditions. In conclusion, this disparity is largely attributable to heterogeneity in GGCMs' responses, as opposed to uncertainty in historical weather forcings, and is responsible for widely divergent impacts of climate on future crop yields.

54 ENVIRONMENTAL SCIENCES↗

Toward Immersed Boundary Simulation of High Reynolds Number Flows

In the immersed boundary (IB) method, the surface of an object is reconstructed with forcing terms in the underlying flow field equations. The surface may split a computational cell removing the constraint of the near wall gridlines to be aligned with the surface. This feature greatly simplifies the grid generation process which is cumbersome and expensive in particular for structured grids and complex geometries. The IB method is ideally suited for Cartesian flow solvers. The flow equations written in Cartesian coordinates appear in a very simple form and several numerical algorithms can be used for an efficient solution of the equations. In addition, the accuracy of numerical algorithms is dependent on the underlying grid and it usually deteriorates when the grid deviates from a Cartesian mesh. The challenge for the IB method lies in the representation of the wall boundaries and in providing an adequate near wall flow field resolution. The issue of enforcing no-slip boundary conditions at the immersed surface has been addressed by several authors by imposing a local reconstruction of the solution. Initial work by Verzicco et al. was based on a simple linear, one-dimensional operator and this approach proved to be accurate for boundaries largely aligned with the grid lines. Majumdar et al. used various multidimensional and high order polynomial interpolations schemes. These high order schemes, however, are keen to introduce wiggles and spurious extrema. Iaccarino & Verzicco and Kalitzin & Iaccarino proposed a tri-linear reconstruction for the velocity components and the turbulent scalars. A modified implementation that has proven to be more robust is reported in this paper. The issue of adequate near wall resolution in a Cartesian framework can initially be addressed by using a non-uniform mesh which is stretched near the surface. In this paper, we investigate an unstructured approach for local grid refinement that utilizes Cartesian mesh features. The computation of high Reynolds number wall bounded flows is particularly challenging as it requires the consideration of thin turbulent boundary layers, i.e. near wall regions with large gradients of the flow field variables. For such flows, the representation of the wall boundary has a large impact on the accuracy of the computation. It is also critical for the robustness and convergence of the flow solver.

Kalitzin, Georgi↗

Estimates of the stratospheric residual circulation using the downward control principle

The transformed Eulerian-mean momentum and continuity equations are used to calculate the residual mean meridional circulations for the lower stratosphere and troposphere. Momentum and temperature fluxes required for the computation are estimated from U.K. Meteorological Office (UKMO) analyzed geopotential heights. National Center for Atmospheric Research (NCAR) CCM2 model output is used to assess the errors associated with the calculation. The model comparisons showed that the method works reasonably well for solstice seasons, but is inadequate for equinox seasons. In addition, it is found that some parameterization of gravity wave drag needs to be included with the planetary wave forcing to accurately estimate the residual mean circulation using this method.

Rosenlof, Karen H.↗

Serial segment method for measuring remnant

For tracking tasks where a sum of sine waves forcing function is used it is often difficult and/or expensive to obtain the pilot's remnant in the vicinity of the sine waves. For the case where each sine wave has at least four times an integer number of cycle per run length, this paper illustrates the serial segments method for measuring remnant power spectral density in a frequency band centered on each sine wave. This method can be implemented on digital, hybrid, or analog Fourier coefficient analyzers, and is particularly advantageous on the latter since properties of Fourier coefficients are exploited to yield both a remnant measure and an improved estimate of the correlated component.

Magdaleno, R. E.↗

What causes the spread of model projections of ocean dynamic sea-level change in response to greenhouse gas forcing?

Abstract Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level ( ζ ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ . The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific.

54 ENVIRONMENTAL SCIENCES↗

Method to Control Crimping Processes Using Ultrasonic Transmission Analysis

A method of crimping wires includes positioning a wire/terminal combination between first and second crimp forming tools. A force is applied to the crimp forming tools to deform the wire/terminal combination. The method further includes measuring ultrasonic energy that is transmitted across the wire/terminal combination as the terminal is being deformed. A rate of change of the magnitude of the ultrasonic energy is also determined as the terminal is being deformed. The crimping process is terminated if the rate of change of the magnitude of the ultrasonic energy falls below a predefined threshold level. Data gathered during the crimping process can also be utilized to determine if a faulty crimp has occurred.

Cramer, Karl Elliott↗

A new strategy to engineer polymer bulk heterojunction solar cells with thick active layers via self-assembly of the tertiary columnar phase

Here, we report that the addition of a non-photoactive tertiary polymer phase in the binary bulk heterojunction (BHJ) polymer solar cell leads to a self-assembled columnar nanostructure, enhancing the charge mobilities and photovoltaic efficiency with surprisingly increased optimal active blend thicknesses over 300 nm, 3–4 times larger than that of the binary counterpart. Using the prototypical poly(3-hexylthiophene) (P3HT):fullerene blend as a model BHJ system, we discover that the inert poly(methyl methacrylate) (PMMA) added in the binary BHJ blend self-assembles into vertical columns, which not only template the phase segregation of electron acceptor fullerenes but also induce the out-of-plane rotation of the edge-on-orientated crystalline P3HT phase. Using complementary interrogation methods including neutron reflectivity, X-ray scattering, atomic force microscopy, transmission electron microscopy, and molecular dynamics simulations, we show that the enhanced charge transport originates from the more randomized molecular stacking of the P3HT phase and the spontaneous segregation of fullerenes at the P3HT/PMMA interface, driven by the high surface tension between the two polymeric components. The results demonstrate a potential method for increasing the thicknesses of high-performance polymer BHJ solar cells with improved photovoltaic efficiency, alleviating the burden of stringently controlling the ultrathin blend thickness during the roll-to-roll-type large-area manufacturing environment.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Steady and Oscillatory, Subsonic and Supersonic, Aerodynamic Pressure and Generalized Forces for Complex Aircraft Configurations and Applications to Flutter

A general method for analyzing aerodynamic flows around complex configurations is presented. By applying the Green function method, a linear integral equation relating the unknown, small perturbation potential on the surface of the body, to the known downwash is obtained. The surfaces of the aircraft, wake and diaphragm (if necessary) are divided into small quadrilateral elements which are approximated with hyperboloidal surfaces. The potential and its normal derivative are assumed to be constant within each element. This yields a set of linear algebraic equations and the coefficients are evaluated analytically. By using Gaussian elimination method, equations are solved for the potentials at the centroids of elements. The pressure coefficient is evaluated by the finite different method; the lift and moment coefficients are evaluated by numerical integration. Numerical results are presented, and applications to flutter are also included.

Chen, L. T.↗

A study on the forced vibration of a Timoshenko beam

By using Galerkin's variational method an approximate solution is developed for a system of two differential equations with linear partial derivatives of the second order. This system of differential equations corresponds to the physical model, known in the literature as the Timoshenko Beam. The results obtained are applied to two particular cases representing respectively: the case of a beam with a rectangular section, with a constant height and a basis with a linear variation; and the case of a beam with a constant basis and a height with cubic variation.

Zainea, B.↗

Calculation of non-stationary aerodynamic forces in the near sonic range

The development of a mathematical method for calculating nonstationary supersonic flow in the near-sonic range is described. A perturbation formula is derived based on the exact stationary values; it is applicable to the equation for potential. The problem can thus be divided into stationary and nonstationary fields. The pressure distribution in an oscillating profile is determined, based on hyperbolic differential equations. It is shown that there are important corollaries concerning the application of linear theory. With suitable extrapolations, linear theory can be used up to about Mach 0.8. Linear theory is not applicable, however, when determining the moment coefficients; for this case, a special technique is described.

Teipel, I.↗

Open-mode delamination stress concentrations in horseshoe and elliptic composite curved bars subjected to end forces

The multilayer theory of anisotropic elasticity and a finite element method were used to analyze the open-mode delamination stress concentrations in horseshoe and elliptic laminated composite curved bars. Two types of laminations, solid laminations and sandwich laminations, were analyzed. It was found that the open-mode delamination stress concentration could be greatly increased in these two types of curved bars by decreasing their aspect ratios. The open-mode delamination stress concentration generated in the solid laminations was found to be far more severe than that generated in the sandwich laminations. The horseshoe curved bar may be used to determine both the open-mode delamination strength of solidly laminated composites and the open-mode debonding strength of sandwiched laminated composites. However, the elliptic curved bar is only good for determining the open-mode delamination strength of solidly laminated composites.

Ko, William L.↗

Evaluation of Load Analysis Methods for NASAs GIII Adaptive Compliant Trailing Edge Project

The Air Force Research Laboratory (AFRL), NASA Armstrong Flight Research Center (AFRC), and FlexSys Inc. (Ann Arbor, Michigan) have collaborated to flight test the Adaptive Compliant Trailing Edge (ACTE) flaps. These flaps were installed on a Gulfstream Aerospace Corporation (GAC) GIII aircraft and tested at AFRC at various deflection angles over a range of flight conditions. External aerodynamic and inertial load analyses were conducted with the intention to ensure that the change in wing loads due to the deployed ACTE flap did not overload the existing baseline GIII wing box structure. The objective of this paper was to substantiate the analysis tools used for predicting wing loads at AFRC. Computational fluid dynamics (CFD) models and distributed mass inertial models were developed for predicting the loads on the wing. The analysis tools included TRANAIR (full potential) and CMARC (panel) models. Aerodynamic pressure data from the analysis codes were validated against static pressure port data collected in-flight. Combined results from the CFD predictions and the inertial load analysis were used to predict the normal force, bending moment, and torque loads on the wing. Wing loads obtained from calibrated strain gages installed on the wing were used for substantiation of the load prediction tools. The load predictions exhibited good agreement compared to the flight load results obtained from calibrated strain gage measurements.

Aerodynamic Loads↗

Kinematics and control of a fully parallel force-reflecting hand controller for manipulator teleoperation

Design of a parallel force-reflecting hand controller that implements a friction- and inertia canceling control loop about the entire mechanism based on wrench sensing in the mechanism handgrip is discussed. Kinematics of the controller under consideration is analyzed and results are presented using a closed-form solution for the inverse kinematics and Newton-Raphson's method for the forward kinematics. Results indicate that the force control scheme based on a handgrip force sensor provides smaller steady-state errors than the scheme without a handigrip sensor.

Bryfogle, Mark D.↗

Nonlinear Aeroelastic Analysis Using a Time-Accurate Navier-Stokes Equations Solver

A method to simulate limit cycle oscillation (LCO) due to control surface freeplay using a modified CFL3D, a time-accurate Navier-Stokes computational fluid dynamics (CFD) analysis code with structural modeling capability, is presented. This approach can be used to analyze aeroelastic response of aircraft with structural behavior characterized by nonlinearity in the force verses displacement curve. A limited validation of the method, using very low Mach number experimental data for a three-degrees-of-freedom (pitch/plunge/flap deflection) airfoil model with flap freeplay, is also presented.

Kuruvila, Geojoe↗

Accurate force field for molybdenum by machine learning large materials data

In this work, we present a highly accurate spectral neighbor analysis potential (SNAP) model for molybdenum (Mo) developed through the rigorous application of machine learning techniques on large materials data sets. Despite Mo's importance as a structural metal, existing force fields for Mo based on the embedded atom and modified embedded atom methods still do not provide satisfactory accuracy on many properties. We will show that by fitting to the energies, forces and stress tensors of a large density functional theory (DFT)-computed dataset on a diverse set of Mo structures, a Mo SNAP model can be developed that achieves close to DFT accuracy in the prediction of a broad range of properties, including energies, forces, stresses, elastic constants, melting point, phonon spectra, surface energies, grain boundary energies, etc. We will outline a systematic model development process, which includes a rigorous approach to structural selection based on principal component analysis, as well as a differential evolution algorithm for optimizing the hyperparameters in the model fitting so that both the model error and the property prediction error can be simultaneously lowered. We expect that this newly developed Mo SNAP model will find broad applications in large-scale, long-time scale simulations.

36 MATERIALS SCIENCE↗

A Transient Response Method for Linear Coupled Substructures

A method is presented for determining the transient response of a discrete coordinate model of a linear structural system composed of substructures. The method is applicable to systems consisting of any number of substructures, both determinate and indeterminate interface boundaries, and any topological arrangement of the substructures. The method is simple to implement from a computational point of view because the equations of motion of each of the substructures are solved independently, and the interface boundary compatibility conditions are enforced at each integration time step by a matrix multiplication. The method is demonstrated for a structural system consisting of two beam segments and acted upon by a time dependent force. The numerical results from the demonstration problem validates the accuracy of the method. The application of this method to structural systems with changing interface boundary conditions between substructures is discussed.

Admire, J. R.↗

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intercomparison of observed cloud radiative forcing - A zonal and global perspective

Cloud radiative forcing determined by the Earth Radiation Budget Experiment (ERBE) scene identification method and two regression methods is compared using identical data sources and the same analysis periods. Results confirm that from the global perspective the increased reflection of SW radiation induced by clouds is greater than outgoing longwave energy trapped by clouds. The estimates of net effect obtained from the three methods are found to be in good agreement (within 3.5 W/sq m). It is concluded that differences in published global mean values of net cloud radiative forcing are mainly due to different data sources and analysis periods. The best estimated annual mean among all previous estimates is considered to be the ERBE measurement (-17.3 W/sq m).

Sohn, Byung-Ju↗