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 73 records · Page 4

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↗

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗

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↗

Performance assessment of near-fault buildings subjected to physics-based simulated earthquake ground motions with fling step

The effects of the co-seismic static offset (known as fling step) and associated velocity pulses on civil structures have been difficult to study because the static offset is typically removed during the processing of earthquake ground motion records. Simulated ground motions contain fling features and require no processing; therefore, they create new opportunities for representing fling features in seismic hazard analysis and assessing their influence on the seismic demands on near-fault structures. We use physics-based fault rupture simulations to study the characteristics of ground motions with fling step and the sensitivity of the near-fault structural demands to strong fling features. We uncover that simulated ground motions with a large fling step tend to have higher spectral intensity than those without a fling step at the same rupture distance, especially at periods longer than 2 s. As a result, the structural demands on flexible buildings tend to be the most sensitive to the fling features. Statistical analysis suggests that the ground motion spectral shape (represented by spectral accelerations at multiple periods) is—in most cases—a sufficient predictor of the structural demands on near-fault low-rise and mid-rise buildings at locations that are susceptible to strong fling effects. Finally, ground motion record selection experiments reveal that representing the spectral shape features at periods that are most relevant to a given structure may be an effective strategy to reduce the bias in the estimated demands on near-fault long-period structures when the available database of records is considered deficient in fling features.

Fling step↗

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.↗

Shapiro steps and stability of skyrmions interacting with alternating anisotropy under the influence of ac and dc drives

Here we use atomistic simulations to examine the sliding dynamics of a skyrmion in a two-dimensional system containing a periodic one-dimensional stripe pattern of variations between low and high values of the perpendicular magnetic anisotropy. The skyrmion changes in size as it crosses the interface between two anisotropy regions. On applying combined dc and ac driving in either parallel or perpendicular directions, we observe a wide variety of Shapiro steps, Shapiro spikes, and phase-locking phenomena. The phase-locked orbits have two-dimensional dynamics due to the gyrotropic or Magnus dynamics of the skyrmions and are distinct from the phase-locked orbits found for strictly overdamped systems. Along a given Shapiro step when the ac drive is perpendicular to the dc drive, the velocity parallel to the ac drive is locked while the velocity in the perpendicular direction increases with increasing drive to form Shapiro spikes. At the transition between adjacent Shapiro steps, the parallel velocity jumps up to the next step value, and the perpendicular velocity drops. The skyrmion Hall angle shows a series of spikes as a function of increasing dc drive, where the jumps correspond to the transition between different phase-locked steps. At high drives, the Shapiro steps and Shapiro spikes are lost. When both the ac and dc drives are parallel to the stripe periodicity direction, Shapiro steps appear, while if the dc drive is parallel to the stripe periodicity direction and the ac drive is perpendicular to the stripe periodicity, then there are only two locked phases, and the skyrmion motion consists of a combination of sliding along the interfaces between the two anisotropy values and jumping across the interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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.↗

On the Error Covariance Correction Step of an ESKF Attitude Update

The attitude states of an error-state Kalman filter (ESKF) behave differently than most other states in the system due to their multiplicative (rather than additive) nature. One way in which they differ is an error covariance correction step after an ESKF error reset, which is not required for, for example, position and velocity states. This covariance correction step is not intuitive, and it has only been recently derived for coordinate transform matrices. The author of this memo, however, found the provided derivation in [1] confusing due to a lack of clarity surrounding the invoked reference frames, and clarity is required as there are at least 4 different ways to parameterize small-angle attitude errors in an ESKF. Furthermore, while reproducing the work, the author of this memo found a more straightforward derivation that provides additional insight into the correction step. This memo offers a derivation of the attitude error covariance correction step of an ESKF, which pays specific attention to the coordinate reference frames.

97 MATHEMATICS AND COMPUTING↗