Search NASA⌕ Search

SEARCH · Search NASA

Results for “complexity analysis”

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 883 records · Page 49

GrAMME: Semisupervised Learning Using Multilayered Graph Attention Models

Modern data analysis pipelines are becoming increasingly complex due to the presence of multiview information sources. While graphs are effective in modeling complex relationships, in many scenarios, a single graph is rarely sufficient to succinctly represent all interactions, and hence, multilayered graphs have become popular. Though this leads to richer representations, extending solutions from the single-graph case is not straightforward. Consequently, there is a strong need for novel solutions to solve classical problems, such as node classification, in the multilayered case. In this article, we consider the problem of semisupervised learning with multilayered graphs. Though deep network embeddings, e.g., DeepWalk, are widely adopted for community discovery, we argue that feature learning with random node attributes, using graph neural networks, can be more effective. To this end, we propose to use attention models for effective feature learning and develop two novel architectures, GrAMME-SG and GrAMME-Fusion, that exploit the interlayer dependences for building multilayered graph embeddings. Using empirical studies on several benchmark data sets, we evaluate the proposed approaches and demonstrate significant performance improvements in comparison with the state-of-the-art network embedding strategies. The results also show that using simple random features is an effective choice, even in cases where explicit node attributes are not available.

97 MATHEMATICS AND COMPUTING↗

URBANopt: An Open-Source Software Development Kit for Community and Urban District Energy Modeling: Preprint

Urban building modeling tools are developing rapidly; these tools use emerging simulation workflows for specific urban environmental design tasks, such as assessing the impacts of energy efficiency technologies at a district scale. However, with the emergence of new environmental design tasks, addressing all possible use cases and tasks is challenging and cannot be covered by a single tool. Urban-scale analysis at this level of complexity often requires linking multiple emerging tools, rather than using a single tool, to adequately evaluate a variety of possible fields in urban environmental design. To achieve this, flexible platforms are needed to support multiple input formats (e.g., geometric and non-geometric building properties), enabling the mapping of such inputs to underlying simulation engines. This paper provides an overview of the open-source URBANopt Software Development Kit (SDK) for modeling high-performance buildings and energy systems at a district scale. URBANopt's flexible SDK is composed of several modules that can be customized to integrate with other tools and generate new workflows to perform urban environmental design tasks, such as capturing interactions between individual buildings, district energy systems, distributed energy resources, and the electric distribution grid. We describe the functionality of the core SDK modules in URBANopt (called Core Gem, GeoJSON Gem, and Scenario Gem) and discuss the flexibility of these modules as a means of integration with a variety of tools. We also document and demonstrate technical details of writing and combining new modules to create customized workflows. Finally, we present a case study that uses the URBANopt SDK to model a hypothetical mixed-use urban project and simulate various scenarios to meet district energy performance goals.

buildings↗

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Lastly, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical analysis of HAADF-STEM images to determine the surface coverage and distribution of immobilized molecular complexes

The surface immobilization of molecular catalysts is attractive because it combines the benefits of homogeneous and heterogeneous catalysis. However, determining the surface coverage and distribution of a molecular catalyst on a solid support is often challenging, inhibiting our ability to design improved catalytic systems. Here, in this work, we demonstrate that the combination of scanning transmission electron microscopy (STEM) and image analysis of the individual positions of heavy atoms in transition metal complexes via a convolutional neural network (CNN) allows statistically robust determination of the surface coverage and distribution of immobilized molecular catalysts. These observations provide information about how changes in the functionalization conditions, attachment group, and structure of the molecular catalyst affect the surface coverage and distribution, providing insight into the chemical mechanism of surface immobilization. The method could be generally valuable for correlating the surface coverage and distribution to the activity, selectivity, and stability of a catalytic system.

HAADF-STEM↗

Evolution of the nanoporous microstructure of sintered Ag at high temperature using in-situ X-ray nanotomography

Silver pastes sintering is a potential candidate for die bonding in power electronic modules. The joints, obtained by sintering, exhibit a significant pore fraction, thus reducing the density of the material compared to bulk silver. The mechanical properties (Young's modulus, yield strength and ultimate tensile stress) are then drastically altered. However, while careful analysis of the nanoporous structure has been reported in 2D, little is known about its quantitative spatial evolution during thermal aging and more specifically during temperature jumps endured by the assembly during operation. High temperature evolutions of the 3D nanoporous structure have been observed in-situ using a heater fitted into the beamline 6-2c at SLAC-SSRL. Segmentation of the porosity and subsequent statistical analysis of the tomographic dataset reveal a complex evolution of the porous nanostructure including growth, separation and coalescence at a constant density. Such an analysis provides insight into the microstructural evolution of sintered nanoporous Ag joints in-service.

36 MATERIALS SCIENCE↗

Fused Global-Local Economic Model Predictive Control for Real-Time Eco-Optimal Control of a Heavy Duty Truck

This paper introduces a fused global-local economic model predictive control (MPC) with varying time steps for maximizing the fuel economy of a heavy-duty truck, where global MPC is performed over a long horizon using cloud-based dynamic programming (DP), and local MPC is performed on-board the vehicle using sequential quadratic programming (SQP). Because of the computational complexity of DP, an in-depth analysis of how the discretization levels, time constant, and horizon length influence the energy savings and computational complexity of online DP has been performed. The fused optimization technique has been evaluated over different route profiles, accounting for traffic, using a medium-fidelity proprietary Volvo model. Results show 4-7% fuel savings, as compared to a baseline model.

33 ADVANCED PROPULSION SYSTEMS↗

Taxogenomic analysis of Pichia senei sp. nov. and new insights into hybridization events in the Pichia cactophila species complex

Three strains of a novel yeast species were isolated from necrotic cactus tissues of Cereus saddianus and Micranthocereus dolichospermaticus and from phytotelmata of Bromelia karatas. DNA sequence analysis of the Internal Transcribed Spacer (ITS) region and D1/D2 domains of the large subunit ribosomal RNA, along with whole genome phylogenomic analysis, showed that this yeast is most closely related to Pichia insulana, Pichia cactophila, and Pichia inconspicua. The new species differs by 10–13 nucleotide substitutions from these species in D1/D2 sequences and exhibits <90% genome-wide average nucleotide identity to them. The name Pichia senei sp. nov. is proposed for the novel species, which is homothallic and produces asci with one to four hat-shaped ascospores. The holotype is CBS 16311 (MycoBank MB 858723). Taxogenomic analyses of the P. cactophila species complex, including P. senei, provide new insights about the hybridizations events that shaped this group. Pichia insulana and P. inconspicua are identified as the parental lineages that originated P. cactophila, and P. senei also appears closely related to one of the progenitors of P. inconspicua. We assess phylogeny, heterozygosity, and ploidy to explore the processes shaping diversity, showing how genomic data support yeast species delimitation and reveal complex hybridization.

59 BASIC BIOLOGICAL SCIENCES↗

Dose Consequence and Probabilistic Risk Assessment Integration into Digital Documented Safety Analysis

This is an intern poster presentation. The current reactor authorization process is complex and error prone. The development of a digital Documented Safety Analysis has been proposed to provide an automated and integrated solution to enhance the design and authorization process. The digital DSA will consist of interlinked models, analyses, and reports, all of which will be updated using automated workflows when design changes are made. This poster examines the integration of the probabilistic risk assessment (PRA) with the transient and dose consequence analyses. A PRA for a generic high temperature gas-cooled reactor (HTGR) has been constructed which will drive the input parameters for a transient analysis model currently being constructed. Dose consequence will be determined using the results of the transient analysis and the Radiation Safety Analysis Computer (RSAC) code. Dose consequence data will then be input back into the PRA to drive design parameters.

42 ENGINEERING↗

Trace Analysis and Reaction Monitoring by Nanophotonic Ionization Mass Spectrometry from Elevated Bowtie and Silicon Nanopost Arrays

We present that silicon nanopost arrays (NAPA) are used in trace analysis by mass spectrometry (MS) because they enable highly efficient ion production from small molecules and thin tissue sections by UV laser desorption ionization (LDI). Such nanophotonic ionization of adsorbates relies on localized interactions between a nanostructured substrate and laser radiation. In LDI from NAPA, only the component of the oscillating electric field vector that is parallel with the posts couples the laser energy into the nanostructure. Enhancements in control over adsorbate ionization and fragmentation are expected if the surface-parallel component can also interact with the nanostructure. Here, an alternative nanophotonic ionization platform is introduced for LDI-MS, the elevated bowtie (EBT) array by adding triangular chromium features on top of silicon post pairs to form bowties. Compared to NAPA, the threshold fluence for ionization from EBT is lower, and at low laser fluences the ionization efficiency is increased by a factor of ≈17. The EBT platform with optimized apex angle exhibits a higher survival yield for molecular ions produced from biomolecules and xenobiotics and allows more control over fragmentation by adjusting the fluence. Finally, these unique nanophotonic ionization attributes are utilized for trace analysis and reaction monitoring in complex biological samples.

36 MATERIALS SCIENCE↗

Structural and functional validation of a highly specific Smurf2 inhibitor

Smurf1 and Smurf2 are two closely related member of the HECT (homologous to E6AP carboxy terminus) E3 ubiquitin ligase family and play important roles in the regulation of various cellular processes. Both were initially identified to regulate transforming growth factor‐β and bone morphogenetic protein signaling pathways through regulating Smad protein stability and are now implicated in various pathological processes. Generally, E3 ligases, of which over 800 exist in humans, are ideal targets for inhibition as they determine substrate specificity; however, there are few inhibitors with the ability to precisely target a particular E3 ligase of interest. In this work, we explored a panel of ubiquitin variants (UbVs) that were previously identified to bind Smurf1 or Smurf2. In vitro binding and ubiquitination assays identified a highly specific Smurf2 inhibitor, UbV S2.4, which was able to inhibit ligase activity with high potency in the low nanomolar range. Orthologous cellular assays further demonstrated high specificity of UbV S2.4 toward Smurf2 and no cross‐reactivity toward Smurf1. Structural analysis of UbV S2.4 in complex with Smurf2 revealed its mechanism of inhibition was through targeting the E2 binding site. In summary, we investigated several protein‐based inhibitors of Smurf1 and Smurf2 and identified a highly specific Smurf2 inhibitor that disrupts the E2–E3 protein interaction interface.

59 BASIC BIOLOGICAL SCIENCES↗

Chaos, complexity, and random matrices

Chaos and complexity entail an entropic and computational obstruction to describing a system, and thus are intrinsically difficult to characterize. In this paper, we consider time evolution by Gaussian Unitary Ensemble (GUE) Hamiltonians and analytically compute out-of-time-ordered correlation functions (OTOCs) and frame potentials to quantify scrambling, Haar-randomness, and circuit complexity. While our random matrix analysis gives a qualitatively correct prediction of the late-time behavior of chaotic systems, we find unphysical behavior at early times including an O(1) scrambling time and the apparent breakdown of spatial and temporal locality. The salient feature of GUE Hamiltonians which gives us computational traction is the Haar-invariance of the ensemble, meaning that the ensemble-averaged dynamics look the same in any basis. Motivated by this property of the GUE, we introduce k-invariance as a precise definition of what it means for the dynamics of a quantum system to be described by random matrix theory. Here,we envision that the dynamical onset of approximate k-invariance will be a useful tool for capturing the transition from early-time chaos, as seen by OTOCs, to late-time chaos, as seen by random matrix theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗

High Strain Response and Dielectric Properties of Bi 1/2 (Na 0.78 K 0.22 ) 1/2 TiO 3 Ceramics Doped with (Fe 1/2 Nb 1/2 ) 4+

Environmentally friendly piezoceramics Bi 1/2 (Na 0.78 K 0.22 ) 1/2 Ti 1-x (Fe 1/2 Nb 1/2 ) x O 3 (BNKT-FN) with FeNb (FN) ratios 0, 0.01, 0.02, 0.03 and 0.04 were synthesized via conventional solid-state reaction (CSSR). X-ray diffraction (XRD) analysis revealed that the FN complexes were completely dissolved in BNKT perovskite and the structural transition occurred from ferroelectric to relaxor state. The average grain sizes for x = 0, 0.01, 0.02, 0.03 and 0.04 were found to be 0.34 μm, 0.35 μm, 0.38 μm, 0.32 μm and 0.39 μm, respectively. It is observed that the grain size increased with the addition of FN, except for x = 0.03. Here, the coercive field (E c ) and remnant polarization (P r ) of the system were decreased with increasing FN concentration and the highest E c of 21.1 kV/cm and P r of 29.2 μC/cm 2 were noted for the base pure sample. In addition, the increase of FN content steadily broadened the dielectric anomaly at maximum dielectric constant temperature (T m ) in the BNKT-FN system. Consequently, unipolar normalized strain (S max /E max ) was found to be enhanced monotonically achieving the highest value of around 654 pm/V (x = 0.02). These properties suggest that the investigated system could be a promising eco-friendly candidate for piezoelectric actuators and many other applications requiring high piezoelectric induced strain.

25 ENERGY STORAGE↗

A new predictive multi-zone model for HCCI engine combustion

Here, this work introduces a new predictive multi-zone model for the description of combustion in Homogeneous Charge Compression Ignition (HCCI) engines. The model exploits the existing OpenSMOKE++ computational suite to handle detailed kinetic mechanisms, providing reliable predictions of the in-cylinder auto-ignition processes. All the elements with a significant impact on the combustion performances and emissions, like turbulence, heat and mass exchanges, crevices, residual burned gases, thermal and feed stratification are taken into account. Compared to other computational approaches, this model improves the description of mixture stratification phenomena by coupling a wall heat transfer model derived from CFD application with a proper turbulence model. Furthermore, the calibration of this multi-zone model requires only three parameters, which can be derived from a non-reactive CFD simulation: these adaptive variables depend only on the engine geometry and remain fixed across a wide range of operating conditions, allowing the prediction of auto-ignition, pressure traces and pollutants. This computational framework enables the use of detail kinetic mechanisms, as well as Rate of Production Analysis (RoPA) and Sensitivity Analysis (SA) to investigate the complex chemistry involved in the auto-ignition and the pollutants formation processes. In the final sections of the paper, these capabilities are demonstrated through the comparison with experimental data.

30 DIRECT ENERGY CONVERSION↗

Experimental and fuel-surrogates modeling study of the high-pressure pyrolysis of specialty cetane number fuels: implications for fall-off in ethylene unimolecular dissociation

Single pulse shock tube experiments were conducted at 50 atm nominal pressure and 4 ms nominal reaction time over a temperature range of 900–1800 K, to study the pyrolysis speciation of a multi-component jet fuel, F-24, and six cetane number (CN) specialty fuels - CN30, CN35, CN40, CN45, CN50, and CN55. Gas chromatography (GC) was used to qualitatively and quantitatively analyze the post shock gases. The relationship between the formation of key pyrolysis species and the chemically controlled combustion propensity as reflected by the cetane number of each fuel was examined. A surrogate-based mechanism from the CRECK Modelling Group and chemical-functional group based optimized surrogates (CFGO) were used to simulate the pyrolysis speciation results. The model was able to capture the chemistry of most species except two important pyrolysis intermediates – ethylene and acetylene. Chemical kinetic analyses were performed to identify the important reactions which affect the chemistry of these species; however, the rate parameters of critical reactions were found to be unsuitable for simulating the present high-pressure studies. Here, to address this unsuitability, a theory-based fall-off analysis for three reactions representing the decomposition of ethylene and subsequent formation of acetylene was performed, and these are included in an updated version of the CRECK mechanism. This update resolves discrepancies between the experimental results and simulations for ethylene and acetylene. Reaction flux analyses using the updated surrogate model were also performed to identify the important reaction pathways responsible for the formation of crucial species and to provide an analysis of the chemistry of complex multi-component fuel systems. The fundamental reactions responsible for driving pyrolysis chemistry were greatly influenced by the chemical functional groups present in these fuels. In addition to updating the rate parameters of specific reactions to improve modeling, this study also emphasizes the effectiveness of the fuel-surrogate approach, where surrogates representing the chemical functional group composition of the parent fuel serve as a valuable tool for predicting the combustion chemistry of novel fuels.

Chemical Kinetics↗

The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient

Numerical schemes for nonlinear parabolic equations based on the harmonic averaging of cell-centered diffusion coefficients break down when some of these coefficients go to zero or their ratio grows. To tackle this problem, we propose new mimetic finite difference schemes that use a staggered discretization of the diffusion coefficient. The primary mimetic operator approximates div (k•); the derived (dual) mimetic operator approximates - ∇(•). The new mimetic schemes preserve symmetry and positive-definiteness of the continuum problem which allows us to use algebraic solvers with optimal complexity. We perform detailed numerical analysis of the new schemes for linear elliptic problems and a specially designed linear parabolic problem that has solution dynamics typical for nonlinear problems. We also show that the new schemes are competitive with the state-of-the-art schemes for steady-state problems but provide much more accurate solution dynamics for the transient problem.

97 MATHEMATICS AND COMPUTING↗

G EANT 4 atomic relaxation data for transfermium nuclei (Z = 101–104)

Advanced theoretical methods can accurately calculate various atomic observables and predict electronic structure. Still, systematic computations of the radiative and non-radiative transition probabilities and energies are missing for the actinides and all the transfermium elements. However, these compilations are needed for comprehensive Monte-Carlo simulations (such as GEANT4) of the radioactive decay of transfermium nuclei. These simulations can forma basis for data analysis of experiments, especially with complex detection setups. Investigation of the transfermium nuclei is crucial for understanding the nature of the nuclear force. In this study, simulations based on data from the Jena Atomic Calculator (JAC) and the data from the Evaluated Atomic Data Library (EADL) present in GEANT4 were found compatible for the three elements Ba(Z = 56), U(Z = 92), and Fm(Z = 100), thus, validating the JAC calculations. For Z> 100, we also found sound agreement between simulations that used data generated with JAC and experimental results involving No(Z = 102) and Rf(Z = 104) isotopes. In conclusion, these results demonstrate that JAC can produce reliable atomic data sets for transfermium elements, which will assist in analyzing nuclear-decay-spectroscopy experiments.

GEANT4↗