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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Relationship between the contact force strength and numerical inaccuracies in piecewise-smooth systems

This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. Finally, this study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.

42 ENGINEERING↗

NASA NTRS2004

Solar Pumped High Power Solid State Laser for Space Applications

Highly coherent laser light provides a nearly optimal means of transmitting power in space. The simplest most direct means of converting sunlight to coherent laser light is a solar pumped laser oscillator. A key need for broadly useful space solar power is a robust solid state laser oscillator capable of operating efficiently in near Earth space at output powers in the multi hundred kilowatt range. The principal challenges in realizing such solar pumped laser oscillators are: (1) the need to remove heat from the solid state laser material without introducing unacceptable thermal shock, thermal lensing, or thermal stress induced birefringence to a degree that improves on current removal rates by several orders of magnitude and (2) to introduce sunlight at an effective concentration (kW/sq cm of laser cross sectional area) that is several orders of magnitude higher than currently available while tolerating a pointing error of the spacecraft of several degrees. We discuss strategies for addressing these challenges. The need to remove the high densities of heat, e.g., 30 kW/cu cm, while keeping the thermal shock, thermal lensing and thermal stress induced birefringence loss sufficiently low is addressed in terms of a novel use of diamond integrated with the laser material, such as Ti:sapphire in a manner such that the waste heat is removed from the laser medium in an axial direction and in the diamond in a radial direction. We discuss means for concentrating sunlight to an effective areal density of the order of 30 kW/sq cm. The method integrates conventional imaging optics, non-imaging optics and nonlinear optics. In effect we use a method that combines some of the methods of optical pumping solid state materials and optical fiber, but also address laser media having areas sufficiently large, e.g., 1 cm diameter to handle the multi-hundred kilowatt level powers needed for space solar power.

Fork, Richard L.↗

NASA NTRS1995

Lagrangian mixed layer modeling of the western equatorial Pacific

Processes that control the upper ocean thermohaline structure in the western equatorial Pacific are examined using a Lagrangian mixed layer model. The one-dimensional bulk mixed layer model of Garwood (1977) is integrated along the trajectories derived from a nonlinear 1 1/2 layer reduced gravity model forced with actual wind fields. The Global Precipitation Climatology Project (GPCP) data are used to estimate surface freshwater fluxes for the mixed layer model. The wind stress data which forced the 1 1/2 layer model are used for the mixed layer model. The model was run for the period 1987-1988. This simple model is able to simulate the isothermal layer below the mixed layer in the western Pacific warm pool and its variation. The subduction mechanism hypothesized by Lukas and Lindstrom (1991) is evident in the model results. During periods of strong South Equatorial Current, the warm and salty mixed layer waters in the central Pacific are subducted below the fresh shallow mixed layer in the western Pacific. However, this subduction mechanism is not evident when upwelling Rossby waves reach the western equatorial Pacific or when a prominent deepening of the mixed layer occurs in the western equatorial Pacific or when a prominent deepening of the mixed layer occurs in the western equatorial Pacific due to episodes of strong wind and light precipitation associated with the El Nino-Southern Oscillation. Comparison of the results between the Lagrangian mixed layer model and a locally forced Eulerian mixed layer model indicated that horizontal advection of salty waters from the central Pacific strongly affects the upper ocean salinity variation in the western Pacific, and that this advection is necessary to maintain the upper ocean thermohaline structure in this region.

Shinoda, Toshiaki↗

NASA NTRS1994

The fine structure of Langmuir waves observed upstream of the bow shock at Venus

Highly structured Langmuir waves, also known as electron plasma oscillations, have been observed in the foreshock of Venus using the plasma wave experiment on the Galileo spacecraft during the gravity assist flyby on February 10, 1990. The Galileo wideband sampling system provides digital electric field waveform measurements at sampling rates up to 201,600 samples per second, much higher than any previous instrument of this type. The main Langmuir wave emission band occurs near the local electron plasma frequency, which was approximately 43 kHz. The Langmuir waves are observed to shift above and below the plasma frequency, sometimes by as much as 20 kHz. The shifts in frequency are closely correlated with the downstream distance from the tangent field line, implying that the shifts are controlled by the electron beam velocity. Considerable fine structure is also evident, with time scales as short as 0.15 milliseconds, corresponding to spatial scales of a few tens of Debye lengths. The frequency spectrum often consists of beat-type waveforms, with beat frequencies ranging from 0.2 to 7 kHz, and in a few cases, isolated wavepackets. The peak electric field strengths are approximately 1 mV/m. These field strengths are too small for strongly nonlinear processes to be important. The beat-type waveforms are suggestive of a parametric decay process.

Hospodarsky, G. B.↗

NASA NTRS1994

Fine structure of Langmuir waves observed upstream of the bow shock at Venus

Highly structured Langmuir waves, also known as electron plasma oscillations, have been observed in the foreshock of Venus using the plasma wave experiment on the Galileo spacecraft during the gravity assist flyby on February 10, 1990. The Galileo wideband sampling system provides digital electric field waveform measurements at sampling rates up to 201,600 samples per second, much higher than any previous instrument of this type. The main Langmuir wave emission band occurs near the local electron plasma frequency, which was approximately 43 kHz. The Langmuir waves are observed to shift above and below the plasma frequency, sometimes by as much as 20 kHz. The shifts in frequency are closely correlated with the downstream distance from the tangent field line, implying that the shifts are controlled by the electron beam velocity. Considerable fine structure is also evident, with timescales as short as 0.15 ms, corresponding to spatial scales of a few tens of Debye lengths. The frequency spectrum often consists of beat-type waveforms, with beat frequencies ranging from 0.2 to 7 kHz, and in a few cases, isolated wave packets. The peak electric field strengths are approximately 1 mV/m. These field strengths are too small for strongly nonlinear processes to be important. The beat-type waveforms are suggestive of a parametric decay process.

Hospodarsky, G. B.↗

NASA NTRS1998

Turbulent Compressible Convection with Rotation: Mean Flows and Differential Rotation - 2

The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Brummell, Nicholas H.↗

NASA NTRS2003

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

Negative Differential Resistance (NDR) frequency conversion with gain

The dependence of the I-V characteristic of the negative differential resistance (NDR) devices on the power level and frequency of the rf input signal has been theoretically analyzed with a modified large- and small-signal nonlinear circuit analysis program. The NDR devices we used in this work include both the tunnel diode (without the antisymmetry in the I-V characteristic) and resonant-tunneling devices (with the antisymmetry in the I-V characteristic). Absolute negative conductance can be found from a zero-biased resonant tunneling device when the applied pump power is within a small range. This study verifies the work of Sollner et al. Variable negative conductances at the fundamental and harmonic frequencies can also be obtained from both the unbiased and biased tunnel diodes. The magnitude of the negative conductances can be adjusted by varying the pump amplitude -- a very useful circuit property. However, the voltage range over which the negative conductance occurs moves towards the more positive side of the voltage axis with increasing frequency. Furthermore, the range of the pumping amplitude to obtain negative conductance varies with the parasitics (resistance and capacitance) of the device. The theoretical observation of the dependence of the I-V characteristic of the NDR devices on the power and frequency of the applied pump signal is supported by the experimental results. In addition, novel functions of a NDR device such as self-oscillating frequency multiplier and mixer with gain have been experimentally demonstrated. The unbiased oscillator have also been successfully realized with a NDR device with an antisymmetrical I-V characteristic. Finally, the applications of these device functions will be discussed.

Hwu, R. J.↗

DOE OSTI2022

Learning Distributed Geometric Koopman Operator for Sparse Networked Dynamical Systems

Koopman operator theory provides an alternative to study nonlinear networked dynamical systems by mapping the state space to an abstract higher dimensional space where the system evolution is linear. Recent works show the application of graph neural networks (GNNs) to learn state to object-centric embeddings and achieve centralized block-wise computation of Koopman operator (KO) under additional assumptions on the underlying agents properties and constraints on the KO structure. However, the computational complexity of learning the Koopman increases exponentially for networked systems where the number of possible system states grows in a combinatorial fashion with the number of nodes. The learning challenge is further amplified for sparse networks by two factors: 1) sample sparsity for learning the Koopman operator in the non-linear space, and 2) the divergence in the dynamics of individual nodes or from one subgraph to another. Our work aims to address these challenge by formulating the representation learning of networked dynamical systems into a multi-agent paradigm and learning the Koopman operator in a distributive manner. The computational as well as performance advantages of distributed Koopman is predominant for sparse networks whereas for fully connected networks, it is shown to coincide with the centralized one. The empirical study on rope system, network of oscillators and a synthetic power system show comparable and superior performance along with computational benefits with the state-of-the-art methods.

Mukherjee, Sayak↗

NASA NTRS1989

Currents between tethered electrodes in a magnetized laboratory plasma

Laboratory experiments on important plasma physics issues of electrodynamic tethers were performed. These included current propagation, formation of wave wings, limits of current collection, nonlinear effects and instabilities, charging phenomena, and characteristics of transmission lines in plasmas. The experiments were conducted in a large afterglow plasma. The current system was established with a small electron-emitting hot cathode tethered to an electron-collecting anode, both movable across the magnetic field and energized by potential difference up to V approx.=100 T(sub e). The total current density in space and time was obtained from complete measurements of the perturbed magnetic field. The fast spacecraft motion was reproduced in the laboratory by moving the tethered electrodes in small increments, applying delayed current pulses, and reconstructing the net field by a linear superposition of locally emitted wavelets. With this technique, the small-amplitude dc current pattern is shown to form whistler wings at each electrode instead of the generally accepted Alfven wings. For the beam electrode, the whistler wing separates from the field-aligned beam which carries no net current. Large amplitude return currents to a stationary anode generate current-driven microinstabilities, parallel electric fields, ion depletions, current disruptions and time-varying electrode charging. At appropriately high potentials and neutral densities, excess neutrals are ionized near the anode. The anode sheath emits high-frequency electron transit-time oscillations at the sheath-plasma resonance. The beam generates Langmuir turbulence, ion sound turbulence, electron heating, space charge fields, and Hall currents. An insulated, perfectly conducting transmission line embedded in the plasma becomes lossy due to excitation of whistler waves and magnetic field diffusion effects. The implications of the laboratory observations on electrodynamic tethers in space are discussed.

Stenzel, R. L.↗

DOE OSTI2021

Driven-dissipative phase transition in a Kerr oscillator: From semiclassical PT symmetry to quantum fluctuations

We study a minimal model that has a driven-dissipative quantum phase transition, namely a Kerr non-linear oscillator subject to driving and dissipation. Using mean-field theory, exact diagonalization, and the Keldysh formalism, we analyze the critical phenomena in this system, showing which aspects can be captured by each approach and how the approaches complement each other. Then critical scaling and finite-size scaling are calculated analytically using the quantum Langevin equation. Furthermore, the physics contained in this simple model is surprisingly rich: it includes a continuous phase transition, Z 2 symmetry breaking, PT symmetry, state squeezing, and critical fluctuations. Due to its simplicity and solvability, this model can serve as a paradigm for exploration of open quantum many-body physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Quantum scars of bosons with correlated hopping

Abstract Recent experiments on Rydberg atom arrays have found evidence of anomalously slow thermalization and persistent density oscillations, which have been interpreted as a many-body analog of the phenomenon of quantum scars. Periodic dynamics and atypical scarred eigenstates originate from a “hard” kinetic constraint: the neighboring Rydberg atoms cannot be simultaneously excited. Here we propose a realization of quantum many-body scars in a 1D bosonic lattice model with a “soft” constraint in the form of density-assisted hopping. We discuss the relation of this model to the standard Bose-Hubbard model and possible experimental realizations using ultracold atoms. We find that this model exhibits similar phenomenology to the Rydberg atom chain, including weakly entangled eigenstates at high energy densities and the presence of a large number of exact zero energy states, with distinct algebraic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Entangled two-photon absorption with Brownian-oscillator fluctuations

Here, we theoretically investigate the two-photon absorption signals of a three-band (g, e, f) system diagonally coupled to an over-damped Brownian oscillator bath, which induces random Gaussian modulations of energy levels with an arbitrary degree of correlation. For fast modulation, extra 2ω eg and 2ω fe peaks may obscure the g–f transitions in the classical two-photon absorption (CTPA) spectra for nearly resonant e states. These peaks arise from one-photon resonant g–e or e–f transitions. In the slow modulation limit, these peaks vanish because of the short tails of the Gaussian line shape. CTPA strongly depends on the correlations between energy fluctuations. In entangled two-photon absorption, the extra peaks are eliminated because of the broad one-photon but narrow two-photon spectrum of the twin photons. The variation of the coherences between f states with the correlation between energy fluctuations is explored.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Producing entangled photon pairs and quantum squeezed states in plasmas

Plasma is capable of mediating the conversion of two pump photons into two different photons through a relativistic four-wave mixing nonlinearity. Spontaneously created photon pairs are emitted at symmetric angles with respect to the colinear pump direction, and the emission rate is largest if they have identical frequency. Thus, two orthogonally polarized pumps can produce polarization-entangled photon pairs through a millimeter-long homogeneous plasma. Here, the noise from Raman scattering can be avoided if the pump detuning differs from twice the plasma frequency. However, pump detuning exactly equal to twice the plasma frequency can significantly enhance the interaction rate, which allows for the production of strong two-mode squeezed states. Remarkably, the amplified noise from Raman scattering are correlated and hence can be suppressed in one of the output quadratures, thereby maintaining the squeezing magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Analytical WKB theory for high-harmonic generation and its application to massive Dirac electrons

Here we propose an analytical approach to high-harmonic generation (HHG) for nonperturbative low-frequency and high-intensity fields based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation. By properly taking into account Stokes phenomena of WKB solutions, we obtain wave functions that systematically include the repetitive dynamics of production and acceleration of electron-hole pairs and quantum interference due to phase accumulation between different pair production times (Stückelberg phase). Using the obtained wave functions without relying on any phenomenological assumptions, we explicitly compute electric current (including intra- and interband contributions) as the source of HHG for a massive Dirac system in (1+1) dimensions under an ac electric field. We demonstrate that the WKB approximation agrees well with numerical results obtained by solving the time-dependent Schrödinger equation and point out that the quantum interference is important in HHG. We also predict in the deep nonperturbative regime that (1) harmonic intensities oscillate with respect to electric-field amplitude E 0 and frequency Ω with a period determined by the Stückelberg phase, (2) the cutoff order of HHG is determined by 2eE 0 ℏΩ 2 , with e being the electron charge, and that (3) noninteger harmonics, controlled by the Stückelberg phase, appear as a transient effect. Our WKB theory is particularly suited for a parameter regime, where the Keldysh parameter γ = (Δ/2)Ω/eE 0 ,with Δ being the gap size, is small. This parameter regime corresponds to intense lasers in the terahertz regime for realistic massive Dirac materials. Our analysis implies that the so-called HHG plateau can be observed at the terahertz frequency within the current technology.

36 MATERIALS SCIENCE↗

NASA NTRS2017

Simulation Study of Landau Damping Near the Persisting to Arrested Transition

A one-dimensional electrostatic filtered Vlasov-Poisson simulation study is discussed. The transition from persisting to arrested Landau damping that is produced by increasing the strength of a sinusoidal perturbation on a background Vlasov-Poisson equilibrium is explored. Emphasis is placed on observed features of the electron phase-space distribution when the perturbation strength is near the transition value. A single ubiquitous waveform is found perturbing the space-averaged phase-space distribution at almost any time in all of the simulations; the sole exception is the saturation stage that can occur at the end of the arrested damping scenario. This waveform contains relatively strong, very narrow structures in velocity bracketing (plus or minus velocity (sub res)) - the velocities at which electrons must move to traverse the dominant field mode wavelength in one of its oscillation periods - and propagating with (plus or minus velocity (sub res)) respectively. Local streams of electrons are found in these structures crossing the resonant velocities from low speed to high speed during Landau damping and from high speed to low speed during Landau growth. At the arrest time, when the field strength is briefly constant, these streams vanish. It is conjectured that the expected transfer of energy between electrons and field during Landau growth or damping has been visualized for the first time. No evidence is found in the phase-space distribution to support recent well-established discoveries of a second-order phase transition in the electric field evolution. While trapping is known to play a role for larger perturbation strengths, it is shown that trapping plays no role at any time in any of the simulations near the transition perturbation strength.

plasma simulation↗

DOE OSTI2022

Dynamics and bifurcations of laminar annular swirling and non-swirling jets

This paper presents bifurcation analyses characterising the nonlinear dynamics of fully developed laminar annular jets with respect to the centrebody diameter ${ {d}}$ , Reynolds number ${ {Re}}$ , and swirl ratio ${ {S}}$ . Similar flows appear in numerous applications and feature a vibrant range of topological and dynamical characteristics associated with phenomena including shear layer separation and vortex breakdown. Our results begin by describing the non-monotonic evolution of the axisymmetric jet's steady topology under varying ${ {S}}$ . In accord with earlier reports, the jet progresses through a sequence of wake, breakdown and wall jet regimes in a qualitatively similar manner across a wide span of ${ {d}}$ and ${ {Re}}$ values. In the wake regime, the non-swirling jet bifurcates to a plane-symmetric, but not axisymmetric, steady flow pattern beyond a ${ {d}}$ -dependent critical ${ {Re}}$ value. With further increase in ${ {Re}}$ , the steady non-swirling jet destabilises subsequently via multiple distinct Hopf bifurcations. Introducing ${ {S}}>0$ to the jet also induces unsteadiness by twisting the singly azimuthally periodic ( $|m|=1$ ) asymmetric wake structure and causing it to precess periodically in time about the central axis. Intermediate swirl stabilises this unsteady dynamics and restores the jet's axisymmetry. This stabilising effect is then reversed in the breakdown regime at higher ${ {S}}$ , where a variety of different $|m|=1$ and $|m|=2$ instabilities bifurcate from the steady flow as ${ {S}}$ is increased. Several instances of hysteresis and subcritical behaviour are reported and discussed, including one that manifests precessing vortex core oscillations.

Douglas, Christopher M. (ORCID:0000000259683315)↗