SEARCH · Engineering Papers
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.
The results of a stagnation-line analysis of the radiative heating of a phenolic-nylon ablator are reported. The analysis includes flowfield coupling with the ablator surface, equilibrium chemistry, multicomponent or binary diffusion and a coupled line and continuum radiation calculations. Users manual and operating instructions for the computer programs listed in the report are included.
In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.
Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.
The specular reflection process is shown to be a key aspect of radiation transfer by plant canopies. Polarization measurements are demonstrated as the tool for determining the specular and diffuse portions of the canopy radiance. The magnitude of the specular fraction of the reflectance is significant compared to the magnitude of the diffuse fraction. Therefore, it is necessary to consider specularly reflected light in developing and evaluating light-canopy interaction models for wheat canopies. Models which assume leaves are diffuse reflectors correctly predict only the diffuse fraction of the canopy reflectance factor. The specular reflectance model, when coupled with a diffuse leaf model, would predict both the specular and diffuse portions of the reflectance factor. The specular model predicts and the data analysis confirms that the single variable, angle of incidence of specularly reflected sunlight on the leaf, explains much of variation in the polarization data as a function of view-illumination directions.
We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.
It has previously been demonstrated that simple coupled reaction-diffusion models can approximate the aging behavior of PMR-15 resin subjected to different oxidative environments. Based on empirically observed phenomena, a model coupling chemical reactions, both thermal and oxidative, with diffusion of oxygen into the material bulk should allow simulation of the aging process. Through preliminary modeling techniques such as this it has become apparent that accurate analytical models cannot be created until the phenomena which cause the aging of these materials are quantified. An experimental program is currently underway to quantify all of the reaction/diffusion related mechanisms involved. The following contains a summary of the experimental data which has been collected through thermogravimetric analyses of neat PMR-15 resin, along with analytical predictions from models based on the empirical data. Thermogravimetric analyses were carried out in a number of different environments - nitrogen, air and oxygen. The nitrogen provides data for the purely thermal degradation mechanisms while those in air provide data for the coupled oxidative-thermal process. The intent here is to effectively subtract the nitrogen atmosphere data (assumed to represent only thermal reactions) from the air and oxygen atmosphere data to back-figure the purely oxidative reactions. Once purely oxidative (concentration dependent) reactions have been quantified it should then be possible to quantify the diffusion of oxygen into the material bulk.
In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.
This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.
Translation and rotation are the two most fundamental forms of diffusion, yet their coupling mechanism is not clear, especially under confinement. Here, we provided evidence of the coupling between rotation and translation using a substituted benzene molecule as an example. A counterintuitive behavior was observed where the movement of the smaller molecule with an asymmetric shape was unexpectedly slower than the larger one with a symmetric shape in confined channels of zeolite. We showed that this diffusion behavior was caused by the presence of the specific and selective interaction of the asymmetric guest with the pores, which increased the local restricted residence time, thus inhibiting the translation under confinement, as further confirmed by dynamic breakthrough curves, uptake measurements, quasi-elastic neutron scattering, and 2 H solid-state NMR techniques. Our work correlated asymmetric rotation and diffusion under a confined environment, which enriched our understanding of the coupling between rotation and translation and could shed light on a fundamental understanding of the diffusion process.
A time dependent model of F region structure decay by 'classical' cross field diffusion and electrical coupling along magnetic field lines to the E region is examined. The temporal behavior of the ion concentration fluctuations is determined by the electric field in the coupled system as well as by the initial perturbation spectra and the E region recombination rate. The formation of image structure in the E region ion concentration affects the lifetime of F layer structure in a scale size dependent way. Once an image is formed, the image amplitude and the driving F region structure amplitude decay at the same rate. At large scale sizes of at least lambda(2pi/k), this rate is proportional to k2 and the ratio of the temperatures in each region. At small scale sizes it depends on the E region recombination rate and the temperatures of the two regions but is only very weakly dependent on k. The background E region concentration determines the wave number beyond which the structure amplitude decay rate is almost independent of its scale size.
The Lead-Zinc binary equilibrium system is currently being investigated. Ground based studies of this system were performed to examine the possibility of obtaining a couple which, after diffusion, could be examined continuously along the diffusion axis by quantitative metallography to determine the extent of diffusion. The specimens were analyzed by X-ray fluorescence in the scanning electron microscope to provide exact information on the chemical composition gradient. Two diffusion experiments were run simultaneously in the multipurpose furnace, each in its own isothermal cavity. Two flight samples, two flight backup samples, and two flight space samples were generated.
The effects of diffusion on the coupling between counterpropagating modes in a two-level medium are studied in the rate equation limit. It is shown that diffusion can greatly decrease saturation gratings in low-pressure gases. Among many other effects which rely on spatial hole burning, this can alter the stability of bidirectional solutions of the third-order Lamb equation model of a ring laser.
A new approach based on the coupling of an Eulerian and a Lagrangian solver, aimed at reducing the numerical diffusion errors of standard Eulerian time-marching finite-volume solvers, is presented. The approach is applied to the computation of the secondary flow in two bent pipes and the flow around a 3D wing. Using convective point markers the Lagrangian approach provides a correction of the basic Eulerian solution. The Eulerian flow in turn integrates in time the Lagrangian state-vector. A comparison of coarse and fine grid Eulerian solutions makes it possible to identify numerical diffusion. It is shown that the Eulerian/Lagrangian approach is an effective method for reducing numerical diffusion errors.
Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.
The electron precipitation bursts have been shown to be a major sink for the radiation belt relativistic electrons. As underlying mechanism of such bursts, we propose particle scattering into the loss cone due to nonlinear resonance interaction between electrons and chorus. Stochastic heating due to the coupling leads to diffusion in pitch angle, and the rate of diffusion would be sufficient to account for the emptying of the Earth's radiation belt over the time of the main phase of geomagnetic storms. The results obtained in the present paper account for a strong energy dependence in the electron precipitation event and the correlation between the energization and loss processes on macroscopic timescales, which is primarily attributed to the cooperative effects of the coupling. This mechanism of chorus scattering should produce pitch angle distributions that are energy-dependent and butterfly-shaped. The calculated timescales and the total energy input to the atmosphere from precipitating relativistic electrons are in reasonable agreement with experimental data.
Discussion of the temperature-density phase anomaly on the basis of a quasi-three-dimensional model in which the thermosphere dynamics associated with wind circulation is considered in a self-consistent form. Included in this analysis are the first three harmonics, which involve nonlinear coupling between diurnal and semidiurnal tides. It is shown that the phase anomaly with exospheric temperature peaks near 1600 LT and mass density peaks between 1400 and 1445 LT can be reproduced in a self-consistent theory without invoking ad hoc assumptions and boundary conditions that would mask the physical processes to be explored. A number of factors and processes are found to contribute to the phase anomaly, including the semidiurnal and particularly the terdiurnal components, heat advection, diffusion, and energy coupling with the lower atmosphere.
A general model of particle acceleration by plasma waves coupled with adiabatic radial diffusion in a planetary magnetic field is developed. The model assumes that a spectrum of lower hybird waves is present to resonantly accelerate ions transverse to the magnetic field. The steady state Green's function for the combined radial diffusion and wave acceleration equation is found in terms of a series expansion. The results provide a rigorous demonstration of how a quasi-Maxwellian distribution function is formed in the absence of particle collisons and elucidate the nature of turbulent heating of magnetospheric plasmas. The solution is applied to the magnetosphere of Neptune for which a number of examples are given illustrating how the spectrum of pickup N(+) ions from Triton evolves.
The results for a stagnation-line analysis of the radiative heating of a phenolic-nylon ablator are presented. The analysis includes flow field coupling with the ablator surface, equilibrium chemistry, a step-function diffusion model and a coupled line and continuum radiation calculation. This report serves as the documentation, i e. users manual and operating instructions for the computer programs listed in the report.