SEARCH · Engineering Papers
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.
The problem of the finite displacement and buckling, of a shallow spherical dome is investigated both theoretically and experimentally. Experimental results seem to indicate that the classical criterion of buckling is applicable to very shallow spherical domes for which the theoretical calculation was made. A transition to energy criterion for higher domes is also indicated.
Leakage and wear are two fundamental problems in all traditional turbine seals that contribute to an engine's inefficiency. The solutions to seal leakage and wear conflict in the conventional design space. Reducing the clearance between the seal and rotating shaft reduces leakage but increases wear because of increased contact incidents. Increasing the clearance to reduce the contact between parts reduces wear but increases parasitic leakage. The goal of this effort is to develop a seal that restricts leakage flow using acoustic pressure while operating in a noncontacting manner, thereby increasing life. In 1996, Dr. Timothy Lucas announced his discovery of a method to produce shock-free high-amplitude pressure waves. For the first time, the formation of large acoustic pressures was possible using dissonant resonators. A pre-prototype acoustic seal developed at the NASA Glenn Research Center exploits this fundamental acoustic discovery: a specially shaped cavity oscillated at the contained fluid's resonant frequency produces high-amplitude acoustic pressure waves of a magnitude approaching those required of today's seals. While the original researchers are continuing their development of acoustic pumps, refrigeration compressors, and electronic thermal management systems using this technology, the goal of researchers at Glenn is to apply these acoustic principles to a revolutionary sealing device. When the acoustic resonator shape is optimized for the sealing device, the flow from a high-pressure cavity to a low-pressure cavity will be restricted by a series of high-amplitude standing pressure waves of higher pressure than the pressure to be sealed. Since the sealing resonator cavity will not touch the adjacent sealing structures, seal wear will be eliminated, improving system life. Under a cooperative agreement between Glenn and the Ohio Aerospace Institute (OAI), an acoustic-based pre-prototype seal was demonstrated for the first time. A pressurized cavity was attached to one end of the resonator while the other end remained open to ambient pressure. Measurements were taken at several values of applied pressure with the assembly stationary, oscillated at an off-resonance frequency, and then oscillated on-resonance. The three cases show that the flow through the conical resonator can be reduced by oscillating the resonator at the resonance frequency of the air contained within the cavity. The results are currently being compared with results obtained from a commercial computational fluid dynamics code. The objective is to improve the design through numerical simulation before fabricating a next-generation prototype sealing device. Future work is aimed at implementing acoustic seal design improvements to further reduce the leakage flow rate through the device and at reducing the device's overall size.
A digital computer program known as SATANS (static and transient analysis, nonlinear, shells) for the geometrically nonlinear static and dynamic response of arbitrarily loaded shells of revolution is presented. Instructions for the preparation of the input data cards and other information necessary for the operation of the program are described in detail and two sample problems are included. The governing partial differential equations are based upon Sanders' nonlinear thin shell theory for the conditions of small strains and moderately small rotations. The governing equations are reduced to uncoupled sets of four linear, second order, partial differential equations in the meridional and time coordinates by expanding the dependent variables in a Fourier sine or cosine series in the circumferential coordinate and treating the nonlinear modal coupling terms as pseudo loads. The derivatives with respect to the meridional coordinate are approximated by central finite differences, and the displacement accelerations are approximated by the implicit Houbolt backward difference scheme with a constant time interval. The boundaries of the shell may be closed, free, fixed, or elastically restrained. The program is coded in the FORTRAN 4 language and is dimensioned to allow a maximum of 10 arbitrary Fourier harmonics and a maximum product of the total number of meridional stations and the total number of Fourier harmonics of 200. The program requires 155,000 bytes of core storage.
Aerospace structures and spacecraft are a complex assemblage of structural components that are subjected to a variety of complex, cyclic, and transient loading conditions. Significant modeling uncertainties are present in these structures, in addition to the inherent randomness of material properties and loads. To properly account for these uncertainties in evaluating and assessing the reliability of these components and structures, probabilistic structural mechanics (PSM) procedures must be used. Much research has focused on basic theory development and the development of approximate analytic solution methods in random vibrations and structural reliability. Practical application of PSM methods was hampered by their computationally intense nature. Solution of PSM problems requires repeated analyses of structures that are often large, and exhibit nonlinear and/or dynamic response behavior. These methods are all inherently parallel and ideally suited to implementation on parallel processing computers. New hardware architectures and innovative control software and solution methodologies are needed to make solution of large scale PSM problems practical.
The stability of the two-dimensional flow induced by the tranverse oscillation of a cylinder in a viscous fluid is investigated in both the linear and weakly nonlinear regime. The major assumption that is made to simplify the problem is that the oscillation frequency is large in which case an unsteady boundary layer is set up on the cylinder. Results are given for cylinders of elliptic cross section and it is found that for any given eccentricity the most dangerous configuration is when the cylinder oscillates parallel to its minor axis. Some discussion of nonlinear effects is also given and for the circular cylinder it is shown that the steady streaming boundary layer of the basic flow is significantly altered by the instability.
The stability of the two-dimensional flow induced by the transverse oscillation of a cylinder in a viscous fluid is investigated in both the linear and weakly nonlinear regime. The major assumption that is made to simplify the problem is that the oscillation frequency is large in which case an unsteady boundary layer is set up on the cylinder. Results are given for cylinders of elliptic cross section and it is found that for any given eccentricity the most dangerous configuration is when the cylinder oscillates parallel to its minor axis. Some discussion of nonlinear effects is also given and for the circular cylinder it is shown that the steady streaming boundary layer of the basic flow is significantly altered by the instability.
We present a numerical study of spatially quasi-periodic gravity-capillary waves of finite depth in both the initial value problem and travelling wave settings. We adopt a quasi-periodic conformal mapping formulation of the Euler equations, where one-dimensional quasi-periodic functions are represented by periodic functions on a higher-dimensional torus. We compute the time evolution of free surface waves in the presence of a background flow and a quasi-periodic bottom boundary and observe the formation of quasi-periodic patterns on the free surface. Two types of quasi-periodic travelling waves are computed: small-amplitude waves bifurcating from the zero-amplitude solution and larger-amplitude waves bifurcating from finite-amplitude periodic travelling waves. We derive weakly nonlinear approximations of the first type and investigate the associated small-divisor problem. We find that waves of the second type exhibit striking nonlinear behaviour, e.g. the peaks and troughs are shifted non-periodically from the corresponding periodic waves due to the activation of quasi-periodic modes.
The formulation for studying nonlinear Alfvenic states, sustainable in Hall Magnetohydrodynamics (HMHD), becomes considerably simpler and more tractable when circularly polarized Beltrami vectors (the eigenstates of linear HMHD) are used as the basis functions. Nonlinear HMHD is, then, reduced to a rather simple looking set of scalar equations from which a model problem of three interacting Beltrami modes is formulated and analytically solved. The triplet interactions span a variety of familiar nonlinear processes leading to a redistribution as well as periodic exchange of energy. The energy exchange processes (whose strength is measured by an energy exchange/depletion time) will, perhaps, play a dominant role in determining the spectral content of an eventual Alfvenic state. All nonlinearities (sensitive functions of the interacting wave vectors) operate at par, and none is dominant over any substantial region of k-space; their intricate interplay prevents a “universal” picture from emerging; few generalizations on the processes that may, for instance, lead to a turbulent state, are possible. However, the theory can definitely claim: (1) the energy tends to flow from lower to higher k and (2) the higher k z (in the direction of the ambient magnetic field) components of a mode with a given k are depleted/oscillate faster—in some cases much faster. It is noteworthy that the mode coupling is the strongest (with the shortest depletion time) when the participating wave vectors are nearly perpendicular; perhaps, an expected consequence of the curl (cross product) nonlinearities. Numerical simulations will be necessary to help create a fully reliable picture.
A numerical procedure has been developed to investigate the nonlinear and strain rate dependent deformation response of polymer matrix composite laminated plates under high strain rate impact loadings. A recently developed strength of materials based micromechanics model, incorporating a set of nonlinear, strain rate dependent constitutive equations for the polymer matrix, is extended to account for the transverse shear effects during impact. Four different assumptions of transverse shear deformation are investigated in order to improve the developed strain rate dependent micromechanics model. The validities of these assumptions are investigated using numerical and theoretical approaches. A method to determine through the thickness strain and transverse Poisson's ratio of the composite is developed. The revised micromechanics model is then implemented into a higher order laminated plate theory which is modified to include the effects of inelastic strains. Parametric studies are conducted to investigate the mechanical response of composite plates under high strain rate loadings. Results show the transverse shear stresses cannot be neglected in the impact problem. A significant level of strain rate dependency and material nonlinearity is found in the deformation response of representative composite specimens.
The Newton Modified Barrier Method (NMBM) is applied to structural optimization problems with large a number of design variables and constraints. This nonlinear mathematical programming algorithm was based on the Modified Barrier Function (MBF) theory and the Newton method for unconstrained optimization. The distinctive feature of the NMBM method is the rate of convergence that is due to the fact that the design remains in the Newton area after each Lagrange multiplier update. This convergence characteristic is illustrated by application to structural problems with a varying number of design variables and constraints. The results are compared with those obtained by optimality criteria (OC) methods and by the ASTROS program.
The familiar suboptimal regulator design approach is recast as a constrained optimization problem and incorporated in a Computer Aided Design (CAD) package where both design objective and constraints are quadratic cost functions. This formulation permits the separate consideration of, for example, model following errors, sensitivity measures and control energy as objectives to be minimized or limits to be observed. Efficient techniques for computing the interrelated cost functions and their gradients are utilized in conjunction with a nonlinear programming algorithm. The effectiveness of the approach and the degree of insight into the problem which it affords is illustrated in a helicopter regulation design example.
The problem of low-speed impact of a one-dimensional sandwich panel by a rigid cylindrical projectile is considered. The core of the sandwich panel is functionally graded such that the density, and hence its stiffness, vary through the thickness. The problem is a combination of static contact problem and dynamic response of the sandwich panel obtained via a simple nonlinear spring-mass model (quasi-static approximation). The variation of core Young's modulus is represented by a polynomial in the thickness coordinate, but the Poisson's ratio is kept constant. The two-dimensional elasticity equations for the plane sandwich structure are solved using a combination of Fourier series and Galerkin method. The contact problem is solved using the assumed contact stress distribution method. For the impact problem we used a simple dynamic model based on quasi-static behavior of the panel - the sandwich beam was modeled as a combination of two springs, a linear spring to account for the global deflection and a nonlinear spring to represent the local indentation effects. Results indicate that the contact stiffness of thc beam with graded core Increases causing the contact stresses and other stress components in the vicinity of contact to increase. However, the values of maximum strains corresponding to the maximum impact load arc reduced considerably due to grading of thc core properties. For a better comparison, the thickness of the functionally graded cores was chosen such that the flexural stiffness was equal to that of a beam with homogeneous core. The results indicate that functionally graded cores can be used effectively to mitigate or completely prevent impact damage in sandwich composites.
The problem of the exitation of obliquely propagating magnetosonic waves which can steepen up (also known as shocklets) is considered. Shocklets have been observed upstream of the Earth's bow shock and at comets Giacobini-Zinner and Grigg-Skjellerup. Linear theory as well as two-dimensional (2-D) hybrid (fluid electrons, particle ions) simulations are used to determine the properties of waves generated by ring-beam velocity distributions in great detail. The effects of both proton and oxygen ring-beams are considered. The study of instabilities excited by a proton ring-beam is relevant to the region upstream of the Earth's bow shock, whereas the oxygen ring-beam corresponds to cometary ions picked up by the solar wind. Linear theory has shown that for a ring-beam, four instabilities are found, one on the nonresonant mode, one on the Alfven mode, and two along the magnetosonic/whistler branch. The relative growth rate of these instabilities is a sensitive function of parameters. Although one of the magnetosonic instabilities has maximum growth along the magnetic field, the other has maximum growth in oblique directions. We have studied the competition of these instabilities in the nonlinear regime using 2-D simulations. As in the linear limit, the nonlinear results are a function of beam density and distribution function. By performing the simulations as both initial value and driven systems, we have found that the outcome of the simulations can vary, suggesting that the latter type simulations is needed to address the observations. A general conclusion of the simulation results is that field-aligned beams do not result in the formation of shocklets, whereas ring-beam distributions can.
An attempt is made to solve the problem of orbital transfer between coplanar elliptical orbits. Pure propulsive transfer is analyzed with the restriction of two impulses for the transfer. The optimal switching conditions are reviewed, and it is shown that the solution is obtained by solving a set of three nonlinear equations for three unknowns. A semianalytical solution is obtained to the problem of planar rotation of an orbit for the pure propulsive maneuver, and it is shown that, for high eccentricity and rotation angle, aeroassisted transfer is a fuel-saving maneuver. It is demonstrated that complete circularization of the intermediate orbit is not necessary in the optimal aeroassisted transfer. An analytical proof is presented, giving an explicit condition for noncircularization. A complete numerical solution is presented for a case of optimal aeroassisted transfer from a low-energy orbit to a high-energy orbit.
This study presents a connected vehicles (CVs)-based traffic signal optimization framework for a coordinated arterial corridor. The signal optimization and coordination problem are first formulated in a centralized scheme as a mixed-integer nonlinear program (MINLP). The optimal phase durations and offsets are solved together by minimizing fuel consumption and travel time considering an individual vehicle’s trajectories. Due to the complexity of the model, we decompose the problem into two levels: an intersection level to optimize phase durations using dynamic programming (DP), and a corridor level to optimize the offsets of all intersections. In order to solve the two-level model, a prediction-based solution technique is developed. The proposed models are tested using traffic simulation under various scenarios. Compared with the traditional actuated signal timing and coordination plan, the signal timing plans generated by solving the MINLP and the two-level model can reasonably improve the signal control performance. When considering varies vehicle types under high demand levels, the proposed two-level model reduced the total system cost by 3.8% comparing to baseline actuated plan. MINLP reduced the system cost by 5.9%. It also suggested that coordination scheme was beneficial to corridors with relatively high demand levels. For intersections with major and minor street, coordination conducted for major street had little impacts on the vehicles at the minor street.
Drones are projected to alter last-mile delivery, but their short travel range is a concern. In this study, we propose a drone delivery network design using automated battery swapping machines (ABSMs) to extend ranges. The design minimizes the long-term delivery costs, including ABSM investment, drone ownership, and cost of the delivery time, and locates ABSMs to serve a set of customers. We build a mixed-integer nonlinear program that captures the nonlinear waiting time of drones at ABSMs. To solve the problem, we create an exact solution algorithm that finds the globally optimal solution using a derivative-supported cutting-plane method. To validate the applicability of our program, we conduct a case study on the Chicago Metropolitan area using cost data from leading ABSM manufacturer and geographical data from the planning and operations language for agent-based regional integrated simulation (more commonly known as POLARIS). A sensitivity analysis identifies that ABSM service times and costs are the key parameters impacting the long-term adoption of drone delivery.
This paper presents the formulation of simultaneous state and parameter estimation problems for flexible structures in terms of least-squares minimization problems. The approach combines an on-line order determination algorithm, with least-squares algorithms for finding estimates of modal approximation functions, modal amplitudes, and modal parameters. The approach combines previous results on separable nonlinear least squares estimation with a regression analysis formulation of the state estimation problem. The technique makes use of sequential Householder transformations. This allows for sequential accumulation of matrices required during the identification process. The technique is used to identify the modal prameters of a flexible beam.
A static flexible fighter aircraft wing configuration is presently addressed by a multilevel optimization technique, based on both a full-potential concept and a rapid structural optimization program, which can be applied to such aircraft-design problems as maneuver load control, aileron reversal, and lift effectiveness. It is found that nonlinearities are important in the design of an aircraft whose flight envelope encompasses the transonic regime, and that the present structural suboptimization produces a significantly lighter wing by reducing ply thicknesses.