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 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Space plasma turbulent dissipation - Reality or myth?

A prevalent approach to understanding magnetospheric dynamics is to combine a hydromagnetic description of the large scale magnetic structure and convection flows with a locally determined anomalous dissipation which develops in boundary layers. Three problems (nose and tail reconnection, auroral field-aligned currents, and diffuse auroral precipitation) are critically examined to test the validity of this theoretical philosophy. Although the expected plasma wave turbulence is observed for each case, the concept of local anomalous dissipation fails to provide an adequate or complete description of the phenomenae.

Coroniti, F. V.↗

NASA NTRS1986

Radiation transport in dust in disk geometry

The main objective of the research program is twofold: (1) to develop a computer code to solve the problem of scattering, absorption and emission of photons by dust grains in a dusty medium with 2 dimensional disk geometry, and (2) to study the various physical and geometrical effects of 2 dimensional radiation transport on the thermal structure and radiation field. These tasks were accomplished and are briefly summarized. The method for solving the radiation transport problem in disk geometry is a generalization of the quasi-diffusion method (QDM) previously developed by the author.

Chun, Ming Leung↗

NASA NTRS1990

Diffuse reflection from a stochastically bounded, semi-infinite medium

In order to determine the diffuse reflection from a medium bounded by a rough surface, the problem of radiative transfer in a boundary layer characterized by a statistical distribution of heights is considered. For the case that the surface is defined by a multivariate normal probability density, the propagation probability for rays traversing the boundary layer is derived and, from that probability, a corresponding radiative transfer equation. A solution of the Eddington (two stream) type is found explicitly, and examples are given. The results should be applicable to reflection from the regoliths of solar system bodies, as well as from a rough ocean surface.

Lumme, K.↗

NASA NTRS1993

The application of preconditioning in viscous flows

The present time-derivative preconditioning algorithm is effective in flow conditions ranging from inviscid to very diffusive flows, as well as low subsonic to supersonic flow velocities. By means of a preconditioning matrix that (1) introduces well-conditioned eigenvalues and (2) avoids nonphysical time reversals for viscous flows, a mechanism is obtained which controls the inviscid and viscous time-step parameters at very diffusive flows. These capabilities are demonstrated for a variety of sample problems; convergence rates of solutions that are indistinguishable from those obtained without preconditioning are shown to be accelerated by as much as two orders of magnitude.

Choi, Y.-H.↗

NASA NTRS1994

A class of the van Leer-type transport schemes and its application to the moisture transport in a general circulation model

A generalized form of the second-order van Leer transport scheme is derived. Several constraints to the implied subgrid linear distribution are discussed. A very simple positive-definite scheme can be derived directly from the generalized form. A monotonic version of the scheme is applied to the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) for the moisture transport calculations, replacing the original fourth-order center-differencing scheme. Comparisons with the original scheme are made in idealized tests as well as in a summer climate simulation using the full GLA GCM. A distinct advantage of the monotonic transport scheme is its ability to transport sharp gradients without producing spurious oscillations and unphysical negative mixing ratio. Within the context of low-resolution climate simulations, the aforementioned characteristics are demonstrated to be very beneficial in regions where cumulus convection is active. The model-produced precipitation pattern using the new transport scheme is more coherently organized both in time and in space, and correlates better with observations. The side effect of the filling algorithm used in conjunction with the original scheme is also discussed, in the context of idealized tests. The major weakness of the proposed transport scheme with a local monotonic constraint is its substantial implicit diffusion at low resolution. Alternative constraints are discussed to counter this problem.

Lin, Shian-Jiann↗

NASA NTRS2022

Multi-Scale Modelling of the Bound Metal Deposition Manufacturing of Ti6Al4V

Nonlinear shrinkage of the metal part during manufacturing by bound metal deposition, both on the ground and under microgravity, is considered. A multi-scale physics-based approach is developed to address the problem. It spans timescales from atomistic dynamics on the order of nanoseconds to full-part shrinkage on the order of hours. This approach enables estimation of the key parameters of the problem, including the widths of grain boundaries, the coefficient of surface diffusion, the initial redistribution of particles during the debinding stage, the evolution of the microstructure from round particles to densely-packed grains, the corresponding changes in the total and chemical free energies, and the sintering stress. The method has been used to predict shrinkage at the levels of two particles, of the filament cross-section, of the sub-model, and of the whole green, brown, and metal parts.

Nonlinear shrinkage↗

DOE OSTI2024

Investigation of Local Parameters of PERC Solar Cells Metallized with Screen Printed Cu-Paste

For 40 TW of PV required to transition our planet to 100% renewables, the silver (Ag) should disappear from PV production. Advantages of copper (Cu) over silver (Ag) include: 1) bulk Cu has a similar conductivity to Ag (1.7 one millionth O-cm for Cu, 1.6 one millionth O-cm for Ag, and 2) Cu is -100 times cheaper than Ag, making it an excellent potential replacement. Problems associated with copper contacts include: 1) easy oxidation, and 2) diffusion into the Si cell and recombination activity. To summarize: 1) Rs map shows some regions with very high Rs, indicating no contact in those areas; 2) the histogram shows the peak Rs -5 ohm.cm2, which explains the high FF loss due to Rs; and 3) DLIT indicates non-uniformity in J01 and J02.

copper↗

DOE OSTI2026

High-Fidelity, Low-Dissipation/Symmetry-Preserving Numerical Scheme for Solving the Euler Equations with Unstructured, Metric-Based Mesh Adaptation

This work presents an overview of a high-fidelity compressible Euler solver that utilizes the continuous Galerkin (CG) method with added artificial numerical diffusion for stabilization to solve a variety of unsteady and steady benchmark inviscid flow problems. This work shows that discretizing the Euler equations with this CG approach and first order basis functions produces a cost-effective stencil as well as simple well-posed boundary conditions. We show through convergence testing with manufactured solutions that the reduced stencil of CG, combined with the low amount of artificial diffusion required when using the stabilization method outlined in this work, leads to stable and highly accurate results for a variety of unsteady and steady applications. When combined with the adaptive mesh refinement approach used for many of the cases in this work, our results show that the flow solver achieves even more accurate results. A variety of inviscid flow cases are presented in this work, including transient 2D cases with complex shock structures and several steady 3D airfoils sections with a constant span.

Doetsch, Kevin [ORNL] (ORCID:0000000267051705)↗

NASA NTRS1975

A nonlinear theory of cosmic ray pitch angle diffusion in homogeneous magnetostatic turbulence

A plasma strong turbulence, weak coupling, theory is applied to the problem of cosmic ray pitch angle scattering in magnetostatic turbulence. The theory used is a rigorous generalization of Weinstock's resonance-broadening theory and contains no ad hoc approximations. A detailed calculation is presented for a model of slab turbulence with an exponential correlation function. The results agree well with numerical simulations. The rigidity dependence of the pitch angle scattering coefficient differs from that found by previous researchers. The differences result from an inadequate treatment of particle trajectories near 90 deg pitch angle in earlier work.

Goldstein, M. L.↗

NASA NTRS1976

A nonlinear theory of cosmic-ray pitch-angle diffusion in homogeneous magnetostatic turbulence

A plasma strong turbulence, weak coupling, theory is applied to the problem of cosmic-ray pitch-angle scattering in magnetostatic turbulence. The theory used is a rigorous generalization of Weinstock's 'resonance broadening' theory and contains no ad hoc approximations. A detailed calculation is presented for a model of 'slab' turbulence with an exponential correlation function. The results agree well with numerical simulations. The rigidity dependence of the pitch-angle scattering coefficient differs from that found by previous researchers. The differences result from an inadequate treatment of particle trajectories near 90 deg pitch angle in earlier work.

Goldstein, M. L.↗

NASA NTRS2002

Mathematical model of diffusion-limited gas bubble dynamics in unstirred tissue with finite volume

Models of gas bubble dynamics for studying decompression sickness have been developed by considering the bubble to be immersed in an extravascular tissue with diffusion-limited gas exchange between the bubble and the surrounding unstirred tissue. In previous versions of this two-region model, the tissue volume must be theoretically infinite, which renders the model inapplicable to analysis of bubble growth in a finite-sized tissue. We herein present a new two-region model that is applicable to problems involving finite tissue volumes. By introducing radial deviations to gas tension in the diffusion region surrounding the bubble, the concentration gradient can be zero at a finite distance from the bubble, thus limiting the tissue volume that participates in bubble–tissue gas exchange. It is shown that these deviations account for the effects of heterogeneous perfusion on gas bubble dynamics, and are required for the tissue volume to be finite. The bubble growth results from a difference between the bubble gas pressure and an average gas tension in the surrounding diffusion region that explicitly depends on gas uptake and release by the bubble. For any given decompression, the diffusion region volume must stay above a certain minimum in order to sustain bubble growth.

NASA Discipline Environmental Health↗

DOE OSTI2014

Steady fall of isothermal, resistive-viscous, compressible fluid across magnetic field

This is a basic MHD study of the steady fall of an infinite, vertical slab of isothermal, resistive-viscous, compressible fluid across a dipped magnetic field in uniform gravity. This double-diffusion steady flow in unbounded space poses a nonlinear but numerically tractable, one-dimensional (1D) free-boundary problem, assuming constant coefficients of resistivity and viscosity. The steady flow is determined by a dimensionless number μ{sub 1} proportional to the triple product of the two diffusion coefficients and the square of the linear total mass. For a sufficiently large μ{sub 1}, the Lorentz, viscous, fluid-pressure, and gravitational forces pack and collimate the fluid into a steady flow of a finite width defined by the two zero-pressure free-boundaries of the slab with vacuum. The viscous force is essential in this collimation effect. The study conjectures that in the regime μ{sub 1}→0, the 1D steady state exists only for μ{sub 1}∈Ω, a spectrum of an infinite number of discrete values, including μ{sub 1} = 0 that corresponds to two steady states, the classical zero-resistivity static slab of Kippenhahn and Schlüter [R. Kippenhahn and A. Schlüter, Z. Astrophys. 43, 36 (1957)] and its recent generalization [B. C. Low et al., Astrophys. J. 755, 34 (2012)] to admit an inviscid resistive flow. The pair of zero-pressure boundaries of each of the μ{sub 1}→0 steady-state slabs are located at infinity. Computational evidence suggests that the Ω steady-states are densely distributed around μ{sub 1} = 0, as an accumulation point, but are sparsely separated by open intervals of μ{sub 1}-values for which the slab must be either time-dependent or spatially multi-dimensional. The widths of these intervals are vanishingly small as μ{sub 1}→0. This topological structure of physical states is similar to that described by Landau and Liftshitz [L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Addison-Wesley, Reading, MA, 1959)] to explain the onset of hydrodynamic turbulence. The implications of this MHD study are discussed, with an interest in the prominences in the solar atmosphere and the interstellar clouds in the Galaxy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the polynomial order used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Verification of the DIF3D Software to Support Fast Reactor Analysis (Rev. 3)

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1972

Applications of NASTRAN to nuclear problems

The extent to which suitable solutions may be obtained for one physics problem and two engineering type problems is traced. NASTRAN appears to be a practical tool to solve one-group steady-state neutron diffusion equations. Transient diffusion analysis may be performed after new levels that allow time-dependent temperature calculations are developed. NASTRAN piecewise linear anlaysis may be applied to solve those plasticity problems for which a smooth stress-strain curve can be used to describe the nonlinear material behavior. The accuracy decreases when sharp transitions in the stress-strain relations are involved. Improved NASTRAN usefulness will be obtained when nonlinear material capabilities are extended to axisymmetric elements and to include provisions for time-dependent material properties and creep analysis. Rigid formats 3 and 5 proved to be very convenient for the buckling and normal-mode analysis of a nuclear fuel element.

Spreeuw, E.↗

DOE OSTI2015

Driven microswimmers on a 2D substrate: A stochastic towed sled model

We investigate, both numerically and analytically, the diffusion properties of a stochastic sled sliding on a substrate, subject to a constant towing force. The problem is motivated by the growing interest in controlling transport of artificial microswimmers in 2D geometries at low Reynolds numbers. We simulated both symmetric and asymmetric towed sleds. Remarkable properties of their mobilities and diffusion constants include sidewise drifts and excess diffusion peaks. We interpret our numerical findings by making use of stochastic approximation techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Summer term 2018 work with IMC radiative transport in binary stochastic media

In the quest for new and improved methods to solve thermal radiation transport in stochastic media, this summer project was to characterize radation transport in a deterministic lattice problem that represents a real, binary stochastic media. A representative geometry for a binary stochastic media consisting of a cube of a high volume fraction media containing constant-diameter non-overlappying spheres randomly distributed within was run previously. Previous thermal radiation wave problems in binary stochastic media displayed three characteristic thermal energy deposition wavefronts in the following order, front of the spheres, bulk material, and the backs of the spheres. Several physical aspects of the representative geometry for a stochastic binary media using an optically thin background with optically thick spheres were analyzed. This analysis includes wave speed propagation using a one dimensional diffusion approximation, energy fluence transmission, and material energy deposition. The approaches to analysis included representing one sphere of the high opacity material in a diffusive background of the low opacity material, as well as repeating spheres in a column con guration with reflective boundary conditions to approximate a lattice. Additionally, both atomic mix and Levermore-Pomraning homogenization schemes were compared against the representative cube geometry as well as the sphere column lattice geometry. The wave speed analysis with a one dimensional Marshak diffusion wave was largely ineffective at determining an effective opacity value to which a problem could be homogenized, likely due to low mean free path resolution within the problem. However, the sphere column lattice geometry showed promising results for approximating the energy fluence transmitted through the problem domain. One concern with the sphere column lattice geometry indicated the high forward-bias effect of fluence reaching the problem boundary due to streaming paths around the spheres in the case of small mean free path resolution of the background material.

36 MATERIALS SCIENCE↗