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,135 records · Page 63

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

Time-Resolved Evolution of the 3D Nanoporous Structure of Sintered Ag by X-Ray Nanotomography: Role of the Interface with a Copper Substrate

Here, the evolution of the nanoporous structure of cylindrical sintered silver samples during high temperature aging, ranging from 200 to 350 °C for 350 min, is studied through in situ computed X-ray tomography. Investigations are carried out for two types of specimens: pure sintered silver and specimens containing a silver-copper interface. It is shown that the overall pore evolution is driven by the evolution of very few large ones. The smaller pores, although being more numerous, do not really evolve before being absorbed by the few bigger ones. In pure silver, pore evolution is driven by diffusion (Ostwald ripening) but the presence of an interface promotes faster growth kinetics until the aging time reaches a threshold value, after which a deviation from Ostwald ripening occurs. The transition is a function of the aging temperature. This behavior is associated with the competition between elastic relaxation and surface energy minimization.

42 ENGINEERING↗

Stochastic and coarse-grained two-dimensional modeling of directional particle movement

In this work, we study the evolution of interacting groups of pedestrians in two-dimensional geometries. We introduce a microscopic stochastic model that includes floor fields modeling the global flow of individual groups as well as local interaction rules. From this microscopic model we derive an analytically-tractable system of conservation laws that governs the evolution of the macroscopic pedestrian densities. Numerical simulations show good agreement between the system of conservation laws and the microscopic model, though the latter is slightly more diffusive. We conclude by deriving second-order corrections to the system of conservation laws.

97 MATHEMATICS AND COMPUTING↗

Characterization of Grain‐Size Distribution, Thermal Conductivity, and Gas Diffusivity in Variably Saturated Binary Sand Mixtures

Core Ideas Parameterized models are relevant to applications where different porous media mixtures are used. Models include a grain‐size distribution function to describe bimodal behavior for binary mixtures. Improved model describes observed thermal conductivity–saturation relations for binary mixtures. Combined Buckingham–Penman model used to describe observed gas diffusivity–air content relations. This work is relevant for proper simulation of mixed porous media. Characterization of differently textured porous materials, as well as different volumetric porous media mixtures, in relation to mass and heat transport is vital for many engineering and research applications. Functional relations describing physical (e.g., grain‐size distribution, total porosity), thermal, and gas diffusion properties of porous media and mixtures are necessary to optimize the design of porous systems that involve heat and gas transport processes. However, only a limited number of studies provide characterization of soil physical, thermal, and gas diffusion properties and the functional relationships of these properties under varying soil water contents, especially for soil mixtures, complicating optimization efforts. To better understand how mixing controls the physical, thermal, and gas diffusion properties of porous media, a set of laboratory experiments was performed using five volumetric mixtures of coarse‐ and fine‐grained sand particles. For each mixture, the grain‐size distribution (GSD), thermal conductivity, and gas diffusivity were obtained and parameterized using existing and suggested parametric models. Results show that the extended, two‐region Rosin–Rammler particle‐size distribution model proposed in this study could successfully describe the bimodal behavior of the GSD of binary mixtures. Further, the modified Côté and Konrad thermal conductivity model adequately described the thermal conductivity–water saturation relations observed in different mixtures. The proposed simple soil‐gas diffusivity descriptive model parameterized the upper limit, average, and lower limit behavior in gas diffusivity–air content relations in apparently texture‐invariant gas diffusivity data. Results further show a close analogy between gas diffusivity and thermal conductivity and their variation with saturation across different binary mixtures. Overall, the results of the study provide useful numerical insight into the physical, thermal, and gas transport characteristics of binary mixtures, with wide implications for future engineering and research applications that involve multicomponent porous systems.

Deepagoda, T. K. K. Chamindu↗

On the calculation of turbulent heat transport downstream from an abrupt pipe expansion

A numerical study is reported of flow and heat transfer in the separated flow region created by an abrupt pipe expansion. Computations employed an adaptation of the TEACH-2E computer program with the standard model of turbulence. Emphasis is given to the simulation, from both a physical and numerical viewpoint, of the region in the immediate vicinity of the wall where turbulent transport gives way to molecular conduction and diffusion. Wall resistance laws or wall functions used to bridge this near-wall region are based on the idea that, beyond the viscous sublayer, the turbulent length scale is universal, increasing linearly with distance from the wall. Predictions of expermental data for a diameter ratio of 0.54 show generally encouraging agreement with experiment. At a diameter of 0.43 different trends are discernible between measurement and calculation though this appears to be due to effects unconnected with the wall region studied.

Chieng, C. C.↗

A numerical model study of long-term planetary wave predictability

A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.

Wobus, R. L.↗

Effects of unsteady free stream velocity and free stream turbulence on stagnation point heat transfer

The combined effects of transient free stream velocity and free stream turbulence on heat transfer at a stagnation point over a cylinder situated in a crossflow are studied. An eddy diffusivity model was formulated and the governing momentum and energy equations are integrated by means of the steepest descent method. The numerical results for the wall shear stress and heat transfer rate are correlated by a turbulence parameter. The wall friction and heat transfer rate increase with increasing free stream turbulence intensity.

Gorla, R. S. R.↗

Combined influence of unsteady free stream velocity and free stream turbulence on stagnation point heat transfer

An analysis is presented to study the combined effects of transient free stream velocity and free stream turbulence on heat transfer at a stagnation point over a cylinder situated in a crossflow. An expression for the eddy diffusivity has been formulated and the governing momentum and energy equations are integrated by means of the steepest descent method. The numerical results for the heat transfer rate are correlated by a turbulence parameter. It has been found that the stagnation point heat transfer rate increases with increasing free stream turbulence intensity.

Gorla, R. S. R.↗

Infrared radiation parameterizations in numerical climate models

This study presents various approaches to parameterizing the broadband transmission functions for utilization in numerical climate models. One-parameter scaling is applied to approximate a nonhomogeneous path with an equivalent homogeneous path, and the diffuse transmittances are either interpolated from precomputed tables or fit by analytical functions. Two-parameter scaling is applied to parameterizing the carbon dioxide and ozone transmission functions in both the lower and middle atmosphere. Parameterizations are given for the nitrous oxide and methane diffuse transmission functions.

Chou, Ming-Dah↗

Navier-Stokes simulation of flow through a highly contoured subsonic diffuser

The flow through a highly offset subsonic diffuser with cross-sectional profiles that varied from rectangular at the duct entrance to circular at the engine face was numerically simulated. A multizonal approach combined with a two-grid topology was used to represent both the internal and external flowfields, and an implicit, approximately-factored, partially flux-split finite-difference algorithm was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computed static pressures along the inlet wall and total pressures on the engine face were compared with experimental data. In addition, the overall flowfield within the duct was examined in detail. Good agreement is shown between experiment and computations, with the limiting factor being the lack of a reliable turbulence model for internal flow problems.

Chyu, W. J.↗

Mixing and mixing enhancement in supersonic reacting flowfields

Propulsion systems planned for use late in this century and beyond will require appropriate physical models for describing supersonic combustion and numerical techniques for solving the model governing equations. A computer program to study these flows is reported which considers the multicomponent diffusion and convection of important chemical species, the finite-rate reaction of these species, and the resulting interaction of the field mechanics and the chemistry. The application of the program to a spatially developing and reacting mixing layer, which serves an an excellent physical model for the mixing and reaction processes that take place in a scramjet combustor, is reported. Several techniques to enhance the fuel-air mixing and growth of that layer and improve its overall combustion efficiency are considered.

Drummond, J. P.↗

Cosmic Rays and Their Radiative Processes in Numerical Cosmology

A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive, acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of 128 (exp 3) cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.

Ryu, Dongsu↗

Effects of Heat Release on a Two-Dimensional Reacting Shear Layer

The diffusion controlled binary reaction between initially segregated reactants in a two-dimensional low Mach number mixing layers is studied via numerical simulation. The stoichiometric ratio of the reactants is chosen to be much larger than one, as is typical of hydrocarbon flames in air. This results is a flame that is offset from the main vortical region of the mixing layer. In agreement with experimental observations, the flame remains surprisingly uncontorted during the flow evolution and is not entrained into the mixing layer. The effect of the heat release of the flame on the evolution of the mixing layer is thus similar to the effect of a difference in free-stream density between the two sides of the layer. The resulting baroclinic torque inhibits the familiar rolup and pairing of mixing layer vortices common in constant density flows. This also contributes to the layers inability to entrain the flame. The increase in viscosity caused by the heating of the flame reduces the effective Reynolds number of the flow. But, contrary to what has commonly been suggested, this is not the major reason for the inhibition of the usual large-scale mixing layer structures.

Moser, Robert↗

Fast Three-Dimensional Method of Modeling Atomic Oxygen Undercutting of Protected Polymers

A method is presented to model atomic oxygen erosion of protected polymers in low Earth orbit (LEO). Undercutting of protected polymers by atomic oxygen occurs in LEO due to the presence of scratch, crack or pin-window defects in the protective coatings. As a means of providing a better understanding of undercutting processes, a fast method of modeling atomic-oxygen undercutting of protected polymers has been developed. Current simulation methods often rely on computationally expensive ray-tracing procedures to track the surface-to-surface movement of individual "atoms." The method introduced in this paper replaces slow individual particle approaches by substituting a model that utilizes both a geometric configuration-factor technique, which governs the diffuse transport of atoms between surfaces, and an efficient telescoping series algorithm, which rapidly integrates the cumulative effects stemming from the numerous atomic oxygen events occurring at the surfaces of an undercut cavity. This new method facilitates the systematic study of three-dimensional undercutting by allowing rapid simulations to be made over a wide range of erosion parameters.

Snyder, Aaron↗

APC: A New Code for Atmospheric Polarization Computations

A new polarized radiative transfer code Atmospheric Polarization Computations (APC) is described. The code is based on separation of the diffuse light field into anisotropic and smooth (regular) parts. The anisotropic part is computed analytically. The smooth regular part is computed numerically using the discrete ordinates method. Vertical stratification of the atmosphere, common types of bidirectional surface reflection and scattering by spherical particles or spheroids are included. A particular consideration is given to computation of the bidirectional polarization distribution function (BPDF) of the waved ocean surface.

surface reflection↗

Distal Ejecta from Lunar Impacts: Extensive Regions of Rocky Deposits

Lunar Reconnaissance Orbiter (LRO) Diviner Radiometer, Mini-RF, and LRO Camera data were used to identify and characterize rocky lunar deposits that appear well separated from any potential source crater. Two regions are described: 1) A approximate 18,000 sq km area with elevated rock abundance and extensive melt ponds and veneers near the antipode of Tycho crater (167.5 deg E, 42.5 deg N). This region has been identified previously, using radar and aging data. 2) A much larger and more diffuse region, covering approximately 730,000 sq km, centered near 310 deg E, 35 deg S, containing elevated rock abundance and numerous granular flow deposits on crater walls. The rock distributions in both regions favor certain slope azimuths over others, indicating a directional component to the formation of these deposits. The spatial distribution of rocks is consistent with the arrival of ejecta from the west and northwest at low angles (approximately 10-30 deg) above the horizon in both regions. The derived age and slope orientations of the deposits indicate that the deposits likely originated as ejecta from the Tycho impact event. Despite their similar origin, the deposits in the two regions show significant differences in the datasets. The Tycho crater antipode deposit covers a smaller area, but the deposits are pervasive and appear to be dominated by impact melts. By contrast, the nearside deposits cover a much larger area and numerous granular flows were triggered. However, the features in this region are less prominent with no evidence for the presence of impact melts. The two regions appear to be surface expressions of a distant impact event that can modify surfaces across wide regions, resulting in a variety of surface morphologies. The Tycho impact event may only be the most recent manifestation of these processes, which likely have played a role in the development of the regolith throughout lunar history

Bandfield, Joshua L.↗

Spinbox: tools for many-body quantum systems in a Monte Carlo context

Spinbox is a piece of software that facilitates quantum mechanical calculations relevant to Monte Carlo simulation of atomic nuclei. At the front lines of research on the nuclear many-body problem are a large number of supercomputer-scale simulation codes. These codes produce valuable results but can be hard to understand, especially for those without intimate knowledge of the relevant theoretical methods. Thus, tools that fill pedagogical roles are extremely valuable. Spinbox makes it easy for one to replicate and analyze the computational processes relevant to a Quantum Monte Carlo (QMC) simulation that may be difficult to understand/debug/analyze due to the scale of the corresponding simulation software. Spinbox is written in Python using other state-of-the-art Python modules for numerical calculations. While a number of Python libraries exist that are suited to general quantum many-body calculations, the motivation of Spinbox is quite particular. In Diffusion Monte Carlo methods (DMC, GFMC, AFDMC), the central calculation is the imaginary-time propagation of individual samples of the many-body wavefunction. Although quantum wavefunctions generally must be described by a probability distribution over a basis, DMC imbues particles (within one sample) with classical spatial coordinates. This method is unusual, so other Python packages are typically not set up to do this easily. Furthermore, the software has built-in options for nuclear systems assuming isospin symmetry, which can be set up with other libraries but is a nontrivial process to do so. Features: - numerical representation of samples of the many-body wavefunctions, including tensor-product states (used in AFDMC) - numerical representation of many-body operators, including tensor-product operators: general, spin, imaginary-time propagation, etc. - the correct associated arithmetic and algebra, implemented as class methods - classes for representing realistic nuclear two- and three-body Hamiltonians (e.g. Argonne V18, Illinois NNN) - large-scale parallel integration over random variables, crucial for the AFDMC method My goal is to make this package open source so that anyone may use it and contribute to it, particularly other researchers doing AFDMC calculations

Fox, Jordan↗

Numerical results for the thermal scattering functions

A recent formulation in radiative transfer defined the thermal scattering functions that characterize radiative transfer from a general plane-parallel finite medium driven solely by an internal distribution of thermal sources. Exiting diffuse intensities are expressed as space convolutions of the thermal scattering functions with any thermal source distribution. A parametric study is presented to obtain the basic structure of these scattering functions. The independent variables of these azimuthally independent functions are the direction cosine and source location, while the parameters are the single-scattering albedo, total optical depth, and the asymmetry factor in the Henyey-Greenstein phase function. The basic functional trends are discussed by using various parametric plots, and selected results are given to allow numerical checks. The computational method is invariant imbedding.

Cogley, A. C.↗