Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion rate”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗
DOE OSTI2024

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

Herein this paper focuses on developing a reduction-based algebraic multigrid (AMG) method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based AMG approach, $\ell \text{AIR}$ (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection-dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion-dominated problems in two or three dimensions. Motivated by the success of $\ell \text{AIR}$ in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion-dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy-minimization AMG methods with the local approximation of ideal operators used in $\ell \text{AIR}$. The resulting constrained $\ell \text{AIR}$ algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that was previously difficult for reduction-based methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1946

Performance of the Modified V-1710-93 Engine-Stage Supercharger with a Constant-Area Vaneless Diffuser

As part of an investigation to increase the power output of the V-1710-93 engine at altitude, the engine-stage supercharger was combined with a constant-area vaneless diffuser designed to improve the performance of the engine-stage supercharger at the rated engine operating point. The performance of the modified supercharger was investigated in a variable-component supercharger test rig and compared with that of the standard supercharger with an 8-vaned diffuser. A separate evaluation of the component efficiencies and a study of the flow characteristics of the modified supercharger was made possible by internal diffuser instrumentation. At the volume flow required by the engine for rated operating conditions, the modified supercharger increased the over-all adiabatic efficiency 0.05 and the over-all pressure coefficient 0.035. Furthermore, the capacity of the engine-stage supercharger was increased by replacing the standard 8-vaned diffuser with the vaneless diffuser. The peak over-all adiabatic efficiency for the modified supercharger, however, was 0.05 to 0.07 lower than that of the standard unit over the range of tip speeds investigated. The improved performance of the modified supercharger at rated engine operating conditions resulted from a shift of the point of peak adiabatic efficiency and pressure coefficient of the standard supercharger to a higher flow. The energy loss through the vaneless diffuser was found to be small. Because of the restricted diffuser diameter, however, diffusion was inadequate, which resulted in a relatively small static-pressure rise through the diffuser, high diffuser-exit velocities, and excessive collector-case losses.

Douglas, John E.↗
DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
NASA NTRS1988

Diffusional transport during the cyclic oxidation of gamma + beta, Ni-Cr-Al(Y, Zr) alloys

The cyclic oxidation behavior of several cast gamma + beta, Ni-Cr-Al(Y, Zr) alloys and one low-pressure plasma spraying gamma + beta, Ni-Co-Cr-Al(Y) alloy was studied. Cyclic oxidation was found to result in a decreasing Al concentration at the oxide-metal interface due to a high rate of Al consumption coupled with oxide scale cracking and spalling. Diffusion paths plotted on the ternary phase diagram showed higher Ni concentrations with increasing cyclic oxidation exposures. The alloy with the highest rate of Al consumption and the highest Al content underwent breakaway oxidation following 500 1-hr cycles at 1200 C.

Nesbitt, J. A.↗
NASA NTRS1977

Prediction of laminar and turbulent boundary layer flow separation in V/STOL engine inlets

The paper provides a theoretical description of the development of the boundary layer on the lip and diffuser surface of a subsonic inlet at arbitrary operating conditions of mass flow rate, freestream velocity and incidence angle. Both laminar separation on the lip and turbulent separation in the diffuser are discussed. The agreement of the theoretical results with model experimental data illustrates the capability of the theory to predict separation. The effects of throat Mach number, inlet size, and surface roughness on boundary-layer development and separation are illustrated.

Chou, D. C.↗
NASA NTRS2022

Biofilm Study Under Simulated Microgravity

The goal of this study was to understand biofilm formation under microgravity (µg), in support of biofilm mitigation efforts in exploration water recovery systems. The technical approach was to conduct a mass transfer and bacterial culture study under both simulated µg and ambient gravity. The aim was to correlate nutrient consumption to gene expression to better understand biofilm formation. A representative species of bacteria that is commonly cultured from the International Space Station (ISS) Water Processor Assembly (WPA) was cultured in a WPA influent water ersatz formulation that is tailored for microbiology studies. A mass transfer rate study was carried out using the ersatz WPA influent water by introducing a water-soluble dye to represent dissolved nutrients and nutrient particles. Imaging of dye diffusion over time allowed for the comparison of mass transport rates under a series of rotation per minute (RPM) speeds for the High Aspect Ratio Vessels (HARVs) on a Rotating Wall Vessel (RWV). This was done to determine the speed that will most accurately simulate the low convective rates experienced under actual µg conditions. Three biological replicates of the Burkholderia contaminans (B. contaminans) microbe were cultured under simulated µg with a rotating (R) control in the horizontal plane at the determined optimal RPM of 15, along with a stationary (S) reference culture. At T=0, and then at T=1,2,3 (in exponential phase) and T=4 (in S phase), the bacterial culture and ersatz were harvested for transcriptomic and nutrient content analysis, respectively. The experimental results illustrated that phosphate is a limiting nutrient in the WPA ersatz formula. Nutrient analysis illustrated that the µg treatment culture took up essential nutrients more rapidly than the R and S control cultures, yet non-essential nutrients remained higher in the µg treatment than in the controls at later timepoints. The rapid uptake and subsequent starvation of phosphate in the culture under µg conditions is further illustrated in the transcriptomic response when compared to that of the R control condition. The subsequent starvation response may serve as one element to explain a moderate enhancement of biofilm formation in the µg treatment. One implication of this work is that biofilm mitigation in the ISS environment could be supported by ensuring a steady flow of water as a vehicle for phosphate within the WPA to avoid complete phosphate consumption, which occurs in times of no flow and leads to undesired biofilm formation.

Aubrie O’Rourke↗
NASA NTRS2017

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗
DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗
DOE OSTI2025

Role of pore dilation in molecular transport through the nuclear pore complex: Insights from polymer scaling theory

The nuclear pore complex (NPC), a channel within the nuclear envelope filled with intrinsically disordered proteins, regulates the transport of macromolecules between the nucleus and the cytoplasm. Recent studies have highlighted the NPC’s ability to adjust its diameter in response to the membrane tension, underscoring the importance of exploring how variations in pore size influence molecular transport through the NPC. In this study, we investigated the relationship between pore size and transport rate and proposed a mathematical model describing this connection. We began by theoretically analyzing how the pore size scales with the characteristic dimensions of the mesh-like structure within the pore. By introducing key assumptions about how the meshwork structure influences molecular diffusion, we derived a mathematical expression for the transport rate based on the size of the pore and the transported molecules. To validate our model, we conducted Brownian dynamics simulations using a coarse-grained representation of the NPC. These simulations, performed across a range of pore sizes, demonstrated strong agreement with our model’s predictions, confirming its accuracy and applicability. Our model is specifically tailored for small-to-medium-sized molecules, approximately 5 nanometers in size, making it relevant to a wide range of transcription factors and signaling molecules. It also extends to molecules with weak and transient interactions with FG-Nups, such as importin-β. By presenting this model formula, our study offers a quantitative framework for analyzing the effects of pore dilation on nucleocytoplasmic transport.

Biochemistry & Molecular Biology↗
NASA NTRS1995

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗
NASA NTRS1977

Photochemistry of the Martian atmosphere

A critical analysis is carried out for models which may be fine tuned to give agreement with observational constraints for O2, CO, and O3 in the Martian lower atmosphere and which are in accord also with upper atmospheric data for O and CO. The models have a number of features in common: dynamic mixing must be exceedingly rapid at heights above 90 km; if heterogeneous processes are not assumed to play a major role in the removal of atmospheric CO and O2, mixing must be rapid at heights between 30 and 40 km. The value implied for the diffusion coefficient in this region depends on the assumptions made regarding the rates for reaction of OH with HO2 to form H2O and the rate of reaction of HO2 with itself to form H2O2. If these reaction rates are taken to have values similar to those used in current models of the earth's stratosphere, the eddy diffusion coefficient at 40 km above the Martian surface should be about 5 by 10 to the seventh power sq cm/sec.

Kong, T. Y.↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tensile behavior of diffusion bonded AA6061 - AA6061 with variation in cooling method

Hot isostatically pressed AA6061 cladding is an important structural component of the high performance, Zr-laminated U-10Mo monolithic fuel system for the application in research and test reactors. In this study, the mechanical behavior of two diffusion bonded aluminum alloy, AA6061, was examined using tensile testing. Solid-to-solid diffusion bonding between two pieces of AA6061 was performed by isothermal annealing at 560 °C for 1.5 h, and diffusion couples were subsequently cooled via three different cooling methods: furnace cooling, air cooling, and water quenching. Dog-bone shaped tensile specimens, with 10 mm in gauge length (with diffusion bonded interface in the middle), and 1.5 × 1.5 mm 2 gauge cross-sections, were fabricated from the diffusion bonded AA6061 by electro-discharge machining. Yield strength (% EL at failure) of furnace cooled, air cooled and water quenched tensile specimens determined was 82 ~ 89 MPa (10 ~ 30%), 112 ~ 116 MPa (10 ~ 14%), and 149 ~ 164 MPa (10 ~ 17%), respectively. This variation in mechanical behavior was examined with cooling-rate dependent, concentrated precipitation of Mg 2 Si at the diffusion bonded interface, with due respect for mechanical properties of the AA6061 alloy that inherently vary as a function of cooling rate from 560 °C. Finite element analysis using ABAQUS was employed to augment experimental findings with the appropriate microstructural constituents and alloy properties. Results suggest that the strength is dominated by matrix/bulk properties of AA6061, while ductility is strongly influenced by the cooling method dependent presence of Mg 2 Si precipitates at the interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Soil porous microstructure control over soil organic matter mobility: A multimethod workflow for understanding chemistry-dependent organic matter binding in soil

Soil organic matter (SOM) has attracted a great deal of interest; particularly for its potential to mitigate human derived CO 2 emissions. Studies have demonstrated that SOM plays a critical role in carbon storage and CO 2 sequestration. However, the sorption properties of SOM, which influence its transport in pore water and stabilization within the soil, remain poorly understood. This study develops a workflow to: (1) examine compound-specific advective and diffusive transport and desorption behaviors, (2) quantify desorption rates through stop-flow and continuous-flow column experiments, and (3) evaluate the impact of soil microporosity on SOM mobility using high-resolution imaging and extractions. Intact core column experiments were conducted on Uncultivated (Natural) and Cultivated soil samples, both were arid soils, collected in Washington State. X-ray computed tomography was employed to measure porosity and pore connectivity, while Fourier-transform ion cyclotron resonance mass spectrometry was used to analyze SOM composition. The findings revealed that cultivation increased total carbon and nitrogen levels due to irrigation and fertilization, enhancing carbon capture potential in arid soils. In contrast, the Natural soil, characterized by higher porosity and connectivity, contained more oxidized carbon. Pore network analysis indicated that soil compaction in the Cultivated soil may lead to longer diffusion pathways, significantly influencing SOM transport and stability.

Hydraulic Properties↗
NASA NTRS1977

Prediction of Laminar and Turbulent Boundary Layer Flow Separation in V/STOL Engine Inlets

A description is presented of the development of the boundary layer on the lip and diffuser surface of a subsonic inlet at arbitrary operating conditions of mass flow rate, free stream velocity and incidence angle. Both laminar separation on the lip and turbulent separation in the diffuser are discussed. The agreement of the theoretical results with model experimental data illustrates the capability of the theory to predict separation. The effects of throat Mach number, inlet size, and surface roughness on boundary layer development and separation are illustrated.

Chou, D. C.↗
NASA NTRS1989

Proposed uses of laser light scattering instruments for polymerization studies

Microgravity offers a unique environment for studying polymer diffusion and polymer polymerization reactions. The absence of convection currents, which are the major mode of mixing at the molecular level on Earth, are eliminated or reduced in the microgravity environment. More importantly, the prediction of unique copolymer composition development in microgravity allows controlled formation of new compositions of matter. The absence of mixing at the molecular level should produce unique short block copolymers available for the first time for comonomer compositions which normally lead to random or long block copolymer under good mixing. The investigation of fundamental polymer diffusion and polymer polymerization processes in microgravity is proposed. This effort will involve fundamental studies of monomer and polymer diffusion; their effects on initiation, propagation, and especially termination kinetics rate constant; and the accurate evaluation of copolymerization reactivity ratios in microgravity. The experimental design is presented for these studies along with an evaluation technique for in situ monitoring of polymer diffusion and polymerization kinetics.

Mathias, Lon J.↗