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,621 records · Page 90

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗

Method for improving the performance of oxidizable ceramic materials in oxidizing environments

Improved adhesion of thermal barrier coatings to nonmetallic substrates using a dense layer of ceramic on an underlying nonmetallic substrate that includes at least one oxidizable component. The improved adhesion occurs because the application of the dense ceramic layer forms a diffusion barrier for oxygen. This diffusion barrier prevents the oxidizable component of the substrate from decomposing. The present invention applies ceramic by a process that deposits a relatively thick and dense ceramic layer on the underlying substrate. The formation of the dense layer of ceramic avoids the problem of void formation associated with ceramic formation by most prior art thermal decomposition processes. The formation of voids has been associated with premature spalling of thermal barrier layers and other protective layers applied to substrates.

Nagaraj, Bangalore A.↗

Comparison of subchannel and averaged channel thermal-hydraulic descriptions on coupled pin-by-pin neutronic calculations

Best-estimate transient analysis of Light Water Reactors (LWRs) requires high accuracy code systems solving neutronics (NK), thermal-hydraulics (TH) and thermal-mechanics to model the core behavior. For long time, nodal diffusion methods for solution of the transport equation using assembly size computational meshes have been applied neglecting the details within the fuel assembly, or using 'pin power reconstruction' techniques to predict those details. Moreover, the multiphysics analysis has been commonly performed coupling those neutronic diffusion methods to thermal-hydraulic system codes using a nodalization of several fuel assemblies per computational node without taking into account the effects that this type of averaging can have on the results obtained. This paper presents the work accomplished with the coupled system COBAYA4/CTF for a preliminary analysis of the effect of the different accuracy levels in the NK and TH domain. COBAYA4 is a neutronic code able to solve the diffusion equation at two resolution levels, nodal and pin by pin (PbP). CTF is a thermo-hydraulic code able to solve the coolant conservation equations also at two levels, channel (fuel assembly averaged solution) and sub-channel (real cooling channel) level. The coupled system COBAYA4/CTF is used through SALOME platform, i.e. SALOME is the software able to map conveniently the NK and TH meshes. This coupling allows combining different resolution levels for NK and TH problems. Therefore, the study of the effect of using different resolutions for the NK and TH nodalizations has been performed applying the system to a fuel assembly with selected boundary conditions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic histories of dust grains in the interstellar medium

The purpose is to study an evolving system of refractory dust grains within the Interstellar Medium (ISM). This is done via a combination of Monte Carlo processes and a system of partial differential equations, where refractory dust grains formed within supernova remnants and ejecta from high mass loss stars are subjected to the processes of sputtering and collisional fragmentation in the diffuse media and accretion within the cold molecular clouds. In order to record chemical detail, the authors take each new particle to consist of a superrefractory core plus a more massive refractory mantle. The particles are allowed to transfer to and fro between the different phases of the interstellar medium (ISM) - on a time scale of 10(exp 8) years - until either the particles are destroyed or the program finishes at a Galaxy time of 6x10(exp 9) years. The resulting chemical and size spectrum(s) are then applied to various astrophysical problems with the following results. For an ISM which has no collisional fragmentation of the dust grains, roughly 10 percent by mass of the most refractory material survives the rigors of the ISM intact, which leaves open the possibility that fossilized isotopically anomalous material may have been present within the primordial solar nebula. Stuctured or layered refractory dust grains within the model cannot explain the observed interstellar depletions of refractory material. Fragmentation due to grain-grain collisions in the diffuse phase plus the accretion of material in the molecular cloud phase can under certain circumstances cause a bimodal distribution in grain size.

Liffman, Kurt↗

Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple

The development of more efficient thermoelectric couple technology capable of operating with high-grade heat sources up to 1,275 K is key to improving the performance of radioisotope thermoelectric generators. Lanthanum telluride La3-xTe4 and 14-1-11 Zintls (Yb14MnSb11) have been identified as very promising materials. The fabrication of advanced high-temperature thermoelectric couples requires the joining of several dissimilar materials, typically including a number of diffusion bonding and brazing steps, to achieve a device capable of operating at elevated temperatures across a large temperature differential (up to 900 K). A thermoelectric couple typically comprises a heat collector/ exchanger, metallic interconnects on both hot and cold sides, n-type and ptype conductivity thermoelectric elements, and cold-side hardware to connect to the cold-side heat rejection and provide electrical connections. Differences in the physical, mechanical, and chemical properties of the materials that make up the thermoelectric couple, especially differences in the coefficients of thermal expansion (CTE), result in undesirable interfacial stresses that can lead to mechanical failure of the device. The problem is further complicated by the fact that the thermoelectric materials under consideration have large CTE values, are brittle, and cracks can propagate through them with minimal resistance. The inherent challenge of bonding brittle, high-thermal-expansion thermoelectric materials to a hot shoe material that is thick enough to carry the requisite electrical current was overcome. A critical advantage over prior art is that this device was constructed using all diffusion bonds and a minimum number of assembly steps. The fabrication process and the materials used are described in the following steps: (1) Applying a thin refractory metal foil to both sides of lanthanum telluride. To fabricate the n-type leg of the advanced thermoelectric couple, the pre-synthesized lanthanum telluride coupon was diffusion bonded to the metal foil using a thin adhesion layer. (2) Repeating a similar process for the 14-1-11 Zintl p-type leg of the advanced thermoelectric couple. (3) Bonding thick CTE-matched metal plates on the metallized lanthanum telluride and Yb14MnSb11 to form the hot and cold sides of the thermoelectric couple. The calculated conversion efficiency of such an advanced couple would be about 10.5 percent, about 35 percent better than heritage radioisotope thermoelectric technology that relies on Si-Ge alloys. In addition, unlike Si-Ge alloys, these materials can be combined with many other thermoelectric materials optimized for operation at lower temperatures to achieve conversion efficiency in excess of 15 percent (a factor of 2 increase over heritage technology).

Ravi, Vilupanur A.↗

Colloidal State Machines as Smart Tracers for Chemical Reactor Analysis

A widely utilized tool in reactor analysis is passive tracers that report the residence time distribution, allowing estimation of the conversion and other properties of the system. Recently, advances in microrobotics have introduced powered and functional entities with sizes comparable to some traditional tracers. This has motivated the concept of Smart Tracers that could record the local chemical concentrations, temperature, or other conditions as they progress through reactors. Herein, the design constraints and advantages of Smart Tracers by simulating their operation in a laminar flow reactor model conducting chemical reactions of various orders are analyzed. It is noted that far fewer particles are necessary to completely map even the most complex concentration gradients compared with their conventional counterparts. Design criteria explored herein include sampling frequency, memory storage capacity, and ensemble number necessary to achieve the required accuracy to inform a reactor model. Cases of severe particle diffusion and sensor noise appear to bind the functional upper limit of such probes and require consideration for future design. The results of the study provide a starting framework for applying the new technology of microrobotics to the broad and impactful set of problems classified as chemical reactor analysis.

97 MATHEMATICS AND COMPUTING↗

Surface transport processes in charged porous media

Surface transport processes are important in chemistry, colloidal sciences, engineering, biology, and geophysics. Natural or externally produced charges on surfaces create electrical double layers (EDLs) at the solid-liquid interface. The existence of the EDLs produces several complex processes including bulk and surface transport of ions. In this work, a model is presented to simulate bulk and transport processes in homogeneous porous media comprising big pores. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A volume averaging technique is used to derive the averaged transport equations in the limit of thin electrical double layers. Description of the EDL between the electrolyte solution and the charged wall is accomplished using the Gouy-Chapman-Stern (GCS) model. The surface transport terms enter into the average equations due to the use of boundary conditions for diffuse interfaces. Two extra surface transports terms appear in the closed average equations. One is a surface diffusion term equivalent to the transport process in non-charged porous media. The second surface transport term is a migration term unique to charged porous media. The effective bulk and transport parameters for isotropic porous media are calculated solving the corresponding closure problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multipoint Correlations in Poisson Media

Multipoint correlations in randomly heterogeneous composite media are determined by the probability that a set of points belong to specific phases. They determine a wide range of macroscopic transport properties such as conductivity, dielectric constant, diffusion coefficient, and transmittance. The Poisson model—a random tesselation of space by hyperplanes—provides realistic descriptions of heterogeneous media in, e.g., radiation-transport applications; yet, until now, it has lacked closed-form expressions for its multipoint correlations. We resolve this problem by presenting an exact solution for the multipoint correlations in the Poisson model. By comparing it to Monte Carlo simulations of four-point correlations in three dimensions, we demonstrate the accuracy of our solution. In conclusion, we visualize the multipoint correlations and discuss their features.

Amorphous materials↗

Multischeme transport with upwinding for Rattlesnake code

This report conclude our efforts at INL to make Rattlesnake, the MOOSE-based radiation multigroup transport application, fully support the multischeme transport capability. Specifically, it details all five types of interface condition with upwinding method for coupling eight available discretization schemes within Rattlesnake. These discretization schemes include SN (discrete ordinates method), PN (spherical harmonics expansion method) and diffusion for angular discretization, and continuous finite element method (CFEM) and discontinuous finite element method (DFEM) for spatial discretization. The other type of method based on Lagrange multipliers for coupling schemes based on CFEM has been studied in the past. Besides simple test problems, a more realistic application of multischeme transport for TREAT (transient reactor test facility) experiment modeling are presented. Several enhancements are reviewed at the end which could potentially be useful in the future.

97 MATHEMATICS AND COMPUTING↗

Major unresolved issues in the formation of chondrules and CAIs

So much has been learned and published about chondrites and so many models have been floated (but neither proven nor refuted) that the subject has become extremely diffuse and amorphous. It is impossible to remain simultaneously aware of all the many constraints, from many different subdisciplines, that bear on nebular processes. This paper attempts to impose some structure on the situation and loosen the interpretational logjam by breaking the chondrite problem into a series of well-defined questions. Most are open questions; those that some might consider closed by now may deserve formal reexamination.

Wood, J. A.↗

On Multi-Dimensional Unstructured Mesh Adaption

Anisotropic unstructured mesh adaption is developed for a truly multi-dimensional upwind fluctuation splitting scheme, as applied to scalar advection-diffusion. The adaption is performed locally using edge swapping, point insertion/deletion, and nodal displacements. Comparisons are made versus the current state of the art for aggressive anisotropic unstructured adaption, which is based on a posteriori error estimates. Demonstration of both schemes to model problems, with features representative of compressible gas dynamics, show the present method to be superior to the a posteriori adaption for linear advection. The performance of the two methods is more similar when applied to nonlinear advection, with a difference in the treatment of shocks. The a posteriori adaption can excessively cluster points to a shock, while the present multi-dimensional scheme tends to merely align with a shock, using fewer nodes. As a consequence of this alignment tendency, an implementation of eigenvalue limiting for the suppression of expansion shocks is developed for the multi-dimensional distribution scheme. The differences in the treatment of shocks by the adaption schemes, along with the inherently low levels of artificial dissipation in the fluctuation splitting solver, suggest the present method is a strong candidate for applications to compressible gas dynamics.

Wood, William A.↗

How Does the Electron Dynamics Affect the Reconnection Rate in a Typical Reconnection Layer?

The question of whether the microscale controls the macroscale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the properties of electron diffusion region dynamics both for guide field and anti-parallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

How Does the Electron Dynamics Affect the Global Reconnection Rate

The question of whether the microscale controls the macroscale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the properties of electron diffusion region dynamics both for guide field and anti-parallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using monotonic basin hopping (MBH) have become well developed during the past decade [1, 2, 3, 4, 5, 6]. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing random variable (RV)s from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by J. Englander [3, 6]) significantly improves monotonic basin hopping (MBH) performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness. Efficiency is finding better solutions in less time. Robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive random walks (RWs) originally developed in the field of statistical physics.

autonomous↗

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using MBH have become well developed during the past decade. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing RVs from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by Englander significantly improves MBH performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness, where efficiency is finding better solutions in less time, and robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive RWs originally developed in the field of statistical physics.

mission design↗