Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 127 records · Page 7

Calculation of three-dimensional compressible turbulent boundary layers

From comparisons of high speed data with low speed closure procedures using variable mean density, there does not appear to be any appreciable influence of compressibility upon turbulent shear stress modeling in compressible turbulent boundary layers, even for extreme cases such as Mach 14 to 20 with a change in density across the layer of up to a factor of 100. Other evidence of apparent lack of compressibility caused new physics which may alter the shear stress for the compressible boundary layer cases including: (1) fluctuation Mach number was generally less than 1; (2) the shear stress distribution through the boundary layer was not a function of Mach number for zero pressure gradient flows; (3) the Morkovin hypothesis was valid up to Mach 5 (based on fluctuation data); (4) profile N power was not a function of Mach number, at least up to Mach 10; and (5) the nondimensional burst period was approximately the same as that for low speed.

Source record↗

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗

Putting Priors in Mixture Density Mercer Kernels

This paper presents a new methodology for automatic knowledge driven data mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly infinite dimensional feature space. We describe a new method called Mixture Density Mercer Kernels to learn kernel function directly from data, rather than using predefined kernels. These data adaptive kernels can en- code prior knowledge in the kernel using a Bayesian formulation, thus allowing for physical information to be encoded in the model. We compare the results with existing algorithms on data from the Sloan Digital Sky Survey (SDSS). The code for these experiments has been generated with the AUTOBAYES tool, which automatically generates efficient and documented C/C++ code from abstract statistical model specifications. The core of the system is a schema library which contains template for learning and knowledge discovery algorithms like different versions of EM, or numeric optimization methods like conjugate gradient methods. The template instantiation is supported by symbolic- algebraic computations, which allows AUTOBAYES to find closed-form solutions and, where possible, to integrate them into the code. The results show that the Mixture Density Mercer-Kernel described here outperforms tree-based classification in distinguishing high-redshift galaxies from low- redshift galaxies by approximately 16% on test data, bagged trees by approximately 7%, and bagged trees built on a much larger sample of data by approximately 2%.

Srivastava, Ashok N.↗

Influence of suspension on the oxidative burst by rat neutrophils

The influence of spaceflight on the oxidative burst of neutrophils is not known. The present study was designed to evaluate the influence of antiorthostatic suspension, a ground-based modeling system designed to simulate certain aspects of weightlessness that occur after spaceflight, on the capacity of rat neutrophils to express the oxidative burst, an important host defense mechanism against microbial pathogens. Rats were suspended in whole body harnesses in the antiorthostatic orientation for a 3- or 7-day period. Control rats were suspended orthostatically or allowed to remain in vivarium cages without the attachment of any suspension materials. After suspension, peripheral blood was harvested and neutrophils were isolated by density gradient centrifugation. The enriched neutrophil preparations were stimulated with N-formyl-methionyl-leucine-phenylalanine and phorbol myristic acid to induce the oxidative burst. It was found that neutrophils isolated from suspended animals released the same levels of superoxide anion as did vivarium control animals that were not suspended, indicating that whole body suspension did not alter this aspect of rat neutrophil function.

Non-NASA Center↗

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

Transport by weak electrostatic density drift fluctuations

Expressions are derived for transport by weak electrostatic microfluctuations driven unstable by a density gradient perpendicular to a uniform magnetic field in a Vlasov plasma. These expressions, which are evaluated for the universal and lower-hybrid density drift instabilities, provide improved physical interpretations for various wave-particle exchange frequencies previously reported. It is demonstrated that wave-particle effects generally reduce the free energy driving these instabilities, and that some forms of steady-state distribution functions are inappropriate for nonlinear theory.

Sanderson, J. J.↗

Nature of convection-stabilized dc arcs in dual-flow nozzle geometry. I - The cold flow field and dc arc characteristics. II - Optical diagnostics and theory

Steady-state low-current air arcs in a dual-flow nozzle system are studied experimentally. The cold flow field with no arc is investigated using a 12.7-mm diameter dual-flow nozzle in a steady-flow facility. Mach number and mass flux distributions are determined for various nozzle-pressure ratios and nozzle-gap spacing. It is found that the shock waves in the converging-diverging nozzles result in a decrease in overal resistance by about 15 percent. Also, Schlieren and differential interferometry techniques are used to visualize the density gradients within the arc plasma and thermal mantle. Both optical techniques reveal a laminar arc structure for a reservoir pressure of 1 atm at various current levels. Experimentally determined axial static pressure and cold-flow mass flux rate distributions and a channel-flow model with constant arc temperatre are used to solve the energy integral for the arc radius as a function of axial distance. The arc electric field strength, voltage, resistance, and power are determined with Ohm's law and the total heat transfer is related to arc power.

Serbetci, Ilter↗

Early Spring Post-Fire Snow Albedo Dynamics in High Latitude Boreal Forests Using Landsat-8 OLI Data

Taking advantage of the improved radiometric resolution of Landsat-8 OLI which, unlike previous Landsat sensors, does not saturate over snow, the progress of fire recovery progress at the landscape scale (less than 100 m) is examined. High quality Landsat-8 albedo retrievals can now capture the true reflective and layered character of snow cover over a full range of land surface conditions and vegetation densities. This new capability particularly improves the assessment of post-fire vegetation dynamics across low- to high-burn severity gradients in Arctic and boreal regions in the early spring, when the albedos during recovery show the greatest variation. We use 30 m resolution Landsat-8 surface reflectances with concurrent coarser resolution (500 m) MODIS high quality full inversion surface Bidirectional Reflectance Distribution Functions (BRDF) products to produce higher resolution values of surface albedo. The high resolution full expression shortwave blue sky albedo product performs well with an overall RMSE of 0.0267 between tower and satellite measures under both snow-free and snow-covered conditions. While the importance of post-fire albedo recovery can be discerned from the MODIS albedo product at regional and global scales, our study addresses the particular importance of early spring post-fire albedo recovery at the landscape scale by considering the significant spatial heterogeneity of burn severity, and the impact of snow on the early spring albedo of various vegetation recovery types. We found that variations in early spring albedo within a single MODIS gridded pixel can be larger than 0.6. Since the frequency and severity of wildfires in Arctic and boreal systems is expected to increase in the coming decades, the dynamics of albedo in response to these rapid surface changes will increasingly impact the energy balance and contribute to other climate processes and physical feedback mechanisms. Surface radiation products derived from Landsat-8 data will thus play an important role in characterizing the carbon cycle and ecosystem processes of high latitude systems.

Landsat-8 snow albedo↗

Anisotropies in the fluxes of Pioneer 10 protons

One-hour-averaged fluxes of 1.8- to 2.15-MeV protons observed by the LET2 detector on Pioneer 10 on the inbound trajectory showed anisotropies attributable to corotation of Jupiter's magnetodisc only when Pioneer was near the dipole equator. Most of the time the anisotropy greatly exceeded the value expected from corotation. Gradients in the distribution function can be used to account for this excess anisotropy, but the amount of gradient required is unacceptably large by 1-2 orders of magnitude. If they were taken as real, these gradients would predict almost complete disappearance of these protons from Jupiter's magnetosphere in a matter of hours. The remedy is to introduce into the model of the distribution function proton flow along field lines away from the equator into both the southern and the northern hemisphere. The parallel flux at the southernmost latitudes reached by Pioneer can reach 25% of the product of proton density and velocity, i.e., 25% of the maximum possible.

Northrop, T. G.↗

Towards Surrogate Modeling of Subgrid Turbulent Transport for 3D Radiative Hydrodynamic Simulations of the Quiet Sun

In this work, we investigate the use of deep learn-ing techniques as surrogate models, to enhance the estimationof effects of subgrid turbulent transport for 3D radiatuve hy-drodynamic simulations of the quiet Sun. We develop two dis-tinct 3D Convolutional Neural Networks (3DCNNs) to capturespatio-temporal dependencies in 3D velocity fields, leveragingdifferent activation functions and architectural designs. Thesemodels integrate both averaged velocity vector components andscalar features such as plasma density to enhance predictionaccuracy. Additionally, a Multilayer Perceptron (MLP) modelis employed to approximate complex nonlinear relationships,offering a comparison in performance between convolutionaland fully connected architectures. Logarithmic transformationis applied to the targets to handle heavily skewed data, im-proving model performance. All models are compared againsta physics-based Gradient Model. Results show that the 3DCNNmodels excel at approximating Reynolds stress tensors, makingthem a candidate for assisting in producing reduced resolutionsimulations, and thereby reducing computational overheadwhile maintaining higher accuracy than the baseline. Thesefindings demonstrate the potential of deep learning, particu-larly CNNs, to advance scalable and accurate simulations ofsolar dynamics, offering a promising alternative to traditionalturbulence models.

Heliophysics↗

Structure and functions of water-membrane interfaces and their role in proto-biological evolution

Among the most important developments in proto-biological evolution was the emergence of membrane-like structures. These are formed by spontaneous association of relatively simple amphiphilic molecules that would have been readily available in the primordial environment. The resulting interfacial regions between water and nonpolar interior of the membrane have several properties which made them uniquely suitable for promoting subsequent evolution. They can (1) selectively attract organic material and mediate its transport, (2) serve as simple catalysts for chemical reactions, and (3) promote the formation of trans-membrane electrical and chemical gradients which could provide energy sources for proto-cells. Understanding the structure of interfaces, their interactions with organic molecules and molecular mechanisms of their functions is an essential step to understanding proto-biological evolution. In our computer simulation studies, we showed that the structure of water at interfaces with nonpolar media is significantly different from that in the bulk. In particular, the average surface dipole density points from the vapor to the liquid. As a result, negative ions can approach the interface more easily than positive ions. Amphiphilic molecules composed of hydrocarbon conjugated rings and polar substituents (e.g., phenol) assume at the interface rigid orientations in which polar groups are buried in water while hydrocarbon parts are located in the nonpolar environment. These orientational differences are of special interest in connection with the ability of some of these molecules to efficiently absorb photons. Flexible molecules with polar substituents often adopt at interfaces conformations different from those in the bulk aquaeous solution and in the gas phase. As a result, in many instances both specificity and kinetics of chemical reactions in which these molecules can participate is modified by the presence of surfaces. Of special interest is the mechanism by which polar molecules are transferred across interface between water and a nonpolar medium. Our recent study showed that simple ionophores bind ions by the same mechanisms as ion channels and carriers from modern cells.

Pohorille, A.↗

A model for simulating random atmospheres as a function of latitude, season, and time

An empirical stochastic computer model was developed with the capability of generating random thermodynamic profiles of the atmosphere below an altitude of 99 km which are characteristic of any given season, latitude, and time of day. Samples of temperature, density, and pressure profiles generated by the model are statistically similar to measured profiles in a data base of over 6000 rocket and high-altitude atmospheric soundings; that is, means and standard deviations of modeled profiles and their vertical gradients are in close agreement with data. Model-generated samples can be used for Monte Carlo simulations of aircraft or spacecraft trajectories to predict or account for the effects on a vehicle's performance of atmospheric variability. Other potential uses for the model are in simulating pollutant dispersion patterns, variations in sound propagation, and other phenomena which are dependent on atmospheric properties, and in developing data-reduction software for satellite monitoring systems.

Campbell, J. W.↗

Low-Earth Orbit Trajectory Optimization in the Presence of Atmospheric Uncertainty

The previous 20 to 25 years have seen a tremendous increase in space exploration, and with that an increase in the level of logistics planning needed to ensure mission success. For spacecraft that are designed to be periodically re-supplied, a key logistics consumable is propellant, as it constitutes the greatest up-mass on re-supply vehicles. A trajectory design strategy is therefore desired that minimizes propellant usage in order to ease the demand for propellant re-supply missions. This thesis develops such a strategy in three stages, and uses the International Space Station (ISS) as its testbed, as no other LEO spacecraft is more challenging from a space logistics standpoint. First, the ISS trajectory planning problem is formulated as a constrained burn optimization problem assuming a deterministic atmosphere. The cost function is total ∆v, with constraints imposed on longitude of ascending viii node (LAN) and semi-major axis (SMA) altitude. Analytic derivatives are constructed for both the cost and constraints, which are necessary given the 6-week to 2-year time frames being considered. A gradient-based optimizer is then utilized to find locally-optimal solutions to real-world ISS trajectory planning problems. Second, atmospheric uncertainty is addressed by constructing a probabilistic model of space weather data using Gaussian Processes (GPs). Bayesian inference is performed using the GP model to generate mean and covariance estimates for space weather predictions, whose pedigree is assessed against test data. The predictions are then mapped into atmospheric density via the analytic Jacchia-Roberts density model, and the effect of space weather uncertainty on orbital lifetime is examined. Third, an ISS burn execution uncertainty model is developed. This model, along with the space weather uncertainty model, are deployed in a linear covariance analysis to ascertain their combined effect on LAN and SMA altitude dispersions. The deterministic constraints from the original problem are re-formulated as stochastic constraints, where now the constraint uncertainty interval is required to fall within specified bounds. An updated optimization framework is constructed using the original ∆v cost function along with the stochastic constraints to solve the trajectory optimization problem under atmospheric uncertainty. Finally, the complete architecture is summarized for deployment in an operational setting.

Trajectory Optimization↗

System noise analysis of the dumbbell tethered satellite for gravity-gradient measurements

An analysis of the dumbbell gravity gradiometer concept for measuring short wavelength variations in the earth's gravity gradient is presented. Variations in the gradient are recorded by measuring tension variations in a vertically stabilized satellite consisting of heavy masses connected by a long wire or rod. Tension noise arises from the excitation of various mechanical oscillations of the system. The principal noise sources that were identified are fluctuations in atmospheric drag heating and drag force resulting from density variations and winds. Approximate analytical expressions are presented for the tension noise as a function of the system design parameters for various possible configurations. Computer simulations using numerical integration were performed to study the tension noise for several sample cases. Three designs consistent with Shuttle launch capabilities are discussed.

Colombo, G.↗

The Influence of Tungsten on the Chemical Composition of a Temporally Evolving Nanostructure of a Model Ni-Al-Cr Superalloy

The influence of W on the temporal evolution of gamma' precipitation toward equilibrium in a model Ni-Al-Cr alloy is investigated by three-dimensional atom-probe (3DAP) microscopy and transmission electron microscopy (TEM). We report on the alloys Ni-10 Al-8.5 Cr (at.%) and Ni-10 Al-8.5 Cr-2 W (at.%), which were aged isothermally in the gamma + gamma' two-phase field at 1073 K, for times ranging from 0.25 to 264 h. Spheroidal-shaped gamma' precipitates, 5-15 nm diameter, form during quenching from above the solvus temperature in both alloys at a high number density (approx. 10(exp 23/cu m). As gamma' precipitates grow with aging at 1073 K, a transition from spheriodal- to cuboidal-shaped precipitates is observed in both alloys. The elemental partitioning and spatially resolved concentration profiles across the gamma' precipitates are obtained as a function of aging time from three-dimensional atom-by-atom reconstructions. Proximity histogram concentration profiles of the quaternary alloy demonstrate that W concentration gradients exist in gamma' precipitates in the as-quenched and 0.25-h aging states, which disappear after 1 h of aging. The diffusion coefficient of W in gamma' is estimated to be 6.2 x 10(exp -20) sq m/s at 1073 K. The W addition decreases the coarsening rate constant, and leads to stronger partitioning of Al to gamma' and Cr to gamma.

Sudbrack, Chantal K.↗

Compressible boundary layer with normal pressure gradients: Quasi-similarity equations - Their properties at the wall and at sharp and blunt leading edges.

The failure of most viscous-inviscid interaction methods at strong interactions is attributed to the presence of a normal pressure gradient. A new theory is proposed for supersonic laminar boundary layers that can generate normal pressure gradients. The Navier-Stokes equations are reexamined by an order of magnitude analysis and all first and second order terms are retained. The approximation is found to be dependent not only on the boundary layer thickness but also on the ratio of the dimensionless viscosity and density. The equations are transformed into two quasi-similar, nonlinear, third order, ordinary integro-differential equations for the velocity and pressure as functions of a single transverse variable. The properties of the equations at the boundaries are discussed.

Seginer, A.↗

Solar Wind Helium Abundance as a Function of Speed and Heliographic Latitude: Variation through a Solar Cycle

We present a study of the variation of the relative abundance of helium to hydrogen in the solar wind as a function of solar wind speed and heliographic latitude over the previous solar cycle. The average values of A(sub He), the ratio of helium to hydrogen number densities, are calculated in 25 speed intervals over 27-day Carrington rotations using Faraday Cup observations from the Wind spacecraft between 1995 and 2005. The higher speed and time resolution of this study compared to an earlier work with the Wind observations has led to the discovery of three new aspects of A(sub He), modulation during solar minimum from mid-1995 to mid-1997. First, we find that for solar wind speeds between 350 and 415 km/s, A(sub He), varies with a clear six-month periodicity, with a minimum value at the heliographic equatorial plane and a typical gradient of 0.01 per degree in latitude. For the slow wind this is a 30% effect. We suggest that the latitudinal gradient may be due to an additional dependence of coronal proton flux on coronal field strength or the stability of coronal loops. Second, once the gradient is subtracted, we find that A(sub He), is a remarkably linear function of solar wind speed. Finally, we identify a vanishing speed, at which A(sub He), is zero, is 259 km/s and note that this speed corresponds to the minimum solar wind speed observed at one AU. The vanishing speed may be related to previous theoretical work in which enhancements of coronal helium lead to stagnation of the escaping proton flux. During solar maximum the A(sub He), dependences on speed and latitude disappear, and we interpret this as evidence of two source regions for slow solar wind in the ecliptic plane, one being the solar minimum streamer belt and the other likely being active regions.

Kasper, J. C.↗

Interchange Method in Compressible Magnetized Couette Flow: Magnetorotational and Magnetoconvective Instabilities

We obtain the general forms of the axisymmetric stability criteria in a magnetized compressible Couette flow using an energy variational principle, the so-called interchange or Chandrasekhar s met hod, which we applied successfully in the incompressible case. This formulation accounts for the simultaneous presence of gravity, rotation, a toroidal magnetic field, a weak axial magnetic field, entropy gradients, and density gradients in the initial equilibrium state. The power of the method lies in its simplicity which allows us to derive extremely compact and physically clear expressions for the relevant stability criteria despite the inclusion of so many physical effects. In the implementation of the method, all the applicable conservation laws are explicitly taken into account during the variations of a quantity with dimensions of energy which we call the free energy function. As in the incompressible case, the presence of an axial field invalidates the conservation laws of angular momentum and azimuthal magnetic flux and introduces instead isorotation and axial current conservation along field lines. Our results are therefore markedly different depending on whether an axial magnetic field is present, and generalize in two simple expressions all previously known, partial stability criteria for the appearance of magnetorotational instability. Furthermore, the coupling between magnetic tension and buoyancy and its influence to the dynamics of nonhomoentropic magnetized flows becomes quite clear from our results. In the limits of plane-parallel atmospheres and homoentropic flows, our formulation easily recovers the stability criteria for suppression of convective and Parker instabilities, as well as some related special cases studied over 40 years ago by Newcomb and Tserkovnikov via laborious variational techniques.

Christodoulou, Dimitris M.↗