Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Characterization of the Head Stabilization Response to a Lateral Perturbation During Walking in Older Adults

A main contributor of fractures in older adults is from a lateral fall. The decline in sensory systems results in difficulty maintaining balance stability. Head stabilization contributes to postural control by serving as a stable platform for the sensory systems. The purpose of this study was to characterize the head stabilization response to a lateral perturbation while walking. A total of 16 healthy older adults, aged 66-81 years, walked across a foam pathway 6 times. One piece of the foam pathway covered a movable platform that translated to the left when the subject stepped on the foam. Three trials were randomized in which the platform shifted. Angular rate sensors placed on the center of mass of the head and trunk collected head and trunk movement in all three planes of motion. The roll plane was analyzed to examine motion in the plane of the perturbation. Subjects stepped onto the platform with the right foot. Recovery step time and distance were recorded. The first trial was analyzed to capture the novelty of the perturbation. Results indicate a significant difference in footfall distance t=0.004, p<0.05, as well as the speed of foot recovery t=0.001, p<0.05, between natural and perturbed walking. Results indicate that the head t=0.005, p<0.05, and trunk t=0.0001, p<0.05, velocities increase during perturbed compared to natural walking. Older adults place their recovery foot down faster when perturbed to re-establish their base of support. Head and trunk segments are less stable and move with greater velocities to reestablish stability when perturbed.

Buccello-Stout, Regina R.↗
NASA NTRS1978

Lean Combustion Limits of a Confined Premixed-Prevaporized Propane Jet

Lean blowout limits were reported for a premixed prevaporized propane jet issuing into a cylindrical combustor. A single hole in a flat plate was used as a flameholder. Flameholders with various hole diameters were used. Jet velocities were varied from 3 to 290 meters per second. The combustor cross sectional area was changed by using different quartz liners of 12.7 and 22.2 millimeters diameters. As a result the combustor Reynolds number varied from 1000 to 9000. Stability was achieved at laminar as well as turbulent conditions. Three zones of flame stability were observed. The blowout equivalence ratio varied with step size and the combustor and jet Reynolds numbers. The combustor inlet mixture temperature was 395 K, and the combustor pressure was 1 atmosphere.

Huck, K. L.↗
NASA NTRS1989

Microwave measurements of ocean currents for the SAXON experiment

Measurements of ocean surface currents carried out during the SAXON experiment using the University of Massachusetts stepped-frequency Delta-K radar (SFDKR) are presented. The measurement site was the Chesapeake light tower situated in the north Atlantic Ocean, 15 km off the coast of Virginia Beach, Virginia. The radar was placed on a platform 26 m above mean sea level. Measurements were made at grazing angles between 3 and 10 deg. A nearly uninterrupted record of ocean current velocity was obtained for a 26-day period. The measurement dependence on various system and hydrodynamic parameters and a comparison with an independent in-situ current measurement (EM current meter) are presented. Experimental results show that the Delta K spectrum with a distinguishable resonant line can be measured under all weather conditions, from calm seas with large swells to the fully developed seas occurring under winds as high as 20 ms-1. The only parameter that needs to be adjusted in real time is the frequency separation, Delta f, which has to be modified according to the varying shape of the ocean wave spectrum. Currents could not be accurately measured for very low wind conditions, even though resonant lines were present in the Delta K spectrum. The SFDKR can accurately measure wind-driven surface currents for wide range of wind conditions.

Popstefanija, Ivan↗
NASA NTRS2000

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗
NASA NTRSDate not supplied

Split-Step Simulations of Sonic Boom Propagation Beyond the Lateral Cutoff in a Turbulent Atmosphere

Recent flight tests during the Quiet Supersonic Flights 2018 (QSF18) study reported sonic booms heard outside of the primary carpet region. In the absence of turbulence, the lateral cutoff region separates the primary sonic boom carpet from the shadow zone, where the sonic boom signal experiences significant attenuation. However, when turbulence is present in the atmospheric boundary layer (ABL), additional scattering of the sonic boom to the shadow zone region occurs. A method is presented for simulating sonic boom propagation in a turbulent atmospheric boundary layer beyond the lateral cutoff region into the shadow zone. A split-step method is used to integrate a partially one-way equation for the acoustic pressure. Inhomogeneous turbulence, representative of the ABL, is generated in the computational domain with a Fourier synthesis approach. Distributions of several loudness metrics in the shadow zone region for a sonic boom N-wave and a shaped boom are examined. Increasing both turbulence root-mean-square velocity and integral length scale are found to increase the average loudness of booms in the shadow zone. (This research is supported by the Commercial Supersonic Technology Project of the National Aeronautics and Space Administration under Grant No. 80NSSC19K1685.)

sonic boom↗
DOE OSTI2022

Detailed biomass fast pyrolysis kinetics integrated to computational fluid dynamic (CFD) and discrete element modeling framework: Predicting product yields at the bench-scale

Fast pyrolysis is an intricate process due to the variability and anisotropy of lignocellulosic biomass and the complicated chemistry and physics during conversion in a bubbling fluidized bed reactor (BFBR). The complexity of biomass fast pyrolysis lends itself well to computational fluid dynamics (CFD) and discrete element (DEM) analysis, which promises to reduce experimental time and its associated cost. This work investigated switchgrass fast pyrolysis simulated by computational fluid dynamics coupled with a discrete element method to track individual reacting biomass particles throughout a bench-scale BFBR reactor. We accounted for the fast pyrolysis chemistry through a comprehensive reaction scheme with secondary cracking reactions. We performed a three-step reduction for secondary cracking reactions to convert the full cracking scheme into a reduced scheme easily incorporated into our model. We assessed the impact of operational conditions on the steady-state yields of liquid bio-oil, non-condensable gases (NCG), at 550 °C over a range of fluidization numbers (2 – 6 Umf), reported as a ratio to the minimum fluidization velocity (Umf). At steady-state, the volatile bio-oil yield had a range of 49.3–50.4 wt%. Levoglucosan was the primary volatile component present with 21 wt% of the bio-oil while water was the second largest with 20 wt%. The reduction of the secondary reaction schemes did not appreciably affect the overall yields of switchgrass pyrolysis compared to the full secondary scheme.

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

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗
NASA NTRS1993

Jet to freestream velocity ratio computations for a jet in a crossflow

The flowfield induced by a single, subsonic jet exhausting perpendicularly from a flat plate into a subsonic crossflow has been numerically investigated. Time-averaged solutions were obtained using the thin-layer Navier-Stokes equations and two overlapping grids. Test cases were chosen to match available experimental data where the jet Mach number was 0.78 and the freestream Mach number was varied to represent effective velocity ratios R from 4 to 12. Comparisons of the pressures induced on the flat plate are presented. The results show that the best agreement is obtained for R of 4 and 6, R of 8 is more difficult to resolve, and the R of 12 solution was not satisfactory. It is anticipated that the resolution could be improved by using either smaller time steps or finer grid spacing as R increases.

Margason, Richard J.↗
NASA NTRS1996

Large-eddy simulation of a backward facing step flow using a least-squares spectral element method

We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Chan, Daniel C.↗
NASA NTRS1982

Steps toward the Hubble constant. VIII - The global value

The problem of determining the global Hubble constant is approached by using two distance indicators whose dispersion in absolute magnitude is small, and hence which are virtually free of the Malmquist bias. New data are presented on the mean absolute-magnitude of the brightest individual red and blue supergiant stars in very nearby galaxies, calibrated relative to the Cepheids. The results are used to calibrate the maximum absolute blue magnitude for two well-observed Type I supernovae (SNe I). The velocity-relative distance diagram for 16 SNe I with known maximum apparent blue magnitudes is then calibrated using distances for IC 4182 and NGC 4214, with the result that the global Hubble constant equals 50 + or - 7 km per second per megaparsec. The result is compared with the known age of globular clusters in the Galaxy to show that the cosmological and the evolutionary time scales now closely agree.

Sandage, A.↗
DOE OSTI2025

Vortical interactions in turbulent thermoacoustic systems

This study examines the dynamics of vortical interactions and their implications for mitigating thermoacoustic instability in a turbulent combustor. The regions of intense vortical interactions are identified as vortical communities in the network space of weighted directed vortical networks constructed from two-dimensional experimental velocity data. One can expect vortical interactions in the combustor to be strongest near the moment of vortex shedding, as the shed vortices gradually weaken due to dissipation while convecting downstream. However, we show that, during the state of thermoacoustic instability, there is a non-trivial consistent phase lag of approximately 52° between the shedding of the coherent structures from the backward-facing step and the time instant when the vortical interactions attain their local maximum value. We explain this phase lag by investigating the correlation between acoustic pressure fluctuations, spatio-temporal dynamics of coherent structures and vortical interactions in the reaction field of the combustor. We also show the aperiodic variation of vortical interactions during the states of combustion noise and aperiodic epochs of intermittency. Furthermore, the spatio-temporal evolution of pairs of vortical communities with the maximum inter-community interactions provides insight into explaining the critical regions detected in the reaction field during the states of intermittency and thermoacoustic instability, also identified in previous studies. As a result, we further show that the most efficient suppression of thermoacoustic instability via air microjet injection is achieved when steady air jets are introduced to disrupt the maximum inter-community interactions present during the state of thermoacoustic instability.

Sahay, Ankit [Indian Institute of Technology Madra↗
DOE OSTI2023

A TOPAS model for lens-based proton radiography

Abstract Objective. Proton Radiography can be used in conjunction with proton therapy for patient positioning, real-time estimates of stopping power, and adaptive therapy in regions with motion. The modeling capability shown here can be used to evaluate lens-based radiography as an instantaneous proton-based radiographic technique. The utilization of user-friendly Monte Carlo program TOPAS enables collaborators and other users to easily conduct medical- and therapy- based simulations of the Los Alamos Neutron Science Center (LANSCE). The resulting transport model is an open-source Monte Carlo package for simulations of proton and heavy ion therapy treatments and concurrent particle imaging. Approach. The four-quadrupole, magnetic lens system of the 800-MeV proton beamline at LANSCE is modeled in TOPAS. Several imaging and contrast objects were modelled to assess transmission at energies from 230–930 MeV and different levels of particle collimation. At different proton energies, the strength of the magnetic field was scaled according to βγ, the inverse product of particle relativistic velocity and particle momentum. Main results. Materials with high atomic number, Z, (gold, gallium, bone-equivalent) generated more contrast than materials with low-Z (water, lung-equivalent, adipose-equivalent). A 5-mrad collimator was beneficial for tissue-to-contrast agent contrast, while a 10-mrad collimator was best to distinguish between different high-Z materials. Assessment with a step-wedge phantom showed water-equivalent path length did not scale directly according to predicted values but could be mapped more accurately with calibration. Poor image quality was observed at low energies (230 MeV), but improved as proton energy increased, with sub-mm resolution at 630 MeV. Significance. Proton radiography becomes viable for shallow bone structures at 330 MeV, and for deeper structures at 630 MeV. Visibility improves with use of high-Z contrast agents. This modality may be particularly viable at carbon therapy centers with accelerators capable of delivering high energy protons and could be performed with carbon therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Multigrid simulation of asymptotic curved-duct flows using a semi-implicit numerical technique

Asymptotic flows inside curved ducts of rectangular as well as polar cross section are analyzed using the Navier-Stokes equations in terms of the axial velocity and vorticity and the cross-flow stream function. Numerical solutions of the three second-order coupled elliptic partial differential equations governing this flow are obtained efficiently using the coupled alternating-direction implicit (ADI) method as well as the multigrid strongly-implicit (SI) scheme. For the flow configuration studied, the ADI method is found to be more sensitive to the time steps used than is the SI scheme. Use of the multigrid-coupled-strongly-implicit (MG-SI) scheme makes it possible to efficiently obtain fine-grid solutions for configurations having strong secondary flow. It is shown that, for this asymptotic curved-duct flow, the similarity parameter of significance is the Dean's number K rather than the Reynolds number Re. Results are obtained for curved ducts with square cross sections for K up to 900, which here corresponds to Re = 9,000 for this internal flow configuration.

Ghia, K. N.↗
NASA NTRS1994

Location and characteristics of the reconnection X-line deduced from low-altitude satellite and radar observations

We present an analysis of a cusp ion step observed between two poleward-moving events of enhanced ionospheric electron temperature. From the computed variation of the reconnection rate and the onset times of the associated ionospheric events, the distance between the satellite and the X-line can be estimated, but with a large uncertainty due to that in the determination of the low-energy cut-off of the ion velocity distribution function, f(E). Nevertheless, analysis of the time series f(t) shows the reconnection site to be on the dayside magnetopause, consistent with the pulsating cusp model, and the best estimate of the X-line location is 13 R(E) from the satellite. The ion precipitation is used to reconstruct the field-parallel part of the Cowley-D ion distribution function injected into the open low latitude boundary layer (LLBL) in the vicinity of the X-line. From this the Alfven speed, plasma density, magnetic field, parallel ion temperature, and flow velocity of the magnetosheath near the X-line can be derived.

Lockwood, M.↗
NASA NTRS2016

The Effect of Acoustic Forcing on Instabilities and Breakdown in Swept-Wing Flow over a Backward-Facing Step

Instability interaction and breakdown were experimentally investigated in the flow over a swept backward-facing step. Acoustic forcing was used to excite the Tollmien-Schlichting (TS) instability and to acquire phase-locked results. The phase-averaged results illustrate the complex nature of the interaction between the TS and stationary cross flow instabilities. The weak stationary cross flow disturbance causes a distortion of the TS wavefront. The breakdown process is characterized by large positive and negative spikes in velocity. The positive spikes occur near the same time and location as the positive part of the TS wave. Higher-order spectral analysis was used to further investigate the nonlinear interactions between the TS instability and the traveling cross flow disturbances. The results reveal that a likely cause for the generation of the spikes corresponds to nonlinear interactions between the TS, traveling cross flow, and stationary cross flow disturbances. The spikes begin at low amplitudes of the unsteady and steady disturbances (2-4% U (sub e) (i.e. boundary layer edge velocity)) but can achieve very large amplitudes (20-30 percent U (sub e) (i.e. boundary layer edge velocity)) that initiate an early, though highly intermittent, breakdown to turbulence.

Eppink, Jenna L.↗
NASA NTRS2009

Cooperative Three-Robot System for Traversing Steep Slopes

Teamed Robots for Exploration and Science in Steep Areas (TRESSA) is a system of three autonomous mobile robots that cooperate with each other to enable scientific exploration of steep terrain (slope angles up to 90 ). Originally intended for use in exploring steep slopes on Mars that are not accessible to lone wheeled robots (Mars Exploration Rovers), TRESSA and systems like TRESSA could also be used on Earth for performing rescues on steep slopes and for exploring steep slopes that are too remote or too dangerous to be explored by humans. TRESSA is modeled on safe human climbing of steep slopes, two key features of which are teamwork and safety tethers. Two of the autonomous robots, denoted Anchorbots, remain at the top of a slope; the third robot, denoted the Cliffbot, traverses the slope. The Cliffbot drives over the cliff edge supported by tethers, which are payed out from the Anchorbots (see figure). The Anchorbots autonomously control the tension in the tethers to counter the gravitational force on the Cliffbot. The tethers are payed out and reeled in as needed, keeping the body of the Cliffbot oriented approximately parallel to the local terrain surface and preventing wheel slip by controlling the speed of descent or ascent, thereby enabling the Cliffbot to drive freely up, down, or across the slope. Due to the interactive nature of the three-robot system, the robots must be very tightly coupled. To provide for this tight coupling, the TRESSA software architecture is built on a combination of (1) the multi-robot layered behavior-coordination architecture reported in "An Architecture for Controlling Multiple Robots" (NPO-30345), NASA Tech Briefs, Vol. 28, No. 10 (October 2004), page 65, and (2) the real-time control architecture reported in "Robot Electronics Architecture" (NPO-41784), NASA Tech Briefs, Vol. 32, No. 1 (January 2008), page 28. The combination architecture makes it possible to keep the three robots synchronized and coordinated, to use data from all three robots for decision- making at each step, and to control the physical connections among the robots. In addition, TRESSA (as in prior systems that have utilized this architecture) , incorporates a capability for deterministic response to unanticipated situations from yet another architecture reported in Control Architecture for Robotic Agent Command and Sensing (NPO-43635), NASA Tech Briefs, Vol. 32, No. 10 (October 2008), page 40. Tether tension control is a major consideration in the design and operation of TRESSA. Tension is measured by force sensors connected to each tether at the Cliffbot. The direction of the tension (both azimuth and elevation) is also measured. The tension controller combines a controller to counter gravitational force and an optional velocity controller that anticipates the motion of the Cliffbot. The gravity controller estimates the slope angle from the inclination of the tethers. This angle and the weight of the Cliffbot determine the total tension needed to counteract the weight of the Cliffbot. The total needed tension is broken into components for each Anchorbot. The difference between this needed tension and the tension measured at the Cliffbot constitutes an error signal that is provided to the gravity controller. The velocity controller computes the tether speed needed to produce the desired motion of the Cliffbot. Another major consideration in the design and operation of TRESSA is detection of faults. Each robot in the TRESSA system monitors its own performance and the performance of its teammates in order to detect any system faults and prevent unsafe conditions. At startup, communication links are tested and if any robot is not communicating, the system refuses to execute any motion commands. Prior to motion, the Anchorbots attempt to set tensions in the tethers at optimal levels for counteracting the weight of the Cliffbot; if either Anchorbot fails to reach its optimal tension level within a specified time, it sends message to the other robots and the commanded motion is not executed. If any mechanical error (e.g., stalling of a motor) is detected, the affected robot sends a message triggering stoppage of the current motion. Lastly, messages are passed among the robots at each time step (10 Hz) to share sensor information during operations. If messages from any robot cease for more than an allowable time interval, the other robots detect the communication loss and initiate stoppage.

Stroupe, Ashley↗
NASA NTRS1982

Utilization of normal mode initial conditions for detecting errors in the dynamics part of primitive equation global models

When a global atmospheric basic state has constant angular velocity and its temperature varies with altitude only, there exist normal mode solutions to the linearized global primitive equations. The use of these normal modes, which have known behavior in time, is superior to the use of the Rossby-Haurwitz wave as initial conditions for detecting errors in the dynamics part of primitive equation global models. With these initial conditions, integration through only one time step is sufficient to detect many formulation and coding errors. Other tests are still required for detecting problems of nonlinear instability and conservation of integral properties, however.

Chao, W. C.↗