Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 829 records · Page 46

Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state

We develop a flexible, relativistically covariant parametrization of the dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code smash, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be used to study critical behavior in dense nuclear matter, both at and away from equilibrium. We also show that two-particle correlations calculated from hadronic transport simulation data follow theoretical expectations based on the second-order cumulant ratio, and constitute a clear signature of the crossover region above the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

36 MATERIALS SCIENCE↗

Band Heater Heat Flux Characterization Using Inverse Heat Conduction Problem Models

Abstract An experimental apparatus was constructed to correlate water flowrate and temperature rise under an external band heater. Due to the physical characteristics of the band heater, its transient heating behavior is unknown. This paper investigates the application of inverse heat conduction problem (IHCP) methods to characterize the heat flux from the band heater. Three experiments with different heating times (5, 10, and 20 s) and no flowrate were conducted to measure the transient temperature under the 400 W band heater. Type-T thermocouples measure surface temperature at the centerline of the band heater. The experimental results are computed with five different heat conduction models. The models are chosen to identify how the heat flux response varies from a simplified to a realistic model. Additionally, the results of the experimental heat flux are compared to the manufacturer band heater data (58.9 kW/m2) for each model. The minimum time needed for the heater to fully energize the system is from 10 to 12 s. The residuals for each model are analyzed and used to evaluate the appropriateness of the five different models. The results show that the use of simpler models can achieve results similar to those of complex models, with less time and computational cost.

Engineering↗

Thermodynamics-informed latent space dynamics identification

This software showcases a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.

Cheung, Siu Wun↗

Advanced Cost-Effective Coal-Fired Rotating Detonation Combustor for High Efficiency Power Generation

Coal-based detonation as a means of heat addition in power generation devices is a revolutionary “out of the box” technology concept, wherein it is theoretically possible to achieve an increase in total pressure across the combustor as opposed to a loss of available energy as is the case in conventional combustion systems. Detonation is a revolutionary technology concept of pressure gain combustion that exploits pressure rise from the combustion process to augment high flow momentum and pressures. Recognizing that pressure gain combustion has become a topic of elevated national interest, high energetic modes (exergy) pressure gain combustion (PGC) systems are required. Pressure gain combustion is a key enabling technology for maintaining technological superiority in the development of advanced power generation systems. Pressure gain combustion is an innovative scheme of turbulent combustion that considerably increases thermodynamic cycle efficiencies (~10-20%). This project extends the State-of-the-Art by providing the first direct experimental and computational data of operability dynamics and performance characterization data and create the first validated Coal-Fired RDC investigations. Coal-based RDC is an innovative revolutionary technology that would result in high efficiency power generation and reduced of harmful emissions. Also, the fundamental and open source nature of the project will benefit a broader scientific and industry audience who will be able to utilize the outcomes of this project for advancing the coal-fired RDC power generation development.

01 COAL, LIGNITE, AND PEAT↗

Introduction to Finite Element Analysis: Abaqus Tutorial [Slides]

Abaqus is a suite of powerful engineering simulation programs, based on the finite element method, that can solve a variety of problems, ranging from simple to complex, including the fields of: Structural mechanics, Structural dynamics, Thermodynamics, Heat transfer, and more

42 ENGINEERING↗

The aluminum I autoionization doublet in the quiet solar spectrum

Observations are presented of the Al I autoionization doublet 1932 A and 1936 A in the quiet solar spectrum, obtained from the NRL slit spectrograph aboard Skylab and from the University of Hawaii Echelle Rocket Spectrograph. The observed profiles are compared with theoretical spectra computed for the Harvard Smithsonian Reference Atmosphere and the Vernazza, Avrett and Loeser (1976) solar models. It is found that nonlocal thermodynamic equilibrium effects are important in the line-formation problem and the synthetic spectra are in good agreeement with the data.

Heasley, J. N.↗

Evaluation of on-board hydrogen storage methods f or high-speed aircraft

Hydrogen is the fuel of choice for hypersonic vehicles. Its main disadvantage is its low liquid and solid density. This increases the vehicle volume and hence the drag losses during atmospheric flight. In addition, the dry mass of the vehicle is larger due to larger vehicle structure and fuel tankage. Therefore it is very desirable to find a fuel system with smaller fuel storage requirements without deteriorating the vehicle performance substantially. To evaluate various candidate fuel systems, they were first screened thermodynamically with respect to their energy content and cooling capacities. To evaluate the vehicle performance with different fuel systems, a simple computer model is developed to compute the vehicle parameters such as the vehicle volume, dry mass, effective specific impulse, and payload capacity. The results indicate that if the payload capacity (or the gross lift-off mass) is the most important criterion, only slush hydrogen and liquid hydrogen - liquid methane gel shows better performance than the liquid hydrogen vehicle. If all the advantages of a smaller vehicle are considered and a more accurate mass analysis can be performed, other systems using endothermic fuels such as cyclohexane, and some boranes may prove to be worthy of further consideration.

Akyurtlu, Ates↗

Multi-Fidelity Simulation of a Turbofan Engine With Results Zoomed Into Mini-Maps for a Zero-D Cycle Simulation

A Zero-D cycle simulation of the GE90-94B high bypass turbofan engine has been achieved utilizing mini-maps generated from a high-fidelity simulation. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled 3D computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the 3D component models are integrated into the cycle model via partial performance maps generated from the CFD flow solutions using one-dimensional mean line turbomachinery programs. This paper highlights the generation of the high-pressure compressor, booster, and fan partial performance maps, as well as turbine maps for the high pressure and low pressure turbine. These are actually "mini-maps" in the sense that they are developed only for a narrow operating range of the component. Results are compared between actual cycle data at a take-off condition and the comparable condition utilizing these mini-maps. The mini-maps are also presented with comparison to actual component data where possible.

Turner, Mark G.↗

Computational and Experimental Investigation of Li-doped Ionic Liquid Electrolytes: [pyr14][tfsi], [pyr13][fsi], and [EMIM][BF4]

We employ molecular dynamics (MD) simulation and experiment to investigate the structure, thermodynamics, and transport of N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N -methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt mole fraction (0.05 xLi+ 0.33) and temperature (298 K T 393 K). Structurally, Li+ is shown to be solvated by three anion neigh- bors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi+ we find the presence of lithium aggregates. Pulsed field gradient spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, diffusion, viscosity, and ionic conductivity) are obtained with MD and show excellent agreement with experiment. While the diffusion exhibits a systematic decrease with increasing xLi+, the contribution of Li+ to ionic conductivity increases until reach- ing a saturation doping level of xLi+ 0.10. Comparatively, the Li+ conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. Our transport results also demonstrate the necessity of long MD simulation runs ( 200 ns) required to converge transport properties at room T. The differences in Li+ transport are reflected in the residence times of Li+ with the anions (Li), which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, to comment on the relative kinetics of Li+ transport in each liquid, we find that while the net motion of Li+ with its solvation shell (vehicular) significantly contributes to net diffusion in all liquids, the importance of transport through anion exchange (hopping) increases at high xLi+ and in liquids with large anions.

electrolytes↗

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical study to assess sulfur hexafluoride as a medium for testing multielement airfoils

A methodology is described for computing viscous flows of air and sulfur hexafluoride (SF6). The basis is an existing flow solver that calculates turbulent flows in two dimensions on unstructured triangular meshes. The solver has been modified to incorporate the thermodynamic model for SF6 and used to calculate the viscous flow over two multielement airfoils that have been tested in a wind tunnel with air as the test medium. Flows of both air and SF6 at a free-stream Mach number of 0.2 and a Reynolds number of 9 x 10(exp 6) are computed for a range of angles of attack corresponding to the wind-tunnel test. The computations are used to investigate the suitability of SF6 as a test medium in wind tunnels and are a follow-on to previous computations for single-element airfoils. Surface-pressure, lift, and drag coefficients are compared with experimental data. The effects of heavy gas on the details of the flow are investigated based on computed boundary-layer and skin-friction data. In general, the predictions in SF6 vary little from those in air. Within the limitations of the computational method, the results presented are sufficiently encouraging to warrant further experiments.

Bonhaus, Daryl L.↗

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

Development of a Gibbs Energy Minimiser for the MOOSE-based Corrosion Modelling App Yellowjacket and Validation of MSTDB

Nuclear materials are highly complex multiscale, multiphysics systems,and an effective prediction of nuclear reactor performance and safety requires simulation capabilities that tightly couple different physical phenomena. The Idaho National Laboratory’s Multiphysics Object Oriented Simulation Environment (MOOSE) provides the computational foundation for performing such simulations. With the move towards advanced reactors, such as the Molten Salt Reactor (MSR), that employ high temperature fluids compared to conventional reactors, corrosion has become a problem of great interest. A new application called Yellowjacket is currently under development to directly couple thermodynamic equilibrium and kinetics with phase field models in order to model corrosion in MSRs. As part of Yellowjacket, a Gibbs energy minimiser is being developed to perform thermochemical equilibrium calculations for a range of different materials, which is currently in its infancy. This report describes the further progress towards the development of Yellowjacket Gibbs energy minimiser. Ontario Tech University is developing a new Gibbs energy minimiser for Yellowjacket which is the primary contribution of this work. The aim to develop a thermochemistry solver for the MOOSE framework following the same development philosophy and using the same tools and libraries. A special focus is on performance, documentation and SQA. Furthermore through a scope extension partway through the fiscal year a thorough assessment of the MSTDB-TC v1.3 was performed at Ontario Tech University in the context of continuous improvement and quality assurance. The objective of the work was to have an arm’s length review of the database to give confidence that the database is performing as it was intended while assessing its current state to give recommendations to future developments. This assessment involved two parts: A) a quantitative assessment, and B) a qualitative assessment. Part A involved developing an automated test-suite that would compute values from the database using Thermochimica with comparisons to experimental measurements for validation purposes, which gives confidence to the database’s stakeholders that its functioning properly. Part B involved reviewing all binary systems in the database and making a qualitative assessment with two performance indicators: comprehensiveness and overall confidence. It is important to note that the models in the database are empirical, which is to say that the quality of any model is highly dependent on the experimental data used to inform its development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-order Accurate Finite-Volume Formulations for the Pressure Gradient Force in Layered Ocean Models

Discretisation of the horizontal pressure gradient force in layered ocean models is a challenging task, with non-trivial interactions between the thermodynamics of the fluid and the geometry of the layers often leading to numerical difficulties. We present two new finite-volume schemes for the pressure gradient operator designed to address these issues. In each case, the horizontal acceleration is computed as an integration of the contact pressure force that acts along the perimeter of an associated momentum control-volume. A pair of new schemes are developed by exploring different control-volume geometries. Non-linearities in the underlying equation-of-state definitions and thermodynamic profiles are treated using a high-order accurate numerical integration framework, designed to preserve hydrostatic balance in a non-linear manner. Numerical experiments show that the new methods achieve high levels of consistency, maintaining hydrostatic and thermobaric equilibrium in the presence of strongly-sloping layer geometries, non-linear equations-of-state and non-uniform vertical stratification profiles. These results suggest that the new pressure gradient formulations may be appropriate for general circulation models that employ hybrid vertical coordinates andor terrain-following representations.

pressure gradient force↗

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗

eQuilibrator 3.0: a database solution for thermodynamic constant estimation

Abstract eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg2+ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms.

59 BASIC BIOLOGICAL SCIENCES↗

Theoretical study of thermodynamic properties and reaction rates of importance in the high-speed research program

One of the primary goals of NASA's high-speed research program is to determine the feasibility of designing an environmentally safe commercial supersonic transport airplane. The largest environmental concern is focused on the amount of ozone destroying nitrogen oxides (NO(x)) that would be injected into the lower stratosphere during the cruise portion of the flight. The limitations placed on NO(x) emission require more than an order of magnitude reduction over current engine designs. To develop strategies to meet this goal requires first gaining a fundamental understanding of the combustion chemistry. To accurately model the combustor requires a computational fluid dynamics approach that includes both turbulence and chemistry. Since many of the important chemical processes in this regime involve highly reactive radicals, an experimental determination of the required thermodynamic data and rate constants is often very difficult. Unlike experimental approaches, theoretical methods are as applicable to highly reactive species as stable ones. Also our approximation of treating the dynamics classically becomes more accurate with increasing temperature. In this article we review recent progress in generating thermodynamic properties and rate constants that are required to understand NO(x) formation in the combustion process. We also describe our one-dimensional modeling efforts to validate an NH3 combustion reaction mechanism. We have been working in collaboration with researchers at LeRC, to ensure that our theoretical work is focused on the most important thermodynamic quantities and rate constants required in the chemical data base.

Langhoff, Stephen↗