Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Degenerate parametric down-conversion facilitated by exciton-plasmon polariton states in a nonlinear plasmonic cavity

Abstract We study the effect of degenerate parametric down-conversion (DPDC) in an ensemble of two-level quantum emitters (QEs) coupled via near-field interactions to a single surface plasmon (SP) mode of a nonlinear plasmonic cavity. For this purpose, we develop a quantum driven-dissipative model capturing non-equilibrium dynamics of the system in which incoherently pumped QEs have transition frequency tuned near the second-harmonic response of the SPs. Considering the strong coupling regime, i.e. the SP-QE interaction rate exceeds system dissipation rates, we find a critical SP-QE coupling attributed to the phase transition between normal and lasing steady states. Examining fluctuations above the system’s steady states, we predict new elementary excitations, namely, the exciton-plasmon polaritons formed by the two-SP quanta and single-exciton states of QEs. The contribution of two-SP quanta results in the linear scaling of the SP-QE interaction rate with the number of QEs,  o , as opposed to the  o -scaling known for the Dicke and Tavis–Cummings models. We further examine how SP-QE interaction scaling affects the polariton dispersions and power spectra in the vicinity of the critical coupling. For this purpose, we compare the calculation results assuming a finite ensemble of QEs and the model thermodynamic limit. The calculated power spectra predict an interplay of coherent photon emission by QEs near the second-harmonic frequency and correlated photon-pair emission at the fundamental frequency by the SPs (i.e. the photonic DPDC effect).

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Quantitative Nonlinear Optical Polarimetry with High Spatial Resolution

Nonlinear optical microscopy such as in the optical second-harmonic generation (SHG) modality has become a popular tool today for probing materials in the physical and biological sciences. While imaging and spectroscopy are widely used in the microscopy mode, nonlinear polarimetry, which can shed light on materials’ symmetry and microstructure, is relatively underdeveloped. This is partly because quantitative analytical modeling of the optical SHG response for anisotropic crystals and films largely assumes low-numerical aperture (NA) focusing of light, where the plane-wave approximation is sufficient. Tight focusing provides unique benefits in revealing out-of-plane polarization responses, which cannot be detected by near-plane-wave illumination at normal incidence. Here, we outline a method for quantitatively analyzing SHG polarimetry measurements obtained under high-NA focusing within a microscope geometry. Experiments and simulations of a variety of standard samples, from single crystals to thin films, are in good agreement, including measured and simulated spatial SHG maps of ferroelectric domains. A solution to the inverse problem is demonstrated, where the spatial distribution of an SHG tensor with unknown tensor coefficient magnitudes is determined by experimentally measured polarimetry. The ability to extract the out-of-plane component of the nonlinear polarization in normal incidence is demonstrated, which can be valuable for high-resolution polarimetry of 2D materials, thin films, heterostructures, and uniaxial crystals with a strong out-of-plane response.

36 MATERIALS SCIENCE↗
NASA NTRS2000

The Investigation of the Effects of Gravity on Single Bubble Sonoluminescence

In single bubble following it's rapid collapse each cycle of oscillation of an ultrasonic field. Since widely varying length and time scales affect the bubble dynamics and optical emission processes, it is difficult to anticipate the importance of the effects of gravity present for observations on earth. Our bubble is driven in an acoustically resonating cavity at it's first harmonic mode. The acoustical radiation pressure (Bjerknes force) will then keep it suspended in the center near the pressure antinode. When driven in a region where the diffusive processes balance the bubble it acts in a nonlinear but regular way, emitting a short (approx. 200ps) burst of light each acoustic cycle. Balancing the Bjerknes force with buoyancy, as in, we can see that the bubble should be displaced from the velocity node approximately 20m at normal gravity. Therefore, water flows past the bubble at the time of collapse. Gravitation also changes the ambient pressure at the bubble's location, as Delta.P = rho.g.h this gives a change of approximately -0.5% in our experiment when going from 1.8g to 0g. Studies of ambient pressure changes were also done in order to assess these effects. Inside a pressure sealed chamber a spherical glass cell is filled with distilled water which has been degassed to 120mmHg. A bubble is then trapped in the center and driven by a piezoelectric transducer at 32.2kHz attached to the side of the cell. An optical system is then set up to take strobbed video images along and light emission data simultaneously. Temperature, pressure, drive voltage, and listener voltage are also monitored. PMT output in Volts The radii of the bubbles for both experiment s are fit using the Rayleigh-Plesset equation and the acoustic drive amplitude and the ambient bubble radius are found. There is little change in the acoustic drive amplitude as we expect, since we are not varying the drive voltage. However. the ambient bubble radius goes up considerably. These changes (increased light output, increased maximum bubble radius, and increased ambient bubble radius) are also observed when the ambient pressure is varied in the laboratory by an amount similar to that due to gravitation. The changes in the ambient bubble radius and light output with a change in ambient pressure are predicted by the "dissociation hypothesis" and have been observed by other groups in the laboratory. It seems clear that buoyancy's effect on light output and bubble radius, are at best on the same order as the effects of ambient pressure.

Dzikowicz, Ben↗
DOE OSTI2023

Isotope effects on energy transport in the core of ASDEX-Upgrade tokamak plasmas: Turbulence measurements and model validation

Design and operation of future tokamak fusion reactors using a deuterium–tritium 50:50 mix requires a solid understanding of how energy confinement properties change with ion mass. This study looks at how turbulence and energy transport change in L-mode plasmas in the ASDEX Upgrade tokamak when changing ion species between hydrogen and deuterium. For this purpose, both experimental turbulence measurements and modeling are employed. Local measurements of ion-scale (with wavevector of fluctuations perpendicular to the B-field k⊥< 2 cm−1, k⊥ρs< 0.2, where ρs is the ion sound Larmor radius using the deuterium ion mass) electron temperature fluctuations have been performed in the outer core (normalized toroidal flux ρTor=0.65−0.8) using a multi-channel correlation electron cyclotron emission diagnostic. Lower root mean square perpendicular fluctuation amplitudes and radial correlation lengths have been measured in hydrogen vs deuterium. Measurements of the cross-phase angle between a normal-incidence reflectometer and an ECE signal were made to infer the cross-phase angle between density and temperature fluctuations. The magnitude of the cross-phase angle was found larger (more out-of-phase) in hydrogen than in deuterium. TRANSP power balance simulations show a larger ion heat flux in hydrogen where the electron-ion heat exchange term is found to play an important role. These experimental observations were used as the basis of a validation study of both quasilinear gyrofluid trapped gyro-Landau fluid-SAT2 and nonlinear gyrokinetic GENE codes. Linear solvers indicate that, at long wavelengths (k⊥ρs<1), energy transport in the deuterium discharge is dominated by a mixed ion-temperature-gradient (ITG) and trapped-electron mode turbulence while in hydrogen transport is exclusively and more strongly driven by ITG turbulence. The Ricci validation metric has been used to quantify the agreement between experiments and simulations taking into account both experimental and simulation uncertainties as well as four different observables across different levels of the primacy hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Flow field and acoustic properties of a Mach number 0.9 jet at a low Reynolds number

An experimental study of the flow field and acoustic properties of a low Reynolds number M = 0.9 jet has been performed in a low pressure anechoic test chamber. The mean flow field was surveyed with a conventional pitot pressure probe and flow fluctuations were detected with a normal hot wire probe. Also, condenser microphone measurements were made in the acoustic field. The flow fluctuations within the jet were found to be dominated initially by a relatively discrete, large-scale, wave-like instability centered around a Strouhal number of 0.44. The axial wavelength of this instability was determined to be 1.45 jet diameters and its azimuthal character includes the n = 0 and n = plus or minus 1 modes. The growth of this instability coupled with its nonlinear breakdown are major contributors to the termination of the potential core region of the jet. The acoustic field of the jet, in contrast to the flow field, has a broad frequency spectrum with a peak amplitude near a Strouhal number of St = 0.2. The results indicate that a nonlinear mechanism involving the large-scale flow instability is responsible for a dominant portion of the noise generated from this jet.

Stromberg, J. L.↗
NASA NTRS1988

Oscillations of drops in zero gravity with weak viscous effects

Nonlinear oscillations and other motions of large axially symmetric liquid drops in zero gravity are studied numerically by a boundary-integral method. The effect of small viscosity is included in the computations by retaining first-order viscous terms in the normal stress boundary condition. This is accomplished by making use of a partial solution of the boundary-layer equations which describe the weak vortical surface layer. Small viscosity is found to have a relatively large effect on resonant mode coupling phenomena.

Lundgren, T. S.↗
DOE OSTI2022

Theory of disordered superconductors with applications to nonlinear current response

I present a review of the theory and basic equations for charge transport in superconducting alloys starting from the Keldysh formulation of the quasiclassical transport equations developed by Eilenberger, Larkin and Ovchinnikov, and Eliashberg. This formulation is the natural extension of Landau’s theory of normal Fermi liquids to the superconducting state of strongly correlated metals. For dirty metals the transport equations reduce to equations for charge diffusion, with the current response given by the Drude conductivity at low temperatures. The extension of the diffusion equation for the charge and current response of a strongly disordered normal metal to the superconducting state yields Usadel’s equations for the nonequilibrium quasiclassical Keldysh propagator. The conditions for the applicability of the Usadel equations are discussed, the pair-breaking effect of disorder on the current response, including the nonlinear current response to an electromagnetic field in the dirty limit, τ &NestedLessLess; ℏ/Δ, are reported. The same nonlinearity is shown to lead to source currents for photon generation and nonlinear Kerr rotation driven by the nonlinear response to excitation of the superconductor by a multi-mode electromagnetic field. The potential relevance of the nonlinear source currents to superconducting radio-frequency cavities as detectors of axion-like dark matter candidates is briefly discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1981

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗
NASA NTRS2013

Direct Numerical Simulation of Transition in a Swept-Wing Boundary Layer

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing derived from a two-dimensional, partial-differential-equation based eigenvalue computation; and the mode selected for forcing corresponds to the most amplified secondary instability mode which, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. Both the growth of the secondary instability wave and the resulting onset of laminar-turbulent transition are captured within the DNS computations. The growth of the secondary instability wave in the DNS solution compares well with linear secondary instability theory when the amplitude is small; the linear growth is followed by a region of reduced growth resulting from nonlinear effects before an explosive onset of laminar breakdown to turbulence. The peak fluctuations are concentrated near the boundary layer edge during the initial stage of transition, but rapidly propagates towards the surface during the process of laminar breakdown. Both time-averaged statistics and flow visualization based on the DNS reveal a sawtooth transition pattern that is analogous to previously documented surface flow visualizations of transition due to stationary crossflow instability. The memory of the stationary crossflow vortex is found to persist through the transition zone and well beyond the location of the maximum skin friction.

Duan, Lian↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2023

Unveiling the Significance of Correlations in K-Space and Configuration Space for Drift Wave Turbulence in Tokamaks

Turbulence and transport phenomena play a crucial role in the confinement and stability of tokamak plasmas. Turbulent fluctuations in certain physical quantities, such as density or temperature fluctuations, can have a wide range of spatial scales, and understanding their correlation length is important for predicting and controlling the behavior of the plasma. The correlation length in the radial direction is identified as the critical length in real space. The dynamics in real space are of significant interest because transport in configuration space is primarily focused on them. When investigating transport caused by the E×B drift, the correlation length in real space represents the size of E×B whirls. It was numerically discovered that in drift wave turbulence, this length is inversely proportional to the normalized mode number of the fastest growing mode relative to the drift frequency. Considerable time was required before a proper analytical derivation of this condition was accomplished. Therefore, a connection has been established between phenomena occurring in real space and those occurring in k-space. Although accompanied by a turbulent spectrum in k-space with a substantial width, transport in real space is uniquely determined by the correlation length, allowing for accurate transport calculations through the dynamics of a single mode. Naturally, the dynamics are subject to nonlinear effects, with resonance broadening in frequency being the most significant nonlinear effect. Thus, mode number space is once again involved. Resonance broadening leads to the detuning of waves from particles, permitting a fluid treatment. It should be emphasized that the consideration here involves the total electric field, including the induction part, which becomes particularly important at higher beta plasmas.

Weiland, Jan (ORCID:0000000253496025)↗
NASA NTRS1997

Analysis of Hydrodynamic (Landau) Instability in Liquid-Propellant Combustion at Normal and Reduced Gravity

The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.

Margolis, Stephen B.↗
NASA NTRS1989

A performance comparison of integration algorithms in simulating flexible structures

Asymptotic formulas for the characteristic root errors as well as transfer function gain and phase errors are presented for a number of traditional and new integration methods. Normalized stability regions in the lambda h plane are compared for the various methods. In particular, it is shown that a modified form of Euler integration with root matching is an especially efficient method for simulating lightly-damped structural modes. The method has been used successfully for structural bending modes in the real-time simulation of missiles. Performance of this algorithm is compared with other special algorithms, including the state-transition method. A predictor-corrector version of the modified Euler algorithm permits it to be extended to the simulation of nonlinear models of the type likely to be obtained when using the discretized structure approach. Performance of the different integration methods is also compared for integration step sizes larger than those for which the asymptotic formulas are valid. It is concluded that many traditional integration methods, such as RD-4, are not competitive in the simulation of lightly damped structures.

Howe, R. M.↗
DOE OSTI2020

Conformation and dynamics of ring polymers under symmetric thin film confinement

Understanding the structure and dynamics of polymers under confinement has been of widespread interest, and one class of polymers that have received comparatively little attention under confinement is that of ring polymers. The properties of non-concatenated ring polymers can also be important in biological fields because ring polymers have been proven to be a good model to study DNA organization in the cell nucleus. From our previous study, linear polymers in a cylindrically confined polymer melt were found to segregate from each other as a result of the strong correlation hole effect that is enhanced by the confining surfaces. By comparison, our subsequent study of linear polymers in confined thin films at similar levels of confinements found only the onset of segregation. Here, we use molecular dynamics simulation to investigate the chain conformations and dynamics of ring polymers under planar (1D) confinement as a function of film thickness. Our results show that conformations of ring polymers are similar to the linear polymers under planar confinement, except that ring polymers are less compressed in the direction normal to the walls. While we find that the correlation hole effect is enhanced under confinement, it is not as pronounced as the linear polymers under 2D confinement. Finally, we show that chain dynamics far above Tg are primarily affected by the friction from walls based on the monomeric friction coefficient we get from the Rouse mode analysis.

36 MATERIALS SCIENCE↗
NASA NTRS2019

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↗
DOE OSTI2020

Dynamic State Estimation for Radial Microgrid Protection

Dynamic state estimation (DSE) is a generalization of differential protection that offers a reduced likelihood of misoperation, particularly in the case of devices with nonlinear characteristics such as transformers which are being energized. It is also useful in cases where distance protection performs poorly, such as transmission lines with series compensation or mutually coupled transmission lines. DSE has been previously applied to microgrid branch protection. This study investigates the use of DSE for protection of radial portions of a microgrid or distribution system. This can be a challenge in microgrids or distribution systems with distributed generation on account of lack of fault current from inverter-interfaced generation, varying fault current between grid- connected and islanded modes, the potential for normally-meshed operation and unbalanced operation due to single-phase loads. Admittance relaying has been investigated as a solution for protection of microgrids, though it has been observed to have issues with grounded-wye connected loads, and additional relaying is necessary to prevent misoperation.

24 POWER TRANSMISSION AND DISTRIBUTION↗