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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Calibrating radar wind profiler reflectivity factor using surface disdrometer observations

Abstract. This study uses surface disdrometer reflectivity factor estimates to calibrate the vertical and off-vertical pointing radar beams produced by an ultra high frequency (UHF) band radar wind profiler (RWP) deployed at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) Central Facility in northern Oklahoma from April 2011 through July 2019. The methodology consists of five steps. First, the recorded Doppler velocity power spectra are adjusted to account for Nyquist velocity aliasing and coherent integration filtering effects. Second, the spectrum moments are calculated. The third step increases the signal-to-noise ratio (SNR) due to inflated noise power estimates during convective rain events that cause SNR to be biased low. The fourth step determines the RWP calibration constant for one radar beam (called the “reference” beam) by comparing uncalibrated RWP reflectivity factors at 500 m above the ground to 1 min resolution surface disdrometer reflectivity factors. The last step uses the calibrated reference beam reflectivity factor to calibrate the other radar beams during precipitation. There are two key findings. The RWP sensitivity decreased by approximately 3 to 4 dB yr−1 as the hardware aged. This drift was slow enough that the reference calibration constant can be estimated over 3-month intervals using episodic rain events. The calibrated moments are available on the DOE ARM data archive, and the Python processing code is available on public repositories.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2001

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

NASA NTRS2003

Step Bunch Evolution on Vicinal Faces of KDP

For in-situ studies of the formation and evolution of step patterns in solution growth, we have assembled an experimental setup based on Michelson interferometry with the growing crystal surface as one of the reflective surfaces. The device allows data collection over a relatively large area (approximately 4 sq. mm) in situ and in real time during growth. The depth resolution is improved over traditional interferometry using phase-shifted images combining by a suitable algorithm. We achieve a depth resolution of approximately 50 Angstroms. Lateral resolution, dependent on the degree of magnification, is around 0.3 to 5 microns. The crystal chosen as a model in this work is potassium dihydrogen phosphate (KDP), the optically non-linear material widely used in frequency doubling applications. Kinetics of KDP crystallization is well studied so that KDP can serve as a benchmark for our investigations. We present quantitative results on the onset, initial stages and development of instabilities in moving step trains on vicinal crystal surfaces at varying supersaturation, flow rate, and flow direction. The kinetics data suggest that at low supersaturations, step bunching is caused by impurity retardation of the steps, while at higher supersaturations, we link the non-linearity during growth to interdependence of the velocity and density of the steps evidenced in independent experiments. The behavior on the surface is very dynamic, small bunches both merge and split from larger bunches as they travel across the facet. We present evidence that despite these dynamics, under steady conditions there exists a limiting value to step bunch height. This height is reached at distances between 600 and 1000 microns from the step source. In our experiments, we observed the retention of this step bunch height limit up to the path of 1500 microns.

Booth, N. A.↗

NASA NTRS1992

LDA measurements in a Mach 2 flow over a rearward facing step with staged transverse injection

Measurements of the mean velocity field and selected turbulence statistics have been obtained via 2D LDA in a Mach 2 flow over a rearward facing step with downstream transverse injection. Axial mean velocity profiles, profiles of the axial and normal RMS levels, and the correlation coefficient between the axial and normal fluctuating components are presented for locations upstream of the step. These data, through comparison with other reported measurements in zero pressure gradient compressible boundary layers indicate that the boundary layers upstream of the step are consistent with 'quasi-equilibrium' turbulent boundary layers with the adiabatic wall boundary condition. Mean velocity field measurements on the symmetry plane of the tunnel are compared to laser induced iodine fluorescence measurements reported in the literature and obtained in the same facility. This comparison demonstrates the quality of the present LDA data set and shows that particle lag is not significant in the LDA measurements despite the complex nature of the downstream flow.

Wang, J.-A.↗

NASA NTRS1973

Error analysis of multi-conic techniques

A general error analysis of three recently developed multi-conic methods of three-body trajectory integration has been carried out. Single-step error functions for position and velocity have been derived as Taylor series in powers of the time step and also in integral form. These error functions are used to investigate the relative accuracy of the three methods in various regions of the earth-moon space and to provide a method of variable step size control for the trajectory integration procedure. Numerical results are used to compare the multi-step performance of the methods for both large and small step sizes.

D'Amario, L. A.↗

NASA NTRS1976

Turbulent atmospheric flow over a backward-facing step

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.

Kaul, U. K.↗

NASA NTRS1981

Approximate factorization for incompressible flow

For computational solution of the incompressible Navier-Stokes equations, the approximate factorization (AF) algorithm is used to solve the vectorized momentum equation in delta form based on the pressure calculated in the previous time step. The newly calculated velocities are substituted into the pressure equation (obtained from a linear combination of the continuity and momentum equation), which is then solved by means of line SOR. Computational results are presented for the NACA 66 sub 3 018 airfoil at Reynolds numbers of 1000 and 40,000 and attack angles of 0 and 6 degrees. Comparison with wind tunnel data for Re = 40,000 indicates good qualitative agreement between measured and calculated pressure distributions. Quantitative agreement is only fair, however, with the calculations somewhat displaced from the measurements. Furthermore, the computed velocity profiles are unrealistically thick around the airfoil, due to the excessive amount of artificial viscosity needed for stability. Based on the performance of the algorithm with regard to stability, it is concluded that AF/SOR is suitable for calculations at Reynolds numbers less than 10,000. Speedwise, the method is faster than point SOR by at least a factor of two.

Bernard, R. S.↗

NASA NTRS2003

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

NASA NTRS1949

The Compressible Flow Past Various Plane Profiles Near Sonic Velocity

In an earlier report UM No.1117 by Gothert,the single-source method was applied to the compressible flow around circles, ellipses, lunes, and around an elongated body of revolution at different Mach numbers and the results compared as far as possible with the calculations by Lamla ad Busemann. Essentially, it was found that with favorable source arrangement the single-source method is in good agreement with the calculations of the same degree of approximation by.Lamla and Busemann. Near sonic velocity the number of steps must be increased considerably in order to sufficiently approximate the adiabatic curve. After exceeding a certain Mach number where local supersonic fields occur already, it was no longer possible, in spite of the substantially increased number of steps, to obtain a systematic solution because the calculation diverged. This result,was interpreted to mean that above this point of divergence the symmetrical type of flow ceases to exist and changes into the unsymmetrical type characterized by compressibility shocks.

Goethert, B.↗

NASA NTRSDate not supplied

Deep Space navigation for the BioSentinel spacecraft science orbit

BioSentinel is an astrobiology small spacecraft mission. The payload consists of two parts, the first has optical and microfluidics sensors, and the second is a Linear Energy Transfer spectrometer that has the objective to measure deep space radiation from events such as coronal mass ejections. The goal of the mission is to observe potential DNA damage due to the radiation in heliocentric space on the living organism Saccharomyces cerevisiae, which is a budding yeast. Two types of this living organism are included in the payload. The first is a natural type that is more radiation tolerant, while the second is a mutant strain that has a deficiency in a gene that allows DNA repair once damage occurs. The impact caused by the radiation on the DNA is compared to an identical sample aboard the International Space Station, as well as another identical sample at a laboratory on the ground. The BioSentinel mission consists of a 6U CubeSat currently ,as of January 2024, active in heliocentric orbit. The spacecraft was launched aboard the first SLS flight as part of the Artemis-I campaign in November 2022. After successful deployment from the launch vehicle, it performed a lunar flyby with an altitude of 406 km. The delta-V imparted by the flyby provided the necessary energy to achieve a heliocentric orbit, in an Earth-trailing pattern. The navigation analysis consisted of a Kalman-filter that utilized data from the Deep Space Network and the ESA Estrack network. All those antennas were needed since the Artemis-1 campaign included the deployment of several other cubesats, therefore the scheduling process required more antenna assets than usual due to simultaneous demands from various missions. The processed tracking data was later also refined with a smoother in order to obtain a more accurate solution. The type of tracking data included TCP, Sequential Range, Doppler and Range formats. The solar radiation pressure coefficient, as well as the delta-V from the deployment and the flyby were modeled to obtain suitable solutions that could decrease the position and velocity uncertainties at several steps along the mission concept of operations. The final product each time resulted in updated ephemeris files that were used by the mission and the antenna networks as the mission progressed. Once in the final science orbit, the utilized antennas are only from the DSN network and the data format is bounded to just TCP. Regular orbit determination is performed, every two weeks. The spacecraft is in a nominal well-known orbit, performing regular operations. This paper includes an analysis of the final science orbit, the techniques and procedures utilized to perform orbit determination and a description of the overall navigation campaign produced during the mission and, more specifically, during the final science operations in Deep Space.

BioSentinel↗

NASA NTRS2011

Real-Time Thermographic-Phosphor-Based Temperature Measurements of Thermal Barrier Coating Surfaces Subjected to a High-Velocity Combustor Burner Environment

Surface temperature measurements were conducted on metallic specimens coated with an yttria-stabilized zirconia (YSZ) thermal barrier coating (TBC) with a YAG:Dy phosphor layer that were subjected to an aggressive high-velocity combustor burner environment. Luminescence-based surface temperature measurements of the same TBC system have previously been demonstrated for specimens subjected to static furnace or laser heating. Surface temperatures were determined from the decay time of the luminescence signal of the YAG:Dy phosphor layer that was excited by a pulsed laser source. However, the furnace and laser heating provides a much more benign environment than that which exists in a turbine engine, where there are additional challenges of a highly radiant background and high velocity gases. As the next step in validating the suitability of luminescence-based temperature measurements for turbine engine environments, new testing was performed where heating was provided by a high-velocity combustor burner rig at Williams International. Real-time surface temperature measurements during burner rig heating were obtained from the decay of the luminescence from the YAG:Dy surface layer. The robustness of several temperature probe designs in the sonic velocity, high radiance flame environment was evaluated. In addition, analysis was performed to show whether the luminescence decay could be satisfactorily extracted from the high radiance background.

Eldridge, Jeffrey I.↗

NASA NTRS1993

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗

NASA NTRS2015

The Effect of Backward-Facing Step Height on Instability Growth and Breakdown in Swept Wing Boundary-Layer Transition

A low-speed experiment was performed on a swept at plate model with an imposed pressure gradient to determine the effect of a backward-facing step on transition in a stationary-cross flow dominated flow. Detailed hot-wire boundary-layer measurements were performed for three backward-facing step heights of approximately 36, 45, and 49% of the boundary-layer thickness at the step. These step heights correspond to a subcritical, nearly-critical, and critical case. Three leading-edge roughness configurations were tested to determine the effect of stationary-cross flow amplitude on transition. The step caused a local increase in amplitude of the stationary cross flow for the two larger step height cases, but farther downstream the amplitude decreased and remained below the baseline amplitude. The smallest step caused a slight local decrease in amplitude of the primary stationary cross flow mode, but the amplitude collapsed back to the baseline case far downstream of the step. The effect of the step on the amplitude of the primary cross flow mode increased with step height, however, the stationary cross flow amplitudes remained low and thus, stationary cross flow was not solely responsible for transition. Unsteady disturbances were present downstream of the step for all three step heights, and the amplitudes increased with increasing step height. The only exception is that the lower frequency (traveling crossflow-like) disturbance was not present in the lowest step height case. Positive and negative spikes in instantaneous velocity began to occur for the two larger step height cases and then grew in number and amplitude downstream of reattachment, eventually leading to transition. The number and amplitude of spikes varied depending on the step height and cross flow amplitude. Despite the low amplitude of the disturbances in the intermediate step height case, breakdown began to occur intermittently and the flow underwent a long transition region.

Eppink, Jenna L.↗

DOE OSTI2023

Wave Detection and Tracking Within a Rotating Detonation Engine Through Object Detection

As the operational time window of experimental rotating detonation engines (RDEs) is expanded and the technology matures toward integration within gas turbines, monitoring techniques must evolve to offer computationally efficient and highly time-resolved diagnostics. In this study, computer vision object detection methodology that seeks to reduce data processing time and calculate wave velocity within drastically reduced time intervals as compared to traditional high-frame-rate RDE images analysis techniques is proposed. The adapted you-only-look-once object detection network is trained to detect individual detonation waves within single down-axis RDE images. The wave location and rotational direction detected within a frame are tracked through a series of high-speed images to calculate the frame-to-frame wave velocity with the time-step resolution of $\mathrm{20 μs}$ across a series of frames. The analysis of the annotation box size and image linearization effects is presented, demonstrating the lowest frame-to-frame velocity total uncertainty of $\mathrm{±3.8\%}$ and the highest classification speed of 9.5 frames per second using linearized images. Linearized images “unwrap” the RDE annulus pixel region to a reduced image size. Here, this new method offers great reductions in data processing times and unsteady detonation behavior insight at intervals more comparable to the timescales of detonation wave interactions via the application of machine learning to experimental RDE data.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1991

Faster Algorithm For Computation Of Incompressible Flow

Improved algorithm yields faster numerical solutions of Navier-Stokes equations of steady or unsteady three-dimensional flow of incompressible fluid. In artificial-compressibility method, unsteady flow treated as incompressible in advancing from one time step to next, but at each time step (or in steady state), fluid treated as having variable compressibility enabling propagation of flow field, and subiterations performed in increments of pseudotime until effects of compressibility subside. Directly couples pressure and velocity fields at same time step and converts elliptic incompressible Navier-Stokes equations to hyperbolic form more amenable to numerical integration.

Rogers, S. E.↗

DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗

NASA NTRS2000

Artificial gravity as a countermeasure in long-duration space flight

Long-duration exposure to weightlessness results in bone demineralization, muscle atrophy, cardiovascular deconditioning, altered sensory-motor control, and central nervous system reorganizations. Exercise countermeasures and body loading methods so far employed have failed to prevent these changes. A human mission to Mars might last 2 or 3 years and without effective countermeasures could result in dangerous levels of bone and muscle loss. Artificial gravity generated by rotation of an entire space vehicle or of an inner chamber could be used to prevent structural changes. Some of the physical characteristics of rotating environments are outlined along with their implications for human performance. Artificial gravity is the centripetal force generated in a rotating vehicle and is proportional to the product of the square of angular velocity and the radius of rotation. Thus, for a particular g-level, there is a tradeoff between velocity of rotation and radius. Increased radius is vastly more expensive to achieve than velocity, so it is important to know the highest rotation rates to which humans can adapt. Early studies suggested that 3 rpm might be the upper limit because movement control and orientation were disrupted at higher velocities and motion sickness and chronic fatigue were persistent problems. Recent studies, however, are showing that, if the terminal velocity is achieved over a series of gradual steps and many body movements are made at each dwell velocity, then full adaptation of head, arm, and leg movements is possible. Rotation rates as high as 7.5-10 rpm are likely feasible. An important feature of the new studies is that they provide compelling evidence that equilibrium point theories of movement control are inadequate. The central principles of equilibrium point theories lead to the equifinality prediction, which is violated by movements made in rotating reference frames. Copyright 2000 Wiley-Liss, Inc.

Review Literature↗

DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗