Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A comparative study of discrete fracture network and equivalent continuum models for simulating flow and transport in the far field of a hypothetical nuclear waste repository in crystalline host rock

One of the major challenges of simulating flow and transport in the far field of a geologic repository in crystalline host rock is related to reproducing the properties of the fracture network over the large volume of rock with sparse fracture characterization data. Various approaches have been developed to simulate flow and transport through the fractured rock. The approaches can be broadly divided into Discrete Fracture Network (DFN) and Equivalent Continuum Model (ECM). The DFN explicitly represents individual fractures, while the ECM uses fracture properties to determine equivalent continuum parameters. In this paper, we compare DFN and ECM in terms of upscaled observed transport properties through generic fracture networks. The major effort was directed on making the DFN and ECM approaches similar in their conceptual representations. This allows for separating differences related to the interpretation of the test conditions and parameters from the differences between the DFN and ECM approaches. The two models are compared using a benchmark test problem that is constructed to represent the far field (1 × 1 × 1 km 3 ) of a hypothetical repository in fractured crystalline rock. The test problem setting uses generic fracture properties that can be expected in crystalline rocks. The models are compared in terms of the: 1) effective permeability of the domain, and 2) nonreactive solute breakthrough curves through the domain. The principal differences between the models are mesh size, network connectivity, matrix diffusion and anisotropy. We demonstrate how these differences affect the flow and transport. Finally, we identify the factors that should be taken in consideration when selecting an approach most suitable for the site-specific conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Positivity enhancements to the P n transport equations

One of the common methods for solving the radiation transport equation is to use a polynomial expansion for the angle variable(s). Here, while each solution technique has its disadvantages, in some instances this method can result in oscillations large enough to give a negative radiation energy density, which is nonphysical. Some authors have recast the closure of the transport equations to guarantee positivity. Unfortunately, this results in messy nonlinear equations. Other authors post process each time step’s solution to improve positivity. In this work, the P n equations are modified to be more diffusive, less oscillatory, more positive, and more accurate. No new solution techniques are required, just the scaling of some terms with constant factors and the addition of constant filter opacities; therefore, the solutions are still linear. This new solution method is called stretched and filtered P n . Test problems are presented in one and two dimensions. Comparisons to the D n and TP n closures are made in one dimension. Modifications to D n and TP n are suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Investigating the Effects of Anisotropic Mass Transport on Dendrite Growth in High Energy Density Lithium Batteries

Dendrite formation on the electrode surface of high energy density lithium (Li) batteries causes safety problems and limits their applications. Suppressing dendrite growth could significantly improve Li battery performance. Dendrite growth and morphology is a function of the mixing in the electrolyte near the anode interface. Most research into dendrites in batteries focuses on dendrite formation in isotropic electrolytes (i.e., electrolytes with isotropic diffusion coefficient). In this work, an anisotropic diffusion reaction model is developed to study the anisotropic mixing effect on dendrite growth in Li batteries. The model uses a Lagrangian particle-based method to model dendrite growth in an anisotropic electrolyte solution. The model is verified by comparing the numerical simulation results with analytical solutions, and its accuracy is shown to be better than previous particle-based anisotropic diffusion models. Several parametric studies of dendrite growth in an anisotropic electrolyte are performed and the results demonstrate the effects of anisotropic transport on dendrite growth and morphology, and show the possible advantages of anisotropic electrolytes for dendrite suppression.

SPH, high energy density lithium batteries, precip↗
NASA NTRS1987

The physics of interstellar shock waves

This review discusses the observations and theoretical models of interstellar shock waves, in both diffuse cloud and molecular cloud environments. It summarizes the relevant gas dynamics, atomic, molecular and grain processes, radiative transfer, and physics of radiative and magnetic precursors in shock models. It then describes the importance of shocks for observations, diagnostics, and global interstellar dynamics. It concludes with current research problems and data needs for atomic, molecular and grain physics.

Shull, J. Michael↗
NASA NTRS1990

A cost-effective strategy for nonoscillatory convection without clipping

Clipping of narrow extrema and distortion of smooth profiles is a well known problem associated with so-called high resolution nonoscillatory convection schemes. A strategy is presented for accurately simulating highly convective flows containing discontinuities such as density fronts or shock waves, without distorting smooth profiles or clipping narrow local extrema. The convection algorithm is based on non-artificially diffusive third-order upwinding in smooth regions, with automatic adaptive stencil expansion to (in principle, arbitrarily) higher order upwinding locally, in regions of rapidly changing gradients. This is highly cost effective because the wider stencil is used only where needed-in isolated narrow regions. A recently developed universal limiter assures sharp monotonic resolution of discontinuities without introducing artificial diffusion or numerical compression. An adaptive discriminator is constructed to distinguish between spurious overshoots and physical peaks; this automatically relaxes the limiter near local turning points, thereby avoiding loss of resolution in narrow extrema. Examples are given for one-dimensional pure convection of scalar profiles at constant velocity.

Leonard, B. P.↗
NASA NTRS1991

Assessment of a prototype sunphotometer for network applications, phase 1

Data taken in Tucson, AZ, for comparing the CE sunphotometer's results against the Univ. of Arizona's no. 3 solar radiometer are in very good agreement. The spectral voltage intercepts obtained on 15 different days taken with the CE instrument during a seven week period are stable and follow the same daily trends as the voltage intercepts for similar wavelengths contained in the Univ. of Arizona's no. 3 radiometer. The spectral optical depths obtained from both instruments are also in very good agreement and follow the same trends. The only problem that exists in the CE instrument does not have a heater to stabilize the detector temperature which is important around 1000 nm. Zenith sky radiance scans are in good agreement, but both instruments are biased by direct solar rays entering the instrument's respective entrance apertures. All diffuse light data shown were, therefore, taken with the instruments shielded from the direct sun. The CE sunphotometer results pertaining to the water vapor wavelength's are also in good agreement with the Univ. of Arizona's no. 3 and should prove useful for water vapor retrieval.

Reagan, J. A.↗
DOE OSTI2019

Hybrid super homogenization and discontinuity factor method for continuous finite element diffusion

This paper presents a novel homogenization equivalence technique aiming to simultaneously leveragethe simplicity of the Super Homogenization (SPH) method to reproduce reference reaction rates andthe ability to preserve reference leakage rates at desired surfaces through Discontinuity Factors (DF).The need for this new class of methods arises from the inability of the current state-of-the-art SPH technology to properly reproduce the reactor multiplication factor for problems with significant leakage. This work shows that this defect lies in the use of normalization factors in the SPH algorithm: while they are necessary with purely reflecting problems to ensure uniqueness of the solution, they introduce homogenization inconsistencies if at least one vacuum boundary condition is present. Two solutions to this problem are presented in this work: (i) to simply remove the normalization factor or (ii) to introduce additional degrees of freedom in the form of DFs in such a way that the normalization factors can still constrain the problem. While the former clearly offers unrivaled simplicity and can lead to very satisfactory results for leakage-dominated cores or if spatial restriction of the SPH regions is applied – as demonstrated for the High Temperature Test Reactor, the latter is more robust and does not require to scale the fission terms with the reference eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

Space Station Freedom combustion research

Extended operations in microgravity, on board spacecraft like Space Station Freedom, provide both unusual opportunities and unusual challenges for combustion science. On the one hand, eliminating the intrusion of buoyancy provides a valuable new perspective for fundamental studies of combustion phenomena. On the other hand, however, the absence of buoyancy creates new hazards of fires and explosions that must be understood to assure safe manned space activities. These considerations - and the relevance of combustion science to problems of pollutants, energy utilization, waste incineration, power and propulsion systems, and fire and explosion hazards, among others - provide strong motivation for microgravity combustion research. The intrusion of buoyancy is a greater impediment to fundamental combustion studies than to most other areas of science. Combustion intrinsically heats gases with the resulting buoyant motion at normal gravity either preventing or vastly complicating measurements. Perversely, this limitation is most evident for fundamental laboratory experiments; few practical combustion phenomena are significantly affected by buoyancy. Thus, we have never observed the most fundamental combustion phenomena - laminar premixed and diffusion flames, heterogeneous flames of particles and surfaces, low-speed turbulent flames, etc. - without substantial buoyant disturbances. This precludes rational merging of theory, where buoyancy is of little interest, and experiments, that always are contaminated by buoyancy, which is the traditional path for developing most areas of science. The current microgravity combustion program seeks to rectify this deficiency using both ground-based and space-based facilities, with experiments involving space-based facilities including: laminar premixed flames, soot processes in laminar jet diffusion flames, structure of laminar and turbulent jet diffusion flames, solid surface combustion, one-dimensional smoldering, ignition and flame spread of liquids, drop combustion, and quenching of panicle-air flames. Unfortunately, the same features that make microgravity attractive for fundamental combustion experiments, introduce new fire and explosion hazards that have no counterpart on earth. For example, microgravity can cause broader flammability limits, novel regimes of flame spread, enhanced effects of flame radiation, slower fire detector response, and enhanced combustion upon injecting fire extinguishing agents, among others. On the other hand, spacecraft provide an opportunity to use 'fire-safe' atmospheres due to their controlled environment. Investigation of these problems is just beginning, with specific fire safety experiments supplementing the space based fundamental experiments listed earlier; thus, much remains to be done to develop an adequate technology base for fire and explosion safety considerations for spacecraft.

Faeth, G. M.↗
DOE OSTI2014

Bridging the Transition from Mesoscale to Microscale Turbulence in Numerical Weather Prediction Models

With a focus towards developing multiscale capabilities in numerical weather prediction models, the specific problem of the transition from the mesoscale to the microscale is investigated. For that purpose, idealized one-way nested mesoscale to large-eddy simulation (LES) experiments were carried out using the Weather Research and Forecasting model framework. It is demonstrated that switching from one-dimensional turbulent diffusion in the mesoscale model to three-dimensional LES mixing does not necessarily result in an instantaneous development of turbulence in the LES domain. However, very large fetches are needed for the natural transition to turbulence to occur. The computational burden imposed by these long fetches necessitates the development of methods to accelerate the generation of turbulence on a nested LES domain forced by a smooth mesoscale inflow. To that end, four new methods based upon finite amplitude perturbations of the potential temperature field along the LES inflow boundaries are developed, and investigated under convective conditions. Each method accelerated the development of turbulence within the LES domain, with two of the methods resulting in a rapid generation of production and inertial range energy content associated to microscales that is consistent with non-nested simulations using periodic boundary conditions. The cell perturbation approach, the simplest and most efficient of the best performing methods, was investigated further under neutral and stable conditions. Successful results were obtained in all the regimes, where satisfactory agreement of mean velocity, variances and turbulent fluxes, as well as velocity and temperature spectra, was achieved with reference non-nested simulations. In contrast, the non-perturbed LES solution exhibited important energy deficits associated to a delayed establishment of fully-developed turbulence. The cell perturbation method has negligible computational cost, significantly accelerates the generation of realistic turbulence, and requires minimal parameter tuning, with the necessary information relatable to mean inflow conditions provided by the mesoscale solution.

17 WIND ENERGY↗
DOE OSTI2018

Residence time distribution in a structured packing unit for monitoring aerosol emissions

Solvent absorption in a structured packed column has been demonstrated as a viable option for post-combustion carbon capture. Amine emissions in the form of aerosol particles pose a major challenge that severely impacts the separation efficiency of the packed column. The aerosol particles quickly achieve the motion of the gas flow. In this context, residence time distribution (RTD) of the gaseous phase will provide insights for managing aerosol particles, thereby controlling amine loss in the packed column. Computational fluid dynamics (CFD) modeling of such columns is a multiscale problem because of a wide range of disparity in length scales. Accordingly, three-dimensional CFD simulations in the representative elementary unit of a structured packing are conducted using OpenFOAM. The RTD is analyzed in terms of residence time distribution function, mean residence time, and variance. Effects of the diffusivity and flow rates on RTD are extensively studied, and the dispersion coefficient is calculated using a one-dimensional axial dispersion model. Furthermore, RTDs computed in a smaller domain are extrapolated to a larger domain using a convolution integral. The computed and corresponding CFD models derived from the convolution integral match well at different flow rates and domain sizes. A comparison of RTD among different packing designs also is presented at a fixed number of repeating units and different flow rates.

structured packings, residence time distribution↗
DOE OSTI2025

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

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

Plasma Science and Innovation Center at Washington, Wisconsin, and Utah State: Final Scientific Report for the University of Wisconsin-Madison

The University of Wisconsin-Madison component of the Plasma Science and Innovation Center (PSI Center) contributed to modeling capabilities and algorithmic efficiency of the Non-Ideal Magnetohydrodynamics with Rotation (NIMROD) Code, which is widely used to model macroscopic dynamics of magnetically confined plasma. It also contributed to the understanding of direct-current (DC) injection of electrical current for initiating and sustaining plasma in three spherical torus experiments: the Helicity Injected Torus-II (HIT-II), the Pegasus Toroidal Experiment, and the National Spherical Torus Experiment (NSTX). The effort was funded through the PSI Center's cooperative agreement with the University of Washington and Utah State University over the period of March 1, 2005 - August 31, 2016. In addition to the computational and physics accomplishments, the Wisconsin effort contributed to the professional education of four graduate students and two postdoctoral research associates. The modeling for HIT-II and Pegasus was directly supported by the cooperative agreement, and contributions to the NSTX modeling were in support of work by Dr. Bickford Hooper, who was funded through a separate grant. Our primary contribution to model development is the implementation of detailed closure relations for collisional plasma. Postdoctoral associate Adam Bayliss implemented the temperature-dependent effects of Braginskii's parallel collisional ion viscosity. As a graduate student, John O'Bryan added runtime options for Braginskii's models and Ji's K2 models of thermal conduction with magnetization effects and thermal equilibration. As a postdoctoral associate, O'Bryan added the magnetization effects for ion viscosity. Another area of model development completed through the PSI-Center is the implementation of Chodura's phenomenological resistivity model. Finally, we investigated and tested linear electron parallel viscosity, leveraged by support from the Center for Extended Magnetohydrodynamic Modeling (CEMM). Work on algorithmic efficiency improved NIMROD's element-based computations. We reordered arrays and eliminated a level of looping for computations over the data points that are used for numerical integration over elements. Moreover, the reordering allows fewer and larger communication calls when using distributed-memory parallel computation, thereby avoiding a data starvation problem that limited parallel scaling over NIMROD's Fourier components for the periodic coordinate. Together with improved parallel preconditioning, work that was supported by CEMM, these developments allowed NIMROD's first scaling to over 10,000 processor cores. Another algorithm improvement supported by the PSI Center is nonlinear numerical diffusivities for implicit advection. We also developed the Stitch code to enhance the flexibility of NIMROD's preprocessing. Our simulations of HIT-II considered conditions with and without fluctuation-induced amplification of poloidal flux, but our validation efforts focused on conditions without amplification. A significant finding is that NIMROD reproduces the dependence of net plasma current as the imposed poloidal flux is varied. The modeling of Pegasus startup from localized DC injectors predicted that development of a tokamak-like configuration occurs through a sequence of current-filament merger events. Comparison of experimentally measured and numerically computed cross-power spectra enhance confidence in NIMROD's simulation of magnetic fluctuations; however, energy confinement remains an open area for further research. Our contributions to the NSTX study include adaptation of the helicity-injection boundary conditions from the HIT-II simulations and support for linear analysis and computation of 3D current-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Robust RF Photocathodes (Phase II Final Report)

This program addressed requirements for photocathodes with high quantum efficiency (QE) and long lifetimes. Currently, photocathodes with QEs above 1% exhibit lifetimes measured in a few hours. This is unacceptable for commercial devices using photocathodes as the electron source. This includes accelerators for industrial, medical, defense, and educational applications as well as emerging applications, such as time resolved electron microscopy. Calabazas Creek Research, Inc. (CCR) and the University of Maryland (UMD) attempted to develop reservoir photocathodes that dramatically increased lifetime and increased the ability to recover from contaminating events, such as loss of system vacuum. The Phase I program demonstrated the viability of the reservoir approach for higher QE cathodes. The Phase II program focused on design, fabrication, and delivery. This included increasing understanding of photocathode physics, investigating design options, addressing performance issues, and developing fabrication and packaging techniques leading to commercial products. The initial task identified specific photocathode materials that provide high QE and are compatible with the reservoir approach. The goal was that these photocathodes be manufacturable in a production environment, understanding that the “environment” is highly specialized. The photocathode geometry was modified to address several photogun performance issues related primarily to emittance. A goal was to reduce or eliminate thermal emittance issues arising from the reservoir heater. The program developed fabrication, handling, and transport techniques and facilities to allow commercial photocathode production for certain materials. Not only was complex chemistry involved, but practical issues of transporting and installing the cathodes into user systems were addressed. A critical goal was demonstration of reservoir photocathodes in operating RF guns, where performance could be directly measured. Success depended on clear demonstration that the cathodes could be produced, delivered, installed, and provide the required performance. The principal effort in this program was to develop reservoir cesium antimonide photocathodes. Cesium antimonide cathodes provide a work function in excess of the 1% goal, while the current generation exhibits lifetimes measured in hours. The reservoir in CCRs’ controlled porosity reservoir (CPR) cathodes provides a mechanism for replacing cesium evaporating from the surface by transporting replacement cesium from the reservoir. This is accomplished by gently heating the reservoir to induce cesium diffusion through micron-scale pores in the reservoir cap and photocathode emitter. Experiments at the UMD determined that cesium diffused from the reservoir would only react with antimony on the surface if the temperature exceeded 300C. At this temperature, however, the surface cesium evaporated faster than it could be replenished from the reservoir. Consequently, the QE decreased each time recessiation was attempted. This appears to be a fundamental problem with this approach. The sections below provide information on accomplishments for each task in the program. Once the issue with cesium antimonide cathodes was identified, CCR refocused the program to begin investigation of graphene coated photocathodes in collaboration with SLAC National Accelerator Laboratory (SLAC). The photocathode test chamber developed in the program was transferred to SLAC for use in the graphene research. This is described below in a task added to the program.

43 PARTICLE ACCELERATORS↗
NASA NTRS1982

Transient thermal stress problem for a circumferentially cracked hollow cylinder

The transient thermal stress problem for a hollow elasticity cylinder containing an internal circumferential edge crack is considered. It is assumed that the problem is axisymmetric with regard to the crack geometry and the loading, and that the inertia effects are negligible. The problem is solved for a cylinder which is suddenly cooled from inside. First the transient temperature and stress distributions in an uncracked cylinder are calculated. By using the equal and opposite of this thermal stress as the crack surface traction in the isothermal cylinder the crack problem is then solved and the stress intensity factor is calculated. The numerical results are obtained as a function of the Fourier number tD/b(2) representing the time for various inner-to-outer radius ratios and relative crack depths, where D and b are respectively the coefficient of diffusivity and the outer radius of the cylinder.

Nied, H. F.↗
NASA NTRS1987

Numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory

A review of numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory is presented. The methods reviewed include explicit, implicit, quasi-linearization, time linearization, operator-splitting, random walk and finite-element techniques and methods of lines. Adaptive and nonadaptive procedures are also reviewed. These techniques are applied first to solve two model problems which have exact traveling wave solutions with which the numerical results can be compared. This comparison is performed in terms of both the wave profile and computed wave speed. It is shown that the computed wave speed is not a good indicator of the accuracy of a particular method. A fourth-order time-linearized, Hermitian compact operator technique is found to be the most accurate method for a variety of time and space sizes.

Ramos, J. I.↗
NASA NTRS1990

A parallel algorithm for nonlinear convection-diffusion equations

A parallel algorithm for the efficient solution of nonlinear time-dependent convection-diffusion equations with small parameter on the diffusion term is presented. The method is based on a physically motivated domain decomposition that is dictated by singular perturbation analysis. The analysis is used to determine regions where certain reduced equations may be solved in place of the full equation. The method is suitable for the solution of problems arising in the simulation of fluid dynamics. Experimental results for a nonlinear equation in two-dimensions are presented.

Scroggs, Jeffrey S.↗
NASA NTRS1991

A computational model for the analysis of diffusive binary alloy solid/liquid phase change energy transport

A model is presented for the prediction of solid/liquid, binary alloy, phase change energy transport. The model incorporates only one-dimensional, diffusive transport of both energy and species but is not inherently limited in any way to such a restriction. The model incorporates a previously developed algorithm for single constituent phase change energy transport and is applied to the energy transport as well as the species transport aspects of the problem. The model is applied to four example cooling problems and is shown to perform exceedingly well from an algorithmic point of view.

Schneider, G. E.↗
NASA NTRS1992

Parabolized reduced Navier-Stokes computational techniques

A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Rubin, Stanley G.↗