Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogenization theory”

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 685 records · Page 38

Structure of turbulence at high shear rate

The structure of homogeneous turbulence subject to high shear rate has been investigated by using three-dimensional, time-dependent numerical simulations of the Navier-Stokes equations. This study indicates that high shear rate alone is sufficient for generation of the streaky structures, and that the presence of a solid boundary is not necessary. Evolution of the statistical correlations is examined to determine the effect of high shear rate on the development of anisotropy in turbulence. It is shown that the streamwise fluctuating motions are enhanced so profoundly that a highly anisotropic turbulence state with a 'one-component' velocity field and 'two-component' vorticity field develops asymptotically as total shear increases. Because of high-shear rate, rapid distortion theory predicts remarkably well the anisotropic behavior of the structural quantities.

Lee, Moon Joo↗

Chapter 21. Atomic-Level Dislocation Dynamics in Irradiated Metals

Primary damage and microstructure evolution in structural nuclear materials operating under conditions of a high flux of energetic atomic particles and high temperature and stress lead to the formation of a high concentration, non-homogeneous distribution of defect clusters in the form of dislocation loops, voids, gas-filled bubbles and radiation-induced precipitates of nanometer scale. They cause changes in many material properties. Being obstacles to dislocation glide, they strongly affect the mechanical properties, with an increase in yield and flow stresses and a reduction in ductility. Atomic-scale computer simulations can provide details of how these effects are influenced by the obstacle structure, applied stress, strain rate and temperature. Processes such as obstacle cutting, transformation, absorption and drag are observed. Some recent results for body-centered and face-centered cubic metals are described in this review and, where appropriate, comparisons are drawn with predictions based on the elasticity theory of crystal defects. Perspectives on how to use this information at higher scales and in particular in mesoscale, dislocation dynamics simulations are also discussed.

Osetskiy, Yury N.↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective field theory for quasicrystals and phasons dynamics

We build an effective field theory (EFT) for quasicrystals -- aperiodic incommensurate lattice structures -- at finite temperature, entirely based on symmetry arguments and a well-define action principle. By means of Schwinger-Keldysh techniques, we derive the full dissipative dynamics of the system and we recover the experimentally observed diffusion-to-propagation crossover of the phason mode. From a symmetry point of view, the diffusive nature of the phason at long wavelengths is due to the fact that the internal translations, or phason shifts, are symmetries of the system with no associated Noether currents. The latter feature is compatible with the EFT description only because of the presence of dissipation (finite temperature) and the lack of periodic order. Finally, we comment on the similarities with certain homogeneous holographic models and we formally derive the universal relation between the pinning frequency of the phonons and the damping and diffusion constant of the phason.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The radial evolution of the amplitudes of 'dissipationless' turbulent solar wind fluctuations

We examine the evolution with heliocentric distance of the amplitude of fluctuations in the interplanetary magnetic and velocity fields assuming a model of homogeneous, steady state turbulence. Both the Alfven ratio and the normalized cross helicity are taken to be constant, compared to other variations, and the turbulence is assumed to be nearly incompressive. The fluctuation amplitudes are found by assuming the heating of the plasma by the fluctuations to be negligible; quasi-steady state turbulence with damping balanced by large-scale energy input will lead to the same conclusions for inertial range fluctuations. While the assumptions of this model contrast sharply with those for purely Alfvenic fluctuations, the radial dependence of the amplitude of the fluctuations for reasonable parameters is very nearly that found from both WKB analysis and recent turbulence modeling. The robustness of this result suggests why some predictions of WKB theory are apparently correct in solar wind conditions where the theory is not expected to be valid.

Verma, M. K.↗

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗

Geochemical comparison of K-T boundaries from the Northern and Southern Hemispheres

Closely spaced (cm-scale) traverses through the K-T boundary at Stevns Klint (Denmark), Woodside Creek (New Zealand) and a new Southern Hemisphere site at Richards Bay (South Africa) were subjected to trace element and isotopic (C, O, Sr) investigation. Intercomparison between these data-sets, and correlation with the broad K-T database available in the literature, indicate that the chemistry of the boundary clays is not globally constant. Variations are more common than similarities, both of absolute concentrations, and interelement ratios. For example, the chondrite normalized platinum-group elements (PGE) patterns of Stevns Klint are not like those of Woodside Creek, with the Pt/Os ratios showing the biggest variation. These differences in PGE patterns are difficult to explain by secondary alteration of a layer that was originally chemically homogeneous, especially for elements of such dubious crustal mobility as Os and Ir. The data also show that enhanced PGE concentrations, with similar trends to those of the boundary layers, occur in the Cretaceous sediments below the actual boundary at Stevns Klint and all three the New Zealand localities. This confirms the observations of others that the geochemistry of the boundary layers apparently does not record a unique component. It is suggested that terrestrial processes, eg. an extended period of Late Cretaceous volcanism can offer a satisfactory explanation for the features of the K-T geochemical anomaly. Such models would probably be more consistent with the observed stepwise, or gradual, palaeontological changes across this boundary, than the instant catastrophe predicated by the impact theory.

Tredous, M.↗

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗

In situ x-ray absorption investigations of a heterogenized molecular catalyst and its interaction with a carbon nanotube support

A highly active heterogenized molecular CO 2 reduction catalyst on a conductive carbon support is investigated to identify if its improved catalytic activity can be attributed to strong electronic interactions between catalyst and support. The molecular structure and electronic character of a [Re +1 (tBu-bpy)(CO) 3 Cl] (tBu-bpy = 4,4'-tert-butyl-2,2'-bipyridine) catalyst deposited on multiwalled carbon nanotubes are characterized using Re L 3 -edge x-ray absorption spectroscopy under electrochemical conditions and compared to the homogeneous catalyst. The Re oxidation state is characterized from the near-edge absorption region, while structural changes of the catalyst are assessed from the extended x-ray absorption fine structure under reducing conditions. Chloride ligand dissociation and a Re-centered reduction are both observed under applied reducing potential. The results confirm weak coupling of [Re(tBu-bpy)(CO) 3 Cl] with the support, since the supported catalyst exhibits the same oxidation changes as the homogeneous case. However, these results do not preclude strong interactions between a reduced catalyst intermediate and the support, preliminarily investigated here using quantum mechanical calculations. Thus, our results suggest that complicated linkage schemes and strong electronic interactions with the initial catalyst species are not required to improve the activity of heterogenized molecular catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Retrieval of Sulfur Dioxide Amount and Altitude from Spaceborne Hyperspectral UV Measurements: Theory and Application

We describe the physical processes by which a vertically localized absorber perturbs the top-of-atmosphere solar backscattered ultraviolet (UV) radiance. The distinct spectral responses to perturbations of an absorber in its column amount and layer altitude provide the basis for a practical satellite retrieval technique, the Extended Iterative Spectral Fitting (EISF) algorithm, for the simultaneous retrieval of these quantities of a SO2 plume. In addition, the EISF retrieval provides an improved UV aerosol index for quantifying the spectral contrast of apparent scene reflectance at the bottom of atmosphere bounded by the surface and/or cloud; hence it can be used for detection of the presence or absence of UV absorbing aerosols. We study the performance and characterize the uncertainties of the EISF algorithm using synthetic backscattered UV radiances, retrievals from which can be compared with those used in the simulation. Our findings indicate that the presence of aerosols (both absorbing and nonabsorbing) does not cause large errors in EISF retrievals under most observing conditions when they are located below the SO2 plume. The EISF retrievals assuming a homogeneous field of view can provide accurate column amounts for inhomogeneous scenes, but they always underestimate the plume altitudes. The EISF algorithm reduces systematic errors present in existing linear retrieval algorithms that use prescribed SO2 plume heights. Applying the EISF algorithm to Ozone Monitoring Instrument satellite observations of the recent Kasatochi volcanic eruption, we demonstrate the successful retrieval of effective plume altitude of volcanic SO2, and we also show the improvement in accuracy in the corresponding SO2 columns.

Yang, Kau↗

Investigation of the dilatational dissipation in compressible homogeneous shear flow

The dilatational dissipation rate within compressible homogeneous turbulent shear flow is studied using data from direct numerical simulations. It is found that the dilatational dissipation rate is mainly associated with large scale acoustic waves. Eddy shocklets are observed; however, they have little contribution to the average dissipation rate. A mechanism for the generation of eddy shocklets is shown to be the focusing of acoustic waves. Turbulence models for the dilatational dissipation rate are compared with data from the simulations. It is found that the formulation of Zeman used by Viegas & Rubesin to calculate a compressible mixing layer agrees well with the numerical simulation results for turbulent Mach numbers less than 0.3. However, it is also found that, for the Mach number range occurring in mixing layers, the model does not accurately represent the theory upon which it is based.

Blaisdell, G. A.↗

IUE observations of a hot DAO white dwarf: Implications for diffusion theory and photospheric stratification

Observations of the DAO white dwarf PG1210+533, including the first high dispersion spectrum of a hybrid H-He object of this nature were obtained by IUE. In contrast with hot DAs in the 50,000 K temperature range, PG1210+533 shows no narrow interstellar-like metal lines, in spite of an optically observed He/H abundance of 0.1. This lack of metal makes accretion from the ISM an unlikely source for the He in the PG1210+533 photosphere. A significant discovery in the high dispersion spectrum is the existence of a sharp, non-LTE like, core seen in the He II 1640 line. Such features are detected in DO white dwarfs. A small aperture SWP low dispersion observation reveals the Lyman alpha profile of PG1210+533 to be surprisingly weak and narrow. Fits of this profile using pure H models yielded a T(eff) = 56,000 K. Fits of the Balmer H gamma profile however, yield T(eff) = 42,300 K and log g = 8.5 + or - 0.5 for the same models. It is unlikely that homogeneously mixed H-He atmospheres can resolve the inconsistency between the Lyman alpha and H gamma features in this star. Stratified models involving thin H photospheres may be necessary to explain these results.

Holberg, J. B.↗

Atomic-resolution imaging as a mechanistic tool for studying single-site heterogeneous catalysis

Heterogeneous catalysts dominate the chemical industry but typically feature diverse, incompletely defined active sites. Thus, describing structure-activity relationships, unlike homogeneous catalysts, remains challenging. In contrast, molecularly defined single-site heterogeneous catalysts (SSHCs), using appropriate tools, are poised to address these challenges and provide new avenues for catalysis research and development. The present study explores eco-friendly H 2 production mediated by discrete MoO 2 sites supported on carbon nanohorns (CNHs) and active for alcohol dehydrogenation. While informative, detailed ensemble EXAFS/XANES/XPS, kinetic measurements, and DFT analysis alone cannot provide a full molecular picture of the reaction pathway. Here, using single-molecule atomic-resolution time-resolved electron microscopy (SMART-EM), we propose the identification of four key catalytic intermediates anchored to CNHs and uncover a new reaction pathway involving alkoxide/hemiacetal equilibration and acetal oligomerization. Furthermore, these intermediates are inferred through a combination of theory and SMART-EM, showcasing the potential of SMART-EM as a complementary tool for exploring mechanistic hypotheses in catalysis.

36 MATERIALS SCIENCE↗

Estimation of bare soil evaporation using multifrequency airborne SAR

It is shown that for homogeneous areas soil moisture can be derived from synthetic aperture radar (SAR) measurements, so that the use of microwave remote sensing can given realistic estimates of energy fluxes if coupled to a simple two-layer model repesenting the soil. The model simulates volumetric water content (Wg) using classical meterological data, provided that some of the soil thermal and hydraulic properties are known. Only four parameters are necessary: mean water content, thermal conductivity and diffusitivity, and soil resistance to evaporation. They may be derived if a minimal number of measured values of Wg and surface layer temperature (Tg) are available together with independent measurements of energy flux to compare with the estimated values. The estimated evaporation is shown to be realistic and in good agreement with drying stage theory in which the transfer of water in the soil is in vapor form.

Soares, Joao V.↗

Revealing the Configuration and Conformation of Surface Organometallic Catalysts with DNP-Enhanced NMR

Although single-site, supported organometallic catalysts were designed with homogeneity in mind, little is strictly known about their atomic-level structure and uniformity. This is in large part due to the inability of conventional characterization tools to provide structural information with adequate resolution and range for this purpose. In this work, we show that dynamic nuclear polarization (DNP)-enhanced solid-state nuclear magnetic resonance (NMR) enables the measurement of distances between the surface and individual carbons, making it possible to orient complexes and molecules on surfaces. We use this approach to determine the orientation and configuration of naturally 13 C-abundant acetate and supported organoiridium(III) pincer complex, both supported on γ-Al 2 O 3 . By combining solid-state NMR and extended X-ray absorption fine structure spectroscopy (EXAFS) experiments with periodic density functional theory (DFT) calculations, we are able to determine the three-dimensional arrangement of ligands around the metal center and thereby identify the structure of the supported complex at a resolution reminiscent of that associated with single-crystal diffractometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution of the age-included nearest pair distribution in disperse multiphase flows

The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermal Neutron Scattering Law Evaluation for Zirconium Carbide and Critical Mass Calculations

Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗