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 55 records · Page 3

High-Frequency Secondary Instabilities Downstream of a Forward-Facing Step

Localized high-frequency velocity fluctuations were measured downstream of a forward-facing step on a 30◦swept flat plate with an imposed pressure gradient. Phase-locked time-resolved particle image velocimetry measurements were performed to determine the source of the fluctuations, which ultimately lead to transition. The steps result in localized regions of reversed flow for a short region downstream of the step. These regions are highly localized due to the impact of the stationary crossflow vortices. The velocity fluctuations typically increase significantly in amplitude shortly downstream of reattachment in these localized regions, indicating that the fluctuations are related to reattachment of the shear layer. The phase-averaged measurements confirm that the unsteadiness in the shear layer results in the shedding of a plethora of vortices, some oriented in the streamwise direction, and others in the spanwise direction. These vortices are impacted by the local properties of the mean flow, such as streamwise vorticity and streamwise velocity gradient. The origin of the unsteadiness in the shear layer is not definitively shown, but evidence suggests it is due to a shear-layer instability. Ultimately, the unsteadiness of the shear layer, and the vortices that are shed as a result of this unsteadiness, are responsible for transition in this scenario.

Jenna L Eppink↗

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗

Separated and Recovering Turbulent Boundary Layer Flow Behind a Backward Facing Step For Different Reynolds Numbers

Experimental results for a two-dimensional separated turbulent boundary layer behind a backward facing step for five different Reynolds numbers are reported. Results are presented in the form of tables, graphs and a floppy disk for an easy access of the data. Reynolds number based on the step height was varied by changing the reference velocity upstream of the step, U(sub o), and the step height, h. Hot-wire measurement techniques were used to measure three Reynolds stresses and four triple-velocity correlations. In addition, surface pressure and skin friction coefficients were measured. All hot-wire measurements were acquired in a measuring domain which excluded recirculating flow region due to the directional insensitivity of hot-wires. The downstream extent of the domain from the step was 51 h for the largest and I 14h for the smallest step height. This significant downstream length permitted extensive study of the flow recovery. Prediction of perturbed flows and their recovery is particularly attractive for popular turbulence models since variations of turbulence length and time scales and flow interactions in different regions are generally inadequately predicted. The data indicate that the flow in the free shear layer region behaves like the plane mixing layer up to about 2/3 of the mean reattachment length when the flow interaction with the wall commences the flow recovery to that of an ordinary turbulent boundary layer structure. These changes of the flow do not occur abruptly with the change of boundary conditions. A reattachment region represents a transitional region where the flow undergoes the most dramatic adjustments to the new boundary conditions. Large eddies, created in the upstream free-shear layer region, are being torn, recirculated, reentrained back into the main stream interacting with the incoming flow structure. It is foreseeable that it is quite difficult to describe the physics of this region in a rational and quantitative manner other than statistical. Downstream of the reattachment point the flow recovers at different rates near the wall, in the newly developing internal boundary layer, and in the outer part of the flow. It appears that Reynolds stresses do not fully recover up to the longest recovery length of 114 h.

Jovic, Srba↗

Morphological Stability of Faceted Interfaces

The major focus of this investigation is to study the fundamentals of layer spreading mechanisms during growth of doped Ge (a facet forming material), and to determine the conditions for morphological instability of vicinal solid-liquid interfaces. The investigation will also lead to the determination of the effect of dopants on the layer growth kinetics, step free energy, and dopant capture by the advancing ledges. The theoretical treatment of growth of faceted interfaces indicates that the kinetics of a step on a growing vicinal interface considerably depends on its angle of inclination, the melt concentration, and characteristics of flow currents in the melt. The morphological stability of the interface also depends on these parameters, as well as on the density and spreading velocity of the steps. However, the treatment of the instability of the interface by the layer growth mechanism is rather difficult because it requires exact knowledge of the thermal and solutal fields, hydrodynamics of the melt, and supercooling at the interface. The results of recent space experiments of the principal investigator involving directional solidification of faceted Bi-Sn alloys have shown that the morphological stability of various crystallographic orientations is significantly affected by the anistropy in interfacial properties of the faceted alloy in general, and the interface kinetics in particular. These findings have also raised many important and fundamental questions, particularly with respect to the behavior of interfacial steps, which need to be addressed via additional groundbased and microgravity experiments. For the present investigation we will use a novel crystal growth technique which provides axial heat flux close to the solid-liquid boundary. The Axial Heat Processing (AHP) technique allows for precise control and determination of the heat and mass transfer close to the crystallization front, and the establishment of a planar interface over the entire cross-section of the growing crystal.

Abbaschian, Reza↗

Influence of a Backward-Facing Step on Swept-Wing Boundary-Layer Transition

Experimental measurements were performed on a swept flat-plate model with an airfoil leading edge and imposed chordwise pressure gradient to determine the effects of a backward-facing step on transition in a low-speed stationary crossflow-dominated boundary layer. Detailed hot-wire measurements were performed for three step heights ranging from 36 to 49% of the boundary-layer thickness at the step and corresponding to subcritical, nearly critical, and critical cases. In general, the step had a small localized effect on the growth of the stationary crossflow vortex, whereas the unsteady disturbance amplitudes increased with increasing step height. Intermittent spikes in instantaneous velocity began to appear for the two larger step heights. A physical explanation was provided for the mechanism leading to transition and the sudden movement in the transition front due to the critical steps. The large localized velocity spikes, which ultimately led to an intermittent breakdown of the boundary layer, were the result of nonlinear interactions of the different types of unsteady instabilities with each other and with the stationary crossflow vortices. Thus, the unsteady disturbances played the most important role in transition, but the stationary crossflow vortices also had a significant role via the modulation and the increased amplitude of the unsteady disturbances.

Eppink, Jenna L.↗

In-Time UAV Flight-Trajectory Estimation and Tracking Using Bayesian Filters

Rapid increase of UAV operation in the next decade in areas of on-demand delivery, medical transportation services, law enforcement, traffic surveillance and several others pose potential risks to the low altitude airspace above densely populated areas. Safety assessment of airspace demands the need for a novel UAV traffic management (UTM) framework for regulation and tracking of the vehicles. Particularly for low-altitude UAV operations, quality of GPS measurements feeding into the UAV is often compromised by loss of communication link caused by presence of trees or tall buildings in proximity to the UAV flight path. Inaccurate GPS locations may yield to unreliable monitoring and inaccurate prognosis of remaining battery life and other safety metrics which rely on future expected trajectory of the UAV. This work therefore proposes a generalized monitoring and prediction methodology for autonomous UAVs using in-time GPS measurements. Firstly, a typical 4D smooth trajectory generation technique from a series of waypoint locations with associated expected times-of-arrival based on B-spline curves is presented. Initial uncertainty in the vehicle's expected cruise velocity is quantified to compute confidence intervals along the entire flight trajectory using error interval propagation approach. Further, the generated planned trajectory is considered as the prior knowledge which is updated during its flight with incoming GPS measurements in order to estimate its current location and corresponding kinematic profiles. Estimation of position is denoted in dicrete state-space representation such that position at a future time step is derived from position and velocity at current time step and expected velocity at the future time step. A linear Bayesian filtering algorithm is employed to efficiently refine position estimation from noisy GPS measurements and update the confidence intervals. Further, a dynamic re-planning strategy is implemented to incorporate unexpected detour or delay scenarios. Finally, critical challenges related to uncertainty quantification in trajectory prognosis for autonomous vehicles are identified, and potential solutions are discussed at the end of the paper. The entire monitoring framework is demonstrated on real UAV flight experiments conducted at the NASA Langley Research Center.

Banerjee, Portia↗

Neural basis for eye velocity generation in the vestibular nuclei of alert monkeys during off-vertical axis rotation

Activity of "vestibular only" (VO) and "vestibular plus saccade" (VPS) units was recorded in the rostral part of the medial vestibular nucleus and caudal part of the superior vestibular nucleus of alert rhesus monkeys. By estimating the "null axes" of recorded units (n = 79), the optimal plane of activation was approximately the mean plane of reciprocal semicircular canals, i.e., lateral canals, left anterior-right posterior (LARP) canals or right anterior-left posterior (RALP) canals. All units were excited by rotation in a direction that excited a corresponding ipsilateral semicircular canal. Thus, they all displayed a "type I" response. With the animal upright, there were rapid changes in firing rates of both VO and VPS units in response to steps of angular velocity about a vertical axis. The units were bidirectionally activated during vestibular nystagmus (VN), horizontal optokinetic nystagmus (OKN), optokinetic after-nystagmus (OKAN) and off-vertical axis rotation (OVAR). The rising and falling time constants of the responses to rotation indicated that they were closely linked to velocity storage. There were differences between VPS and VO neurons in that activity of VO units followed the expected time course in response to a stimulus even during periods of drowsiness, when eye velocity was reduced. Firing rates of VPS units, on the other hand, were significantly reduced in the drowsy state. Lateral canal-related units had average firing rates that were linearly related to the bias or steady state level of horizontal eye velocity during OVAR over a range of +/- 60 deg/s. These units could be further divided into two classes according to whether they were modulated during OVAR. Non-modulated units (n = 5) were VO types and all modulated units (n = 5) were VPS types. There was no significant difference between the bias level sensitivities relative to eye velocity of the units with and without modulation (P > 0.05). The modulated units had no sustained change in firing rate in response to static head tilts and their phases relative to head position varied from unit to unit. The phase did not appear to be linked to the modulation of horizontal eye velocity during OVAR. The sensitivities of unit activity to eye velocity were similar during all stimulus modalities despite the different gains of eye velocity vs stimulus velocity during VN, OKN and OVAR. Therefore, VO and VPS units are likely to carry an eye velocity signal related to velocity storage.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

On spectral multigrid methods for the time-dependent Navier-Stokes equations

A splitting scheme is proposed for the numerical solution of the time-dependent, incompressible Navier-Stokes equations by spectral methods. A staggered grid is used for the pressure, improved intermediate boundary conditions are employed in the split step for the velocity, and spectral multigrid techniques are used for the solution of the implicit equations.

Zang, T. A.↗

On spectral multigrid methods for the time-dependent Navier-Stokes equations

A splitting scheme is proposed for the numerical solution of the time-dependent, incompressible Navier-Stokes equations by spectral methods. A staggered grid is used for the pressure, improved intermediate boundary conditions are employed in the split step for the velocity, and spectral multigrid techniques are used for the solution of the implicit equations.

Zang, T. A.↗

Numerical simulation of the flow field and fuel sprays in an IC engine

A two-dimensional model for axisymmetric piston-cylinder configurations is developed to study the flow field in two-stroke direct-injection Diesel engines under motored conditions. The model accounts for turbulence by a two-equation model for the turbulence kinetic energy and its rate of dissipation. A discrete droplet model is used to simulate the fuel spray, and the effects of the gas phase turbulence on the droplets is considered. It is shown that a fluctuating velocity can be added to the mean droplet velocity every time step if the step is small enough. Good agreement with experimental data is found for a range of ambient pressures in Diesel engine-type microenvironments. The effects of the intake swirl angle in the spray penetration, vaporization, and mixing in a uniflow-scavenged two-stroke Diesel engine are analyzed. It is found that the swirl increases the gas phase turbulence levels and the rates of vaporization.

Nguyen, H. L.↗

The Fractional Step Method Applied to Simulations of Natural Convective Flows

This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The development of the momentum and continuity equations for natural convection in a fluid, a permeable medium, and in a binary alloy undergoing directional solidification will be presented. Finally, results for natural convection in a pure liquid, natural convection in a medium with a constant permeability, and for directional solidification will be presented.

Westra, Douglas G.↗

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗

Simulation of three-dimensional separated flows with vortex filaments

A Lagrangian vorticity method for numerical simulation of three-dimensional separated flows with vortex filaments, about solid bodies, is discussed. The method is an extension of an earlier one for three-dimensional rotational flows away from solid boundaries e.g., jets, vortex rings, and aircraft trailing vortices after initial rollup. In this method an harmonic contribution to the velocity field is computed at each time step to ensure tangency of the velocity field at the boundary. The mechanics of the boundary layer are approximated with sufficient accuracy so that the separation lines are located correctly on the surface of the body.

Leonard, A.↗

Turbulent solution of the Navier-Stokes equations for an inhomogenous developing shear layer

To study the nonlinear physics of inhomogeneous turbulent shear flow, the unaveraged Navier-Stokes equations are solved numerically. For initial conditions a three-dimensional cosine velocity fluctuation and a mean-velocity profile with a step are used. Although the initial conditions are nonrandom. The flow soon becomes turbulent. Concentrated turbulent energy develops near the plane where the mean velocity gradient is initially infinite. The terms in the one-point correlation equation for turbulent energy, including those for the diffusion and production of turbulence, are calculated, the diffusion terms tend to make the turbulence more homogeneous.

Deissler, R. G.↗

Seismic velocity structure of the lunar mantle

The recently completed set of seismic arrival times from the Apollo lunar seismic network are inverted to estimate the average seismic velocities in three sections of the lunar mantle: two for the upper mantle and one for the middle mantle. The method used is a variation of the linearized least squares inversion where the inversion is accomplished in steps. The estimated average velocities in the upper mantle decrease from Vp = 7.74 km/sec and Vs = 4.49 km/sec in the section above 270-km depth to Vp = 7.46 km/sec and Vs = 4.25 km/sec in the section between 270- and 500-km depth, confirming the earlier finding of negative gradients based on seismic amplitude variations. The average velocities in the middle mantle between the depths of 500 km and 1000 km of Vp = 8.26 km/sec and Vs = 4.65 km/sec are significantly higher than those in the upper mantle, contradicting earlier estimates based on more limited data. The higher velocities may suggest initial melting of the moon down to at least 1000-km depth.

Nakamura, Y.↗

Direction-specific adaptation effects acquired in a slow rotation room

Thirty-eight subjects were required to execute 120 head movements in a slow rotation room at each 1-rpm increase in velocity of the room between 0 and 6 rpm and, after a single-step gradual return to zero velocity, execute 120 head movements either immediately after the return or after delay periods varying from 1 to 24 hours unless, at any time, more than mild symptoms of motion sickness were elicited. A second stress profile differed by the sequential addition of an incremental adaptation schedule in which the direction of rotation was reversed. The experimental findings demonstrated the acquisition of direction-specific adaptation effects that underwent spontaneous decay with a short time constant (hours). Speculations are presented which could account for the simultaneous acquisition of short-term and long-term adaptation effects. The findings support the theory that motion sickness, although a consequence of vestibular stimulation, has its immediate origin in nonvestibular systems, implying a faculative or temporary linkage between the vestibular and nonvestibular systems.

Graybiel, A.↗

Direction-specific adaptation effects acquired in a slow rotation room.

Thirty-eight subjects were required to execute 120 head movements in a slow rotation room at each 1-rpm increase in velocity of the room between 0 and 6 rpm and, after a single-step gradual return to zero velocity, execute 120 head movements either immediately after return or after delay periods varying from 1 to 24 hours unless, at any time, more than mild symptoms of motion sickness were elicited. A second stress profile differed by the sequential addition of an incremental adaptation schedule (identical to the first) in which the direction of rotation was reversed. The experimental findings demonstrated the acquisition of direction-specific adaptation effects that underwent spontaneous decay with a short time constant (hours). With their disappearance a nondirection-specific adaptation was revealed with a long-time constant (days). Speculations are presented which could account for the simultaneous acquisition of short-term and long-term adaptation effects.

Graybiel, A.↗

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗