Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - DRAG, PROFILE - BICKNELL”

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 487 records · Page 27

Forecasts on the Dark Matter Density Profiles of Dwarf Spheroidal Galaxies with Current and Future Kinematic Observations

We forecast parameter uncertainties on the mass profile of a typical Milky Way dwarf spheroidal galaxy (dSph) using the spherical Jeans equation and Fisher matrix formalism. For a Draco-like system we show that radial velocity measurements for 1000 individual stars can constrain the mass contained within the effective radius of a dSph to within 5%. This is consistent with constraints extracted from current observational data. We compare two systems, a cusp and core, and demonstrate that a minimum sample of 100,000 (10,000) stars with both radial and proper motions measurements is required to disentangle their inner slopes at the 2σ (1σ) level. If using the log-slope measured at the half-light radius as a proxy for differentiating between a core or cusp slope, only 1000 line-of-sight and proper motions measurements are required; however, we show this choice of radius does not always unambiguously differentiate between core and cusped profiles. Once observational errors are below half the value of the intrinsic dispersion, improving the observational precision yields little change in the density profile uncertainties. The choice of priors in our profile shape analysis plays a crucial role when the number of stars in a system is less than 100 but does not affect the resulting uncertainties for larger kinematic samples. Our predicted 2D confidence regions agree well with those from a full likelihood analysis run on a mock kinematic data set taken from the Gaia Challenge, validating our Fisher predictions. Our methodology is flexible, allowing us to predict density profile uncertainties for a wide range of current and future kinematic data sets.

79 ASTRONOMY AND ASTROPHYSICS↗

Generalized wave envelope analysis of sound propagation in ducts with stepped noise source profiles and variable axial impedance

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution tends not to oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude. Example solutions indicate that stepped noise source profiles have much higher attenuation than plane waves in a uniform impedance liner. Also, multiple stepped impedance liners are shown to have higher attenuation than uniform ducts if the impedances are chosen properly. For optimum noise reduction with axial variations in impedance, the numerical analysis indicates that for a plane wave input the resistance should be near zero at the entrance of a suppressor duct, while the reactance should be near the optimum value associated with the least-attenuated mode in a uniform duct.

Baumeister, K. J.↗

Turbulent solution of the Navier-Stokes equations for an inhomogenous developing shear layer

To study the nonlinear physics of inhomogeneous turbulent shear flow, the unaveraged Navier-Stokes equations are solved numerically. For initial conditions a three-dimensional cosine velocity fluctuation and a mean-velocity profile with a step are used. Although the initial conditions are nonrandom. The flow soon becomes turbulent. Concentrated turbulent energy develops near the plane where the mean velocity gradient is initially infinite. The terms in the one-point correlation equation for turbulent energy, including those for the diffusion and production of turbulence, are calculated, the diffusion terms tend to make the turbulence more homogeneous.

Deissler, R. G.↗

Model for rotating and nonuniform planetary exospheres.

The model neutral exosphere with a uniformly rotating exobase is generalized by allowing variations in exobase density and temperature which characterize the thermosphere just below the base. The corresponding velocity distribution function, satisfying the collisionless Boltzmann equation, is constructed and used to form a general expression for the velocity moments. Resulting density profiles of rotating exospheres with nonuniform densities and temperatures on the exobase are compared with corresponding nonrotating exospheres. Density enhancements due to rotation are found to be greatest above regions of exobase density or temperature minima. Equatorial density enhancements of terrestrial hydrogen, resulting from rotation, are estimated to be 15%-17% at altitudes of 10-20 earth radii. Corresponding increases in terrestrial helium are 30%-50% on the equator at altitudes of 0.7 to 1 earth radii even when there is a polar density bulge in the barosphere.

Hartle, R. E.↗

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resilient Operating Constraints for Power Distribution Systems under Setpoint Attacks

Integration and operation of distributed generation (DG) and energy storage (ES) in power distribution systems are enabled by communication networks and embedded sensor and control devices that increase the vulnerability of the systems to cyber-threats, broadening the attack surface and making adversary actions more unpredictable. This paper proposes a methodology that uses ellipsoidal approximations to quantify the potential damage caused by successful attacks that affect, directly or indirectly, the desired operation setpoints and may drive the power distribution operation to unsafe states by violating the limits of voltage or line flows. More specifically, a new methodology is introduced to find the optimal non-symmetric operating constraints that can be imposed to each DG and ES in order to guarantee that the power distribution system is resilient to any malicious setpoints. The proposed method takes as inputs the system topology, DG and ES capabilities, and load limits to solve a convex optimization problem formulated using linear matrix inequalities (LMIs) and the power flow equations. The proposed solution is agnostic to the attacker's action or load profile and it does not require any assumption about the location or means of the attack. The numerical results on a test distribution feeder with several DG and ES illustrate how the proposed resilient operating constraints guarantee the security of the power distribution system under setpoint attacks.

Giraldo, Jairo↗

Performance-Portable GPU Acceleration of the EFIT Tokamak Plasma Equilibrium Reconstruction Code

This paper presents the steps followed to GPU-offload parts of the core solver of EFIT-AI, an equilibrium reconstruction code suitable for tokamak experiments and burning plasmas. For this work, we will focus on the fitting procedure that consists of a Grad–Shafranov (GS) equation inverse solver that calculates equilibrium reconstructions on a grid. We will show profiling results of the original code (CPU-baseline), as well as the directives used to GPU-offload the most time-consuming function, initially to compare OpenACC and OpenMP on NVIDIA and AMD GPUs and later on to assess OpenMP performance portability on NVIDIA, AMD and Intel GPUs. We will make a performance comparison for different spatial grid sizes and show the speedup achieved on NVIDIA A100 (Perlmutter-NERSC), AMD MI250X (Frontier-OLCF) and Intel PVC GPUs (Sunspot-ALCF). Finally, we will draw some conclusions and recommendations to achieve high-performance portability for an equilibrium reconstruction code on the new HPC architectures

GPU↗

Theoretical equations of state

Theoretical equations of state in geophysics, considering systematics approach to laboratory data, seismic velocity profiles, finite strain and atomistic approach

Sammis, C. G.↗

Carbon vaporization into a nonequilibrium, stagnation-point boundary layer

The heat transfer to the stagnation point of an ablating carbonaceous heat shield, where both the gas-phase boundary layer and the heterogeneous surface reactions are not in chemical equilibrium, is examined. Specifically, the nonequilibrium changes in the mass fraction profiles of carbon species calculated for frozen flow are studied. A set of equations describing the steady-state, nonequilibrium laminar boundary layer in the axisymmetric stagnation region, over an ablating graphite surface, is solved, with allowance for the effects of finite rate of carbon vaporization.

Suzuki, T.↗

Prediction of boundary-layer characteristics of an oscillating airfoil

The evolution of unsteady boundary layers on oscillating airfoils is investigated by solving the governing equations by the Characteristic Box scheme. The difficulties associated with computing the first profile on a given time line, and the velocity profiles with partial flow reversal are solved. A sample calculation is performed for an external velocity distribution typical of those found near the leading edge of thin airfoils. The viability of the calculation procedure is demonstrated.

Cebeci, T.↗

Flight in low-level wind shear

Results of studies of wind shear hazard to aircraft operation are summarized. Existing wind shear profiles currently used in computer and flight simulator studies are reviewed. The governing equations of motion for an aircraft are derived incorporating the variable wind effects. Quantitative discussions of the effects of wind shear on aircraft performance are presented. These are followed by a review of mathematical solutions to both the linear and nonlinear forms of the governing equations. Solutions with and without control laws are presented. The application of detailed analysis to develop warning and detection systems based on Doppler radar measuring wind speed along the flight path is given. A number of flight path deterioration parameters are defined and evaluated. Comparison of computer-predicted flight paths with those measured in a manned flight simulator is made. Some proposed airborne and ground-based wind shear hazard warning and detection systems are reviewed. The advantages and disadvantages of both types of systems are discussed.

Frost, W.↗

Prediction of boundary-layer characteristics of an oscillating airfoil

The evolution of unsteady boundary layers on oscillating airfoils is investigated by solving the governing equations by the Characteristic Box scheme. The difficulties associated with computing the first profile on a given time line, and the velocity profiles with partial flow reversal are solved. A sample calculation is performed for an external velocity distribution typical of those found near the leading edge of thin airfoils. The viability of the calculation procedure is demonstrated. Previously announced in STAR as N81-28391

Cebeci, T.↗

Evaluation and development of new VAS remote sensing algorithms

An algorithm, developed to permit the simultaneous retrieval of temperature/moisture profiles and surface skin temperature by direct analytical solution of the radiative transfer equation, alleviates the problem associated with the interdependencies of water vapor retrieval on temperature retrieval and their dependencies on surface emissions. Simultaneous solution of all these quantities is achieved in one calculation using the available radiance observations. Since only a single matrix inversion is required for the specification of all parameters, the solution is computationally efficient. Ancillary observations of temperature and/or moisture from surface sensors on aircraft can be readily incorporated into the solution. Simulation tests of the method indicate improved performance over the previous iterative technique, particularly for the lower troposphere and for water vapor.

Smith, W. L.↗

Electrotransport and diffusivity of molybdenum, rhenium, tungsten, and zirconium in beta-thorium

The electric mobilities, diffusivities, and effective valences were determined for molybdenum, rhenium, tungsten, and zirconium in beta-thorium. All four solutes migrated in the same direction as the electron flow. Rhenium and molybdenum were found to be very mobile, with tungsten somewhat slower. Zirconium was found to move at a rate near that of the self-diffusion of beta-thorium, viz., about 10 to the -11th sq m/s at 1500 C. The electromigration velocities showed a similar trend. A comparison was made between experimental data obtained by scanning laser mass spectrometry and theoretical transport equations for two purification experiments. Good agreement was obtained with both the concentration profile predicted by DeGroot and the purification ratio predicted by Verhoeven.

Schmidt, F. A.↗

Infrared observations of the Uranian system

Initial results are presented from analyses of IR interferometer spectrometer data collected by Voyager 2 during passes of Uranus and the moons Miranda and Ariel. The data covered reflected solar radiation in the visible and near-IR and thermal emissions between 25 and 50 microns. An equatorial atmospheric He mole fraction of 0.10-0.20 and a mass fraction of 0.16-0.34 was observed for Uranus. Para-hydrogen was also a significant fraction Uranus CH4 atmosphere. Average polar and equatorial spectra were employed to generate vertical temperature profiles between 60-900 mbar, with the finding that the temperatures at the equator and the poles are markedly similar between 400-900 mbar. The effective temperature of Uranus had a calculated upper limit of about 59.4 K, while the subsolar temperatures of Miranda and Arial were around 86 and 84 K, respectively. The albedos of the two satellites indicated surface microstructures composed of isotopically scattering grains.

Hanel, R.↗

Marine gravity of the Southern Ocean and Antarctic margin from Geosat

Geosat altimeter data, collected from an orbit with a ground rack that repeated every 17 days and overlayed one of the 17-day Seasat ground tracks, were used to map the gravity field of the Southern Ocean and the continental margin of Antarctica. The combination of ascending an descending profiles produced a typical Geosat ground track spacing of 70 km at the equator, with the best coverage occurring between the latitudes of 60 and 72 deg in both the Northern and Southern hemispheres. The new data reveal many previously uncharted seamounts and fracture zones in the extreme Southern Ocean areas adjacent to Antarctica, showing the detailed gravity signatures of the passive and active continental margins of Antarctica. Seven large age-offset fracture zones apparent in the Geosat data record the early breakup of Gondwana.

Sandwell, David T.↗