Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complexity 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 1,693 records · Page 94

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗

A Perspective on Multistep Pathways of Nucleation

While views on crystal nucleation throughout the 20th century were dominated by classical nucleation theory (CNT) in which the ordered nucleus forms directly from monomeric species, observations over the past two decades across scientific fields have revealed a plethora of multi-step nucleation processes, often involving more complex species. These processes have come to be known as non-classical, reflecting distinctions from the expectations of CNT. Here I consider the various classes of non-classical nucleation pathways and examine the extent to which they can be easily placed within the framework of CNT modified to include well-understood complexities in the free energy landscape and/or kinetic factors. I highlight recent examples that illustrate these pathways before turning to processes that cannot readily be included in this framework. I conclude by describing future research directions needed to fill the knowledge gaps that prevent the development of quantitative models for multi-step nucleation processes

De Yoreo, James J.↗

Crystal structure of deracoxib, C 17 H 14 F 3 N 3 O 3 S

The crystal structure of deracoxib has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Deracoxib crystallizes in space group Pbca (#61) with a = 9.68338(11), b = 9.50690(5), c = 38.2934(4) Å, V = 3525.25(3) Å 3 , and Z = 8. The molecules stack in layers parallel to the ab-plane. N–H∙∙∙O hydrogen bonds link the molecules along the b-axis, in chains with the graph set C1,1(4), as well as more-complex patterns. N–H∙∙∙N hydrogen bonds link the layers. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.

36 MATERIALS SCIENCE↗

Transfer Learning Meets Embedded Correlated Wavefunction Theory for Chemically Accurate Molecular Simulations: Application to Calcium Carbonate Ion Pairing

Achieving chemical accuracy for molecular simulations remains a central challenge in computational chemistry. Here, we present an embedded correlated wavefunction transfer learning (ECW-TL) framework for accurately simulating molecular dynamics in the condensed phase. ECW-TL incorporates high-level electron exchange and correlation effects in ECW theory while preserving the training and computational efficiency of machine-learned interatomic potentials. We demonstrate the framework on Ca 2+ –CO 3 2– ion pairing in aqueous solution, a key process underlying CO 2 mineralization in seawater. As proof of principle, we first show that fine-tuning a DFT-revPBE-D3(BJ) baseline model with embedded-DFT-SCAN data reproduces the DFT-SCAN free-energy surface within 1 kcal/mol across all solvation states. Extending the framework to embedded MP2 and localized natural-orbital CCSD(T) further refines the free-energy profile, revealing the crucial role of exact electron exchange and correlation in determining ion-pair stability and structure. The computed ion-pair association free energy is in quantitative agreement with experimental measurements, further validating the accuracy of the ECW-TL framework. ECW-TL thus provides a general, data-efficient route for transferring CW accuracy to efficient simulations of complex aqueous and interfacial chemical processes.

cluster chemistry↗

The origin of continental crust: Outlines of a general theory

The lower continental crust, formerly very poorly understood, has recently been investigated by various geological and geophysical techniques that are beginning to yield a generally agreed on though still vague model (Lowman, 1984). As typified by at least some exposed high grade terranes, such as the Scottish Scourian complex, the lower crust in areas not affected by Phanerozoic orogeny or crustal extension appears to consist of gently dipping granulite gneisses of intermediate bulk composition, formed from partly or largely supracrustal precursors. This model, to the degree that it is correct, has important implications for early crustal genesis and the origin of continental crust in general. Most important, it implies that except for areas of major overthrusting (which may of course be considerable) normal superposition relations prevail, and that since even the oldest exposed rocks are underlain by tens of kilometers of sial, true primordial crust may still survive in the lower crustal levels (of. Phinney, 1981).

Lowman, P. D., Jr.↗

Error and Uncertainty Quantification in the Numerical Simulation of Complex Fluid Flows

The failure of numerical simulation to predict physical reality is often a direct consequence of the compounding effects of numerical error arising from finite-dimensional approximation and physical model uncertainty resulting from inexact knowledge and/or statistical representation. In this topical lecture, we briefly review systematic theories for quantifying numerical errors and restricted forms of model uncertainty occurring in simulations of fluid flow. A goal of this lecture is to elucidate both positive and negative aspects of applying these theories to practical fluid flow problems. Finite-element and finite-volume calculations of subsonic and hypersonic fluid flow are presented to contrast the differing roles of numerical error and model uncertainty. for these problems.

Barth, Timothy J.↗

Deep level defects and cation sublattice disorder in ZnGeN 2

The III-nitrides have revolutionized lighting technology and power electronics. Expanding the nitride semiconductor family to include the heterovalent ternary nitrides opens up new and exciting opportunities for device design that may help overcome some of the limitations of the binary nitrides as well as new applications. However, the more complex cation sublattice also gives rise to new complexities with both native and extrinsic defects, such as disorder on the cation sublattice. Here, depth-resolved cathodoluminescence spectroscopy and surface photovoltage spectroscopy measurements of defect energy levels in ZnGeN 2 combined with transmission microscopy and crystal lattice diffraction reveal optical signatures of mid-gap states that can be associated with cation sublattice disorder. The energies of these characteristic optical signatures in ZnGeN 2 thin films grown by metal-organic chemical vapor deposition are in good agreement with energy levels predicted by density functional theory for specific thermodynamically stable native point defects. To test the dependence of these defect signatures on crystal atomic composition, we correlated spatially-resolved optical and atomic composition measurements using spatially resolved x-ray photoelectron spectroscopy with systematically varied growth conditions on the same ZnGeN 2 films. The resultant elemental maps versus defect spectral energies and intensities suggest that cation antisite complexes (Zn Ge -Ge Zn ) versus isolated native point defects form preferentially and introduce mid-gap states that dominate electron-hole pair recombination. Complexing of Zn Ge and Ge Zn antisites manifests as disorder in the cation sub-lattice and leads to the formation of wurtzitic ZnGeN 2 as indicated by transmission electron microscopy diffraction patterns and x-ray diffraction (XRD) reciprocal space maps. Furthermore, these findings emphasize the importance of growth conditions to control cation place exchange.

42 ENGINEERING↗

Space-Split Algorithm for Sensitivity Analysis of Discrete Chaotic Systems With Multidimensional Unstable Manifolds

Accurate approximations of the change of a system's output and its statistics with respect to the input are highly desired in computational dynamics. Ruelle's linear response theory provides breakthrough mathematical machinery for computing the linear response of chaotic dynamical systems. In this paper, we propose an algorithm for sensitivity analysis of discrete chaos with an arbitrary number of positive Lyapunov exponents. We combine the concept of perturbation space-splitting, which regularizes Ruelle's original expression, together with measure-based parameterization of the expanding subspace. We use these tools to rigorously derive trajectory-following recursive relations that converge exponentially fast, and construct a memory-efficient Monte Carlo scheme for derivatives of the output statistics. Thanks to the regularization and lack of simplifying assumptions on the system's behavior, our method is immune to the common problems of other popular methods such as the exploding tangent solutions and unphysical shadowing directions. Here, we provide a ready-to-use algorithm, analyze its complexity, and demonstrate several numerical examples of sensitivity computation using physically-inspired low-dimensional systems.

97 MATHEMATICS AND COMPUTING↗

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leveraging Intramolecular π-Stacking to Access an Exceptionally Long-Lived 3 MC Excited State in an Fe(II) Carbene Complex

The ability to manipulate excited-state decay cascades using molecular structure is essential to the application of abundant-metal photosensitizers and chromophores. Ligand design has yielded some spectacular results elongating charge-transfer excited state lifetimes of Fe(II) coordination complexes, but triplet metal-centered ( 3 MC) excited states - recently demonstrated to be critical to the photoactivity of isoelectronic Co(III) polypyridyls - have to date remained elusive, with temporally isolable examples limited to the picosecond regime. Here, with this report, we show how strong-field donors and intramolecular π-stacking can conspire to stabilize a long-lived 3 MC excited state for a remarkable 4.1 ± 0.3 ns in fluid solution at ambient temperature. Analysis of variable-temperature time-resolved absorption data using theoretical models ranging from Arrhenius to semiclassical Marcus theory, combined with computational modeling and X-ray crystallography, reveal a Jahn−Teller stabilized excited state with a high activation barrier for ground-state recovery. The net result is a chromophore with a 3 MC excited-state lifetime that is orders of magnitude longer than anything yet observed for an Fe(II) complex.

carbene compounds↗

Cellular Phosphate Sensing and Anion Binding by an Azacrown‐Calixpyrrole Hybrid

A hybrid receptor-sensor for anions originating from the merging of positively charged ammonium moieties for electrostatic attraction/stronger binding of azacrowns with directionality of calixpyrrole hydrogen bond donors for selectivity is investigated. As demonstrated this hybrid receptor-sensor shows a remarkable selectivity for orthophosphate even in the presence of other phosphates and anions found in cellular materials (K assoc H 2 PO 4 − >H 2 P 2 O 7 2− >AMP − ≫ADP 2− or ATP 3− over halides, nitrate, or hydrogen sulfate; all Na + salts in water) but also cellular polyphosphate or phospholipids. This selectivity is harnessed in a real-time monitoring of cell lysis by lysozyme, which releases orthophosphate and other phosphates and anions from the cells. This sensitive (LOD 0.4 μM) fluorescence-based microscale method compares favorably with the state-of-the-art techniques but can easily be practiced in a high-throughput screening (HTS) manner. The anion binding and selectivity in aqueous solutions were investigated by NMR and put in context with phosphate binding of the parent calix[4]pyrrole. The microscopic understanding of anion binding by the hybrid receptor was then obtained from a combination of density functional theory (DFT), classical molecular dynamics (MD) with explicit water solvation, and ab initio MD (AIMD) simulations. Correlating the NMR and fluorescence binding data with studies of solvation of the receptor, phosphate anion, and the resulting complex confirms the binding is largely driven by entropic component (TΔS) associated with receptor and anion desolvation.

Anions↗

Unique Aluminum Clusters Stabilized by Cation‐Ligand Cooperativity

Abstract Two novel cyclic Al 3 clusters, Li 2 [Al 3 (PR 2 ) 6 ] ⋅ 2(Et 2 O) [R=Ph ( 1 ), Cy ( 2 )] were synthesized through salt metathesis of metastable AlCl ⋅ (Et 2 O) n solutions with LiPR 2 . Both complexes were characterized by single crystal X‐Ray diffraction and their solid state and electronic structures are compared to previously reported cyclic Al 3 clusters. Compounds 1 and 2 were further characterized using density functional theory (DFT) methods to understand the nature of the free electron and the three‐center two‐electron Al−Al bonds. Calculations revealed that both clusters are influenced by pseudo Jahn‐Teller distortion and are likely in equilibrium between two low energy bonding states. The solid‐state structures of 1 and 2 are unique among Al 3 clusters, while the electronic structure is similar to the previously reported [Si t 4 Al 3 ]. In aluminum chemistry, 1 and 2 are rare examples of clusters containing isoelectronic cores in different ligand environments. The effects of the phosphide ligand sphere and strongly coordinated lithium atoms are discussed in detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Horizon constraints on holographic Green’s functions

We explore a new class of general properties of thermal holographic Green’s functions that can be deduced from the near-horizon behaviour of classical perturbations in asymptotically anti-de Sitter spacetimes. We show that at negative imaginary Matsubara frequencies and appropriate complex values of the wavenumber the retarded Green’s functions of generic operators are not uniquely defined, due to the lack of a unique ingoing solution for the bulk perturbations. From a boundary perspective these ‘pole-skipping’ points correspond to locations in the complex frequency and momentum planes at which a line of poles of the retarded Green’s function intersects with a line of zeroes. As a consequence the dispersion relations of collective modes in the boundary theory at energy scales ω ~ T are directly constrained by the bulk dynamics near the black-brane horizon. For the case of conserved U(1) current and energy-momentum tensor operators we give examples where the dispersion relations of hydrodynamic modes pass through a succession of poleskipping points as real wavenumber is increased. We discuss implications of our results for transport, hydrodynamics and quantum chaos in holographic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Discrete element modeling of granular hopper flow of irregular-shaped deformable particles

Many natural and engineered granular materials have relatively deformable particles. Besides particle size and shape, particle deformability is another salient factor that significantly impacts the material’s flow behavior. Here, in this work, the flow of irregular-shaped deformable particles in a wedge-shaped hopper is investigated using discrete element simulations. A bonded-sphere model is developed to simultaneously capture irregular particle shapes and particle-wise deformations (e.g., compression, deflection, and distortion). Quantitative analysis of the effects of irregular shapes and particle deformations shows that the increase in particle stiffness tends to increase initial packing porosity and decrease the flow rate in the hopper. Rigid particles tend to have clogging issues, whereas deformable particles have less chance to, indicating particle deformation reduces the critical bridging width in the hopper flow. Detailed analysis of stress fields is also conducted to provide insights into the mechanism of particle flow and clogging. Stresses and discharge rates calculated from numerical simulations are compared and show good agreement with Walker’s theory and the extended Beverloo formula. Simulations with various particle shape combinations are also performed and show that the initial packing porosity decreases with an increasing percentage of fibers while the discharge rate has a complex dependency on particle shapes.

09 BIOMASS FUELS↗

Integrating computation and experiment to investigate photoelectrodes for solar water splitting at the microscopic scale.

CONSPECTUS: Photoelectrochemical water-splitting is a promising and sustainable way to store the energy of the sun in chemical bonds and use it to produce hydrogen gas, a clean fuel. The key components in photoelectrochemical cells (PECs) are photoelectrodes, including a photocathode that reduces water to hydrogen gas and a photoanode that oxidizes water to oxygen gas. Materials used in photoelectrodes for PECs must effectively absorb sunlight, yield photogenerated carriers, and exhibit electronic properties that enable the efficient shuttling of carriers to the surface to participate in relevant water-splitting reactions. Discovering and understanding the key characteristics of optimal photoelectrode materials is paramount to the realization of PEC technologies. Oxide-based photoelectrodes can satisfy many of these materials requirements, including stability in aqueous environments, band edges with reasonable alignment with the redox potentials for water splitting, and ease of synthesis. However, oxide photoelectrodes generally suffer from poor charge transport properties and considerable bulk electron-hole separation, and they have relatively large band gaps. Numerous strategies have been proposed to improve these aspects and understand how these improvements are reflected in the photoelectrochemical performance. Unfortunately, the structural and compositional complexity of multinary oxides accompanied by the inherent complexity of photoelectrochemical processes makes it challenging to understand the individual effects of composition, structure, and defects in the bulk and on the surface on a material's photoelectrochemical properties. The integration of experiment and theory has great potential to increase our atomic-level understanding of structure-composition-property relationships in oxide photoelectrodes. In this Account, we describe how integrating experiment and theory is beneficial for achieving scientific insights at the microscopic scale. We highlight studies focused on understanding the role of (i) bulk composition via solid-state solutions, intercalation, and comparison with isoelectronic compounds, (ii) dopants for both the anion and cation and their interactions with oxygen vacancies, and (iii) surface/interface structure in the photocurrent generation and photoelectrochemical performance in oxide photoelectrodes. In each instance, we outline strategies and considerations for integrating experiment and theory and describe how this integration led to valuable insights and new directions in uncovering structure-composition-property relationships. Our aim is to demonstrate the unique value of combining experiment and theory in studying photoelectrodes and to encourage the continued effort to bring experiment and theory in closer step with each other.

Wang, Wennie↗

First-Principles Insights into Proton-Coupled Electron Transfer versus Hydrogen Evolution Reaction Selectivity from a Base-Appended Cobaltocene Mediator

Performing selective proton-coupled electron transfer (PCET) to substrates such as N 2 , CO 2 , and unsaturated organic molecules under electrochemical conditions requires the suppression of the competing hydrogen evolution reaction (HER). To address this challenge, our laboratory previously demonstrated a PCET mediator strategy using a dimethylaniline-appended cobaltocene complex, [(CpCoCp NMe2 )H] + , which performs selective reductive chemistry while suppressing the HER. However, the origin of the suppressed, yet still observable, HER has not been thoroughly established. In this work, we perform density functional theory (DFT) calculations to elucidate the HER mechanism involving this redox mediator and to provide atomistic insights into the bifurcation between the PCET and HER pathways. We find that protonation of the aniline moiety to form [CpCoCp NMe2H ] + is more favorable, both kinetically and thermodynamically, than formation of the ring-protonated species [(CpCo(Cp-H) NMe2 )] + . Furthermore, PCET to acetophenone is energetically more favorable via [CpCoCp NMe2H ] + than via [(CpCo(Cp-H) NMe2 )] +1/0 . In contrast, the most favorable HER pathway involves the ring-protonated Co(I) species. These results offer mechanistic insights into HER versus PCET bifurcation and establish guiding principles for designing PCET mediators for selective electroreductive transformations.

evolution reactions↗