Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Velocity Model”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2017

CFD-DEM modeling the effect of column size and bed height on minimum fluidization velocity in micro fluidized beds with Geldart B particles

The fluidization behavior of Geldart B particles in micro fluidized beds is investigated numerically using Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) available in the open-source Multiphase Flow with Interphase eXchanges (MFIX) code. The effects of different bed inner diameters (D) of 8 mm, 12 mm, 16 mm and various initial static bed heights (H) were examined. It is found that both decreasing the column diameter and increasing the bed height in a micro fluidized bed increases the minimum fluidization velocity (Umf). The observed overshoot in pressure drop that occurs before the onset of fluidization decreases in magnitude with increasing column diameter, however there is less sensitivity to bed height. Overall, the numerical results agree qualitatively with existing theoretical correlations and experimental studies. The simulations show that both column diameter and particle-wall friction contribute to the variation in minimum fluidization velocity. Finally, these two factors are coupled and hard to separate. The detailed influences of wall friction on minimum fluidization velocity are then investigated for a prescribed column diameter of 8 mm by varying the wall friction from 0 to 0.4.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Information-theoretic astrophysical uncertainties in the effective theory of dark matter direct detection

The impact of astrophysical uncertainties in direct detection searches can vary significantly across particle dark matter models and detector targets, due to the different velocity and momentum dependencies of the scattering cross section. We address these uncertainties for all operators of the nonrelativistic effective field theory of dark-matter/nucleon interactions, making use of the Kullback-Leibler (KL) information divergence to measure the deviation of the true dark matter velocity distribution from the Maxwell-Boltzmann form. This approach quantifies how astrophysical uncertainties affect each operator in the effective theory, without assuming any specific functional form for the velocity distribution. While for some operators the uncertainties are smaller than 1 order of magnitude for entropically motivated deviations from the Maxwell-Boltzmann form, for other operators, these uncertainties can be as large as three orders of magnitude near threshold. Furthermore, we identify the dependence of the scattering rate for various operators of the effective theory with different velocity-weighted moments of the velocity distribution, functionally analogous to the mean, variance, or skewness. This provides new analytic insight into which features of the velocity distribution are most relevant to detect a given particle dark matter model. Our technique is general and could be applied to a broader class of physics problems where a physical observable depends on the statistical moments of an uncertain theoretical distribution.

Herrera, Gonzalo [MIT, MKI; Harvard U.; Virginia T↗
DOE OSTI2017

Increase in the Amplitude of Line-of-sight Velocities of the Small-scale Motions in a Solar Filament before Eruption

We present a study on the evolution of the small-scale velocity field in a solar filament as it approaches the eruption. The observation was carried out by the Solar Dynamics Doppler Imager (SDDI) that was newly installed on the Solar Magnetic Activity Research Telescope at Hida Observatory. The SDDI obtains a narrowband full-disk image of the Sun at 73 channels from H α − 9.0 Å to H α + 9.0 Å, allowing us to study the line-of-sight (LOS) velocity of the filament before and during the eruption. The observed filament is a quiescent filament that erupted on 2016 November 5. We derived the LOS velocity at each pixel in the filament using the Becker’s cloud model, and made the histograms of the LOS velocity at each time. The standard deviation of the LOS velocity distribution can be regarded as a measure for the amplitude of the small-scale motion in the filament. We found that the standard deviation on the previous day of the eruption was mostly constant around 2–3 km s{sup −1}, and it slightly increased to 3–4 km s{sup −1} on the day of the eruption. It shows a further increase, with a rate of 1.1 m s{sup −2}, about three hours before eruption, and another increase, with a rate of 2.8 m s{sup −2}, about an hour before eruption. From this result we suggest that the increase in the amplitude of the small-scale motions in a filament can be regarded as a precursor of the eruption.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

A relationship between peak temperature drop and velocity differential in a microburst

Results from numerical microburst simulations using the Terminal Area Simulation System (Proctor, 1987) are used to develop a relationship between wind velocity differential and peak temperature drop. The numerical model and the relationships derived from the model are described. The relationship between peak temperature drop and differential wind velocity is shown to be valid during microburst development, for all precipitation shaft intensities and diameters. It is found that the relationship is not valid for low-reflectivity microburst events or in the presence of ground-based stable layers. The use of the relationship in IR wind shear detection systems is considered.

Proctor, Fred H.↗
NASA NTRS1994

Transversely accelerated ions in the topside ionosphere

Data from the rocket campaigns Mechanism in the Auroral Region for Ion Energization (MARIE) and TOpside Probe of the Auroral Zone (TOPAZ) III, within regions of low-altitude transversely accelerated ions, are interpreted to explain the acceleration of the ions. Using the Monte Carlo kinetic technique to evaluate the ion heating produced by the simultaneously observed lower hybrid waves, we find that their observed electric field amplitudes are sufficient to explain the observed ion energies in the MARIE event. Much of the uncertainty in evaluating the efficiency of a plasma wave induced particle heating process which is dependent on a velocity resonance comes from the lack of information on the phase velocities of the waves. In the case of the MARIE observations, our modeling efforts show that features in the ion velocity distribution are consistent with the wave phase velocities inferred from interferometer measurements of wavelengths. The lower hybrid waves with which low-altitude transversely accelerated ions are associated are frequently observed to be concentrated in small-scale wave packets called 'spikelets'. We demonstrate through the scaling of the size of these wave packets that they are consistent with the theory of lower hybrid collapse. Using the Monte Carlo technique, we find that if the lower hybrid field energy is concentrated in these wave packets, it is still adequate to accelerate the ionospheric ions to the observed energies.

Retterer, John M.↗
NASA NTRS2017

Planetary Probe Entry Atmosphere Estimation Using Synthetic Air Data System

This paper develops an atmospheric state estimator based on inertial acceleration and angular rate measurements combined with an assumed vehicle aerodynamic model. The approach utilizes the full navigation state of the vehicle (position, velocity, and attitude) to recast the vehicle aerodynamic model to be a function solely of the atmospheric state (density, pressure, and winds). Force and moment measurements are based on vehicle sensed accelerations and angular rates. These measurements are combined with an aerodynamic model and a Kalman-Schmidt filter to estimate the atmospheric conditions. The new method is applied to data from the Mars Science Laboratory mission, which landed the Curiosity rover on the surface of Mars in August 2012. The results of the new estimation algorithm are compared with results from a Flush Air Data Sensing algorithm based on onboard pressure measurements on the vehicle forebody. The comparison indicates that the new proposed estimation method provides estimates consistent with the air data measurements, without the use of pressure measurements. Implications for future missions such as the Mars 2020 entry capsule are described.

Karlgaard, Chris↗
NASA NTRS2017

The Velocity and Density Distribution of Earth-Intersecting Meteoroids: Implications for Environment Models

Meteoroids are known to damage spacecraft: they can crater or puncture components, disturb a spacecraft's attitude, and potentially create secondary electrical effects. Because the damage done depends on the speed, size, density, and direction of the impactor, accurate environment models are critical for mitigating meteoroid-related risks. Yet because meteoroid properties are derived from indirect observations such as meteors and impact craters, many characteristics of the meteoroid environment are uncertain. In this work, we present recent improvements to the meteoroid speed and density distributions. Our speed distribution is derived from observations made by the Canadian Meteor Orbit Radar. These observations are de-biased using modern descriptions of the ionization efficiency. Our approach yields a slower meteoroid population than previous analyses (see Fig. 1 for an example) and we compute the uncertainties associated with our derived distribution. We adopt a higher fidelity density distribution than that used by many older models. In our distribution, meteoroids with TJ less than 2 are assigned to a low-density population, while those with TJ greater than 2 have higher densities (see Fig. 2). This division and the distributions themselves are derived from the densities reported by Kikwaya et al. These changes have implications for the environment: for instance, the helion/antihelion sporadic sources have lower speeds than the apex and toroidal sources and originate from high-T(sub J) parent bodies. Our on-average slower and denser distributions thus imply that the helion and antihelion sources dominate the meteoroid environment even more completely than previously thought. Finally, for a given near-Earth meteoroid cratering rate, a slower meteoroid population produces a comparatively higher rate of satellite attitude disturbances.

Moorhead, A. V.↗
DOE OSTI2016

Development of an Integrated Hydrological Model for Simulation of Surface Runoff and Stream Flow in Tims Branch Watershed, Savannah River Site, South Carolina - 16203

This research is part of continued efforts to develop an integrated transport model to predict the distribution of tin within the overland flow in Tims Branch Watershed (TBW) at Savannah River Site (SRS), South Carolina. Tin was introduced into TBW during the application of an innovative mercury remediation technology implemented by the U.S. Department of Energy's Office of Environmental Management, which involved the injection of stannous (tin) chloride into mercury contaminated groundwater. Understanding the fate and transport of tin and its compounds is of primary importance due to their potential impact on the environment. In this study, we developed an overland hydrology model (MIKE SHE) capable of simulating surface hydrology throughout the TBW. The modeling application used historical precipitation, groundwater levels, geological data, and river discharges that were retrieved from government databases and input to the model. The model was developed to simulate spatio-temporal distribution of flow discharges, flow duration, and water levels in the TBW. The preliminary simulation results indicated that the model was capable of predicting flow depth and velocity within the study area during extreme climate events. The developed hydrology model will later on be coupled with transport model to simulate tin distribution within Tims Branch stream during seasonal storm events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Modeling Dense Particle Flow in Multistage and Obstructed Flow Receivers Using High Fidelity Simulations

Particles are a leading contender for next-generation, concentrating solar power technologies, and the design of the particle receiver is critical to minimize the levelized cost of electricity. Falling particle receivers (FPRs) are a viable receiver concept, but many new designs feature complex particle obstructions that include dense discrete phase flows. This creates additional challenges for modeling as particle-to-particle interactions (i.e., collisions) and particle drag become more complex. To improve upon existing modeling strategies, a CFD-DEM simulation capability was created by coupling two independent codes: Sierra/Fuego and LAMMPS. A suitable receiver model was then defined using a traditional continuum-based model for the air and a granular model for the particle curtain. A sensitivity study was executed using this model to determine the relevance of different granular model inputs on important quantities of interest in obstructed flow FPRs: the particle velocity and curtain opacity. The study showed that the granular model inputs had little effect on the particle velocity magnitude and curtain opacity after an obstruction.

Mills, Brantley↗
DOE OSTI2016

Evaluation of interfacial area transport models for horizontal bubbly flow

This paper presents evaluation results of interfacial area transport equation (IATE) for horizontal air-water bubbly flow. The steady-state one-dimensional, one-group IATE for adiabatic air-water bubbly flow developed by Talley [1] is employed in the current evaluation study. The models evaluated include: (a) drift-flux-based closure relations for void-weighted bubble velocity and void fraction; (b) frictional pressure drop model based on the Lockhart-Martinelli approach; (c) bubble interaction mechanisms; and (d) covariance associated with asymmetric bubble distribution in horizontal two-phase flow. In order to evaluate the models, a local two-phase flow database obtained by the four-sensor conductivity probe in fourteen different test conditions is employed, which covers a wide range of flow conditions in horizontal bubbly flow regime, with superficial liquid and gas velocities ranging from 3.50-6.00 m/s and 0.09- 1.00 m/s, respectively. It is demonstrated that both the void weighted bubble velocity and void fraction are predicted well with an average percent difference of ±2.7%. The fictional pressure loss can be predicted with an average percent difference of ±1.12%. The covariance parameters calculated from the modeling approach generally compare well with those calculated from experimental data within ±20% difference. It is found that the IATE generally predicts the measured interfacial area concentration with an average percent difference of ±5.7%. Furthermore, the effect of liquid velocity on individual bubble interaction mechanisms is investigated. It is found that both the contributions of random collision and turbulent impact decrease as liquid velocity decreases. (authors)

42 ENGINEERING↗
NASA NTRS1995

Plate Motion and Crustal Deformation Estimated with Geodetic Data from the Global Positioning System

We use geodetic data taken over four years with the Global Positioning System (GPS) to estimate: (1) motion between six major plates and (2) motion relative to these plates of ten sites in plate boundary zones. The degree of consistency between geodetic velocities and rigid plates requires the (one-dimensional) standard errors in horizontal velocities to be approx. 2 mm/yr. Each of the 15 angular velocities describing motion between plate pairs that we estimate with GPS differs insignificantly from the corresponding angular velocity in global plate motion model NUVEL-1A, which averages motion over the past 3 m.y. The motion of the Pacific plate relative to both the Eurasian and North American plates is observed to be faster than predicted by NUVEL-1A, supporting the inference from Very Long B ase- line Interferometry (VLBI) that motion of the Pacific plate has speed up over the past few m.y. The Eurasia-North America pole of rotation is estimated to be north of NUVEL-1A, consistent with the independent hypothesis that the pole has recently migrated northward across northeast Asia to near the Lena River delta. Victoria, which lies above the main thrust at the Cascadia subduction zone, moves relative to the interior of the overriding plate at 30% of the velocity of the subducting plate, reinforcing the conclusion that the thrust there is locked beneath the continental shelf and slope.

Argus, Donald F.↗
NASA NTRS1985

An experimental and computational study of turbulent boundary layers in oscillating flows

An experimental-computational study of the behavior of turbulent boundary layers for oscillating air flows over a plane surface with a small favorable mean pressure gradient is described. Experimental data were obtained for boundary layers generated on the test section wall of a new type of facility that produces a mean flow velocity with a superposed sinusoidal component over a wide frequency range. Experimental results for unsteady velocity profiles and turbulence intensity obtained from hot wire measurements at a fixed stream-wise position on the surface are presented for frequencies ranging from 5 to 21 Hz. Numerical predictions were obtained using an unsteady boundary layer computational code and the Cebeci-Smith and Glushko turbulence models. Predicted quantities related to unsteady velocity profiles exhibit fair agreement with experiment when the Cebeci-Smith model is used.

Cook, W. J.↗
DOE OSTI2017

A Semi-analytical Model for Wind-fed Black Hole High-mass X-Ray Binaries: State Transition Triggered by Magnetic Fields from the Companion Star

We propose a mechanism of state transition in wind-fed black hole (BH) binaries (high-mass X-ray binaries) such as Cyg X-1 and LMC X-1. Modeling a line-driven stellar wind from the companion by two-dimensional hydrodynamical calculations, we investigate the processes of wind capture by, and accretion onto, the BH. We assume that the wind acceleration is terminated at the He ii ionization front because ions responsible for line-driven acceleration are ionized within the front, i.e., the He iii region. It is found that the mass accretion rate inferred from the luminosity is remarkably smaller than the capture rate. Considering the difference, we construct a model for the state transition based on the accretion flow being controlled by magnetorotational instability. The outer flow is torus-like, and plays an important role to trigger the transition. The model can explain why state transition does occur in Cyg X-1, while not in LMC X-1. Cyg X-1 exhibits a relatively low luminosity, and then the He ii ionization front is located and can move between the companion and BH, depending on its ionizing photon flux. On the other hand, LMC X-1 exhibits too high luminosity for the front to move considerably; the front is too close to the companion atmosphere. The model also predicts that each state of high-soft or low-hard would last fairly long because the luminosity depends weakly on the wind velocity. In the context of the model, the state transition is triggered by a fluctuation of the magnetic field when its amplitude becomes comparable to the field strength in the torus-like outer flow.

79 ASTRONOMY AND ASTROPHYSICS↗