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 703 records · Page 39

Interchange Reconnection Dynamics in A Solar Coronal Pseudo-Streamer

Context. The generation of the slow solar wind remains an open problem in heliophysics. One of the current theories among those aimed at explaining the injection of coronal plasma in the interplanetary medium is based on interchange reconnection. It assumes that the exchange of magnetic connectivity between closed and open fields allows the injection of coronal plasma in the interplanetary medium to travel along the newly reconnected open field. However, the exact mechanism underlying this effect is still poorly understood. Aims. Our objective is to study this scenario in a particular magnetic structure of the solar corona: a pseudo-streamer. This topological structure lies at the interface between open and closed magnetic field and is thought to be involved in the generation of the slow solar wind. Methods. We performed innovative 3D magnetohydrodynamic (MHD) simulations of the solar corona with a pseudo-streamer, using the Adaptively Refined MHD Solver (ARMS). By perturbing the quasi-steady ambient state with a simple photospheric, large-scale velocity flow, we were able to generate a complex dynamics of the open-and-closed boundary of the pseudo-streamer. We studied the evolution of the connectivity of numerous field lines to understand its precise dynamics. Results. We witnessed different scenarios of opening of the magnetic field initially closed under the pseudo-streamer: one-step interchange reconnection dynamics, along with more complex scenarios, including a coupling between pseudo-streamer and helmet streamer, as well as back-and-forth reconnections between open and closed connectivity domains. Finally, our analysis revealed large-scale motions of a newly opened magnetic field high in the corona that may be explained by slipping reconnection. Conclusions. By introducing a new analysis method for the magnetic connectivity evolution based on distinct closed-field domains, this study provides an understanding of the precise dynamics underway during the opening of a closed field, which enables the injection of closed-field, coronal plasma in the interplanetary medium. Further studies shall provide synthetic observations for these diverse outgoing flows, which could be measured by Parker Solar Probe and Solar Orbiter.

magnetic fields↗

Sample interchange of MST radar data from the Urbana radar

As a first step in interchange of data from the Urbana mesosphere-stratosphere-troposphere (MST) radar, a sample tape has been prepared in 9-track 1600-bpi IBM format. It includes all Urbana data for April 1978 (the first month of operation of the radar). The 300-ft tape contains 260 h of typical mesospheric power and line-of-sight velocity data.

Bowhill, S. A.↗

Numerical investigation of chemically reacting flows in ramjet dump combustors

The time-dependent Navier-Stokes equations, including second-order turbulence model, are numerically integrated by using four-stage Runge-Kutta scheme to predict the steady-state supersonic flow structures in ramjet dump combustors. The formulation is derived for reacting flows with finite-rate chemistry. In the present study, it is firstly attempted to assess the accuracy of existing high-order turbulence model in supersonic flows. The comparison shows reasonable agreement between calculated and measured data in terms of velocity distributions. It is indicated that a modified constant C-mu for calculating turbulent eddy viscosity is needed in the supersonic flow regime and the adaptive meshing is preferred to capture the recirculation zone. In the reacting flow calculation, the results from a test case of hydrogen and air combustion at premixed conditon show that the rearward facing step is able to increase flow residence time and stabilize the flame in supersonic flows.

Hsieh, Kwang-Chung↗

Energy considerations in computational aeroacoustics

A finite-volume multistage time-stepping Euler code is used to investigate the use of CFD algorithms for the direct calculation of acoustics. The 2D compressible inviscid flow about an accelerating or decelerating circular cylinder is used as a model problem. The time evolution of the energy transfer from the cylinder to the fluid, as the cylinder is moved from rest to some nonnegligible velocity, is clearly seen. By examining the temporal and spatial characteristics of the numerical solution, a distinction can be made between the propagating acoustic energy, the convecting energy associated with the entropy change in the fluid, and the energy contained in the local aerodynamic field. Systematic variation of the cylinder acceleration shows that the radiated acoustic energy depends strongly upon the rate of acceleration or deceleration. The computational grid has a large effect on the ratio of acoustic energy to nonphysical entropy associated energy, while the role of the explicit artificial viscosity seems to be of second order. The entropy term was nearly negligible in all cases the cylinder was started slowly.

Brentner, Kenneth S.↗

Fast Plasma Instrument for MMS: Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. The Dual Electron Spectrometer (DES) of the Fast Plasma Instrument (FPI) for MMS meets these demanding requirements by acquiring the electron velocity distribution functions (VDFs) for the full sky with high-resolution angular measurements every 30 ms. This will provide unprecedented access to electron scale dynamics within the reconnection diffusion region. The DES consists of eight half-top-hat energy analyzers. Each analyzer has a 6 deg. x 11.25 deg. Full-sky coverage is achieved by electrostatically stepping the FOV of each of the eight sensors through four discrete deflection look directions. Data compression and burst memory management will provide approximately 30 minutes of high time resolution data during each orbit of the four MMS spacecraft. Each spacecraft will intelligently downlink the data sequences that contain the greatest amount of temporal structure. Here we present the results of a simulation of the DES analyzer measurements, data compression and decompression, as well as ground-based analysis using as a seed re-processed Cluster/PEACE electron measurements. The Cluster/PEACE electron measurements have been reprocessed through virtual DES analyzers with their proper geometrical, energy, and timing scale factors and re-mapped via interpolation to the DES angular and energy phase-space sampling measurements. The results of the simulated DES measurements are analyzed and the full moments of the simulated VDFs are compared with those obtained from the Cluster/PEACE spectrometer using a standard quadrature moment, a newly implemented spectral spherical harmonic method, and a singular value decomposition method. Our preliminary moment calculations show a remarkable agreement within the uncertainties of the measurements, with the results obtained by the Cluster/PEACE electron spectrometers. The data analyzed was selected because it represented a potential reconnection event as currently published.

Figueroa-Vinas, Adolfo↗

Helioseismic Measurements of Convective Power in Solar Cycle 24

Constraining the parameters under which convection in the solar interior operates has important implications for describing how energy is transported by plasma motions, and various physical models have been employed to provide some theoretical estimates on the expected power spectrum. Past attempts to measure the convective power distributed among large spatial scales have, however, found differing and incompatible values. Here, we present measurements of the convective power spectrum in the upper convection zone for Carrington rotations in Solar Cycle 24 obtained from the helioseismic signal corresponding to East-West flows. We also perform calibration on synthetic data using the global acoustic GALE code to make assessments of the flow velocities at various length scales without the need for performing inversions. This allows us to derive the convective power spectrum from the flow maps and to compare with the results from the raw travel times. These results are compared against predictions of convective power in global models of convection produced by the EULAG code. Finally, we show how the steps in our analysis procedure (for example data segmentation, filtering, etc.) affect our estimates of the convective power by comparing with the synthetic data from the GALE code.

Heliophysics↗

A method for estimating light quenching in inorganic scintillator detectors for radioactive ion beam experiments

In recent experiments, inorganic scintillators have been used to study the decays of exotic nuclei, providing an alternative to silicon detectors and enabling measurements that were previously impossible. However, proper use of these materials requires us to understand and quantify the scintillation process, specifically in response to very heavy nuclei. Here, in this work, we show a simplified method based on the models of Birks (1951) and Meyer and Murray (1962) to parametrize the light output of inorganic scintillators in response to beams of energetic heavy ions over a broad range of energies. We test the accuracy of our parametrization approach by calculating light output and quenching factors for various ions and comparing them with experimental data from Lutetium Yttrium Orthosilicate (LYSO:Ce), a common inorganic scintillator. The Meyer–Murray model suggests that, for sufficiently heavy ions at high energies, the majority of the light output is associated with the creation of delta electrons, which are induced by the passage of the beam through the material. These delta electrons dramatically impact the response of detection systems when subject to ions with velocities typical of beams in modern fragmentation facilities. To illustrate this, we also present a qualitative estimate of the effects of delta rays on overall light output using the Birks–Meyer–Murray parametrization. The approach presented herein will serve as a basic framework for further, more rigorous studies of scintillator response to heavy ions. This work is a crucial first step in planning future experiments where energetic exotic nuclei are interacting with scintillator detectors.

Heavy ion↗

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The effect of shearing strain-rate on the ultimate shearing resistance of clay

The shearing resistance of cohesive soils subjected to a high rate of shearing strain was investigated. A fast step-loading torque apparatus was used to induce a state of pure shear in a hollow cylindrical soil specimen. The relationship between shearing resistance and rate of shear deformation was established for various soil densities expressed in terms of initial void ratio or water content. For rate of shearing deformation studies to date, the shearing resistance increases initially with shearing velocity but subsequently reaches a terminal value as the shearing velocity increases. The terminal shearing resistance is also found to increase as the density of the soil increases. The results are useful in the rheological study of clay. It is particularly important for mobility problems of soil runways, since the soil resistance is found to be sensitive to the rate of shearing.

Cheng, R. Y. K.↗

A study of the second and third order closure models of turbulence for prediction of separated shear flows

The hybrid model of the Reynolds-stress turbulence closure is tested for the computation of the flows over a step and disk. Here it is attempted to improve the redistributive action of the turbulence energy among the Reynolds stresses. By evaluating the existing models for the pressure-strain correlation, better coefficients are obtained for the prediction of separating shear flows. Furthermore, the diffusion rate of the Reynolds stresses is reevaluated adopting several algebraic correlations for the triple-velocity products. The models of Cormack et al., Daly-Harlow, Hanjalic-Launder, and Shir were tested for the reattaching shear flows. It was generally observed that all these algebraic models give considerably low values of the triple-velocity products. This is attributed to the fact that none of the algebraic models can take the convective effect of the triple-velocity products into account in the separating shear flows, thus resulting in much lower diffusion rate than Reynolds stresses. In order to improve the evaluation of these quantities correction factors are introduced based on the comparison with some experimental data.

Amano, R. S.↗

Calibration of a Spacecraft Gyro Quadruplet

This work presents a new approach to gyro calibration where, in addition to being used for computing attitude that is needed in the calibration process, the gyro outputs are also used as measurements in a Kalman filter. Gyro calibration as well as calibration of other instruments occurs in two steps. In the first step, the instrument error parameters are estimated. During the second stage, those errors are continuously removed from the gyro readings. In the classical approach to gyro calibration, the gyro outputs are used to maintain or compute body orientation rather than being used as measurements in the context of filtering. In inertial navigation, for example, gyro errors cause erroneous computation of velocity and position, and then when the latter are compared to measured velocity and position, a great portion of the computed velocity and position errors can be determined. The latter errors are then fed into a Kalman filter (KF) that uses the INS error model to infer on the gyro errors. Similarly, when applying the classical approach to spacecraft (SC) altitude determination, the gyro outputs are used to compute the attitude and then attitude measurements are used to determine the attitude errors, which again using a KF, indicates what the gyro errors are. In the approach adopted in this work the gyro outputs are used as angular rate measurements and are compared to estimated angular rate measurements. However, this approach requires the knowledge of the angular rate. In the past, the estimated angular rate was computed in a rather simplistic way assuming that the rate was constant. In the present work, the estimated angular rate is derived using a KF whose input can be any kind of attitude measurement, therefore the angular rate experienced by the SC can be continuously changing, and yet a good estimate of the rate, necessary for calibration, can be obtained.

Harman, Richard R.↗

Plasma Photonics: Manipulating Light Using Plasmas (Full Technical Final Report)

In this LDRD, we made significant advances on multiple fronts on the investigation and design of plasma-based optical structures to manipulate light at extreme intensities. We have unveiled new fundamental properties of laser-plasma interaction that were leveraged to introduce transformative applications. We first investigated new ways to generate high-amplitude plasma-based optical structures, based on the auto-resonance mechanism in order to produce refractive structures in plasmas. New insights into the interaction of such plasma structures with light waves allowed us to propose new methods for polarization manipulation and control for high-intensity lasers in the nonlinear regime, i.e. allowing the manipulation of an intense laser via its interaction with a less intense “control” laser beam in plasmas. Another application that we demonstrated experimentally is the capability to manipulate the velocity of light waves in plasmas by tailoring its refractive index; we reported the first demonstration of “slow light” and supra-luminal light propagation in plasmas. We also developed a new technological concept that would allow the probing of the type of plasma structures investigated in such applications, using high-bandwidth laser probing [Ludwig Phys. Plasmas 2019]. Finally, first steps were taken towards the demonstration of a plasma-based amplifier and compressor for the National Ignition Facility, in order to achieve powers equal of greater than those delivered by the Advanced Radiographic Capability (ARC). This technique uses a plasma-based grating acting like a compressor and amplifier, which could potentially compress and amplify an 88 ps, 56 J beam from the National Ignition Facility to <10 ps and hundreds of J.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of unsteady airfoil loads with and without flap deflection at -90 degrees incidence

A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This unique method provides for the direct solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body-fitted computational mesh incorporating a staggered grid method. The vorticity is determined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and for the conservation of mass at the mesh-cell centers. The method provides for the direct solution of the flow field and satisfies the conservation of mass to machine zero at each time-step. The results of the present analysis and experimental results obtained for a XV-15 airfoil are compared. The comparisons indicate that the calculated drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results. Comparisons of the numerical results of the present method for several airfoils demonstrate the significant influence of airfoil curvature and flap deflection on the predicted download.

Stremel, Paul M.↗

Validation of a Computational Fluid Dynamics (CFD) Code for Supersonic Axisymmetric Base Flow

The ability to accurately and efficiently calculate the flow structure in the base region of bodies of revolution in supersonic flight is a significant step in CFD code validation for applications ranging from base heating for rockets to drag for protectives. The FDNS code is used to compute such a flow and the results are compared to benchmark quality experimental data. Flowfield calculations are presented for a cylindrical afterbody at M = 2.46 and angle of attack a = O. Grid independent solutions are compared to mean velocity profiles in the separated wake area and downstream of the reattachment point. Additionally, quantities such as turbulent kinetic energy and shear layer growth rates are compared to the data. Finally, the computed base pressures are compared to the measured values. An effort is made to elucidate the role of turbulence models in the flowfield predictions. The level of turbulent eddy viscosity, and its origin, are used to contrast the various turbulence models and compare the results to the experimental data.

Tucker, P. Kevin↗

Rotor blade boundary layer measurement hardware feasibility demonstration

A traverse mechanism which allows the measurement of the three dimensional boundary layers on a helicopter rotor blade has been built and tested on a full scale rotor to full scale conditions producing centrifugal accelerations in excess of 400 g and Mach numbers of 0.6 and above. Boundary layer velocity profiles have been measured over a range of rotor speeds and blade collective pitch angles. A pressure scanning switch and transducer were also tested on the full scale rotor and found to be insensitive to centrifugal effects within the normal main rotor operating range. The demonstration of the capability to measure boundary layer behavior on helicopter rotor blades represents the first step toward obtaining, in the rotating system, data of a quality comparable to that already existing for flows in the fixed system.

Clark, D. R.↗

Retarding potential analyzer for the Pioneer-Venus Orbiter Mission

The retarding potential analyzer on the Pioneer-Venus Orbiter Mission has been designed to measure most of the thermal plasma parameters within and near the Venusian ionosphere. Parameters include total ion concentration, concentrations of the more abundant ions, ion temperatures, ion drift velocity, electron temperature, and low-energy (0-50 eV) electron distribution function. To accomplish these measurements on a spinning vehicle with a small telemetry bit rate, several functions, including decision functions not previously used in RPA's, have been developed and incorporated into this instrument. The more significant functions include automatic electrometer ranging with background current compensation; digital, quadratic retarding potential step generation for the ion and low-energy electron scans; a current sampling interval of 2 ms throughout all scans; digital logic inflection point detection and data selection; and automatic ram direction detection. Extensive numerical simulation and plasma chamber tests have been conducted to verify adequacy of the design for the Pioneer Mission.

Knudsen, W. C.↗

Vapor concentration measurement with photothermal deflectometry

Theoretical and experimental results for using the photothermal deflection technique to measure vapor species concentration, while minimizing the disturbance of the transport (material) parameters due to vapor heating, are developed and described. In contrast to common practice, the above constraints require using a pump-beam duty cycle of less than 50 percent. The theoretical description of the shortened heating time is based on a step-function formulation of the pumping cycle. The results are obtained as closed-form solutions of the energy equation for many chopping cycles until steady state is reached, by use of a Green's-function method. The Euler formulation of the Fermat principle is used to calculate the deflection angle. The equations are expanded to include the effects of vapor velocity on both the temperature and temperature gradient profiles. The effects of finite (unfocused) pump and probe beams and thermal (Soret) diffusion are also accounted for. Excellent agreement between theory and experiment is obtained.

Banish, R. Michael↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗