Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “integral solution”

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 127 records · Page 7

Simplified combustion noise theory yielding a prediction of fluctuating pressure level

The first order equations for the conservation of mass and momentum in differential form are combined for an ideal gas to yield a single second order partial differential equation in one dimension and time. Small perturbation analysis is applied. A Fourier transformation is performed that results in a second order, constant coefficient, nonhomogeneous equation. The driving function is taken to be the source of combustion noise. A simplified model describing the energy addition via the combustion process gives the required source information for substitution in the driving function. This enables the particular integral solution of the nonhomogeneous equation to be found. This solution multiplied by the acoustic pressure efficiency predicts the acoustic pressure spectrum measured in turbine engine combustors. The prediction was compared with the overall sound pressure levels measured in a CF6-50 turbofan engine combustor and found to be in excellent agreement.

Huff, R. G.↗

Finite element-integral acoustic simulation of JT15D turbofan engine

An iterative finite element integral technique is used to predict the sound field radiated from the JT15D turbofan inlet. The sound field is divided into two regions: the sound field within and near the inlet which is computed using the finite element method and the radiation field beyond the inlet which is calculated using an integral solution technique. The velocity potential formulation of the acoustic wave equation was employed in the program. For some single mode JT15D data, the theory and experiment are in good agreement for the far field radiation pattern as well as suppressor attenuation. Also, the computer program is used to simulate flight effects that cannot be performed on a ground static test stand.

Baumeister, K. J.↗

An iwatsubo-based solution for labyrinth seals - comparison with experimental results

The basic equations are derived for compressible flow in a labyrinth seal. The flow is assumed to be completely turbulent in the circumferential direction where the friction factor is determined by the Blasius relation. Linearized zeroth and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equation. The first-order equations are solved by a separation of variables solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are compared to published test results.

Childs, D. W.↗

Quantitative non-destructive testing

The work undertaken during this period included two primary efforts. The first is a continuation of theoretical development from the previous year of models and data analyses for NDE using the Optical Thermal Infra-Red Measurement System (OPTITHIRMS) system, which involves heat injection with a laser and observation of the resulting thermal pattern with an infrared imaging system. The second is an investigation into the use of the thermoelastic effect as an effective tool for NDE. As in the past, the effort is aimed towards NDE techniques applicable to composite materials in structural applications. The theoretical development described produced several models of temperature patterns over several geometries and material types. Agreement between model data and temperature observations was obtained. A model study with one of these models investigated some fundamental difficulties with the proposed method (the primitive equation method) for obtaining diffusivity values in plates of thickness and supplied guidelines for avoiding these difficulties. A wide range of computing speeds was found among the various models, with a one-dimensional model based on Laplace's integral solution being both very fast and very accurate.

Welch, C. S.↗

A Comparison of Experimental and Theoretical Results for Labyrinth Gas Seals

The basic equations are derived for a two control volume model for compressible flow in a labyrinth seal. The flow is assumed to be completely turbulent and isoenergetic. The wall friction factors are determined using the Blasius formula. Jet flow theory is used for the calculation of the recirculation velocity in the cavity. Linearized zeroth and first order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth order pressure distribution is found by satisfying the leakage equation. The circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variable solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are compared to experimental test results.

Scharrer, Joseph Kirk↗

Case history of FAA/SRI wind shear models

In order to understand the development of the FAA/SRI wind fields, it is important to understand the operating philosophy of the FAA's Wind Shear Program Office. The goal of the office was to ensure an integrated solution to the wind shear problem which addressed three area: ground based equipment and coordination; airborne systems and procedures; and weather prediction. This triply addressed goal was central to the development of the wind fields. The primary user of the wind shear modeling during the FAA's program was airborne simulation. The project requirement was to use wind shear models that resulted from accidents so that effective procedures and/or equipment could be found for hazardous wind shear encounters. The wind shear model development is discussed in detail.

Schlickenmaier, Herbert↗

Dynamic delamination buckling in composite laminates under impact loading: Computational simulation

A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.

Grady, Joseph E.↗

Numerical optimization of conical flow waveriders including detailed viscous effects

A family of optimized hypersonic waveriders is generated and studied wherein detailed viscous effects are included within the optimization process itself. This is in contrast to previous optimized waverider work, wherein purely inviscid flow is used to obtain the waverider shapes. For the present waveriders, the undersurface is a streamsurface of an inviscid conical flowfield, the upper surface is a streamsurface of the inviscid flow over a tapered cylinder (calculated by the axisymmetric method of characteristics), and the viscous effects are treated by integral solutions of the boundary layer equations. Transition from laminar to turbulent flow is included within the viscous calculations. The optimization is carried out using a nonlinear simplex method. The resulting family of viscous hypersonic waveriders yields predicted high values of lift/drag, high enough to break the L/D barrier based on experience with other hypersonic configurations. Moreover, the numerical optimization process for the viscous waveriders results in distinctly different shapes compared to previous work with inviscid-designed waveriders. Also, the fine details of the viscous solution, such as how the shear stress is distributed over the surface, and the location of transition, are crucial to the details of the resulting waverider geometry. Finally, the moment coefficient variations and heat transfer distributions associated with the viscous optimized waveriders are studied.

Bowcutt, Kevin G.↗

Theory versus experiment for the rotordynamic coefficients of labyrinth gas seals. I - A two control volume model

The basic equations are derived for a two-control-volume model for compressible flow in a labyrinth seal. The recirculation velocity in the cavity is incorporated into the model for the first time. The flow is assumed to be completely turbulent and isoenergetic. The wall friction factors are determined using the Blasius formula. Jet flow theory is used for the calculation of the recirculation velocity in the cavity. Linearized zeroth- and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equations. The first-order equations are solved by a separation of variable solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients.

Scharrer, J. K.↗

Rotordynamic coefficients for stepped labyrinth gas seals

The basic equations are derived for compressible flow in a stepped labyrinth gas seal. The flow is assumed to be completely turbulent in the circumferential direction where the friction factor is determined by the Blasius relation. Linearized zeroth and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variables solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are presented in the form of a parametric study, since there are no known experimental data for the rotordynamic coefficients of stepped labyrinth gas seals. The parametric study investigates the relative rotordynamic stability of convergent, straight and divergent stepped labyrinth gas seals. The results show that, generally, the divergent seal is more stable, rotordynamically, than the straight or convergent seals. The results also show that the teeth-on-stator seals are not always more stable, rotordynamically, then the teeth-on-rotor seals as was shown by experiment by Childs and Scharrer (1986b) for a 15 tooth seal.

Scharrer, Joseph K.↗

Dynamics of coorbital satellite rings

The dynamical behavior of a coorbital satellite ring is studied for N = 2-9 satellites, in terms of a simplified dynamical description, in which the motion is reduced to the separation angles between satellites. The number and stability of different kinds of stationary configurations is explored, revealing that equally spaced rings are not stable against small perturbations for N not greater than 6, while for N = 2-8 there exists another, stable compact solution. Integrations of exact equations confirm these results. Moreover, the systems are found to display chaotic characteristics for a certain range of energy. The behavior can be interpreted in terms of maximum velocity curves, defining the allowed region of the motion in the phase space.

Salo, H.↗

Overview of a flywheel stack energy storage system

The concept of storing electrical energy in rotating flywheels provides an attractive substitute to batteries. To realize these advantages the critical technologies of rotor design, composite materials, magnetic suspension, and high efficiency motor/generators are reviewed in this paper. The magnetically suspended flywheel energy storage system, currently under development at the University of Maryland, consisting of a family of interference assembled rings, is presented as an integrated solution for energy storage.

Kirk, James A.↗

Dynamic delamination buckling in composite laminates under impact loading - Computational simulation

A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.

Grady, Joseph E.↗

Dynamics of spacecraft with deploying flexible appendages

A computational formulation for the dynamic analysis of spacecraft with deploying appendages is presented. The appendage model is based on a geometrically nonlinear beam formulation which accurately accounts for large rotational and large deformation motions. A moving finite element reference grid is incorporated within the nonlinear beam formulation to model the deployment motion. Hamilton's Law is used to formulate the general equations of motion, and a transient integration solution procedure is derived from a space-time finite element discretization of the Hamiltonian variational statement. Computational results of the methodology are presented for a classical gravity gradient stabilized satellite configuration.

Downer, Janice D.↗

A regularization method for extrapolation of solar potential magnetic fields

The mathematical basis of a Tikhonov regularization method for extrapolating the chromospheric-coronal magnetic field using photospheric vector magnetograms is discussed. The basic techniques show that the Cauchy initial value problem can be formulated for potential magnetic fields. The potential field analysis considers a set of linear, elliptic partial differential equations. It is found that, by introducing an appropriate smoothing of the initial data of the Cauchy potential problem, an approximate Fourier integral solution is found, and an upper bound to the error in the solution is derived. This specific regularization technique, which is a function of magnetograph measurement sensitivities, provides a method to extrapolate the potential magnetic field above an active region into the chromosphere and low corona.

Gary, G. A.↗

Fluctuations at the blue edge of saturated wind lines in IUE spectra of O-type stars

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-Lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles and then derived from formal solution integration using this source function. The more appropriate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10 percent or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Owocki, Stanley P.↗

Electrohydrodynamic deformation and interaction of a pair of emulsion drops

The response of a pair of emulsion drops to the imposition of a uniform electric field is examined. The case studied is that of equal-sized drops whose line of centers is parallel to the axis of the applied field. A new boundary integral solution to the governing equations of the leaky dielectric model is developed; the formulation accounts for the electrostatic and hydrodynamic interactions between the drops, as well as their deformations. Numerical calculations show that, after an initial transient during which the drops primarily deform, the pair drift slowly together due to their electrostatic interactions.

Baygents, James C.↗

On the synthesis of resonance lines in dynamical models of structured hot-star winds

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles are then derived from formal solution integration using this source function. Two more approximate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10% or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Puls, J.↗