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 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Primordial black holes arise when the inflaton falls

Primordial Black Holes (PBHs) have entered the forefront of theoretical cosmology, due their potential role in phenomena ranging from gravitational waves, to dark matter, to galaxy formation. While producing PBHs from inflationary fluctuations naively would seem to require a large deceleration of the inflaton from its velocity at the horizon exit of CMB scales, in this work we demonstrate that an acceleration from a relatively small downward step in the potential that is transited in much less than an e-fold amplifies fluctuations as well. Depending on the location of the step, such PBHs could explain dark matter or the black holes detected by the gravitational wave interferometers. Furthermore, the perturbation enhancement has a natural interpretation as particle production due to the non-adiabatic transition associated with the step.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Electrolyte-Assisted Hydrogen Cycling in Lithium and Sodium Alanates at Low Pressures and Temperatures

An investigation of electrolyte-assisted hydrogen storage reactions in complex aluminum hydrides (LiAlH4 and NaAlH4) reveals significantly reduced reaction times for hydrogen desorption and uptake in the presence of an electrolyte. LiAlH4 evolves ~7.8 wt% H2 over ~3 h in the presence of a Li-KBH4 eutectic at 130 °C compared to ~25 h for the same material without the electrolyte. Similarly, NaAlH4 exhibits 4.8 wt% H2 evolution over ~4 h in the presence of a diglyme electrolyte at 150 °C compared to 4.4 wt% in ~15 h for the same material without the electrolyte. These reduced reaction times are composed of two effects, an increase in reaction rates and a change in the reaction kinetics. While typical solid state dehydrogenation reactions exhibit kinetics with rates that continuously decrease with the extent of reaction, we find that the addition of an electrolyte results in rates that are relatively constant over the full desorption window. Fitting the kinetics to an Avrami-Erofe’ev model supports these observations. The desorption rate coefficients increase in the presence of an electrolyte, suggesting an increase in the velocities of the reactant-product interfaces. In addition, including an electrolyte increases the growth parameters, primarily for the second desorption steps, resulting in the observed relatively constant reaction rates. Similar effects occur upon hydrogen uptake in NaH/Al where the presence of an electrolyte enables hydrogenation under more practical low temperature (75 °C) and pressure (50 bar H2) conditions.

25 ENERGY STORAGE↗

NASA NTRS1973

Observations of collisionless drift waves.

Velocity distribution functions of both the electron and the proton component of the magnetospheric plasma are measured in 62 logarithmically spaced steps from 50 eV to 50 keV by an instrument carried upon the ATS-5 spacecraft in circular synchronous orbit (6.6 earth radii). Oscillations which appear to be collisionless drift waves having periods of several minutes or longer are found in this data. Linearized wave theory can be used to predict the theoretical perturbation distribution function in terms of the perturbing fields and the equilibrium distribution function. The functional dependence of the theoretical perturbation distribution function upon particle energy can then be compared with that of the measured perturbation.

Laquey, R. E.↗

NASA NTRS1976

A stop-restart solid propellant study with salt quench

Experiments were conducted to gain insight into the unsatisfactory performance of the salt quench system of solid propellants in earlier studies. Nine open-air salt spray tests were conducted and high-speed cinematographic coverage was obtained of the events. It is shown that the salt spray by the detonator is generally a two-step process yielding two different fractions. The first fraction consists of finely powdered salt and moves practically unidirectionally at a high velocity (thousand of feet per second) while the second fraction consists of coarse particles and moves randomly at a low velocity (a few feet per second). Further investigation is required to verify the speculation that a lower quench charge ratio (weight of salt/propellant burning area) than previously employed may lead to an efficient quench

Kumar, R. N.↗

NASA NTRS1982

Plastic deformation and wear process at a surface during unlubricated sliding

The plastic deformation and wear of a 304 stainless steel surface sliding against an aluminum oxide rider with a spherical surface (the radius of curvature: 1.3 cm) were observed by using scanning electron and optical microscopes. Experiments were conducted in a vacuum of one million Pa and in an environment of fifty thousandth Pa of chlorine gas at 25 C. The load was 500 grams and the sliding velocity was 0.5 centimeter per second. The deformed surface layer which accumulates and develops successively is left behind the rider, and step shaped proturbances are developed even after single pass sliding under both environmental conditions. A fully developed surface layer is gradually torn off leaving a characteristic pattern. The mechanism for tearing away of the surface layer from the contact area and sliding track contour is explained assuming the simplified process of material removal based on the adhesion theory for the wear of materials.

Yamamoto, T.↗

NASA NTRS1983

Plastic deformation and wear process at a surface during unlubricated sliding

Yamamoto, T.↗

NASA NTRS1993

Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.

Palmer, P. L.↗

DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2011

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗

NASA NTRS1998

Prediction of Unsteady Transitional Layers in Turbomachinery Using Navier Stokes Equations

The objective of the research reported in this presentation is to develop computational techniques for the prediction of unsteady transitional flows associated with the rotor stator interaction in turbomachinery. Three low-Reynolds number turbulence models are incorporated in two unsteady Navier-Stokes codes (one is pressure based and the other is time marching with Runge-Kutta time stepping) and evaluated for accuracy in predicting the onset and the end of unsteady transitional patches due to wake passing. The best model is then used for modification and improvement for the leading edge effect. An existing steady Navier-Stokes code was modified to include pseudo-time stepping, which provided acceleration from 5 to 25 times that of the original code. A systematic validation procedure was implemented to assess the effects of the grid, artificial dissipation, physical, and the pseudo-time step for an accurate prediction of transitional flows resulting from the rotor-stator interaction. The ability of the Navier-Stokes code to predict the unsteady transitional flow on a turbomachinery blade is demonstrated. The unsteady pressure and velocity fields are in good agreement with the experimental data and the prediction from the Euler/boundary layer approach. The numerical solver was able to capture all zones (wake induced transitional strip, wake induced turbulent strip, calmed region, etc.) associated with wake induced transition in a compressor cascade. Another significant step is the assessment of k-epsilon turbulence models, including the leading edge modifications. Best results were obtained from the FLB model. The LB model predicted earlier inception of the transition and shorter transition length. Modification of the k-epsilon model was found to be essential for an accurate prediction of the unsteady transitional flow in a compressor cascade. The CH model failed to predict the unsteady transitional flow. Predicted boundary layer was turbulent from the leading edge, even with the modification of the k-epsilon model near the stagnation point. A comparison of the instantaneous shape factor, the skin friction coefficient, and the momentum thickness indicates that the Navier-Stokes predictions are reasonably good. Interaction between the upstream wake and stator wake results in shedding of unsteady vortices from the trailing edge and increased dissipation in the stator wake and, as a consequence, increased rate of decay of the stator wake. The procedure developed at Penn State should aid the designers of turbomachinery in allowing for the unsteady transitional flows due to rotor-stator interaction.

Lakshminarayana, B.↗

DOE OSTI2022

Three-dimensional realizations of flood flow in large-scale rivers using the neural fuzzy-based machine-learning algorithms

Machine learning methods have been extensively used to study the dynamics of complex fluid flows. One such algorithm, known as adaptive neural fuzzy inference system (ANFIS), can generate data-driven predictions for flow fields, but has not been applied to natural geophysical flows in large-scale rivers. Herein, we demonstrate the potential of ANFIS to produce three-dimensional (3D) realizations of the instantaneous flood flow field in several large-scale, virtual meandering rivers. The 3D dynamics of flood flow in large-scale rivers were obtained using large-eddy simulation (LES). The LES results, i.e., the 3D velocity components, were employed to train the learnable coefficients of an ANFIS. Further, the trained ANFIS, along with a few time-steps of LES results (precursor data) were then used to produce 3D realizations of flood flow fields in large-scale rivers with geometries other than the one the ANFIS was trained with. We also used the trained ANFIS to generate 3D realizations of river flow at a discharge other than that the ANFIS was trained with. The flow field results obtained from ANFIS were validated using separate LES runs to assess the accuracy of the 3D instantaneous realizations of the machine learning algorithm. An error analysis was conducted to quantify the discrepancies among the ANFIS and LES results for various flood flow predictions in large-scale rivers.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗

DOE OSTI2021

A pion-argon cross section measurement in the ProtoDUNE-SP experiment with cosmogenic muon

Neutrinos are tiny mysterious fundamental particles with small cross sections. Through neutrino physics, scientists across the world are trying to answer many intriguing questions about nature such as the dominance of matter over antimatter, CP violation in the lepton sector, number of supernovas in the early universe, etc. Detection of neutrinos requires massive particle detectors and intense neutrino beam owing to their small cross section. Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment that is planned to start taking data beginning in 2026. DUNE will consist of 4 massive detectors, the first of which will be using single-phase liquid argon time projection chamber (LArTPC) technology. The ProtoDUNE-SP experiment is a prototype of the DUNE built at the CERN neutrino platform and uses the same detector technology that will be used in DUNE first module. The ProtoDUNE-SP experiment collected months of test beam and cosmic ray data beginning in September 2018. It was built to provide a testbed for the installation of detector parts for DUNE, showing long-term stability of the detector, understanding detector response for different test beam particles (including protons, pions, electrons, kaons, muons), and measurement of hadron-argon cross sections. When a particle passes through LArTPC electron-ion pairs are produced. To reconstruct the position and energy of a particle passing through the medium knowledge of ionization electron drift velocity is essential. The electron drift velocity is distorted by an excess positive charge built up in the detector, known as space charge. This study discusses a novel technique for measuring the ionization electron drift velocity using cosmic-ray muons. The technique uses tracks that travel the entire drift distance of the TPC for drift velocity determination. Secondly, the study discusses a method for converting the charge deposited into energy. The method is carried out in two step s. In th e first step detector response for energetic cosmic ray muons crossing the entire the TPC is used to make the charge deposition uniform throughout the TPC, and in the second step stopping cosmic-ray muons are used for determining the energy scale. Finally, the study discusses a pion-argon cross section measurement based on reweighting of Monte Carlo simulations using J. Calcutt's Geant4Reweight framework. Neutrinos cannot be directly detected; they are identified based on the interaction products. Pions are a common interaction product in a neutrino interaction. For precise modeling of neutrino event generators, it is essential to understand the pion-argon interaction. Pion-argon cross section measurement serves as an important input for neutrino interaction models. The results of the pion-argon total reaction cross section using the Geant4 reweighting technique are found to be in good agreement with Geant4 predictions. The many studies carried out in the ProtoDUNE-SP experi ment wil l be useful for current and future neutrino experiments using LArTPC technology including ICARUS, MicroBooNE, DUNE

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1976

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse to fine. The coarser grids are used to diminish the magnitude of the smooth part of the residuals. The method was applied to the solution of the transonic small disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is studied with meshes of both constant and variable step size.

South, J. C., Jr.↗

NASA NTRS1981

A modified dodge algorithm for the parabolized Navier-Stokes equations and compressible duct flows

A revised version of a split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three-dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard successive overrelaxation iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition.

Cooke, C. H.↗

NASA NTRS1981

A modified Dodge algorithm for the parabolized Navier-Stokes equation and compressible duct flows

A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitive agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions.

Cooke, C. H.↗

NASA NTRS1982

Plastic deformation at surface during unlubricated sliding

The plastic deformation and wear of 304 stainless-steel surface slid against an aluminum oxide rider were observed by using a scanning electron microscope and an optical microscope. Experiments were conducted in a vacuum of 0.000001 Pa and in an environment of 0.0005 Pa chlorine gas at 25 C. The load was 500 grams and the sliding velocity was 0.5 centimeter per second. The deformed surface layer which accumulates and develops successively is left behind the rider, and step-shaped protuberances are developed even after single pass sliding under both environmental conditions. A fully developed surface layer is gradually torn off leaving a characteristic pattern. These observations result from both adhesion and an adhesive wear mechanism.

Yamamoto, T.↗

NASA NTRS1983

A modified Dodge algorithm for the parabolized Navier-Stokes equations and compressible duct flows

A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitative agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions. Previously announced in STAR as N82-16363

Cooke, C. H.↗