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

Light Water Reactor Sustainability Program Modeling of Cu Precipitate Contributions to Reactor Pressure Vessel Steel Microstructure Evolution and Embrittlement

This report summarizes the lower-length-scale effort during FY 2017 in developing simulation capabilities for the mesoscale: models for microstructure evolution and plasticity in reactor pressure vessel steels. During operation, reactor pressure vessels undergo hardening and embrittlement caused by irradiation induced defect accumulation and irradiation enhanced solute precipitation. Both defect production and solute precipitation initiate on the atomic scale and manifest their eventual effects as degradation in the engineering scale properties. To predict the engineering scale property degradation, multiscale modeling and simulation capabilities are needed to better understand the microstructure evolution and to connect the microstructure feature evolution to the engineering scale material properties. In this report the development of coupling across different length scale models to understand the multi-physics property dedgration problem is discussed. The development of coupling a lattice kinetic Monte Carlo code to cluster dynamics to capture irradiation enhanced diffusion in the formation of defects is summarized. The implementation of a cluster dynamics model with Grizzly to predict Cu percipitate formation is discussed. A summary of the coupling of crystal plasticity and cluster dynamics to predict the effect on reactor pressure steel hardening is given.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Cosmic-Ray Diffusion Coefficients throughout the Inner Heliosphere from a Global Solar Wind Simulation

We use a three-dimensional magnetohydrodynamic simulation of the solar wind to calculate cosmic-ray diffusion coefficients throughout the inner heliosphere (2 R{sub ⊙}−3 au). The simulation resolves large-scale solar wind flow, which is coupled to small-scale fluctuations through a turbulence model. Simulation results specify background solar wind fields and turbulence parameters, which are used to compute diffusion coefficients and study their behavior in the inner heliosphere. The parallel mean free path (mfp) is evaluated using quasi-linear theory, while the perpendicular mfp is determined from nonlinear guiding center theory with the random ballistic interpretation. Several runs examine varying turbulent energy and different solar source dipole tilts. We find that for most of the inner heliosphere, the radial mfp is dominated by diffusion parallel to the mean magnetic field; the parallel mfp remains at least an order of magnitude larger than the perpendicular mfp, except in the heliospheric current sheet, where the perpendicular mfp may be a few times larger than the parallel mfp. In the ecliptic region, the perpendicular mfp may influence the radial mfp at heliocentric distances larger than 1.5 au; our estimations of the parallel mfp in the ecliptic region at 1 au agree well with the Palmer “consensus” range of 0.08–0.3 au. Solar activity increases perpendicular diffusion and reduces parallel diffusion. The parallel mfp mostly varies with rigidity (P) as P{sup .33}, and the perpendicular mfp is weakly dependent on P. The mfps are weakly influenced by the choice of long-wavelength power spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Vertical distribution of vibrational energy of molecular nitrogen in a stable auroral red arc and its effect on ionospheric electron densities

Previous solutions of the problem of the distribution of vibrationally excited molecular nitrogen in the thermosphere have either assumed a Boltzmann distribution and considered diffusion as one of the loss processes or solved for the energy level populations and neglected diffusion. Both of the previous approaches are combined by solving the time dependent continuity equations, including the diffusion process, for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to a stable auroral red arc. The primary source of molecular nitrogen excitation was subexcitation, and inelastic collisions between thermal electrons and molecular nitrogen. The reaction rates for this process were calculated from published cross section calculations. The loss processes for vibrational energy were electron and atomic oxygen quenching and vibrational energy exchange. The coupled sets of nonlinear, partial differential equations were solved numerically by employing finite difference equations.

Newton, G. P.↗
NASA NTRS1977

Further investigation of the spontaneous and evoked activity of the primary neurons of statoreceptors (and other receptors) of the labyrinth of the bullfrog before, during and after an extended period of weightlessness, including alternative intervals of artificial gravity

Vestibular neuron activity was examined by studying nerve stimulation and evoked response. A cooling element, applied to the nerve consisted of a silver hook through which a coolant fluid flowed. Temperature changes were recorded via microtermistors on an eight channel brush recorder, together with response. Diffusion of the cooling effect was measured, recovery time was assessed, and the nerve was then studied hystologically and ultrastructurally. Problems in frog preparation were discussed along with problems in maintaining healthy specimens and bacteria controlled aquaria.

Source record↗
NASA NTRS1982

The use of optimization techniques to design controlled diffusion compressor blading

A method for automating compressor blade design using numerical optimization, and applied to the design of a controlled diffusion stator blade row is presented. A general purpose optimization procedure is employed, based on conjugate directions for locally unconstrained problems and on feasible directions for locally constrained problems. Coupled to the optimizer is an analysis package consisting of three analysis programs which calculate blade geometry, inviscid flow, and blade surface boundary layers. The optimizing concepts and selection of design objective and constraints are described. The procedure for automating the design of a two dimensional blade section is discussed, and design results are presented.

Sanger, N. L.↗
NASA NTRS1981

Silicon solar cell process development, fabrication, and analysis

Work has progressed in fabrication and characterization of solar cells from ubiquitous crystallization process (UCP) wafers and LASS ribbons. Gettering tests applied to UCP wafers made little change on their performance compared with corresponding baseline data. Advanced processes such as shallow junction (SJ), back surface field (BSF), and multilayer antireflection (MLAR) were also applied. While BSF by Al paste had shunting problems, cells with SJ and BSF by evaporated Al, and MLAR did achieve 14.1% AMI on UCP silicon. The study of LASS material was very preliminary. Only a few cells with SJ, BSR, (no BSF) and MLAR were completed due to mechanical yield problems after lapping the material. Average efficiency was 10.7% AMI with 13.4% AMI for CZ controls. Relatively high minority carrier diffusion lengths were obtained. The lower than expected Jsc could be partially explained by low active area due to irregular sizes.

Yoo, H. I.↗
NASA NTRS1981

Monte-Carlo model of pitch-angle scattering in solar cosmic ray events

Monte Carlo simulations of the propagation of solar cosmic-rays in interplanetary space are reported, including the effects of pitch-angle scattering and adiabatic focusing. Time intensity profiles agree well with the corresponding spatial diffusion approximation for models close to those used by Gombosi and Owens (1980). Monte Carlo simulations of the same model problem yield results that disagree with those of Gombosi and Owens, and the most likely possible source of error in the Monte Carlo simulation is considered to be in the region near the 90 deg pitch angle where the finite step size introduces certain inconsistencies.

Palmer, I. D.↗
NASA NTRS1987

Particle sizing in rocket motor studies utilizing hologram image processing

A technique of obtaining particle size information from holograms of combustion products is described. The holograms are obtained with a pulsed ruby laser through windows in a combustion chamber. The reconstruction is done with a krypton laser with the real image being viewed through a microscope. The particle size information is measured with a Quantimet 720 image processing system which can discriminate various features and perform measurements of the portions of interest in the image. Various problems that arise in the technique are discussed, especially those that are a consequence of the speckle due to the diffuse illumination used in the recording process.

Netzer, David↗
NASA NTRS1987

Development of a unique laboratory standard: Indium gallium arsenide detector for the 500-1700 nm spectral region

A planar (5 mm diameter) indium gallium arsenide detector having a high (greater than 50 pct) quantum efficiency from the visible into the infrared spectrum (500 to 1700 nm) was fabricated. Quantum efficiencies as high as 37 pct at 510 nm, 58 pct at 820 nm and 62 pct at 1300 nm and 1550 nm were measured. A planar InP/InGaAs detector structure was also fabricated using vapor phase epitaxy to grow device structures with 0, 0.2, 0.4 and 0.6 micrometer thick InP caps. Quantum efficiency was studied as a function of cap thickness. Conventional detector structures were also used by completely etching off the InP cap after zinc diffusion. Calibrated quantum efficiencies were measured. Best results were obtained with devices whose caps were completely removed by etching. Certain problems still remain with these detectors including non-uniform shunt resistance, reproducibility, contact resistance and narrow band anti-reflection coatings.

Source record↗
NASA NTRS1992

Model capabilities, 2-D

Two dimensional (2-D) atmospheric models provide results for altitude versus latitude as a function of time and are developed primarily for two reasons: to help understand atmospheric occurrences and to give assessments and/or make predictions of future changes in the atmosphere. Historically, the formulation of transport in 2-D models has been a difficult problem. Most current 2-D models have a transport that is either an Eulerian mean circulation with large stratospheric eddy diffusion or a residual (diabatic or Lagragian) mean circulation which typically is accompanied with small stratospheric eddy diffusion. Because of the assumption of zonal averaging, 2-D models are primarily useful in making predictions of atmospheric changes of time scales longer than a season. Although decadel atmospheric changes may be reasonably well represented with a 2-D model, the year to year changes which result from interannual transport differences, stratospheric warmings, semiannual oscillations, or quasi-biennial oscillations may not be well represented in the stratosphere and troposphere.

Jackman, Charles H.↗
NASA NTRS1995

Effects of confinement on partially premixed flames

Partially premixed combustion is an intermediate regime between the limiting cases of premixed and nonpremixed combustion. Although combustion problems are generally approached from one of these two limiting cases, there are many practical situations where flames cannot be considered as purely premixed or nonpremixed, and thus the partially premixed approach must be used. In partially premixed combustion, mechanisms from the premixed and nonpremixed regimes can coexist, and as a result some interesting new phenomena can arise. One such phenomenon is the flame stabilization in laminar mixing layers by triple flames. One of the first observations of triple flames was made by Phillips (1965), who investigated a triple flame propagating in a methane mixing layer. Kioni et al. (1993) also examined triple flames both experimentally and numerically. There have also been numerous analytical studies on the shape and propagation of triple flames under various assumptions by Dold (1989), Dold et al. (1991), and Hartley and Dold (1991). In terms of modeling, Muller et al. (1994) have combined the flamelet formulations for premixed and nonpremixed combustion in order to treat lifted diffusion flames. One common feature in the analytical and numerical studies mentioned above is the assumption of zero heat release, which is necessary to make the problem tractable. The effect of heat release on triple flames was investigated by Ruetsch et al. (1995), where for the unconfined case, flame speeds larger than their premixed counterparts were found. One of the most important practical situations in which these conditions arise is in lifted turbulent jet diffusion flames. At a critical velocity the burning zone of a fuel jet lifts off from the nozzle, moves to increasing distances as the jet velocity increases, and finally blows off. The mechanisms that control these phenomena, i.e. that determine the stability of these flames, are still not understood. In addition to regions where diffusion flame stabilization takes place, partially premixed conditions also exist during the ignition process in nonpremixed systems. Numerical simulations by Reveillon et al. (1994) of the ignition process in a weakly stirred mixture of fuel and oxidizer show that triple flames propagate along lines of stoichiometric mixture fraction throughout the fluid. In addition, Peters (1994) notes that NO(x) emissions are likely to be large in such transient cases, and therefore an understanding of triple flames can provide information concerning pollutant formation. This study extends the work previously done and examines the effects of lateral confinement on partially premixed flames. Once again, we study both the flame structure and propagation.

Ruetsch, G. R.↗
NASA NTRS2007

27 Years of Satellite Ozone Data: Merging of Data Records from Multiple Instruments to Observe Global Trends and Recovery

Satellite measurements provide a unique global view of the stratospheric ozone layer. The perspective from satellites allowed for the early mapping of the extent of the phenomenon that became known as the ozone hole. The use of the satellite data for global trends outside of the ozone hole confronts the problem of the possible drift of the calibration of the instrument. The TOMS and SBUV instruments on Nimbus 7 lasted for more than a decade. During that time, the diffuser plate used to reflect sunlight into the measurement degraded (darkened) and the instruments each had a number of events that made calibration determination difficult. Initially the TOMS data were used for global trends by adjusting the overall calibration to agree with a set of ground-based measurement stations. But this was unsatisfactory because the record was not independent of those ground measurements and problems were found in many of the ground stations by using TOMS as a transfer standard. After many years of dedicated work, the TOMS/SBUV team learned how to correct for instrument drift, remove the interfering effects of aerosols, and establish instrument-to-instrument calibrations resulting in a long-term record that can be used for accurate trend and recovery determination. The global view of the satellites allows for determination not only of temporal change in ozone, but spatial fingerprints that allow more confidence in assigning cause to observed changes.

Stolarski, Richard S.↗
NASA NTRS1975

On the use of a finite difference method for solving anisotropic scattering problems

A new method of solving the radiative transfer equation is developed in which the scattering and absorption coefficients may have arbitrary variations with depth, and in which both internal (thermal) emission and incident radiation are allowed. Specular and diffuse reflection at both boundaries also is taken into account. The method begins by forming a paired set of coupled first-order differential equations for the symmetric and antisymmetric parts of the radiation field after writing the scattering integral as a numerical quadrature. These differential equations are broken into finite difference form, in which the symmetric and antisymmetric parts of the radiation field are found on alternate grid points. Numerical results for a number of test problems are shown, demonstrating that the method is very fast, that it returns specific intensities and fluxes that are accurate to at least a percent, and that it can be applied to optically thick problems.

Barkstrom, B. R.↗
NASA NTRS1987

Information retrieval from holographic interferograms: Fundamentals and problems

Holographic interferograms can contain large amounts of information about flow and temperature fields. Their information content can be very high because they can be viewed from many different directions. This multidirectionality, and fringe localization add to the information contained in the fringe pattern if diffuse illumination is used. Additional information, and increased accuracy can be obtained through the use of dual reference wave holography to add reference fringes or to effect discrete phase shift or hetrodyne interferometry. Automated analysis of fringes is possible if interferograms are of simple structure and good quality. However, in practice a large number of practical problems can arise, so that a difficult image processing task results.

Vest, Charles M.↗
NASA NTRS1975

Viscous flow solutions with a cubic spline approximation

A cubic spline approximation is used for the solution of several problems in fluid mechanics. This procedure provides a high degree of accuracy even with a nonuniform mesh, and leads to a more accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several typical integration schemes are presented. For two-dimensional flows a spline-alternating-direction-implicit (SADI) method is evaluated. The spline procedure is assessed and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗
NASA NTRS1997

Diffuse-Interface Methods in Fluid Mechanics

The authors review the development of diffuse-interface models of hydrodynamics and their application to a wide variety of interfacial phenomena. The authors discuss the issues involved in formulating diffuse-interface models for single-component and binary fluids. Recent applications and computations using these models are discussed in each case. Further, the authors address issues including sharp-interface analyses that relate these models to the classical free-boundary problem, related computational approaches to describe interfacial phenomena, and related approaches describing fully-miscible fluids.

SURFACE DIFFUSION↗
DOE OSTI2017

An Integrated Approach to Parameter Learning in Infinite-Dimensional Space

The availability of sophisticated modern physics codes has greatly extended the ability of domain scientists to understand the processes underlying their observations of complicated processes, but it has also introduced the curse of dimensionality via the many user-set parameters available to tune. Many of these parameters are naturally expressed as functional data, such as initial temperature distributions, equations of state, and controls. Thus, when attempting to find parameters that match observed data, being able to navigate parameter-space becomes highly non-trivial, especially considering that accurate simulations can be expensive both in terms of time and money. Existing solutions include batch-parallel simulations, high-dimensional, derivative-free optimization, and expert guessing, all of which make some contribution to solving the problem but do not completely resolve the issue. In this work, we explore the possibility of coupling together all three of the techniques just described by designing user-guided, batch-parallel optimization schemes. Our motivating example is a neutron diffusion partial differential equation where the time-varying multiplication factor serves as the unknown control parameter to be learned. We find that a simple, batch-parallelizable, random-walk scheme is able to make some progress on the problem but does not by itself produce satisfactory results. After reducing the dimensionality of the problem using functional principal component analysis (fPCA), we are able to track the progress of the solver in a visually simple way as well as viewing the associated principle components. This allows a human to make reasonable guesses about which points in the state space the random walker should try next. Thus, by combining the random walker's ability to find descent directions with the human's understanding of the underlying physics, it is possible to use expensive simulations more efficiently and more quickly arrive at the desired parameter set.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗