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,657 records · Page 92

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗

Excitation of Soliton-Type Waves in Crystals of the A{sub 3}B Stoichiometry

Using the method of molecular dynamics and taking Ni{sub 3}Al and Pt{sub 3}Al as examples, crystals of the A{sub 3}B composition are considered for the possibility of excitation of soliton-type waves in them. The potentials obtained by the embedded-atom method were used to describe interatomic interactions. It is shown that the harmonic external stimulus can excite waves of the soliton type in a Pt{sub 3}Al crystal, but not in Ni{sub 3}Al. Such compression–expansion waves are generated because of excitation of discrete breathers with soft nonlinearity that cannot exist in a Ni{sub 3}Al crystal near the affected region. The detected waves are capable of propagating to thousands of nanometers along the Pt{sub 3}Al crystal without losses of integrity and speed. The shape of the obtained wave corresponds to the kink solution of the sine-Gordon equation. The aggregate amount of energy transferred by a wave is determined by the number of rows of atoms involved in fluctuations; this may involve dozens and hundreds of electron volts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electrochemical Synthesis of Monodisperse Platinum-Cobalt Nanocrystals

The synthesis of platinum-cobalt nanocrystals (NCs) using colloidal solvothermal techniques is well understood. However, for monodisperse NCs to form, high temperatures and environmentally detrimental solvents are needed. We report a room temperature, aqueous method of platinum-cobalt NC synthesis using electrochemical reduction as the driving force for nucleation and growth. It is found that colloidal NCs will form in both the presence and absence of surfactant. Additionally, we report a monodisperse electrochemical deposition of NCs utilising a transparent conducting oxide electrode. The methods developed here will allow for a synthetic method to produce nanocatalysts with minimal environmental impact and should be readily applicable to other NC systems, including single- and multi-component alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimation of Directional Stability Derivatives at Small Angles and Subsonic and Supersonic Speeds

Methods are presented for estimating the directional stability derivative increments contributed by the stabilizing surfaces of subsonic and supersonic aircraft. These methods are strictly applicable at zero angle of attack and small angles of sideslip. The procedure of totaling the incremental coefficients to obtain an estimation of the total empennage side-force and yawing-moment coefficient derivatives is also shown, together with numerical examples. A correlation is presented between estimated and experimental incremental coefficients which indicates that the methods of this report generally estimate the increment of side force gained by the addition of a panel to within +/-10 percent of the experimental value while the yawing-moment increment is generally estimated to within +/-20 percent. This is true for both subsonic and supersonic Mach numbers. An example application of the methods to one of the problems in directional stability, that of minimizing the effect of Mach number on the side-force coefficient derivative of the empennage, is discussed.

Goodwin, Frederick K.↗

On-the-fly active learning of interpretable Bayesian force fields for atomistic rare events

Abstract Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.

36 MATERIALS SCIENCE↗

Performance Assessment of Fluidic Oscillators Tested on the NASA Hump Model

Flow separation control over a wall-mounted hump model was studied experimentally to assess the performance of fluidic oscillators (sweeping jet actuators). An array of fluidic oscillators was used to control flow separation. The results showed that the fluidic oscillators were able to achieve substantial control over the separated flow by increasing the upstream suction pressure and downstream pressure recovery. Using the data available in the literature, the performance of the fluidic oscillators was compared to other active flow control (AFC) methods such as steady blowing, steady suction, and zero-net-mass-flux (ZNMF) actuators. Several integral parameters, such as the inviscid flow comparison coefficient, pressure drag coefficient, and modified normal force coefficient, were used as quality metrics in the performance comparison of the AFC methods. These quality metrics indicated the superiority of the steady suction method, especially at lower excitation amplitudes that is followed by the fluidic oscillators, steady blowing, and the ZNMF actuators, respectively. An aerodynamic figure of merit (AFM) was also constructed using the integral parameters and AFC power usage. The AFM results revealed that, for this study, steady suction was the most efficient AFC method at lower excitation amplitudes. The steady suction loses its efficiency as the excitation amplitude increases, and the fluidic oscillators become the most efficient AFC method. Both the steady suction and the fluidic oscillators have an AFM > 1 for the range tested in this study, indicating that they provide a net benefit when the AFC power consumption is also considered. On the other hand, both the steady blowing and ZNMF actuators were found to be inefficient AFC methods (AFM < 1) for the current configuration. Although they improved the flow field by controlling flow separation, the power requirement was more than their benefit.

Fluidic oscillators↗

Embedded function methods for compressible high speed turbulent flow

Fundamental issues relating to compressible turbulent flow are addressed. The focus has been on developing methods and testing concepts for attached flows rather than trying to force a conventional law of the wall into a zone of backflow. Although the dynamics of the near-wall flow in an attached turbulent boundary layer are relatively well documented, the dynamical features of a zone of reversed turbulent flow are not, nor are they well understood. Incompressibility introduces effects and issues that have been dealt with only marginally in the literature, therefore, the present work has been focussed on attached high-speed flows. The wall function method has been extended up through the supersonic to hypersonic speeds. Algorithms have been successfully introduced into the code that calculates the flow all the way to the wall, and testing is being carried out for progressively more complex flow situations.

Walker, J. D. A.↗

Rapid Freeform Sheet Metal Forming: Technology Development and System Verification

The objective of this project is to develop a transformational RApid Freeform sheet metal Forming Technology (RAFFT) in an industrial environment, which has the potential to increase manufacturing energy efficiency up to ten times, at a fraction of the cost of conventional technologies. The RAFFT technology is a flexible and energy-efficient process that eliminates the need for having geometry-specific forming dies. The innovation lies in the idea of using the energy resource at the local deformation area which provides greater formability, process control, and process flexibility relative to traditional methods. Double-Sided Incremental Forming (DSIF), the core technology in RAFFT, is a new concept for sheet metal forming. A blank sheet is clamped around its periphery and gradually deformed into a complex 3D freeform part by two strategically aligned stylus-type tools that follow a pre-described toolpath. The two tools, one on each side of the blank, can form a part with sharp features for both concave and convex shapes. Since deformation happens locally, the forming force at any instant is significantly decreased when compared to traditional methods. The key advantages of DSIF are its high process flexibility, high energy-efficiency, low capital investment, and the elimination of the need for massive amounts of die casting and machining. Additionally, the enhanced formability and process flexibility of DSIF can open up design spaces and result in greater weight savings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Computer-aided roll pass design in rolling of airfoil shapes

This paper describes two computer-aided design (CAD) programs developed for modeling the shape rolling process for airfoil sections. The first program, SHPROL, uses a modular upper-bound method of analysis and predicts the lateral spread, elongation, and roll torque. The second program, ROLPAS, predicts the stresses, roll separating force, the roll torque and the details of metal flow by simulating the rolling process, using the slab method of analysis. ROLPAS is an interactive program; it offers graphic display capabilities and allows the user to interact with the computer via a keyboard, CRT, and a light pen. The accuracy of the computerized models was evaluated by (a) rolling a selected airfoil shape at room temperature from 1018 steel and isothermally at high temperature from Ti-6Al-4V, and (b) comparing the experimental results with computer predictions. The comparisons indicated that the CAD systems, described here, are useful for practical engineering purposes and can be utilized in roll pass design and analysis for airfoil and similar shapes.

Akgerman, N.↗

Dynamic analysis of nonlinear rotor-housing systems

Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine (SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the nonlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increment. The method is applied to a nonlinear generic model of the high pressure oxygen turbopump (HPOTP). As compared to the fourth order Runge-Kutta numerical integration methods, the convolution approach proved to be more accurate and more highly efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic responses fo the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-rotor models to their coordinates at the bearing clearances. Recommendations are included for further development of the method, for extending the analysis to aperiodic and chaotic regimes and for conducting critical parameteric studies of the nonlinear response of the current SSME turbopumps.

Noah, Sherif T.↗

Evaluation of Methods for In-Situ Calibration of Field-Deployable Microphone Phased Arrays

Current field-deployable microphone phased arrays for aeroacoustic flight testing require the placement of hundreds of individual sensors over a large area. Depending on the duration of the test campaign, the microphones may be required to stay deployed at the testing site for weeks or even months. This presents a challenge in regards to tracking the response (i.e., sensitivity) of the individual sensors as a function of time in order to evaluate the health of the array. To address this challenge, two different methods for in-situ tracking of microphone responses are described. The first relies on the use of an aerial sound source attached as a payload on a hovering small Unmanned Aerial System (sUAS) vehicle. The second relies on the use of individually excited ground-based sound sources strategically placed throughout the array pattern. Testing of the two methods was performed in microphone array deployments conducted at Fort A.P. Hill in 2015 and at Edwards Air Force Base in 2016. The results indicate that the drift in individual sensor responses can be tracked reasonably well using both methods. Thus, in-situ response tracking methods are useful as a diagnostic tool for monitoring the health of a phased array during long duration deployments.

Humphreys, William M.↗

Drop Ejection From an Oscillating Rod

The dynamics of a drop of a Newtonian liquid that is pendant from or sessile on a solid rod that is forced to undergo time-periodic oscillations along its axis is studied theoretically. The free boundary problem governing the time evolution of the shape of the drop and the flow field inside it is solved by a method of lines using a finite element algorithm incorporating an adaptive mesh. When the forcing amplitude is small, the drop approaches a limit cycle at large times and undergoes steady oscillations thereafter. However, drop breakup is the consequence if the forcing amplitude exceeds a critical value. Over a wide range of amplitudes above this critical value, drop ejection from the rod occurs during the second oscillation period from the commencement of rod motion. Remarkably, the shape of the interface at breakup and the volume of the primary drop formed are insensitive to changes in forcing amplitude. The interface shape at times close to and at breakup is a multi-valued function of distance measured along the rod axis and hence cannot be described by recently popularized one-dimensional approximations. The computations show that drop ejection occurs without the formation of a long neck. Therefore, this method of drop formation holds promise of preventing formation of undesirable satellite droplets.

Wilkes, E. D.↗

Finite element analyses of single particle crushing tests incorporating computed tomography imaging and damage mechanics

Fracture of individual grains within a granular assembly affect the strength and deformation behavior of granular materials under large stress states. Interparticle forces, leading to contact stresses and ultimately fracture initiation, are influenced by particle morphology. A numerical method that incorporates both particle fracture and morphology can provide a more accurate model of a granular material’s macro-scale response. This research investigates an approach to modeling fracture initiation and propagation within individual grains, utilizing high resolution X-ray computed tomography (CT) imaging to consider grain morphology and an explicit finite element code optimized for high performance computing which incorporates damage mechanics. Using precisely measured force-displacement curves from single particle crushing tests of manufactured quartz spheres, the method is validated, and the effects of fracture properties on grain fragmentation are determined. This approach is then used to simulate single particle crushing of Ottawa sand, with grain morphology incorporated through discretized CT images. Here the effect of grain orientation on the force and applied displacement at catastrophic splitting is also explored. The proposed discrete particle based finite element framework can be applied to modeling granular assemblies, considering the effect of individual particle shape and fracture on the assembly’s deformation response through well calibrated numerical simulations.

42 ENGINEERING↗

ISS Payload Operations Training During the COVID-19 Pandemic: Impacts and Solutions

Introduction: The onset and the protracted duration of the COVID-19 pandemic brought a dramatic and rapid transformation to almost every aspect of humanity in 2020. The world’s space agencies and their missions were not immune to the unpredictable and wide-sweeping changes. One discipline principally affected was mission operations and the various groups supporting that critical function. Mission support teams, especially for complex and crewed missions like the International Space Station (ISS) were forced to rapidly rethink how and where control center staff performed their vital work. Operations training – an essential element to mission ops, had unique hurdles to overcome. Operations training is responsible for preparing astronaut crews for their scientific missions, providing recurring proficiency and currency material to flight controllers, as well as ensuring new team members are ready to join their colleagues on console. Every element of training was impacted by COVID restrictions. From orientation and introductory classes for new controllers, training for scientists and payload developers, to simulations and advanced lessons, critical On the Job Training (OJT) and final evaluations; all aspects faced new challenges. Operations trainers at NASA’s Marshall Space Flight Center (MSFC) in Huntsville, Alabama, were forced to rethink, retool, and become more efficient with trainees and training resources in order to continue supporting payload operations on the ISS. Impacts and Barriers to Operations Training: The pandemic started in the United States in early March 2020. Immediately, NASA mandated that the support for ISS real-time operations was critical. As a result, physical access to the ISS Payload Operations Center was restricted to only those people supporting console operations. Thus, use of on-site NASA facilities and any trainer to trainee face-to-face interaction was not authorized. Trainers and trainees had to quickly shift to 100% remote learning. In the short term, this was not a problem since NASA had already prepared Information Technology (IT) resources that facilitated remote access to mission resources like classes, documents and teleconferences. However, instructors quickly discovered the classes that they were accustomed to delivering in a large classroom environment did not translate well to remote teaching. The lack of non-verbal communication with the class became a major barrier to instruction. Another significant hurdle to operations training was the mandate that all simulations could only be held remotely. The logistics of planning and executing even small simulations proved to be challenging, almost impossible, due to IT technology and public internet limitations. With simulations essentially stopped, as well as the restrictions on most OJT training, trainees were essentially stopped in their advancement towards certification. An additional barrier that had to be overcome was the formal documentation of the curriculum objectives. Even though NASA had been deliberately transitioning to electronic documents for years, some organizations were still relying on “pen and ink” documents to formalize training completion. The logistics of implementing secure and controlled e-signature documents in a standardized, NASA approved format had to be addressed. Lesson Learned and Efficiencies: Once the various limitations were identified, ISS trainers at MSFC discussed, proposed, and prioritized new options and techniques. Transitioning to all electronic learning materials and certifications was a relevantly easy fix. Teaching courses to large and diverse groups took additional shifts in the NASA training paradigm. Likewise, methods and procedures for preparing astronauts for their missions were forced to be revised. Simulation supervisors and their teams were also forced to find more efficient techniques to provide realistic training experience. Communication and approvals from management was essential. In every case, the ISS payload operations instructors found novel and efficient solutions to all functions listed above. This paper will discuss the factors and solutions payloads operations trainers found to keep scientific research on the ISS flying forward to mission success.

ISS Operations↗

Future climate emulations using quantile regressions on large ensembles

The study of climate change and its impacts depends on generating projections of future temperature and other climate variables. For detailed studies, these projections usually require some combination of numerical simulation and observations, given that simulations of even the current climate do not perfectly reproduce local conditions. We present a methodology for generating future climate projections that takes advantage of the emergence of climate model ensembles, whose large amounts of data allow for detailed modeling of the probability distribution of temperature or other climate variables. The procedure gives us estimated changes in model distributions that are then applied to observations to yield projections that preserve the spatiotemporal dependence in the observations. We use quantile regression to estimate a discrete set of quantiles of daily temperature as a function of seasonality and long-term change, with smooth spline functions of season, long-term trends, and their interactions used as basis functions for the quantile regression. A particular innovation is that more extreme quantiles are modeled as exceedances above less extreme quantiles in a nested fashion, so that the complexity of the model for exceedances decreases the further out into the tail of the distribution one goes. We apply this method to two large ensembles of model runs using the same forcing scenario, both based on versions of the Community Earth System Model (CESM), run at different resolutions. The approach generates observation-based future simulations with no processing or modeling of the observed climate needed other than a simple linear rescaling. The resulting quantile maps illuminate substantial differences between the climate model ensembles, including differences in warming in the Pacific Northwest that are particularly large in the lower quantiles during winter. We show how the availability of two ensembles allows the efficacy of the method to be tested with a “perfect model” approach, in which we estimate transformations using the lower-resolution ensemble and then apply the estimated transformations to single runs from the high-resolution ensemble. Finally, we describe and implement a simple method for adjusting a transformation estimated from a large ensemble of one climate model using only a single run of a second, but hopefully more realistic, climate model.

54 ENVIRONMENTAL SCIENCES↗

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

Efficient Integration of Algebraic Constraint for Exponential Time Integration

This work is a continuation of a previous project where the high efficiency of exponential integration methods for DRMHD systems was demonstrated. Often algebraic constraints must also be enforced for DRMHD systems of interest. Straightforward application of exponential methods to DRMHD equations with constraints leads to prohibitively computationally expensive methods. In this work, we propose new exponential schemes that allow the constraints to be removed from evaluation of exponential matrix functions which drastically reduces computational complexity by eliminating the need to perform computations enforcing constraints in exponential calculations while still preserving a high order of accuracy and allowing for a large time step even when the problem is stiff. This idea is similar to the W-methods and is achieved by carefully designing a method with the desired order of accuracy, even with an incomplete Jacobian matrix used as an argument of exponential-like functions. The constraints are accounted for by including them in the evaluation of the right-hand-side forcing function of the spatially discretized system. We study performance of the new methods on test problems and outline future research directions that this work opens.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗