Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Principles of robust innate immune response to viral infections: a multiplex network analysis

The human innate immune response, particularly the type-I interferon (IFN) response, is highly robust and effective first line of defense against virus invasion. IFN molecules are produced and secreted frominfected cells upon virus infection and recognition. They then act as signaling/communicationmolecules to activate an antiviral response in neighboring cells so that those cells become refractory to infection. Previous experimental studies have identified the detailed molecular mechanisms for the IFN signaling and response. However, the principles underlying how host cells use IFN to communicate with each other to collectively and robustly halt an infection is not understood. Here we take a multiplex network modeling approach to provide a theoretical framework to identify key factors that determine the effectiveness of the IFN response against virus infection of a host. In this approach, we consider the virus spread among host cells and the interferon signaling to protect host cells as a competition process on a two-layer multiplex network. We focused on two types of network topology, i.e., the Erdős-Rényi (ER) network and the Geometric Random (GR) network, which represent the scenarios when infection of cells is mostly well mixed (e.g., in the blood) and when infection is spatially segregated (e.g., in tissues), respectively.We show that in general, the IFN response works effectively to stop viral infection when virus infection spreads spatially (a most likely scenario for initial virus infection of a host at the peripheral tissue). Importantly, we show that the effectiveness of the IFN response is robust against large variations in the distance of IFN diffusion as long as IFNs diffuse faster than viruses and they can effectively induce antiviral responses in susceptible host cells. This suggests that the effectiveness of the IFN response is insensitive to the specific arrangement of host cells in peripheral tissues. Thus, our work provides a quantitative explanation of why the IFN response can serve an effective and robust response in different tissue types to a wide range of viral infections of a host.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Effect of Ceramic Rolling and Annealing on Mechanical Properties of AlCoCrFeNi{sub 2.1} Eutectic High-Entropy Alloys

It was firstly proved that AlCoCrFeNi{sub 2.1} could optimize the strength-ductility combination by ceramic rolling at room temperature. The technique greatly enhanced the strength of the materials at the price of slight ductility, which is different from the previous studies. When the thickness reduction was 23.06%, the yield strength and compressive strength were 1611.42 and 2212.99 MPa, respectively, with the plastic strain of 28.36%. Such an excellent combination of strength and ductility was attributed to the distribution of the residual stress, the high hydrostatic stress and the shear bands in ceramic rolling process. The plastic strain of the material increased to 38.54%, and the strength reduced to 1900.73 MPa after annealing at 800 °C for 1 h. A lot of like-flower patterns of recrystallized ultrafine grain were observed after annealing at 800 °C due to the sluggish diffusion in HEAs and the effect of ceramic rolling. In addition, after annealing at 800 °C, the fracture mode changed to shearing fracture of microvoid aggregation from cleavage/quasi-cleavage fracture of as-rolled sample due to the deposition of the second phase in the tested sample.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Effect of Long-Term 1093 K Exposure to Air or Vacuum on the Structure of Several Wrought Superalloys

Long-term 1093 K heat treatments of three commercial superalloy sheet materials were undertaken in air -IAI and vacuum. With either exposure, significant precipitation of second phases occurred in the Co-base Haynes(R) Alloy 188 (HA 188) and the Ni-base Haynes(R) Alloy 230 (HA 230); however, much less precipitation was found in the exposed Ni-base alloy Inconel(R) 617 (IN 617). Although some grain growth occurred in HA 198, no changes in the grain size of either HA 230 or IN 617 were observed after 22,500 h at temperature. Oxidation during air heat treatments led to weight gain due to the formation of chromia + spinel scales and surface-connected grain boundary pits/oxides in all three superalloys. Both the weight gain and depth of intergranular attack were dependent on the square root of time, which is indicative of diffusion-controlled phenomena. Because many alloy samples had neighbors in close proximity, most vacuum heat treated specimens did not suffer significant loss of volatile elements. However, some exposed samples were subjected to unrestricted vacuum heat treatments, allowing estimates of volatilization to be made. Based on the data for HA 188, the weight loss during 1093 K vacuum exposure was diffusion controlled once the inhibiting effects of surface films on the as-received alloys were broken down.

J.D. Whittenberger↗
DOE OSTI2018

Effect of Heat Treatment on Madagascar Dravite Tourmaline: UV-Visible and Diffuse Reflectance Infrared Spectroscopic Characterization

The color change and chemical behavior of Madagascar dravite when subjected to heating in air atmosphere were studied by energy-dispersive X-ray fluorescence spectrometry, UV-visible (UV-vis) spectroscopy and diffuse reflectance infrared Fourier transform (DRIFT) spectroscopy. Heat treatment was performed with temperatures of 500, 600, and 700°C for 2 h. Colors and color differences were measured and evaluated using CIELAB color measurement. The origin of color and its change after treatment are assessed via UV-vis spectra. The diffuse reflectance infrared spectra of the dravites revealed major absorption bands of OH, BO{sub 3}, Si{sub 6}O{sub 18}, and R–OH (where R = Fe, Mg). The variation of infrared spectra with increasing temperature indicates the deformation of R–OH bonds in an octahedron and, further, the deformation of BO{sub 3} and Si{sub 6}O{sub 18}. The results of color change may be utilized to enhance the color or clarity in dravite tourmaline and for generating satisfactory color.

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

Effects of irradiation on the interface between U-Mo and zirconium diffusion barrier

Irradiated fuel plates were characterized by microscopy that focused on the interface between U–Mo and Zr. Before irradiation, there were three major sub-layers identified in the U–Mo/Zr interface, namely, UZr 2 , Mo 2 Zr, and U with low Mo. The typical total thickness of this U–Mo/Zr interaction is 2–3 microns. The UZr 2 sub-layer formed during fuel fabrication remains stable after irradiation, without large bubbles/porosity accumulation. However, this sub-layer becomes increasingly discontinuous as burnup increases. The low-Mo sub-layer exhibits numerous sub-micron bubbles/porosity at low burnup. Larger, interconnected porosity in this sub-layer was observed in a medium-burnup fuel specimen. However, at higher burnup, regions with the extra-large bubbles/porosity (i.e., larger than 5 ?m) were observed at least 5 microns away from the original location of this sub-layer. The mechanism for the formation of the extra-large bubbles/porosity is still unclear at this time. In general, the U–Mo/Zr interface in monolithic U–Mo fuels is relatively stable after irradiation. As a result, no large detrimental defects, such as large interfacial bubbles or cracks/delamination, were observed in the fuel plates characterized.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1981

Double-beam optical method and apparatus for measuring thermal diffusivity and other molecular dynamic processes in utilizing the transient thermal lens effect

A sample material was irradiated by relatively high power, short pulses from a dye laser. Energy from the pulses was absorbed by the sample material, thereby forming a thermal lens in the area of absorption. The pulse repetition rate was chosen so that the thermal lens is substantially dissipated by the time the next pulse reaches the sample material. A probe light beam, which in a specific embodiment is a relatively low power, continuous wave (cw) laser beam, irradiated the thermal lens formed in the sample material. The intensity characteristics of the probe light beam subsequent to irradiation of the thermal lens is related to changes in the refractive index of the sample material as the thermal lens is formed and dissipated.

Moacanin, J.↗
NASA NTRS2018

Effect of the Reconnection Electric Field on Electron Distribution Functions in the Diffusion Region of Magnetotail Reconnection

Electron distribution functions in the electron diffusion region during symmetric magnetic reconnection are investigated by means of theory and fully kinetic simulations. Crescent‐like striations are formed in distribution functions in the velocity plane perpendicular to the magnetic field. Using an analytical current sheet, we solve the equation of motion for electrons, and derive the shape of a crescent distribution, as a function of the distance from the neutral line, field gradients, and the reconnection electric field. Each crescent is tilted in the velocity plane because of the acceleration by the reconnection electric field, and multiple stripes appear due to multiple meandering bounces. Applying the theory to distribution functions observed in Earth's magnetotail, we deduce the amplitude of the reconnection electric field.

Bessho, N.↗
DOE OSTI2019

Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The DSA scheme has been implemented in PROTEUS to speed up SN transport calculations. The existing DSA scheme employs a consistent finite element formulation on the same fine mesh structure as the SN solution While the use of consistent discretizations makes the DSA method effective, it does not necessarily make it efficient because the numerical solution of a discretized elliptic diffusion problem itself can be costly, particularly when the problem size becomes large. We developed a new discontinuous Galerkin (DG) discretization of the diffusion equation, called DG-DSA, which can effectively and efficiently accelerate the SN transport iterations. As compared with the previous work, the novelty of our method is that the diffusion equation is solved on a coarse-mesh grid using the DG methods, and the DG diffusion discretization incorporates local hp adaptation, i.e., local adaptation of mesh size and/or polynomial degree, based on local total cross section (or optical thickness). Therefore, the resulting number of degrees of freedom (DOF) of the DG discretization is much less than the conventional consistent DSA discretizations, and thus DG-DSA can achieve significant improvement in computational efficiency. We implemented this scheme in PROTEUS-SN. In addition, we developed a new stabilization scheme named linear prolongation CMFD (lpCMFD). A novel feature of this scheme is that the conventional flat flux ratio–based scaling approach is replaced with a linear interpolation of the scalar flux differences at the coarse-mesh cell edges between the neutron transport and CMFD calculations. Fourier convergence analysis and numerical results show that lpCMFD is unconditionally stable and more effective than other CMFD based on acceleration schemes such as pCMFD, odCMFD. We also developed a new nonlinear diffusion acceleration scheme for solving neutron transport equations. This scheme, called LR-NDA, employs a local refinement approach on the framework of CMFD by solving a local boundary value problem of the scalar flux on the coarse-mesh structure to replace the piecewise constant scalar flux obtained by CMFD. The refined flux is then used to update the scalar flux in the neutron transport source iteration. It has been demonstrated that that LR-NDA is much more effective and stable than CMFD for a wide range of optical thicknesses. LR-NDA is a local adaptive method, which means LR-NDA does not necessarily require local refinement for all the coarse-mesh cells on the problem domain, i.e., it can be used only for relatively optically thick regions where the standard CMFD scheme would encounter the convergence problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination

Here in this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Immobile Pore-Water Storage Enhancement and Retardation of Gas Transport in Fractured Rock

The effect of immobile pore-water on gas transport in fractured rock has implications for numerical modeling of soil vapor extraction, methane leakage, gaseous CO 2 leakage from sequestration operations, radionuclide transport from underground nuclear explosions, and nuclear waste disposal. While the ability for immobile pore-water storage to effect gas transport has been recognized in the past, the details and specific scenarios leading to enhanced, retarded, or unaffected gas transport have not been explored. In this paper, we performed numerical investigations into the enhancement and retardation of gas transport due to immobile pore-water storage in order to identify implications for gas transport applications. To do this, we developed a numerical approach to model gas transport with a single-phase flow solution coupled to the advection–dispersion equation modified to account for immobile pore-water storage. Other than the immobility of pore water, the formulation contains all other physics included in two-phase formulations (advective and diffusive gas transport in fractures and rock matrix and dissolution in immobile pore water). The assumption of immobile pore water is valid here since for many applications involving transport of soluble gases in fractured rock, the rate of aqueous transport is insignificant compared to gas transport. We verify our modeling approach with analytical solutions of: (1) 1D gas diffusion, (2) 1D gas advection, (3) barometric pumping of a fracture, and (4) gas transport with uniform fracture flow and diffusion into the matrix. We account for pore-water storage in our model by implementing a kinetic formulation of gas dissolution wherein the dissolved (aqueous) phase is considered an immobile constituent. Using this formulation, we model the effect of dissolution rate and saturation on the retardation of gas transport during pure diffusion and pure advection. We also demonstrate that although it is commonly believed that pore-water storage will always enhance gas transport in fractures during oscillatory flow (e.g., during reversing pressure gradients such as barometric pumping cycles), our simulations indicate that this may not always be the case. Our numerical investigations indicate that scenarios with lower effective diffusion coefficients ( &#x2A85;<!-- <≈ --> 10 &#x2212;<!-- - --> 5 &#xA0;m 2 / s ) and lower dissolution coefficients (i.e., the dissolution diffusion coefficient) ( &#x2A85;<!-- <≈ --> 10 &#x2212;<!-- - --> 11 &#xA0;m 2 / s ) will result in enhanced gas transport. Other combinations of gaseous diffusion and dissolution coefficients result in delayed gas transport or insignificant effect on gas transport. Based on these results, and given the fact that the free-air diffusion coefficients for many gases are slightly above 10 &#x2212;<!-- - --> 5 &#xA0;m 2 / s (near the boundary of enhanced/non-enhanced gas transport), tortuosity would have to be significant for gas transport enhancement to occur. Similarly, if we consider the free-water diffusion coefficient of gases (around 10 &#x2212;<!-- - --> 9 &#xA0;m 2 / s for many gases of interest) to be a maximum bound for the effective dissolution coefficients (i.e., the rate of diffusion into the water away from the air/water interface will limit the dissolution rate), the effective dissolution coefficient would have to be at least 2 orders of magnitude less for enhanced gas transport to occur. Otherwise, pore-water storage will delay or have negligible effect on the gas transport. The results also indicate that the assumption of instantaneous equilibrium, often invoked in numerical codes for dissolution processes (i.e., dissolution coefficient is effectively infinite), will fail to accurately capture important details of soluble gas transport. Finally, the implication of our results is that gas breakthrough times in fractured rock during oscillatory flow may be highly dependent on the effective gaseous diffusion and dissolution coefficients.

58 GEOSCIENCES↗
DOE OSTI2021

Short-range thermal magnon diffusion in magnetic garnet

Here, using the spin Seebeck effect (SSE), we study the propagation distance of thermally induced spin currents inside a magnetic insulator thin film in the short-range regime. We disambiguate spin currents driven by temperature and chemical potential gradients by comparing the SSE signal before and after adding a heat-sinking capping layer on the same device. We report that the measured spin decay behavior near the heat source is well accounted for by a diffusion model where the magnon diffusion length is in submicron range, in other words, two orders of magnitude smaller than previous estimates inferred from the long-range behavior. Our results highlight the caveat in applying a diffusive theory to describe thermally generated magnon transport, where a single decay length may not capture the behavior on all length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Size dependent behavior of Fe 3 O 4 crystals during electrochemical (de)lithiation: an in situ X-ray diffraction, ex situ X-ray absorption spectroscopy, transmission electron microscopy and theoretical investigation

Here, the iron oxide magnetite, Fe 3 O 4 , is a promising conversion type lithium ion battery anode material due to its high natural abundance, low cost and high theoretical capacity. While the close packing of ions in the inverse spinel structure of Fe 3 O 4 enables high energy density, it also limits the kinetics of lithium ion diffusion in the material. Nanosizing of Fe 3 O 4 to reduce the diffusion path length is an effective strategy for overcoming this issue and results in improved rate capability. However, the impact of nanosizing on the multiple structural transformations that occur during the electrochemical (de)lithiation reaction in Fe 3 O 4 is poorly understood. In this study, the influence of crystallite size on the lithiation-conversion mechanisms in Fe 3 O 4 is investigated using complementary X-ray techniques along with transmission electron microscopy (TEM) and continuum level simulations on electrodes of two different Fe 3 O 4 crystallite sizes. In situ X-ray diffraction (XRD) measurements were utilized to track the changes to the crystalline phases during (de)lithiation. X-ray absorption spectroscopy (XAS) measurements at multiple points during the (de)lithiation processes provided local electronic and atomic structural information. Tracking the crystalline and nanocrystalline phases during the first (de)lithiation provides experimental evidence that (1) the lithiation mechanism is non-uniform and dependent on crystallite size, where increased Li + diffusion length in larger crystals results in conversion to Fe 0 metal while insertion of Li + into spinel-Fe 3 O 4 is still occurring, and (2) the disorder and size of the Fe metal domains formed when either material is fully lithiated impacts the homogeneity of the FeO phase formed during the subsequent delithiation.

25 ENERGY STORAGE↗
NASA NTRS1982

Diffusive-convective physical vapor transport of PbTe from a Te-rich solid source

Crystal growth of PbTe by physical vapor transport (sublimation) in a closed ampoule is governed by the vapor species in thermal equilibrium with the solid compound. Deviations from stoichiometry in the source material cause diffusion limitation of the transport rate, which can be modified by natural (gravity-driven) convection. Mass-transport experiments have been performed using Te-rich material wherein sublimation rates have been measured in order to study the effects of natural convection in diffusion-limited vapor transport. Linear velocities for both crystal growth and evaporation (back sublimation) have been measured for transport in the direction of gravity, horizontally, and opposite to gravity. The experimental results are discussed in terms of both the one-dimensional diffusive-advective model and current, more sophisticated theory which includes natural convection. There is some evidence that convection effects from radial temperature gradients and solutal density gradients have been observed.

Zoutendyk, J.↗
NASA NTRS1992

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗
NASA NTRS1993

Diffusion of Na and K in the uppermost regolith of Mercury

We examine the conditions under which ions impacting the Hermean surface can act as a regional source of enhanced atmospheric column through ion implantation and subsequent release at the surface. At most latitudes, energetic (a few keV) Na or K ions which impact the nightside surface, are released quickly (well before noon) upon warming of the surface. We show that the relative sunrise/sunset difference produced by ion implantation is ((n)(sub SR) - (n)(sub SS)/(n)(sub Ave) = f, where f is the fraction of the photo-ions recycled, (n)(sub SR) is the average zenith column above the sunrise portion of the illuminated hemisphere, (n)(sub SS) is the average zenith column in the sunset portion, and (n)(sub Ave) is the average column over the sunlit disk. Thus, to produce a large sunrise/sunset difference via ion implantation and subsequent release requires efficient (close to total) recycling. We show that the most extensive set of available data reduced to Na column abundance does not show any sunrise enhancement. We argue that the K data do not permit an unambiguous interpretation in favor of sunrise/sunset differences. We further find that if an efficient surface loss process for the alkali is not operating after sunrise, the initial Na and K distribution will relax into the bulk of the solid. Preserved abundance gradients of Na and K in lunar glasses that the Arrhenius coefficients for impact glasses are likely more modest than those for laboratory glasses derived from rock by nonimpact processes, but the effect of the more modest diffusion rates is only to delay the efficient loss of Na by a few Earth days. We argue that implantation can lead to observable regional increases in the observed Na or K column densities only if it occurs at very high latitudes, where diffusion is slow. It is typically lower energy ions which impact at high latitude and these are both more numerous than the high energy ions and possess smaller average penetration depths; thus there are additional reasons to favor a high-latitude locus for any possible prompt return related increases in zenith column. We find that we cannot rule out sputtering as a source process, as the low-energy ions are efficient sputterers. Finally, we argue that the observed Na/K ratio in the atmopshere may be the expression of their different asymptotic rates of loss from the interiors of the regolith grains. If so, this allows us to fix the importance of sputtering and photon stimulated desorption relative to impact vaporization.

Killen, R. M.↗