Nonlinear Systems- The Parameter Analysis and Design.
Book on nonlinear systems analysis and design covering parameter mapping, symmetrical, transient and forced oscillations and stability analysis
Siljak, D. D.↗
SEARCH · Engineering Papers
Results for “nonlinear oscillation”
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.
A macroscopic plasma Lagrangian and its application to wave interactions and resonances
The derivation of a macroscopic plasma Lagrangian is considered, along with its application to the description of nonlinear three-wave interaction in a homogeneous plasma and linear resonance oscillations in a inhomogeneous plasma. One approach to obtain the Lagrangian is via the inverse problem of the calculus of variations for arbitrary first and second order quasilinear partial differential systems. Necessary and sufficient conditions for the given equations to be Euler-Lagrange equations of a Lagrangian are obtained. These conditions are then used to determine the transformations that convert some classes of non-Euler-Lagrange equations to Euler-Lagrange equation form. The Lagrangians for a linear resistive transmission line and a linear warm collisional plasma are derived as examples. Using energy considerations, the correct macroscopic plasma Lagrangian is shown to differ from the velocity-integrated low Lagrangian by a macroscopic potential energy that equals twice the particle thermal kinetic energy plus the energy lost by heat conduction.
Peng, Y. K. M.↗
A numerical study of the stability of radiative shocks
Attention is given to the oscillatory instability of optically thin radiative shocks in time-dependent numerical calculations of accretion flows onto degenerate dwarfs. The present nonlinear calculations yield good quantitative agreement with the linear results obtained for oscillation frequencies, damping rates, and critical alpha-values. The fundamental mode and the first overtone in the shock radius and luminosity variations can be clearly identified, and evidence is sometimes seen for the second overtone. Time-dependent calculations are also performed which include additional physics relevant to degenerate dwarf accretion, such as electron thermal conduction, unequal electron and ion temperatures, Compton cooling, and relativistic corrections to the bremsstrahlung cooling law. All oscillatory modes are found to be damped, and hence stable, in the case of a 1-solar mass white dwarf accreting in spherical symmetry.
Imamura, J. N.↗
HEMT Frequency Doubler with Output at 300 GHz
An active frequency doubler in the form of an InP-based monolithic microwave integrated circuit (MMIC) containing a high-electron-mobility transistor (HEMT) has been demonstrated in operation at output frequencies in the vicinity of 300 GHz. This is the highestfrequency HEMT doubler reported to date, the next-highest-frequency active HEMT doubler having been previously reported to operate at 180 GHz. While the output power of this frequency doubler is less than that of a typical Schottky diode, this frequency doubler is considered an intermediate product of a continuing effort to realize the potential of active HEMT frequency doublers to operate with conversion efficiencies greater than those of passive diode frequency doublers. An additional incentive for developing active HEMT frequency doublers lies in the fact that they can be integrated with amplifiers, oscillators, and other circuitry on MMIC chips. The circuitry of the doubler MMIC (see Figure 1) features grounded coplanar waveguides. Air bridges and vias are used to make contact with the ground plane. The HEMT is biased for Class-A operation (in which current is conducted throughout each cycle of oscillation), which would ordinarily be better suited to linear amplification than to frequency doubling. Ordinarily, class-B operation (in which current is conducted during about half of each cycle of oscillation) would be more suitable for frequency doubling because of the essential nonlinearity of partial-cycle conduction. The reason for the unusual choice of class A was that computational simulations had shown that in this case, the efficiency in class B would be less than in class A. The input matching circuit of this doubler includes transmission lines that afford a good impedance match at the fundamental frequency, plus an open stub to prevent leakage of the second harmonic through the input terminals. The output circuit was designed to suppress the fundamental while providing a good match for the second harmonic.
Samoska, Lorene↗
Highly nondegenerate all-resonant optical parametric oscillator
We show that a nondegenerate multi-frequency parametric oscillator has different properties compared with an usual three-wave parametric oscillator. We consider, as an example, a scheme of a resonant cw monolithic microwave-optical parametric oscillator based on high-Q whispering gallery modes excited in a nonlinear dielectric cavity. In this paper we recall main properties of usual three-mode OPO, discuss four-mode OPO, study conditions of generation of a frequency comb in all-resonant OPO with two-frequency optical pumping, and discuss stability of a cavity with moving mirrors.
whispering↗
Dynamic Instability of Solidification Fronts and Formation of Transverse Layered Structures
A plane front of crystallization in a unidirectional process of solidification of a binary melt or a solution is shown to be unstable with respect to virtual perturbations of the front velocity within some region of physical and operating parameters. This dynamic instability has nothing in common with morphological instability of the front also possible under certain conditions. It gives rise to the origination of a self-oscillating regime of solidification. The oscillations of the front velocity are accompanied by oscillations of other relevant variables, such as the melt concentration and temperature at the front. What is most important is that this provides for a natural explanation of phenomena of the formation of ordered striations commonly observed when producing alloys and withdrawing crystals. A theory of nonlinear instability of the solidification front is developed which also allows for properties of both self-oscillating regimes of solidification and transverse layered structures it allows to be described in an explicit form, without introducing any additional hypotheses and empirical parameters. Apart from other things, the properties are strongly dependent on convective transport originated in consequence of a difference between the densities of the liquid melt and solid alloy. The theory results in a remarkably good agreement with experimental evidence on crystal withdrawal obtained in normal gravity and under space flight conditions. The problem of interplay between modes of the dynamic and morphological types of instability is also discussed. There are parametric regions in which either only one of these types occurs or both of them arise simultaneously.
Buyevich, Yu A.↗
Nonlinear Analysis of a Two-Parachute System Undergoing Pendulum Motion
Motion resembling that of a pendulum undergoing large-amplitude limit cycle oscillation was observed during a series of flight tests of an unoccupied Orion Capsule Parachute Assembly System (CPAS) comprised of two parachutes and a capsule payload. Large excursions away from vertical by the capsule could cause it to strike the ground or ocean at a large angle with respect to vertical, or at a large horizontal speed. These conditions are undesirable because they would endanger the occupants of the capsule in an actual mission. A simplified planar dynamics model in conjunction with a nonlinear normal force coefficient vs. angle of attack model serves as the basis of an analytical investigation of the fundamental dynamics of this pendulum motion. Output error methodology from system identification theory was used to identify the parameters of the nonlinear aerodynamics model. The identified model yielded excellent comparison with portions of flight test data where the pendulum motion occurred. Due to the inherent nonlinear nature of the pendulum motion limit cycle, traditional nonlinear analysis techniques were applied to gain further insight into the system. Lyapunov’s direct method provided mathematical proof in the absolute stability of the pendulum mode. Describing Function method was used to predict the amplitude and frequency of the limit cycle oscillation. Finally, phase plane analysis allowed easy visualization on the size and shape of the limit cycle with respect to variations in key aerodynamic parameters.
Pei, Jing↗
Characterizing Observed Limit Cycles in the Cassini Main Engine Guidance Control System
The Cassini spacecraft dynamics-related telemetry during long Main Engine (ME) burns has indicated the presence of stable limit cycles between 0.03-0.04 Hz frequencies. These stable limit cycles cause the spacecraft to possess non-zero oscillating rates for extended periods of time. This indicates that the linear ME guidance control system does not model the complete dynamics of the spacecraft. In this study, we propose that the observed limit cycles in the spacecraft dynamics telemetry appear from a stable interaction between the unmodeled nonlinear elements in the ME guidance control system. Many nonlinearities in the control system emerge from translating the linear engine gimbal actuator (EGA) motion into a spacecraft rotation. One such nonlinearity comes from the gear backlash in the EGA system, which is the focus of this paper. The limit cycle characteristics and behavior can be predicted by modeling this gear backlash nonlinear element via a describing function and studying the interaction of this describing function with the overall dynamics of the spacecraft. The linear ME guidance controller and gear backlash nonlinearity are modeled analytically. The frequency, magnitude, and nature of the limit cycle are obtained from the frequency response of the ME guidance controller and nonlinear element. In addition, the ME guidance controller along with the nonlinearity is simulated. The simulation response contains a limit cycle with similar characterstics as predicted analytically: 0.03-0.04 Hz frequency and stable, sustained oscillations. The analytical and simulated limit cycle responses are compared to the flight telemetry for long burns such as the Saturn Orbit Insertion and Main Engine Orbit Trim Maneuvers. The analytical and simulated limit cycle characteristics compare well with the actual observed limit cycles in the flight telemetry. Both have frequencies between 0.03-0.04 Hz and stable oscillations. This work shows that the stable limit cycles occur due to the interaction between the unmodeled nonlinear elements and linear ME guidance controller.
stable limit cycle↗
A general numerical analysis program for the superconducting quasiparticle mixer
A user-oriented computer program SISCAP (SIS Computer Analysis Program) for analyzing SIS mixers is described. The program allows arbitrary impedance terminations to be specified at all LO harmonics and sideband frequencies. It is therefore able to treat a much more general class of SIS mixers than the widely used three-frequency analysis, for which the harmonics are assumed to be short-circuited. An additional program, GETCHI, provides the necessary input data to program SISCAP. The SISCAP program performs a nonlinear analysis to determine the SIS junction voltage waveform produced by the local oscillator. The quantum theory of mixing is used in its most general form, treating the large signal properties of the mixer in the time domain. A small signal linear analysis is then used to find the conversion loss and port impedances. The noise analysis includes thermal noise from the termination resistances and shot noise from the periodic LO current. Quantum noise is not considered. Many aspects of the program have been adequately verified and found accurate.
Hicks, R. G.↗
The Stability of Radiatively Cooling Jets. 2: Nonlinear Evolution
We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.
Stone, James M.↗
The collective excitation of g-modes in the sun
Oscillations of the solar interior (linear g-modes) may be strongly driven by the collective influence of all the modes upon the nuclear reactions in the core. This heretofore neglected effect could couple the modes, reduce the effective amplitudes near the center, and spatially concentrate most of the oscillation energy into just a portion of the radiative interior. If operating at sufficient strength, this can reverse the conventional conclusion, drawn from single mode calculations, that almost all solar g-modes are damped. Furthermore, it would put the theory into rough harmony with three otherwise troubling observations: (1) the 'low' neutrino flux measured by Davis (1978), (2) the high correspondence found by Wolff (1976) between recurrence periods in solar activity and the rotational beat periods of g-modes, and (3) the fluctuations in the sun's diameter which imply g-mode activity at high angular harmonics (Hill and Caudell 1979). A nonlinear expression is derived for the local rate of work done on an array of oscillation modes by the nuclear reactions. Three additional tests of the model are suggested.
Wolff, C. L.↗
Can El Nino Amplify the Solar Forcing of Climate?
ENSO (El Nino and the Southern Oscillation) is considered as a stochastic driver that excites the atmospheric anomaly states, such as Pacific North American pattern. This can make the 11 year solar activity forcing feasible to climate through stochastic resonance -- a phenomenon that amplifies a weak input to a nonlinear bistable system by the assistance of noise.
ENSO El Nino Oscillation climate stochastic resona↗
Higher-Order Spectral Analysis of a Nonlinear Pitch and Plunge Apparatus
Simulated aeroelastic responses of a nonlinear pitch and plunge apparatus are analyzed using various statistical signal processing techniques including higher-order spectral methods. A MATLAB version of the Nonlinear Aeroelastic Testbed Apparatus (NATA) at the Texas A&M University is used to generate various aeroelastic response data including limit cycle oscillations (LCO). Traditional and higher-order spectral (HOS) methods are applied to the simulated aeroelastic responses. Higher-order spectral methods are used to identify critical signatures that indicate the transition from linear to nonlinear (LCO) aeroelastic behavior.
Silva, Walter A.↗
Intra-Seasonal Monthly Oscillations in Stratospheric NCEP Data and Model Results
Intra-seasonal oscillations (ISO) are observed in the zonal-mean of mesospheric wind and temperature measurements-and the numerical spectral model (NSM) generates such oscillations. Relatively large temperature ISO are evident also in stratospheric CPC (NCEP) data at high latitudes, where the NSM produces amplitudes around 3 K at 30 km. Analyzing the NCEP data for the years 1996-2006, we find in Fourier spectra signatures of oscillations with periods between 1.7 and 3 months. With statistical confidence levels exceeding 70%, the spectral features are induced by nonlinear interactions involving the annual and semi-annual variations. The synthesized data show for the 10-year average that the temperature ISO peak in winter, having amplitudes close to 4 K. The synthesized complete spectrum for periods around 2 months produces oscillations, varying from year to year, which can reach peak amplitudes of 15 and 5 K respectively at northern and southern polar latitudes.
Mayr, H. G.↗
A 200 GHz tripler using a single barrier varactor
The GaAs Schottky varactor diode is the nonlinear device most commonly used for submillimeter wave harmonic generation. Output power adequate to serve as a local oscillator source for SIS tunnel junctions has been demonstrated with whisker-contacted GaAs Schottky varactor multipliers in waveguide mounts up to about 800 GHz. In this paper, we present results for a tripler to 200 GHz using a new multiplier device, the single barrier varactor (SBV). This new varactor has potential advantages such as stronger nonlinearities or special symmetry, which make it attractive for submillimeter wave frequency multiplication. The performance of a tripler using a SBV over a output frequency range from 186 to 207 GHz has been measured in a crossed waveguide mount. The theoretical performance of the device has been calculated using large signal analysis. A comparison of theoretical and measured results and a discussion of various losses in the mount and the varactor have also been presented.
Choudhury, Debabani↗
A 200 GHz tripler using single barrier varactor
The GaAs Schottky varactor diode is the nonlinear device most commonly used for submillimeter wave harmonic generation. Output power adequate to serve as a local oscillator source for SIS tunnel junctions has been demonstrated with whisker-contacted GaAs Schottky varactor multipliers in waveguide mounts up to about 800 GHz. In this paper, we present results for a tripler to 200 GHz using a new multiplier device, the single barrier varactor (SBV). This new varactor has a potential advantages such as stronger nonlinearities or special symmetry, which make it attractive for submillimeter wave frequency multiplication. The performance of a tripler using a SBV over a output frequency range from 186 to 207 GHz has been measured in a crossed waveguide mount. The theoretical performance of the device has been calculated using large signal analysis. A comparison of theoretical and measured results and a discussion of various losses in the mount and the varactor have also been presented.
Choudhury, Debabani↗
Observations and 3D MHD Modeling of Small-scale Solar Prominence Oscillations
Small scaleoscillations were recently observed in prominences with high resolution Hinode/SOT and IRIS SJI in unprecedented detailand identified as fast magnetosonic waves. The high-resolution observations provide evidence of flows of cool material, fluid instabilities such as Kelvin-Helmholtz and Rayleigh-Taylor, and the spectroscopic information provides evidence of Doppler shifts, heating, and cooling of the prominence material. Recently nonlinear fast magnetosonic waves were observed withperiods in the range of 5–11 minutes and wavelengths in the plane of the sky (POS) of ~2000kmand flows, and the flows in narrow threads withPOS speed ∼16–46 km/s. Thenonlinear fast magnetosonic waves and flows were modeled using 2.5D MHD (Ofman & Kucera 2020). Here, we extend the studyof the nonlinear fast magnetosonic waves using more realistic 3D MHD model. The model includes the 3D structure of the density and the magnetic fieldas well as mode coupling anddemonstrate the propagating density compressions associated with the nonlinear fast magnetosonicwaves, as well as the guided modesin the pillar. The 3DMHD modelingimprove understanding ofthe POSeffects in limb prominence oscillations observations.The results of the modeling are useful for coronal seismology of the prominence structure.
Leon Ofman↗
A comparison of the bounded derivative and the normal-mode initialization methods using real data
Application of the bounded-derivative and normal-mode methods to a simple linear barotropic model at a typical middle latitude shows that the two methods lead to identical constraints up to a certain degree of approximation. Beyond this accuracy the two methods may differ from each other. When applied to a global nonlinear barotropic model using real data, again the two methods lead to similar balanced initial states. The gravity oscillations in the unbalanced height field, which have amplitudes of up to 60 m with a dominant periodicity of about 5 to 6 h, are practically eliminated by both initialization methods. The rotational wind component is smooth even for the unbalanced initial state. The small-scale spatial features of the irrotational wind component are drastically reduced by initialization. Both the nonlinear normal-mode and the bounded-derivative initialization methods yield similar divergence fields centered around the areas of highest orography. The comparison shows that there is no significant loss of information in the mass and momentum fields, despite the fact that the bounded-derivative method employs only the original, rotational wind component to construct a balanced initial state compared to the normal-mode method, which, in addition, makes use of the unbalanced divergent wind and height fields.
Semazzi, F. H. M.↗