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 919 records · Page 51

Regenerative life support system research

Sections on modeling, experimental activities during the grant period, and topics under consideration for the future are contained. The sessions contain discussions of: four concurrent modeling approaches that were being integrated near the end of the period (knowledge-based modeling support infrastructure and data base management, object-oriented steady state simulations for three concepts, steady state mass-balance engineering tradeoff studies, and object-oriented time-step, quasidynamic simulations of generic concepts); interdisciplinary research activities, beginning with a discussion of RECON lab development and use, and followed with discussions of waste processing research, algae studies and subsystem modeling, low pressure growth testing of plants, subsystem modeling of plants, control of plant growth using lighting and CO2 supply as variables, search for and development of lunar soil simulants, preliminary design parameters for a lunar base life support system, and research considerations for food processing in space; and appendix materials, including a discussion of the CELSS Conference, detailed analytical equations for mass-balance modeling, plant modeling equations, and parametric data on existing life support systems for use in modeling.

Source record↗

Preconditioned methods for solving the incompressible and low speed compressible equations

The introduction of artificial time derivatives which permit faster convergence to the steady state is used to precondition steady state incompressible equation-solving acceleration methods. An optimal value is obtained for beta, given constant alpha, where alpha and beta are arbitrary functions in the general preconditioning stage. A symmetric hyperbolic system is formed by the resulting incompressible equations. Compressible equations in conservation form with slow flow are also considered.

Turkel, Eli↗

TOPEX Microwave Radiometer - Thermal design verification test and analytical model validation

The testing of the TOPEX Microwave Radiometer (TMR) is described in terms of hardware development based on the modeling and thermal vacuum testing conducted. The TMR and the vacuum-test facility are described, and the thermal verification test includes a hot steady-state segment, a cold steady-state segment, and a cold survival mode segment totalling 65 hours. A graphic description is given of the test history which is related temperature tracking, and two multinode TMR test-chamber models are compared to the test results. Large discrepancies between the test data and the model predictions are attributed to contact conductance, effective emittance from the multilayer insulation, and heat leaks related to deviations from the flight configuration. The TMR thermal testing/modeling effort is shown to provide technical corrections for the procedure outlined, and the need for validating predictive models is underscored.

Lin, Edward I.↗

Spectral solution of inviscid supersonic flows over wedges and axisymmetric cones

A shock-fitted multidomain spectral collocation method is used to solve both steady and unsteady inviscid supersonic flows over bodies. New aspects of the method include two subdomain interface types and a zonal solution procedure to get efficient convergence to steady-state in supersonic regions. Steady-state examples include flow over a sharp cone, a hyperbolic cone and a hyperbolic wedge. For the cone, the exact flow solution is used to show that the method is spectrally accurate. As an example of an unsteady flow, a calculation of the interaction of a free-stream hot ring with the flow over a sharp cone is presented.

Kopriva, David A.↗

A Heuristic Method for Determining the Necessary Time Duration of Electron Beam Tests ESD of Spacecraft Dielectrics

Electrostatic discharge or ESD can pose a significant risk to spacecraft in many space environments. Laboratory electron beam facilities can be used to test the performance of candidate spacecraft dielectrics. However, limited resources necessitate accelerated testing. The aim of this work is the development of a criterion for determining when an ESD test has run for sufficient time to capture representative ESD behavior. Such a criterion has the potential of saving hours of personnel and facility time per test. A comparison of the distributions of ESD event magnitudes from consecutive segments of an ESD test can be used to determine when the test has reached a quasi-steady state. Once this quasi-steady state has been observed the test may be truncated without a significant reduction in test fidelity.

Chinn, James↗

Transient Weighted Moving-Average Model of Photovoltaic Module Back-Surface Temperature

Accurate modeling of photovoltaic (PV) performance requires the precise calculation of module temperature. Currently, most temperature models rely on steady-state assumptions that do not account for the transient climatic conditions and thermal mass of the module. On the other hand, complex physics-based transient models are computationally expensive and difficult to parameterize. In order to address this, a new approach to transient thermal modeling was developed, in which the steady-state predictions from previous timesteps are weighted and averaged to accurately predict the module temperature at finer time scales. This model is informed by 3-D finite-element analyses, which are used to calculate the effect of wind speed and module unit mass on module temperature. The model, in application, serves as an added filter over existing steady-state models that smooths out erroneous values that are a result of intermittency in solar resource. Validation of this moving-average model has shown that it can improve the overall PV energy performance model accuracy by as much as 0.58% over steady-state models based on mean absolute error improvements and can significantly reduce the variability between the model predictions and measured temperature times series data.

14 SOLAR ENERGY↗

Frequency response in short thermocouple wires

Theoretical expressions are derived for the steady state frequency response of a thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for a nonuniform wire with unequal material properties and wire diameters across the junction. The amplitude ratio at low frequency omega approaches 0 agrees with the results of Scadron and Warshawsky (1952) for a steady state temperature distribution. Moreover, the frequency response for a nonuniform wire in the limit of infinite length l approaches infinity is shown to reduce to a simple expression that is analogous to the classic first order solution for a thermocouple wire with uniform properties. Theoretical expressions are also derived for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire and a two material wire with unequal material properties across the junction. For the case of a one material supported wire, an exact solution is derived which compares favorably with an approximate expression that only matches temperatures at the support junction. Moreover, for the case of a two material supported wire, an analytical expression is derived that closely correlates numerical results. Experimental measurements are made for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire (type K) and a two material wire (type T) with unequal material properties across the junction. The data for the amplitude ratio and phase angle are correlated to within 10 pct. with the theoretical predictions of Forney and Fralick (1991). This is accomplished by choosing a natural frequency omega sub n for the wire data to correlate the first order response at large gas temperature frequencies. It is found that a large bead size, however, will increase the amplitude ratio at low frequencies but decrease the natural frequency of the wire. The phase angle data are also distorted for imperfect junctions.

Forney, L. J.↗

An Experimental Evaluation of the Performance of Two Combination Pitot Pressure Probes

Experimental tests have been completed which recorded the ability of two combination steady state and high response time varying Pitot probe designs to accurately measure steady stagnation pressure at a single location in a flow field. Tests were conducted of double-barreled and coannular Prati probes in a 3.5 in. diameter free jet probe calibration facility from Mach 0.1 to 0.9. Geometric symmetry and pitch (-40 deg to 40 deg) and yaw (0 deg to 40 deg) angle actuation were used to fully evaluate the probes. These tests revealed that the double-barreled configuration induced error in its steady state measurement at zero incidence that increased consistently with jet Mach number to 1.1 percent at Mach 0.9. For all Mach numbers, the double-barreled probe nulled at a pitch angle of approximately 7.0 deg and provided inconsistent measurements when yawed. The double-barreled probe provided adequate measurements via both its steady state and high response tubes (within +/- 0.15 percent accuracy) over unacceptable ranges of biased pitch and inconsistent yaw angles which varied with Mach number. By comparison, the coannular probe provided accurate measurements (at zero incidence) for all jet Mach numbers as well as over a flow angularity range which varied from +/- 26.0 deg at Mach 0.3 deg to +/- 14.0 deg at Mach 0.9. Based on these results, the Prati probe is established as the preferred design. Further experimental tests are recommended to document the frequency response characteristics of the Prati probe.

Arend, David J.↗

Task scheduling in dataflow computer architectures

Dataflow computers provide a platform for the solution of a large class of computational problems, which includes digital signal processing and image processing. Many typical applications are represented by a set of tasks which can be repetitively executed in parallel as specified by an associated dataflow graph. Research in this area aims to model these architectures, develop scheduling procedures, and predict the transient and steady state performance. Researchers at NASA have created a model and developed associated software tools which are capable of analyzing a dataflow graph and predicting its runtime performance under various resource and timing constraints. These models and tools were extended and used in this work. Experiments using these tools revealed certain properties of such graphs that require further study. Specifically, the transient behavior at the beginning of the execution of a graph can have a significant effect on the steady state performance. Transformation and retiming of the application algorithm and its initial conditions can produce a different transient behavior and consequently different steady state performance. The effect of such transformations on the resource requirements or under resource constraints requires extensive study. Task scheduling to obtain maximum performance (based on user-defined criteria), or to satisfy a set of resource constraints, can also be significantly affected by a transformation of the application algorithm. Since task scheduling is performed by heuristic algorithms, further research is needed to determine if new scheduling heuristics can be developed that can exploit such transformations. This work has provided the initial development for further long-term research efforts. A simulation tool was completed to provide insight into the transient and steady state execution of a dataflow graph. A set of scheduling algorithms was completed which can operate in conjunction with the modeling and performance tools previously developed. Initial studies on the performance of these algorithms were done to examine the effects of application algorithm transformations as measured by such quantities as number of processors, time between outputs, time between input and output, communication time, and memory size.

Katsinis, Constantine↗

Analog Techniques for Measuring the Frequency Response of Linear Physical Systems Excited by Frequency-Sweep Inputs

Data-reduction methods using general-purpose analog computer equipment and compatible testing techniques for determining the frequency response of linear physical systems are examined. The techniques considered may be classed as steady state or transient depending on the method of excitation. The relative merits of periodic, slow sweep, and transient (rapid sweep) forcing functions are discussed and applications are given that relate to dynamic-response tests of aeroelastic systems. Two frequency-sweep-input methods are considered in detail. one case the sweep rate is sufficiently slow that the response is approximately the same as that for steady-state conditions. With this input the frequency response can be evaluated and displayed in real time while the test is in progress. Errors due to treating sweep data as steady state can be eliminated, when desired, by reanalyzing tape-recorded time histories of the input and output as transient rather than as periodic data. In the second method the frequency-response function is deter- mined from the system's transient response to a very rapid sweep input. The purpose of frequency sweep in this case is to provide sufficient harmonic content in the input to overcome noise while keeping the test time as short as possible. tests and limited flight-test data presented herein, it appears that a transient-type rapid-sweep forcing function offers a considerable saving in test time while preserving the accuracy possible with steady-state sinusoidal inputs.

Wilmer H. Reed III↗

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗

Development of OPTIMOD - The new Space Shuttle Orbiter/Payload Thermal Integration Model

The development of the Orbiter/Payload Thermal Integration Model (OPTIMOD), a simplified Space Shuttle Orbiter thermal math model to be used for payload thermal integration analysis, is described. The Orbiter lower fuselage math model has been reduced to the minimum level of detail consistent with accuracy requirements. Instead, an 'environment radiation sensor', a simple closed-end cylinder oriented in a fixed position with respect to the Orbiter, is introduced. Equations are presented which relate the cylinder surface heating rates to external fuselage steady-state temperatures. These equations were derived from a regression correlation analysis of steady-state Orbiter thermal solutions for 2372 points on orbit, representing a complete range of vehicle-earth-sun orientations. An exponential equation is developed to relate the steady-state data to time-variant temperatures for the particular transient attitude timeline being simulated. The accuracy of OPTIMOD is shown to be excellent in the payload bay and to compare favorably with previous models in Orbiter elements outside the bay.

Russell, David J.↗

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The flammability boundary of a solid is U-shaped with the ordinate as ambient oxygen percentage (or total pressure) and the abscissa represents flow velocity, gravity level, stretch rate, etc. The left side of the flammability boundary is the quenching branch, and the right side is the blow-off branch. Quenching is due to a weak flame whose heat release from combustion cannot overcome the excessive heat losses from, for example, radiation. Blow-off is due to a short gas residence time that is smaller than the chemical reaction time (traditionally referred to as the small Damköhler number effect). For thick solids, the heat conduction into the solid interior has a long transient after ignition before the solid temperature reaches the steady state. During this slow transient, there is excessive (i.e. more than that at the steady state) gas phase heat conduction into the solid interior. From the viewpoint of the quasi-static gaseous flame, this is an additional heat loss that will affect the gas flame's flammability. In this work, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is extended to account for the excessive heat loss into the solid beyond the steady burning state. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper to account for the transient solid effect. For each Ψ, the gas phase is considered reaching a quasi-steady state. Extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The U-shaped flammability boundary is now a series of boundaries with Ψ as a parameter. Flame structures under the same environment condition but with different Ψ are compared. For a given oxygen ambient, the extinction boundary Ψ vs stretch rate is an inverted U-shape. The boundary consists of a blowoff branch and a radiative quenching branch. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, species concentration, and reaction rate are presented along the extinction boundaries. Heat balance analysis is performed for both the gas phase and the gas-solid interface. Along the extinction boundary, the relative importance of conduction, convection, radiation, and chemical reaction are analyzed in detail.

Stagnation point diffusion flame↗

Development of BEM for ceramic composites

BEST-CMS (boundary element solution technology - composite modeling system) is an advanced engineering system for the micro-analysis of fiber composite structures. BEST-CMS is based upon the boundary element program BEST3D which was developed for NASA by Pratt and Whitney Aircraft and the State University of New York at Buffalo under contract NAS3-23697. BEST-CMS presently has the capabilities for elastostatic analysis, steady-state and transient heat transfer analysis, steady-state and transient concurrent thermoelastic analysis, and elastoplastic and creep analysis. The fibers are assumed to be perfectly bonded to the composite matrix, or in the case of static or steady-state analysis, the fibers may be assumed to have spring connections, thermal resistance, and/or frictional sliding between the fibers and the composite matrix. The primary objective of this user's manual is to provide an overview of all BEST-CMS capabilities, along with detailed descriptions of the input data requirements. In the next chapter, a brief review of the theoretical background is presented for each analysis category. Then, chapter three discusses the key aspects of the numerical implementation, while chapter four provides a tutorial for the beginning BEST-CMS user. The heart of the manual, however, is in chapter five, where a complete description of all data input items is provided. Within this chapter, the individual entries are grouped on a functional basis for a more coherent presentation. Chapter six includes sample problems and should be of considerable assistance to the novice. Chapter seven includes capsules of a number of fiber-composite analysis problems that have been solved using BEST-CMS. This chapter is primarily descriptive in nature and is intended merely to illustrate the level of analysis that is possible within the present BEST-CMS system. Chapter eight contains a detail description of the BEST-CMS Neutral File which is helpful in writing an interface between BEST-CMS and any graphic post-processor program. Finally, all pertinent references are listed in chapter nine.

Henry, D. P.↗

Further Developments of BEM for Micro and Macromechanical Analyses of Composites: Boundary Element Software Technology-Composite User's Manual

BEST-CMS (Boundary Element Solution Technology - Composite Modeling System) is an advanced engineering system for the micro-analysis of fiber composite structures. BEST-CMS is based upon the boundary element program BEST3D which was developed for NASA by Pratt and Whitney Aircraft and the State University of New York at Buffalo under contract NAS3-23697. BEST-CMS presently has the capabilities for elastostatic analysis, steady-state and transient heat transfer analysis, steady-state and transient concurrent thermoelastic analysis and elastoplastic and creep analysis. The fibers are assumed to be perfectly bonded to the composite matrix, or in the case of static or steady-state analysis, the fibers may be assumed to have spring connections, thermal resistance, and/or frictional sliding between the fibers and the composite matrix. The primary objective of this User's Manual is to provide an overview of all BEST-CMS capabilities, along with detailed descriptions of the input data requirements. A brief review of the theoretical background is presented for each analysis category. Then, Chapter 3 discusses the key aspects of the numerical implementation, while Chapter 4 provides a tutorial for the beginning BEST-CMS user. The heart of the manual, however, is in Chapter 5, where a complete description of all data input items is provided. Within this chapter, the individual entries are grouped on a functional basis for a more coherent presentation. Chapter 6 includes sample problems and should be of considerable assistance to the novice. Chapter 7 includes capsules of a number of fiber-composite analysis problems that have been solved using BEST-CMS. This chapter is primarily descriptive in nature and is intended merely to illustrate the level of analysis that is possible within the present BEST-CMS system. Chapter 8 contains a detailed description of the BEST-CMS Neutral File which is helpful in writing an interface between BEST- CMS and any graphic post-processor program. Finally, all pertinent references are listed in Chapter 9.

Banerjee, P. K.↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

CFD-Based Frequency Domain Method for Dynamic Stability Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a dynamic stability estimation technique obtained from high-fidelity CFD simulations of the Mach 0.8 Transonic Truss-Braced Wing (TTBW). A series of unsteady RANS CFD simulations in FUN3D is performed on the TTBW in pitch and plunge oscillations at various reduced frequencies. The time-domain data are transformed into the frequency-domain data by Fourier series. Transfer functions of the dynamic stability derivatives are then estimated by a frequency-domain regression. The dynamic stability derivatives with respect to the angle of attack are determined by the regression of the unsteady aerodynamic coefficients for the plunge motion. The dynamic stability derivatives with respect to the pitch rate are deter-mined by the regression of the differential unsteady aerodynamic coefficients for the pitch motion upon the removal of the angle of attack contribution by the plunge motion. The steady-state stability derivatives are then compared to the results obtained from a stability analysis code VSPAERO as well as steady-state FUN3D simulations. The comparison of the steady-state stability derivatives shows excellent agreement.

Dynamic Stability Derivatives↗