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 1,009 records · Page 56

Nonlinear effects in a plain journal bearing. I - Analytical study. II - Results

In the first part of this work, a numerical model is presented which couples the variable-property Reynolds equation with a rotor-dynamics model for the calculation of a plain journal bearing's nonlinear characteristics when working with a cryogenic fluid, LOX. The effects of load on the linear/nonlinear plain journal bearing characteristics are analyzed and presented in a parametric form. The second part of this work presents numerical results obtained for specific parametric-study input variables (lubricant inlet temperature, external load, angular rotational speed, and axial misalignment). Attention is given to the interrelations between pressure profiles and bearing linear and nonlinear characteristics.

Choy, F. K.↗

On the solution of a lubrication problem with particulate solids

The lubrication characteristic of a fluid with solid particles is studied using the continuum theory of mixtures. The governing equations are formulated and appropriate boundary conditions are introduced for an arbitrary-shaped lubricant film thickness. As a special case, closed-form analytical perturbation solutions for pressure and shear stress are obtained for a mixture of a conventional oil and solid particles with small values of solid-volume fraction sheared in the clearance space of an infinitely long slider bearing. It is found that when compared with a pure fluid, the mixture of the fluid and solid generates a higher pressure and therefore a higher load-carrying capacity with the added advantage of a reduction in the coefficient of friction.

Dai, F.↗

The Lateral Stability of Equal-flanged Aluminum-alloy I-beams Subjected to Pure Bending

Equal-flange beams of a special extruded I-section of 27ST aluminum alloy were tested in pure bending. Complete end fixity was not attained. Loading was continued until a definite maximum value had been reached. Tensile tests were made on specimens cut from the flanges and the web of each beam. Compressive stress-strain characteristics were determined by pack compression tests on specimens cut from the flanges. Values computed from an equation previously suggested by one of the authors for the critical stress at which such beams become unstable were found to be in good agreement with values computed from experimentally determined critically bending moments.

Dumont, C↗

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗

Inertial Mass Viewed as Reaction of the Vacuum to Accelerated Motion

Preliminary analysis of the momentum flux (or of the Poynting vector) of the classical electromagnetic version of the quantum vacuum consisting of zero-point radiation impinging on accelerated objects as viewed by an inertial observer suggests that the resistance to acceleration attributed to inertia may be a force of opposition originating in the vacuum. This analysis avoids the ad hoc modeling of particle-field interaction dynamics used previously by Haisck Rueda and Puthoff (1994) to derive a similar result. This present approach is not dependent upon what happens at the particle point but on how an external observer assesses the kinematical characteristics of the zero-point radiation impinging on the accelerated object. A relativistic form of the equation of motion results from the present analysis.

Rueda, Alfonso↗

Development of the Electron Drift Instrument (EDI) for Cluster

The Electron Drift Instrument (EDI) is a new technique for measuring electric fields in space by detecting the effect on weak beams of test electrons. This U.S. portions of the technique, flight hardware, and flight software were developed for the Cluster mission under this contract. Dr. Goetz Paschmann of the Max Planck Institute in Garching, Germany, was the Principle Investigator for Cluster EDI. Hardware for Cluster was developed in the U.S. at the University of New Hampshire, Lockheed Palo Alto Research Laboratory, and University of California, San Diego. The Cluster satellites carrying the original EDI instruments were lost in the catastrophic launch failure of first flight of the Arianne-V rocket in 1996. Following that loss, NASA and ESA approved a rebuild of the Cluster mission, for which all four satellites were successfully launched in the Summer of 2000. Limited operations of EDI were also obtained on the Equator-S satellite, which was launched in December, 1997. A satellite failure caused a loss of the Equator-S mission after only 5 months, but these operations were extremely valuable in learning about the characteristics and operations of the complex EDI instrument. The Cluster mission, satellites, and instruments underwent an extensive on-orbit commissioning phase in the Fall of 2000, carrying over through January 2001. During this period all elements of the instruments were checked and careful measurements of inter-experiments interferences were made. EDI is currently working exceptionally well in orbit. Initial results verify that all aspects of the instrument are working as planned, and returning highly valuable scientific information. The first two papers describing EDI on-orbit results have been submitted for publication in April, 2001. The principles of the EDI technique, and its implementation on Cluster are described in two papers by Paschmann et al., attached as Appendices A and B. The EDI presentation at the formal Cluster Commissioning Review, held at ESA Headquarters in Paris, is attached as Appendix C.

Quinn, Jack↗

A Transport Model for Non-Local Heating of Electrons in ICP Reactors

A new model has been developed for non-local heating of electrons in ICP reactors, based on a hydrodynamic approach. The model has been derived using the electron momentum conservation in azimuthal direction with electromagnetic and frictional forces respectively as driving force and damper of harmonic oscillatory motion of electrons. The resulting transport equations include the convection of azimuthal electron momentum in radial and axial directions, thereby accounting for the non-local effects. The azimuthal velocity of electrons and the resulting electrical current are coupled to the Maxwell's relations, thus forming a self-consistent model for non-local heating. This model is being implemented along with a set of Navier-Stokes equations for plasma dynamics and gas flow to simulate low-pressure (few mTorr's) ICP discharges. Characteristics of nitrogen plasma in a TCP 300mm etch reactor is being studied. The results will be compared against the available Langmuir probe measurements.

Chang, C. H.↗

The Acoustic Analogy and Alternative Theories for Jet Noise Prediction

This paper describes several methods for the prediction of jet noise. All but one of the noise prediction schemes are based on Lighthill's or Lilley's acoustic analogy while the other is the jet noise generation model recently proposed by Tam and Auriault. In all the approaches some assumptions must be made concerning the statistical properties of the turbulent sources. In each case the characteristic scales of the turbulence are obtained from a solution of the Reynolds-averaged Navier Stokes equation using a k - epsilon turbulence model. It is shown that, for the same level of empiricism, Tam and Auriault's model yields better agreement with experimental noise measurements than the acoustic analogy. It is then shown that this result is not because of some fundamental flaw in the acoustic analogy approach: but, is associated with the assumptions made in the approximation of the turbulent source statistics. If consistent assumptions are made, both the acoustic analogy and Tam and Auriault's model yield identical noise predictions. The paper concludes with a proposal for an acoustic analogy that provides a clearer identification of the equivalent source mechanisms and a discussion of noise prediction issues that remain to be resolved.

Morris, Philip J.↗

Right and left ventricular volumes in vitro by a new nongeometric method

We present an evaluation of a new nongeometric technique for calculating right and left ventricular volumes. This method calculates ventricular chamber volumes from multiple cross-sectional echocardiographic views taken from a single point as the echo beam is tilted progressively through the ventricle. Right and left ventricular volumes are calculated from both the approximate short axis and approximate apical position on 20 in vitro human hearts and compared with the actual chamber volumes. The results for both ventricles from both positions are excellent. Correlation coefficients are > 0.95 for all positions; the standard errors are in the range of 5 to 7 mL and the slopes and intercepts for the regression lines are not significantly different from 1 and 0, respectively (except for the left ventricular short-axis intercept). For all positions, approximately 6 to 8 views are needed for peak accuracy (7.5 degrees to 10 degrees separation). This approach offers several advantages. No geometric assumptions about ventricular shape are made. All images are acquired from a single point (or window), and the digitized points can be used to make a three-dimensional reconstruction of the ventricle. Also, during the calculations a volume distribution curve for the ventricle is produced. The shape of this curve can be characteristic for certain situations (ie, right ventricle, short axis) and can be used to make new simple equations for calculating volume. We conclude that this is an accurate nongeometric method for determining both right and left ventricular volumes in vitro.

NASA Discipline Cardiopulmonary↗

The Acoustic Analogy and Alternative Theories for Jet Noise Prediction

This paper describes several methods for the prediction of jet noise. All but one of the noise prediction schemes are based on Lighthill's or Lilley's acoustic analogy while the other is the jet noise generation model recently proposed by Tam and Auriault. In all the approaches some assumptions must be made concerning the statistical properties of the turbulent sources. In each case the characteristic scales of the turbulence are obtained from a solution of the Reynolds-averaged Navier Stokes equation using a k-epsilon turbulence model. It is shown that, for the same level of empiricism, Tam and Auriault's model yields better agreement with experimental noise measurements than the acoustic analogy. It is then shown that this result is not because of some fundamental flaw in the acoustic analogy approach: but, is associated with the assumptions made in the approximation of the turbulent source statistics. If consistent assumptions are made, both the acoustic analogy and Tam and Auriault's model yield identical noise predictions. The paper concludes with a proposal for an acoustic analogy that provides a clearer identification of the equivalent source mechanisms and a discussion of noise prediction issues that remain to be resolved.

Morris, Philip J.↗

Flutter Analysis for Turbomachinery Using Volterra Series

The objective of this paper is to describe an accurate and efficient reduced order modeling method for aeroelastic (AE) analysis and for determining the flutter boundary. Without losing accuracy, we develop a reduced order model based on the Volterra series to achieve significant savings in computational cost. The aerodynamic force is provided by a high-fidelity solution from the Reynolds-averaged Navier-Stokes (RANS) equations; the structural mode shapes are determined from the finite element analysis. The fluid-structure coupling is then modeled by the state-space formulation with the structural displacement as input and the aerodynamic force as output, which in turn acts as an external force to the aeroelastic displacement equation for providing the structural deformation. NASA's rotor 67 blade is used to study its aeroelastic characteristics under the designated operating condition. First, the CFD results are validated against measured data available for the steady state condition. Then, the accuracy of the developed reduced order model is compared with the full-order solutions. Finally the aeroelastic solutions of the blade are computed and a flutter boundary is identified, suggesting that the rotor, with the material property chosen for the study, is structurally stable at the operating condition, free of encountering flutter.

computational fluid dynamics↗

Implementation and Use of a Convergence Test for Iterative Wind Tunnel Balance Load Predictions

A modified and extended version of a convergence test for wind tunnel strain-gage balance load iterations was implemented. The test uses an upper bound of the Lipschitz constant to assess convergence characteristics of balance load predictions if the Iterative Method is applied. Convergence is expected within the use envelope of the balance whenever this upper bound is less than the threshold of one. It is explained in great detail how the convergence test can be applied to the two load iteration equation types that are currently being used in the aerospace testing community. In addition, the application of the test to balances with bi-directional output characteristics is discussed. It is also shown how numerical differentiation can be used in order to obtain partial derivatives that are needed for the calculation of the Lipschitz constant. Finally, machine calibration data of NASA's MC60E six-component force balance is selected to demonstrate both implementation and use of the convergence test.

calibration analysis↗

Aeroelastic Computations of a Compressor Stage Using the Harmonic Balance Method

The aeroelastic characteristics of a compressor stage were analyzed using a computational fluid dynamic (CFD) solver that uses the harmonic balance method to solve the governing equations. The three dimensional solver models the unsteady flow field due to blade vibration using the Reynolds-Averaged Navier-Stokes equations. The formulation enables the study of the effect of blade row interaction through the inclusion of coupling modes between blade rows. It also enables the study of nonlinear effects of high amplitude blade vibration by the inclusion of higher harmonics of the fundamental blade vibration frequency. In the present work, the solver is applied to study in detail the aeroelastic characteristics of a transonic compressor stage. Various parameters were included in the study: number of coupling modes, blade row axial spacing, and operating speeds. Only the first vibration mode is considered with amplitude of oscillation in the linear range. Both aeroelastic stability (flutter) of rotor blade and unsteady loading on the stator are calculated. The study showed that for the stage considered, the rotor aerodynamic damping is not influenced by the presence of the stator even when the axial spacing is reduced by nearly 25 percent. However, the study showed that blade row interaction effects become important for the unsteady loading on the stator when the axial spacing is reduced by the same amount.

Reddy, T. S. R.↗

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗

Identification of an object by input and output spectral characteristics

The problem is discussed of identification of a linear object of known structure, the movement of which is described by a system of differential equations of the type y = Ay + Bu, where y is an n-dimensional output vector, u is an m-dimensional vector of stationary, random disturbances (inputs), A and B are matrices of unknown parameters of the dimension, n x n and n x m, respectively. The spectral and reciprocal spectral densities of the inputs and outputs are used as the initial information on the object.

Redko, S. F.↗

Regression analysis of traction characteristics of traction fluids

Traction data for Santotrac 50 and TDF-88 over a wide range of operating conditions were analyzed. An eight term correlation equation to predict the maximum traction coefficient and a six term correlation equation to predict the initial slope of the traction curve were developed. The slope correlation was corrected for size effect considering the compliance of the disks. The effects of different operating conditions on the traction performance of each traction fluid were studied. Both fluids exhibited a loss in traction with increases in spin, but the losses with the TDF-88 fluid were not as severe as those with Santotrac 50. Overall, both fluids exhibited similar performance, showing an increase in traction with contact pressure up to about 2.0 GPa, and a reduction in traction with higher surface speeds up to about 100 m/sec. The apparent stiffness of the traction contact, that is, film disk combination, increases with contact pressure and decreases with speed.

Loewenthal, S. H.↗

Interactions among soil texture, pore structure, and labile carbon influence soil carbon gains

Perennial vegetation with high plant diversity, e.g., restored prairie, is known for stimulation of soil carbon (C) gains, due in part to enhanced formation of pore structure beneficial for long-term C storage. However, the prevalence of this phenomenon across soils of different types remains poorly understood. The aim of the study was to assess the associations between pore structure, soil C, and their differences in monoculture switchgrass and polyculture restored prairie vegetation across a wide range of soils dominating the Upper Midwest of the USA. Six experimental sites were sampled, representing three soil types with texture ranging from sandy to silt loams. The two vegetation systems studied at each site were (i) monoculture switchgrass (Panicum virgatum L.), and (ii) polyculture restored prairie, also containing switchgrass as one of its species. X-ray computed micro-tomography (µCT) was employed to analyze soil pore structure. Structural equation modeling and multiple path analyses were used to assess direct and indirect effects of soil texture and pore characteristics on microbial biomass C (MBC), particulate organic matter (POM), dissolved organic C (DOC), short-term respiration (CO 2 ), and, ultimately, soil organic C (SOC). Across studied sites, prairie increased fractions of medium (50–150 µm Ø) pores by 11–45 %, SOC by 3–69 %, and MBC by 18–59 % (except for one site). The greater were the prairie-induced increases in the medium pore volumes, the greater were the prairie-induced SOC gains. Greater C losses via CO 2 and DOC contributed to slower C accumulation in the prairie soil. We surmise that the interactive feedback loop relating medium pores and soil C acts across a wide range of soil textures and is an important mechanism through which perennial vegetation with high plant diversity, such as restored prairie, promotes rapid SOC gains.

54 ENVIRONMENTAL SCIENCES↗