Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinematics and dynamics”

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 685 records · Page 38

Kinematic Dependence of the Dilution Factor in the SpinQuest Experiment E1039 at Fermilab

Thirty years after the "Spin Crisis" began, there is still uncertainty surrounding the spin structure of the nucleon. SpinQuest Experiment E1039 at FermiLab will attempt to probe the internal transverse dynamics of the nucleon via Drell-Yan scattering of a 120 GeV proton beam from a polarized proton target to determine if the orbital angular momentum of the sea-quarks contributes to the overall spin of the nucleon. In order to understand a leading contribution to the overall systematic uncertainty in this experiment, a study was performed to calculate and evaluate the uncertainty in the dilution factor of the experimental target. This thesis catalogs that effort, highlights where it overlaps with other sources of systematic error, and suggests future work aimed at reducing the overall systematic error of the experiment.

Conover, Arthur↗

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

The Role of Musculoskeletal Dynamics and Neuromuscular Control in Stress Development in Bone

The role of forces produced by the musculotendon units in the stress development of the long bones during gait has not been fully analyzed. It is well known that the musculotendons act as actuators producing the joint torques which drive the body. Although the joint torques required to perform certain motor tasks can be recovered through a kinematic analysis, it remains a difficult problem to determine the actual forces produced by each muscle that resulted in these torques. As a consequence, few studies have focused on the role of individual muscles in the development of stress in the bone. This study takes a control theoretic approach to the problem. A seven-link, eight degrees of freedom model of the body is controlled by various muscle groups on each leg to simulate gait. The simulations incorporate Hill-type models of muscles with activation and contraction dynamics controlled through neural inputs. This direct approach allows one to know the exact muscle forces exerted by each musculotendon throughout the gait cycle as well the joint torques and reaction forces at the ankle and knee. Stress and strain computed by finite element analysis on skeletal members will be related to these derived loading conditions. Thus the role of musculoskeletal dynamics and neuromuscular control in the stress development of the tibia during gait can be analyzed.

DeWoody, Yssa↗

Do steady fast magnetic dynamos exist?

This paper considers the question of the existense of a steady fast kinematic magnetic dynamo for a conducting fluid with a steady velocity field and vanishingly small electrical resistivity. The analysis of examples of steady dynamos, found by considering the zero-resistivity dynamics, indicated that, for sufficiently small resistivity, dynamo action can indeed occur in steady smooth three-dimensional chaotic fluid flows and that fast dynamos should consequently be a typical occurrence for such flows.

Finn, John M.↗

The Migdal effect in semiconductors for dark matter with masses below ~ 100 MeV

Dark matter scattering off a nucleus has a small probability of inducing an observable ionization through the inelastic excitation of an electron, called the Migdal effect. We use an effective field theory to extend the computation of the Migdal effect in semiconductors to regions of small momentum transfer to the nucleus, where the final state of the nucleus is no longer well described by a plane wave. Our analytical result can be fully quantified by the measurable dynamic structure factor of the semiconductor, which accounts for the vibrational degrees of freedom (phonons) in a crystal. We show that, due to the sum rules obeyed by the structure factor, the inclusive Migdal rate and the shape of the electron recoil spectrum is well captured by approximating the nuclei in the crystal as free ions; however, the exclusive differential rate with respect to energy depositions to the crystal depends on the phonon dynamics encoded in the dynamic structure function of the specific material. Our results now allow the Migdal effect in semiconductors to be evaluated even for the lightest dark matter candidates (m χ ≳ 1 MeV) that can kinematically excite electrons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The oblique impingement of an axisymmetric jet

The mechanics of the oblique impingement of an axisymmetric jet on a plane surface are examined in detail. The stagnation point is discussed. A schematic drawing of the problem and coordinate system used to describe the flow field are given. The kinematic features of the flow above the plate are examined in the context of the conservation of mass, the vorticity of the jet, and the vorticity introduced by the jetplate interaction. The dynamic features of the flow are examined in terms of the surface pressure distribution and the cause-effect relationships which exist between the pressure and velocity/vorticity distributions. Flow calculations performed are given. The investigation is relevant to the flow resulting from the interaction of the propulsion jet with the main airfoil (STOL aircraft), and is appropriate to an over- or under- wing configuration.

Foss, J. F.↗

The Structure of Active Galactic Nuclei

We are continuing our systematic investigation of the nuclear structure of nearby active galactic nuclei (AGN). Upon completion, our study will characterize hypothetical constructs such as narrow-line clouds, obscuring tori, nuclear gas disks. and central black holes with physical measurements for a complete sample of nearby AGN. The major scientific goals of our program are: (1) the morphology of the NLR; (2) the physical conditions and dynamics of individual clouds in the NLR; (3) the structure and physical conditions of the warm reflecting gas; (4) the structure of the obscuring torus; (5) the population and morphology of nuclear disks/tori in AGN; (6) the physical conditions in nuclear disks; and (7) the masses of central black holes in AGN. We will use the Hubble Space Telescope (HST) to obtain high-resolution images and spatially resolved spectra. Far-UV spectroscopy of emission and absorption in the nuclear regions using HST/FOS and the Hopkins Ultraviolet Telescope (HUT) will help establish physical conditions in the absorbing and emitting gas. By correlating the dynamics and physical conditions of the gas with the morphology revealed through our imaging program, we will be able to examine mechanisms for fueling the central engine and transporting angular momentum. The kinematics of the nuclear gas disks may enable us to measure the mass of the central black hole. Contemporaneous X-ray observations using ASCA will further constrain the ionization structure of any absorbing material. Monitoring of variability in the UV and X-ray absorption will be used to determine the location of the absorbing gas, possibly in the outflowing warm reflecting gas, or the broad-line region, or the atmosphere of the obscuring torus. Supporting ground-based observations in the optical, near-IR, imaging polarimetry, and the radio will complete our picture of the nuclear structures. With a comprehensive survey of these characteristics in a complete sample of nearby AGN, our conclusions should be more reliably extended to AGN as a class.

Kriss, Gerard A.↗

Turbofan forced mixer-nozzle internal flowfield. Volume 2: Computational fluid dynamic predictions

A general program was conducted to develop and assess a computational method for predicting the flow properties in a turbofan forced mixed duct. The detail assessment of the resulting computer code is presented. It was found that the code provided excellent predictions of the kinematics of the mixing process throughout the entire length of the mixer nozzle. The thermal mixing process between the hot core and cold fan flows was found to be well represented in the low speed portion of the flowfield.

Werle, M. J.↗

Aeroelastic characteristics of the AH-64 bearingless tail rotor

The results of a wind tunnel test program to determine the performance loads and dynamic characteristics of the Composite Flexbeam Tail Rotor (CFTR) for the AH-64 Advanced Attack Helicopter are reported. The CFTR uses an elastomeric shear attachment of the flexbeam to the hub to provide soft-inplane S-mode and stiff-inplane C-mode configuration. The properties of the elastomer were selected for proper frequency placement and scale damping of the inplane S-mode. Kinematic pitch-lag coupling was introduced to provide the first cyclic inplane C-mode damping at high collective pitch. The CFTR was tested in a wind tunnel over the full slideslip envelop of the AH-64. It is found that the rotor was aeroelastically stable throughout the complete collective pitch range and up to rotor speeds of 1403 rpm. The dynamic characteristics of the rotor were found to be satisfactory at all pitch angles and rotor speeds of the tunnel tests. The design characteristics of the rotor which permit the high performance characteristics are discussed. Several schematic drawings and photographs of the rotor are provided.

Banerjee, D.↗

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗

Multivariable control of the Space Shuttle Remote Manipulator System using linearization by state feedback

This paper demonstrates an approach to nonlinear control system design that uses linearization by state feedback to allow faster maneuvering of payloads by the Shuttle Remote Manipulator System (SRMS). A nonlinear feedback law is defined to cancel the nonlinear plant dynamics so that a linear controller can be designed for the SRMS. First a nonlinear design model was generated via SIMULINK. This design model included nonlinear arm dynamics derived from the Lagrangian approach, linearized servo model, and linearized gearbox model. The current SRMS position hold controller was implemented on this system. Next, a trajectory was defined using a rigid body kinematics SRMS tool, KRMS. The maneuver was simulated. Finally, higher bandwidth controllers were developed. Results of the new controllers were compared with the existing SRMS automatic control modes for the Space Station Freedom Mission Build 4 Payload extended on the SRMS.

Gettman, Chang-Ching L.↗

Solving challenges in QCD 3D partonic distributions

The superior luminosity and high polarization of CEBAF, combined with the high resolution and excellent particle identification capabilities of its detectors, as well as the ability for multidimensional and multiparticle detection using polarized targets, make Jefferson Lab uniquely positioned to disentangle the genuine intrinsic transverse structure of hadrons encoded in 3D partonic distributions—particularly in the kinematic regime dominated by valence quarks. The inclusion of Jefferson Lab data on semi-inclusive and hard exclusive hadron production has the potential to dramatically advance our understanding of non-perturbative QCD dynamics. Although a wealth of data from various Jefferson Lab experiments is already available, its incorporation into phenomenological studies has been slow, and a significant portion of data from other leptoproduction experiments is still missing from global fits. In this contribution, we discuss the existing challenges and outline a path forward for improving the analysis of low center-of-mass electroproduction experiments in general, and Jefferson Lab data in particular.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stability and slosh wave for rotating bubble due to gravity-jitters under microgravity

Time-dependent fluid behaviors in microgravity environments are examined with particular attention given to the dynamics of bubble configurations. The initial bubble profiles, computed from the steady-state formulation, are used as the initial input for the time dependent computation. The parameters considered are liquid density and its kinematic viscosity, gas density and its kinematic viscosity, surface tension coefficient, angular velocity, and gravity environment. This computer algorithm can be used to simulate the fluid behavior in a microgravity environment, in particular the excitation of slosh waves due to different frequencies of gravity jitters. Results show that lower-frequency gravity jitters excite slosh waves with a higher ratio of maximum amplitude to wave length than that of the slosh waves generated by the higher frequency gravity jitters.

Hung, R. J.↗

Mobility analysis, simulation, and scale model testing for the design of wheeled planetary rovers

The use of computer based techniques to model and simulate wheeled rovers on rough natural terrains is considered. Physical models of a prototype vehicle can be used to test the correlation of the simulations in scaled testing. The computer approaches include a quasi-static planar or two dimensional analysis and design tool based on the traction necessary for the vehicle to have imminent mobility. The computer program modeled a six by six wheel drive vehicle of original kinematic configuration, called the Rocker Bogie. The Rocker Bogie was optimized using the quasi-static software with respect to its articulation parameters prior to fabrication of a prototype. In another approach used, the dynamics of the Rocker Bogie vehicle in 3-D space was modeled on an engineering workstation using commercial software. The model included the complex and nonlinear interaction of the tire and terrain. The results of the investigation yielded numerical and graphical results of the rover traversing rough terrain on the earth, moon, and Mars. In addition, animations of the rover excursions were also generated. A prototype vehicle was then used in a series of testbed and field experiments. Correspondence was then established between the computer models and the physical model. The results indicated the utility of the quasi-static tool for configurational design, as well as the predictive ability of the 3-D simulation to model the dynamic behavior of the vehicle over short traverses.

Lindemann, Randel A.↗

STEREO/SECCHI Stereoscopic Observations Constraining the Initiation of Polar Coronal Jets

We report on the first stereoscopic observations of polar coronal jets made by the EUVI/SECCHI imagers on board the twin STEREO spacecraft. The significantly separated viewpoints (approximately 11 degrees ) allowed us to infer the 3D dynamics and morphology of a well-defined EUV coronal jet for the first time. Triangulations of the jet's location in simultaneous image pairs led to the true 3D position and thereby its kinematics. Initially the jet ascends slowly at approximately equal to 10-20 kilometers per second and then, after an apparent 'jump' takes place, it accelerates impulsively to velocities exceeding 300 kilometers per second with accelerations exceeding the solar gravity. Helical structure is the most important geometrical feature of the jet which shows evidence of untwisting. The jet structure appears strikingly different from each of the two STEREO viewpoints: face-on in the one viewpoint and edge-on in the other. This provides conclusive evidence that the observed helical structure is real and is not resulting from possible projection effects of single viewpoint observations. The clear demonstration of twisted structure in polar jets compares favorably with synthetic images from a recent MHD simulation of jets invoking magnetic untwisting as their driving mechanism. Therefore, the latter can be considered as a viable mechanism for the initiation of polar jets.

Patsourakos, S.↗

Comparing Theoretically Scaled Biomechanical Models

BACKGROUND An investigation into incorporating space suit aspects required a 50th-percentile male model which was not part of the existing dataset of models. Previously, theoretical models for 5th-percentile female and 95th-percentile male were created utilizing a scaled test subject as close as possible to the target height and weight. Scaling factors were created by taking the height ratio and applying it uniformly to the model body segments, then fine tuning to ensure the theoretical model’s height is as expected. This study was initiated to examine the existing method of isometrically scaling in OpenSim [1,2] and to create alternative methods which do not rely on an existing subject being close in height and weight to the theoretical model of interest. Generating a theoretical biomechanical model provides additional abilities without the reliance on available OpenSim models or real subjects. Gained abilities include creating different percentile models and attaining representative anthropometry for any target height/weight, such as targeting specific crew populations or gaps within the current dataset. METHODS AND RESULTS This study includes 3 methods for generating the theoretical model; utilizing the modified unscaled OpenSim Full Body Rajagopal Model (FBRM) [3,4], a dataset of scaled OpenSim models, and existing scaling factors from Dumas et al. [5]. All methods use the Anthropometric Survey of US Army Personnel (ANSUR II) [6] as an input of necessary anthropometric measurements. The following measurements are retrieved from the ANSUR II collection: mass, stature, cervicale height, acromial height, axilla height, waist height, trochanterion height, lateral femoral epicondyle height, lateral malleolus height, acromion-radiale length, radiale-stylion length, palm length, ball of foot length, bicristal breadth, bimalleolar breadth. Some of the measurements are direct segment lengths and others are utilized to derive segment lengths. Height, weight, and age ranges are the required inputs to parse the collection for a mean value of the measurements. A simple iterative process may be necessary for the output of mean mass and stature consistent with the theoretical model of interest. In some cases, one may specify exact values instead of a range but that is dependent upon whether those exact values pertain to a single subject from the ANSUR II collection. Equations from Dumas et al. were used to calculate scaling factors for the two methods that utilize OpenSim scaled and unscaled models. In the third method, Dumas’ scaling factors were used directly instead of generating our own. After scaling factors are achieved, they are applied along with the mass and segment lengths retrieved from ANSUR II in order to calculate the segment mass, Center of Mass (CoM), mass Moment of Inertia (MoI), and the joint location in the parent frame. The theoretical models from the different methods were compared to the previous existing method of isometrically scaling in OpenSim, as well as the standards found in NASA STD-3000 [7] and the NASA Human Integration Design Handbook [8]. The whole-body center of mass is the main comparison performed between the models and the standards. The body segment mass properties were also compared when possible. In some cases, the standards do not provide complete data, or the center of mass location is not provided with respect to the appropriate joint coordinate system. For a 50th-percentile male subject, results from isometric scaling in OpenSim and the method utilizing the unscaled FBRM were compared to the standards in NASA STD-3000. According to the standard, the whole-body vertical CoM location was taken from the head vertex to the CoM. The standards provide 80.2 cm for the CoM vertical component whereas isometric scaling in OpenSim and utilizing the unscaled FBRM provide 82.2 and 82.3 cm, respectively. Through definition, the whole-body CoM in the medial/lateral direction agreed across the models. Investigation is on-going to clearly identify the reference point pertaining to dorsal/ventral whole-body CoM position in standards and relate that to the theoretical models. REFERENCES [1] Delp, S.L., et al., “OpenSim: Open-source Software to Create and Analyze Dynamic Simulations of Movement”, IEEE Transactions on Biomedical Engineering, (2007). [2] Seth, A., Hicks J.L., Uchida, T.K., Habib, A., Dembia, C.L., Dunne, J.J., Ong, C.F., DeMers, M.S., Rajagopal, A., Millard, M., Hamner, S.R., Arnold, E.M., Yong, J.R., Lakshmikanth, S.K., Sherman,n M.A., Delp, S.L. OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. Plos Computational Biology, 14(7). (2018) [3] R. K. Huffman, W. K. Thompson, C. A. Gallo, L. J. Quiocho, “Improvement of Scaling and Inverse Kinematic Results with Additional Upper Body Joints Added to Opensim Rajagopal Model”, NASA Human Research Program Investigator’s Workshop, (2019). [4] Huffman R.K., Thompson W., Gallo C., “Modified OpenSim Rajagopal Full Body Model”, New Technology Report, MSC-26872-1. (August 2020). [5] Dumas R., Cheze L. and Verriest J., 2007. “Adjustments to McConville et al. and Young et al. body segment inertial parameters”. Journal of Biomechanics 40, pp. 543-553 [6] Gordon C.C., 2012 “Anthropometric Survey of U.S. Army Personnel: Methods and Summary Statistics”, US Army Natick Soldier RD&E Center. [7] NASA STD-3000/REV-B, 1995. “The Man-System Integration Standards” [8] NASA/SP-2010-3407/REV-1, 2014. “The Human Integration Design Handbook (HIDH)”

F N Matari↗

First inverse kinematics measurement of resonances in Be 7 ( α , γ ) C 11 relevant to neutrino-driven wind nucleosynthesis using DRAGON

A possible mechanism to explain the origin of the light $\textit{p}$ nuclei in the Galaxy is the nucleosynthesis in the proton-rich neutrino-driven wind ejecta of core-collapse supernovas via the $\textit{νp}$ process. However, this production scenario is very sensitive to the underlying supernova dynamics and the nuclear physics input. As far as the nuclear uncertainties are concerned, the breakout from the $\textit{pp}$ chains via the 7 Be(α,γ) 11 C reaction has been identified as an important link which can influence the nuclear flow and, therefore, the efficiency of the $\textit{νp}$ process. However, its reaction rate is poorly known over the relevant temperature range, $\textit{T}$ = 1.5–3 GK. We report on the first direct measurement of two resonances of the 7 Be(α,γ) 11 C reaction with previously unknown strengths using an intense radioactive 7 Be beam from the Isotope Separator and Accelerator (ISAC-I) Center facility and the DRAGON recoil separator in inverse kinematics. In this work, we have decreased the 7 Be(α,γ) 11 C reaction rate uncertainty to ≈9.4–10.7% over the relevant temperature region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗