Analysis of Transonic Limit Cycle Oscillations under Uncertainty

For the computation of limit cycle oscillations (LCO) at transonic speeds, CFD is required to capture the nonlinear flow features present. The Harmonic Balance method provides an effective means for the computation of LCOs and this paper exploits its efficiency to investigate the impact of variability (both structural a nd aerodynamic) on the aeroelastic behaviour of a 2 dof aerofoil. A Harmonic Balance inviscid CFD solver is coupled with the structural equations and is validated against time marching analyses. Polynomial chaos expansions are employed for the stochastic investiga tion as a faster alternative to Monte Carlo analysis. Adaptive sampling is employed when discontinuities are present. Uncertainties in aerodynamic parameters are looked at first followed by the inclusion of structural variability. Results show the nonlinear effect of Mach number and it’s interaction with the structural parameters on supercritical LCOs. The bifurcation boundaries are well captured by the polynomial chaos.

Prediction of Limit Cycle Oscillations Using an Implicit Aeroelastic-Harmonic Balance Method

This work proposes a extends a novel approach to compute tran sonic Limit Cycle Oscillations using high fidelity analysis. CFD based Harmonic Balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a methodology to determine the unknown frequency of oscillations using an implicit for- mulation of the HB method to accurately capture Limit Cycle Oscillations (LCOs); this is achieved by defining a frequency updating procedure based on a coupled CFD/CSD Harmonic Balance formulation to find the LCO condition. A pitch/plunge aerofoil and respective linear structural models is used to exercise the new method. Results show consistent agreement between the proposed and time-marching methods for both LCO amplitude and frequency.

Investigation of Limit Cycle Oscillations Using Aeroelastic-Harmonic Balance Method

This work investigates limit cycle oscillations in the transonic regime. A novel approach to predict Limit Cycle Oscillations using high fidelity analysis is exploited to accelerate calculations. The method used is an Aeroeasltic Harmonic Balance approach, which has been proven to be efficient and able to predict periodic phenomena. The behaviour of limit cycle oscillations is analysed using uncertainty quantification tools based on polynomial chaos expansions. To improve the efficiency of the sampling process for the polynomial-chaos expansions an adaptive sampling procedure is used. These methods are exercised using two problems: a pitch/plunge aerofoil and a delta-wing. Results indicate that Mach n. variability is determinant to the amplitude of the LCO for the 2D test case, whereas for the wing case analysed here, variability in the Mach n. has an almost negligible influence in amplitude variation and the LCO frequency variability has an almost linear relation with Mach number. Further test cases are required to understand the generality of these results.

Prediction of limit cycle oscillations under uncertainty using a Harmonic Balance method

The Harmonic Balance method is an attractive solution for computing periodic responses and can be an alternative to time domain methods, at a reduced computational cost. The current paper investigates using a Harmonic Balance method for simulating limit cycle oscillations under uncertainty. The Harmonic Balance method is used in conjunction with a non-intrusive polynomial-chaos approach to propagate variability and is validated against Monte Carlo analysis. Results show the potential of the approach for a range of nonlinear dynamical systems, including a full wing configuration exhibiting supercritical and subcritical bifurcations, at a fraction of the cost of performing time domain simulations.

Non-Markovian qubit dynamics in a circuit-QED setup

We consider a circuit-QED setup that allows the induction and control of non-Markovian dynamics of a qubit. Non-Markovianity is enforced over the qubit by means of its direct coupling to a bosonic mode which is controllably coupled to another qubit-mode system. We show that this configuration can be achieved in a circuit-QED setup consisting of two initially independent superconducting circuits, each formed by one charge qubit and one transmission-line resonator, which are put in interaction by coupling the resonators to a current-biased Josephson junction. We solve this problem exactly and then proceed with a thorough investigation of the emergent non-Markovianity in the dynamics of the qubits. Our study might serve the context for the first experimental assessment of non-Markovianity in a multielement solid-state device.

Turbulent Fluctuations in G-band and K-line Intensities Observed with the Rapid Oscillations in the Solar Atmosphere (ROSA) Instrument

Using the Rapid Oscillation in the Solar Atmosphere (ROSA) instrument at the Dunn Solar Telescope we have found that the spectra of fluctuations of the G-band (cadence 1.05 s) and Ca II K-line (cadence 4.2 s) intensities show correlated fluctuations above white noise out to frequencies beyond 300 mHz and up to 70 mHz, respectively. The noise-corrected G-band spectrum presents a scaling range (Ultra High Frequency “UHF”) for f = 25-100 mHz, with an exponent consistent with the presence of turbulent motions. The UHF power, is concentrated at the locations of magnetic bright points in the intergranular lanes, it is highly intermittent in time and characterized by a positive kurtosis κ. Combining values of G-band and K-line intensities, the UHF power, and κ, reveals two distinct “states” of the internetwork solar atmosphere. State 1, with κ ≍ 6, which includes almost all the data, is characterized by low intensities and low UHF power. State 2, with κ ≍ 3, including a very small fraction of the data, is characterized by high intensities and high UHF power. Superposed epoch analysis shows that for State 1, the K-line intensity presents 3.5 min chromospheric oscillations with maxima occurring 21 s after G-band intensity maxima implying a 150-210 km effective height difference. For State 2, the G-band and K-line intensity maxima are simultaneous, suggesting that in the highly magnetized environment sites of G-band and K-line emission may be spatially close together. Analysis of observations obtained with Hinode/SOT confirm a scaling range in the G-band spectrum up to 53 mHz also consistent with turbulent motions as well as the identification of two distinct states in terms of the H-line intensity and G-band power as functions of G-band intensity.

Rapid Oscillations in the Solar Atmosphere: Spectra and Physical Effects

High-frequency fluctuations are observed with the Rapid Oscillations in the Solar Atmosphere (ROSA) instrument (Jess et al. 2010, Solar Phys, 261, 363) at the Dunn Solar Telescope. This can produce simultaneous observations in up to six channels, at different heights in the photosphere and chromosphere, at an unprecedentedly high cadence of 0.5 seconds, and at a spatial resolution of 100 km after photometrically correct speckle reconstruction. Here we concentrate on observations at two levels. The first is in the G-band of the CH radical at 4305.5Å, bandpass 9.2Å, with height of formation z <250 km at a cadence of 0.525 sec corresponding to Nyquist frequency 950 mHz. The second is in the Ca II K-line core at 3933.7Å, bandpass 1.0Å, with height of formation z <1300 km, and cadence 4.2 sec giving Nyquist frequency 120 mHz. The data span 53 min, and the maximum field of view is 45 Mm. The data were taken on 28 May 2009 in internetwork and network near disk center. Using both Fourier and Morlet wavelet methods we find evidence in the G-band spectra for intensity fluctuations above noise out to frequencies f >> 100 mHz. The K-line signal is noisier and is seen only for f <50 mHz. With wavelet techniques we find that G-band spectral power with 20 <f <100 mHz is clearly concentrated in the intergranular lanes and especially at the locations of magnetic elements indicated by G-band bright points. This wavelet power is highly intermittent in time. By cross-correlating the data we find that pulses of high-frequency G-band power in the photosphere tend to be followed by increases in K-line emission in the chromosphere with a time lag of about 2 min.

Prediction of Transonic Limit-Cycle Oscillations Using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic limit-cycle oscillations using high-fidelity analysis. Computational-Fluid-Dynamics based harmonic balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling harmonic balance methods to accurately capture limit-cycle oscillations; this is achieved by defining a frequency-updating procedure based on a coupled computational-fluid-dynamics/computational-structural-dynamics harmonic balance formulation to find the limit-cycle oscillation condition. A pitch/plunge airfoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both limit-cycle oscillation amplitude and frequency while producing at least a one-order-of-magnitude reduction in computational time.

Wave Damping Observed in Upwardly Propagating Sausage-Mode Oscillations Contained within a Magnetic Pore

We present observational evidence of compressible MHD wave modes propagating from the solar photosphere through to the base of the transition region in a solar magnetic pore. High cadence images were obtained simultaneously across four wavelength bands using the Dunn Solar Telescope. Employing Fourier and wavelet techniques, sausage-mode oscillations displaying significant power were detected in both intensity and area fluctuations. The intensity and area fluctuations exhibit a range of periods from 181 to 412 s, with an average period∼290 s, consistent with the global p-mode spectrum. Intensity and area oscillations present in adjacent band passes were found to be out of phase with one another, displaying phase angles of 6.°12, 5.°82,and 15.°97 between the 4170 Å continuum–G-band,G-band–Na i D1, and Na i D1–Ca ii K heights, respectively, reiterating the presence of upwardly propagating sausage-mode waves. A phase relationship of ∼0° between same-bandpass emission and area perturbations of the pore best categorizes the waves as belonging to the “slow” regime of a dispersion diagram. Theoretical calculations reveal that the waves are surface modes, with initial photospheric energies in excess of 35,000 Wm‑2. The wave energetics indicate a substantial decrease in energy with atmospheric height, confirming that magnetic pores are able to transport waves that exhibit appreciable energy damping, which may release considerable energy into the local chromospheric plasma.

Reducing Parametric Uncertainty in Limit Cycle Oscillations

The assimilation of discrete higher fidelity data points with model predictions can be used to achieve a reduction in the uncertainty of the model input parameters which generate accurate predictions. The problem investigated here involves the prediction of limit-cycle oscillations using a High-Dimensional Harmonic Balance method (HDHB). The efficiency of the HDHB method is exploited to enable calibration of structural input parameters using a Bayesian inference technique. Markov-chain Monte Carlo is employed to sample the posterior distributions. Parameter estimation is carried out on both a pitch/plunge aerofoil and Goland wing configuration. In both cases significant refinement was achieved in the distribution of possible structural parameters allowing better predictions of their
true deterministic values.

Closed circuit rebreathing to achieve inert gas wash-in for multiple breath wash-out

Multiple breath wash-out (MBW) testing requires prior wash-in of inert tracer gas. Wash-in efficiency can be enhanced by a rebreathing tracer in a closed circuit. Previous attempts to deploy this did not account for the impact of CO2 accumulation on patients and were unsuccessful. We hypothesised that an effective rebreathe wash-in could be delivered and it would not alter wash-out parameters. Computer modelling was used to assess the impact of the rebreathe method on wash-in efficiency. Clinical testing of open and closed circuit wash-in–wash-out was performed in healthy controls and adult patients with cystic fibrosis (CF) using a circuit with an effective CO2 scrubber and a refined wash-in protocol. Wash-in efficiency was enhanced by rebreathing. There was no difference in mean lung clearance index between the two wash-in methods for controls (6.5 versus 6.4; p=0.2, n=12) or patients with CF (10.9 versus 10.8; p=0.2, n=19). Test time was reduced by rebreathe wash-in (156 versus 230 s for CF patients, p<0.001) and both methods were well tolerated. End wash-in CO2 was maintained below 2% in most cases. Rebreathe–wash-in is a promising development that, when correctly deployed, reduces wash-in time and facilitates portable MBW testing. For mild CF, wash-out outcomes are equivalent to an open circuit.

Equivalent Circuit Model of Twisted Split Ring Frequency Selective Surfaces

This paper proposes a wideband equivalent circuit model for a twisted split ring frequency selective surface (FSS). Such surfaces can be used for modelling and design of polarisation sensitive surfaces such as circularly polarized selective surfaces as well as structures with asymmetric transmission. The proposed model is based extraction of equivalent circuit parameters from a single split ring (SRR) FSS and magnetic coupling from periodic eigenmode analysis of the coupled SRR. The resulting equivalent circuit model demonstrates excellent agreement with full-wave simulations.