Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 433 records · Page 24

High-rate deposition of ultrathick bismuth films by direct-current magnetron sputtering

Bismuth coatings deposited at room temperature by physical vapor deposition typically feature irregularly-shaped grains and a high density of voids, thereby preventing their use in many potential applications. Here, we demonstrate the deposition of low-stress, high-density 40-μm-thick Bi films by systematically studying the effect of the deposition rate and substrate tilt on the properties of Bi films deposited at room temperature by direct-current magnetron sputtering. Furthermore, results show that, with increasing deposition rate, the morphology, structure, and properties of Bi films greatly improve, while an increase in substrate tilt leads to a monotonic decrease in film thickness, density, and electrical conductivity.

Bismuth↗

Dynamics and bifurcations of laminar annular swirling and non-swirling jets

This paper presents bifurcation analyses characterising the nonlinear dynamics of fully developed laminar annular jets with respect to the centrebody diameter ${ {d}}$ , Reynolds number ${ {Re}}$ , and swirl ratio ${ {S}}$ . Similar flows appear in numerous applications and feature a vibrant range of topological and dynamical characteristics associated with phenomena including shear layer separation and vortex breakdown. Our results begin by describing the non-monotonic evolution of the axisymmetric jet's steady topology under varying ${ {S}}$ . In accord with earlier reports, the jet progresses through a sequence of wake, breakdown and wall jet regimes in a qualitatively similar manner across a wide span of ${ {d}}$ and ${ {Re}}$ values. In the wake regime, the non-swirling jet bifurcates to a plane-symmetric, but not axisymmetric, steady flow pattern beyond a ${ {d}}$ -dependent critical ${ {Re}}$ value. With further increase in ${ {Re}}$ , the steady non-swirling jet destabilises subsequently via multiple distinct Hopf bifurcations. Introducing ${ {S}}>0$ to the jet also induces unsteadiness by twisting the singly azimuthally periodic ( $|m|=1$ ) asymmetric wake structure and causing it to precess periodically in time about the central axis. Intermediate swirl stabilises this unsteady dynamics and restores the jet's axisymmetry. This stabilising effect is then reversed in the breakdown regime at higher ${ {S}}$ , where a variety of different $|m|=1$ and $|m|=2$ instabilities bifurcate from the steady flow as ${ {S}}$ is increased. Several instances of hysteresis and subcritical behaviour are reported and discussed, including one that manifests precessing vortex core oscillations.

Douglas, Christopher M. (ORCID:0000000259683315)↗

Effects of thermophoresis on high-pressure binary-species boundary layers with uniform and non-uniform compositions

Direct numerical simulations of high-pressure binary-species temporal boundary layers are performed to investigate the flow physics for three situations: (1) uniform and equal composition, (2) uniform but unequal compositions and (3) non-uniform composition. Both colder- and hotter-wall situations compared with the free stream are simulated. The working fluid is a nitrogen/methane mixture. The analysis is performed at a case-specific self-similar state. Even when the initial composition is uniform, the methane mean mass fraction decreases near the colder wall, whereas it increases near the hotter wall and the mass fraction fluctuates in the entire boundary layer. Analysis of the species-mass diffusion balance and flow structures reveal that both mass-fraction variation and fluctuations are induced by the Soret effect. When the initial composition is non-uniform and the wall is colder, the methane mean mass fraction monotonically increases from the wall akin to its initial profile. However, when the wall is hotter the mean mass fraction decreases near the wall in contrast to its initial profile, a fact traced through the species-mass diffusion balance to the Soret effect being large and enriching methane near the wall. In contrast, the direction of the Soret flux is opposite for the colder wall, thus keeping the methane concentration small. Although the initial magnitude of the difference between the wall and free-stream temperature is the same in all cases, the situation is not symmetric between colder- and hotter-wall cases; the flow structure exhibits much smaller scales when the wall is hotter than when the wall is colder.

Mechanics↗

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗

Plasma potential shaping using end-electrodes in the Large Plasma Device

We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗

Comparison of the Forward and Reverse Photocycle Dynamics of Two Highly Similar Canonical Red/Green Cyanobacteriochromes Reveals Unexpected Differences

Cyanobacteriochromes (CBCRs) are cyanobacterial photoreceptors that exhibit photochromism between two states: a thermally stable dark-adapted state and a metastable light-adapted state with bound linear tetrapyrrole (bilin) chromophores possessing 15Z and 15E configurations, respectively. The photodynamics of canonical red/green CBCRs have been extensively studied; however, the time scales of their excited-state lifetimes and subsequent ground-state evolution rates widely differ and, at present, remain difficult to predict. Here, we compare the photodynamics of two closely related red/green CBCRs that have substantial sequence identity (~68%) and similar chromophore environments: AnPixJg2 from Anabaena sp. PCC 7120 and NpR6012g4 from Nostoc punctiforme. Using broadband transient absorption spectroscopy on the primary (125 fs to 7 ns) and secondary (7 ns to 10 ms) time scales together with global analysis modeling, our studies revealed that AnPixJg2 and NpR6012g4 have comparable quantum yields for initiating the forward ( 15Z P r → 15E P g ) and reverse ( 15E P g → 15Z P r ) reactions, which proceed through monotonic and nonmonotonic mechanisms, respectively. In addition to small discrepancies in the kinetics, the secondary reverse dynamics resolved unique features for each domain: intermediate shunts in NpR6012g4 and a Meta-Gf intermediate red-shifted from the 15Z P r photoproduct in AnPixJg2. Altogether, this study supports the conclusion that sequence similarity is a useful criterion for predicting pathways of the light-induced evolution and quantum yield of generating primary intermediate Φ p within subfamilies of CBCRs, but more studies are still needed to develop a comprehensive molecular level understanding of these processes.

59 BASIC BIOLOGICAL SCIENCES↗

Pressure-Induced Superconductivity in the Wide-Band-Gap Semiconductor Cu 2 Br 2 Se 6 with a Robust Framework

We report pressure-induced superconductivity in a ternary and nonmagnetic Cu-containing semiconductor, Cu 2 Br 2 Se 6 , with a wide band gap of 1.89 eV, in which the Cu and Br atoms generate infinite 2, helical chains along the c-axis and are linked by the cyclohexane-like Se 6 rings to form a three-dimensional framework. We find that this framework is remarkably robust under compression, and the ambient-pressure phase survives at least to our experimental limit of 32.1 GPa. Concurrent semiconductor-to-metal transition and superconductivity are observed above 21.0 GPa. The superconducting temperature monotonically increases from 4.0 to 6.7 K at 40.0 GPa. Furthermore, first-principles calculations show that the emergence of superconductivity is associated with the formation of weak multicentered bonds that involve the increase in coordination of the Cu atoms and a subset of the Se atoms. The observation of superconductivity in this type of nonmagnetic transition-metal-based material will inspire the exploration of related new superconductors under pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Influence of Annealing and Composition on the Crystal Structure of Mixed-Halide, Ruddlesden–Popper Perovskites

Mixed-halide two-dimensional (2D) hybrid organic–inorganic perovskites offer an important opportunity to control the band gap for applications in optoelectronic devices. In this study, we focus on phenethylammonium lead halide [(PEA) 2 Pb(I 1– x Br x ) 4 ] films, the pure iodide form of which is one of the most widely studied optically active 2D perovskite systems. Resonant infrared, matrix-assisted pulsed laser evaporation is used to grow films to explore the effects of post-growth annealing and mixed-halide composition. The composition, crystal structure, and optical properties are studied for as-grown and annealed films of (PEA) 2 Pb(I 1– x Br x ) 4 for x = 0, 0.25, 0.5, 0.75, and 1. First-principles calculations are used in conjunction with the experimental data to explain the mixed-halide behavior that does not trend monotonically with the bromide content. Important results of this work are as follows: (i) X-ray diffraction reveals evidence for halide phase separation around x = 0.25, consistent with first-principles calculations, whereas no phase separation is observed for x = 0.5 and above and (ii) a unique photoluminescence (PL) peak splitting is observed for x = 0.75, a composition for which no compositional phase separation is observed. The PL splitting is tentatively explained by the coexistence of two distinct types of halide short-range ordering at x = 0.75. Overall, this study demonstrates that mixed-halide n = 1 Ruddlesden–Popper perovskites are not simple random alloys but that instead, they display distinct sites and ordering preferences of the different halide anions. These preferences are critical to understand and rationally tune the properties of the materials.

36 MATERIALS SCIENCE↗

Correlated Ion Transport Governed by Dynamic Local Structure in High Concentration and Localized High Concentration Electrolytes

Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.

Mohanakrishnan, Rohith Srinivaas [University of Ca↗

A Simple Criterion for Feasibility of Heat Integration between Distillation Streams Based on Relative Volatilities

In a multicomponent distillation configuration, there are numerous sources and sinks of heat, and a potential way to reduce the heat duty requirement is to perform heat integration. Unfortunately, an algorithmic search of the optimal heat integration opportunities is intractable when the required temperatures of intermediate streams are computed via complex models. Instead, in this work, we introduce pressure-scaled pseudo relative volatility, a new metric to compare stream temperatures. Here, we justify the use of pseudo relative volatility by proving that this variable is a monotonically increasing function of the liquid fraction in a saturated mixture stream. Using this metric, we derive a shortcut criterion to check the feasibility of various heat integration opportunities, such as thermal coupling via heat transfer (TCH). The advantage of this approach is that it circumvents the need for explicit temperatures and instead relies on composition, component relative volatilities, and pressure—quantities that are readily available in shortcut models for optimization of distillation configurations. Leveraging this fact, we propose a new optimization framework that identifies feasible TCHs that we consider within the formulation while minimizing the total heat duty of a distillation configuration. We demonstrate, on a few examples, that our formulation can identify heat duty efficient configurations, some of which are multieffect configurations. Using this methodology, we discover configurations that are not only simpler than the fully thermally coupled (FTC) configuration but also have a much lower heat duty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields↗

The many-body expansion for aqueous systems revisited: 1. Water-water interactions

We revisit the Many-Body Expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the Basis Set Superposition Error (BSSE) correction, on the various terms for selected sizes of water clusters up to n = 21. The analysis is performed at the second order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality for the (H2O)n, n = 7, 10, 13, 16 and 21 clusters, for which we report either the complete MBE (n = 7, 10) or the one through the 6-body (n = 13) and the 5-body terms (n = 16, 21). Our results suggest that any sizeable contributions to the total cluster binding energy arising from the 5-body and larger terms are solely an artifact of the size of the basis set. Indeed, all terms above the 4-body converge monotonically to practically zero at the Complete Basis Set (CBS) limit. The magnitudes of the 5-body and above terms are accurately reproduced even with the smaller basis set of the series (aug-cc-pVDZ) once the BSSE correction is taken into account. The same level of theory (MP2/aug-cc-pVDZ, BSSEcorrected) also accurately reproduces the MP2/CBS values of the 3- and 4-body terms. The contribution of electron correlation to the 3- and 4-body terms is quite small, so that neglecting the correlation contribution in all terms above the 3-body results in an error of ??(0.1%). The BSSE correction to the largest 2-body term in the MBE was accurately estimated from the function ?? [1 + erf(-?? ??)], which is proportional to the common (overlapping) area between two Gaussian distributions whose centers are separated by R, with the constants ??, ?? fitted to the calculated BSSE corrections for the individual 2-body terms of the clusters with each basis set. Our results unquestionably demonstrate that the MBE for water-water interactions converges to the 4-body term, since any finite terms above the 4-body are artifacts of the size of the basis set and can be safely truncated at the 4-body term when either a very large basis set is used or BSSE corrections are taken into account even with the aug-cc-pVDZ basis set. We expect these findings to have important consequences in the pursuit of accurate many-body molecular dynamics simulations for aqueous systems.

Heindel, Joseph↗

Effects of Field Strength and Silver Nanowire Size on Plasmon-Enhanced N 2 Dissociation

Dissociation of the nitrogen molecule via plasmon-enhanced catalysis using noble metal nanoparticles has been investigated both experimentally and computationally in recent years. However, the mechanism of plasmon-enhanced nitrogen dissociation is still not very clear. In this work, we apply theoretical approaches to examine the dissociation of a nitrogen molecule on atomically-thin Ag n nanowires (n = 6, 8, 10, 12) and a Ag 19 + nanorod. Ehrenfest dynamics provides information about the motion of nuclei during the dynamics process and real-time TDDFT calculations show the electronic transitions and population of electrons over the first 10s of fs time scale. The activation and dissociation of nitrogen is typically enhanced when the electric field strength increases. However, the enhancement is not always monotonic with field strength. As the length of the Ag wire increases, nitrogen is typically easier to dissociate and thus requires lower field strengths, even though the plasmon frequency is lower. The Ag 19 + nanorod leads to faster dissociation of N 2 than the atomically-thin nanowires. Altogether, our detailed study yields insights into the mechanisms involved in plasmon-enhanced N 2 dissociation, as well as provides information about factors that can be used to improve adsorbate activation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗

Hydrogen Bond Network Disruption by Hydration Layers in Water Solutions with Salt and Hydrogen-Bonding Polymers (PEO)

A mean field theory model describing the interaction of ion hydration layers with the network of hydrogen bonds of both water and the nonionic polymer poly(ethylene oxide) (PEO) is presented. The predictions of the model for types and statistics of hydrogen bonds, the number of water molecules bound to PEO, or their dependence on temperature are successfully verified from all-atom simulations at different NaCl and PEO concentrations. Furthermore, our simulations show that the binding of cations to PEO increases monotonically with salt concentration, in agreement with recent experimental results, through a mechanism in which the sum of the number of bound water and cations is independent of salt concentration. Here, the model introduced is general and can describe any salt or hydrogen-bond-forming polymer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetics of Shear-Induced Structural Ordering in Dense Colloids

The macroscopic rheological response of a colloidal solution is highly correlated with the local microscopic structure, as revealed by an in situ Rheo-SAXS experiment with a high temporal resolution. Oscillatory shear can induce a strain-controlled ordering-to-disorder transition, resulting in a shear-thickening process that is different from the normal shear-thickening behavior that is driven by hydrodynamics and particle friction. Here, we reveal that there is a complex time-dependent kinetics toward structural ordering under different applied strains. When the strain amplitude reaches a critical value that starts to induce disordering in the system, the pathway toward the dynamic equilibrium can also become highly non-monotonic. Within the same oscillatory cycle, there is a strong correlation of ordering with different phases of the oscillation, with the system oscillating between two dynamic metastable states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗