Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex computing”

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 649 records · Page 36

How generic is iodide-tagging photoelectron spectroscopy: An extended investigation on the Gly· X - (Gly = glycine, X = Cl or Br) complexes

We report a joint negative ion photoelectron spectroscopy (NIPES) and quantum chemical computational study on glycine-chloride/bromide complexes (denoted Gly·X - , X = Cl/Br) in close comparison to the previously studied Gly·I - cluster ion. Combining experimental NIPE spectra and theoretical calculations, various Gly·X - complexes were found to adopt the same types of low-lying isomers, albeit with different relative energies. Despite more congested spectral profiles for Gly·Cl - and Gly·Br - , spectral assignments were accomplished with the guidance of the knowledge learned from Gly·I - , where a larger spin–orbit splitting of iodine afforded well-resolved, recognizable spectral peaks. Three canonical plus one zwitterionic isomer for Gly·Cl - and four canonical conformers for Gly·Br - were experimentally identified and characterized in contrast to the five canonical ones observed for Gly·I - under similar experimental conditions. Taken together, this study investigates both genericity and variations in binding patterns for the complexes composed of glycine and various halides, demonstrating that iodide-tagging is an effective spectroscopic means to unravel diverse ion-molecule binding motifs for cluster anions with congested spectral bands by substituting the respective ion with iodide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rigorous computation of short-range order unifies its controversial effects in complex concentrated alloys

Direct experimental observations of chemical short-range order (SRO) in complex concentrated alloys (CCAs) have triggered high interest. However, the reported effects of SRO on yield stresses are controversial, and their atomic-scale mechanisms are elusive, which limits our ability to utilize SRO in alloy design. Here we tackle this challenge using an advanced computational approach that rigorously takes into account the critical lattice distortion in CCAs and further verify our theoretical predictions with experiments. We show that the CoCrNi model alloy has a narrow temperature window around 670 °C for SRO formation. This explains why the mechanical effect of SRO is observed in some experiments but not in others. Here, we propose an effective alloy-doping method to control SRO and reveal atomic-bonding types that dominate SRO formation for different alloys. The strategies and insights generally apply to a broad spectrum of alloys, laying the foundation for designing advanced alloys by manipulating their SRO.

36 MATERIALS SCIENCE↗

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

Negative complexity of formation: the compact dimensions strike back

We show that the vacuum-subtracted maximal volume, the proposed holographic dual to complexity of formation, can be negative when contributions from compact directions are included. We construct explicit solutions with arbitrarily negative complexity of formation in asymptotically AdS 4 × S 7 SUGRA. These examples rely critically on the compact directions, specifically the fact that the full eleven-dimensional spacetime is not asymptotically AdS 11 . While there is some ambiguity in the extension of the holographic complexity proposal to the compact directions, we show that the two natural candidates can both have arbitrarily negative complexity of formation in SUGRA solutions. We further find examples in which complexity can even decrease at late times, including cases of both single-sided geometries and two-sided wormholes. In particular, we construct a cosmological wormhole with simultaneously negative and decreasing complexity of formation (as computed by volume) at late times. We find a distinguished role for relevant primaries in these constructions and comment on possible interpretations.

2D Gravity↗

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗

Clustering and Cliques in Preferential Attachment Random Graphs with Edge Insertion

In this paper, we investigate the global clustering coefficient (a.k.a transitivity) and clique number of graphs generated by a preferential attachment random graph model with an additional feature of allowing edge connections between existing vertices. Specifically, at each time step t, either a new vertex is added with probability f(t), or an edge is added between two existing vertices with probability 1 – f(t). We establish concentration inequalities for the global clustering and clique number of the resulting graphs under the assumption that f(t) is a regularly varying function at infinity with index of regular variation –$\gamma$, where $\gamma$ $\in$ [0, 1). Finally, we also demonstrate an inverse relation between these two statistics: the clique number is essentially the reciprocal of the global clustering coefficient.

97 MATHEMATICS AND COMPUTING↗

Effects of Aluminum and Molybdenum on the Phase Stability of Iron-Chromium Alloys: A First-Principles Study

The interaction between solute atoms is critical to the thermodynamic behavior of Fe-Cr alloys, but the effects of non-dilute Al and Mo on the Fe-Cr phase stability and vacancy formation energy are not clearly understood. In this study, density functional theory, cluster expansion, and Monte Carlo simulation are used to predict the effects of ternary solute elements on the thermodynamic properties in multicomponent Fe-Cr alloys. The machine learning regression approach is applied to train and construct energy models that accurately describe the complex chemical interactions. The computational outcomes include the prediction of the partial ternary phase diagram, mixing enthalpy, and vacancy formation energy for different compositions. The phase boundary calculation predicts a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from a’ precipitates. The mixing enthalpy calculation shows that Mo may also reduce the Cr solubility in bcc Fe. Additionally, the simulations show that non-dilute Cr decreases the vacancy formation energy in bcc Fe, while adding Al results in a less significant effect. Finally, the results demonstrate important applications of using machine learning energy models to study model or commercial alloys with multicomponent solute species and point defects.

36 MATERIALS SCIENCE↗

Neural message-passing for objective-based uncertainty quantification and optimal experimental design

Various real-world scientific applications involve the mathematical modeling of complex uncertain systems with numerous unknown parameters. Accurate parameter estimation is often practically infeasible in such systems, as the available training data may be insufficient and the cost of acquiring additional data may be high. In such cases, based on a Bayesian paradigm, we can design robust operators retaining the best overall performance across all possible models and design optimal experiments that can effectively reduce uncertainty to enhance the performance of such operators maximally. While objective-based uncertainty quantification (objective-UQ) based on MOCU (mean objective cost of uncertainty) provides an effective means for quantifying uncertainty in complex systems, the high computational cost of estimating MOCU has been a challenge in applying it to real-world scientific/engineering problems. In this work, we propose a novel scheme to reduce the computational cost for objective-UQ via MOCU based on a data-driven approach. We adopt a neural message-passing model for surrogate modeling, incorporating a novel axiomatic constraint loss that penalizes an increase in the estimated system uncertainty. As an illustrative example, we consider the optimal experimental design (OED) problem for uncertain Kuramoto models, where the goal is to predict the experiments that can most effectively enhance robust synchronization performance through uncertainty reduction. We show that our proposed approach can accelerate MOCU-based OED by four to five orders of magnitude, without any visible performance loss compared to the state-of-the-art. The proposed approach applies to general OED tasks, beyond the Kuramoto model.

97 MATHEMATICS AND COMPUTING↗

Surrogate modeling of Monte Carlo radiation transport with convolutional neural networks for shielding optimization

Here, we present a machine learning (ML)-based surrogate model using convolutional neural networks (CNN) designed to emulate the attenuation of neutron fields as they pass through various shielding materials. This model can compute the outgoing neutron flux almost instantaneously and achieves reasonable accuracy compared to traditional Monte Carlo (MC)-based codes, which are computationally intensive. This emulator alleviates the complexity of neutron radiation transport through shielding materials by reducing the dimensionality and enables shielding optimization for a known radiation environment. This optimization process, which would have taken an unrealistic timeline due to several complex radiation transport simulations, can now be achieved in minutes, thus increasing computational capabilities in radiation shielding assessment. We demonstrate the applications of this emulator in computing effective dose rates and optimizing shielding solutions for a heavy-ion accelerator facility, such as the Facility for Rare Isotope Beams, where secondary neutrons produced via beam interactions dominate the radiation environment.

accelerator shielding↗

Urban building energy modeling (UBEM) tools: A state-of-the-art review of bottom-up physics-based approaches

Regulations corroborate the importance of retrofitting existing building stocks or constructing new energy-efficient districts. There is, thus, a need for modeling tools to evaluate energy scenarios to better manage and design cities, and numerous methodologies and tools have been developed. Among them, Urban Building Energy Modelling (UBEM) tools allow the energy simulation of buildings at large scales. Choosing an appropriate UBEM tool, balancing the level of complexity, accuracy, usability, and computing needs, remains a challenge for users. The review focuses on the main bottom-up physics-based UBEM tools, comparing them from a user-oriented perspective. Five categories are used: (i) the required inputs, (ii) the reported outputs, (iii) the exploited workflow, (iv) the applicability of each tool, and (v) the potential users. Moreover, a critical discussion is proposed, focusing on interests and trends in research and development. The results highlighted major differences between UBEM tools that must be considered to choose the proper one for an application. Finally, barriers of adoption of UBEM tools include the needs of a standardized ontology, a common three-dimensional city model, a standard procedure to collect data, and a standard set of test cases. This feeds into future development of UBEM tools to support cities’ sustainability goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Resonant X-ray Diffraction Reveals the Location of Counterions in Doped Organic Mixed Ionic Conductors

Organic mixed ionic–electronic conductors (OMIECs) have the potential to enable diverse new technologies, ranging from novel in situ biosensors to flexible energy storage devices and neuromorphic computing platforms. However, their complex behavior in functional films involving electrolyte-induced swelling, ion ingress, and electrochemical doping inhibits rational material design. Of critical importance is an understanding of the specific location of the ions in the volumetrically doped material, yet this information is not readily available. Here, in this report, we present the use of grazing- incidence resonant X-ray diffraction (RXRD, also known as anomalous diffraction) at S and Cl K-edges to determine the structure of a doped, prototypical, semicrystalline polymer OMIEC based on oligo(ethylene glycol) substitution of regioregular polythiophene. The RXRD measurement provides two key insights. Quantitative analysis of the RXRD allows the determination of the position of the ion relative to the polymer backbone in the crystalline regions. We find that the anion is relatively distant from the backbone, nearer to the lamella mid-plane naively in conflict with expected Coulombic attraction between the ion and the doped polymer polaron. Comparison of RXRD to Cl – fluorescence (total Cl – ) allows determination of the relative order of doping between the crystalline and amorphous regions. We find preferential doping of the crystalline regions. Both insights, the preferential doping of crystals at low potential and the specific location of the counterion with respect to the polymer backbone, are critical to developing a microscopic understanding of transport in OMIECs.

36 MATERIALS SCIENCE↗

Modeling Oxidative Dehydrogenation of Propane with Supported Vanadia Catalysts Using Multireference Methods

The oxidative dehydrogenation of propane over supported vanadium oxide catalysts poses significant computational challenges due to complex electronic structure changes along the reaction coordinate, driven primarily by changes in the oxidation states of vanadium. To address these challenges, we systematically test quantum chemical methods, including multireference (MR) approaches, domain-based local pair natural orbital coupled cluster theory (DLPNO-CCSD(T)), and density functional theory (DFT). The initial C–H bond-breaking transition state requires MR treatment due to its multireference character, while subsequent steps permit efficient single-reference calculations. For the rate-limiting C–H activation step mediated by the vanadyl moiety, complete 1 active space second-order perturbation theory (CASPT2) yields an apparent activation barrier (E app 600K ) of 138 kJ/mol, consistent with experimental values (134 ± 4 kJ/mol; Gruene et al. Catal. Today 2010, 157, 137). In contrast, DLPNO-CCSD(T) overestimates this barrier (198 kJ/mol), whereas DFT predictions span 125–150 kJ/mol, depending on the functional. Our multireference investigation of this transition metal oxide-catalyzed process demonstrates that an active space that incorporates the C–H σ and V=O σ/π bonding orbitals, oxygen lone pairs, and their antibonding counterparts adequately captures electronic structure changes along the chemical transformation. Furthermore, these findings provide a general strategy for active space selection in transition metal oxide-catalyzed C/O–H bond activation reactions. The reference dataset from this work, which includes MR calculations with manually selected active spaces for all intermediates and transition states in the propane ODH reaction network, will serve as a benchmark for automating active space selection in similar systems.

Catalysts↗

Dinuclear Gold(I) Complexes Bearing Alkyl-Bridged Bis(N-heterocyclic carbene) Ligands as Catalysts for Carboxylative Cyclization of Propargylamine: Synthesis, Structure, and Kinetic and Mechanistic Comparison to the Mononuclear Complex [Au(IPr)Cl]

Eight new dinuclear gold(I) complexes, [Au 2 (L)X 2 ] (1–8), were synthesized using a straightforward synthetic procedure under very mild conditions. The complexes have been characterized by NMR spectroscopy, elemental analysis, and single-crystal X-ray structure analysis. Their catalytic activity was investigated in the carboxylative cyclization of propargylamine (PPA). A superior performance in comparison to [Au(IPr)Cl] (9) was obtained for complexes 1 and 2 having an eight-methylene bridge connecting two NHCs with an arene bearing an isopropyl substituent for X = Cl, Br. This prompted more detailed kinetic and mechanistic studies by FTIR comparing dinuclear complex 2 of X = Cl to complex 9. Fortuitously the FTIR studies allowed monitoring of the formation of the products carbamic acid (CA) and carbamate salt (CS), as well as a key cyclized intermediate first discovered by Ikariya. These data allow additional insight into the mechanism as well as the central role which may be played by Au(I) carbamate formation as a higher energy resting state present in the catalytic cycle. In conclusion, the crystal structures of four of the new complexes and a detailed computational study relevant to the role of carbamic acid (CA) and carbamates in the catalytic cycle are also reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ALD-Derived WO 3– x Leads to Nearly Wake-Up-Free Ferroelectric Hf 0.5 Zr 0.5 O 2 at Elevated Temperatures

Breaking the memory wall in advanced computing architectures will require complex 3D integration of emerging memory materials such as ferroelectrics─either within the back-end-of-line (BEOL) of CMOS front-end processes or through advanced 3D packaging technologies. Achieving this integration demands that memory materials exhibit high thermal resilience, with the capability to operate reliably at elevated temperatures, such as 125°C, due to the substantial heat generated by front-end transistors. However, silicon-compatible HfO 2 -based ferroelectrics tend to exhibit antiferroelectric-like behavior in this temperature range, accompanied by a more pronounced wake-up effect, posing significant challenges to their thermal reliability. Here, we report that by introducing a thin tungsten oxide (WO 3–x ) layer─known as an oxygen reservoir─and carefully tuning its oxygen content, ultrathin Hf 0.5 Zr 0.5 O 2 (5 nm) films can be made robust against the ferroelectric-to-antiferroelectric transition at elevated temperatures. This approach not only minimizes polarization loss in the pristine state but also effectively suppresses the wake-up effect, reducing the required wake-up cycles from 10 5 to only 10 at 125°C, a qualifying temperature for back-end memory integrated with front-end logic, as defined by the JEDEC standard. First-principles density functional theory (DFT) calculations reveal that WO 3 enhances the stability of the ferroelectric orthorhombic phase (o-phase) at elevated temperatures by increasing the tetragonal-to-orthorhombic phase energy gap and promoting favorable phonon mode evolution, thereby supporting o-phase formation under both thermodynamic and kinetic constraints.

36 MATERIALS SCIENCE↗

How Equivalent Are Equivalent Porous Media?

Abstract Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

58 GEOSCIENCES↗

A data-driven approach to modeling cancer cell mechanics during microcirculatory transport

In order to understand the effect of cellular level features on the transport of circulating cancer cells in the microcirculation, there has been an increasing reliance on high-resolution in silico models. Accurate simulation of cancer cells flowing with blood cells requires resolving cellular-scale interactions in 3D, which is a significant computational undertaking warranting a cancer cell model that is both computationally efficient yet sufficiently complex to capture relevant behavior. Given that the characteristics of metastatic spread are known to depend on cancer type, it is crucial to account for mechanistic behavior representative of a specific cancer’s cells. To address this gap, in the present work we develop and validate a means by which an efficient and popular membrane model-based approach can be used to simulate deformable cancer cells and reproduce experimental data from specific cell lines. Here, cells are modeled using the immersed boundary method (IBM) within a lattice Boltzmann method (LBM) fluid solver, and the finite element method (FEM) is used to model cell membrane resistance to deformation. Through detailed comparisons with experiments, we (i) validate this model to represent cancer cells undergoing large deformation, (ii) outline a systematic approach to parameterize different cell lines to optimally fit experimental data over a range of deformations, and (iii) provide new insight into nucleated vs. non-nucleated cell models and their ability to match experiments. While many works have used the membrane-model based method employed here to model generic cancer cells, no quantitative comparisons with experiments exist in the literature for specific cell lines undergoing large deformation. Here, we describe a phenomenological, data-driven approach that can not only yield good agreement for large deformations, but explicitly detail how it can be used to represent different cancer cell lines. This model is readily incorporated into cell-resolved hemodynamic transport simulations, and thus offers significant potential to complement experiments towards providing new insights into various aspects of cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗

Preliminary Introduction of a Free Vortex Wake Method Into OpenFAST

Increasing wind turbine size can have substantial economic benefits. However, increasing wind turbine turbine size requires the use of models that accurately capture the aerodynamic response to large deflections and nonstraight blade geometries. Additionally, the installation of offshore floating turbines requires turbine models that can handle large rotor motion. Free vortex wake (FVW) methods are capable of modeling this type of complex physics while remaining computationally tractable to perform the many simulations necessary for turbine design. In this work, a FVW model is added to the National Renewable Energy Laboratory engineering tool OpenFAST to allow for the aerodynamic modeling of highly flexible turbines along with the aero-hydro-servo-elastic response capabilities of OpenFAST. Code implementation is not complete and thus no results are presented. However, areas of potential modeling improvement are discussed, such as modifications to the circulation calculation and the need for different FVW input model specifications when high-thrust cases are being modeled.

17 WIND ENERGY↗