SEARCH · Engineering Papers
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.
Powerful, machine-independent program calculates theoretical thermodynamic properties of chemical systems. Aids in design of compressors, turbines, engines, heat exchangers, and chemical processing equipment.
This program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form (written in FORTRAN 2 for the IBM 7094 computer). A synchronized material property concept utilizing incremental time steps and the finite element matrix displacement approach has been adopted for the current analysis. A special option enables employment of constant time steps in the logarithmic scale, thereby reducing computational efforts resulting from accumulative material memory effects. A wide variety of structures with elastic or viscoelastic material properties can be analyzed by VISCEL. The program is written in FORTRAN 5 language for the Univac 1108 computer operating under the EXEC 8 system. Dynamic storage allocation is automatically effected by the program, and the user may request up to 195K core memory in a 260K Univac 1108/EXEC 8 machine. The physical program VISCEL, consisting of about 7200 instructions, has four distinct links (segments), and the compiled program occupies a maximum of about 11700 words decimal of core storage.
Computer program for calculating inviscid one-dimensional equilibrium, frozen and nonequilibrium nozzle expansion of propellant exhaust mixtures
Computer program calculates the inviscid one-dimensional equilibrium, frozen, and nonequilibrium nozzle expansion of propellant exhaust mixtures containing these six elements - carbon, hydrogen, oxygen, nitrogen, fluorine, and chlorine plus either aluminum, beryllium, boron or lithium. This program will perform calculations for contoured and conical nozzles.
Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.
This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.
We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.
We develop a non-cooperative game-theoretic model for the problem of graph minor-embedding to show that optimal compiling of adiabatic quantum programs in the sense of Nash equilibrium is possible.
This project extends and leverages the work completed under funding from ARPA-E and other sources by exploring aspects of the Cryogenic Carbon Capture™ (CCC) process and issues discovered during those projects, but that were not within the scope or budget of previous projects or that became apparent during field tests. Our successful completion of the tasks outlined in this project has increased the reliability, efficiency, and scalability of the CCC process and we are better prepared for pilot-scale testing. The first phase of this project improved key areas of the process through iterative design and experiment, culminating with recommendations for improvements to the existing skid-scale system developed under non-federal funding. The second phase of this project implemented the improved unit operations into the Cryogenic Carbon Capture External Cooling Loop (CCC ECL™) skid and field tests the improved system at a coal-fired power plant. The specific areas of work in this project are: Budget Period 1: 1) Selection of the most energy-efficient drying system for ambient-pressure flue gas. Sustainable Energy Solutions, LLC (SES) has successfully demonstrated a dryer that removes bulk water using a low-temperature brine to remove about 99% of the water followed by a mol sieve 3A desiccant dryer to remove the remaining 1% water; 2) Selection and demonstration of efficient solution to fouling due to dissolved CO 2 . SES has tested several commercial heat exchangers that purport to be able to handle fouling issues similar to those experienced in the contact liquid chiller. Replicated tests using the CCC ECL™ skid and dimpled brazed plate heat exchangers demonstrate that quick throttling of refrigerant combined with a burst of clean gas mitigate the fouling caused by CO 2 solubility in contact liquid with only minor increases in parasitic load. SES models indicate that parasitic load would increase by approximately 0.0125 GJ/tonne, an order of magnitude less than what was determined to be the maximum allowable additional parasitic load caused by this system; 3) Selection and demonstration of efficient and scalable solid–liquid separation system. The screw press technology met or exceeded the 70% solids loading benchmark during several test runs exceeding 24 hours on the CCC ECL™ skid. The average solids loading during the entire 50-hour run was 67%. The energy demand for this system is less than 0.1 GJ e /tonne of CO 2 processed; 4) Optimize design and demonstration of patented desublimating heat exchangers. We designed and built a hybrid spray tower that was tested over many hundreds of hours and a 50-hour continuous test run was completed capturing over 0.5 tonne CO 2 /day at >90% capture; (5) Plan for safe dilution of oxygen-depleted, cleaned flue-gas stream. The AERSCREEN simulator developed by the EPA modeled the dispersion of oxygen-depleted flue gas from the CCC process at a commercial plant. The results indicate that this plan effectively mitigates the flue gas hazards; (6) Simulation of and experimental work related to multi-pollutant capture. SES programmed an upgrade to its solid/vapor/liquid equilibrium (SVLE) algorithm. This enables simultaneous modeling of multiple pollutants (i.e., CO 2 , SO x , NO x , Hg, HCl). Budget Period 2; (7) Modify the CCC ECL™ skid for field testing and complete shakedown testing. The plans for modifying the CCC ECL™ skid were finalized and put into place. Shakedown testing of the updated CCC ECL™ skid occurred at SES using closed-loop mode (i.e., using a synthetic flue gas consisting of only N 2 and CO 2 ) and at the Hunter Power Plant using flue gas from the coal-fired boiler. Shakedown testing took longer than expected due to delays in receiving equipment and unexpected equipment performance in the field; (8) Operate the CCC ECL™ skid in the field. Longer testing of the CCC ECL™ skid occurred at the Hunter Power Plant. Multiple tests exceeding 35 hours were completed, with a total of over 600 hours of CO 2 capture completed using flue gas at the Hunter Power Plant with an average CO 2 capture rate of over 91% during those hours. Extensive notes and documentation were acquired informing improved methods for controls, startup, operation, and shut down; (9) Determine figures of merit and quantify process improvement. SES defined figures of merit for unit operation performance and overall CCC ECL™ process performance. All figures of merit were achieved or surpassed, except for 500 continuous hours of demonstration at the Hunter Power Plant; (10) Complete an Environmental Health and Safety Risk assessment by a third party. Tri-State completed an Environmental Health and Safety Risk assessment; (11) Update techno-economic analysis with third-party validation. The techno-economic analysis has been completed. It confirms the potential large reduction in energy and cost of the CCC process compared with leading alternatives and that the estimates and estimation techniques are reliable.
Program performs calculations such as chemical equilibrium for assigned thermodynamic states, theoretical rocket performance for both equilibrium and frozen compositions during expansion, incident and reflected shock properties, and Chapman-Jouget detonation properties. Features include simplicity of input and storage of all thermodynamic and transport property data on master tape.
A second order numerical method employing reference plane characteristics has been developed for the calculation of geometrically complex three dimensional nozzle-exhaust flow fields, heretofore uncalculable by existing methods. The nozzles may have irregular cross sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. The nozzles may have highly nonuniform entrance conditions, the medium considered being an equilibrium hydrogen-air mixture. The program calculates and carries along the underexpansion shock and contact as discrete discontinuity surfaces, for a nonuniform vehicle external flow.
The properties of clustering alloy systems and the manner in which they are influenced by rapid quenching from a containerless undercooled melt are discussed. It was postulated that rapid quenching under such conditions would result in highly disordered metastable alloys, and furthermore, that alloys in such conditions would possess physical properties characteristically different from those of alloys in the annealed equilibrium state. The scope of the program is essentially to gauge the influence of containerless undercooling on the submicrostructure of clustering-type alloys, using certain physical properties as diagnostic tools. Microstructures and macrostructures were to be examined using optical- and scanning-electron microscopy.
Surfactant-induced Marangoni effects strongly alter the stresses exerted along fluid particle interfaces. In low gravity processes, these stresses can dictate the system behavior. The dependence of Marangoni effects on surfactant physical chemistry is not understood, severely impacting our ability to predict and control fluid particle flows. A droplet in an extensional flow allows the controlled study of stretching and deforming interfaces. The deformations of the drop allow both Marangoni stresses, which resist tangential shear, and Marangoni elasticities, which resist surface dilatation, to develop. This flow presents an ideal model system for studying these effects. Prior surfactant-related work in this flow considered a linear dependence of the surface tension on the surface concentration, valid only at dilute surface concentrations, or a non-linear framework at concentrations sufficiently dilute that the linear approximation was valid. The linear framework becomes inadequate for several reasons. The finite dimensions of surfactant molecules must be taken into account with a model that includes surfaces saturation. Nonideal interactions between adsorbed surfactant molecules alter the partitioning of surfactant between the bulk and the interface, the dynamics of surfactant adsorptive/desorptive exchange, and the sensitivity of the surface tension to adsorbed surfactant. For example, cohesion between hydrocarbon chains favors strong adsorption. Cohesion also slows the rate of desorption from interfaces, and decreases the sensitivity of the surface tension to adsorbed surfactant. Strong cohesive interactions result in first order surface phase changes with a plateau in the surface tension vs surface concentration. Within this surface concentration range, the surface tension is decoupled from surface concentration gradients. We are engaged in the study of the role of surfactant physical chemistry in determining the Marangoni stresses on a drop in an extensional flow in a numerical and experimental program. Using surfactants whose dynamics and equilibrium behavior have been characterized in our laboratory, drop deformation will be studied in ground-based experiment. In an accompanying numerical study, predictive drop deformations will be determined based on the isotherm and equation of state determined in our laboratory. This work will improve our abilities to predict and control all fluid particle flows.
The main aim of this research program is to determine the ionization equilibrium and abundances in quasar outflows. Especially in the broad absorption line QSO PG 0946+301. We find that the outflow's metalicity is consistent with being solar, while the abundance ratio of phosphorus to other metals is at least ten times solar. These findings are based on diagnostics that are not sensitive to saturation and partial covering effects in the BALs (Broad Adsorption Lines), which considerably weakened previous claims for enhanced metalicity. Ample evidence for these effects is seen in the spectrum.
A computer program for calculating quasi-one-dimensional gas flow in axisymmetric and two-dimensional nozzles and rectangular channels is presented. Flow is assumed to start from a state of thermochemical equilibrium at a high temperature in an upstream reservoir. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. Electronic nonequilibrium effects can be included using a two-temperature model. An approximate laminar boundary layer calculation is given for the shear and heat flux on the nozzle wall. Boundary layer displacement effects on the inviscid flow are considered also. Chemical equilibrium and transport property calculations are provided by subroutines. The code contains precoded thermochemical, chemical kinetic, and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It provides calculations of the stagnation conditions on axisymmetric or two-dimensional models, and of the conditions on the flat surface of a blunt wedge. The primary purpose of the code is to describe the flow conditions and test conditions in electric arc heated wind tunnels.
The NATA code is a computer program for calculating quasi-one-dimensional gas flow in axisymmetric nozzles and rectangular channels, primarily to describe conditions in electric archeated wind tunnels. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. The shear and heat flux on the nozzle wall are calculated and boundary layer displacement effects on the inviscid flow are taken into account. The program contains compiled-in thermochemical, chemical kinetic and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It calculates stagnation conditions on axisymmetric or two-dimensional models and conditions on the flat surface of a blunt wedge. Included in the report are: definitions of the inputs and outputs; precoded data on gas models, reactions, thermodynamic and transport properties of species, and nozzle geometries; explanations of diagnostic outputs and code abort conditions; test problems; and a user's manual for an auxiliary program (NOZFIT) used to set up analytical curvefits to nozzle profiles.
The Python-based program is a laboratory automation tool designed to control and monitor electrochemical systems. The tool was developed for capacitive deionization (CDI) experiments, but it can be used for any system that requires controlled voltage or current segments and multi-parameter monitoring. The program integrates hardware components to run user-defined experimental parameters, providing operational control of a programmable power supply, peristaltic pump, and data acquisition devices. Currently, the program is structured with a workflow that includes an initialization (or pre-run) phase, a main loop, and a post-experiment stabilization (or post-run) phase. The initialization phase prepares and stabilizes the cell, ensuring that the electrodes and solution reach a baseline state before the experiment begins. The main loop consists of multiple voltage segments that repeat, controlling the experiment while recording key parameters such as time, voltage, current, pH, and conductivity. Finally, the post-experiment stabilization phase allows the system to stabilize after the experiment, returning the cell and solution to equilibrium conditions before ending the sequence. The program is designed with four variations, each tailored to different experimental needs. All variations include both the initialization and post-experiment stabilization stages, which run for a set amount of time, voltage, current, and flow rate before and after the main experiment block. The main loop runs for a set number of cycles, as defined by the user input, and each cycle is composed of 2 or 4 segments. The 4 program variations are described as follows: Program 1: The main program includes 2 segments. Each segment is defined to have a set duration, flow rate, voltage, and current. This program measures conductivity, flow rate, voltage, and current. Program 2: The main program expands Program 1 to include 4 segments. Each segment has a specified duration, flow rate, voltage, and current. Like Program 1, it measures conductivity, flow rate, voltage, and current. Program 3: The main program consists of 2 segments, each defined by time, flow rate, voltage, and current. In addition to conductivity, flow rate, voltage, and current, Program 3 collects pH and temperature data through a 4-channel data acquisition device. Program 4: This program independently controls two channels of a multi-channel power supply simultaneously. While conductivity can only be measured for one cell at a time, the dual-channel control makes it possible to operate two cells simultaneously under different voltage/current conditions. The main program includes 2 segments.For each program, all measurements are automatically logged and integrated into a single Excel output file. Data are displayed in numerical format and plotted, both in real time, to track system performance. A key feature of the program is its ability to synchronize all outputs so that every measurement shares a single timestamp, ensuring accurate alignment of voltage, current, pH, conductivity, and pH data.By combining hardware control, real-time monitoring, and unified data collection, this program significantly reduces manual workload and minimizes errors, making it a reliable platform for researchers, engineers, and laboratory technicians conducting CDI experiments, among other electrochemical tests.
The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.