Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toolkit”

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 973 records · Page 54

The abelian confinement mechanism revisited: New aspects of the Georgi–Glashow model

The confinement problem remains one of the most difficult problems in theoretical physics. An important step toward the solution of this problem is Polyakov’s work on abelian confinement. The Georgi–Glashow model is a natural testing ground for this mechanism which has been surprising us by its richness and wide applicability. In this work, we shed light on two new aspects of this model in 2+1 D. First, we develop a many-body description of the effective degrees of freedom. Namely, we consider a non-relativistic gas of W-bosons in the background of monopole–instanton plasma. Many-body treatment is a standard toolkit in condensed matter physics. However, we add a new twist by supplying the monopole–instantons as external background field. Using this construction along with a mean-field approximation, we calculate the form of the potential between two electric probes as a function of their separation. This potential is expressed in terms of the Meijer-G function which interpolates between logarithmic and linear behavior at small and large distances, respectively. Second, we develop a systematic approach to integrate out the effect of the W-bosons at finite temperature in the range 0≤T<M{sub W}, where M{sub W} is the W-boson mass, starting from the full relativistic partition function of the Georgi–Glashow model. Using a heat kernel expansion that takes into account the non-trivial thermal holonomy, we show that the partition function describes a three-dimensional two-component Coulomb gas. We repeat our analysis using the many-body description which yields the same result and provides a check on our formalism. At temperatures close to the deconfinement temperature, the gas becomes essentially two-dimensional recovering the partition function of the dual sine-Gordon model that was considered in a previous work. -- Highlights: • We consider the problem of abelian confinement in the Georgi–Glashow model from a new perspective. • We develop a many-body description of the degrees of freedom of this model in the background of monopole–instanton plasma. • We use the many-body formalism to find an analytic expression of the potential between two probes at zero temperature. • We also use a systematic approach to integrate out the W-bosons starting from the full relativistic partition function. • This results in a three-dimensional two-component Coulomb gas with long range and Aharonov–Bohm phase interaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Hamiltonian approach to Thermodynamics

In the present work we develop a strictly Hamiltonian approach to Thermodynamics. A thermodynamic description based on symplectic geometry is introduced, where all thermodynamic processes can be described within the framework of Analytic Mechanics. Our proposal is constructed on top of a usual symplectic manifold, where phase space is even dimensional and one has well-defined Poisson brackets. The main idea is the introduction of an extended phase space where thermodynamic equations of state are realized as constraints. We are then able to apply the canonical transformation toolkit to thermodynamic problems. Throughout this development, Dirac’s theory of constrained systems is extensively used. To illustrate the formalism, we consider paradigmatic examples, namely, the ideal, van der Waals and Clausius gases. - Highlights: • A strictly Hamiltonian approach to Thermodynamics is proposed. • Dirac’s theory of constrained systems is extensively used. • Thermodynamic equations of state are realized as constraints. • Thermodynamic potentials are related by canonical transformations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electric load shape benchmarking for small- and medium-sized commercial buildings

Small- and medium-sized commercial buildings owners and utility managers often look for opportunities for energy cost savings through energy efficiency and energy waste minimization. However, they currently lack easy access to low-cost tools that help interpret the massive amount of data needed to improve understanding of their energy use behaviors. Benchmarking is one of the techniques used in energy audits to identify which buildings are priorities for an energy analysis. Traditional energy performance indicators, such as the energy use intensity (annual energy per unit of floor area), consider only the total annual energy consumption, lacking consideration of the fluctuation of energy use behavior over time, which reveals the time of use information and represents distinct energy use behaviors during different time spans. To fill the gap, this study developed a general statistical method using 24-hour electric load shape benchmarking to compare a building or business/tenant space against peers. Specifically, the study developed new forms of benchmarking metrics and data analysis methods to infer the energy performance of a building based on its load shape. We first performed a data experiment with collected smart meter data using over 2,000 small- and medium-sized businesses in California. We then conducted a cluster analysis of the source data, and determined and interpreted the load shape features and parameters with peer group analysis. Finally, we implemented the load shape benchmarking feature in an open-access web-based toolkit (the Commercial Building Energy Saver) to provide straightforward and practical recommendations to users. The analysis techniques were generic and flexible for future datasets of other building types and in other utility territories.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A novel stochastic modeling method to simulate cooling loads in residential districts

District cooling systems are widely used in urban residential communities in China. Most of such systems are oversized, which leads to wasted investment, low operational efficiency and, thus, waste of energy. The accurate prediction of district cooling loads that can support the rightsizing of cooling plant equipment remains a challenge. This study develops a novel stochastic modeling method that consists of (1) six prototype house models representing most apartments in a district, (2) occupant behavior models of residential buildings reflecting their spatial and temporal diversity as well as their complexity based on a large-scale residential survey in China, and (3) a stochastic sampling process to represent all apartments and occupants in the district. The stochastic method was applied to a case study using the Designer's Simulation Toolkit (DeST) to simulate the cooling loads of a residential district in Wuhan, China. The simulation results agreed well with the measured data based on five performance metrics representing the aggregated cooling consumption, the peak cooling loads, the spatial load distribution, the temporal load distribution and the load profiles. Two prevalent simulation methods were also employed to simulate the district cooling loads. Here, the results showed that oversimplified assumptions about occupant behavior could lead to significant overestimation of the peak cooling load and the total cooling loads in the district. Future work will aim to simplify the workflow and data requirements of the stochastic method for its application, and to explore its use in predicting district heating loads and in commercial or mixed-use districts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Extending Discrete Exterior Calculus to a Fractional Derivative

Fractional partial differential equations (FDEs) are used to describe phenomena that involve a “non-local” or “long-range” interaction of some kind. Accurate and practical numerical approximation of their solutions is challenging due to the dense matrices arising from standard discretization procedures. In this paper, we begin to extend the well-established computational toolkit of Discrete Exterior Calculus (DEC) to the fractional setting, focusing on proper discretization of the fractional derivative. We define a Caputo-like fractional discrete derivative, in terms of the standard discrete exterior derivative operator from DEC, weighted by a measure of distance between $p$-simplices in a simplicial complex. We discuss key theoretical properties of the fractional discrete derivative and compare it to the continuous fractional derivative via a series of numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Solution-based thermodynamic modeling of the Ni-Al-Mo system using first-principles calculations

A solution-based thermodynamic description of the ternary Ni–Al–Mo system is developed here, incorporating first-principles calculations and reported modeling of the binary Ni–Al, Ni–Mo and Al–Mo systems. To search for the configurations with the lowest energies of the N phase, the Alloy Theoretic Automated Toolkit (ATAT) was employed and combined with VASP. The liquid, bcc and γ-fcc phases are modeled as random atomic solutions, and the γ'-Ni3Al phase is modeled by describing the ordering within the fcc structure using two sublattices, summarized as (Al,Mo,Ni)0.75(Al,Mo,Ni)0.25. Thus, γ-fcc and γ'-Ni3Al are modeled with a single Gibbs free energy function with appropriate treatment of the chemical ordering contribution. In addition, notable improvements are the following: first, the ternary effects of Mo and Al in the B2-NiAl and D0a-Ni3Mo phases, respectively, are considered; second, the N-NiAl8Mo3 phase is described as a solid solution using a three-sublattice model; third, the X-Ni14Al75Mo11 phase is treated as a stoichiometric compound. Model parameters are evaluated using first-principles calculations of zero-Kelvin formation enthalpies and reported experimental data. In comparison with the enthalpies of formation for the compounds ψ-AlMo, θ-Al8Mo3 and B2-NiAl, the first-principles results indicate that the N-NiAl8Mo3 phase, which is stable at high temperatures, decomposes into other phases at low temperature. Resulting phase equilibria are summarized in the form of isothermal sections and liquidus projections. To clearly identify the relationship between the γ-fcc and γ'-Ni3Al phases in the ternary Ni–Al–Mo system, the specific γ-fcc and γ'-Ni3Al phase fields are plotted in x(Al)–x(Mo)–T space for a temperature range 1200–1800 K.

36 MATERIALS SCIENCE↗

Studying Two-Dimensional Zeolites with the Tools of Surface Science: MFI Nanosheets on Au(111)

While surface science has provided fundamental insights into a variety a materials, the most used catalysts in the industry, namely zeolites, still remain a challenge. The recent preparation of two-dimensional versions of MFI zeolite frameworks and the possibility of their deposition on electrically conductive supports provides for the first time a viable strategy to perform detailed studies on industrially relevant zeolites using the vast toolkit of surface science. In this work we demonstrate the use of infrared reflection absorption spectroscopy (IRRAS) and synchrotron-based x-ray photoelectron spectroscopy (XPS) to study these materials. Furthermore, polarization modulation IRRAS is used to study the adsorption of methanol and its effect in phonon vibrations of the zeolite framework. The possibility of using surface science methods, in particular under ambient pressure conditions, for the study of well-defined zeolites and other microporous structures opens new avenues to understand structural and mechanistic aspects of these materials as catalysts, adsorbents and molecular sieves.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum

Budding yeasts (subphylum Saccharomycotina) are found in every biome and are as genetically diverse as plants or animals. To understand budding yeast evolution, we analyzed the genomes of 332 yeast species, including 220 newly sequenced ones, which represent nearly one-third of all known budding yeast diversity. In this study, we establish a robust genus-level phylogeny comprising 12 major clades, infer the timescale of diversification from the Devonian period to the present, quantify horizontal gene transfer (HGT), and reconstruct the evolution of 45 metabolic traits and the metabolic toolkit of the budding yeast common ancestor (BYCA). We infer that BYCA was metabolically complex and chronicle the tempo and mode of genomic and phenotypic evolution across the subphylum, which is characterized by very low HGT levels and widespread losses of traits and the genes that control them. More generally, our results argue that reductive evolution is a major mode of evolutionary diversification.

59 BASIC BIOLOGICAL SCIENCES↗

CRISPR-Cas9 Circular Permutants as Programmable Scaffolds for Genome Modification

The ability to engineer natural proteins is pivotal to a future, pragmatic biology. CRISPR proteins have revolutionized genome modification, yet the CRISPR-Cas9 scaffold is not ideal for fusions or activation by cellular triggers. Here in this paper, we show that a topological rearrangement of Cas9 using circular permutation provides an advanced platform for RNA-guided genome modification and protection. Through systematic interrogation, we find that protein termini can be positioned adjacent to bound DNA, offering a straightforward mechanism for strategically fusing functional domains. Additionally, circular permutation enabled protease-sensing Cas9s (ProCas9s), a unique class of single-molecule effectors possessing programmable inputs and outputs. ProCas9s can sense a wide range of proteases, and we demonstrate that ProCas9 can orchestrate a cellular response to pathogen-associated protease activity. Together, these results provide a toolkit of safer and more efficient genome-modifying enzymes and molecular recorders for the advancement of precision genome engineering in research, agriculture, and biomedicine.

59 BASIC BIOLOGICAL SCIENCES↗

helixCAM: A platform for programmable cellular assembly in bacteria and human cells

Interactions between cells are indispensable for signaling and creating structure. The ability to direct precise cell-cell interactions would be powerful for engineering tissues, understanding signaling pathways, and directing immune cell targeting. In humans, intercellular interactions are mediated by cell adhesion molecules (CAMs). However, endogenous CAMs are natively expressed by many cells and tend to have cross-reactivity, making them unsuitable for programming specific interactions. Here, we showcase “helixCAM,” a platform for engineering synthetic CAMs by presenting coiled-coil peptides on the cell surface. helixCAMs were able to create specific cell-cell interactions and direct patterned aggregate formation in bacteria and human cells. Based on coiled-coil interaction principles, we built a set of rationally designed helixCAM libraries, which led to the discovery of additional high-performance helixCAM pairs. Finally, we applied this helixCAM toolkit for various multicellular engineering applications, such as spherical layering, adherent cell targeting, and surface patterning.

59 BASIC BIOLOGICAL SCIENCES↗

Alkali-activated materials

This paper, which forms part of the UNEP White Papers series on Eco-Efficient Cements, provides a brief discussion of the class of cementing materials known as ‘alkali-activated binders’, which are identified to have potential for utilization as a key component of a sustainable future global construction materials industry. These cements are not expected to offer a like-for-like replacement of Portland cement across its full range of applications, for reasons related to supply chain limitations, practical challenges in some modes of application, and the need for careful control of formulation and curing. However, when produced using locally-available raw materials, with well-formulated mix designs (including in particular consideration of the environmental footprint of the alkaline activator) and production under adequate levels of quality control, alkali-activated binders are potentially an important and cost-effective component of the future toolkit of sustainable construction materials.

36 MATERIALS SCIENCE↗

Linear model decision trees as surrogates in optimization of engineering applications

Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. Importantly, we observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).

42 ENGINEERING↗

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗

Scale and the representation of human agency in the modeling of agroecosystems

Human agency is an essential determinant of the dynamics of agroecosystems. However, the manner in which agency is represented within different approaches to agroecosystem modeling is largely contingent on the scales of analysis and the conceptualization of the system of interest. While appropriate at times, narrow conceptualizations of agroecosystems can preclude consideration for how agency manifests at different scales, thereby marginalizing processes, feedbacks, and constraints that would otherwise affect model results. Modifications to the existing modeling toolkit may therefore enable more holistic representations of human agency. Model integration can assist with the development of multi-scale agroecosystem modeling frameworks that capture different aspects of agency. In addition, expanding the use of socioeconomic scenarios and stakeholder participation can assist in explicitly defining context-dependent elements of scale and agency. Finally, such approaches, however, should be accompanied by greater recognition of the meta agency of model users and the need for more critical evaluation of model selection and application.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Scalable Metropolis Monte Carlo for simulation of hard shapes

We design and implement a scalable hard particle Monte Carlo simulation toolkit (HPMC), and release it open source as part of HOOMD-blue. HPMC runs in parallel on many CPUs and many GPUs using domain decomposition. We employ BVH trees instead of cell lists on the CPU for fast performance, especially with large particle size disparity, and optimize inner loops with SIMD vector intrinsics on the CPU. Our GPU kernel proposes many trial moves in parallel on a checkerboard and uses a block-level queue to redistribute work among threads and avoid divergence. HPMC supports a wide variety of shape classes, including spheres/disks, unions of spheres, convex polygons, convex spheropolygons, concave polygons, ellipsoids/ellipses, convex polyhedra, convex spheropolyhedra, spheres cut by planes, and concave polyhedra. NVT and NPT ensembles can be run in 2D or 3D triclinic boxes. Additional integration schemes permit Frenkel–Ladd free energy computations and implicit depletant simulations. In a benchmark system of a fluid of 4096 pentagons, HPMC performs 10 million sweeps in 10 min on 96 CPU cores on XSEDE Comet. The same simulation would take 7.6 h in serial. HPMC also scales to large system sizes, and the same benchmark with 16.8 million particles runs in 1.4 h on 2048 GPUs on OLCF Titan.

Computer Science↗

eqtools: Modular, extensible, open-source, cross-machine Python tools for working with magnetic equilibria

As plasma physics research for fusion energy transitions to an increasing emphasis on cross-machine collaboration and numerical simulation, it becomes increasingly important that portable tools be developed to enable data from diverse sources to be analyzed in a consistent manner. This paper presents eqtools, a modular, extensible, open-source toolkit implemented in the Python programming language for handling magnetic equilibria and associated data from tokamaks. eqtools provides a single interface for working with magnetic equilibrium data, both for handling derived quantities and mapping between coordinate systems, extensible to function with data from different experiments, data formats, and magnetic reconstruction codes, replacing the diverse, non-portable solutions currently in use. Moreover, while the open-source Python programming language offers a number of advantages as a scripting language for research purposes, the lack of basic tokamak-specific functionality has impeded the adoption of the language for regular use. Implementing equilibrium-mapping tools in Python removes a substantial barrier to new development in and porting legacy code into Python. In this paper, we introduce the design of the eqtools package and detail the workflow for usage and expansion to additional devices. Additionally, the implementation of a novel three-dimensional spline solution (in two spatial dimensions and in time) is also detailed. Finally, verification and benchmarking for accuracy and speed against existing tools are detailed. Wider deployment of these tools will enable efficient sharing of data and software between institutions and machines as well as self-consistent analysis of the shared data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solution of the Skyrme-Hartree–Fock–Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis. (VIII) HFODD (v2.73y): A new version of the program

Here, we describe the new version (v2.73y) of the code hfodd which solves the nuclear Skyrme Hartree–Fock or Skyrme Hartree–Fock–Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we have implemented the following new features: (i) full proton–neutron mixing in the particle–hole channel for Skyrme functionals, (ii) the Gogny force in both particle–hole and particle–particle channels, (iii) linear multi-constraint method at finite temperature, (iv) fission toolkit including the constraint on the number of particles in the neck between two fragments, calculation of the interaction energy between fragments, and calculation of the nuclear and Coulomb energy of each fragment, (v) the new version 200d of the code hfbtho, together with an enhanced interface between HFBTHO and HFODD, (vi) parallel capabilities, significantly extended by adding several restart options for large-scale jobs, (vii) the Lipkin translational energy correction method with pairing, (viii) higher-order Lipkin particle-number corrections, (ix) interface to a program plotting single-particle energies or Routhians, (x) strong-force isospin-symmetry-breaking terms, and (xi) the Augmented Lagrangian Method for calculations with 3D constraints on angular momentum and isospin. Finally, an important bug related to the calculation of the entropy at finite temperature and several other little significant errors of the previous published version were corrected.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PANNA: Properties from Artificial Neural Network Architectures

We report prediction of material properties from first principles is often a computationally expensive task. Recently, artificial neural networks and other machine learning approaches have been successfully employed to obtain accurate models at a low computational cost by leveraging existing example data. Here, we present a software package “Properties from Artificial Neural Network Architectures” (PANNA) that provides a comprehensive toolkit for creating neural network models for atomistic systems following the Behler–Parrinello topology. Besides the core routines for neural network training, it includes data parser, descriptor builder for Behler–Parrinello class of symmetry functions and force-field generator suitable for integration within molecular dynamics packages. PANNA offers a variety of activation and cost functions, regularization methods, as well as the possibility of using fully-connected networks with custom size for each atomic species. PANNA benefits from the optimization and hardware-flexibility of the underlying TensorFlow engine which allows it to be used on multiple CPU/GPU/TPU systems, making it possible to develop and optimize neural network models based on large datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗