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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Analytic gradients for coupled-cluster energies that include noniterative connected triple excitations - Application to cis- and trans-HONO

An efficient formulation of the analytic energy gradient for the single and double excitation coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations is presented. The formulation has a small computational cost, and the algebraic manipulations may be applied generally to the analytic gradient of Moller-Plesset perturbation theory energies. The new formulation has been implemented in an efficient set of programs that utilize highly vectorized algorithms and has been used to investigate the equilibrium structures, harmonic vibrational frequencies, IR intensities, and energy separation of cis- and trans-HONO.

Lee, Timothy J.↗

DOE OSTI2024

Denoising and Extension of Response Functions in the Time Domain

Response functions of quantum systems, such as electron Green’s functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Further, response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗

NASA NTRS2003

Computational Aeroheating Predictions for Mars Lander Configurations

The proposed Mars Science Laboratory (MSL) mission is intended to deliver a large rover to the Martian surface within 10 km of the target site. This paper presents computational fluid dynamics (CFD) predictions of forebody heating rates for two MSL entry configurations with fixed aerodynamic trim tabs. Results are compared to heating on a 70-deg sphere-cone reference geometry. All three heatshield geometries are designed to trim hypersonically at a 16 deg angle of attack in order to generate the lift-to-drag ratio (LID) required for precision landing. Comparisons between CFD and tunnel data are generally in good agreement for each configuration, but the computations predict more flow separation and higher heating on a trim tab inclined 10 deg relative to the surface. CFD solutions at flight conditions were obtained using an 8-species Mars gas in chemical and thermal non-equilibrium. Laminar and Baldwin-Lomax solutions were used to estimate the effects of the trim tabs and turbulence on heating. A tab extending smoothly from the heatshield flank is not predicted to increase laminar or turbulent heating rates above the reference levels. Laminar heating on a tab deflected 10 deg from the conical heatshield is influenced by flow separation and is up to 35% above the baseline heating rate. The turbulent solution on the inclined tab configuration predicts attached flow and a 43% heating increase above the reference level.

Edquist, Karl T.↗

NASA NTRS2011

Technique for Calculating Solution Derivatives With Respect to Geometry Parameters in a CFD Code

A solution has been developed to the challenges of computation of derivatives with respect to geometry, which is not straightforward because these are not typically direct inputs to the computational fluid dynamics (CFD) solver. To overcome these issues, a procedure has been devised that can be used without having access to the mesh generator, while still being applicable to all types of meshes. The basic approach is inspired by the mesh motion algorithms used to deform the interior mesh nodes in a smooth manner when the surface nodes, for example, are in a fluid structure interaction problem. The general idea is to model the mesh edges and nodes as constituting a spring-mass system. Changes to boundary node locations are propagated to interior nodes by allowing them to assume their new equilibrium positions, for instance, one where the forces on each node are in balance. The main advantage of the technique is that it is independent of the volumetric mesh generator, and can be applied to structured, unstructured, single- and multi-block meshes. It essentially reduces the problem down to defining the surface mesh node derivatives with respect to the geometry parameters of interest. For analytical geometries, this is quite straightforward. In the more general case, one would need to be able to interrogate the underlying parametric CAD (computer aided design) model and to evaluate the derivatives either analytically, or by a finite difference technique. Because the technique is based on a partial differential equation (PDE), it is applicable not only to forward mode problems (where derivatives of all the output quantities are computed with respect to a single input), but it could also be extended to the adjoint problem, either by using an analytical adjoint of the PDE or a discrete analog.

Mathur, Sanjay↗

DOE OSTI2024

Integrating subsystem embedding subalgebras and coupled cluster Green’s function: a theoretical foundation for quantum embedding in excitation manifold

Here, in this study, we introduce a novel approach to coupled-cluster Green's function (CCGF) embedding by seamlessly integrating conventional CCGF theory with the state-of-the-art sub-system embedding sub-algebras coupled cluster (SES-CC) formalism. This integration focuses primarily on delineating the characteristics of the sub-system and the corresponding segments of the Green's function, defined explicitly by active orbitals. Crucially, our work involves the adaptation of the SES-CC paradigm, addressing the left eigenvalue problem through a distinct form of Hamiltonian similarity transformation. This advancement not only facilitates a comprehensive representation of the interaction between the embedded sub-system and its surrounding environment but also paves the way for the quantum mechanical description of multiple embedded domains, particularly by employing the emergent quantum flow algorithms. Our theoretical underpinnings further set the stage for a generalization to multiple embedded sub-systems. This expansion holds significant promise for the exploration and application of non-equilibrium quantum systems, enhancing the understanding of system–environment interactions. In doing so, the research underscores the potential of SES-CC embedding within the realm of quantum computations and multi-scale simulations, promising a good balance between accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Chemically Reacting One-Dimensional Gas-Particle Flows

The governing equations for the one-dimensional flow of a gas-particle system are discussed. Gas-particle effects are coupled via the system momentum and energy equations with the gas assumed to be chemically frozen or in chemical equilibrium. A computer code for calculating the one-dimensional flow of a gas-particle system is discussed and a user's input guide presented. The computer code provides for the expansion of the gas-particle system from a specified starting velocity and nozzle inlet geometry. Though general in nature, the final output of the code is a startline for initiating the solution of a supersonic gas-particle system in rocket nozzles. The startline includes gasdynamic data defining gaseous startline points from the nozzle centerline to the nozzle wall and particle properties at points along the gaseous startline.

Tevepaugh, J. A.↗

NASA NTRS1980

Application of finite-element-techniques to the interaction of conduction and radiation in an absorbing, scattering and emitting medium

In this paper, the authors demonstrate that a Galerkin finite element method of analysis, utilizing isoparametric elements, offers a viable means of solving continuum thermal radiation problems with conduction in a participating medium. The participating medium was considered to be a gray radiation medium exhibiting isotropic absorption, emission, and scattering characteristics and optical properties that are independent of temperature. The medium was considered to be bounded by infinite parallel opaque, gray surfaces with diffuse emission and reflection characteristics. In solving this problem, a finite element formulation was developed to describe a system in radiative equilibrium. Then the results of this first analysis were linked with a second finite element model which incorporated conduction into the analysis. The results of this study were found to be in good agreement with existing published data. The model offers the following advantageous features: geometric generality, a computational algorithm which is 'convenient' and 'computable', and a functional basis for extension of the radiation model to higher order approximation.

Wu, S. T.↗

NASA NTRS1995

Structural design using equilibrium programming formulations

Solutions to increasingly larger structural optimization problems are desired. However, computational resources are strained to meet this need. New methods will be required to solve increasingly larger problems. The present approaches to solving large-scale problems involve approximations for the constraints of structural optimization problems and/or decomposition of the problem into multiple subproblems that can be solved in parallel. An area of game theory, equilibrium programming (also known as noncooperative game theory), can be used to unify these existing approaches from a theoretical point of view (considering the existence and optimality of solutions), and be used as a framework for the development of new methods for solving large-scale optimization problems. Equilibrium programming theory is described, and existing design techniques such as fully stressed design and constraint approximations are shown to fit within its framework. Two new structural design formulations are also derived. The first new formulation is another approximation technique which is a general updating scheme for the sensitivity derivatives of design constraints. The second new formulation uses a substructure-based decomposition of the structure for analysis and sensitivity calculations. Significant computational benefits of the new formulations compared with a conventional method are demonstrated.

Scotti, Stephen J.↗

NASA NTRS1999

Non-Invasive Investigation of Bone Adaptation in Humans to Mechanical Loading

Experimental studies have identified peak cyclic forces, number of loading cycles, and loading rate as contributors to the regulation of bone metabolism. We have proposed a theoretical model that relates bone density to a mechanical stimulus derived from average daily cumulative peak cyclic 'effective' tissue stresses. In order to develop a non-invasive experimental model to test the theoretical model we need to: (1) monitor daily cumulative loading on a bone, (2) compute the internal stress state(s) resulting from the imposed loading, and (3) image volumetric bone density accurately, precisely, and reproducibly within small contiguous volumes throughout the bone. We have chosen the calcaneus (heel) as an experimental model bone site because it is loaded by ligament, tendon and joint contact forces in equilibrium with daily ground reaction forces that we can measure; it is a peripheral bone site and therefore more easily and accurately imaged with computed tomography; it is composed primarily of cancellous bone; and it is a relevant site for monitoring bone loss and adaptation in astronauts and the general population. This paper presents an overview of our recent advances in the areas of monitoring daily ground reaction forces, biomechanical modeling of the forces on the calcaneus during gait, mathematical modeling of calcaneal bone adaptation in response to cumulative daily activity, accurate and precise imaging of the calcaneus with quantitative computed tomography (QCT), and application to long duration space flight.

Whalen, R.↗

DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Quantum algorithms from fluctuation theorems: Thermal-state preparation

Fluctuation theorems provide a correspondence between properties of quantum systems in thermal equilibrium and a work distribution arising in a non-equilibrium process that connects two quantum systems with Hamiltonians H 0 and H 1 = H 0 + V. Building upon these theorems, we present a quantum algorithm to prepare a purification of the thermal state of H 1 at inverse temperature β ≥ 0 starting from a purification of the thermal state of H 0 . The complexity of the quantum algorithm, given by the number of uses of certain unitaries, is $\mathcal{O}$ (e β(ΔA - w l )/2 ), where ΔA is the free-energy difference between H 1 and H 0 , and w l is a work cutoff that depends on the properties of the work distribution and the approximation error ϵ > 0. If the non-equilibrium process is trivial, this complexity is exponential i β∥V∥, where ∥V∥ is the spectral norm of V. This represents a significant improvement of prior quantum algorithms that have complexity exponential in β∥H 1 ∥ in the regime where ∥V∥$\ll$ ∥H 1 ∥. The dependence of the complexity in ϵ varies according to the structure of the quantum systems. It can be exponential in 1/ϵ in general, but we show it to be sublinear in 1/ϵ if H 0 and H 1 commute, or polynomial in 1/ϵ if H 0 and H 1 are local spin systems. The possibility of applying a unitary that drives the system out of equilibrium allows one to increase the value of w l and improve the complexity even further. To this end, we analyze the complexity for preparing the thermal state of the transverse field Ising model using different non-equilibrium unitary processes and see significant complexity improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Multi-cycle reload analysis of a long cycle gas-cooled fast modular reactor

There is currently significant interest in deploying HALEU-fueled fast reactors, including the General Atomics (GA) Fast Modular Reactor (FMR). Such reactors can achieve very long fuel cycles, but with multi-batch loading will take decades to reach equilibrium. This motivates design and analysis of both the initial core and multi-cycle reload, which is typically performed using fast-running, deterministic fast reactor codes such as the Argonne Reactor Computation (ARC) codes. In this paper, multicycle reload of the GA FMR is analyzed using the ARC codes. The GA FMR utilizes 19.75 % enriched fuel in a 16 year cycle with a three-batch strategy, with twice-burned fuel placed on the core periphery. The GA FMR has a softened neutron spectrum due to reflecting elements in the core, so the neutronic solution is first benchmarked against the OpenMC Monte Carlo code. Discrepancy on k eff is 400–600 pcm, likely due to the softened neutron spectrum, heterogeneous fuel assembly design and central reflector. However, the rms discrepancy on the assembly power distribution is only 0.6 %, despite the presence of the central reflector. A reload strategy is devised for the first three cycles of such a reactor, ultimately spanning the first 45–48 years of its operation. The fresh core uses 19.75 %, 19.25 % and 16.75 % enriched fuel in place of fresh, once-burned and twice-burned and is then subsequently refueled with only 19.75 % enriched fuel. The cycle length is varied over 3 cycles of operation to balance fuel utilization and reactor availability, specifically with use of an extended 18-year Cycle 1, followed by a shortened 11-year Cycle 2. Cycle 3 is close to the target 16-year length. Finally, placing twice burned assemblies next to the GA FMR central reflector can reduce power peaking by 3 %, at the expense of slightly reducing the cycle length.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment (FY2021--Part 1)

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The tasks completed for this report are as follows: First, Task 1d: The net radiation transfer method was implemented into MOOSE for modeling reactor cavity cooling system geometries. RCCS models for two experiments were created: (1) Natural Convection Shutdown Heat Removal Test Facility (NSTF) experiment R022, and (2) HTTR VCS mockup. For validation, computed temperature distributions were compared to measured temperatures. Next, Task 4c: An algorithm for computing the pebble bed reactor equilibrium core isotopic com-position was developed and an initial version is implemented into the reactor multi-physics code Griffin. Initial results for a simplified axisymmetric pebble bed reactor are presented. Finally, Task 7: generation of a reference plant model for molten salt cooled pebble bed reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2015

Spin-Up and Tuning of the Global Carbon Cycle Model Inside the GISS ModelE2 GCM

Planetary carbon cycle involves multiple phenomena, acting at variety of temporal and spacial scales. The typical times range from minutes for leaf stomata physiology to centuries for passive soil carbon pools and deep ocean layers. So, finding a satisfactory equilibrium state becomes a challenging and computationally expensive task. Here we present the spin-up processes for different configurations of the GISS Carbon Cycle model from the model forced with MODIS observed Leaf Area Index (LAI) and prescribed ocean to the prognostic LAI and to the model fully coupled to the dynamic ocean and ocean biology. We investigate the time it takes the model to reach the equilibrium and discuss the ways to speed up this process. NASA Goddard Institute for Space Studies General Circulation Model (GISS ModelE2) is currently equipped with all major algorithms necessary for the simulation of the Global Carbon Cycle. The terrestrial part is presented by Ent Terrestrial Biosphere Model (Ent TBM), which includes leaf biophysics, prognostic phenology and soil biogeochemistry module (based on Carnegie-Ames-Stanford model). The ocean part is based on the NASA Ocean Biogeochemistry Model (NOBM). The transport of atmospheric CO2 is performed by the atmospheric part of ModelE2, which employs quadratic upstream algorithm for this purpose.

Climate model↗

NASA NTRS1964

Discussion of the General Electric proposal to develop a computer program to select thermal coatings for passive temperature control of satellites

The General Electric Company, Missiles and Space Division, submitted/a proposal to NASA for the development of an IBM 7094 computer program which would select the exterior surface coatings for passively controlling spacecraft temperatures. G. E. claims that the "trial and error" procedures currently used can be accomplished more rationally and can therefore be programmed for a digital computer. In the ASME paper, 63-HT-41, which was presented at the ASME-AIChE Heat Transfer Conference at Boston in August, 1963, Costello, Harper, and Cline, described the procedures that have been used at G. E. to develop a coating selection program subject to the following restrictions: 1. Steady-state conditions prevail, 2. Heat transfer occurs by radiation only, 3. Temperatures are optimized at only one interior point in the spacecraft, an d 4. Only the solar absorptance of the external coatings is varied to optimize temperature. The emittance must remain constant at initially specified values. General Electric proposes to develop a generalized program in three steps: 1. The program would vary both solar absorptance and hemispherical emittance to obtain the optimum coating patterns; 2. Temperatures would be optimized at more than one interior point; and 3. The equations would be modified to account for both conduction and radiation heat transfer. In the development of the general program, the scope would be restricted to steady-state heat transfer. Since the thermal designs of most spacecraft are based primarily on nearly equilibrium conditions, the proposed program could have wide application. An obvious extension of the proposed program would be to account for transient temperatures.

IBM 7090 COMPUTER↗

NASA NTRS1996

A Parametric Study of a Plug Nozzle, Using the Liquid Propellant Program (LPP) Code

The Liquid Propellant Program (LPP) computer code is a super-set of the industry standard Two Dimensional Kinetics (TDK) computer code. The TDK code uses a two dimensional method of characteristics solution with fully coupled finite rate kinetics for axially symmetric nozzles. The chemical reactions are modeled with a generalized reaction package that include three dimensional body efficiencies and four reaction rate forms. The code performs optional solutions for frozen or equilibrium flow. TDK evaluates discrete shocks, both attached or induced. The Transonic module models variable mixture ratio profiles from the combustion chamber injector. The Mass Addition Boundary Layer module (MABL) calculates the boundary parameters with the same chemistry options, and includes transpiration or tangential slot injection of gas at the wall. The LPP upgrades include: planar nozzle, scarfed nozzles, plug nozzles, and scramjet nozzle configurations. The code evaluates both upper and lower wall flow simulation, and includes the interaction with the external flow. The MABL module evaluates equilibrium radiation heat transfer for both upper and lower walls. In addition, LPP code models combustion effects due to injector inefficiencies with the Spray Combustion Analysis Program (SCAP) module. The LPP package provides extensive post plotting capabilities for flow visualization. The LPP is sufficiently fast and robust to provide performance predictions for extensive parametric studies and sufficiently accurate to provide flow field and performance solutions for detailed studies.

Dunn, Stuart S.↗

DOE OSTI2023

Mechanical theory of nonequilibrium coexistence and motility-induced phase separation

Nonequilibrium phase transitions are routinely observed in both natural and synthetic systems. The ubiquity of these transitions highlights the conspicuous absence of a general theory of phase coexistence that is broadly applicable to both nonequilibrium and equilibrium systems. Here, in this study, we present a general mechanical theory for phase separation rooted in ideas explored nearly a half-century ago in the study of inhomogeneous fluids. The core idea is that the mechanical forces within the interface separating two coexisting phases uniquely determine coexistence criteria, regardless of whether a system is in equilibrium or not. We demonstrate the power and utility of this theory by applying it to active Brownian particles, predicting a quantitative phase diagram for motility-induced phase separation in both two and three dimensions. This formulation additionally allows for the prediction of novel interfacial phenomena, such as an increasing interface width while moving deeper into the two-phase region, a uniquely nonequilibrium effect confirmed by computer simulations. The self-consistent determination of bulk phase behavior and interfacial phenomena offered by this mechanical perspective provide a concrete path forward toward a general theory for nonequilibrium phase transitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗