Search NASA⌕ Search

SEARCH · Search NASA

Results for “solution transfer”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗
NASA NTRS1999

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗
DOE OSTI2022

Conservative and Accurate Solution Transfer Between High-Order and Low-Order Refined Finite Element Spaces

Here, in this paper, we introduce general transfer operators between high-order and low-order refined finite element spaces that can be used to couple high-order and low-order simulations. Under natural restrictions on the low-order refined space we prove that both the high-to-low-order and low-to-high-order linear mappings are conservative, constant preserving, and high-order accurate. While the proof holds for affine geometries, numerical experiments indicate that the results hold for more general curved meshes. We present several numerical results confirming our analysis and demonstrate the utility of the new mappings in the context of adaptive mesh refinement and conservative multidiscretization coupling.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Radiative transfer solutions for cylindrical coordinates with emitting, absorbing, and anisotropic scattering medium

The radiative transfer equation is solved by the S-N discrete ordinates method in one-dimensional radial and two-dimensional r-z coordinates systems. The walls of the enclosures are diffuse, and the participating medium absorbs, emits and anisotropically scatters the radiative energy. Diffuse wall incidence, isothermal medium emission and collimated incidence problems are considered. Effects of the scattering phase functions on average incident radiation and net radiative heat fluxes are studied. In addition, the effects of scattering albedo, optical thickness and the wall emissivity are briefly discussed.

Jendoubi, Slah↗
NASA NTRS1999

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗
NASA NTRS1992

Comparisons of cylindrical and rectangular media radiative transfer solutions using the S-N discrete ordinates method

The application of the S-N discrete ordinates method is compared for r-z and x-y problems with similar boundary conditions. The complexities of the r-z solution come from the additional angular derivative, which makes this into a virtual three-dimensional problem. Effects of spatial differencing schemes, phase functions, and grid spacing are considered and compared for these x-y and r-z problems.

Jendoubi, Slah↗
DOE OSTI2024

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗
DOE OSTI2020

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Inverted Region in Bimolecular Electron Transfer in Solution Enabled by Delocalization

Rate constants for bimolecular electron transfer (ET) increased with driving force, -ΔG°, reached a plateau, and then decreased in an inverted region. This rate data was described well by electron transfer theory subject to a diffusion-controlled limit. These were for ET from radical anions of polydecylthiophene (P3DT) to a series of acceptors in THF solution. When the donor was the smaller anion of quaterthiophene (T4 •- ) the inverted region was much less prominent, and still less so for when the donor was the anion of bithiophene (T2 •- ). Description of the data using ET theory identifies smaller electronic couplings for the highly-delocalized P3DT anions as enabling the inverted behavior: The presence of a Marcus inverted region is a consequence of delocalized electronic states. The results of this study further imply that electronic couplings smaller than usually found for molecules in contact could boost efficiency of energy storage by electron transfer, and identifies size-mismatch as an important concept in control of electronic couplings.

14 SOLAR ENERGY↗
DOE OSTI2021

Development of Approaches to Model Excited State Charge and Energy Transfer in Solution

The development of next generation energy conversion and catalytic systems requires fundamental understanding of the interplay between photo-excitation and the resulting proton and electron transfer processes that occur in solution. To address these challenges requires accurate and efficient methods to compute ground and excited states, as well as the ability to treat the dynamics of energy transfer in the presence of solvent fluctuations. Our team is developing accurate and efficient theoretical models for solution phase reactions. These developments provide an improved understanding of photo-initiated excitations, electron transfer, and proton coupled electron transfer processes. Our research team brings together expertise in electron and nuclear dynamics (classical and quantum), electronic structure, and analysis of solvation networks. We are developing and validating techniques for modeling condensed phase charge and energy transfer and to simulate the methods to monitor them using optical and vibrational spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1992

A simple analytical function for bidirectional reflectance

A simple bidirectional reflectance model based on physical scattering laws is developed to calculate the bidirectional reflectance of a wide variety of surfaces. Chandrasekhar's radiative transfer solution, obtained for the anisotropically scattering semi-infinite medium, is used to compute the multiple-scattered radiances. However, since the radiative transfer solution does not take architectural effects into consideration, we adopt Hapke's (1986) approach and add an empirical term to explain the hot-spot phenomenon and use Cox and Munk's (1954) formulation to take into account the specular reflection. The physical parameters of the model are retrieved from bidirectional reflectance measurements. For each surface, only one set of model parameters is needed for application to all illumination and viewing geometries. The validity of the model is established by comparing the computed and measured reflectances for sets of viewing and illumination angles that were not included in the inversion algorithm. Good agreement is shown between the model-computed and the observed reflectances for dense prairie vegetation canopies, a sparse desert scrub community, a plowed agricultural field, and an alkali flat.

Ahmad, Suraiya P.↗
NASA NTRS1994

Some Problems and Solutions in Transferring Ecosystem Simulation Codes to Supercomputers

Many computer codes for the simulation of ecological systems have been developed in the last twenty-five years. This development took place initially on main-frame computers, then mini-computers, and more recently, on micro-computers and workstations. Recent recognition of ecosystem science as a High Performance Computing and Communications Program Grand Challenge area emphasizes supercomputers (both parallel and distributed systems) as the next set of tools for ecological simulation. Transferring ecosystem simulation codes to such systems is not a matter of simply compiling and executing existing code on the supercomputer since there are significant differences in the system architectures of sequential, scalar computers and parallel and/or vector supercomputers. To more appropriately match the application to the architecture (necessary to achieve reasonable performance), the parallelism (if it exists) of the original application must be exploited. We discuss our work in transferring a general grassland simulation model (developed on a VAX in the FORTRAN computer programming language) to a Cray Y-MP. We show the Cray shared-memory vector-architecture, and discuss our rationale for selecting the Cray. We describe porting the model to the Cray and executing and verifying a baseline version, and we discuss the changes we made to exploit the parallelism in the application and to improve code execution. As a result, the Cray executed the model 30 times faster than the VAX 11/785 and 10 times faster than a Sun 4 workstation. We achieved an additional speed-up of approximately 30 percent over the original Cray run by using the compiler's vectorizing capabilities and the machine's ability to put subroutines and functions "in-line" in the code. With the modifications, the code still runs at only about 5% of the Cray's peak speed because it makes ineffective use of the vector processing capabilities of the Cray. We conclude with a discussion and future plans.

Skiles, J. W.↗