Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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
NASA NTRS2010

Modeling of Failure for Analysis of Triaxial Braided Carbon Fiber Composites

In the development of advanced aircraft-engine fan cases and containment systems, composite materials are beginning to be used due to their low weight and high strength. The design of these structures must include the capability of withstanding impact loads from a released fan blade. Relatively complex triaxially braided fiber architectures have been found to yield the best performance for the fan cases. To properly work with and design these structures, robust analytical tools are required that can be used in the design process. A new analytical approach models triaxially braided carbon fiber composite materials within the environment of a transient dynamic finite-element code, specifically the commercially available transient dynamic finite-element code LS-DYNA. The geometry of the braided composites is approximated by a series of parallel laminated composites. The composite is modeled by using shell finite elements. The material property data are computed by examining test data from static tests on braided composites, where optical strain measurement techniques are used to examine the local strain variations within the material. These local strain data from the braided composite tests are used along with a judicious application of composite micromechanics- based methods to compute the stiffness properties of an equivalent unidirectional laminated composite required for the shell elements. The local strain data from the braided composite tests are also applied to back out strength and failure properties of the equivalent unidirectional composite. The properties utilized are geared towards the application of a continuum damage mechanics-based composite constitutive model available within LS-DYNA. The developed model can be applied to conduct impact simulations of structures composed of triaxially braided composites. The advantage of this technology is that it facilitates the analysis of the deformation and damage response of a triaxially braided polymer matrix composite within the environment of a transient dynamic finite-element code such as LS-DYNA in a manner which accounts for the local physical mechanisms but is still computationally efficient. This methodology is tightly coupled to experimental tests on the braided composite, which ensures that the material properties have physical significance. Aerospace or automotive companies interested in using triaxially braided composites in their structures, particularly for impact or crash applications, would find the technology useful. By the development of improved design tools, the amount of very expensive impact testing that will need to be performed can be significantly reduced.

Goldberg, Robert K.↗
DOE OSTI2023

Magnetic frustration driven by conduction carrier blocking in Nd 2 Co 0.85 Si 2.88

The intermetallic compound Nd 2 Co 0.85 Si 2.88 having a triangular lattice could be synthesized in single phase only with defect crystal structure. Investigation through different experimental techniques indicate the presence of two magnetic transitions in the system. As verified experimentally and theoretically, the high-temperature transition $T_H$ ~ 140 K is associated with the development of ferromagnetic interaction between itinerant Co moments, whereas the low-temperature transition at $T_L$ ~ 6.5 K is due to the coupling among Nd–4$\textit{f}$ and Co–3$\textit{d}$ moments, which is antiferromagnetic in nature. Detailed studies of temperature-dependent dc magnetic susceptibility, field dependence of isothermal magnetization, nonequilibrium dynamical behavior, viz., magnetic relaxation, aging effect, magnetic-memory effect, and temperature dependence of heat capacity, along with density functional theory (DFT) calculations, suggest that the ground state is magnetically frustrated spin glass in nature, having competing magnetic interactions of equivalent energies. DFT results further reveal that the 3$\textit{d}$/5$\textit{d}$-conduction carriers are blocked in the system and act as a barrier for the 4$\textit{f}$–4$\textit{f}$ RKKY interactions, resulting in spin frustration. Presence of vacancy defects in the crystal are also conducive to the spin frustration. This is an unique mechanism of magnetic frustration, not emphasized so far in any of the ternary $R_2TX_3$ ($\textit{R}$ = rare earth, $\textit{T}$ = transition elements, and $\textit{X}$ = Si, Ge, In) type compounds. Due to the competing character of the itinerant 3$\textit{d}$ and localized 4$\textit{f}$ moments, the compound exhibits anomalous field dependence of magnetic coercivity. Further, the system also exhibits a considerable magnetic entropy change of –Δ$S_M$ ~ 13.3 J/kg K with a relative cooling power (RCP) of 220 J/kg and adiabatic temperature change $ΔT_{\text{ad}}$ of 6 K for magnetic field change of 70 kOe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗
DOE OSTI2020

Design and analysis of a wake steering controller with wind direction variability

Wind farm control strategies are being developed to mitigate wake losses in wind farms, increasing energy production. Wake steering is a type of wind farm control in which a wind turbine's yaw position is misaligned from the wind direction, causing its wake to deflect away from downstream turbines. Current modeling tools used to optimize and estimate energy gains from wake steering are designed to represent wakes for fixed wind directions. However, wake steering controllers must operate in dynamic wind conditions and a turbine's yaw position cannot perfectly track changing wind directions. Research has been conducted on robust wake steering control optimized for variable wind directions. In this paper, the design and analysis of a wake steering controller with wind direction variability is presented for a two-turbine array using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented wake model. First, the authors propose a method for modeling the turbulent and low-frequency components of the wind direction, where the slowly varying wind direction serves as the relevant input to the wake model. Next, we explain a procedure for finding optimal yaw offsets for dynamic wind conditions considering both wind direction and yaw position uncertainty. We then performed simulations with the optimal yaw offsets applied using a realistic yaw offset controller in conjunction with a baseline yaw controller, showing good agreement with the predicted energy gain using the probabilistic model. Using the Gaussian wake model in FLORIS as an example, we compared the performance of yaw offset controllers optimized for static and dynamic wind conditions for different turbine spacings and turbulence intensity values, assuming uniformly distributed wind directions. For a spacing of five rotor diameters and a turbulence intensity of 10 %, robust yaw offsets optimized for variable wind directions yielded an energy gain equivalent to 3.24 % of wake losses recovered, compared to 1.42 % of wake losses recovered with yaw offsets optimized for static wind directions. In general, accounting for wind direction variability in the yaw offset optimization process was found to improve energy production more as the separation distance increased, whereas the relative improvement remained roughly the same for the range of turbulence intensity values considered.

17 WIND ENERGY↗
NASA NTRS1979

Parameter identification applied to analytic hingeless rotor modeling

It is known that dynamic rotor inflow has a substantial effect on rotor dynamic loads. Despite the complexity of the unsteady flow problem, simple analytical models can be made useful by identifying their parameters from transient response tests without performing flow measurements. Two analytical inflow models are studied: the first is based on an equivalent blade Lock number, the second is based on a time delayed unsteady momentum inflow. In preparation for the experimental data analysis, identifications from simulated test data and an eigenvalue analysis are performed. The experimental results show that the first analytical inflow model is accurate for rotor advance ratios of 0.4 and above. For lower advance ratios, the second inflow model provides better accuracy. Prediction studies with experimental data not used for the identification are performed to determine the accuracy of the mathematical models.

Banerjee, D.↗
NASA NTRS1980

Transmission line model of the interaction of a long metal wire with the ionosphere

An equivalent electric circuit model is used to study the electrodynamic interactions of long orbiting metallic tethers with the ionospheric plasma and, in particular, to derive current and potential profiles along bare metallic tethers. In contrast with other models, this approach is dynamic, enabling both the transient behavior of the wire and its final equilibrium state to be derived. A comparison with the results of other models indicates the advantage of the present approach, especially in those cases where the internal resistance of the tether plays a major role in determining the current and potential distributions.

Arnold, D. A.↗
NASA NTRS1991

Modeling of the slewing control of a flexible structure

This paper presents a formulation for the modeling of a single-link flexible beam that is undergoing slewing motion at an actively controlled pinned end, where the other end of the beam is free. This pinned end, or slewing axis, is of fixed orientation such that the beam rotates in the horizontal plane. A geared dc electric motor is connected to the beam at the slewing axis. A position and a velocity sensor are placed at this hinged location, and their proportional feedback provides a position control system about the slewing axis. The motor characteristics, gear ratio, and the position feedback constant determine an equivalent rotational spring constant, often called the servo stiffness. This paper generalizes the structure's boundary conditions at the slewing axis to include the effects of the servo system. It is shown that the clamped-free beam assumption for the dynamics of the structure is a valid assumption if the ratio of the servo stiffness to beam flexibility is high. However, for moderate or low values of this ratio, the clamped-free beam leads to erroneous system models so that it becomes necessary to consider the effects of the driving servo on the dynamics of the flexible beam.

Garcia, Ephrahim↗
NASA NTRS2010

Modification of a Macromechanical Finite-Element Based Model for Impact Analysis of Triaxially-Braided Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. For the current analytical approach, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The constitutive model requires stiffness and strength properties of an equivalent unidirectional composite. Simplified micromechanics methods are used to determine the equivalent stiffness properties, and results from coupon level tests on the braided composite are utilized to back out the required strength properties. Simulations of quasi-static coupon tests of several representative braided composites are conducted to demonstrate the correlation of the model. Impact simulations of a represented braided composites are conducted to demonstrate the capability of the model to predict the penetration velocity and damage patterns obtained experimentally.

Goldberg, Robert K.↗
NASA NTRS2018

Data Assimilation of Terrestrial Water Storage to Adjust Precipitation Fluxes

The Gravity Recovery and Climate Experiment (GRACE) mission has provided unprecedented observations of terrestrial water storage (TWS) dynamics at basin to continental scales. TWS is defined as the sum of groundwater, soil moisture, snow, surface water, ice and biomass water. Data assimilation of GRACE TWS observations has been shown to improve simulation of groundwater, streamflow, and snow water equivalent, and has also proven useful for drought monitoring and identifying human impacts on the water cycle. From a modeling perspective, the TWS components are defined as "prognostic hydrological states". Existing GRACE data assimilation schemes update these prognostic states directly. In this work, we propose an alternate approach in which precipitation fluxes are adjusted in order to achieve the desired change in the hydrological prognostic states. Limitations of such an approach include the assumption that all errors in TWS originate from errors in precipitation. Nonetheless, benefits comprise (1) the water balance is maintained, as opposed to having to add increments to the water budget components, (2) the model automatically determines how to distribute the updates among the TWS prognostic states, and (3) it is not necessary to know the exact time of the observation TWS, because the TWS change timing is determined by the precipitation forcing.

Girotto, Manuela↗
DOE OSTI2021

Dynamic process intensification: Fundamentals and implementation to ternary distillation

Distillation remains a key technology for separating liquid mixtures. Its versatility comes with the disadvantage of high energy consumption. We previously used empirical arguments to introduce dynamic process intensification (DPI) as a strategy for improving the energy efficiency of binary distillation. In this article, we focus on ternary distillation; we begin by providing a rigorous basis for DPI, then formulate the problem of identifying the operating states for DPI as a nonlinear optimization problem. Via an extensive case study considering a hydrocarbon mixture, we demonstrate that reboiler energy use can be reduced by more than 2.3% relative to an equivalent column operating at steady state. Here, we prove that this result is due to the fact that DPI aims to satisfy product flow and quality constraints on average in time, affording broader opportunities for optimizing column economics than steady-state operation, where constraints are met strictly at all times.

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

Power System Network Reduction for Power Hardware-in-the-Loop Simulation: Preprint

This paper proposes a single-port equivalent and a two-port equivalent to respectively reduce single-port and two-port areas in a large power network. Parameters of the reduced systems are rigorously derived, which guarantees that the electrical quantities at the port remain unchanged over the reduction, including voltage magnitude and phase, active and reactive power injections into the area to be reduced. The proposed techniques are applied to reduce a practical Maui grid, where the total numbers of buses, lines and transformers are respectively reduced from 212, 106 and 108 to 45, 30 and 13. Dynamic behaviors between the full model and the reduced model are compared in detail to illustrate the efficacy and accuracy of the proposed network reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2024

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2024

Greybox Thermal Parameter Identification of Electric Machine Stators

The parameters of electric machine thermal equivalent circuit networks are difficult to predict due to material and manufacturing uncertainties. In this paper, a Greybox system identification approach is used to identify parameters of electric machine stator lumped parameter thermal networks (LPTNs). LPTNs provide a low order, computationally efficient, dynamic model of temperatures at specific locations. Second and third order LPTN model structures are defined as state space equations with stator thermal parameters to be identified. To test the Greybox electric machine stator thermal system identification, five stator motorette prototypes were constructed with controlled variations in slot fill and slot liner thickness. The variation in the motorette thermal parameters and thermal time constants are detected using the Greybox identification. Special attention is given to the impact of sampling rate and Greybox data record length on parameter estimation accuracy.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1990

Modeling of vortex-induced oscillations based on indicial response approach

The indicial response approach, a modeling approach originally used for studying nonlinear problems in flight dynamics is applied to the study of vortex-induced oscillations phenomena. The indicial response of the velocity field is derived for the problem studied with emphasis on physical postulates involved. A full account of fluid dynamics effects is taken by considering the incompressible Navier-Stokes equations. The theory is applied to the particular case of flow past a cylinder in periodically forced motion to derive some salient features. The indicial response approach is shown to be equivalent to a currently popular approach based on the use of the amplitude equation. Jump and hysteresis phenomena that experiments indicate occur within the lock-in regime (where the velocity field oscillates at the forcing frequency) are captured by the approach.

Truong, Khiem Van↗
NASA NTRS1988

Applications of Digital Image Processing 11

A new technique, digital image velocimetry, is proposed for the measurement of instantaneous velocity fields of time dependent flows. A time sequence of single-exposure images of seed particles are captured with a high-speed camera, and a finite number of the single-exposure images are sampled within a prescribed period in time. The sampled images are then digitized on an image processor, enhanced, and superimposed to construct an image which is equivalent to a multiple exposure image used in both laser speckle velocimetry and particle image velocimetry. The superimposed image and a single-exposure Image are digitally Fourier transformed for extraction of information on the velocity field. A great enhancement of the dynamic range of the velocity measurement is accomplished through the new technique by manipulating the Fourier transform of both the single-exposure image and the superimposed image. Also the direction of the velocity vector is unequivocally determined. With the use of a high-speed video camera, the whole process from image acquisition to velocity determination can be carried out electronically; thus this technique can be developed into a real-time capability.

Cho, Y. -C.↗
NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗
NASA NTRS2002

Vestibular convergence patterns in vestibular nuclei neurons of alert primates

Sensory signal convergence is a fundamental and important aspect of brain function. Such convergence may often involve complex multidimensional interactions as those proposed for the processing of otolith and semicircular canal (SCC) information for the detection of translational head movements and the effective discrimination from physically congruent gravity signals. In the present study, we have examined the responses of primate rostral vestibular nuclei (VN) neurons that do not exhibit any eye movement-related activity using 0.5-Hz translational and three-dimensional (3D) rotational motion. Three distinct neural populations were identified. Approximately one-fourth of the cells exclusively encoded rotational movements (canal-only neurons) and were unresponsive to translation. The canal-only central neurons encoded head rotation in SCC coordinates, exhibited little orthogonal canal convergence, and were characterized with significantly higher sensitivities to rotation as compared to primary SCC afferents. Another fourth of the neurons modulated their firing rates during translation (otolith-only cells). During rotations, these neurons only responded when the axis of rotation was earth-horizontal and the head was changing orientation relative to gravity. The remaining one-half of VN neurons were sensitive to both rotations and translations (otolith + canal neurons). Unlike primary otolith afferents, however, central neurons often exhibited significant spatiotemporal (noncosine) tuning properties and a wide variety of response dynamics to translation. To characterize the pattern of SCC inputs to otolith + canal neurons, their rotational maximum sensitivity vectors were computed using exclusively responses during earth-vertical axis rotations (EVA). Maximum sensitivity vectors were distributed throughout the 3D space, suggesting strong convergence from multiple SCCs. These neurons were also tested with earth-horizontal axis rotations (EHA), which would activate both vertical canals and otolith organs. However, the recorded responses could not be predicted from a linear combination of EVA rotational and translational responses. In contrast, one-third of the neurons responded similarly during EVA and EHA rotations, although a significant response modulation was present during translation. Thus this subpopulation of otolith + canal cells, which included neurons with either high- or low-pass dynamics to translation, appear to selectively ignore the component of otolith-selective activation that is due to changes in the orientation of the head relative to gravity. Thus contrary to primary otolith afferents and otolith-only central neurons that respond equivalently to tilts relative to gravity and translational movements, approximately one-third of the otolith + canal cells seem to encode a true estimate of the translational component of the imposed passive head and body movement.

Non-NASA Center↗