Search NASA⌕ Search

SEARCH · Search NASA

Results for “force vector”

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 325 records · Page 18

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

Computation of solar magnetic fields from photospheric observations

The observational difficulties of obtaining the magnetic field distribution in the chromosphere and corona of the sun has led to methods of extending photospheric magnetic mesurements into the solar atmosphere by mathematical procedures. A new approach to this problem presented here is that a constant alpha force-free field can be uniquely determined from the tangential components of the measured photospheric flux alone. The vector magnetographs now provide measurements of both the solar photospheric tangential and the longitudinal magnetic field. This paper presents derivations for the computation of the solar magnetic field from these type of measurements. The fields considered are assumed to be a constant alpha force-free fields or equivalent, producing vanishing Lorentz forces. Consequently, magnetic field lines and currents are related by a constant and hence show an identical distribution. The magnetic field above simple solar regions are described from the solution of the field equations.

Hannakam, L.↗

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗

Thrust-induced effects on low-speed aerodynamics of fighter aircraft

Results of NASA Langley has conducted wind-tunnel investigations of several fighter configurations conducted to determine the effects of both thrust vectoring and spanwise blowing are reviewed. A recent joint NASA/Grumman Aerospace Corporation/U.S. Air Force Wright Aeronautical Laboratory wind-tunnel investigation was conducted to examine the effects of spanwise blowing on the trailing-edge flap system. This application contrasts with the more familiar method of spanwise blowing near the wing leading edge. Another joint program among NASA/McDonnell Aircraft Company/U.S. Air Force Wright Aeronautical Laboratory investigated the effects of reverse thrust on the low-speed aerodynamics of an F-15 configuration. The F-15 model was fitted with a rotating van thrust reverser concept which could simulate both in-flight reversing for approach and landing or full reversing for ground roll reduction. The significant results of these two joint programs are reported.

Banks, D. W.↗

Implementation and Verification of Cable Bending Stiffness in MoorDyn: Preprint

Cable bending stiffness capabilities have been added to MoorDyn to enable modeling of dynamic power cables. The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. MoorDyn is a lumped-mass mooring dynamics model that is used for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim). Its existing line elasticity formulation considers axial stiffness only. To properly capture the dynamics of power cables, a cable bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent line elements. The resulting bending moment is realized by applying forces on adjacent nodes, leaving the underlying lumped-mass formulation unchanged. The bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗

Controller Design for EMA in TVC Incorporating Force Feedback

The objective of this research was to develop control schemes and control design procedures for electromechanical actuators (EMA) in thrust vector control (TVC) applications. For a variety of reasons, there is a tendency within the aerospace community to use electromechanical actuators in applications where hydraulics have traditionally been employed. TVC of rocket engines is one such application. However, there is considerable research, development, and testing to be done before EMA will be accepted by the community at large for these types of applications. Besides the development of design procedures for the basic position controller, two major concerns are dealt with in this research by incorporating force feedback: 1) the effects of resonance on the performance of EMA-TVC-rocket-engine systems, and 2) the effects of engine start transients on EMA. This report only highlights the major contributions of this research.

Schinstock, Dale E.↗

Warm gas TVC design study

A warm gas thrust vector control system was studied to optimize the injection geometry for a specific engine configuration, and an injection valve was designed capable of meeting the base line requirements. To optimize injection geometry, studies were made to determine the performance effects of varying injection location, angle, port size, and port configuration. Having minimized the injection flow rate required, a warm gas valve was designed to handle the required flow. A direct drive hydraulic servovalve capable of operating with highly contaminated hydraulic fluid was designed. The valve is sized to flow 15 gpm at 3000 psia and the direct drive feature is capable of applying a spool force of 200 pounds. The baseline requirements are the development of 6 deg of thrust vector control utilizing 2000 F (total temperature) gas for 180 seconds on a 1.37 million pound thrust engine burning LOX and RP-1 at a chamber pressure of 250 psia with a 155 inch long conical nozzle having a 68 inch diameter throat and a 153 inch diameter exit.

Moorhead, S. B., Jr.↗

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗

Testing of turbulent seals for rotordynamic coefficients

A facility has been developed for dynamic testing of straight and convergent tapered seals with the capability of measuring the radial and tangential force components which result from a circular centered orbit. The test apparatus causes the seal journal to execute small-eccentricity centered circular orbits within its clearance circle. Dynamic measurements are made and recorded of the seal displacement-vector components, and of the pressure field. The pressure field is integrated to yield seal tangential and radial reaction-force components. Representative test data are provided and discussed for straight seals.

Childs, D. W.↗

Kinematics and aerodynamics of the velocity vector roll

The velocity vector roll is an angular rotation of an aircraft about its instantaneous velocity vector, constrained to be performed at constant angle-of-attack (AOA), no sideslip, and constant velocity. Consideration of the aerodynamic force equations leads to requirements for body-axis yawing and pitching rotations that satisfy these constraints. Here, the body axis rotations, and the constraints, are used in the moment equations to determine the aerodynamic moments required to perform the velocity vector roll. For representative tactical aircraft, the conditions for maximum pitching moment are a function of orientation, occurring at about 90 deg of bank in a level trajectory. Maximum required pitching moment occurs at peak roll rate, and is achieved at AOA above 45 deg. The conditions for maximum rolling moment depend on the value of the roll mode time constant. For a small time constant (fast response) the maximum rolling moment occurs at maximum roll acceleration and zero AOA, largely independent of aircraft orientation; for a large time constant, maximum required rolling moment occurs at maximum roll rate, at maximum AOA, and at 180 deg of bank in level flight. Maximum yawing moment occurs at maximum roll acceleration, maximum AOA, and is largely independent of airplane orientation.

Durham, Wayne C.↗

Spherical balloon response to three-dimensional time-dependent flows

The concept of the Lagrangian displacement of a balloon is introduced. It is shown that the general balloon response problem is extremely complicated because the wind-forcing functions in the balloon equations of motion are functions of the wind velocity vector and its Eulerian first derivatives evaluated at the location of the balloon. The linear perturbation equations for a spherical balloon are derived by perturbing the components of velocity of the balloon about a terminal velocity state which is in equilibrium with a space-time invariant mean horizontal flow. The atmospheric flow is also perturbed such that the resulting equations can be used to analyze the responses of spherical balloons to three-dimensional time-dependent flows. The wind field is represented in terms of a four-fold Fourier integral that involves three orthogonal wave numbers and a frequency, while the balloon components of velocity are represented as Fourier integrals involving a frequency which, in turn, is a function of the wind field wave numbers and frequency and the unperturbed flow components of velocity.

Fichtl, G. H.↗

Finite amplitude Alfven waves in a multi-ion plasma - Propagation, acceleration, and heating

An expression is derived for the wave action flux of finite-amplitude Alfven waves in a multi-ion plasma. The expression is valid in the presence of dissipative forces and permits an arbitrary angle between the average magnetic field and the wave vector. Applying the conservation of wave action and the first law of thermodynamics yields, for a multi-ion plasma, an expression for the spatial evolution of Alfven wave amplitude in the absence of dissipation. It also gives the relationship between the wave amplitude and the dissipative heating, as well as an expression for the acceleration of an ion species by finite-amplitude Alfven waves. It is pointed out that the acceleration comprises a nondissipative wave pressure that is identical to that derived previously under more restrictive conditions and a new term giving the acceleration that must accompany dissipative heating. The results are discussed in the context of the observations of heavy ions in the solar wind.

Isenberg, P. A.↗

The thermosphere as a sink of magnetospheric energy - A review of recent observations of dynamics

It is pointed out that the past few years have seen an unprecedented influx of new experimental information on the dynamics of the neutral upper atmosphere of the earth. Vector wind measurements provide new information for studies of the thermospheric response to magnetospheric forcing. This response occurs through the medium of convecting ionospheric ions set into motion by electric fields of magnetospheric origin. The ultimate sink for much of the energy and momentum coming from the magnetosphere is the neutral thermosphere whose dynamics have, in the past, received far less attention than their ionospheric counterpart because of basic experimental limitations. In this paper, a review is provided of the progress made in the last few years on the basis of the Dynamics Explorer neutral wind observations, taking into account the coupling between the magnetosphere and the thermosphere via the ionosphere.

Killeen, T. L.↗

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

Integral Velocity Sampling for Unsteady Rotor Models on Cartesian Meshes

An unsteady body-force rotor model utilizing integral velocity sampling is developed for Cartesian meshes with a focus on accurate performance predictions. In this approach, the rotors are modeled as body-force source terms, and integral sampling is used to define the local freestream velocity vector. The method is verified in two dimensions for airfoils, providing accurate angle of attack and freestream Mach number estimates even on relatively coarse meshes, and demonstrating excellent mesh convergence across a variety of source region thicknesses. Extensions of the model to three dimensions for both fixed and rotating wings are considered, with numerical experiments demonstrating the need for a spanwise-continuous representation of the induced velocity or angle of attack. The model provides reasonable spanwise lift distributions without any explicit tip-loss model. Comparisons to high-fidelity simulations and experimental data show excellent agreement for predictions of an isolated DJI Phantom 3 propeller in hover. Finally, simulations of the full aircraft in edgewise forward flight also compare well with wind tunnel data across a variety of pitch angles and propeller speeds.

ARMD↗

Quaternion normalization in additive EKF for spacecraft attitude determination

This work introduces, examines, and compares several quaternion normalization algorithms, which are shown to be an effective stage in the application of the additive extended Kalman filter (EKF) to spacecraft attitude determination, which is based on vector measurements. Two new normalization schemes are introduced. They are compared with one another and with the known brute force normalization scheme, and their efficiency is examined. Simulated satellite data are used to demonstrate the performance of all three schemes. A fourth scheme is suggested for future research. Although the schemes were tested for spacecraft attitude determination, the conclusions are general and hold for attitude determination of any three dimensional body when based on vector measurements, and use an additive EKF for estimation, and the quaternion for specifying the attitude.

Bar-Itzhack, I. Y.↗