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 217 records · Page 12

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.↗

Rapid method for computing the mechanical resonances of irregular objects

A solid object's geometry, density, and elastic moduli completely determine its spectrum of normal modes. Solving the inverse problem—determining a material's elastic moduli given a set of resonance frequencies and sample geometry—relies on the ability to compute resonance spectra accurately and efficiently. Established methods for calculating these spectra are either fast but limited to simple geometries, or are applicable to arbitrarily shaped samples at the cost of being prohibitively slow. Here, we describe a method to rapidly compute the normal modes of irregularly shaped objects using entirely open-source software. Our method's accuracy compares favorably with existing methods for simple geometries and shows a significant improvement in speed over existing methods for irregular geometries.

Acoustics↗

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.↗

A framework to evaluate IMEX schemes for atmospheric models

We present a new evaluation framework for implicit and explicit (IMEX) Runge–Kutta time-stepping schemes. The new framework uses a linearized nonhydrostatic system of normal modes. We utilize the framework to investigate the stability of IMEX methods and their dispersion and dissipation of gravity, Rossby, and acoustic waves. We test the new framework on a variety of IMEX schemes and use it to develop and analyze a set of second-order low-storage IMEX Runge–Kutta methods with a high Courant–Friedrichs–Lewy (CFL) number. We show that the new framework is more selective than the 2-D acoustic system previously used in the literature. Schemes that are stable for the 2-D acoustic system are not stable for the system of normal modes.

54 ENVIRONMENTAL SCIENCES↗

Fluid-Dynamic Mechanisms Underlying Wind Turbine Wake Control with Strouhal-Timed Actuation

A reduction in wake effects in large wind farms through wake-aware control has considerable potential to improve farm efficiency. This work examines the success of several emerging, empirically derived control methods that modify wind turbine wakes (i.e., the pulse method, helix method, and related methods) based on Strouhal numbers on the O(0.3). Drawing on previous work in the literature for jet and bluff-body flows, the analyses leverage the normal-mode representation of wake instabilities to characterize the large-scale wake meandering observed in actuated wakes. Idealized large-eddy simulations (LES) using an actuator-line representation of the turbine blades indicate that the n=0 and ±1 modes, which correspond to the pulse and helix forcing strategies, respectively, have faster initial growth rates than higher-order modes, suggesting these lower-order modes are more appropriate for wake control. Exciting these lower-order modes with periodic pitching of the blades produces increased modal growth, higher entrainment into the wake, and faster wake recovery. Modal energy gain and the entrainment rate both increase with streamwise distance from the rotor until the intermediate wake. This suggests that the wake meandering dynamics, which share close ties with the relatively well-characterized meandering dynamics in jet and bluff-body flows, are an essential component of the success of wind turbine wake control methods. A spatial linear stability analysis is also performed on the wake flows and yields insights on the modal evolution. In the context of the normal-mode representation of wake instabilities, these findings represent the first literature examining the characteristics of the wake meandering stemming from intentional Strouhal-timed wake actuation, and they help guide the ongoing work to understand the fluid-dynamic origins of the success of the pulse, helix, and related methods.

17 WIND ENERGY↗

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.↗

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.↗

A Unified Development of Basis Reduction Methods for Rotor Blade Analysis

The axial foreshortening effect plays a key role in rotor blade dynamics, but approximating it accurately in reduced basis models has long posed a difficult problem for analysts. Recently, though, several methods have been shown to be effective in obtaining accurate,reduced basis models for rotor blades. These methods are the axial elongation method,the mixed finite element method, and the nonlinear normal mode method. The main objective of this paper is to demonstrate the close relationships among these methods, which are seemingly disparate at first glance. First, the difficulties inherent in obtaining reduced basis models of rotor blades are illustrated by examining the modal reduction accuracy of several blade analysis formulations. It is shown that classical, displacement-based finite elements are ill-suited for rotor blade analysis because they can't accurately represent the axial strain in modal space, and that this problem may be solved by employing the axial force as a variable in the analysis. It is shown that the mixed finite element method is a convenient means for accomplishing this, and the derivation of a mixed finite element for rotor blade analysis is outlined. A shortcoming of the mixed finite element method is that is that it increases the number of variables in the analysis. It is demonstrated that this problem may be rectified by solving for the axial displacements in terms of the axial forces and the bending displacements. Effectively, this procedure constitutes a generalization of the widely used axial elongation method to blades of arbitrary topology. The procedure is developed first for a single element, and then extended to an arbitrary assemblage of elements of arbitrary type. Finally, it is shown that the generalized axial elongation method is essentially an approximate solution for an invariant manifold that can be used as the basis for a nonlinear normal mode.

Ruzicka, Gene C.↗

Suhl instabilities in nanoscopic spheroids

We simulate the magnetization dynamics of a permalloy spheroid of nanoscopic size in zero external field, such that both dipolar and exchange interactions are important. Low excitation power is used to obtain the frequencies and mode patterns of many normal modes. At higher power, non-linear three and four mode couplings between magnons carrying orbital angular momentum are observed to give rise to Suhl instabilities. Suhl’s analysis of the selection rules governing the allowed processes is extended to initial states other than uniform precession. These rules are studied and confirmed by the simulations. Both down- and up-conversion are seen as well as three and four-mode processes. As a result, general trends are inferred for preferred instabilities among those that are allowed, although the thresholds for some instabilities appear to be very high.

3D nanomagnetism↗

A substructure technique for dynamics of flexible mechanical systems with contact-impact

A substructure synthesis method is proposed to account for contact-impact effects in flexible components of mechanical systems. Components that may come into contact is divided into substructures, on each of which local deformation modes are defined to described deformation fields of components. Constraint modes and fixed interface normal modes are used to account for elastic deformation within each substructure. A constraint addition-deletion technique is used to account for contact between impacting bodies. Lagrange multipliers associated with the constraints, which represent constraint forces, are used to determine separation of contacting nodes. Use of the method is illustrated on problems of longitudinal and transverse impact of bodies.

Wu, S.-C.↗

Diskoseismology: Probing accretion disks. I - Trapped adiabatic oscillations

The normal modes of acoustic oscillations within thin accretion disks which are terminated by an innermost stable orbit around a slowly rotating black hole or weakly magnetized compact neutron star are analyzed. The dominant relativistic effects which allow modes to be trapped within the inner region of the disk are approximated via a modified Newtonian potential. A general formalism is developed for investigating the adiabatic oscillations of arbitrary unperturbed disk models. The generic behavior is explored by way of an expansion of the Lagrangian displacement about the plane of symmetry and by assuming separable solutions with the same radial wavelength for the horizontal and vertical perturbations. The lowest eigenfrequencies and eigenfunctions of a particular set of radial and quadrupole modes which have minimum motion normal for the plane are obtained. These modes correspond to the standard dispersion relation of disk theory.

Nowak, Michael A.↗

A self-consistent linear-mode model of stellar convection

A normal-mode expansion of the linearized fluid equations in terms of small subset of spherical harmonics can provide a foundation for a physically motivated, self-consistent description of a solar-type convection zone. In the absence of dissipation, a second-order differential equation governs the radial dependence of the modes, so that interpretation of the effects on convection quantities of the normal-form 'potential well' is straightforward. The philosophy is quite different from the more recent work of Narasimha and Antia (1982): all envelopes presented here differ substantially from MLT envelopes, and therefore, from theirs, which are constructed to be consistent with MLT. The amplitude of all modes is set by a Kelvin-Helmholtz-('shear'-) instability argument unrelated to solar observations, with the result that the convection description may be considered to arise from 'first-hueristic-principles'. The thermodynamics modelled vaguely resemble the sun's, and more vigorously convective envelopes show some phenomena qualitatively like solar observations (e.g., atmospheric velocity spectra).

Macauslan, J.↗

Improved orthogonality check for measured modes

A method is proposed for performing an orthogonality check for normal modes derived from ground vibration testing. The method utilizes partitioned mass and stiffness matrices for a linear undamped representation of a structure. The normalization of the modes by the proposed method inherently includes the effects of significant displacements which were not measured; and the method may allow the use of fewer measurement points than would be necessary with the conventional method.

Berman, A.↗