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 1,387 records · Page 77

Isopycnal eddy mixing across the Kuroshio Extension: Stable versus unstable states in an eddying model

The Kuroshio Extension (KE) jet transitions between stable and unstable states on interannual time scales. Cross-jet eddy mixing in the two states is contrasted in the KE region (28°-240°N; 125°-165°E), using a global eddying 0.1° configuration of the Parallel Ocean Program with online numerical particles. The 4 year period chosen (June 1994 to May 1998) covers a full cycle of the stable state, unstable state and the transition period. Large values of cross-jet eddy diffusivities within the KE jet are concentrated in the upper 1000 m. In the upper ocean, elevated cross-jet mixing within the KE jet is mainly concentrated in the downstream part of the KE jet, where the jet is weak but eddy activity is strong. The simulated time-mean KE jet is more intense and extends further east in the stable state than in the unstable state. Consequently, strong cross-jet mixing within the KE jet is located west of 150°E during June 1996 to May 1997 (a typical unstable state), but east of 150°E during June 1995 to May 1996 (a typical stable state). However, average mixing within the KE jet is indistinguishable in the typical stable and unstable states. In the deep ocean, mixing is strongly influenced by topography, and thus their horizontal structures have less inter-annual variability than in the upper ocean. One caveat is that results here cover one representative cycle of the two states. To obtain the climate mean mixing structures for the stable or unstable state, one would need numerical output covering a period much longer than 4 years.

58 GEOSCIENCES↗

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗

A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 1: The 1-D case

In this first part of two papers, we extend the C-method developed in for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities in one space dimension. For gas dynamics, the C-method couples the Euler equations to a scalar reaction-diffusion equation, whose solution C serves as a space-time smooth artificial viscosity indicator. The purpose of this paper is the development of a high-order numerical algorithm for shock-wall collision and bounce-back. Specifically, we generalize the original C-method by adding a new collision indicator, which naturally activates during shock-wall collision. Additionally, we implement a new high-frequency wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplified WENO-based discretization scheme. We show that our scheme improves the order of accuracy of our WENO algorithm, handles extremely strong discontinuities (ranging up to nine orders of magnitude), allows for shock collision and bounce back, and removes high frequency noise. The causes of the well-known “wall heating” phenomenon are discussed, and we demonstrate that this particular pathology can be effectively treated in the framework of the C-method. This method is generalized to two space dimensions in the second part of this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Realizability-preserving DG-IMEX method for the two-moment model of fermion transport

Building on the framework of Zhang & Shu [1], [2], we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function f . The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky & Bludman [3] and Banach & Larecki [4]. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Pauli's exclusion principle, which demands a bounded distribution function (i.e., f ϵ [0,1]). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. The IMEX scheme is formally only first-order accurate, but works well in the diffusion limit, and — without interactions with the background — reduces to the optimal second-order strong stability-preserving explicit Runge-Kutta scheme of Shu & Osher [5]. Numerical results demonstrate the realizability-preserving properties of the scheme. As a result, we also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.

97 MATHEMATICS AND COMPUTING↗

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermal transport exceeding bulk heat conduction due to nonthermal micro/nanoscale phonon populations

While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.

42 ENGINEERING↗

Fluctuating hydrodynamics of electrolytes at electroneutral scales

At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral, and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit, and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross-diffusion and non-ideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species, and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preliminary Reporting of Thermophysical Property Measurements for the Ghareb Formation

Accurate knowledge of thermophysical properties of rock is vital to develop meaningful models of high level nuclear waste emplacement scenarios. The Israel Atomic Energy Commission is considering storing high level nuclear waste in the Ghareb formation, a porous kerogen bearing chalk. Sandia is supporting this effort with an evolving lab- based geomechanics testing program. We have completed measurements of thermal properties up to 275C and room temperature hydrostatic compaction measurements. We report thermal conductivity, thermal diffusivity, specific heat, and mass loss from our thermal measurements, and we report bulk moduli and porosity loss from our compaction measurements. These values are crucial for the numerical models to simulate heat transfer and formation compressibility around a heat generating repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ice Caps and Ice Belts: The Effects of Obliquity on Ice−Albedo Feedback

Planetary obliquity determines the meridional distribution of the annual mean insolation. For obliquity exceeding 55°, the weakest insolation occurs at the equator. Stable partial snow and ice cover on such a planet would be in the form of a belt about the equator rather than polar caps. An analytical model of planetary climate is used to investigate the stability of ice caps and ice belts over the widest possible range of parameters. The model is a non-dimensional diffusive Energy Balance Model, representing insolation, heat transport, and ice−albedo feedback on a spherical planet. A complete analytical solution for any obliquity is given and validated against numerical solutions of a seasonal model in the “deep-water” regime of weak seasonal ice line migration. Multiple equilibria and unstable transitions between climate states (ice-free, Snowball, or ice cap/belt) are found over wide swaths of parameter space, including a “Large Ice-Belt Instability” and “Small Ice-Belt Instability” at high obliquity. The Snowball catastrophe is avoided at weak radiative forcing in two different scenarios: weak albedo feedback and inefficient heat transport (favoring stable partial ice cover), or efficient transport at high obliquity (favoring ice-free conditions). From speculative assumptions about distributions of planetary parameters, three-fourths to four-fifths of all planets with stable partial ice cover should be in the form of Earth-like polar caps.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗