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 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Boundary estimation problems arising in thermal tomography

Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in a output least squares approach. Parameter estimation techniques based on the method of mappings are discussed and approximation schemes are developed based on a finite element Galerkin approach. Theoretical convergence results for computational techniques are given and the results are applied to experimental data for the identification of flaws in the thermal testing of materials.

Banks, H. T.↗

NASA NTRS1997

The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems

In this paper, we study the Local Discontinuous Galerkin methods for nonlinear, time-dependent convection-diffusion systems. These methods are an extension of the Runge-Kutta Discontinuous Galerkin methods for purely hyperbolic systems to convection-diffusion systems and share with those methods their high parallelizability, their high-order formal accuracy, and their easy handling of complicated geometries, for convection dominated problems. It is proven that for scalar equations, the Local Discontinuous Galerkin methods are L(sup 2)-stable in the nonlinear case. Moreover, in the linear case, it is shown that if polynomials of degree k are used, the methods are k-th order accurate for general triangulations; although this order of convergence is suboptimal, it is sharp for the LDG methods. Preliminary numerical examples displaying the performance of the method are shown.

Cockburn, Bernardo↗

NASA NTRS2010

Meridional Distribution of CH3C2H and C4H2 in Saturn's Stratosphere from CIRS/Cassini Limb and Nadir Observations

Limb and nadir spectra acquired by Cassini/CIRS (Composite InfraRed Spectrometer) are analyzed in order to derive, for the first time, the meridional variations of diacetylene (C4H2) and methylacetylene (CH3C2H) mixing ratios in Saturn's stratosphere, from 5 hPa up to 0.05 hPa and 80 deg S to 45 deg N. We find that the C4H2 and CH3C2H meridional distributions mimic that of acetylene (C2H2), exhibiting small-scale variations that are not present in photochemical model predictions. The most striking feature of the meridional distribution of both molecules is an asymmetry between mid-southern and mid-northern latitudes. The mid-southern latitudes are found depleted in hydrocarbons relative to their northern counterparts. In contrast, photochemical models predict similar abundances at north and south mid-latitudes. We favor a dynamical explanation for this asymmetry, with upwelling in the south and downwelling in the north, the latter coinciding with the region undergoing ring shadowing. The depletion in hydrocarbons at mid-southern latitudes could also result from chemical reactions with oxygen-bearing molecules. Poleward of 60 deg S, at 0.1 and 0.05 hPa, we find that the CH3C2H and C4H2 abundances increase dramatically. This behavior is in sharp contradiction with photochemical model predictions, which exhibit a strong decrease towards the south pole. Several processes could explain our observations, such as subsidence, a large vertical eddy diffusion coefficient at high altitudes, auroral chemistry that enhances CH3C2H and C4H2 production, or shielding from photolysis by aerosols or molecules produced from auroral chemistry. However, problems remain with all these hypotheses, including the lack of similar behavior at lower altitudes. Our derived mean mixing ratios at 0.5 hPa of (2.4 +/- 0.3) 10(exp -10) for C4H2 and of (1.1 +/- 0.3) 10(exp -9) for CH3C2H are compatible with the analysis of global-average ISO observations performed by Moses et al. Finally, we provide values for the ratios [CH3C2H]/[C2H2] and [C4H2]/[C2H2] that can constrain the coupled chemistry of these hydrocarbons.

Guerlet, Sandrine↗

NASA NTRS1960

Comparison of Turbulent Heat-Transfer Results for Uniform Wall Heat Flux and Uniform Wall Temperature

The purpose of this note is to examine in a more precise way how the Nusselt numbers for turbulent heat transfer in both the fully developed and thermal entrance regions of a circular tube are affected by two different wall boundary conditions. The comparisons are made for: (a) Uniform wall temperature (UWT); and (b) uniform wall heat flux (UHF). Several papers which have been concerned with the turbulent thermal entrance region problem are given. 1 Although these analyses have all utilized an eigenvalue formulation for the thermal entrance region there were differences in the choices of eddy diffusivity expressions, velocity distributions, and methods for carrying out the numerical solutions. These differences were also found in the fully developed analyses. Hence when making a comparison of the analytical results for uniform wall temperature and uniform wall heat flux, it was not known if differences in the Nusselt numbers could be wholly attributed to the difference in wall boundary conditions, since all the analytical results were not obtained in a consistent way. To have results which could be directly compared, computations were carried out for the uniform wall temperature case, using the same eddy diffusivity, velocity distribution, and digital computer program employed for uniform wall heat flux. In addition, the previous work was extended to a lower Reynolds number range so that comparisons could be made over a wide range of both Reynolds and Prandtl numbers.

Siegel, R.↗

NASA NTRSDate not supplied

The Calculation of Diffusion Properties for Open Shell Molecules

- Navier-Stokes requires as input various collision integrals for gases - Diffusion, - Viscosity - Thermal conductivity - Thermal diffusion ... - Levin, Stallcop and Partridge state-of-the-art - Non-Relativistic Born-Oppenheimer Approximation (freeze nuclei first, let them move second) - Solve electron problem: Compute realistic Potential Curves (PC) from first principles, - Solve nuclei problem: Compute scattering properties for each PC using sophisticated semi-classical method - Compute cross sections - Compute appropriate thermal collision integrals

STMD↗

DOE OSTI2016

A Variational Nodal Approach to 2D/1D Pin Resolved Neutron Transport: I Diffusion Theory

The variational nodal method (VNM) as implemented in the VARIANT code has been widely employed to fast reactor cores, but much less so for thermal reactors. In what follows we report the initial phase of development of the PANX code (Purdue - Argonne - Northwestern - Xi'an) for three-dimensional pin resolved PWR transport calculations. In two dimensions, prototype modifications of VARIANT using finite elements to obtain heterogeneous nodes and quasi-reflected interface conditions to reduce the high angular-order interface conditions, have been successfully employed to the 2D C5G7 benchmark. In three dimensions, however, the polynomial interface conditions employed in VARIANT smear fuel and coolant at each axial interface, leading to lack of convergence as the axial mesh is refined. This problem is also encountered in the 2D/1D approximations in MOC based methods. The intent of this work is to eliminate these interface homogenization effects from VNM calculations. For clarity, in this work we examine these phenomena in diffusion theory. The next stage of development will incorporate spherical harmonic trial functions to achieve more accurate transport solutions. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

The role of native defects and their diffusivity in CdTe

CdTe is a leading thin-film technology for solar-cell applications. Basic features such as the optimum band gap of 1.5 eV and simple manufacturing process have been driving the research on improving the current record efficiency of 22% for CdTe solar cells toward the theoretical limit of 32%. Low p-type doping and short carrier lifetimes in the CdTe absorber have been suggested as the main limiting factors and attributed to defect compensation and carrier recombination. Understanding the role of point defects and impurities in CdTe is crucial to solving these problems. Despite many years of research, the sources of compensation and the microscopic recombination mechanisms are still unclear. By using hybrid density functional calculations with spin-orbit coupling and large supercells, we investigate the electronic properties, the formation energies, and the diffusivities of the native defects in CdTe. As possible sources of compensation for p-type conductivity, we discuss the structural and electronic properties of Cd interstitial, Te vacancies, and TeCd antisites. The stability of each defect is discussed in terms of both the defect formation energy and the migration barrier. We also address possible mechanisms that lead to non-radiative recombination.

14 SOLAR ENERGY↗

NASA NTRS2010

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗

NASA NTRS1993

Effects of Buoyancy on Laminar, Transitional, and Turbulent Gas Jet Diffusion Flames

Gas jet diffusion flames have been a subject of research for many years. However, a better understanding of the physical and chemical phenomena occurring in these flames is still needed, and, while the effects of gravity on the burning process have been observed, the basic mechanisms responsible for these changes have yet to be determined. The fundamental mechanisms that control the combustion process are in general coupled and quite complicated. These include mixing, radiation, kinetics, soot formation and disposition, inertia, diffusion, and viscous effects. In order to understand the mechanisms controlling a fire, laboratory-scale laminar and turbulent gas-jet diffusion flames have been extensively studied, which have provided important information in relation to the physico-chemical processes occurring in flames. However, turbulent flames are not fully understood and their understanding requires more fundamental studies of laminar diffusion flames in which the interplay of transport phenomena and chemical kinetics is more tractable. But even this basic, relatively simple flame is not completely characterized in relation to soot formation, radiation, diffusion, and kinetics. Therefore, gaining an understanding of laminar flames is essential to the understanding of turbulent flames, and particularly fires, in which the same basic phenomena occur. In order to improve and verify the theoretical models essential to the interpretation of data, the complexity and degree of coupling of the controlling mechanisms must be reduced. If gravity is isolated, the complication of buoyancy-induced convection would be removed from the problem. In addition, buoyant convection in normal gravity masks the effects of other controlling parameters on the flame. Therefore, the combination of normal-gravity and microgravity data would provide the information, both theoretical and experimental, to improve our understanding of diffusion flames in general, and the effects of gravity on the burning process in particular.

Bahadori, M. Yousef↗

DOE OSTI2018

Lumping Techniques for Mixed Finite Element Diffusion Discretizations

We are interested in positivity preserving spatial discretizations applicable to arbitrary order mixed finite element spatial discretizations of the radiative diffusion equations. Though not guaranteed to yield strict positivity, matrix lumping is appealing as a method of increasing solution positivity because unlike ad hoc fixups or strictly non-negative, non-linear solution representations, matrix lumping keeps the mixed finite element diffusion equations as a linear system of equations. Self-lumping schemes have previously been shown to be a generalization of more traditional matrix lumping techniques. Self-lumping schemes restrict integration quadrature points to the interpolation points used to define interpolatory finite element basis functions. As shown previously, the positivity and order of accuracy of self-lumping schemes are directly related to the choice of finite element interpolation point. Here, we develop self-lumping schemes for arbitrary order, mixed finite element spatial discretizations of the diffusion equation in Cartesian geometry. After introducing several possible self-lumping schemes appropriate for mixed finite elements, we test the ability of each scheme to increase solution positivity relative to an unlumped discretization using a simple test problem designed to induce negative solutions, evaluating efficacy as a function of finite element order and finite element interpolation point type. The method of manufactured solutions is then used to evaluate order of convergence for promising self-lumping schemes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2017

Implementation of CMFD Acceleration Scheme in PROTEUS-MOC

The PROTEUS-MOC code is a three-dimensional (3D) transport code, which is being developed at Argonne National Laboratory under the DOE-NE's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. It provides a fully consistent 3D discretization of the transport equation by combining the two-dimensional (2D) method of characteristics (MOC) for the planar direction with the discontinuous Galerkin method for the treatment of the axial variable. PROTEUS-MOC uses the unstructured finite element method (FEM) and thus allows the explicit representation of any extruded geometries. It also has a capability of massively parallel computation in energy group, angle, plane, and space using a top-level GMRES linear algebra solver. For parallel computation in space, it uses an unstructured domain decomposition by adapting METIS. The computational burden for direct heterogeneous core calculations based on unstructured geometries is excessive for practical reactor problems, and hence it is essential to implement an efficient acceleration scheme in PROTEUS-MOC. Recently, the coarse mesh finite difference (CMFD) method, which was originally applied to the nodal diffusion methods, has been successfully applied to the MOC transport solvers, resulting in a remarkable reduction in the number of fission source iterations. In addition, a modified CMFD method has been developed to overcome the limitation of the conventional CMFD that is applicable only to neutronics solvers with a conservative discretization. This CMFD method was derived to be consistent with the FEM, which is based on a nonconservative discretization. In this method, the neutron balance in each coarse mesh is enforced by introducing a pseudo absorption cross section (PAXS), and the well-established alternating solution process of CMFD and transport calculations is employed to accelerate source convergence. This modified CMFD acceleration scheme was adopted in PROTEUS-MOC to accelerate the fission source iteration. Pin-cell and assembly based CMFD schemes have been implemented by defining individual pin cells or assemblies as a coarse mesh. For this, an algorithm to map the finite elements to coarse meshes was developed, and the routines were developed to determine the coupling coefficients and PAXS of the CMFD equations from the transport solution. Preliminary performance tests were carried out using the C5G7 benchmark problem. This paper presents the main features of this CMFD acceleration scheme and the preliminary test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Multi-scale topology optimization of multi-material structures with controllable geometric complexity — Applications to heat transfer problems

This study presents a topology optimization method to design assemblies of periodic cellular materials with controllable geometric complexity. The framework is based on a novel multi-scale and multi-material design model in which the structure, the layout of the material subdomains, and their micro-structures are optimized concurrently. To allow a tight control over its geometrical complexity, the layout at the macro-scale is described by level-set fields, parameterized by geometric primitives. The micro-scale geometry is represented through a density approach. A nonlinear programming algorithm drives the optimization process using design sensitivities computed by the discrete adjoint method. The proposed design framework is studied with heat transfer problems. Practical design problems, such as a heat sink and a thermal storage unit with phase change are discussed. The macro-scale analysis model relies on a generalized transient diffusion equation. At the micro-scale, homogenization is used to compute equivalent material properties. The numerical examples in this study show that the optimized multi-scale multi-material layouts outperform the corresponding mono-scale structures while maintaining geometric simplicity at the macro-scale level.

42 ENGINEERING↗

DOE OSTI2015

Theory of activated penetrant diffusion in viscous fluids and colloidal suspensions

We heuristically formulate a microscopic, force level, self-consistent nonlinear Langevin equation theory for activated barrier hopping and non-hydrodynamic diffusion of a hard sphere penetrant in very dense hard sphere fluid matrices. Penetrant dynamics is controlled by a rich competition between force relaxation due to penetrant self-motion and collective matrix structural (alpha) relaxation. In the absence of penetrant-matrix attraction, three activated dynamical regimes are predicted as a function of penetrant-matrix size ratio which are physically distinguished by penetrant jump distance and the nature of matrix motion required to facilitate its hopping. The penetrant diffusion constant decreases the fastest with size ratio for relatively small penetrants where the matrix effectively acts as a vibrating amorphous solid. Increasing penetrant-matrix attraction strength reduces penetrant diffusivity due to physical bonding. For size ratios approaching unity, a distinct dynamical regime emerges associated with strong slaving of penetrant hopping to matrix structural relaxation. A crossover regime at intermediate penetrant-matrix size ratio connects the two limiting behaviors for hard penetrants, but essentially disappears if there are strong attractions with the matrix. Activated penetrant diffusivity decreases strongly with matrix volume fraction in a manner that intensifies as the size ratio increases. We propose and implement a quasi-universal approach for activated diffusion of a rigid atomic/molecular penetrant in a supercooled liquid based on a mapping between the hard sphere system and thermal liquids. Calculations for specific systems agree reasonably well with experiments over a wide range of temperature, covering more than 10 orders of magnitude of variation of the penetrant diffusion constant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Preferential heating of oxygen 5{sup +} ions by finite-amplitude oblique Alfvén waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5{sup +} ions by large-scale finite-amplitude Alfvén waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5{sup +} ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfvén-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles θ ≤ 30°. The obliquely propagating Alfvén pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

43 PARTICLE ACCELERATORS↗

NASA NTRS1973

Analysis of gas flow through a multilayer insulation system.

A theoretical investigation of gas flow inside a multilayer insulation system has been made for the case of the broadside pumping process. A set of simultaneous first-order differential equations for the temperature and pressure of the gas mixture was obtained by considering the diffusion mechanism of the gas molecules through the perforations on the insulation layers. A modified Runge-Kutta method was used for numerical experiment. The numerical stability problem was investigated. It has been shown that when the relaxation time is small compared with the time period over which the gas properties change appreciably, the set of differential equations can be replaced by a set of algebraic equations for solution. Numerical examples were given and comparison with experimental data was made.

Lin, J. T.↗

NASA NTRS1973

Thermal engineering research

Shih, C. C.↗

NASA NTRS1976

Design, development, and test of a laser velocimeter for high speed turbomachinery

The paper describes the design, development, and test of a laser Doppler velocimeter (LDV) to map the velocity fields in a small, high speed, 8:1 pressure ratio, 0.91 kg/s (2 lbm/s), centrifugal compressor. This instrument is being used as a diagnostic tool for the study of the basic fluid dynamics of the inducer, impeller and the diffuser regions of high pressure ratio, small, centrifugal compressors. The LDV instrumentation has been optimized to permit the measurement of instantaneous velocities up to approximately 500 m/s, measured in absolute coordinates, within a rotating compressor impeller. Velocities of the same magnitude can also be measured in the two-dimensional radial plane of the diffuser. It is shown that LDV can make nondisturbing gas velocity measurements and surmount the difficulties found in applying other types of aerodynamic instrumentation to this measurement problem.

Runstadler, P. W., Jr.↗