Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chemically controlled pattern formation in self-oscillating elastic shells

Significance Mechanochemical processes drive various biological functions including morphogenesis and morphological changes in responsive synthetic hydrogels. By allowing reactive components to function in an elastic closed shell, autonomous chemical reactions generate microcompartments capable of expanding or contracting or induce buckled patterns. This capability can drive specific functions such as enhancing the rate of diffusion of compartmentalized chemicals or releasing cargos at specific rates. Here, we demonstrate the coupling of mechanical response to autonomous chemical reactions in elastic shells which swell and deswell to generate specific dynamic patterns and morphologies. The reverse feedback from the mechanical input triggers the chemical reaction. Our work inspires the design of responsive shells to enhance the functionality of microcompartments and nanoreactors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1972

Instabilities and turbulence in highly ionized plasmas in a magnetic field

Physical mechanisms were considered which are responsible for plasma turbulence and the establishment of necessary conditions for energy exchange and transfer in the frequency spectrum. In addition, work was performed to better understand the drift instability in the highly inhomogeneous Rensselaer arc, and methods to suppress this instability using feedback stabilization techniques. Correlation techniques were refined to study plasma turbulence, the diffusion wave technique for monitoring cross-field diffusion was extended to include regimes of high turbulence levels, and a technique for coupling stabilizing RF power to the Rensselaer arc was developed.

Jennings, W. C.↗
NASA NTRS1974

Thermal coupling measurement method

Heat flow from an embedded heated wire responds to a change in the ambient environment. The wire is part of a self-balancing bridge system, and heat flow is measured directly in watts. Steady-state and transient thermal coupling can be measured directly and is an indication of the thermal resistance and diffusivity for the system under study. The method is applied to an aerospace electroexplosive component.

Rosenthal, L. A.↗
NASA NTRS1988

A space-marching method for the computation of viscous internal flows

A space-marching method has been developed to compute 3-D viscous flows in internal geometries. The Navier-Stokes equations have been posed as an initial-value problem by neglecting the effects of streamwise diffusion and treating the streamwise pressure gradient as a known source term. The fully coupled system of equations has been solved by a noniterative algorithm at each streamwise step of the computation. A low Mach number formulation of the equations has been used to compute incompressible flow fields. A computer program has been written to implement all aspects of the space-marching algorithm. The program is modular and is easily adapted to the widely varying geometries of internal flows. The space-marching algorithm has been tested by computing simple flows with known analytical solutions. The method has been used to predict complex 3-D turbulent flows. The algorithm is stable and very economical. A single sweep of the flow field by the space-marching method is approximately equivalent to one time-step of the time-marching method.

Govindan, T. R.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 2: User's guide

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the User's Guide, and describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. The governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models are described in detail.

Towne, Charles E.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This User's Guide describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗
NASA NTRS2012

VIIRS On-Orbit Optical Anomaly - Investigation, Analysis, Root Cause Determination and Lessons Learned

A gradual, but persistent, decrease in the optical throughput was detected during the early commissioning phase for the Suomi National Polar-Orbiting Partnership (SNPP) Visible Infrared Imager Radiometer Suite (VIIRS) Near Infrared (NIR) bands. Its initial rate and unknown cause were coincidently coupled with a decrease in sensitivity in the same spectral wavelength of the Solar Diffuser Stability Monitor (SDSM) raising concerns about contamination or the possibility of a system-level satellite problem. An anomaly team was formed to investigate and provide recommendations before commissioning could resume. With few hard facts in hand, there was much speculation about possible causes and consequences of the degradation. Two different causes were determined as will be explained in this paper. This paper will describe the build and test history of VIIRS, why there were no indicators, even with hindsight, of an on-orbit problem, the appearance of the on-orbit anomaly, the initial work attempting to understand and determine the cause, the discovery of the root cause and what Test-As-You-Fly (TAYF) activities, can be done in the future to greatly reduce the likelihood of similar optical anomalies. These TAYF activities are captured in the lessons learned section of this paper.

Iona, Glenn↗
DOE OSTI2021

Measurement techniques of thermoelectric devices and modules

Thermoelectric (TE) energy conversion devices are designed based on measured material transport properties. Although an ideal device efficiency can be calculated using the figure of merit, zT, the actual efficiency of a thermoelectric module is smaller and must be measured. TE device performance testing techniques and procedures are reviewed with a focus on energy conversion efficiency. One of the areas that attract less attention in the literature is the transition from basic thermoelectric materials research to the production of modules and devices. The effects of interface layers, diffusion barriers, thermal insulations, and heat sinks are reviewed. Finally, device performance from single couples to standard modules and long-term reliability are discussed.

Wang, Hsin↗
DOE OSTI2022

Graphite Oxidation Modeling

Oxidation modeling importance ASME consideration Effect on graphite properties Model Overview Numerical Model Oxidation Diffusion Reaction Kinetics Model Parameterization Example problems Editing the model Simple problem example Conclusions Coupling to stress analysis Oxidation result analysis

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
DOE OSTI2021

Strong Local Evaporative Cooling Over Land Due to Atmospheric Aerosols

Abstract Aerosols can enhance terrestrial productivity through increased absorption of solar radiation by the shaded portion of the plant canopy—the diffuse radiation fertilization effect. Although this process can, in principle, alter surface evaporation due to the coupling between plant water loss and carbon uptake, with the potential to change the surface temperature, aerosol‐climate interactions have been traditionally viewed in light of the radiative effects within the atmosphere. Here, we develop a modeling framework that combines global atmosphere and land model simulations with a conceptual diagnostic tool to investigate these interactions from a surface energy budget perspective. Aerosols increase the terrestrial evaporative fraction, or the portion of net incoming energy consumed by evaporation, by over 4% globally and as much as ∼40% regionally. The main mechanism for this is the increase in energy allocation from sensible to latent heat due to global dimming (reduction in global shortwave radiation) and slightly augmented by diffuse radiation fertilization. In regions with moderately dense vegetation (leaf area index >2), the local surface cooling response to aerosols is dominated by this evaporative pathway, not the reduction in incident radiation. Diffuse radiation fertilization alone has a stronger impact on gross primary productivity (+2.18 Pg C y −1 or +1.8%) than on land evaporation (+0.18 W m −2 or +0.48%) and surface temperature (−0.01 K). Our results suggest that it is important for land surface models to distinguish between quantity (change in total magnitude) and quality (change in diffuse fraction) of radiative forcing for properly simulating surface climate.

Chakraborty, TC↗
DOE OSTI2023

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Elucidation of Active Sites in Aldol Condensation of Acetone over Single-Facet Dominant Anatase TiO 2 (101) and (001) Catalysts

Aldol condensations of carbonyl compounds for C–C bond formation are a very important class of reactions in organic synthesis and upgrading of biomass-derived feedstocks. However, the atomic level understanding of reaction mechanisms and structure-activity correlation on widely used transition metal oxide catalysts are limited due to the high degree of struc-tural heterogeneity of catalysts such as commercial TiO 2 powders. Here, we provide a deep understanding of the reaction mechanisms, kinetics, and structure–function relationships for vapor phase acetone aldol condensation through the con-trolled synthesis of two catalysts with high surface areas and clean, dominant facets, coupled with detailed characterization and kinetic studies that are further assisted by density functional theory (DFT) calculations. Temperature-dependent diffuse reflectance infrared Fourier transform spectroscopy showed the existence of abundant acetone bonded to surface hydroxyl groups (acetone-O s H) and acetone bonded to Lewis acid sites (acetone-Ti 5c ) on surface of both {101} and {001} facet dom-inant TiO 2 . Intermolecular C–C coupling of enolate intermediate from acetone-Ti 5c and a vicinal acetone-O s H is a kinetically relevant step, which is consistent with kinetic and isotopic studies and DFT calculations. The {001} facet showed a lower apparent activation energy (or higher reactivity) than the {101} facet. This is likely caused by the weaker Lewis acid and Brønsted base strengths of the {001} facet which favors the reprotonation–desorption of the coupled intermediate, making the C–C coupling step more exothermic on the {001} facet and resulting in an earlier transition state with a lower activation barrier. It is also possible that the {001} facet has a smoother surface configuration and less steric hindrance during intermo-lecular C–C bond formation than the {101} facet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗