Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 559 records · Page 31

Turbulent flame propagation and combustion in spark ignition engines

Pressure measurements synchronized with high-speed motion-picture records of flame propagation have been made in a transparent-piston engine. The data show that the initial expansion speed of the flame front is close to that of a laminar flame. As the flame expands, its speed rapidly accelerates to a quasi-steady value comparable with that of the turbulent velocity fluctuations in the unburned gas. During the quasi-steady propagation phase, a significant fraction of the gas behind the visible front is unburned. Final burnout of the charge may be approximated by an exponential decay in time. The data have been analyzed in a model-independent way to obtain a set of empirical equations for calculating mass burning rates in spark-ignition engines. The burning equations contain three parameters: the laminar burning speed, a characteristic speed (uT), and a characteristic length (lT). The laminar burning speed is known from laboratory measurements. Tentative correlations relating uT and lT to engine geometry and operating variables have been derived from the engine data.

Beretta, G. P.↗

Static and dynamic characteristics of parallel-grooved seals

Presented is an analytical method to determine static and dynamic characteristics of annular parallel-grooved seals. The governing equations were derived by using the turbulent lubrication theory based on the law of fluid friction. Linear zero- and first-order perturbation equations of the governing equations were developed, and these equations were analytically investigated to obtain the reaction force of the seals. An analysis is presented that calculates the leakage flow rate, the torque loss, and the rotordynamic coefficients for parallel-grooved seals. To demonstrate this analysis, we show the effect of changing number of stages, land and groove width, and inlet swirl on stability of the boiler feed water pump seals. Generally, as the number of stages increased or the grooves became wider, the leakage flow rate and rotor-dynamic coefficients decreased and the torque loss increased.

Iwatsubo, Takuzo↗

The use of a numerical filter to correct airborne temperature measurements for the effects of sensor lag

A numerical filter for transforming measured temperature signals into a close approximation of the actual temperature signal is described. The filter is derived by minimizing the mean-square error of the system, and assuming a knowledge of the characteristics of the sensing element and its housing. The equation representing the frequency-response function of the numerical filter is given. Input and output spectra for a filter applied to a case with negligible noise and a noise level of 1.5 percent of the total power in the input spectrum are analyzed, and the numerical weights for these two cases are calculated. Phase angle and gain for the entire system are examined. It is noted that the filter can enhance spectral components as high as 8 Hz with little phase and gain degradation over the bandwidth.

Ritter, John A.↗

Life prediction of thermomechanical fatigue using total strain version of strainrange partitioning (SRP): A proposal

A method is proposed (without experimental verification) for extending the total strain version of Strainrange Partitioning (TS-SRP) to predict the lives of thermomechanical fatigue (TMF) cycles. The principal feature of TS SRP is the determination of the time-temperature-waveshape dependent elastic strainrange versus life lines that are added subsequently to the classical inelastic strainrange versus life lines to form the total strainrange versus life relations. The procedure is based on a derived relation between failure and flow behavior. Failure behavior is represented by conventional SRP inelastic strainrange versus cyclic life relations, while flow behavior is captured in terms of the cyclic stress-strain response characteristics. Stress-strain response is calculated from simple equations developed from approximations to more complex cyclic constitutive models. For applications to TMF life prediction, a new testing technique, bithermal cycling, is proposed as a means for generating the inelastic strainrange versus life relations. Flow relations for use in predicting TMF lives would normally be obtained from approximations to complex thermomechanical constitutive models. Bithermal flow testing is also proposed as an alternative to thermomechanical flow testing at low strainranges where the hysteresis loop is difficult to analyze.

Saltsman, James F.↗

Development of low Reynolds number two equation turbulence models for predicting external heat transfer on turbine blades

A research effort was underway to study the use of two equation low Reynolds number turbulence models in predicting gas side heat transfer on turbine blades. The major objectives of this work are basicly threefold: study the predictive capabilities of two equation low Reynolds number turbulence models under the conditions characteristic of modern gas turbine blades; explore potential improvements to the models themselves as well as to the specification of initial conditions; and provide a comparison of the predictions of these models with the experimental data from a broad range of recently available turbine cascade experiments. The problems associated with predicting the boundary layer transition from laminar to turbulent flow are emphasized, as this may be the most serious deficiency of current modeling techniques. The results and conclusions of the first two phases are briefly described.

Patankar, Suhas V.↗

CAP-TSD: A program for unsteady transonic analysis of realistic aircraft configurations

The development of a new transonic code to predict unsteady flows about realistic aircraft configurations are described. An approximate factorization algorithm for solution of the unsteady transonic small disturbance equation is first described. Because of the superior stability characteristics of the AF algorithm, a new transonic aeroelasticity code was developed which is described in some detail. The new code was very easy to modify to include the additional aircraft components, so in a very short period of time the code was developed to treat complete aircraft configurations. Finally, applications are presented which demonstrate many of the geometry capabilities of the new code.

Batina, John T.↗

Assessment of impingement effects in the isentropic core of a small satellite control thruster plume

Two computational techniques commonly employed in the calculation of rocket and thruster expansion plumes are assessed. These are the method of characteristics (MOC), which is derived from the continuum Euler equations, and the direct simulation Monte Carlo (DSMC) method, which adopts a discrete particle approach. These techniques vary both in the computational expense and in the accuracy and detail of the solutions that they provide, depending upon the regime of application. The assessment is made with reference to the plume expanding from a small monopropellant hydrazine thruster and concentrates on the isentropic core of the jet for the flow regime lying between the continuum and free molecular limits. It is found that the more numerically intensive DSMC method offers the better correspondence to the available experimental data. In addition, large differences in typical impingement effects such as drag force and heat transfer are found at the free molecular limit of the plume expansion for the two predictive techniques. It is concluded that accurate estimation of impingement potential may only be achieved through application of the discrete particle method.

Boyd, I. D.↗

Ozone measurements from the NOAA-9 and the Nimbus-7 satellites: Implications of short and long term variabilities

An overview is given of the measurements of total ozne and ozone profiles by the SBUV/2 instrument on the NOAA-9 spacecraft relative to similar measurements from the SBUV and TOMS instruments on Nimbus-7. It is shown that during the three year period from March 14, 1985, to February 28, 1988, when these data sets overlap, there have been significant changes in the calibrations of the three instruments which may be attributed to the drift of the NOSS-9 orbit to later equator crossing times (for SBUV/2). These changes in instrument characteristics have affected the absolute values of the trends derived from the three instruments, but their geophysical characteristics and response to short term variations are accurate and correlate well among the three instruments. For example, the total column ozone measured by the three instruments shows excellent agreement with respect to its day to day, seasonal, and latitudinal variabilities. At high latitudes, the day to day fluctuations in total ozone show a strong positive correlation with temperature in the lower stratosphere, as one might expect from the dynamical coupling of the two parameters at these latitudes.

Chandra, S.↗

Hypersonic waverider design from given shock waves

An attempt is made to generate waverider flows from given shock wave geometries using two approaches. In the first approach, axisymmetric flows are used to construct more general flows based on the concept of osculating cones. Conical waverider design can thus be extended to yield results also for shocks forming a slope surface. The second approach involves solving the ill-posed problem of prescribing a shock wave and finding the flow field behind it. A new numerical marching technique with some features of characteristic cross marching is used to solve the Euler equation. The code selection and some test cases are discussed.

Sobieczky, H.↗

Computational dispersion properties of horizontal staggered grids for atmospheric and ocean models

The computational dispersion properties of horizontally and time-horizontally staggered grids utilizing corresponding centered-difference techniques for approximation of the adjustment, or gravity wave equations, are examined in terms of their group velocity characteristics. Results are acquired for oceanic and atmospheric models, the former being characterized by a much smaller Rossby radius of deformation. For all grids considered additional filtering is required to control and even eliminate waves with poor computational dispersion characteristics. Computational dispersion properties along with other computational characteristics and requirements give some guidance for an optimal selection of an appropriate grid for an ocean or atmospheric model.

Fox-Rabinovitz, Michael S.↗

Analysis and design of transonic airfoils using streamwise coordinates

A new approach is developed for analysis and design of transonic airfoils. A set of full potential equivalent equations in von Mises coordinates is formulated from the Euler equations under the irrotationality and isentropic assumptions. This set is composed of a main equation for the main variable, y, and a secondary equations for the secondary variable, R. The main equation is solved by type dependent differencing combined with a shock point operator. The secondary equation is solved by marching from a non-characteristic boundary. Sample computations on NACA 0012 and biconvex airfoils show that, for the analysis problem, the present approach achieves good agreement with experimental C sub p distributions. For the design problem, the approach leads to a simple numerical algorithm in which the airfoil contour is calculated as part of the flow field solution.

Barron, R. M.↗

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗

Ionospheric plasma flow about a system of electrically biased flat plates

The steady state interaction of two electrically biased parallel plates immersed in a flowing plasma characteristic of low earth orbit is studied numerically. Fluid equations are developed to describe the motion of the cold positively charged plasma ions, and are solved using finite-differences in two dimensions on a Cartesian grid. The behavior of the plasma electrons is assumed to be described by the Maxwell-Boltzmann distribution. Results are compared to an analytical and a particle simulation technique for a simplified flow geometry consisting of a single semi-infinite negatively biased plate. Comparison of the extent of the electrical disturbance into the flowing plasma and the magnitude of the current collected by the plate is very good. The interaction of two equally biased parallel plates is studied as a function of applied potential. The separation distance at which the current collected by either plate decreases by five and twenty percent is determined as a function of applied potential. The percent decreases were based on a non-interacting case. The decrease in overall current is caused by a decrease in ionic density in the region between the plates. As the separation between the plates decreases, the plates collect the ions at a faster rate than they are supplied to the middle region by the oncoming plasma flow. The docking of spacecraft in orbit is simulated by moving two plates of unequal potential toward one another in a quasi-static manner. One plate is held at a large negative potential while the other floats electrically in the resulting potential field. It is found that the floating plate does not charge continuously negative as it approaches the other more negatively biased plate. Instead, it charges more and then less negative as ionic current decreases and then increases respectively upon approach. When the two plates come into contact, it is expected that the electrically floating plate will charge rapidly negative to a potential near that of the other plate.

Herr, Joel L.↗

Application of spring tabs to elevator controls

Equations are presented for calculating the stick-force characteristics obtained with a spring-tab type of elevator control. The main problems encountered in the design of a satisfactory elevator spring tab are to provide stick forces in the desired range, to maintain the force per g sufficiently constant throughout the speed range, to avoid undesirable "feel" of the control in ground handling or in flight at low airspeeds, and to prevent flutter. Examples are presented to show the design features of spring tabs required to solve these problems for airplanes of various sizes.

Phillips, William H↗

Boundary element analysis on vector and parallel computers

Boundary element analysis (BEA) can be characterized as a numerical technique that generally shifts the computational burden in the analysis toward numerical integration and the solution of nonsymmetric and either dense or blocked sparse systems of algebraic equations. Researchers have explored the concept that the fundamental characteristics of BEA can be exploited to generate effective implementations on vector and parallel computers. In this paper, the results of some of these investigations are discussed. The performance of overall algorithms for BEA on vector supercomputers, massively data parallel single instruction multiple data (SIMD), and relatively fine grained distributed memory multiple instruction multiple data (MIMD) computer systems is described. Some general trends and conclusions are discussed, along with indications of future developments that may prove fruitful in this regard.

Kane, J. H.↗

Theory of 2 omega(sub pe) radiation induced by the bow shock

A new radiation emission mechanism is proposed to explain electomagnetic radiation observed at twice the electron plasma frequency, 2 omega(sub pe), in the upstream region of the Earth's bow shock. This radiation had its origin at the electron foreshock boundary where energetic electron beams and intense narrow-band Langmiur waves are observed. The proposed emission mechanism results from the interaction of the electron beam and Langmuir waves that are backscattered off thermal ions. This interaction is described by a nonlinear dispersion equation which incorporates an effect owing to electron trajectory modulation by the backscattered Langmuir waves. Subsequent analysis of the dispersion equation reveals two important consequences. First, a long-wavelength electrostatic quasi-mode with frequency at 2 omega(sub pe) is excited, and second, the quasi-mode and the electomagnetic mode are nonlinearly coupled. The implication is that, when the excited 2 omega(sub pe) quasi-mode propagates in an inhomgeneous medium with slightly decreasing density, the quasi-mode can be converted directly into an electromagnetic mode. Hense the electomagnetic radiation at twice the plasma frequency is generated. Numerical solutions of the dispersion equation with the choice of parameters that describe physical characteristics of the electron foreshock are presented, which illustrates the viability of the new mechanism.

Yoon, Peter H.↗

A Three-Dimensional Linearized Unsteady Euler Analysis for Turbomachinery Blade Rows

A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic response characteristics of axial-flow turbomachinery blading. The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. In addition, a numerical model for linearized inviscid unsteady flow, which is based upon an existing nonlinear, implicit, wave-split, finite volume analysis, is described. These aerodynamic and numerical models have been implemented into an unsteady flow code, called LINFLUX. A preliminary version of the LINFLUX code is applied herein to selected, benchmark three-dimensional, subsonic, unsteady flows, to illustrate its current capabilities and to uncover existing problems and deficiencies. The numerical results indicate that good progress has been made toward developing a reliable and useful three-dimensional prediction capability. However, some problems, associated with the implementation of an unsteady displacement field and numerical errors near solid boundaries, still exist. Also, accurate far-field conditions must be incorporated into the FINFLUX analysis, so that this analysis can be applied to unsteady flows driven be external aerodynamic excitations.

Montgomery, Matthew D.↗