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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.

Maximum Entropy↗
NASA NTRS2000

Aerothermodynamic Analysis of Stardust Sample Return Capsule with Coupled Radiation and Ablation

An aerothermodynamic analysis of the forebody aeroshell of the Stardust Sample Return Capsule is carried out by using the axisymmetric viscous shock-layer equations with and without fully coupled radiation and ablation. Formulation of the viscous shock-layer equations with shoulder radius as the length scale and implementation of the Vigneron pressure condition allow resolution of the flowfield over the shoulder. With a predominantly supersonic outflow over the shoulder, a globally iterated solution or viscous shock-layer equations can be obtained. The stagnation-point results are obtained along a specified trajectory, whereas detailed calculations along the body are provided at the peak-heating point. The equilibrium calculations with ablation injection are the focus of the present study because of the lack of a general chemical nonequilibrium analysis that accounts for both surface and flowfield effect. The equilibrium calculations also provide a simple way to conserve surface (and flowfield) elemental composition for the current small ablation injection rates, where the surface elemental composition is a mixture of freestream and ablator elements. Therefore, the coupled laminar and turbulent flow solutions with radiation and ablation are obtained by using the equilibrium flow chemistry, whereas a nonequilibrium chemistry model is used for solutions without ablation and turbulence. Various computed results are compared with those obtained by the other researchers.

Gupta, Roop N.↗
NASA NTRS1995

A family of dynamic models for large-eddy simulation

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Carati, D.↗
NASA NTRS1989

Plotting component maps in the Navy/NASA Engine Program (NNEP): A method and its usage

The Navy/NASA Engine Program (NNEP) and the new extended version which handles chemical equilibrium (NNEPEQ) are very general cycle analysis codes that have been used extensively to calculate design and off-design performance of a wide range of turbine engine cycles and configurations. Component maps are used to obtain the off-design engine performance and a matched engine cycle. This paper describes a method of plotting the scaled NNEP compressor and turbine maps as the user runs the NNEP code as well as plotting the operating line defined by all the cases that were computed in that particular NNEP run. Afterwards, an example demonstrates the use of this capability to help analyze an engine cycle model and then make improvements to that cycle.

Plencner, Robert M.↗
NASA NTRS2012

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗
DOE OSTI2021

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
DOE OSTI2026

Understanding and Controlling Photoexcited Molecules in Complex Environments (Final Technical Report)

Predicting the outcome of a photochemical reaction is complicated by factors outside the scope of tradition theoretical chemistry. On such small scales, fluctuations of the environment render outcomes stochastic, nonadiabatic dynamics blur the separation between nuclear and electronic motion, and ultrafast relaxation renders typical assumptions of equilibrium inappropriate. This project aimed to elucidate and control nonadiabatic molecular dynamics in condensed phases. We developed accurate, efficient simulation tools to study irreversible molecular processes computational across a range of scales. We applied these methodological advances to study electron transfer under exotic conditions where spin and topology opened new pathways for reactivity. We also brought these ideas to bare on light harvesting systems, where energy and charge transport are coupled, and where long lifetimes can be deleterious to function. These advances helped to establish general principles for chemical efficiency and guide the design of nanoscale energy materials, molecular machines, and related technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Capturing the run-in of a pebble-bed reactor by using thermal feedback and high-fidelity neutronics simulations

Modeling the run-in of a pebble-bed reactor (PBR) can be challenging as a result of changes in the power, fuel type, and temperatures that occur throughout the run-in period. Previous work utilized high-fidelity neutronics simulations or lower-fidelity coupled neutronics/thermal-hydraulics models to capture the general characteristics of the run-in process. Here, the present work employs high-fidelity neutronics simulations (using Serpent) coupled with thermal-hydraulics simulations (using Griffin–Pronghorn) to capture the thermal feedback present during the run-in and approach to equilibrium for a PBR. Incorporating thermal feedback enables important distinctions to be made about conditions occurring inside the core, as the power distribution, discharge burnup, and isotopic compositions are all affected by the temperature distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A parametric approach to identify synergistic domains of process intensification for reactive separation

Process intensification aims to combine multiple tasks within multi-functional units to drastically improve economic, energy or sustainability metrics of a chemical process. Limited work exists to systematically identify the synergistic domains where intensification outperforms its nonintensified counterpart. In this work, we computationally derive the synergistic domains of a reactive separation system. Specifically, we first postulate general models for both intensified and nonintensified systems. We use these models to generate data to train a ReLU-type artifical neural network (ANN). Further, the trained ReLU-NN model is formulated as a multi-parametric mixed-integer linear program (mp-MILP), and the critical regions of this mp-MILP define the synergistic feasible domains of intensification. We have derived these synergistic domains of vapor–liquid equilibrium (VLE)-based reactive separation for several industrial applications. These synergistic domains enable quick screening of properties that favor intensification.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗
DOE OSTI2024

Quantum Thermodynamics of Nonequilibrium Processes in Lattice Gauge Theories

A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Using Multiplicative Algebraic Reconstruction Techniques (MART) to Derive Instrument Requirements for Computed Tomography Imaging Spectrographs

Computed Tomography Imaging Spectrographs (CTISs), which are generally slitless or large-aperture spectrographs that observe a wide, dispersed field of view in multiple dispersion angles or diffraction orders, provide a unique opportunity to capture spectral information over an extended source, such as the Sun, but require software reconstruction techniques to be fully utilized. A useful data product that can be recovered from CTIS observations is intensity maps of solar features in single spectral lines. These intensity maps can then be used to determine temperature, density, abundance, and equilibrium properties of the emitting plasma. We apply the multiplicative algebraic reconstruction technique (MART) to example data from a variety of CTIS configurations to determine the capability and limitations of the method to return spectrally-pure maps. By completing this study, we aim to establish a path to derive requirements for CTIS instruments that rely on reconstruction techniques to meet their science objectives.

Amy Winebarger↗
NASA NTRS1982

Numerical techniques in radiative heat transfer for general, scattering, plane-parallel media

The study of radiative heat transfer with scattering usually leads to the solution of singular Fredholm integral equations. The present paper presents an accurate and efficient numerical method to solve certain integral equations that govern radiative equilibrium problems in plane-parallel geometry for both grey and nongrey, anisotropically scattering media. In particular, the nongrey problem is represented by a spectral integral of a system of nonlinear integral equations in space, which has not been solved previously. The numerical technique is constructed to handle this unique nongrey governing equation as well as the difficulties caused by singular kernels. Example problems are solved and the method's accuracy and computational speed are analyzed.

Sharma, A.↗
DOE OSTI2024

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Here, a new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Disentangling the impact of quasiparticles and two-level systems on the statistics of superconducting-qubit lifetime

Temporal fluctuations in the superconducting qubit lifetime, T 1 , present additional challenges in the pursuit of fault-tolerant quantum computing. Although the exact mechanisms remain unclear, T 1 fluctuations are generally attributed to strong coupling between the qubit and a few near-resonant two-level systems (TLSs), which can exchange energy with an ensemble of thermally fluctuating two-level fluctuators (TLFs) at low frequencies. Here, we report T 1 measurements of qubits with varying geometrical footprints and surface dielectrics as a function of temperature. By analyzing the noise spectrum of the qubit depolarization rate, Γ 1 = 1 / T 1 , we disentangle the contributions of TLSs, nonequilibrium quasiparticles (QPs), and equilibrium (thermally excited) QPs to the variance in Γ 1 . We find that the Γ 1 variance in qubits with smaller footprints is more susceptible to QP and TLS fluctuations than that in larger-footprint qubits. Furthermore, the QP-induced variances in all qubits align with the theoretical framework of QP diffusion and fluctuation. These findings offer valuable insights for future qubit design and engineering optimization.

Zhu, Shaojiang [Fermilab] (ORCID:0000000293180092)↗
NASA NTRS1995

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.↗
DOE OSTI2021

A first-principles study of hydrogen surface coverage on δ-Pu (100), (111), and (110) surfaces

Hydriding corrosion of plutonium leads to surface cracking, pitting, and ultimately structural failure. Laboratory experiments demonstrate that hydriding begins on the surface or near the subsurface of plutonium. However, there has not yet been a systematic evaluation of hydrogen surface coverage on plutonium. In this work, we compute the surface energies of the low facet surfaces of face-centered cubic δ-Pu. The adsorption free energies of expected hydrogen structures at low and high coverage are presented along with the likely progression for filling sites as the H 2 partial pressure increases. Implications for near-equilibrium pressure hydride nucleation and non-equilibrium millibar pressure hydriding are discussed.

36 MATERIALS SCIENCE↗