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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

High efficiency GaAs solar cells

The present status of the GaAlAs/GaAs heteroface solar cell program is reported. Studies have been concentrated on GaAlAs/GaAs heteroface solar cells; however, some research has been conducted on thin junction, diffused GaAs solar cells. Emphasis has been on obtaining high efficiency (18% to 20%) GaAs solar cells. Two problems that have limited the efficiency of GaAs solar cells are the high recombination velocity of carriers near the surface and the low minority carrier diffusion length in n-GaAs.

Walker, G. H.↗
NASA NTRS1998

Re-examination of a Classic Experiment to Measure the Positronium-Helium Cross Section

In 1975, before the advent of positronium beams, a clever experiment was carried out in an attempt to measure low-energy Ps-atom scattering cross-sections, especially that of Ps-He. A series of parallel thin plates was placed in the vessel containing the sample of gas, and positrons were allowed to enter the chamber from a radioactive source. The time spectrum of annihilation radiation was observed in the traditional way, and the pickoff annihilation due to thermalized Ps striking one of the plates was to be the special feature of the experiment. Using a diffusion analysis the authors derived cross- sections for several gases, and for helium the value obtained was sigma = 0.0166 pi(alpha)(sub o, sup 2). Even then this value was thought to be very much too small, while recent measurements and some calculations suggest a more reasonable value would be sigma = 3 to 13 pi(alpha)(sub o,sup 2). It has always been puzzling that an apparently well-designed experiment should give such an unsatisfactory result. We have now re-examined the analysis of the data with some interesting consequences which may explain the discrepancy. Two main observations support our re-analysis. First, we note that the mean free path for Ps-He scattering would be quite long if the cross-section were as small as that quoted above; the diffusion method itself would seem to be questionable. For the larger values, however, there would be no such problem. Second, it was assumed that when the annihilation rate had settled down and was following an exponential decay curve the asymptotic solution of the diffusion equation had been reached. We find, instead, that a superposition of exponentials can accurately represent the decay curve, consistent with the higher cross-section value.

Drachman, Richard J.↗
DOE OSTI2022

Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver

SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Complex statistics and diffusion in nonlinear disordered particle chains

We investigate dynamically and statistically diffusive motion in a Klein-Gordon particle chain in the presence of disorder. In particular, we examine a low energy (subdiffusive) and a higher energy (self-trapping) case and verify that subdiffusive spreading is always observed. We then carry out a statistical analysis of the motion, in both cases, in the sense of the Central Limit Theorem and present evidence of different chaos behaviors, for various groups of particles. Integrating the equations of motion for times as long as 10{sup 9}, our probability distribution functions always tend to Gaussians and show that the dynamics does not relax onto a quasi-periodic Kolmogorov-Arnold-Moser torus and that diffusion continues to spread chaotically for arbitrarily long times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On the utility of non thermal plasmas

Nonthermal plasmas have properties which can differ substantially even from thermal plasmas of the same density and energy. In practical applications, many of these differences have utility that can be exploited, often with waves. The goal of this thesis is to explore three such applications of nonthermal plasma for the purposes of wave amplification and, ultimately, controlled nuclear fusion. First, we discuss the outstanding potential for a shaped plasma technology that relies on ionizing a beam of aerosol particles. This ionization results in a Coulomb explosion, yielding a highly nonthermal distribution for the plasma in the target volume. Homogeneous, high aspect ratio plasma couplers such as these could be used in a Raman amplification scheme to attain the next generation of laser intensities. Such a powerful laser could be used to study radiation reaction and vacuum pair production, among many other applications. We describe the conditions under which sufficiently dense aerosol beams can be prepared and provide some experimental validation of these claims. Second, we survey the prospects for controlled p-11B fusion in both magnetic and inertial confinement schemes. Whereas most studies have focused on mitigating the bremsstrahlung loss channel, we focus on nonthermal manipulations of the bulk plasma which could improve the reactivity. In particular, we present upper bounds on the p-11B reactivity using alpha channeling to efficiently heat the proton distribution in both schemes. Finally, we place a sharp limit on the energy that can be extracted from a non-thermal particle distribution using waves. In the quasilinear regime of wave-particle interaction, waves rearrange the densities of states in phase space in a diffusive fashion. A salient consequence is that only certain, generally nonthermal, plasma states can be reached. We characterize this restricted state space and thereby reduce the problem of maximizing energy extraction to linear programming. Varying the connectedness of the phase space so that e.g. particles may only be diffused to parts of phase space close in energy, we describe the implications for the state space and the solution of the optimization problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Overcoming a hemihedral twinning problem in tetrahydrofolate-dependent O -demethylase crystals by the microseeding method

A tetrahydrofolate-dependent O -demethylase, LigM, from Sphingobium sp. SYK-6 was crystallized by the hanging-drop vapour-diffusion method. However, the obtained P 3 1 21 or P 3 2 21 crystals, which diffracted to 2.5–3.3 Å resolution, were hemihedrally twinned. To overcome the twinning problem, microseeding using P 3 1 21/ P 3 2 21 crystals as microseeds was performed with optimization of the reservoir conditions. As a result, another crystal form was obtained. The newly obtained crystal diffracted to 2.5–3.0 Å resolution and belonged to space group P 2 1 2 1 2, with unit-cell parameters a = 102.0, b = 117.3, c = 128.1 Å. The P 2 1 2 1 2 crystals diffracted to better than 2.0 Å resolution after optimizing the cryoconditions. Phasing using the single anomalous diffraction method was successful at 3.0 Å resolution with a Pt-derivative crystal. This experience suggested that microseeding is an effective method to overcome the twinning problem, even when twinned crystals are utilized as microseeds.

cryo-protectant↗
DOE OSTI2016

Nonlinear processes upon two-photon interband picosecond excitation of PbWO{sub 4} crystal

A new experimental method is proposed to study the dynamics of nonlinear processes occurring upon two-photon interband picosecond excitation of a lead tungstate crystal and upon its excitation by cw probe radiation in a temporal range from several nanoseconds to several seconds. The method is applied to the case of crystal excitation by a sequence of 25 high-power picosecond pulses with a wavelength of 523.5 nm and 633-nm cw probe radiation. Measuring the probe beam transmittance during crystal excitation, one can investigate the influence of two-photon interband absorption and the thermal nonlinearity of the refractive index on the dynamics of nonlinear processes in a wide range of times (from several nanoseconds to several seconds). The time resolution of the measuring system makes it possible to distinguish fast and slow nonlinear processes of electronic or thermal nature, including the generation of a thermal lens and thermal diffusion. An alternative method is proposed to study the dynamics of induced absorption transformation and, therefore, the dynamics of the development of nonlinear rocesses upon degenerate two-photon excitation of the crystal in the absence of external probe radiation. (nonlinear-optical phenomena)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Some analytical results for radiative transfer in thick atmospheres

Singular eigenmode expansions are a convenient analytical tool with which to study problems of monochromatic radiative transfer in thick or semi-infinite atmospheres. Some closed-form solutions are presented for anisotropic scattering, with the neglect of polarization effects. A basic ingredient for applications to the semi-infinite medium is Chandrasekhar's H-function, which is best defined through the Wiener-Hopf factorization. The solutions of the Milne and albedo problems are discussed in order that the method of matched asymptotic approximations may be used to describe the solution for a thick atmosphere adjacent to a diffusely reflecting ground. Expressions for the emerging distributions are quoted. A possible extension of the general scheme to problems involving polarization is indicated.

Kuscer, I.↗
DOE OSTI2015

On the structure of quasi-stationary laser ablation fronts in strongly radiating plasmas

The effect of strong thermal radiation on the structure of quasi-stationary laser ablation fronts is investigated under the assumption that all the laser flux is absorbed at the critical surface. Special attention is paid to adequate formulation of the boundary-value problem for a steady-state planar ablation flow. The dependence of the laser-to-x-ray conversion efficiency ϕ{sub r} on the laser intensity I{sub L} and wavelength λ{sub L} is analyzed within the non-equilibrium diffusion approximation for radiation transfer. The scaling of the main ablation parameters with I{sub L} and λ{sub L} in the strongly radiative regime 1−ϕ{sub r}≪1 is derived. It is demonstrated that strongly radiating ablation fronts develop a characteristic extended cushion of “radiation-soaked” plasma between the condensed ablated material and the critical surface, which can efficiently suppress perturbations from the instabilities at the critical surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Preliminary Monte Carlo and Thermal Hydraulic Analysis using a Hybrid ETF-Corrected-Diffusion Prediction Block

This paper builds upon previous work to accelerate the Picard iteration (PI) method typically applied for coupled Monte Carlo-Thermal hydraulic (MC-TH) solutions. Previously, the use of the generalized transfer functions (GTFs) to predict variation in macroscopic cross sections following a perturbation in TH properties was demonstrated for a subset of simple 3D problems. In addition, the reduced-order transport prediction block relied on the first order perturbation (FOP) method, which was shown to have computational overheads. Recent work replaced the FOP block with a 1-group nodal diffusion solver to eliminate these overheads. While the use of diffusion is desirable for large-scale problems, the new solver introduces significant homogenization error. This work aims to address this issue by using the Jacobian-Free Newton Krylov (JFNK) method to generate a set of super homogenization (SPH) factors to improve the accuracy of the diffusion solution. The SPH factors will be used in conjunction with an improved cross section prediction method – the expanded transfer function (ETF) method – to produce a highly accurate flux prediction for an axial 1D boiling water reactor (BWR) pincell following a large perturbation in moderator density. The ETF-corrected diffusion (ETF-CD) block is shown to be highly accurate for the 1D test case. Future work will investigate the accuracy of the method for a realistic 3D pressurized water reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Diffusion Synthetic Acceleration for Heterogeneous Domains, Compatible with Voids

A standard approach to solving the S N transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (a) preconditioned source iteration converges rapidly and (b) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using O(1) AMG iterations, convergence of the transport iteration typically requires 2 to 5× fewer iterations than current state-of-the-art “full” DSA, and the proposed method is trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3% of the mesh elements with five to ten AMG iterations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1981

Diffusion in single-phase binary alloys

DBAS 1 computer program provides analyst with simple algorithms for exact rapid solutions of systems with planar, cylindrical, or spherical interfaces. Conventional solutions are complex and present convergence problems. Two algorithm types are figured for each geometry; one converges rapidly for short and the other for long diffusion times. DBAS 1 is written in FORTRAN IV for batch execution.

Tenney, D. R.↗
DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗
DOE OSTI2016

Diffuse coronae in cosmological simulations of Milky Way-sized galaxies

We investigate the properties of halo gas using three cosmological “zoom-in” simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot (T > 10{sup 6} K) halo gas is ∼1% of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot (T=10{sup 5}−10{sup 6} K) gas that extends to the virial radius. Simulated halos agree well outside 20–30 kpc with the β-model of Miller and Bregman based on O vii absorption and O viii emission measurements. Warm-hot and hot gas contribute up to 80% of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for L{sup ∗} galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just 20%–30% of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the “missing baryon problem.” The recovered baryon fraction within 3 virial radii is 90% of the universal value. With a characteristic density of ∼10{sup −4} cm{sup −3} at 50–80 kpc, diffuse coronae meet the requirement for fast and complete ram-pressure stripping of the gas reservoirs in dwarf galaxy satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Risk-Averse Control of Fractional Diffusion with Uncertain Exponent

In this paper, we introduce and analyze a new class of optimal control problems constrained by elliptic equations with uncertain fractional exponents. We utilize risk measures to formulate the resulting optimization problem. We develop a functional analytic framework, study the existence of solution and rigorously derive the first-order optimality conditions. Additionally, we employ a sample-based approximation for the uncertain exponent and the finite element method to discretize in space. Further, we prove the rate of convergence for the optimal risk neutral controls when using quadrature approximation for the uncertain exponent and conclude with illustrative examples.

finite element method↗
NASA NTRS2003

Load Diffusion in Composite and Smart Structures

The research carried out here builds on our previous NASA supported research on the general topic of edge effects and load diffusion in composite structures. Further fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for the multi-functional large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Some specific problems recently considered were those of end effects in smart materials and structures, study of the stress response of pressurized linear piezoelectric cylinders for both static and steady rotating configurations, an analysis of the effect of pre-stressing and pre-polarization on the decay of end effects in piezoelectric solids and investigation of constitutive models for hardening rubber-like materials. Our goal in the study of load diffusion is the development of readily applicable results for the decay lengths in terms of non-dimensional material and geometric parameters. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses.

Horgan, Cornelius O.↗
NASA NTRS2003

Load Diffusion in Composite and Smart Structures

The research carried out here builds on our previous NASA supported research on the general topic of edge effects and load diffusion in composite structures. Further fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for the multi-functional large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Some specific problems recently considered were those of end effects in smart materials and structures, study of the stress response of pressurized linear piezoelectric cylinders for both static and steady rotating configurations, an analysis of the effect of pre-stressing and pre-polarization on the decay of end effects in piezoelectric solids and investigation of constitutive models for hardening rubber-like materials. Our goal in the study of load diffusion is the development of readily applicable results for the decay lengths in terms of non-dimensional material and geometric parameters. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses. The decay behavior of stresses and other field quantities provides a significant aid towards this process. The analysis is also amenable to parameter study with a large parameter space and should be useful in structural tailoring studies. Special purpose analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and in assessing results from general purpose finite element analyses. For example, a rational basis is needed in choosing where to use three-dimensional to two-dimensional transition finite elements in analyzing stiffened plates and shells. The decay behavior of stresses and other field quantities furnished by this research provides a significant aid towards this element transition issue. A priori knowledge of the extent of boundary-layers induced by edge effects is also useful in determination of the instrumentation location in structural verification tests or in material characterization tests.

Horgan, C. O.↗