Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 757 records · Page 42

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum exitation gap, gmin = O(n2(sup -n/2)), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, Vadius↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum excitation gap. g min, = O(n 2(exp -n/2), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to 'the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, V. N.↗

Fermi Gamma-Ray Space Telescope: Science Highlights for the First 8 Months

The Fermi Gamma-ray Space Telescope was launched on June 11, 2008 and since August 2008 has successfully been conducting routine science observations of high energy phenomena in the gamma-ray sky. A number of exciting discoveries have been made during its first year of operation, including blazar flares, high-energy gamma-ray bursts, and numerous new,gamma-ray sources of different types, among them pulsars and Active Galactic Nuclei (AGN). fermi-LAT also performed accurate mea.<;urement of the diffuse gamma-radiation which clarifies the Ge V excess reported by EGRET almost 10 years ago, high precision measurement of the high energy electron spectrum, and other observations. An overview of the observatory status and recent results as of April 30, 2009, are presented. Key words: gamma-ray astronomy, cosmic rays, gamma-ray burst, pulsar, blazar. diffuse gamma-radiation

Moiseev, Alexander↗

Numerical experiments with a stochastic zonal climate model

A zonally averaged energy balance climate model is developed to generate zonal temperature variability through fluctuating meridional energy transports. Stochastic transport fluctuations are included in the model by multiplying the eddy diffusion coefficients by Gaussian random deviates. For a variability of eddy coefficients of 50 percent, the model is found to generate an interannual temperature variability of 0.03 K for the global temperature, and 0.04 and 0.05 K for the Northern and Southern Hemispheres, respectively. It is shown that the temperature variability level generated is linearly related to the transport variability level introduced. An increase in the level of model-generated temperature variability and a change in the shape of variance spectra of temperature anomaly time series are obtained by switching from the multiplicative noise model to an additive noise model. The results of these model studies are compared with a time series of central England temperatures as well as GCM generated climate variability.

Schneider, S. H.↗

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗

Role of Hydrogen and Structure on CMSX-2

The goal was to characterize and understand the detailed effect of hydrogen on the tensile and fatigue behavior of single crystals of CMSX-2. Specimens were solutionized, air cooled, and then aged in one of two ways. The resultant macro and microstructure included a pronounced dendritic structure with numerous low angle boundaries in the interdendritic region and on the average a high level of porosity with a large variation in size and distribution. Hydrogen permeability and diffusivity in this material was found to be very low at room temperature. Hydrogen concentration studies on thin samples yielded values significantly higher than found for similar charging conditions in nickel and stainless steels. The tensile mechanical properties were found to depend on hydrogen concentration. SEM observations of the brittle subsurface zone revealed a strong correlation between crack initiation sites and voids just beneath the surface. The hollow and solid samples were compared to better understand the influence of increased charged volume on the fatigue behavior.

Baker, C. L.↗

The global evolution of the primordial solar nebula

Complete radial, time-dependent calculations of the structure and evolution of the primordial solar nebula during the viscous diffusion stage are presented. The viscous stress is derived from analytic one-zone models of the vertical nebular structure based on detailed grain opacities. Comparisons with full numerical integrations indicate that the effective viscous alpha parameter is about 0.01. The evolution time of a minimum mass nebula is one-million yr or less. The flow pattern of fluid elements in the disk is examined and the implications the results have on the theory of the formation of the solar system are discussed.

Ruden, S. P.↗

Theory of Influence of a Low-Volatility, Soluble Impurity on Spherically-Symmetric Combustion of Fuel Droplets

Analyses are given for the evolution of liquid-phase mass fraction profiles and temperature during combustion of fuel droplets composed of binary miscible mixtures of low-volatility and high-volatility constituents, with small initial mass fractions of the low-volatility material, or slow introduction of low-volatility material into the liquid phase by absorption from the gas phase. The gas phase is assumed to remain quasi-steady and the liquid temperature spatially uniform. The ratio of the liquid-phase diffusion coefficient to the initial burning-rate constant is treated as a small parameter. Asymptotic analyses in this parameter are developed, and the conservation equations are integrated numerically in obtaining descriptions of the combustion history. It is shown that at the surface of the droplet a boundary layer arises in which the mass fraction of the low-volatility component increases with time. The results are used to explain qualitatively some observed conditions of flame contraction and liquid disruption in droplet combustion.

Shaw, B. D.↗

Estimation of the processes controlling variability in phytoplankton pigment distributions on the southeastern U.S. continental shelf

Five coastal-zone-color-scanner images from the southeastern U.S. continental shelf are combined with concurrent moored current meter measurements to assess the processes controlling the variability in chlorophyll concentration and distribution in this region. An equation governing the space and time distribution of a nonconservative quantity such as chlorophyll is used in the calculations. The terms of the equation, estimated from observations, show that advective, diffusive, and local processes contribute to the plankton distributions and vary with time and location. The results from this calculation are compared with similar results obtained using a numerical physical-biological model with circulation fields derived from an optimal interpolation of the current meter observations and it is concluded that the two approaches produce different estimates of the processes controlling phytoplankton variability.

Mcclain, Charles R.↗

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Dynamical Evolution of Ring-Satellite Systems

The goal of this research was to understand dynamical processes related to the evolution of size distribution of particles in planetary rings and application of theoretical results to explain features in the present rings of giant planets. We studied velocity evolution and accretion rates of ring particles in the Roche zone. We developed a new numerical code for the evolution of ring particle size distribution, which takes into account the above results for particle velocity evolution and accretion rates. We also studied radial diffusion rate of ring particles due to inelastic collisions and gravitational encounters. Many of these results can be also applied to dynamical evolution of a planetesimal disk. Finally, we studied rotation rates of moonlets and particles in planetary rings, which would influence the accretional evolution of these bodies. We describe our key accomplishments during the past three years in more detail in the following.

Ohtsuki, Keiji↗

Ring Current Dynamic in the Presence of EMIC Waves

The effect of EMIC waves, generated by a positive ion temperature anisotropy on Earth s RC ions is one of the best known examples of wave-particle interaction in the magnetosphere and the most controversial mechanism of RC losses. Under certain conditions, relativistic electrons with energy 21 MeV can be removed from the outer RB by EMIC wave scattering during a magnetic storm much faster than by any other loss mechanisms. That is why the calculation of EMIC waves is a very critical part of the NASA LWS program. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch- angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) LHW are included as an additional factor that contributes to saturation process of EMIC waves; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗

Non-Additive Voltametric Currents From a Mixture of Two, Three and Four Redox-Active Compounds and Electroanalytical Implications

We have published recently the effect of dissimilar diffusion coefficients on the size of the voltammetric waves from a mixture of two redox-active compounds. Similarly, at the potential range where three redox-active species, decamethylferrocene (dMeFc), ferrocene (Fc) and N-methylphenothiazine (MePTZ), are oxidized simultaneously with rates controlled by linear diffusion, electrogenerated radicals diffusing outwards from the electrode react with the original species diffusing towards the electrode from the bulk; thus, Fc(+) reacts with dMeFc producing Fc and dMeFc(+), while MePTZ(+) reacts both with dMeFc producing MePTZ and dMeFc(+), and with Fc producing MePTZ and Fc(+). These reactions replace dMeFc with Fc at the second plateau, and both dMeFc and Fc with MePTZ at the third plateau. Since the diffusion coefficients of the three species are not equal, the mass-transfer limited currents of the second and the third oxidation wave plateaus change by approx. 10%. Numerical simulations of the experimental voltamograms support this mechanism. Similar results were also obtained for a mixture of four redoxactive compounds. The implications of this non-additive nature of currents on: (a) the use of internal voltammetric standards for quantitative analysis of a mixture of redox-active compounds; and, (b) the half wave potentials (E1/2) of the 2nd, 3rd and 4th waves for qualitative analysis, will be discussed.

Dass, Amala↗

Ring Current-Electromagnetic Ion Cyclotron Waves Coupling

The effect of Electromagnetic Ion Cyclotron (EMIC) waves, generated by ion temperature anisotropy in Earth s ring current (RC), is the best known example of wave- particle interaction in the magnetosphere. Also, there is much controversy over the importance of EMIC waves on RC depletion. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt (RB) by EMIC wave scattering during a magnetic storm. That is why the calculation of EMIC waves must be a very critical part of the space weather studies. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch-angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) the EMIC waves saturation due to their modulation instability and LHW generation are included as an additional factor that contributes to this process; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗

Inter-species molecular attraction effect in the development of a two-species mixing layer

Simulations of fluid mixing occurring at high-pressures mandate the use of real-fluid equations of state to capture the effect of boh intermolecular repulsive forces and attractive forces which are neglected in perfect gases. Dedicated mixing rules account for these forces by providing a generic mapping of the mixture state-space onto a pure fluid domain. However, the specific inter-species molecular forces depend on the particular characteristics of the molecular pairs under consideration, and typically require the knowledge of binary interaction coefficients k_{αβ}. The impact of k_{αβ} on a flow field is as of now unclear; further, k_{αβ} is unknown for many molecular pairs and when unknown, set to be null. This study addresses the impact of k_{αβ} on the temporal evolution of a real-fluid mixing layer and uses Direct Numerical Simulation (DNS) to explore this impact. The results show differences between augmented attraction (k_{αβ}<0) and diminished attraction (k_{αβ}>0), affecting both fluid dynamics and diffusion processes. Specifically, in a binary mixing layer, it is observed that augmented attraction leads to delayed transition and mixing layer growth, manifested in the momentum layer thickness. The fractal dimension, which is a manifestation of interface corrugation, is calculated. These results provide an hitherto undocumented mechanism affecting the development of mixing layers.

Banuti, Daniel↗

Solutions of the Navier-Stokes equations for vortex breakdown

A numerical investigation of vortex breakdown has been undertaken in an attempt to understand its properties, and the mechanisms responsible for it. Solutions of the full steady axisymmetric Navier-Stokes equations for breakdown in an unconfined viscous vortex have been obtained for core Reynolds numbers up to 200, for a two-parameter family of assumed upstream velocity distributions. Diffusion and convection of vorticity away from the vortex core, and the strong coupling between the circumferential and axial velocity fields in highly-swirling flows, are shown to lead to stagnation and reversal of the axial flow near the axis. The various theories of vortex breakdown are considered in light of the present numerical solutions.

Grabowski, W. J.↗

A new methodology to determine kinetic parameters for one- and two-step chemical models

In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.

Mantel, T.↗