Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiff 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 487 records · Page 27

Dynamics of tilting proprotor aircraft in cruise flight

A nine degree-of-freedom theoretical model is developed for investigations of the dynamics of a proprotor operating in high inflow axial flight on a cantilever wing. The basic characteristics of the rotor high inflow aerodynamics and the resulting rotor aeroelastic behavior are discussed. The problems of classical whirl flutter, the two-bladed rotor, and the influence of the proprotor on the stability derivatives of the aircraft are treated briefly. The influence of various elements of the theoretical model is discussed, including the modeling used for the blade and wing aerodynamics, and the influence of the rotor lag degree of freedom. The results from tests of two full-scale proprotors - a gimballed, stiff-inplane rotor and a hingeless, soft-inplane rotor - are presented; comparisons with the theoretical results show good correlation.

Johnson, W.↗

A pressurized cylindrical shell with a fixed end which contains an axial part-through or through crack

A cylindrical shell having a very stiff and plate or a flange is considered. It is assumed that near the end the cylinder contains an axial flaw which may be modeled as a part through surface crack or a through crack. The effect of the end constraining on the stress intensity factor which is the main fracture mechanics parameter is studied. The applied loads acting on the cylinder are assumed to be axisymmetric. Thus the crack problem under consideration is symmetric with respect to the plane of the crack and consequently only the Mode 1 stress intensity factors are nonzero. With this limitation, the general perturbation problem for a cylinder with a built in end containing an axial crack is considered. Reissner's shell theory is used to formulate the problem. The part through crack problem is treated by using a line spring model. In the case of a crack tip terminating at the fixed end it is shown that the integral equations of the shell problem has the same generalized Cauchy kernel as the corresponding plane stress elasticity problem.

Yahsi, O. S.↗

The response of Helmholtz resonators to external excitation. I - Single resonators

The response of single two- and three-dimensional Helmholtz resonators subject to external excitation by a plane acoustic wave is studied. The inviscid linearized problem is solved by the matched-asymptotic-expansion technique in the low-frequency limit, i.e. when the characteristic neck dimension is small compared with the acoustic wavelength. The scale of the cavity is chosen such as to tune the system to the lowest Helmholtz mode. The results are compared with the classical lumped-element model, and refined expressions are derived for the effective mass (added length) and stiffness (effective cavity volume) of the resonator.

Monkewitz, P. A.↗

Model reduction in the physical coordinate system

In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.

Yae, K. Harold↗

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 2: Solutions and numerical results

This paper discusses an incremental-iterative nonlinear solution technique for solving the nonlinear finite element equations of the superconducting state of a superconductor. The untreated equations are highly ill-conditioned and are impossible to solve within the typical 16-place double precision supplied by most computers. A combination of matrix scaling and mesh grading techniques is used to reduce the condition number of the tangent stiffness matrix and increase the accuracy of the current carrying boundary layer representation. Numerical results for a one-dimensional model of a time-independent superconductor treated by the Ginzburg-Landau model are presented and discussed. The computed solutions clearly display the Meissner effect of magnetic field expulsion from the central region of the superconductor. These results are compared to the physics of a low-viscosity fluid problem. From this analogy, a physical argument is advanced about the macroscopic behavior of superconductors.

Schuler, James J.↗

Metal Foam Analysis: Improving Sandwich Structure Technology for Engine Fan and Propeller Blades

The Life Prediction Branch of the NASA Glenn Research Center is searching for ways to construct aircraft and rotorcraft engine fan and propeller blades that are lighter and less costly. One possible design is to create a sandwich structure composed of two metal faces sheets and a metal foam core. The face sheets would carry the bending loads and the foam core would have to resist the transverse shear loads. Metal foam is ideal because of its low density and energy absorption capabilities, making the structure lighter, yet still stiff. The material chosen for the face sheets and core was 17-4PH stainless steel, which is easy to make and has appealing mechanical properties. This material can be made inexpensively compared to titanium and polymer matrix composites, the two current fan blade alternatives. Initial tests were performed on design models, including vibration and stress analysis. These tests revealed that the design is competitive with existing designs; however, some problems were apparent that must be addressed before it can be implemented in new technology. The foam did not hold up as well as expected under stress. This could be due to a number of issues, but was most likely a result of a large number of pores within the steel that weakened the structure. The brazing between the face sheets and the foam was also identified as a concern. The braze did not hold up well under shear stress causing the foam to break away from the face sheets. My role in this project was to analyze different options for improving the design. I primarily spent my time examining various foam samples, created with different sintering conditions, to see which exhibited the most favorable characteristics for our purpose. Methods of analysis that I employed included examining strut integrity under a microscope, counting the number of cells per inch, measuring the density, testing the microhardness, and testing the strength under compression. Shear testing will also be done to examine the strengths of different types of brazes.

Fedor, Jessica L.↗

Development of the Aquarius Antenna Deployment Mechanisms and Spring/Damper Actuator

The Aquarius Instrument s large radar reflector dish needed to be stowed for launch, and then deployed on-orbit. The Deployment Subsystem consisted of a cantilevered boom structure and two single-axis hinge mechanisms to accurately deploy and position the reflector dish relative to the radar feed horns. The cantilevered design demanded high stiffness and accuracy from the deployment mechanism at the root of the boom. A preload-generating end-of-travel latch was also required. To largely eliminate the need for control systems, each deployment mechanism was actuated by a passive spring motor with viscous-fluid damping. Tough requirements and adaptation of a heritage actuator to the new application resulted in numerous challenges. Fabrication, assembly, and testing encountered additional problems, though ultimately the system was demonstrated very successfully. This paper revisits the development to highlight which design concepts worked and the many important lessons learned.

Johnson, Joel A.↗

A parameter identification method for the rotordynamic coefficients of a high Reynolds number hydrostatic bearing

In identifying the rotordynamic coefficients of a high-Reynolds-number hydrostatic bearing, fluid-flow induced forces present a unique problem, in that they provide an unmeasureable and uncontrollable excitation to the bearing. An analysis method is developed that effectively eliminates the effects of fluid-flow induced excitation on the estimation of the bearing rotordynamic coefficients, by using power spectral densities. In addition to the theoretical development, the method is verified experimentally by single-frequency testing, and repeatability tests. Results obtained for a bearing are the twelve rotordynamic coefficients (stiffness, damping, and inertia coefficients) as functions of eccentricity ratio, speed, and supply pressure.

Rouvas, C.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Studies in tilt rotor VTOL aircraft aeroelasticity, volume 2

Two methods for natural mode vibration analysis are discussed. The first consists of a direct approach based on a finite element representation of the complete structure as an entity. The mass and stiffness matrices for the complete structure are assembled by properly combining the mass and stiffness matrices of the individual elements into which the structure has been divided. The second approach is that of component mode synthesis. This method is based on the concept of synthesizing the natural modes of the complete structure from modes of conveniently difined substructures, or components, into which the structure has been partitioned. In this way the expedient of reducing the system degrees of freedom, and thus the size of the eigenvalue problem, can be introduced by partial modal synthesis.

Kvaternik, R. G.↗

Flow-induced vibrations of the SSME LOX inlet tee vanes

This paper describes the results of a program initiated to identify the cause and the mechanism of 4-kHz vibrations of the Space Shuttle main engine's (SSME's) LOX inlet tee, which were revealed in 17 percent of the SSMEs during ground tests. The results of the engines' inspection, structural analyses, computational fluid dynamics studies, literature reviews, expert consultations, and laboratory testings identified the probable cause of the 4-kHz vibrations being related to the fact that the walls of some tees were thinner than in the mojority of tees, resulting in less stiff structures. The mechanism of the vibrations was identified as vortex shedding in combination with a structural coupling with the thrust cone, which increased the vane mobility, which, in turn, aggravated the fluid coupling with the thrust cone. As a result, the vane mobility was increased, aggravating the fluid coupling with the vanes. An internal modification was developed and was successfully demonstrated on an engine that had previously exhibited a 4-kHz problem.

O'Connor, G. M.↗

A numerical method for interface problems in elastodynamics

The numerical implementation of a formulation for a class of interface problems in elastodynamics is discussed. This formulation combines the use of the finite element and boundary integral methods to represent the interior and the exteriro regions, respectively. In particular, the response of a semicylindrical alluvial valley in a homogeneous halfspace to incident antiplane SH waves is considered to determine the accuracy and convergence of the numerical procedure. Numerical results are obtained from several combinations of the incidence angle, frequency of excitation, and relative stiffness between the inclusion and the surrounding halfspace. The results tend to confirm the theoretical estimates that the convergence is of the order H(2) for the piecewise linear elements used. It was also observed that the accuracy descreases as the frequency of excitation increases or as the relative stiffness of the inclusion decreases.

Mcghee, D. S.↗

Overcoming Geometry-Induced Stiffness with IMplicit-Explicit (IMEX) Runge-Kutta Algorithms on Unstructured Grids with Applications to CEM, CFD, and CAA

My goal is to develop and implement efficient, accurate, and robust Implicit-Explicit Runge-Kutta (IMEX RK) methods [9] for overcoming geometry-induced stiffness with applications to computational electromagnetics (CEM), computational fluid dynamics (CFD) and computational aeroacoustics (CAA). IMEX algorithms solve the non-stiff portions of the domain using explicit methods, and isolate and solve the more expensive stiff portions using implicit methods. Current algorithms in CEM can only simulate purely harmonic (up to lOGHz plane wave) EM scattering by fighter aircraft, which are assumed to be pure metallic shells, and cannot handle the inclusion of coatings, penetration into and radiation out of the aircraft. Efficient MEX RK methods could potentially increase current CEM capabilities by 1-2 orders of magnitude, allowing scientists and engineers to attack more challenging and realistic problems.

Kanevsky, Alex↗

Stability analysis of position and force control problems for robot arms

Stability issues involving the control of a robot arm under the influence of external forces are discussed. Several different scenarios are considered: position control with the external force as an unmodeled disturbance, compliant control for a bounded external force in some subspace, and compliant control for a force due to the interaction with an environment whose dynamical behavior can be modeled. In each of these cases, a stability analysis using the Lyapunov method is presented. An explanation of instability is put forth for the case in which the environment has flexibility and the gains are inappropriately chosen. When the environment is stiff in the force control subspace, robust stability can be achieved with the integral force feedback.

Wen, John T.↗

A pressurized cylindrical shell with a fixed end which contains an axial part-through or through crack

In this paper a cylindrical shell having a very stiff end plate or a flange is considered. It is assumed that near the end the cylinder contains an axial flow which may be modeled as a part-through surface crack or through crack. The primary objective is to study the effect of the end constraining on the stress intensity factor which is the main fracture mechanics parameter. The applied loads acting on the cylinder are assumed to be axisymmetric. Thus the crack problem under consideration is symmetric with respect to the plane of the crack and consequently only the mode I stress intensity factors are nonzero. With this limitation, the general perturbation problem for a cylinder with a built-in end containing an axial crack is considered. Reissner's shell theory is used to formulate the problem. The part-through crack problem is treated by using a line-spring model. In the case of a crack tip terminating at the fixed end it is shown that the integral equation of the shell problem has the same generalized Cauchy kernel as the corresponding plane stress elasticity problem. Even though the problem is formulated for a general surface crack profile and arbitrary crack surface tractions, the numerical results are obtained only for a semielliptic part-through axial crack located at the inside or outside surface of the cylinder and for internal pressure acting on the cylinder. The stress intensity factors are calculated and presented for a relatively wide range of dimensionless length parameters of the problem.

Yahsi, O. S.↗

Investigation of supersonic chemically reacting and radiating channel flow

The 2-D time-dependent Navier-Stokes equations are used to investigate supersonic flows undergoing finite rate chemical reaction and radiation interaction for a hydrogen-air system. The explicit multistage finite volume technique of Jameson is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The multidimensional radiative transfer equations for a nongray model are provided for a general configuration and then reduced for a planar geometry. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the participating species. The supersonic inviscid and viscous, nonreacting flows are solved by employing the finite volume technique of Jameson and the unsplit finite difference scheme of MacCormack. The specified problem considered is of the flow in a channel with a 10 deg compression-expansion ramp. The calculated results are compared with those of an upwind scheme. The problem of chemically reacting and radiating flows are solved for the flow of premixed hydrogen-air through a channel with parallel boundaries, and a channel with a compression corner. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the entire flow field.

Mani, Mortaza↗

Time Frequency Analysis of Spacecraft Propellant Tank Spinning Slosh

Many spacecraft are designed to spin about an axis along the flight path as a means of stabilizing the attitude of the spacecraft via gyroscopic stiffness. Because of the assembly requirements of the spacecraft and the launch vehicle, these spacecraft often spin about an axis corresponding to a minor moment of inertia. In such a case, any perturbation of the spin axis will cause sloshing motions in the liquid propellant tanks that will eventually dissipate enough kinetic energy to cause the spin axis nutation (wobble) to grow further. This spinning slosh and resultant nutation growth is a primary design problem of spinning spacecraft and one that is not easily solved by analysis or simulation only. Testing remains the surest way to address spacecraft nutation growth. This paper describes a test method and data analysis technique that reveal the resonant frequency and damping behavior of liquid motions in a spinning tank. Slosh resonant frequency and damping characteristics are necessary inputs to any accurate numerical dynamic simulation of the spacecraft.

Green, Steven T.↗

The effects of aeroelastic deformation on the unaugmented stopped-rotor dynamics of an X-Wing aircraft

A new design concept in the development of VTOL aircraft with high forward flight speed capability is that of the X-Wing, a stiff, bearingless helicopter rotor system which can be stopped in flight and the blades used as two forward-swept and two aft-swept wings. Because of the usual configuration in the fixed-wing mode, there is a high potential for aeroelastic divergence or flutter and coupling of blade vibration modes with rigid-body modes. An aeroelastic stability analysis of an X-Wing configuration aircraft was undertaken to determine if these problems could exist. This paper reports on the results of dynamic stability analyses in the lateral and longitudinal directions including the vehicle rigid-body and flexible modes. A static aeroelastic analysis using the normal vibration mode equations of motion was performed to determine the cause of a loss of longitudinal static margin with increasing airspeed. This loss of static margin was found to be due to aeroelastic washin of the forward-swept blades and washout of the aft-swept blades moving the aircraft aerodynamic center forward of the center of gravity. This phenomenon is likely to be generic to X-Wing aircraft.

Gilbert, Michael G.↗