Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗
DOE OSTI2023

Machine learning-based ethylene and carbon monoxide estimation, real-time optimization, and multivariable feedback control of an experimental electrochemical reactor

Electrochemical reduction of CO 2 gas is a novel CO 2 utilization technique that has the potential to mitigate the global climate crisis caused by anthropogenic CO 2 emissions, and enable the large-scale storage of energy generated from renewable sources in the form of carbon-based chemicals and fuels. However, due to the complexity of the electrochemical reactions, the explicit first-principles models for CO2 reduction are not available yet, and there has been a limited effort to develop process modeling, optimization and control of CO 2 electrochemical reactors. To this end, a rotating cylinder electrode (RCE) reactor has been constructed at UCLA to understand the mass transfer and reaction kinetics effects separately on the productivity. In the RCE reactor, the applied potential strongly influences the reaction energetics and the electrode rotation speed affects the hydrodynamic boundary layer and modifies the film mass transfer coefficient, which involves convective and diffusive transport. Further, the present work aims to develop a multi-input multi-output (MIMO) control scheme for the RCE reactor that integrates techniques from artificial and recurrent neural network modeling, nonlinear optimization, and process controller design. Specifically, production rates of two products from the experimental reactor, ethylene and carbon monoxide, are controlled by manipulating two inputs, applied potential and catalyst rotation speed. Process dynamics and controllability are analyzed, a feedback control strategy is designed and the controllers are tuned accordingly. The experimental electrochemical cell is employed to gather data for process modeling and implement the multivariable control system. Finally, the experimental results are presented which demonstrate excellent closed-loop performance by the control system and regulation of the outputs at three different set-points including an economically-optimal set-point.

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

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cosmic-ray transport in inhomogeneous media

A theory of cosmic-ray transport in multiphase diffusive media is developed, with the specific application to cases in which the cosmic-ray diffusion coefficient has large spatial fluctuations that may be inherently multiscale. We demonstrate that the resulting transport of cosmic rays is diffusive in the long-time limit, with an average diffusion coefficient equal to the harmonic mean of the spatially varying diffusion coefficient. Thus, cosmic-ray transport is dominated by areas of low diffusion even if these areas occupy a relatively small, but not infinitesimal, fraction of the volume. On intermediate time-scales, the cosmic rays experience transient effective subdiffusion, as a result of low-diffusion regions interrupting long flights through high-diffusion regions. In the simplified case of a two-phase medium, we show that the extent and extremity of the subdiffusivity of cosmic-ray transport is controlled by the spectral exponent of the distribution of patch sizes of each of the phases. We finally show that, despite strongly influencing the confinement times, the multiphase medium is only capable of altering the energy dependence of cosmic-ray transport when there is a moderate (but not excessive) level of perpendicular diffusion across magnetic-field lines.

Astronomy and Astrophysics↗
NASA NTRS2024

Toward 2D Dynamo Models Calibrated by Global 3D Relativistic Accretion Disk Simulations

Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.

Matthew D Duez↗
DOE OSTI2021

Solvation Structure and Dynamics of Mg(TFSI) 2 Aqueous Electrolyte

Using ab initio molecular dynamics (AIMD) simulations, classical molecular dynamics (CMD) simulations, small-angle X-ray scattering (SAXS), and pulsed-field gradient nuclear magnetic resonance (PFG-NMR), the solvation structure and ion dynamics of magnesium bis(trifluoromethanesulfonyl)imide (Mg(TFSI) 2 ) aqueous electrolyte at 1, 2, and 3 m concentrations are investigated. From AIMD and CMD simulations, the first solvation shell of an Mg 2+ ion is found to be composed of six water molecules in an octahedral configuration and the solvation shell is rather rigid. The TFSI - ions prefer to stay in the second solvation shell and beyond. Meanwhile, the comparable diffusion coefficients of positive and negative ions in Mg(TFSI) 2 aqueous electrolytes have been observed, which is mainly due to the formation of the stable [Mg(H 2 O) 6 ] 2+ complex, and, as a result, the increased effective Mg ion size. Finally, the calculated correlated transference numbers are lower than the uncorrelated ones even at the low concentration of 2 and 3 m, suggesting the enhanced correlations between ions in the multivalent electrolytes. This work provides a molecular-level understanding of how the solvation structure and multivalency of the ion affect the dynamics and transport properties of the multivalent electrolyte, providing insight for rational designs of electrolytes for improved ion transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen diffusion induced dislocation transformations in a nickel superalloy

The diffusion of hydrogen in metals and alloys induces embrittlement that can adversely affect the structural properties. We examine the adsorption and diffusion of hydrogen in Inconel-718 (IN-718), and scrutinize the ensuing effects on the dislocation behavior in the alloy to elucidate the fundamental mechanisms of hydrogen-microstructure interactions from classical molecular simulations. Hydrogen adsorption increases with time until the surface saturates, while hydrogen diffusion exhibits strong temperature dependence, with diffusion coefficients converging above 1300 K regardless of the initial hydrogen concentration in the alloy. The diffusion in IN-718 is significantly sluggish than in pure Ni, Fe, or Cr, and is strongly impacted by hydrogen concentrations, resulting in an order of magnitude higher diffusion coefficient for hydrogen (10-14 m2/s relative to 10-15 m2/s) at high concentrations, especially below 600 K. Hydrogen diffusion coefficient varies from 10-12 to 10-15 m2/s in IN-718 depending on temperature (500–1400 K). More critically, our results reveal that increasing hydrogen concentration induces microstructural changes in the alloy, transforming perfect dislocations into stair-rods and Shockley partials, with higher temperatures favoring the latter. The results are significant for hydrogen fuel applications to gain insights into the materials chemistry for designing safer and more efficient propulsion systems, particularly in high-performance environments related to controlled hydrogen combustion applications.

Banerjee, Tanumoy↗
DOE OSTI2026

Charge Carrier and Spin Diffusion in a Polycrystalline and Single Crystal Lead Halide Perovskite Semiconductor

The directed movement of spin-bearing excitations in semiconductors can enable many potential schemes for spintronic applications. Understanding the mechanism of spin motion, as opposed to charge carrier or exciton motion, may require specialized techniques and sample conditions. Here we employ the noncontact and time-resolved technique of light-induced transient grating spectroscopy (LITG) to measure both carrier and spin motion in a perovskite semiconductor with varying crystallinity. The carrier motion aligns with expectations from past studies undertaken at low fluence, revealing an ambipolar diffusion coefficient on the order of 1 cm 2 /sec and diffusion length of roughly 1.5 μm for single crystals that is strongly attenuated as crystallite size decreases. The spin diffusion is measured with cross-polarized excitation beams and uncovers an intensity-dependent coefficient rising above 100 cm 2 /sec. The fast room-temperature spin relaxation limits the spin diffusion length, but the remarkable speedup of spin over carrier diffusion suggests a mechanism involving a combination of exchange-mediated and doping effects that enhance spin transport.

14 SOLAR ENERGY↗
NASA NTRS2023

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗
NASA NTRS2000

Effect of Film-Hole Shape on Turbine Blade Film Cooling Performance

The detailed heat transfer coefficient and film cooling effectiveness distributions as well as tile detailed coolant jet temperature profiles on the suction side of a gas turbine blade A,ere measured using a transient liquid crystal image method and a traversing cold wire and a traversing thermocouple probe, respectively. The blade has only one row of film holes near the gill hole portion on the suction side of the blade. The hole geometries studied include standard cylindrical holes and holes with diffuser shaped exit portion (i.e. fanshaped holes and laidback fanshaped holes). Tests were performed on a five-blade linear cascade in a low-speed wind tunnel. The mainstream Reynolds number based on cascade exit velocity was 5.3 x 10(exp 5). Upstream unsteady wakes were simulated using a spoke-wheel type wake generator. The wake Strouhal number was kept at 0 or 0.1. Coolant blowing ratio was varied from 0.4 to 1.2. Results show that both expanded holes have significantly improved thermal protection over the surface downstream of the ejection location, particularly at high blowing ratios. However, the expanded hole injections induce earlier boundary layer transition to turbulence and enhance heat transfer coefficients at the latter part of the blade suction surface. In general, the unsteady wake tends to reduce film cooling effectiveness.

Han, J. C.↗
NASA NTRS1986

Losses in Channels with Increased External Turbulence

An approximate method for determining the effect of the level of turbulence on the aerodynamic characteristics of convergent and diffuser channels is examined. A momentum equation for the boundary layer is in the method, introducing external flow turbulence on the basis of experimental values of the coefficient of friction and the form factor. It is found that at significant levels of external turbulence, losses must be considered not only in the boundary layer but also in the central region of the channel.

Zaryankin, A. Y.↗
NASA NTRS1990

Solar Response to Luminosity Variations

The connection between solar luminosity and magnetic fields is now well-established. Magnetic fields under the guise of sunspots and faculae enhance or suppress heat transfer through the solar surface, leading to changes in the total solar luminosity. This raises the question of the effect that such surface heat transfer perturbations have on the internal structure of the sun. The problem has been considered previously by Foukal and Spruit. Here, researchers generalize the calculation of Spruit, removing the assumption of a constant heat diffusivity coefficient by treating the full mixing length heat transfer expression. Further, they treat the surface conditions in a simpler manner, and show that the previous conclusions of Foukal and Spruit are unaffected by these modifications. The model shows that following the application of a step function emissivity change: a fraction 1 - D(sub 0) of the luminosity change relaxes away after approx. 50 days. This corresponds to the thermal diffusion time across the convection zone, adjusting to a value in correspondence with the surface change. In other words, the whole convection zone feels the perturbation on this timescale. The remaining fraction relaxes away on a timescale of 10 to the 5th power years, corresponding to the convective layer radiating away enough energy so that it can adjust to its new adiabat. These are the same results arrived at by Spruit and Foukal. For variations of sigma on timescales of 10 to 200 years, then, the only important relaxation is the 50 day one. If the amplitude of this relaxation is small, the luminosity follows the sigma variation.

Arendt, S.↗
NASA NTRS1987

Influence of oscillation-induced diffusion on heat transfer in a uniformly heated channel

An analysis of the effect of flow oscillations on laminar flow heat transfer in a channel with uniform heat addition is presented. It is shown that the effect of flow oscillations will be to reduce the channel heat transfer coefficient. This effect is due to the fact that the heat addition along the channel wall produces an increasing fluid temperature along the channel length. The flow oscillations interacting with this positive temperature gradient will induce a heat flow back toward the channel inlet. This will tend to inhibit the heat transfer process and will raise the wall temperature required to transfer away a given amount of heat at the channel wall.

Siegel, R.↗
NASA NTRSDate not supplied

Numerical Investigation of Film Coefficient Engineering Methodology for Dissociated, Chemically Reacting Boundary Layers

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flow field for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for their application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. The sensitivity of the engineering model is assessed based on the impact of the CFD assumptions as an effort to make it available to design engineers across industry.

numerical boundary layer↗
DOE OSTI2024

Resolving electron and hole transport properties in semiconductor materials by constant light-induced magneto transport

The knowledge of minority and majority charge carrier properties enables controlling the performance of solar cells, transistors, detectors, sensors, and LEDs. Here, we developed the constant light induced magneto transport method which resolves electron and hole mobility, lifetime, diffusion coefficient and length, and quasi-Fermi level splitting. We demonstrate the implication of the constant light induced magneto transport for silicon and metal halide perovskite films. We resolve the transport properties of electrons and holes predicting the material’s effectiveness for solar cell application without making the full device. The accessibility of fourteen material parameters paves the way for in-depth exploration of causal mechanisms limiting the efficiency and functionality of material structures. To demonstrate broad applicability, we further characterized twelve materials with drift mobilities spanning from 10 –3 to 10 3 cm 2 V –1 s –1 and lifetimes varying between 10 –9 and 10 –3 seconds. The universality of our method its potential to advance optoelectronic devices in various technological fields.

14 SOLAR ENERGY↗
NASA NTRS1982

Adsorption and excess fission xenon

The adsorption of Xe and Kr on lunar soil 10084 was measured by a method that employs only very low fractions of monolayer coverage. Results are presented as parameters for calculation of the Henry constant for adsorption as a function of temperature. The adsorption potentials are about 3 kcal/mole for Kr and 5 kcal/mole for Xe; heating the sample in vacuum increased the Xe potential to nearly 7 kcal/mole. Henry constants at the characteristic lunar temperature are about 0.3 cu cm STP/g-atm. These data were applied to consider whether adsorption is important in producing the excess fission Xe effect characteristic of highland breccias. Sorption equilibrium with a transient lunar atmosphere vented fission Xe produces concentrations seven orders of magnitude lower than observed concentrations. Higher concentrations result because of the resistance of the regolith to upward diffusion of Xe. A diffusion coefficient of 0.26 sq cm/sec is estimated for this process.

Podosek, F. A.↗
NASA NTRS1992

Flow-induced gelation of living (micellar) polymers

The effect of shear velocity gradients on the size (L) of rodlike micelles in dilute and semidilute solution is considered. A kinetic equation is introduced for the time-dependent concentration of aggregates of length L, consisting of 'bimolecular' combination processes L + L-prime yield (L + L-prime) and unimolecular fragmentations L yield L + (L - L-prime). The former are described by a generalization (from spheres to rods) of the Smoluchowski mechanism for shear-induced coalesence of emulsions, and the latter by incorporating the tension-deformation effects due to flow. Steady-state solutions to the kinetic equation are obtained, with the corresponding mean micellar size evaluated as a function of the Peclet number P (i.e., the dimensionless ratio of the flow rate and the rotational diffusion coefficient). For sufficiently dilute solutions, only a weak dependence of the micellar size on P is found. In the semidilute regime, however, an apparent divergence in the micellar size at P of about 1 suggests a flow-induced first-order gelation phenomenon.

Bruinsma, Robijn↗