Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient problem”

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 577 records · Page 32

Multiple-grid convergence acceleration of viscous and inviscid flow computations

A multiple-grid algorithm for use in efficiently obtaining steady solution to the Euler and Navier-Stokes equations is presented. The convergence of a simple, explicit fine-grid solution procedure is accelerated on a sequence of successively coarser grids by a coarse-grid information propagation method which rapidly eliminates transients from the computational domain. This use of multiple-gridding to increase the convergence rate results is substantially reduced work requirements for the numerical solution of a wide range of flow problems. Computational results are presented for subsonic and transonic inviscid flows and for laminar and turbulent, attached and separated, subsonic viscous flows. Work reduction factors as large as eight, in comparison to the basic fine-grid algorithm, were obtained. Possibilities for further performance improvement are discussed. Previously announced in STAR as N83-21847

Johnson, G. M.↗

Effect of a buried-wire gage on the separation bubble Numerical study

A numerically simulated buried-wire separation gage is investigated with emphasis on its effect on the separation bubble. The conjugated problem of a supersonic, time-dependent, two-dimensional flowfield above a conductive solid wall with an embedded heat source is solved using implicit finite difference algorithms. Steady-state and transient cases were computed for different locations of the heat source within the bubble. Results show that by using a steady heat source, the flow direction near the wall can be detected, without distorting the flowfield, only if the source is located in regions where the bubble is thick (i.e., not too close to the separation). The flow direction near separation can be detected by using a temperature pulse at the solid/fluid interface with insignificant distortion of the flowfield.

Degani, D.↗

X-ray imaging of flare loops and coronal arches

HXIS (Hard X-ray Imaging Spectrometer) results related to preflare, flare-decay, and post-flare phases in the solar corona are summarized. Particularly discussed are miniflares, long-lived X-ray enhancements along filament channels, flare precursors (some signifying the flare position, others possibly related to the onset of mass ejections), growth of post-flare loops, images of sites of field-line reconnections late in the flare development, giant post-flare coronal arches with energy contents comparable to coronal transients, and hot loops interconnecting active regions. The extremely low HXIS background, in particular, made it possible to discover many new, previously unknown features in the solar corona. Also emphasized are some problems that have appeared during the evaluation of HXIS data, for the benefit of those who plan future hard X-ray imaging experiments.

Svestka, Z.↗

Heat Transfer from Radiatively Heated Material in a Low Reynolds Number Microgravity Environment

A mathematical model of the transient three-dimensional heat transfer between a slowly moving ambient gas stream and a thermally thick or thin flat surface heated by external radiation in a microgravity environment is presented. The problem is motivated in part by fire safety issues in spacecraft. The gas phase is represented by variable property convection-diffusion energy and mass conservation equations valid at low Reynolds numbers. The absence of gravity and low Reynolds number together permit the flow to be represented by a self-consistent velocity potential determined by the ambient velocity and the thermal expansion in the gas. The solid exchanges energy with the gas by conduction/convection and with the surroundings by surface absorption and re-emission of radiation. Heat conduction in the solid is assumed to be one dimensional at each point on the surface as a consequence of the limited times (of order of 10 seconds) of interest in these simulations. Despite the apparent simplicity of the model, the results show a complex thermally induced flow near the heated surface. The thermal exchange between the gas and solid produces an outward sourcelike flow upstream of the center of the irradiated area and a sinklike flow downstream. The responses of the temperature fields and the associated flows to changes in the intensity of the external radiation and the ambient velocity are discussed.

Yamashita, H.↗

Separating Gravitational Wave Signals from Instrument Artifacts

Central to the gravitational wave detection problem is the challenge of separating features in the data produced by astrophysical sources from features produced by the detector. Matched filtering provides an optimal solution for Gaussian noise, but in practice, transient noise excursions or "glitches" complicate the analysis. Detector diagnostics and coincidence tests can be used to veto many glitches which may otherwise be misinterpreted as gravitational wave signals. The glitches that remain can lead to long tails in the matched filter search statistics and drive up the detection threshold. Here we describe a Bayesian approach that incorporates a more realistic model for the instrument noise allowing for fluctuating noise levels that vary independently across frequency bands, and deterministic "glitch fitting" using wavelets as "glitch templates", the number of which is determined by a trans-dimensional Markov chain Monte Carlo algorithm. We demonstrate the method's effectiveness on simulated data containing low amplitude gravitational wave signals from inspiraling binary black hole systems, and simulated non-stationary and non-Gaussian noise comprised of a Gaussian component with the standard LIGO/Virgo spectrum, and injected glitches of various amplitude, prevalence, and variety. Glitch fitting allows us to detect significantly weaker signals than standard techniques.

Littenberg, Tyson B.↗

ANALYSIS OF QUARTZ AND TEFLON SHIELDS FOR A PARTICULAR RE-ENTRY MISSION

The transient performance of ablation type heat protection shields is treated herein for the surface of a vehicle returning from outer space to the earth. The vehicle weighs 8640 kg, has a ballistic factor of 500 lb/ft2, re-enters with a speed of 11 km/sec at 120 km altitude, has a lift over drag ratio of -0.5, and is subjected to a maximum deceleration of 7.7 times the gravity constant. By use of well-known equations for the heat transfer and the mass transfer at a heated surface, a numerical calculation method is derived which, for the investigated ablation processes, yields exact transient solutions of the fundamental system of partial differential equations. The method is applied to various quartz shields and to one Teflon shield, which all evaporate so readily under the conditions of the problem at hand that practically no flow of molten shield material exists. The solutions also show comparatively small temperature changes parallel to the surface. The results show that the nose of the vehicle is cooled predominantly by the evaporation of the quartz or the Teflon; the rest of the vehicle's surface is cooled by radiation of the quartz or evaporation of the Teflon. The large mass transfer effects on the nose of the vehicle are detrimental since the resulting low surface temperatures prevent the radiative heat transfer out of the shield, which does not involve any mass loss, from being the desirable governing cooling factor.

ABLATING MATERIAL↗

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗

Observations on the burning of droplets in the absence of buoyancy

?The first stages of an experimental study on droplet burning at zero gravity are described. Motivations for the work are to investigate transient burning and effects of finite-rate chemistry, in a geometrical configuration sufficiently simple for theoretical analyses to be performed in detail and checked against experiment. Two possible approaches to the experimental problem are to establish burning in a freely falling chamber in a drop tower and to perform the experiments in a research facility in space, such as spacelab. The work reported here concerns the drop-tower experiments. The design of these experiments is discussed, and preliminary results of the experimental program are given. Critical problems in drop-tower testing of this kind, especially timing problems that arise in this particular experiment, are emphasized.

Knight, B.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

Experimental and Modeling Study of the Burning of an Ethanol Droplet in Microgravity

The microgravity ethanol droplet combustion experiments were performed aboard the STS-94/MSL-1 Shuttle mission within the Fiber-Supported Droplet Combustion-2 (FSDC-2) program. The burning histories and flame standoffs for pure ethanol and ethanol/water droplets were obtained from the images recorded with two 8 mm videocameras. The obtained results show that average gasification rate is related to the initial droplet size in a manner similar to n-alkanes and methanol and consistent with the results of Hara and Kumagai and the data taken recently in the NASA-Lewis 2.2 s droptower. A transient, moving finite-element chemically reacting flow model applied previously to sphero-symmetric combustion of methanol, methanol/water, n-alkane, and n-alkane binary mixture droplets was adopted for the problem of ethanol droplet combustion. The model includes detailed description of gas-phase reaction chemistry and transport, a simplified description of liquid phase transport, and non-luminous radiative heat transfer. Gas-phase chemistry was described with the detailed reaction mechanism of Norton and Dryer, which consists of 142 reversible elementary reactions of 33 species. Another recently published reaction mechanism of high-temperature ethanol oxidation was also considered. The model predictions were found to compare favorably with the experimental data. The model analysis also indicates that water condensation in the case of ethanol has smaller effect on average droplet gasification rate as compared with previously studied methanol cases. This effect is explained by non-ideal (azeotropic) behavior of binary ethanol-water mixtures. Further analysis of computational results and ethanol droplet radiative extinction behavior will be discussed.

Kazakov, Andrei↗

Flux-corrected transport techniques for transient calculations of strongly shocked flows

New flux-corrected transport algorithms are described for solving generalized continuity equations. These techniques were developed by requiring that the finite difference formulas used ensure positivity for an initially positive convected quantity. Thus FCT is particularly valuable for fluid-like problems with strong gradients or shocks. Repeated application of the same subroutine to mass, momentum, and energy conservation equations gives a simple solution of the coupled time-dependent equations of ideal compressible fluid dynamics without introducing an artificial viscosity. FCT algorithms span Eulerian, sliding-rezone, and Lagrangian finite difference grids in several coordinate systems. The latest FCT techniques are fully vectorized for parallel/pipeline processing.

Boris, J. P.↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Dynamics of high-bypass-engine thrust reversal using a variable-pitch fan

The test program demonstrated that successful and rapid forward-to reverse-thrust transients can be performed without any significant engine operational limitations for fan blade pitch changes through either feather pitch or flat pitch. For through-feather-pitch operation with a flight inlet, fan stall problems were encountered, and a fan blade overshoot technique was used to establish reverse thrust.

Schaefer, J. W.↗

Recent advances in thermostructural finite element analysis

A comparison study is described which has indicated that the finite element and lumped parameter methods are comparable for steady-state conduction/convection analyses. The finite element model was advantageous because of the ease of model verification with computer graphics. Some problems encountered in constructing compatible thermal and structural models have indicated the need for further methodology developments to automate the transition between thermal and structural analyses. The need for more efficient finite element nonlinear transient algorithms has been cited. With further development in the thermal analysis area, the finite element method offers high potential for an integrated thermal/structural analysis capability.

Thornton, E. A.↗

Recent advances in thermal-structural analysis and design

A continuing effort to develop a comprehensive capability for thermal-structural analysis and automated design (sizing) is reviewed. A principal role in the activity is played by the finite element program SPAR which contains both an efficient structural and thermal analysis capability. The benefit of having thermal and structural analyses in the same finite element program is illustrated by the application of SPAR to design calculations for the National Transonic Facility - a cryogenic wind tunnel under construction at Langley. Some experience with large-scale thermal structural analysis problems, particularly the space shuttle orbiter, has led to the identification of some analysis needs. Those needs include automated model generation and data output for lumped parameter thermal analysis, faster solution methods for nonlinear transient heat transfer, automated interpolation of temperature data from a thermal finite element model to a dissimilar structural finite model, and automated techniques to identify the times at which the critical combinations of transient heating and loads occur on a structure.

Adelman, H. M.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation is developed for sound propagation in a two-dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic-difference equations and boundary conditions are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

Sensitivity analysis of discrete systems

A selection is presented of methods used in the calculation of sensitivity derivatives of behavior and response of discretely modeled structure. Attention is given to methods for the selection of finite-difference step sizes for structural constraint derivatives, and techniques for finite difference derivatives in iteratively-solved problems. Also noted are methods for sensitivity of linear static response, and methods for second derivatives of structural response, derivatives of stresses, and nonlinear response; an emerging area of interest concerns the derivatives of transient response.

Adelman, Howard M.↗