Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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 937 records · Page 52

Vacillation, sudden warmings and potential enstrophy balance in the stratosphere

The balance of potential enstrophy and its relationship to vacillation cycles and the sudden warming is studied for a beta-channel model of the stratosphere. It is shown that the mean flow cannot be steady in the presence of large-amplitude quasi-geostrophic waves (approximately 1-0.25 geopotential kilometers /gpkm/) when any dissipation is present, and the maximum wave amplitude allowed is approximately 2 gpkm. If wave forcing (transience plus dissipation) is artificially maintained, the mean flow decelerates slowly at first then explosively as the potential vorticity gradient of the basic state is wiped out over the channel. This process is called wave saturation. The initial phase of the explosive deceleration resembles both the observed and modeled mean flow evolution during a sudden stratospheric warming. A simple vacillation model based upon these ideas shows remarkable similarity to the results of Holton and Mass (1976) and Davies (1981).

Schoeberl, M. R.↗

Simulation of ablation in Earth atmospheric entry

The process of ablation for Earth atmospheric entry is simulated using a computational approach that allows thermo-chemical nonequilibrium of the flow field and ablation gases. The heat pulse into the heat shield is modeled. The flowfield and graphite heat shield are coupled through surface mass and energy balances. The surface thermochemistry involves the oxidation of graphite and allows for catalytic recombination of diatomic oxygen. Steady-state simulations are performed on a one meter nose radius sphere at an altitude of 65/km and at freestream velocities of 8 km/s and 10 km/s. A transient simulation is performed at 65 km altitude and a freestream velocity of 10 km/s.

Keenan, James A.↗

A SIMULINK environment for flight dynamics and control analysis: Application to the DHC-2 Beaver. Part 1: Implementation of a model library in SIMULINK. Part 2: Nonlinear analysis of the Beaver autopilot

The design of advanced Automatic Aircraft Control Systems (AACS's) can be improved upon considerably if the designer can access all models and tools required for control system design and analysis through a graphical user-interface, from within one software environment. This MSc-thesis presents the first step in the development of such an environment, which is currently being done at the Section for Stability and Control of Delft University of Technology, Faculty of Aerospace Engineering. The environment is implemented within the commercially available software package MATLAB/SIMULINK. The report consists of two parts. Part I gives a detailed description of the AACS design environment. The heart of this environment is formed by the SIMULINK implementation of a nonlinear aircraft model in block-diagram format. The model has been worked out for the old laboratory aircraft of the Faculty, the De Havilland DHC-2 'Beaver', but due to its modular structure, it can easily be adapted for other aircraft. Part I also describes MATLAB programs which can be applied for finding steady-state trimmed-flight conditions and for linearization of the aircraft model, and it shows how the built-in simulation routines of SIMULINK have been used for open-loop analysis of the aircraft dynamics. Apart from the implementation of the models and tools, a thorough treatment of the theoretical backgrounds is presented. Part II of this report presents a part of an autopilot design process for the 'Beaver' aircraft, which clearly demonstrates the power and flexibility of the AACS design environment from part I. Evaluations of all longitudinal and lateral control laws by means of nonlinear simulations are treated in detail. The AACS design environment from part I proved to be a very useful tool for designing the control laws of the 'Beaver' autopilot within a very tight time-schedule. The autopilot design process itself will be used as a guideline for future AACS research at the Faculty of Aerospace Engineering. Flight tests of the 'Beaver' autopilot, done after evaluating the control laws in the SIMULINK package, proved to be quite successful. In the future, the AACS design package will evolve into a standardized, integrated design environment which can be applied to virtually any type of aircraft. The AACS design cycle will be shortened further by developing tools for automatically porting control laws from the MATLAB/SIMULINK environment to a piloted real-time flight simulator and the Flight Control Computers of the aircraft.

Flight Control System Design↗

Modeling of Biological Fluid Systems

This study presents a detailed methodology for modeling and validating three- dimensional fluid flow in realistic reconstructed vasculature geometry. Optical scanned data from a sectioned stenotic lumen having -65% area constriction within an arterial vessel of 9 mm diameter is processed and converted into a finite-element mesh. Mass and momentum transfer equations within the reconstructed stenotic vessel with a-flow rate of 200 ml/rnin is solved using the Galarkin Finite Element formulation. The calculated velocities at several axial locations along the length of the vessel are compared with the Particle Image Velocimetry (PIV) measurements. The computation results are in agreement with the initial steady state PIV measurements. In addition, initial results from pulsatile flow computations at different time points of a pulse cycle will be presented.

Banerjee, Rupak↗

Community Noise Impact of Urban Air Mobility Vehicle Operations

Advanced air mobility (AAM) missions, carried out by electrically driven air vehicles, are characterized by ranges of less than about 300-500 nm (about 500-900 km) and include both rural and urban operations. The missions may include public transportation, cargo delivery, air taxi, and emergency response. While the urban air mobility (UAM) subset of AAM is projected to have high economic benefit, it is also the most difficult to develop because it must overcome many barriers, including those associated with the airspace system, safety, and community noise. This presentation focuses on the utilization of UAM source noise data for assessment of community noise impact. Although a limited number of acoustic flight measurement campaigns have been made to characterize the source noise of prototype and preproduction UAM aircraft, prediction-based approaches are primarily considered herein. Following conceptual design, in which the vehicle is appropriately sized for its intended mission, a comprehensive analysis must be performed for a range of operating conditions spanning the flight envelope to determine the corresponding configurations of the vehicle, that is, the trimmed states. For each trimmed state, the noise produced by each source, for example, steady and unsteady rotor noise, may be computed and so-called source noise (hemi)spheres generated. These source data may subsequently be used in various community noise impact analyses. Several use cases are presented including those supporting noise certification and those for auralizations that may, in turn, be used as part of a perception-influenced design process. Land use planning tools for generating noise exposure maps, including those using simulation and integrated modeling approaches, are also presented. Finally, use cases supporting development of low noise flight operations, including an acoustic flight simulator and acoustically aware flight control, are considered.

aircraft community noise↗

A Theory of Control for a Class of Electronic Power Processing Systems: Energy-Storage DC-To-DC Converters

An analytically derived approach to the control of energy-storage dc-to-dc converters, which enables improved system performance and an extensive understanding of the manner in which this improved performance is accomplished, is presented. The control approach is derived from a state-plane analysis of dc-to-dc converter power stages which enables a graphical visualization of the movement of the system state during both steady state and transient operation. This graphical representation of the behavior of dc-to-dc converter systems yields considerable qualitative insight into the cause and effect relationships which exist between various commonly used converter control functions and the system performance which results from them.

Burns, W. W., III↗

Laser processes in plasmadynamic flows

Radiation emitted from the exhaust flow of a multi-megawatt, quasi-steady magnetoplasmadynamic accelerator has been examined photoelectrically for argon mass flows from 6 to 32 g/sec and arc currents up to 25 kA. Vacuum ultraviolet radiation, indicative of resonance transitions, has been detected photoelectrically at various axial locations downstream of the arc chamber, consistent with a relaxation length calculated using a collisional-radiative recombination model. Photoelectric measurements of particular transitions in the neutral and first ionic state of argon show that populations of excited electronic states of AII are in excess of that which would be expected on the basis of Saha-Boltzmann equilibrium at the local electron temperature and number density. Optical depth measurements indicate that for a current of 4 kA and argon mass flow of 12 g/sec, a population inversion exists between the upper and lower states of the 0.4880 micron argon ion transition.

Campbell, E. M.↗

H 2 O-assisted O 2 reduction by H 2 on Pt and PtAu bimetallic nanoparticles: Influences of composition and reactant coverages on kinetic regimes, rates, and selectivities

Hydrogen peroxide (H 2 O 2 ) can replace hazardous oxidants in industrial processes but is currently too expensive for many such applications. While direct synthesis of H 2 O 2 (H 2 + O 2 → H 2 O 2 ) may reduce costs in comparison to incumbent technology, current catalysts lack the requisite stability and selectivity. Here, we examine the direct synthesis of H 2 O 2 on bimetallic Pt 1 Au x (0 ≤ x ≤ 230) and Pt catalysts at steady-state in pure water and relate kinetic parameters for H 2 O 2 and H 2 O formation to possible active site structures informed by complementary characterization methods. X-ray photoelectron spectra show significant Pt surface enrichment compared to the bulk composition. Analysis of infrared spectra of mixed monolayers of 12 CO* and 13 CO* indicate that Pt and Au form substitutional surface alloys. The Pt 1 Au x nanoparticles with the greatest mole fractions of Au predominantly expose Pt monomers (i.e., isolated Pt atoms), yet Pt atoms exposed upon all these nanoparticles possess electronic structures distinct from bulk Pt. Despite these differences, rate measurements are consistent with product formation through proton-electron transfer pathways for all Pt 1 Au x catalysts. In situ XAS indicate that Pt remains metallic during H 2 O 2 synthesis. Under the most oxidizing conditions, selectivities toward H 2 O 2 increase strongly with the Au to Pt ratio from 2% for monometallic Pt to 85% for Pt 1 Au 170 . However, selectivities are similar among all catalysts within reducing conditions. Comparisons of apparent activation enthalpies for the formation of H 2 O 2 and H 2 O across these catalysts and the range of conditions suggest that Pt monomers within Au provide the greatest selectivities for H 2 O 2 formation, because these active sites present high barriers for O-O bond rupture. Further, selectivities decrease with increasing ratios of H 2 to O 2 pressures, because Pt atoms aggregate and form oligomers that readily dissociate dioxygen intermediates. The combined use of spectroscopy, kinetics, and concepts employed in reaching these conclusions take inspiration from the legacy of Prof. Michel Boudart, and specifically his elegant methods for interrogating bimetallic catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unlocking bulk and surface oxygen transport properties of mixed oxide-ion and electron conducting membranes with combined oxygen permeation cell and oxygen probe method

Surface exchange kinetics and bulk diffusion of oxygen are of paramount importance to the activity of oxygen electrocatalysis and performance of electrochemical devices such as fuel cell, metal-air batteries, and oxygen separation membranes. Conventional approaches to obtaining these transport properties are often limited to single property under a specific non-operation related condition. Here we use a combined oxygen permeation cell and oxygen probe methodology to simultaneously attain rates of oxygen surface exchange and bulk conductivity/chemical diffusivity of three representative mixed oxide-ion and electron conductors, namely SrCo 0.9 Ta 0.1 O 3-δ (SCT), La 0.6 Sr 0.4 CoO 3-δ (LSC) and La 0.6 Sr 0.4 FeO 3-δ (LSF), operated under a steady-state oxygen flux. The results explicitly show that SCT exhibit the highest oxide-ion conductivity/chemical diffusivity, fastest rates of surface oxygen exchange kinetics, thus promising to be the best oxygen electrocatalyst. Here, we have also mapped out the distribution of oxygen chemical potential gradient across the membranes and applied B-transport number concept to illustrate the rate-limiting steps in the overall oxygen permeation process.

42 ENGINEERING↗

Monte Carlo shock-like solutions to the Boltzmann equation with collective scattering

The results of Monte Carlo simulations of steady state shocks generated by a collision operator that isotropizes the particles by means of elastic scattering in some locally defined frame of reference are presented. The simulations include both the back reaction of accelerated particles on the inflowing plasma and the free escape of high-energy particles from finite shocks. Energetic particles are found to be naturally extracted out of the background plasma by the shock process with an efficiency in good quantitative agreement with an earlier analytic approximation (Eichler, 1983 and 1984) and observations (Gosling et al., 1981) of the entire particle spectrum at a quasi-parallel interplanetary shock. The analytic approximation, which allows a self-consistent determination of the effective adiabatic index of the shocked gas, is used to calculate the overall acceleration efficiency and particle spectrum for cases where ultrarelativistic energies are obtained. It is found that shocks of the strength necessary to produce galactic cosmic rays put approximately 15 percent of the shock energy into relativistic particles.

Ellison, D. C.↗

Solutal convection during directional solidification

During directional solidification of a binary alloy at constant velocity, buoyancy-driven fluid flow may occur due to the solute gradients generated by the solidification process. Numerical calculations of the solute and fluid flow fields in the melt have been carried out using finite differences in a two-dimensional, time-dependent model that assumes a planar crystal-melt interface and allows time-dependent gravitational accelerations. The container walls are rigid and perfectly insulating for solute. For constant vertical gravitational accelerations, as the solutal Rayleigh number is varied, multiple steady states and time-dependent states may occur. The bifurcation from the quiescent state may be subcritical or transcritical, depending on the aspect ratio of the container. Calculations have also been performed for a gravitational acceleration that is assumed to be uniform in magnitude with its direction rotating uniformly. Numerical results have been obtained for a Schmidt number of 10 and a gravitational acceleration of 0.0001 G. The maximum variation in the solute concentration at the crystal-melt interface is calculated for various values of the rotation rate of the gravitational acceleration.

Mcfadden, G. B.↗

The behavior of beams of relativistic non-thermal electrons under the influence of collisions and synchrotron losses

For many astrophysical situations, such as in solar flares or cosmic gamma-ray bursts, continuum gamma rays with energies up to hundreds of MeV were observed, and can be interpreted to be due to bremsstrahlung radiation by relativistic electrons. The region of acceleration for these particles is not necessarily the same as the region in which the radiation is produced, and the effects of the transport of the electrons must be included in the general problem. Hence it is necessary to solve the kinetic equation for relativistic electrons, including all the interactions and loss mechanisms relevant at such energies. The resulting kinetic equation for non-thermal electrons, including the effects of Coulomb collisions and losses due to synchrotron emission, was solved analytically in some simple limiting cases, and numerically for the general cases including constant and varying background plasma density and magnetic field. New approximate analytic solutions are presented for collision dominated cases, for small pitch angles and all energies, synchrotron dominated cases, both steady-state and time dependent, for all pitch angles and energies, and for cases when both synchrotron and collisional energy losses are important, but for relativistic electrons. These analytic solutions are compared to the full numerical results in the proper limits. These results will be useful for calculation of spectra and angular distribution of the radiation (x rays, gamma-rays, and microwaves) emitted via synchrotron or bremsstrahlung processes by the electrons. These properties and their relevance to observations will be observed in subsequent papers.

Mctiernan, James M.↗

The behavior of beams of relativistic nonthermal electrons under the influence of collisions and synchrotron losses

For many astrophysical situations, such as in solar flares or cosmic gamma-ray bursts, continuum gamma rays with energies up to hundreds of MeV were observed, and can be interpreted to be due to bremsstrahlung radiation by relativistic electrons. The region of acceleration for these particles is not necessarily the same as the region in which the radiation is produced, and the effects of the transport of the electrons must be included in the general problem. Hence it is necessary to solve the kinetic equation for relativistic electrons, including all the interactions and loss mechanisms relevant at such energies. The resulting kinetic equation for non-thermal electrons, including the effects of Coulomb collisions and losses due to synchrotron emission, was solved analytically in some simple limiting cases, and numerically for the general cases including constant and varying background plasma density and magnetic field. New approximate analytic solutions are presented for collision dominated cases, for small pitch angles and all energies, synchrotron dominated cases, both steady-state and time dependent, for all pitch angles and energies, and for cases when both synchrotron and collisional energy losses are important, but for relativistic electrons. These analytic solutions are compared to the full numerical results in the proper limits. These results will be useful for calculation of spectra and angular distribution of the radiation (X rays, gamma-rays, and microwaves) emitted via synchrotron or bremsstrahlung processes by the electrons. These properties and their relevance to observations will be observed in subsequent papers.

Mctiernan, James M.↗

Simplified Analysis of Pulse Detonation Rocket Engine Blowdown Gasdynamics and Performance

Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional designs, but represent a challenging modellng task. A simplified model for an idealized, straight-tube, single-shot PDRE blowdown process and thrust determination is described and implemented. In order to form an assessment of the accuracy of the model, the flowfield time history is compared to experimental data from Stanford University. Parametric Studies of the effect of mixture stoichiometry, initial fill temperature, and blowdown pressure ratio on the performance of a PDRE are performed using the model. PDRE performance is also compared with a conventional steady-state rocket engine over a range of pressure ratios using similar gasdynamic assumptions.

Morris, C. I.↗

Verification of a two-dimensional infiltration model for the resin transfer molding process

A two-dimensional finite element model for the infiltration of a dry textile preform by an injected resin was verified. The model, which is based on the finite element/control volume technique, determines the total infiltration time and the pressure increase at the mold inlet associated with the RTM process. Important input data for the model are the compaction and permeability behavior of the preform along with the kinetic and rheological behavior of the resin. The compaction behavior for several textile preforms was determined by experimental methods. A power law regression model was used to relate fiber volume fraction to the applied compaction pressure. Results showed a large increase in fiber volume fraction with the initial application of pressure. However, as the maximum fiber volume fraction was approached, the amount of compaction pressure required to decrease the porosity of the preform rapidly increased. Similarly, a power law regression model was used to relate permeability to the fiber volume fraction of the preform. Two methods were used to measure the permeability of the textile preform. The first, known as the steady state method, measures the permeability of a saturated preform under constant flow rate conditions. The second, denoted the advancing front method, determines the permeability of a dry preform to an infiltrating fluid. Water, corn oil, and an epoxy resin, Epon 815, were used to determine the effect of fluid type and viscosity on the steady state permeability behavior of the preform. Permeability values measured with the different fluids showed that fluid viscosity had no influence on the permeability behavior of 162 E-glass and TTI IM7/8HS preforms. Permeabilities measured from steady state and advancing front experiments for the warp direction of 162 E-glass fabric were similar. This behavior was noticed for tests conducted with corn oil and Epon 815. Comparable behavior was observed for the warp direction of the TTI IM7/8HS preform and corn oil. Mold filling and flow visualization experiments were performed to verify the analytical computer model. Frequency dependent electromagnetic sensors were used to monitor the resin flow front as a function of time. For the flow visualization tests, a video camera and high resolution tape recorder were used to record the experimental flow fronts. Comparisons between experimental and model predicted flow fronts agreed well for all tests. For the mold filling tests conducted at constant flow rate injection, the model was able to accurately predict the pressure increase at the mold inlet during the infiltration process. A kinetics model developed to predict the degree of cure as a function of time for the injected resin accurately calculated the increase in the degree of cure during the subsequent cure cycle.

Hammond, Vincent H.↗

Assessing the Consequences of Postclosure Criticality in Spent Nuclear Fuel

The U.S. Department of Energy is funding research into studying the consequences of postclosure criticality on the performance of a generic repository by (1) identifying the features, events, and processes (FEPs) that need to be considered in such an analysis, (2) developing the tools needed to model the relevant FEPs in a postclosure performance assessment, and (3) conducting analyses both with and without the occurrence of a postclosure criticality and comparing the results. Here, this paper describes progress in this area of research and presents the results to date of analyzing the consequences of a postulated steady-state criticality in a hypothetical saturated shale repository. Preliminary results indicate that postclosure criticality would not affect repository performance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power densities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. The purpose of this work is to verify the power distribution results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

97 MATHEMATICS AND COMPUTING↗

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power desnities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. The purpose of this work is to verify the power distribution results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗