Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state”

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 1,423 records · Page 79

Measurement of transient strain and surface temperature on simulated turbine blades using noncontacting techniques

Noncontacting techniques were used to measure strain and temperature in thermally cycled simulated turbine blades. An electro-optical extensometer was used to measure the displacement between parallel targets mounted on the leading edge of the blades throughout a complete heating and cooling cycle. An infrared photographic pyrometry method was used to measure blade steady state surface temperature. The blade was cyclically heated and cooled by moving it into and out of a Mach 1 hot-gas stream. Transient leading edge strain and steady state surface temperature distributions are presented for blades of three different configurations.

Calfo, F. D.↗

Modelling the growth of triglycine sulphate crystals in Spacelab 3

Two triglycine sulphate crystals were grown from an aqueous solution in Spacelab 3 aboard a Space Shuttle. Using a diffusion coefficient of 0.00002 sq cm/s, a computerized simulation gave reasonable agreement between experimental and theoretical crystal sizes and interferometric lines in the solution near the growing crystal. This diffusion coefficient is larger than most measured values, possibly due to fluctuating accelerations on the order of .001 g (Earth's gravity). The average acceleration was estimated to be less than .000001 g. At this level, buoyancy driven convection is predicted to add approx. 20 percent to the steady state growth rate. Only very slight distortion of the interferometric lines was observed at the end of a 33 hr run. It is suggested that the time to reach steady state convective transport may be inversely proportional to g at low g, so that the full effect of convection was not realized in these experiments.

Yoo, Hak-Do↗

Nonlinear Asymptotic Integration Algorithms for One-dimensional Autonomous Dissipative First-order Odes

Nonlinear asymptotic integrators are applied to one-dimensional, nonlinear, autonomous, dissipative, ordinary differential equations. These integrators, including a one-step explicit, a one-step implicit, and a one- and two-step midpoint algorithm, are designed to follow the asymptotic behavior of a system approaching a steady state. The methods require that the differential equation be written in a particular asymptotic form. This is always possible for a one-dimensional equation with a globally asymptotic steady state. In this case, conditions are obtained to guarantee that the implicit algorithms are well defined. Further conditions are determined for the implicit methods to be contractive. These methods are all first order accurate, while under certain conditions the midpoint algorithms may also become second order accurate. The stability of each method is investigated and an estimate of the local error is provided.

Haslach, Henry W., Jr.↗

Disturbance torque rejection properties of the NASA/JPL 70-meter antenna axis servos

Analytic methods for evaluating pointing errors caused by external disturbance torques are developed and applied to determine the effects of representative values of wind and friction torque. The expressions relating pointing errors to disturbance torques are shown to be strongly dependent upon the state estimator parameters, as well as upon the state feedback gain and the flow versus pressure characteristics of the hydraulic system. Under certain conditions, when control is derived from an uncorrected estimate of integral position error, the desired type 2 servo properties are not realized and finite steady-state position errors result. Methods for reducing these errors to negligible proportions through the proper selection of control gain and estimator correction parameters are demonstrated. The steady-state error produced by a disturbance torque is found to be directly proportional to the hydraulic internal leakage. This property can be exploited to provide a convenient method of determining system leakage from field measurements of estimator error, axis rate, and hydraulic differential pressure.

Hill, R. E.↗

Dynamic process intensification: Fundamentals and implementation to ternary distillation

Distillation remains a key technology for separating liquid mixtures. Its versatility comes with the disadvantage of high energy consumption. We previously used empirical arguments to introduce dynamic process intensification (DPI) as a strategy for improving the energy efficiency of binary distillation. In this article, we focus on ternary distillation; we begin by providing a rigorous basis for DPI, then formulate the problem of identifying the operating states for DPI as a nonlinear optimization problem. Via an extensive case study considering a hydrocarbon mixture, we demonstrate that reboiler energy use can be reduced by more than 2.3% relative to an equivalent column operating at steady state. Here, we prove that this result is due to the fact that DPI aims to satisfy product flow and quality constraints on average in time, affording broader opportunities for optimizing column economics than steady-state operation, where constraints are met strictly at all times.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Steady and Unsteady Nozzle Simulations Using the Conservation Element and Solution Element Method

This paper presents results from computational fluid dynamic (CFD) simulations of a three-stream plug nozzle. Time-accurate, Euler, quasi-1D and 2D-axisymmetric simulations were performed as part of an effort to provide a CFD-based approach to modeling nozzle dynamics. The CFD code used for the simulations is based on the space-time Conservation Element and Solution Element (CESE) method. Steady-state results were validated using the Wind-US code and a code utilizing the MacCormack method while the unsteady results were partially validated via an aeroacoustic benchmark problem. The CESE steady-state flow field solutions showed excellent agreement with solutions derived from the other methods and codes while preliminary unsteady results for the three-stream plug nozzle are also shown. Additionally, a study was performed to explore the sensitivity of gross thrust computations to the control surface definition. The results showed that most of the sensitivity while computing the gross thrust is attributed to the control surface stencil resolution and choice of stencil end points and not to the control surface definition itself.Finally, comparisons between the quasi-1D and 2D-axisymetric solutions were performed in order to gain insight on whether a quasi-1D solution can capture the steady and unsteady nozzle phenomena without the cost of a 2D-axisymmetric simulation. Initial results show that while the quasi-1D solutions are similar to the 2D-axisymmetric solutions, the inability of the quasi-1D simulations to predict two dimensional phenomena limits its accuracy.

Nozzle gross thrust calculations↗

Analytical description of charged particles transport along arbitrary guiding-field configurations

A new description of focused transport was developed which is valid for arbitrary spatial dependences of both the scattering mean free path and the focusing length. In particular, it describes the supercoherent transition from the coherent modes that occurs when focusing is strong compared to scattering to the diffusive mode that occurs when focusing is weak. Because it specifies the pitch angle distribution in terms of relatively simple functions, the new formulation facilitates not only the detailed interpretation of solar particle events, but also the qualitative understanding of focused transport. In the steady state, it predicts that the angular distribution consists of a highly collimated 'strahl' which is superimposed on an anisotropic 'halo'. This picture provides new insight into observed configurations of steady-state particle fluxes.

Earl, J. A.↗

Factors and processes governing the C-14 content of carbonate in desert soils

A model is presented describing the factors and processes which determine the measured C-14 ages of soil calcium carbonate. Pedogenic carbonate forms in isotopic equilium with soil CO2. Carbon dioxide in soils is a mixture of CO2 derived from two biological sources: respiration by living plant roots and respiration of microorganisms decomposing soil humus. The relative proportion of these two CO2 sources can greatly affect the initial C-14 content of pedogenic carbonate: the greater the contribution of humus-derived CO2, the greater the initial C-14 age of the carbonate mineral. For any given mixture of CO2 sources, the steady-state (14)CO2 distribution vs. soil depth can be described by a production/diffusion model. As a soil ages, the C-14 age of soil humus increases, as does the steady-state C-14 age of soil CO2 and the initial C-14 age of any pedogenic carbonate which forms. The mean C-14 age of a complete pedogenic carbonate coating or nodule will underestimate the true age of the soil carbonate. This discrepancy increases the older a soil becomes. Partial removal of outer (and younger) carbonate coatings greatly improves the relationship between measured C-14 age and true age. Although the production/diffusion model qualitatively explains the C-14 age of pedogenic carbonate vs. soil depth in many soils, other factors, such as climate change, may contribute to the observed trends, particularily in soils older than the Holocene.

Amundson, Ronald↗

In situ measurements of the NO2/NO ratio for testing atmospheric photochemical models

Simultaneous in situ measurements of NO2, NO, O3, ClO, pressure and temperature have been made for the first time, presenting a unique opportunity to test our current understanding of the photochemistry of the lower stratospere. Data were collected from several flights of the ER-2 aircraft at mid-latitudes in May 1993 during NASA's Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE). The daytime ratio of NO2/NO remains fairly constant at 19 km with a typical value of 0.68 and standard deviation of +/- 17. The ratio observations are compared with simple steady-state calculations based on laboratory-measured reaction rates and modeled NO2 photolysis rates. At each measurement point the daytime NO2/NO with its measurements uncertainty overlap the results of steady-state caculations and associated uncertainty. Possible sources of error are examined in both model and measurements. It is shown that more accurate laboratory determinations of the NO + 03 reaction rate and of the NO2 cross-sections in the 200-220 K temperature range characteristic of the lower stratosphere would allow for a more robust test of our knowledge of NO(X) phtochemistry by reducing significant sources if uncertainties in the interpretation of statospheric measurements.

Jaegle, L.↗

Structural Evaluation of a Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump Turbine Blade

Thermal and structural finite-element analyses were performed on the first high pressure fuel turbopump turbine blade of the space shuttle main engine (SSME). A two-dimensional (2-D) finite-element model of the blade and firtree disk attachment was analyzed using the general purpose MARC (finite-element) code. The loading history applied is a typical test stand engine cycle mission, which consists of a startup condition with two thermal spikes, a steady state and a shutdown transient. The blade material is a directionally solidified (DS) Mar-M 246 alloy, the blade rotor is forged with waspalloy material. Thermal responses under steady-state and transient conditions were calculated. The stresses and strains under the influence of mechanical and thermal loadings were also determined. The critical regions that exhibited high stresses and severe localized plastic deformation were the blade-rotor gaps.

Abdul-Aziz, Ali↗

In Situ Measurements of the NO2/NO Ratio for Testing Atmospheric Photochemical Models

Simultaneous in situ measurements of NO2, NO, O3, ClO, pressure and temperature have been made for the first time, presenting a unique opportunity to test our current understanding of the photochemistry of the lower stratosphere. Data were collected from several flights of the ER-2 aircraft at mid-latitudes in May 1993 during NASA's Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE). The daytime ratio of NO2/NO remains fairly constant at 19 km with a typical value of 0.68 and standard deviation of +/- 0.17. The ratio observations are compared with simple steady-state calculations based on laboratory-measured reaction rates and modeled NO2 photolysis rates. At each measurement point the daytime NO2/NO with its measurement uncertainty overlap the results of steady-state calculations and associated uncertainty. However, over all the ER-2 flights examined, the model systematically overestimates the ratio by 40% on average. Possible sources of error are examined in both model and measurements. It is shown that more accurate laboratory determinations of the NO + O3 reaction rate and of the NO2 cross-sections in the 200-220 K temperature range characteristic of the lower stratosphere would allow for a more robust test of our knowledge of NO(x) photochemistry by reducing significant sources of uncertainties in the interpretation of stratospheric measurements. The present measurements are compared with earlier observations of the ratio at higher altitudes.

Jaegle, L.↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the POLAR satellite is studied by using a 3-dimensional Particle-In-Cell (PIC) code. When the satellite passes through the region of low density plasma, the satellite charges to positive potentials as high as 4050Volts due to the photoelectrons emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, an ion-rich Xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at about 2Volts. Accordingly, in our 3-dimensional PIC simulation, we considered that the potential of the satellite is 2Volts as a fixed bias. Considering the relatively high density of the Xenon plasma in the sheath (approx. 10 - 10(exp 3)/cc), the ambient plasma of low density (less than 1/cc) is neglected. In the simulations, the electric fields and plasma dynamics are calculated self-consistently. We found that an "Apple" shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission, a high positive potential hill develops. Near the Thermal Ion Detection Experiment (TIDE) detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations, it takes about a hundred electron gyroperiods for the sheath to reach a quasi-steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. Using the steady state sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Leung, Wing C.↗

Thermal Conductivity of Ceramic Thermal Barrier and Environmental Barrier Coating Materials

Thermal barrier and environmental barrier coatings (TBC's and EBC's) have been developed to protect metallic and Si-based ceramic components in gas turbine engines from high temperature attack. Zirconia-yttria based oxides and (Ba,Sr)Al2Si2O8(BSAS)/mullite based silicates have been used as the coating materials. In this study, thermal conductivity values of zirconia-yttria- and BSAS/mullite-based coating materials were determined at high temperatures using a steady-state laser heat flux technique. During the laser conductivity test, the specimen surface was heated by delivering uniformly distributed heat flux from a high power laser. One-dimensional steady-state heating was achieved by using thin disk specimen configuration (25.4 mm diam and 2 to 4 mm thickness) and the appropriate backside air-cooling. The temperature gradient across the specimen thickness was carefully measured by two surface and backside pyrometers. The thermal conductivity values were thus determined as a function of temperature based on the 1-D heat transfer equation. The radiation heat loss and laser absorption corrections of the materials were considered in the conductivity measurements. The effects of specimen porosity and sintering on measured conductivity values were also evaluated.

Zhu, Dong-Ming↗

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Quantification of the importance of the modeling fidelity in the Best Estimate Plus Uncertainty analysis

In this work, the multi-physics coupling of CTF-PARCS for both steady-state and transient simulations is evaluated focusing on the quantification of the importance of the modeling fidelity when applied to the Best Estimate Plus Uncertainty safety analysis. For the steady-state, a Picard iteration approach is adopted in the coupled code while for the transient an explicit approach has been developed as a first step. The CTF-PARCS coupling is applied to the REA at HFP BOC scenario proposed in the exercise I-2c of LWR-UAM benchmark Phase III. Two different modeling choices are investigated for the nominal calculations. The first consists of including the 3D burnup information in CTF, and the second of modeling the spacers grids in both PARCS and CTF. It is observed that due to the degradation of the fuel thermal conductivity, the 3D burnup inclusion in CTF leads to significant higher centerline temperatures than the base case, with a difference up to 300 K. The spacer grids modeling shows the expected power and temperature decrease in the axial positions of the spacer grids. For the uncertainty quantification study, the multi-physics modeling including the burnup information in CTF is used. The results show that the inserted reactivity has a relative standard deviation of 13.6% mainly attributed to the nuclear data. The maximum linear power has a relative standard deviation of 16.2% due to the nuclear data as well, while the maximum centerline temperature has a relative standard deviation of 3.5% primarily caused by the fuel thermal conductivity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗

Cavity Swelling in Irradiated Materials

Cavity swelling, due to three-dimensional clustering of vacancies, is an undesirable isotropic volume expansion of materials under irradiation. This phenomenon occurs above stage-III recovery temperature where vacancies are mobile, between ~0.3–0.6 Tm where Tm is the absolute melting temperature. The primary mechanism for cavity swelling is “dislocation bias”, i.e., preferential absorption of interstitials by network dislocations, causing a vacancy supersaturation. Cavity swelling is highly temperature dependent and typically follows a bell-shaped curve with peak swelling occurring at an intermediate temperature. With increasing temperature, cavity size increases, and number density decreases logarithmically. With increasing dose, cavity swelling increases, but follows a hockey-stick type pattern, characterized by a low-swelling transient incubation period and then a steady-state regime. The incubation period, dominant in cavity nucleation, is dependent upon material’s chemistry and irradiation conditions such as dose rate/temperature. Cavity swelling increases linearly with dose in the steady-state regime, that is independent of chemistry and dose rate. This regime is dominated by cavity growth. Transmutation gases like He enhance cavity nucleation such that an early onset of cavity formation occurs. However, swelling magnitudes can be either increased or suppressed when He is present. Results in austenitic steels and Cu indicate that cavity swelling peaks at intermediate He/dpa ratios of ~10 appm He/dpa. The most effective method to mitigate cavity swelling is by increasing the fixed point-defect sink density, by promoting precipitation or dispersion strengthening.

Bhattacharya, Arunodaya↗

Air quality and comfort constrained energy efficient operation of multi-zone buildings

Maintaining indoor air quality (IAQ) through effective ventilation is essential for the well-being and productivity of building occupants. Control strategies aimed at improving the efficiency of heating, ventilation and air conditioning (HVAC) systems must jointly determine ventilation and heating and cooling processes. Here, in this paper, we study the problem of minimizing the energy consumption of the HVAC system in a multi-zone building, while meeting thermal comfort and IAQ requirements. We first perform a steady state analysis of the zonal carbon dioxide (CO 2 ) concentration and the temperature dynamics. The resulting expressions are convex in the zonal mass flow rates and zonal temperatures. Guided by the steady state solutions for meeting the thermal comfort constraints, we develop two control policies for improving the energy efficiency of building HVAC systems while jointly satisfying indoor temperature and IAQ constraints. We compare the performance of our proposed approaches with those of multiple baseline approaches which implement separate regimes for controlling zonal temperature and IAQ for a typical work-day in a multi-zone campus building. We have evaluated the performance of our proposed approaches under varying levels of flexibility in zonal temperatures. We have shown that zonal temperature flexibility can result in energy savings up to 32% (for the same control strategies) as compared to the case where no such flexibility is permitted. Our proposed approaches were seen to offer potential savings of nearly 29% compared to the baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗