Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “third harmonic generation”

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.

49 records · Page 3

Observations of Mercury's magnetic field

Magnetic field data obtained by Mariner 10 during the third and final encounter with the planet Mercury on 16 March 1975 were studied. A well developed bow shock and modest magnetosphere, previously observed at first encounter on 29 March 1974, were again observed. In addition, a much stronger magnetic field near closest approach, 400 gamma versus 98 gamma, was observed at an altitude of 327 km and approximately 70 deg north Mercurian latitude. Spherical harmonic analysis of the data provide an estimate of the centered planetary magnetic dipole of 4.7 x 10 to the 22nd power Gauss/cu cm with the axis tilted 12 deg to the rotation axis and in the same sense as Earth's. The interplanetary field was sufficiently different between first and third encounters that in addition to the very large field magnitude observed, it argues strongly against a complex induction process generating the observed planetary field. While a possibility exists that Mercury possesses a remanent field due to magnetization early in its formation, a present day active dynamo seems to be a more likely candidate for its origin.

Ness, N. F.↗

Improving Power Density of Free-Piston Stirling Engines

Analyses and experiments demonstrate the potential benefits of optimizing piston and displacer motion in a free-piston Stirling Engine. Isothermal analysis shows the theoretical limits of power density improvement due to ideal motion in ideal Stirling engines. More realistic models based on nodal analysis show that ideal piston and displacer waveforms are not optimal, often producing less power than engines that use sinusoidal piston and displacer motion. Constrained optimization using nodal analysis predicts that Stirling engine power density can be increased by as much as 58 percent using optimized higher harmonic piston and displacer motion. An experiment is conducted in which an engine designed for sinusoidal motion is forced to operate with both second and third harmonics, resulting in a piston power increase of as much as 14 percent. Analytical predictions are compared to experimental data and show close agreement with indirect thermodynamic power calculations, but poor agreement with direct electrical power measurements.

aerospace engineering↗

Third and higher harmonic plasma emission due to Raman scattering

The theory for third and higher harmonic plasma emission by the weak turbulence (or random phase) process L + T prime to T (where L denotes a Langmuir wave, and T and T prime denote transverse waves) is developed. Kinematic constraints on the characteristics and growth lengths of waves participating in the wave processes are identified. The cases of L waves produced either directly by a streaming instability, or by the decay L to L prime + S (S is an ion sound wave) of L waves generated by a streaming instability, are considered. Limits on the brightness temperature of the radiation are determined, and expressions for the growth rate and path-integrated wave temperatures are derived.

Cairns, Iver H.↗

Direct numerical simulations of a spatially developing plane wake

In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Maekawa, Hiroshi↗

Unsteady transonic airfoil computation using implicit Euler scheme on body-fixed grid

The unsteady Euler equations have been derived for the flow relative motion with respect to a frame of reference that is rigidly attached to the moving airfoil. The grid is generated once by an elliptic solver without a need for dynamic grid computation. An implicit factored finite-volume scheme has been developed and implemented through a fully vectorized computer program. Implicit second-order and explicit second and fourth-order dissipations are added to the scheme. The boundary conditions are explicitly satisfied. The scheme is applied to steady and unsteady transonic airfoil flows and the results are in good agreement with the experimental data. For forced harmonic airfoil motions, periodic solutions are achieved within the third cycle of oscillation.

Kandil, Osama A.↗

Optimum wavelengths for two color ranging

The range uncertainties associated with the refractive atmosphere can be mitigated by the technique of two color, or dual wavelength, ranging. The precision of the differential time of flight (DTOF) measurement depends on the atmospheric dispersion between the two wavelengths, the received pulsewidths and photoelectron counts, and on the amount of temporal averaging. In general, the transmitted wavelengths are not independently chosen but instead are generated via nonlinear optics techniques (harmonic crystals, Raman scattering, etc.) which also determine their relative pulsewidths. The mean received photoelectrons at each wavelength are calculated via the familiar radar link equation which contains several wavelength dependent parameters. By collecting the various wavelength dependent terms, one can define a wavelength figure of merit for a two color laser ranging system. In this paper, we apply the wavelength figure of merit to the case of an extremely clear atmosphere and draw several conclusions regarding the relative merits of fundamental-second harmonic, fundamental-third harmonic, second-third harmonic, and Raman two color systems. We find that, in spite of the larger dispersion between wavelengths, fundamental-third harmonic systems have the lowest figure of merit due to a combination of poor detector performance at the fundamental and poor atmospheric transmission at the third harmonic. The fundamental-second harmonic systems (approximately 700 nm and 350 nm) have the highest figure of merit, but second-third harmonic systems, using fundamental transmitters near 1000 nm, are a close second. Raman-shifted transmitters appear to offer no advantage over harmonic systems because of the relatively small wavelength separation that can be achieved in light gases such as hydrogen and the lack of good ultrashort pulse transmitters with an optimum fundamental wavelength near 400 nm.

Degnan, John J.↗

On the phase-locked polarization variations in Cygnus X-1

The paper presents data on the accumulated V-band orbital phase curves of polarization in Cygnus X-1 which disagree with the analysis of Simmons, Aspin, and Brown (1980). No curves which conform to the canonical model were observed, and a noncanonical structure is indicated by F tests on third, fourth, and fifth harmonics, which have a probability of 0.85 as a group. The fifth harmonic has a probability of 0.99; in contrast, the fundamental harmonic has a probability of only 0.6, which is below the detection threshold. Therefore, major conclusions of Simmons et al based on a nonsymmetrical canonical model assuming that fundamental harmonics are present, but harmonics with n greater than 2 are absent, are unacceptable. Analysis of the author's data by Simmons et al is criticized on the grounds that he did not take into account the cumulative aspects of the data and did not recognize the need to average out the random noise generated by Cygnus X-1 over many orbital cycles.

Kemp, J. C.↗

Enhancements to the SSME transfer function modeling code

This report details the results of a one year effort by Ohio University to apply the transfer function modeling and analysis tools developed under NASA Grant NAG8-167 (Irwin, 1992), (Bartholomew, 1992) to attempt the generation of Space Shuttle Main Engine High Pressure Turbopump transfer functions from time domain data. In addition, new enhancements to the transfer function modeling codes which enhance the code functionality are presented, along with some ideas for improved modeling methods and future work. Section 2 contains a review of the analytical background used to generate transfer functions with the SSME transfer function modeling software. Section 2.1 presents the 'ratio method' developed for obtaining models of systems that are subject to single unmeasured excitation sources and have two or more measured output signals. Since most of the models developed during the investigation use the Eigensystem Realization Algorithm (ERA) for model generation, Section 2.2 presents an introduction of ERA, and Section 2.3 describes how it can be used to model spectral quantities. Section 2.4 details the Residue Identification Algorithm (RID) including the use of Constrained Least Squares (CLS) and Total Least Squares (TLS). Most of this information can be found in the report (and is repeated for convenience). Section 3 chronicles the effort of applying the SSME transfer function modeling codes to the a51p394.dat and a51p1294.dat time data files to generate transfer functions from the unmeasured input to the 129.4 degree sensor output. Included are transfer function modeling attempts using five methods. The first method is a direct application of the SSME codes to the data files and the second method uses the underlying trends in the spectral density estimates to form transfer function models with less clustering of poles and zeros than the models obtained by the direct method. In the third approach, the time data is low pass filtered prior to the modeling process in an effort to filter out high frequency characteristics. The fourth method removes the presumed system excitation and its harmonics in order to investigate the effects of the excitation on the modeling process. The fifth method is an attempt to apply constrained RID to obtain better transfer functions through more accurate modeling over certain frequency ranges. Section 4 presents some new C main files which were created to round out the functionality of the existing SSME transfer function modeling code. It is now possible to go from time data to transfer function models using only the C codes; it is not necessary to rely on external software. The new C main files and instructions for their use are included. Section 5 presents current and future enhancements to the XPLOT graphics program which was delivered with the initial software. Several new features which have been added to the program are detailed in the first part of this section. The remainder of Section 5 then lists some possible features which may be added in the future. Section 6 contains the conclusion section of this report. Section 6.1 is an overview of the work including a summary and observations relating to finding transfer functions with the SSME code. Section 6.2 contains information relating to future work on the project.

Irwin, R. Dennis↗

Growth of single crystals of organic salts with large second-order optical nonlinearities by solution processes for devices

Data obtained from the electric field induced second harmonic generation (EFISH) and Kurtz Powder Methods will be provided to MSFC for further refinement of their method. A theoretical model for predicting the second-order nonlinearities of organic salts is being worked on. Another task is the synthesis of a number of salts with various counterions. Several salts with promising SHG activities and new salts will be tested for the presence of two crystalline forms. The materials will be recrystallized from dry and wet solvents and compared for SHG efficiency. Salts that have a high SHG efficiency and no tendency to form hydrates will be documented. The synthesis of these materials are included in this report. A third task involves method to aid in the growth of large, high quality single crystals by solution processes. These crystals will be characterized for their applicability in the fabrication of devices that will be incorporated into optical computers in future programs. Single crystals of optimum quality may be obtained by crystal growth in low-gravity. The final task is the design of a temperature lowering single crystal growth apparatus for ground based work. At least one prototype will be built.

Leslie, Thomas M.↗

Elastic Constants of Solids and Fluids with Initial Pressure via a Unified Approach Based on Equations-of-State

The second and third-order Brugger elastic constants are obtained for liquids and ideal gases having an initial hydrostatic pressure p(sub 1). For liquids the second-order elastic constants are C(sub 11) = A + p(sub 1), C(sub 12) = A -- p(sub 1), and the third-order constants are C(sub 111) = --(B + 5A + 3p(sub 1)), C(sub 112) = --(B + A -- p(sub 1)), and C(sub 123) = A -- B -- p1, where A and B are the Beyer expansion coefficients in the liquid equation of state. For ideal gases the second order constants are C(sub 11) = p(sub 1)gamma + p9sub 1), C(sub 12) = p(sub 1)gamma -- p(sub 1), and the third-order constants are C(sub 111) = p(sub 1)(gamma(2) + 4gamma + 3), C(sub 112) = --p(sub 1)(gamma(2) -- 1), and C(sub 123) = --p(sub 1) (gamma(2) -- 2gamma + 1), where gamma is the ratio of specific heats. The inequality of C(sub 11) and C(sub 12) results in a nonzero shear constant C(sub 44) = (1/2)(C(sub 11) C(sub 12)) = p(sub 1) for both liquids and gases. For water at standard temperature and pressure the ratio of terms p1/A contributing to the second-order constants is approximately 4.3 x 10(-5). For atmospheric gases the ratio of corresponding terms is approximately 0.7. Analytical expressions that include initial stresses are derived for the material 'nonlinearity parameters' associated with harmonic generation and acoustoelasticity for fluids and solids of arbitrary crystal symmetry. The expressions are used to validate the relationships for the elastic constants of fluids.

Cantrell, John H.↗

Gravitational model improvement at the Goddard Space Flight Center

Major new computations of terrestrial gravitational field models were performed by the Geodynamics Branch of Goddard Space Flight Center (GSFC). This development has incorporated the present state of the art results in satellite geodesy and have relied upon a more consistent set of reference constants than was heretofore utilized in GSFC's GEM models. The solutions are complete in spherical harmonic coefficients out to degree 50 for the gravity field parameters. These models include adjustment for a subset of 66 ocean tidal coefficients for the long wavelength components of 12 major ocean tides. This tidal adjustment was made in the presence of 550 other fixed ocean tidal terms representing 32 major and minor ocean tides and the Wahr frequency dependent solid earth tidal model. In addition 5-day averaged values for Earth rotation and polar motion were derived for the time period of 1980 onward. Two types of models were computed. These are satellite only models relying exclusively on tracking data and combination models which have incorporated satellite altimetry and surface gravity data. The satellite observational data base consists of over 1100 orbital arcs of data on 31 satellites. A large percentage of these observations were provided by third generation laser stations (less than 5 cm). A calibration of the model accuracy of the GEM-T2 satellite only solution indicated that it was a significant improvement over previous models based solely upon tracking data. The rms geoid error for this field is 110 cm to degree and order 36. This is a major advancement over GEM-T1 whose errors were estimated to be 160 cm. An error propagation using the covariances of the GEM-T2 model for the TOPEX radial orbit component indicates that the rms radial errors are expected to be 12 cm. The combination solution, PGS-3337, is a preliminary effort leading to the development of GEM-T3. PGS-3337 has incorporated global sets of surface gravity data and the Seasat altimetry to produce a model complete to (50,50). A solution for the dynamic ocean topography to degree and order 10 was included as part of this adjustment.

Marsh, J. G.↗

Avalanche-diode oscillator circuit with tuning at multiple frequencies

Detailed theoretical analysis of three different modes or types of high efficiency oscillation in a PIN diode are presented. For the TRAPATT mode in a PIN diode, it is shown that a traveling avalanche zone is not necessary to generate a dense trapped plasma. An economical computer program for TRAPATT oscillations in a PIN diode is described. Typical results of diode power, dc-to-RF conversion efficiency, and required circuit impedances are presented for several different current waveforms. A semianalytical solution for a second type of high efficiency mode in a PIN diode is derived assuming a rectangular current waveform. A quasi-static approximation is employed to derive a semianalytical solution for the voltage across a PIN diode in a third mode, where avalanching occurs during a major portion of a half cycle. Calculations for this mode indicate that the power increases proportionally to the magnitude of the drive current with a small decrease in efficiency relative to the ordinary TRAPATT mode. An analytical solution is also given for a PIN diode, where it is assumed that the ionization coefficient is a step function. It is shown that the step-ionization approximation permits one to draw possible patterns of avalanche region in the depletion layer as a function of time. A rule governing admissible patterns is derived and an example solution given for one admissible pattern. Preliminary experimental results on the high-efficiency oscillations are presented and discussed. Two different experimental circuits, which used channel-dropping filters to provide independent harmonic tuning, are described. Simpler circuits used to produce high-efficiency oscillations are discussed. Results of experiments using inexpensive Fairchild FD300 diodes are given.

Parker, D.↗

Modeling Wave Driven Non-linear Flow Oscillations: The Terrestrial QBO and a Solar Analog

The Quasi Biennial Oscillation (QBO) of the zonal circulation observed in the terrestrial atmosphere at low latitudes is driven by wave mean flow interaction as was demonstrated first by Lindzen and Holton (1968), shown in a laboratory experiment by Plumb and McEwan (1978), and modeled by others (e.g., Plumb, Dunkerton). Although influenced by the seasonal cycle of solar forcing, the QBO, in principle, represents a nonlinear flow oscillation that can be maintained by a steady source of upward propagating waves. The wave driven non-linearity is of third or odd order in the flow velocity, which regenerates the fundamental harmonic itself to keep the oscillation going - the fluid dynamical analog of the displacement mechanism in the mechanical clock. Applying Hines' Doppler Spread Parameterization (DSP) for gravity waves (GW), we discuss with a global-scale spectral model numerical experiments that elucidate some properties of the QBO and its possible effects on the climatology of the atmosphere. Depending on the period of the QBO, wave filtering can cause interaction with the seasonal variations to produce pronounced oscillations with beat periods around 10 years. Since the seasonal cycle and its variability influence the period of the QBO, it may also be a potent conduit of solar activity variations to lower altitudes. Analogous to the terrestrial QBO, we propose that a flow oscillation may account for the 22-year periodicity of the solar magnetic cycle, potentially answering Dicke (1978) who asked, "Is there a chronometer hidden deep inside the Sun?" The oscillation would occur below the convection region, where gravity waves can propagate. Employing a simplified, analytic model, Hines' DSP is applied to estimate the flow oscillation. Depending on the adopted horizontal wavelengths of GW's, wave amplitudes less than 10 m/s can be made to produce oscillating zonal flows of about 20 m/s that should be large enough to generate a significant oscillation in the magnetic field. For the large length scales of the Sun, the flow cycle period tends to be very long. The period, however, can be made to be 22 years, provided the buoyancy frequency (stability) is sufficiently small, thus placing the proposed flow near the base of the convection zone where a dynamo is now believed to operate.

Mayr, Hans G.↗