Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal modes”

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 181 records · Page 10

The acoustic field of a point source in a uniform boundary layer over an impedance plane

The acoustic field of a point source in a boundary layer above an impedance plane is investigated anatytically using Obukhov quasi-potential functions, extending the normal-mode theory of Chunchuzov (1984) to account for the effects of finite ground-plane impedance and source height. The solution is found to be asymptotic to the surface-wave term studies by Wenzel (1974) in the limit of vanishing wind speed, suggesting that normal-mode theory can be used to model the effects of an atmospheric boundary layer on infrasonic sound radiation. Model predictions are derived for noise-generation data obtained by Willshire (1985) at the Medicine Bow wind-turbine facility. Long-range downwind propagation is found to behave as a cylindrical wave, with attention proportional to the wind speed, the boundary-layer displacement thickness, the real part of the ground admittance, and the square of the frequency.

Zorumski, W. E.↗

Effects of joint rate and displacement constraints on stability regions

A block diagram showing plant dynamics in normal mode coordinates, a linear feedback controller with command limits, and actuator dynamics with rate and displacement limits is given. The objective is to examine the effects of joint rate and displacement saturation limits on the stability regions. An unstable short period dynamic model is transformed into normal mode coordinates and is augmented with actuator dynamics. A linear feedback controller is used to provide closed-loop stability. The stability with constrained actuator rate limits under varying bandwidth, displacement and command limits is examined.

Shrivastava, P. C.↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

Interpretation of dispersion relations for bounded systems.

Treatment of the problem of constructing normal modes for an arbitrarily bounded system from roots of the linear dispersion relation D(omega, k) = 0 for the corresponding infinite or periodically bounded system. For a system described by continuous macroscopic variables, and of general cylindrical form, each transverse eigenmode gives rise to a set of axial normal modes constructed from a pair of dominant roots of D = 0 satisfying the boundary conditions which are characterized by complex reflection coefficients for the dominant waves. The implications of the results for the interpretation of experiments on plasma waves and instabilities on finite cylinders are discussed, with particular reference to the effects of end-plate damping and axial current on Q-machines.

Rognlien, T. D.↗

Propulsion Aeroelastic Analysis Developed for Flutter and Forced Response

The NASA Glenn Research Center at Lewis Field develops new technologies to increase the fuel efficiency of aircraft engines, improve the safety of engine operation, reduce emissions, and reduce engine noise. With the development of new designs for fans, compressors, and turbines to achieve these goals, the basic aeroelastic requirements are that there should be no flutter (self-excited vibrations) or high resonant blade stresses (due to forced response) in the operating regime. Therefore, an accurate prediction and analysis capability is required to verify the aeroelastic soundness of the designs. Such a three-dimensional viscous propulsion aeroelastic analysis capability has been developed at Glenn with support from the Advanced Subsonic Technology (AST) program. This newly developed aeroelastic analysis capability is based on TURBO, a threedimensional unsteady aerodynamic Reynolds-averaged Navier-Stokes turbomachinery code developed previously under a grant from Glenn. TURBO can model the viscous flow effects that play an important role in certain aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), flutter in the presence of shock and boundary-layer interaction, and forced response due to wakes and shock impingement. In aeroelastic analysis, the structural dynamics representation of the blades is based on normal modes. A finite-element analysis code is used to calculate these in-vacuum vibration modes and the associated natural frequencies. In an aeroelastic analysis using the TURBO code, flutter and forced response are modeled as being uncoupled. To calculate if a blade row will flutter, one prescribes the motion of the blade to be a harmonic vibration in a specified in-vacuum normal mode. An aeroelastic analysis preprocessor is used to generate the displacement field required for the analysis. The work done by aerodynamic forces on the vibrating blade during a cycle of vibration is calculated. If this work is positive, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the blade s oscillation amplitude. The forced-response excitations on a blade row are calculated by modeling the flow through two adjacent blade rows using the TURBO code. The blades are assumed to be rigid. As an option, a single blade row can be modeled with the upstream blade row influence represented by a time-varying disturbance (gust) at the inlet boundary. The unsteady forces on a blade row from such analyses are used in a structural analysis along with the blade structural dynamics characteristics and aerodynamic damping associated with blade vibration to calculate the resulting dynamic stresses on the blade.

Bakhle, Milind A.↗

The initial growth of disturbances in a baroclinic flow

The growth of perturbations in a baroclinic flow is examined as an initial value problem. Although the long time asymptotic behavior is dominated by discrete exponentially growing normal modes when they exist, these do not form a complete set and initial intensification is shown to be dependent on the continuous spectrum. The vertical structure of perturbations emerges as an important influence on initial growth, and physically realistic disturbances are shown to grow to amplitudes where nonlinear effects are important before obtaining normal mode form. Connection is made with the work of Arnol'd (1965) and Blumen (1968) and the numerical experiments of Simmons and Hoskins (1979). Application of these results to cyclogenesis in geographically fixed areas is suggested and implied constraints on numerical models discussed.

Farrell, B. F.↗

A mechanism for magnetospheric substorms

Energy-principle analysis performed on two-dimensional, self-consistent solutions for magnetospheric convection indicates that the magnetosphere is unstable to isobaric (yet still frozen-in) fluctuations of plasma-sheet flux tubes. Normally, pdV work associated with compression maintains stability of the inward/outward oscillating normal mode. However, if Earth's ionosphere can provide sufficient mass flux, isobaric expansion of flux tubes can occur. The growth of a field-aligned potential drop in the near-Earth, midnight portion of the plasma sheet, associated with upward field-aligned currents responsible for the Harang discontinuity, redistributes plasma along field lines in a manner that destabilizes the normal mode. The growth of this unstable mode results in an out-of-equilibrium situation near the inner edge. When this occurs over a downtail extent comparable to the half-thickness of the plasma sheet, collapse ensues and forces thinning of the plasma sheet whereby conditions favorable to reconnection occur. This scenario for substorm onset is consistent with observed upward fluxes of ions, parallel potential drops, and observations of substorm onset. These observations include near Earth onset, pseudobreakups, the substorm current wedge, and local variations of plasma-sheet thickness.

Erickson, G. M.↗

The Effects of Magnetic-field Geometry on Longitudinal Oscillaitons of Solar Prominences

We investigate the influence of the geometry of the solar filament magnetic structure on the large-amplitude longitudinal oscillations. A representative filament flux tube is modeled as composed of a cool thread centered in a dipped part with hot coronal regions on either side.We have found the normal modes of the system and establish that the observed longitudinal oscillations are well described with the fundamental mode. For small and intermediate curvature radii and moderate to large density contrast between the prominence and the corona, the main restoring force is the solar gravity. In this full wave description of the oscillation a simple expression for the oscillation frequencies is derived in which the pressure-driven term introduces a small correction. We have also found that the normal modes are almost independent of the geometry of the hot regions of the tube. We conclude that observed large-amplitude longitudinal oscillations are driven by the projected gravity along the flux tubes and are strongly influenced by the curvature of the dips of the magnetic field in which the threads reside.

Luna, M.↗

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗

Development of mid-infrared solid state lasers for spaceborne lidar

The laser performance of a Ho(3+):Tm(3+):Cr(3+):YAG crystal was measured under flashlamp pumping at various operating temperatures. The normal mode laser thresholds of a Ho(3+)(0.45 at. percent):Tm(3+)(2.5 at. percent):Cr(3+)(1.5 at. percent ):YAG crystal were found to range from 26 to 50 J between 120 and 200 K with slope efficiencies up to 0.36 percent with a 60 percent reflective output mirror. From the Q-switched operations a slope efficiency corresponding to 90 percent of the normal mode operation was observed. Laser wavelengths were measured for various operating conditions and fluorescence spectra were obtained at various temperatures in order to help understand the dynamic energy processes among the Ho(3+), Tm(3+), and Cr(3+) ions. A pulse forming network for a flashlamp pumped Cr:GSAG laser, which is to be used as a high power laser diode simulator in rare earth laser pumping, was completed and tested. The network provided critically damped, 1 ms FWHM, square pulses with a rise time of about 160 micrometers at an input electrical energy of 300 J.

Whitney, Donald A.↗

Local Helioseismology: A New Window into the Sun

In recent years the focus of helioseismology has shifted from the classical measurements of global properties using normal mode frequencies to possibly a more difficult problem of measuring the three-dimensional, time-dependent structure of the solar interior, now called local helioseismology. Problems to be studed include the subsurface structure of observed surface features such as sunspots, convective flow patterns, and the magnetic field of the 1 1-year sunspot cycle. A variety of techniques have been developed to attack these problems, including a 'local' normal mode method, a method of measuring travel times between &stant locations, and a method for holographically reconstructing the acoustic signal in the solar interior. Recent results from the field will be reviewed.

Duvall, Thomas↗

Wave-interactions in supersonic and hypersonic flows

Work completed under the current grant comprises the start of a theoretical and computational attack on the subharmonic route to secondary instabilities in compressible flows. The total flow field in this problem is made up of the following components: (1) a steady streamwise mean boundary layer flow which depends only on the normal space component y; (2) a two-dimensional time dependent T-S wave which moves with wavespeed c and has no spanwise dependence; and (3) a fully three-dimensional, time dependent T-S wave whose streamwise wavenumber is half of the streamwise wavenumber associated with the two-dimensional T-S wave in b. If a frame of reference is adopted which moves with the wavespeed c of the 2-D T-S wave, the time dependence of this portion of the flow can be eliminated. The effective steady mean flow in this problem is now the sum of the original parallel steady mean flow and the initial 2-D T-S instability. Dependence on the streamwise coordinate x in this mean flow can be extracted by assuming normal mode expansions involving complex exponentials and the streamwise wavenumber a. However, it is important to note that, because this is a wave-wave interaction problem, unlike the usual linear instability case, both the complex exponential, and its complex conjugate, must be retained in describing the 2-D T-S wave. The role of the perturbation to the steady mean flow is now played by the 3-D time dependent T-S wave. In treating this wave, normal modes in the streamwise and spanwise directions and time may be used. Consistent with the subharmonic nature of this transition route, the streamwise wavenumber is a/2, and complex conjugates of the complex exponential must be employed. This is not the case with the modes giving z and t dependence with wavespeed o and spanwise wavenumber B as the effective mean flow quantities are independent of z and their time dependence is accounted for by the moving frame of reference. Consequently, the wave-wave interaction which will produce mean flow modification occurs only through the streamwise exponentials.

Lakin, William D.↗

Development of mid-infrared solid state lasers for spaceborne lidar

Laser performance of Ho(3+):Tm(3+):Cr(3+):YAG crystals was investigated under both Cr:GSAG laser and flashlamp pumping. A flashlamp pumped Cr:GSAG laser was built to simulate high power quasi-CW laser diode pumping of a 2.1 micrometer holmium laser. The 2.1 micrometer output laser energy exceeded more than 14 mJ, the highest value reported to date under laser pumping near 785 nm. This was obtained in a pulse length of nearly 650 microsec from a 3 x 3 mm Ho:Tm:Cr:YAG rod by using the flashlamp-pumped Cr:CSAG laser as a pumping source at the diode laser wavelength, 785 micrometers. In addition, Ho:Tm:Cr:YAG crystals with various Tm(3+) concentrations were evaluated for flashlamp-pumped normal mode and Q-switched 2.1 micrometer laser operations under a wide variety of experimental conditions in order to understand internal dynamic processes among the ions and to determine an optimum lasing condition. An increase of the laser slope efficiency was observed with the increase of the Tm(3+) concentration from 2.5 to 4.5 atomic percent. The thermal dependence of the laser performance was also investigated. Q-switched laser output energies corresponding to nearly 100 percent of the normal-mode laser energies were obtained in a strong single spike of 200 ns pulse length by optimizing the opening time of a lithium niobate Q-switch.

Whitney, Donald A.↗

Development of mid-infrared solid state lasers for spaceborne lidar

Researchers investigated laser performance of Ho(3+):Tm(3+):Cr(3+):YAG crystals under both Cr:GSAG laser and flashlamp pumping. A flashlamp pumped Cr:GSAG laser was built to simulate high power quasi-CW laser diode pumping of a 2.1 micron holmium laser. The 2.1 micron output laser energy exceeded more than 14 mJ, the highest value reported to date under laser pumping near 785 nm. This was obtained in a pulse length of nearly 650 microns from a 3 x 3 mm Ho:Tm:Cr:YAG rod by using the flashlamp-pumped Cr:GSAG laser as a pumping source at the diode laser wavelength, 785 microns. In addition, Ho:Tm:Cr:YAG crystals with various Tm(3+) concentrations have been evaluated for flashlamp-pumped normal mode and Q-switched 2.1 micron laser operations under a wide variety of experimental conditions in order to understand internal dynamic processes among the ions and to determine an optimum lasing condition. An increase of the laser slope efficiency was observed with the increase of the Tm(3+) concentration from 2.5 atomic percent to 4.5 atomic percent. The thermal dependence of the laser performance was also investigated. Q-switched laser output energies corresponding to nearly 100 percent of the normal-mode laser energies were obtained in a strong single spike of 200 ns pulse length by optimizing the opening time of a lithium niobate Q-switch.

Whitney, Donald A.↗

Flutter of a uniform wing with an arbitrarily placed mass according to a differential-equation analysis and a comparison with experiment

A method is presented for the calculation of the flutter speed of a uniform wing carrying an arbitrarily placed concentrated mass. The method, an extension of recently published work by Goland and Luke, involves the solution of the differential equations of motion of the wing at flutter speed and therefore does not require the assumption of specific normal modes of vibration. The order of the flutter determinant to be solved by this method depends upon the order of the system of differential equations and not upon the number of modes of vibration involved. The differential equations are solved by operational methods, and a brief discussion of operational methods as applied to boundary-value problems is included in one of two appendixes. A comparison is made with experiment for a wing with a large eccentrically mounted weight and good agreement is obtained. Sample calculations are presented to illustrate the method; and curves of amplitudes of displacement, torque, and shear for a particular case are compared with corresponding curves computed from the first uncoupled normal modes.

Runyan, Harry L↗

The wave pattern produced by a point source on a rotating disk

It is pointed out that the boundary layer on a rotating disk is important in stability theory because it provides a particularly simple way to study the important phenomenon of crossflow instability. This type of instability is responsible for early transition on sweptback wings. Mack and Kendall (1983) have studied the wave patterns formed by harmonic point sources in a Blasius boundary layer on the basis that the source uniformly excites all oblique normal modes of the source frequency. The calculation procedure for planar boundary layers was modified to fit the different geometry of the rotating disk and the lack of an axis of symmetry. Calculations were performed of the wave pattern produced by a zero-frequency point source located at the Reynolds number of the artificial roughness element in an experiment conducted by Wilkinson and Malik (1983). The results provided in the present investigation confirm that the experimental wave pattern is a superposition of the complete azimuthal wavenumber spectrum of zero-frequency normal modes with uniform initial amplitude and phase.

Mack, L. M.↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗