Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computable general equilibrium”

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 109 records · Page 6

Watch and learn—a generalized approach for transferrable learning in deep neural networks via physical principles

Transfer learning refers to the use of knowledge gained while solving a machine learning task and applying it to the solution of a closely related problem. Such an approach has enabled scientific breakthroughs in computer vision and natural language processing where the weights learned in state-of-the-art models can be used to initialize models for other tasks which dramatically improve their performance and save computational time. Here we demonstrate an unsupervised learning approach augmented with basic physical principles that achieves fully transferrable learning for problems in statistical physics across different physical regimes. By coupling a sequence model based on a recurrent neural network to an extensive deep neural network, we are able to learn the equilibrium probability distributions and inter-particle interaction models of classical statistical mechanical systems. Our approach, distribution-consistent learning, DCL, is a general strategy that works for a variety of canonical statistical mechanical models (Ising and Potts) as well as disordered interaction potentials. Using data collected from a single set of observation conditions, DCL successfully extrapolates across all temperatures, thermodynamic phases, and can be applied to different length-scales. This constitutes a fully transferrable physics-based learning in a generalizable approach.

97 MATHEMATICS AND COMPUTING↗

Theoretical studies of X-ray emission from supernova remnants and other non-equilibrium astrophysical plasmas

To improve the accuracy of thermal X-ray emission spectrum predictions, new computations of electron excitation rates, dielectric recombination rates, and updated codes were applied to existing data on supernova remnants and O star winds. Because of its general predominance along with low ionization and recombination rates, Fe(XVII) was analyzed with the inclusion of resonances to the excitation rate data. Corrections to the existing dielectric recombination rates based on detailed Fe (XVII) calculations and the more approximate Quantum Defect Method were accomplished. by using Distorted Wave ionization rates and reliable excitation rates, models for X-ray, ultraviolet, and optical emission line intensities of non-radiative shock wave can be constructed that supply data on shock velocity, preshock ionization states, and elemental abundances. The radiative driving force that drives the shock wave is calculated through line opacities of the shocked gas. X-ray emission spectra are also calculated for O star winds.

Dupree, A. K.↗

Computing Properties Of Chemical Mixtures At Equilibrium

Scientists and engineers need data on chemical equilibrium compositions to calculate theoretical thermodynamic properties of chemical systems. Information essential in design and analysis of such equipment as compressors, turbines, nozzles, engines, shock tubes, heat exchangers, and chemical-processing equipment. CET93 is general program that calculates chemical equilibrium compositions and properties of mixtures for any chemical system for which thermodynamic data are available. Includes thermodynamic data for more than 1,300 gaseous and condensed species and thermal-transport data for 151 gases. Written in FORTRAN 77.

Mcbride, B. J.↗

Statistical mechanics of light elements at high pressure. VI - Liquid-state calculations with Thomas-Fermi-Dirac theory

A model free energy is developed for hydrogen-helium mixtures based on solid-state Thomas-Fermi-Dirac calculations at pressures relevant to the interiors of giant planets. Using a model potential similar to that for a two-component plasma, effective charges for the nuclei (which are in general smaller than the actual charges because of screening effects) are parameterized, being constrained by calculations at a number of densities, compositions, and lattice structures. These model potentials are then used to compute the equilibrium properties of H-He fluids using a charged hard-sphere model. The results find critical temperatures of about 0 K, 500 K, and 1500 K, for pressures of 10, 100, and 1000 Mbar, respectively. These phase separation temperatures are considerably lower (approximately 6,000-10,000 K) than those found from calculations using free electron perturbation theory, and suggest that H-He solutions should be stable against phase separation in the metallic zones of Jupiter and Saturn.

Macfarlane, J. J.↗

Hierarchic Extensions in the Static and Dynamic Analysis of Elastic Beams

Approximate solutions of static and dynamic beam problems by the p-version of the finite element method are investigated. Within a hierarchy of engineering beam idealizations, rigorous formulations of the strain and kinetic energies for straight and circular beam elements are presented. These formulations include rotating coordinate system effects and geometric nonlinearities to allow for the evaluation of vertical axis wind turbines, the motivating problem for this research. Hierarchic finite element spaces, based on extensions of the polynomial orders used to approximate the displacement variables, are constructed. The developed models are implemented into a general purpose computer program for evaluation. Quality control procedures are examined for a diverse set of sample problems. These procedures include estimating discretization errors in energy norm and natural frequencies, performing static and dynamic equilibrium checks, observing convergence for qualities of interest, and comparison with more exacting theories and experimental data. It is demonstrated that p-extensions produce exponential rates of convergence in the approximation of strain energy and natural frequencies for the class of problems investigated.

Watson, Robert A.↗

On the maximum energy release in flux-rope models of eruptive flares

We determine the photospheric boundary conditions which maximize the magnetic energy released by a loss of ideal-magnetohydrodynamic (MHD) equilibrium in two-dimensional flux-rope models. In these models a loss of equilibrium causes a transition of the flux rope to a lower magnetic energy state at a higher altitude. During the transition a vertical current sheet forms below the flux rope, and reconnection in this current sheet releases additional energy. Here we compute how much energy is released by the loss of equilibrium relative to the total energy release. When the flux-rope radius is small compared to its height, it is possible to obtain general solutions of the Grad-Shafranov equation for a wide range of boundary conditions. Variational principles can then be used to find the particular boundary condition which maximizes the magnetic energy released for a given class of conditions. We apply this procedure to a class of models known as cusp-type catastrophes, and we find that the maximum energy released by the loss of equilibrium is 20.8% of the total energy release for any model in this class. If the additional restriction is imposed that the photospheric magnetic field forms a simple arcade in the absence of coronal currents, then the maximum energy release reduces to 8.6%

Forbes, T. G.↗

Computer program for calculation of complex chemical equilibrium compositions and applications. Part 1: Analysis

This report presents the latest in a number of versions of chemical equilibrium and applications programs developed at the NASA Lewis Research Center over more than 40 years. These programs have changed over the years to include additional features and improved calculation techniques and to take advantage of constantly improving computer capabilities. The minimization-of-free-energy approach to chemical equilibrium calculations has been used in all versions of the program since 1967. The two principal purposes of this report are presented in two parts. The first purpose, which is accomplished here in part 1, is to present in detail a number of topics of general interest in complex equilibrium calculations. These topics include mathematical analyses and techniques for obtaining chemical equilibrium; formulas for obtaining thermodynamic and transport mixture properties and thermodynamic derivatives; criteria for inclusion of condensed phases; calculations at a triple point; inclusion of ionized species; and various applications, such as constant-pressure or constant-volume combustion, rocket performance based on either a finite- or infinite-chamber-area model, shock wave calculations, and Chapman-Jouguet detonations. The second purpose of this report, to facilitate the use of the computer code, is accomplished in part 2, entitled 'Users Manual and Program Description'. Various aspects of the computer code are discussed, and a number of examples are given to illustrate its versatility.

Gordon, Sanford↗

Generalized computer-aided discrete time domain modeling and analysis of dc-dc converters

A generalized discrete time domain modeling and analysis technique is presented for all types of switching regulators using any type of duty-cycle controller, and operating in both continuous and discontinuous inductor current. State space techniques are employed to derive an equivalent nonlinear discrete time model that describes the converter exactly. The system is linearized about its equilibrium state to obtain a linear discrete time model for small signal performance evaluations, such as stability, audiosusceptibility and transient response. The analysis makes extensive use of the digital computer as an analytical tool. It is universal, exact and easy to use.

Lee, F. C.↗

General instability analysis

HOLBOAT computer program provides instability analysis of inhomogeneous, anisotropic, right-circular cylinder or segment under combined loading. Program is based on Kirchhoff-Love hypothesis, general anisotropic constitutive equations, and Flugge's differential equations of equilibrium.

Bianca, C. J.↗

The Riemann disks. I - Equilibrium and secular evolution

The two-dimensional analogs to Riemann S-type ellipsoids are considered as simple fluid models for galactic bars. Their equilibrium is discussed and their evolution under the effects of viscosity, gravitational radiation, and external torques from other components of the galaxy is computed. As a bar loses angular momentum due to external torques, it generally becomes less elongated, and its figure speed may either increase or decrease.

Weinberg, M. D.↗

Neural network based fast prediction of β N limits in HL-2M

Artificial neural networks (NNs) are trained, based on the numerical database, to predict the no-wall and ideal-wall β N limits, due to onset of the n = 1 ( n is the toroidal mode number) ideal external kink instability, for the HL-2M tokamak. The database is constructed by toroidal computations utilizing both the equilibrium code CHEASE (Lütjens et al 1992 Comput. Phys. Commun. 69 287) and the stability code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The stability results show that (1) the plasma elongation generally enhances both β N limits, for either positive or negative triangularity plasmas; (2) the effect is more pronounced for positive triangularity plasmas; (3) the computed no-wall β N limit linearly scales with the plasma internal inductance, with the proportionality coefficient ranging between 1 and 5 for HL-2M; (4) the no-wall limit substantially decreases with increasing pressure peaking factor. Furthermore, both the NN model and the convolutional neural network (CNN) model are trained and tested, producing consistent results. The trained NNs predict both the no-wall and ideal-wall limits with as high as 95% accuracy, compared to those directly computed by the stability code. Additional test cases, produced by the Tokamak Simulation Code (Jardin et al 1993 Nucl. Fusion 33 371), also show reasonable performance of the trained NNs, with the relative error being within 10%. The constructed database provides effective references for the future HL-2M operations. The trained NNs can be used as a real-time monitor for disruption prevention in the HL-2M experiments, or serve as part of the integrated modeling tools for ideal kink stability analysis.

Physics↗

Application of an invariant second-order closure model to compressible turbulent shear layers

A second-order closure model for two-dimensional, compressible shear flows is investigated using the invariant modeling technique developed by Donaldson. The invariant model parameters were originally selected by comparison of model predictions with critical experimental data on basic incompressible flows. Additional modeling for compressible flows has been introduced. A number of different shear flows such as the free shear layer, flat plate boundary layer and a simulation of the mixing region of a chemical laser are computed. The model is consistent with first-order closure turbulence models for equilibrium flows and is further capable of predicting nonequilibrium flows that cannot be correctly solved by eddy viscosity models. The results are in generally good agreement with experimental measurements but suggest the need for inclusion of a Mach number dependent model of the pressure diffusion terms in order to adequately represent high speed compressible flows.

Varma, A. K.↗

Integrated Force Method for Indeterminate Structures

Two methods of solving indeterminate structural-mechanics problems have been developed as products of research on the theory of strain compatibility. In these methods, stresses are considered to be the primary unknowns (in contrast to strains and displacements being considered as the primary unknowns in some prior methods). One of these methods, denoted the integrated force method (IFM), makes it possible to compute stresses, strains, and displacements with high fidelity by use of modest finite-element models that entail relatively small amounts of computation. The other method, denoted the completed Beltrami Mitchell formulation (CBMF), enables direct determination of stresses in an elastic continuum with general boundary conditions, without the need to first calculate displacements as in traditional methods. The equilibrium equation, the compatibility condition, and the material law are the three fundamental concepts of the theory of structures. For almost 150 years, it has been commonly supposed that the theory is complete. However, until now, the understanding of the compatibility condition remained incomplete, and the compatibility condition was confused with the continuity condition. Furthermore, the compatibility condition as applied to structures in its previous incomplete form was inconsistent with the strain formulation in elasticity.

Hopkins, Dale A.↗

On principles of emergent organization

After more than a century of concerted effort, physics still lacks basic principles of spontaneous organization. To appreciate why, we first state the problem, outline historical approaches, and survey the present state of the physics of self-organization. This frames the particular challenges arising from mathematical intractability and the resulting need for computational approaches, as well as those arising from a chronic failure to define structure. Then an overview of two modern mathematical formulations of organization—intrinsic computation and evolution operators—lays out a way to overcome these challenges. Together, the vantage point they afford shows how to account for the emergence of structured states via a statistical mechanics of systems arbitrarily far from equilibrium. The result is a constructive path forward to principles of organization that builds on mathematical identification of structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗

On symmetry-resolved generalized entropies

Symmetry-resolved entanglement, capturing the refined structure of quantum entanglement in systems with global symmetries, has attracted a lot of attention recently. In this manuscript, introducing the notion of symmetry-resolved generalized entropies, we aim to develop a computational framework suitable for the study of excited state symmetry-resolved entanglement as well as the dynamical evolution of symmetry-resolved entanglement in symmetry-preserving out-of-equilibrium settings. We illustrate our framework using the example of (1+1)-d free massless compact boson theory, and benchmark our results using lattice computation in the XX chain. As a byproduct, our computational framework also provides access to the probability distribution of the symmetry charge contained within a subsystem and the corresponding full counting statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A generalized chemistry version of SPARK

An extension of the reacting H2-air computer code SPARK is presented, which enables the code to be used on any reacting flow problem. Routines are developed calculating in a general fashion, the reaction rates, and chemical Jacobians of any reacting system. In addition, an equilibrium routine is added so that the code will have frozen, finite rate, and equilibrium capabilities. The reaction rate for the species is determined from the law of mass action using Arrhenius expressions for the rate constants. The Jacobian routines are determined by numerically or analytically differentiating the law of mass action for each species. The equilibrium routine is based on a Gibbs free energy minimization routine. The routines are written in FORTRAN 77, with special consideration given to vectorization. Run times for the generalized routines are generally 20 percent slower than reaction specific routines. The numerical efficiency of the generalized analytical Jacobian, however, is nearly 300 percent better than the reaction specific numerical Jacobian used in SPARK.

Carpenter, Mark H.↗