Multifunction tests of a frequency domain based flutter suppression system

The process is described of analysis, design, digital implementation, and subsonic testing of an active control flutter suppression system for a full span, free-to-roll wind tunnel model of an advanced fighter concept. The design technique uses a frequency domain representation of the plant and used optimization techniques to generate a robust multi input/multi output controller. During testing in a fixed-in-roll configuration, simultaneous suppression of both symmetric and antisymmetric flutter was successfully shown. For a free-to-roll configuration, symmetric flutter was suppressed to the limit of the tunnel test envelope. During aggressive rolling maneuvers above the open-loop flutter boundary, simultaneous flutter suppression and maneuver load control were demonstrated. Finally, the flutter damping controller was reoptimized overnight during the test using combined experimental and analytical frequency domain data, resulting in improved stability robustness

Do Finite-Size Lyapunov Exponents Detect Coherent Structures?

Ridges of the Finite-Size Lyapunov Exponent (FSLE) field have been used as indicators of hyperbolic Lagrangian Coherent Structures (LCSs). A rigorous mathematical link between the FSLE and LCSs, however, has been missing. Here we prove that an FSLE ridge satisfying certain conditions does signal a nearby ridge of some Finite-Time Lyapunov Exponent (FTLE) field, which in turn indicates a hyperbolic LCS under further conditions. Other FSLE ridges violating our conditions, however, are seen to be false positives for LCSs. We also find further limitations of the FSLE in Lagrangian coherence detection, including ill-posedness, artificial jump-discontinuities, and sensitivity with respect to the computational time step.Comment: 22 pages, 7 figures, v3: corrects the z-axis labels of Fig. 2 (left) that appears in the version published in Chao

Treatment of input uncertainty in hydrologic modeling: Doing hydrology backward with Markov chain Monte Carlo simulation

There is increasing consensus in the hydrologic literature that an appropriate framework for streamflow forecasting and simulation should include explicit recognition of forcing and parameter and model structural error. This paper presents a novel Markov chain Monte Carlo (MCMC) sampler, entitled differential evolution adaptive Metropolis (DREAM), that is especially designed to efficiently estimate the posterior probability density function of hydrologic model parameters in complex, high-dimensional sampling problems. This MCMC scheme adaptively updates the scale and orientation of the proposal distribution during sampling and maintains detailed balance and ergodicity. It is then demonstrated how DREAM can be used to analyze forcing data error during watershed model calibration using a five-parameter rainfall-runoff model with streamflow data from two different catchments. Explicit treatment of precipitation error during hydrologic model calibration not only results in prediction uncertainty bounds that are more appropriate but also significantly alters the posterior distribution of the watershed model parameters. This has significant implications for regionalization studies. The approach also provides important new ways to estimate areal average watershed precipitation, information that is of utmost importance for testing hydrologic theory, diagnosing structural errors in models, and appropriately benchmarking rainfall measurement devices

On Norm-Based Estimations for Domains of Attraction in Nonlinear Time-Delay Systems

For nonlinear time-delay systems, domains of attraction are rarely studied despite their importance for technological applications. The present paper provides methodological hints for the determination of an upper bound on the radius of attraction by numerical means. Thereby, the respective Banach space for initial functions has to be selected and primary initial functions have to be chosen. The latter are used in time-forward simulations to determine a first upper bound on the radius of attraction. Thereafter, this upper bound is refined by secondary initial functions, which result a posteriori from the preceding simulations. Additionally, a bifurcation analysis should be undertaken. This analysis results in a possible improvement of the previous estimation. An example of a time-delayed swing equation demonstrates the various aspects.Comment: 33 pages, 8 figures, "This is a pre-print of an article published in 'Nonlinear Dynamics'. The final authenticated version is available online at https://doi.org/10.1007/s11071-020-05620-8

A bounded upwinding scheme for computing convection-dominated transport problems

A practical high resolution upwind differencing scheme for the numerical solution of convection-dominated transport problems is presented. The scheme is based on TVD and CBC stability criteria and is implemented in the context of the finite difference methodology. The performance of the scheme is investigated by solving the 1D/2D scalar advection equations, 1D inviscid Burgers’ equation, 1D scalar convection–diffusion equation, 1D/2D compressible Euler’s equations, and 2D incompressible Navier–Stokes equations. The numerical results displayed good agreement with other existing numerical and experimental data

A fast Bayesian approach to discrete object detection in astronomical datasets - PowellSnakes I

A new fast Bayesian approach is introduced for the detection of discrete objects immersed in a diffuse background. This new method, called PowellSnakes, speeds up traditional Bayesian techniques by: i) replacing the standard form of the likelihood for the parameters characterizing the discrete objects by an alternative exact form that is much quicker to evaluate; ii) using a simultaneous multiple minimization code based on Powell's direction set algorithm to locate rapidly the local maxima in the posterior; and iii) deciding whether each located posterior peak corresponds to a real object by performing a Bayesian model selection using an approximate evidence value based on a local Gaussian approximation to the peak. The construction of this Gaussian approximation also provides the covariance matrix of the uncertainties in the derived parameter values for the object in question. This new approach provides a speed up in performance by a factor of `hundreds' as compared to existing Bayesian source extraction methods that use MCMC to explore the parameter space, such as that presented by Hobson & McLachlan. We illustrate the capabilities of the method by applying to some simplified toy models. Furthermore PowellSnakes has the advantage of consistently defining the threshold for acceptance/rejection based on priors which cannot be said of the frequentist methods. We present here the first implementation of this technique (Version-I). Further improvements to this implementation are currently under investigation and will be published shortly. The application of the method to realistic simulated Planck observations will be presented in a forthcoming publication.Comment: 30 pages, 15 figures, revised version with minor changes, accepted for publication in MNRA