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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Flow Mechanisms and Diffusion Combustion of Turbulent Jets

The problem of flow and combustion of turbulent jets of fuel gas in the external medium of an oxidant (air) is solved with regard to the existence of the actual boundary of the turbulent flow region of a jet. Based on the ideas of the friction force of the external flow acting on the boundary of a jet, the entrainment equation for the external medium is derived that closes the system of equations of motion of turbulent jets. The physical meaning of the dissipation rate of the turbulent energy of a jet is interpreted as the work of the friction force. To describe the combustion kinetics, the limit of instantaneous reactions corresponding to the diffusion combustion mode is used. Calculations of the effective reaction rates for reactants and the volumes occupied by them are based on the representation of a turbulent medium as an aggregation of independent turbulent particles—vortices—whose random contacts lead to the mixing and combustion of reacting substances [31]. The concomitant phenomena of flow and combustion are analyzed, including radiation effects. In particular, it is shown that the apparent increase in the combustion temperature with increasing Reynolds number is in fact attributed to the relative decrease of thermal radiation losses. Qualitative agreement is obtained between the results of the theoretical calculations of the length of a combustion torch and experimental data.

42 ENGINEERING↗

NASA NTRS2013

A Patch Density Recommendation based on Convergence Studies for Vehicle Panel Vibration Response resulting from Excitation by a Diffuse Acoustic Field

Fluid structural interaction problems that estimate panel vibration from an applied pressure field excitation are quite dependent on the spatial correlation of the pressure field. There is a danger of either over estimating a low frequency response or under predicting broad band panel response in the more modally dense bands if the pressure field spatial correlation is not accounted for adequately. Even when the analyst elects to use a fitted function for the spatial correlation an error may be introduced if the choice of patch density is not fine enough to represent the more continuous spatial correlation function throughout the intended frequency range of interest. Both qualitative and quantitative illustrations evaluating the adequacy of different patch density assumptions to approximate the fitted spatial correlation function are provided. The actual response of a typical vehicle panel system is then evaluated in a convergence study where the patch density assumptions are varied over the same finite element model. The convergence study results are presented illustrating the impact resulting from a poor choice of patch density. The fitted correlation function used in this study represents a Diffuse Acoustic Field (DAF) excitation of the panel to produce vibration response.

Smith, Andrew↗

NASA NTRS1975

Study of diffusion coefficients of glasses under zero-G

Diffusion studies of the glass forming ion are examined in zero-g environments and diffusion data obtained from these experiments are unique because of earth based experimental problems. The choice of system for diffusion studies is discussed along with the lab processing. The space processing is described consisting of a heating cycle designed to maximize time exposed to the diffusion temperature without exposing the low viscosity melt to gravitational forces.

Kinser, D. L.↗

NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

NASA NTRS2001

Structure and Early Soot Oxidation Properties of Laminar Diffusion Flames

Soot is an important unsolved problem of combustion science because it is present in most hydrocarbon-fueled flames and current understanding of the reactive and physical properties of soot in flame environments is limited. This lack of understanding affects progress toward developing reliable predictions of flame radiation properties, reliable predictions of flame pollutant emission properties and reliable methods of computational combustion, among others. Motivated by these observations, the present investigation extended past studies of soot formation in this laboratory, to consider soot oxidation in laminar diffusion flames using similar methods. Early work showed that O2 was responsible for soot oxidation in high temperature O2-rich environments. Subsequent work in high temperature flame environments having small O2 concentrations, however, showed that soot oxidation rates substantially exceeded estimates based on the classical O2 oxidation rates of Nagle and Strickland-Constable and suggests that radicals such as O and OH might be strong contributors to soot oxidation for such conditions. Neoh et al. subsequently made observations in premixed flames, supported by later work, that showed that OH was responsible for soot oxidation at these conditions with a very reasonable collision efficiency of 0.13. Subsequent studies in diffusion flames, however, were not in agreement with the premixed flame studies: they agreed that OH played a dominant role in soot oxidation in flames, but found collision efficiencies that varied with flame conditions and were not in good agreement with each other or with Neoh et al. One explanation for these discrepancies is that optical scattering and extinction properties were used to infer soot structure properties for the studies that have not been very successful for representing the optical properties of soot. Whatever the source of the problem, however, these differences among observations of soot oxidation in premixed and diffusion flames clearly must be resolved. Motivated by these findings, the present study undertook measurements of soot and flame properties within the soot oxidation region of some typical laminar diffusion flames and exploited the new measurements to identify soot oxidation mechanisms for these conditions. Present considerations were limited to the early stages of soot oxidation (carbon consumption less than 70%) where reactions at the surface of primary soot particles dominate the process, rather than the later stages when particle porosity and internal particle oxidation become important as discussed by Neoh et al.

El-Leathy, A. M.↗

DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗

NASA NTRS1990

Sunlight Diffusing Tent for Lunar Worksite

The purpose is to provide a solution to problems astronauts encounter with sunlight on the lunar surface. Due to the absence of an atmosphere the Moon is subjected to intense sunlight creating problems with color and contrast. This problem can be overcome by providing a way to reduce intensity and diffuse the light in a working environment. The solution to the problem utilizes an umbrella, tent-like structure covered with a diffusing material. The design takes into account structural materials, stresses, fabrics, and deployment.

Burleson, Blair↗

NASA NTRS1975

Steady-state solution of the semi-empirical diffusion equation for area sources

The problem investigated can be solved exactly in a simple manner if the equations are written in terms of a similarity variable. The exact solution is used to explore two questions of interest in the modelling of urban air pollution, taking into account the distribution of surface concentration downwind of an area source and the distribution of concentration with height.

Lebedeff, S. A.↗

NASA NTRS1987

The origin of the diffuse X-ray background

The Einstein Deep Survey fields are used to reexamine the problem of the origin of the extragalactic diffuse X-ray background. An approach based on fluctuation analysis is adopted, and a sophisticated Monte Carlo scheme is used to compare the X-ray images with a variety of models representing point and diffuse source contributions. Results are presented for two classes of models: those which determine the flux cutoff beyond which pointlike sources fail to contribute significantly if the log N-log S slope is fixed at the measured value of -1.5, and those which set limits on the slope of the log N-log S relation below its last measured point at the Deep Survey threshold. In both cases, it is found that the known population of X-ray-emitting quasars and any reasonable extrapolations therefrom do not comprise the entire X-ray background.

Hamilton, T. T.↗

DOE OSTI2016

Shape and Topology Optimization in Stokes Flow with a Phase Field Approach

In this paper we introduce a new formulation for shape optimization problems in fluids in a diffuse interface setting that can in particular handle topological changes. By adding the Ginzburg–Landau energy as a regularization to the objective functional and relaxing the non-permeability outside the fluid region by introducing a porous medium approach we hence obtain a phase field problem where the existence of a minimizer can be guaranteed. This problem is additionally related to a sharp interface problem, where the permeability of the non-fluid region is zero. In both the sharp and the diffuse interface setting we can derive necessary optimality conditions using only the natural regularity of the minimizers. We also pass to the limit in the first order conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1973

Space simulation ultimate pressure lowered two decades by removal of diffusion pump oil contaminants during operation

The complex problem why large space simulation chambers do not realize the true ultimate vacuum was investigated. Some contaminating factors affecting diffusion pump performance have been identified and some advances in vacuum/distillation/fractionation technology have been achieved which resulted in a two decade or more lower ultimate pressure. Data are presented to show the overall or individual contaminating effect of commonly used phthalate ester plasticizers of 390 to 530 molecular weight on diffusion pump performance. Methods for removing contaminants from diffusion pump silicone oil during operation and reclaiming contaminated oil by high vacuum molecular distillation are described.

Buggele, A. E.↗

DOE OSTI2016

Using Directional Diffusion Coefficients for Nonlinear Diffusion Acceleration of the First Order S{sub N} Equations in Near-Void Regions

The common definition of the diffusion coefficient as the inverse of three times the transport cross section is not compatible with voids. Morel introduced a non-local tensor diffusion coefficient that remains finite in voids. It can be obtained by solving an auxiliary transport problem without scattering or fission. Larsen and Trahan successfully applied this diffusion coefficient for enhancing the accuracy of diffusion solutions of very high temperature reactor (VHTR) problems that feature large, optically thin channels in the z-direction. It was demonstrated that a significant reduction of error can be achieved in particular in the optically thin region. Along the same line of thought, non-local diffusion tensors have been applied to modeling the TREAT reactor confirming the findings of Larsen and Trahan. Previous work of the authors have introduced a flexible Nonlinear-Diffusion Acceleration (NDA) method for the first order S{sub N} equations discretized with the discontinuous finite element method (DFEM). This NDA method uses a scalar diffusion coefficient in the low-order system that is obtained as the flux weighted average of the inverse transport cross section. Hence, it suffers from very large and potentially unbounded diffusion coefficients in the low order problem. However, it was noted that the choice of the diffusion coefficient does not influence consistency of the method at convergence and hence the diffusion coefficient is essentially a free parameter. The choice of the diffusion coefficient does, however, affect the convergence behavior of the nonlinear diffusion iterations. Within this work we use Morel's non-local diffusion coefficient in the aforementioned NDA formulation in lieu of the flux weighted inverse of three times the transport cross section. The goal of this paper is to demonstrate that significant enhancement of the spectral properties of NDA can be achieved in near void regions. For testing the spectral properties of the NDA with non-local diffusion coefficients, the periodical horizontal interface problem is used. This problem consists of alternating stripes of optically thin and thick materials both of which feature scattering ratios close to unity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1978

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗

DOE OSTI2016

Using Directional Diffusion Coefficients for Nonlinear Diffusion Acceleration of the First Order SN Equations in Near-Void Regions

The common definition of the diffusion coeffcient as the inverse of three times the transport cross section is not compat- ible with voids. Morel introduced a non-local tensor diffusion coeffcient that remains finite in voids[1]. It can be obtained by solving an auxiliary transport problem without scattering or fission. Larsen and Trahan successfully applied this diffusion coeffcient for enhancing the accuracy of diffusion solutions of very high temperature reactor (VHTR) problems that feature large, optically thin channels in the z-direction [2]. It is demonstrated that a significant reduction of error can be achieved in particular in the optically thin region. Along the same line of thought, non-local diffusion tensors are applied modeling the TREAT reactor confirming the findings of Larsen and Trahan [3]. Previous work of the authors have introduced a flexible Nonlinear-Diffusion Acceleration (NDA) method for the first order S N equations discretized with the discontinuous finite element method (DFEM), [4], [5], [6]. This NDA method uses a scalar diffusion coeffcient in the low-order system that is obtained as the flux weighted average of the inverse transport cross section. Hence, it su?ers from very large and potentially unbounded diffusion coeffcients in the low order problem. However, it was noted that the choice of the diffusion coeffcient does not influence consistency of the method at convergence and hence the di?usion coeffcient is essentially a free parameter. The choice of the di?usion coeffcient does, however, affect the convergence behavior of the nonlinear di?usion iterations. Within this work we use Morel’s non-local di?usion coef- ficient in the aforementioned NDA formulation in lieu of the flux weighted inverse of three times the transport cross section. The goal of this paper is to demonstrate that significant en- hancement of the spectral properties of NDA can be achieved in near void regions. For testing the spectral properties of the NDA with non-local diffusion coeffcients, the periodical horizontal interface problem is used [7]. This problem consists of alternating stripes of optically thin and thick materials both of which feature scattering ratios close to unity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2015

A robust and accurate numerical method for transcritical turbulent flows at supercritical pressure with an arbitrary equation of state

This paper addresses issues in high-fidelity numerical simulations of transcritical turbulent flows at supercritical pressure. The proposed strategy builds on a tabulated look-up table method based on REFPROP database for an accurate estimation of non-linear behaviors of thermodynamic and fluid transport properties at the transcritical conditions. Based on the look-up table method we propose a numerical method that satisfies high-order spatial accuracy, spurious-oscillation-free property, and capability of capturing the abrupt variation in thermodynamic properties across the transcritical contact surface. The method introduces artificial mass diffusivity to the continuity and momentum equations in a physically-consistent manner in order to capture the steep transcritical thermodynamic variations robustly while maintaining spurious-oscillation-free property in the velocity field. The pressure evolution equation is derived from the full compressible Navier–Stokes equations and solved instead of solving the total energy equation to achieve the spurious pressure oscillation free property with an arbitrary equation of state including the present look-up table method. Flow problems with and without physical diffusion are employed for the numerical tests to validate the robustness, accuracy, and consistency of the proposed approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Fourier analysis of finite element preconditioned collocation schemes

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Deville, Michel O.↗

DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗