Two-Stage Extended Kalman Filters with Derivative-Free Local Linearizations

This paper proposes a derivative-free two-stage extended Kalman filter (2-EKF) especially suited for state and parameter identification of mechanical oscillators under Gaussian white noise. Two sources of modeling uncertainties are considered: (1) errors in linearization, and (2) an inadequate system model. The state vector is presently composed of the original dynamical/parameter states plus the so-called bias states accounting for the unmodeled dynamics. An extended Kalman estimation concept is applied within a framework predicated on explicit and derivative-free local linearizations (DLL) of nonlinear drift terms in the governing stochastic differential equations (SDEs). The original and bias states are estimated by two separate filters; the bias filter improves the estimates of the original states. Measurements are artificially generated by corrupting the numerical solutions of the SDEs with noise through an implicit form of a higher-order linearization. Numerical illustrations are provided for a few single- and multidegree-of-freedom nonlinear oscillators, demonstrating the remarkable promise that 2-EKF holds over its more conventional EKF-based counterparts. DOI: 10.1061/(ASCE)EM.1943-7889.0000255. (C) 2011 American Society of Civil Engineers.

Tracking analysis of an LMS decision feedback equalizer for a wireless channel

We consider a time varying wireless fading channel, equalized by an LMS Decision Feedback equalizer (DFE). We study how well this equalizer tracks the optimal MMSEDFE (Wiener) equalizer. We model the channel by an Autoregressive (AR) process. Then the LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs (ordinary differential equations). Using these ODEs, we show via some examples that the LMS equalizer moves close to the instantaneous Wiener filter after initial transience. We also compare the LMS equalizer with the instantaneous optimal DFE (the commonly used Wiener filter) designed assuming perfect previous decisions and computed using perfect channel estimate (we will call it as IDFE). We show that the LMS equalizer outperforms the IDFE almost all the time after initial transience.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

Simulation of RNA silencing pathway for time-dependent transgene transcription rate

The synthesis of dsRNA is analyzed using a pathway model with amplifications caused by the aberrant RNAs. The transgene influx rate is assumed time-decaying considering the fact that the number of transgenes can not be infinite. The dynamics of the transgene induced RNA silencing is investigated using a system of coupled nonautonomous ordinary nonlinear differential equations which describe the model phenomenologically. The silencing phenomena are detected after a period of transcription. Important contributions of certain parameters are discussed with several numerical examples.

Exact controllability of generalized Hammerstein type equation

In this article, we study the exact controllability of an abstract model described by the controlled generalized Hammerstein type integral equation $$ x(t) = int_0^t h(t,s)u(s)ds+ int_0^t k(t,s,x)f(s,x(s))ds, quad 0 leq t leq T less than infty, $$ where, the state $x(t)$ lies in a Hilbert space $H$ and control $u(t)$ lies another Hilbert space $V$ for each time $t in I=[0,T]$, $T$ greater than 0. We establish the controllability result under suitable assumptions on $h, k$ and $f$ using the monotone operator theory.

Noise reduction with perforated three-duct muffler components

This paper describes the authors’ distributed parameter approach for derivation of closed-form expressions for the four-pole parameters of the perforated three-duct muffler components. In this method, three simultaneous second-order partial differential equations are first reduced to a set of six first-order ordinary differential equations. These equations are then uncoupled by means of a modal matrix. The resulting 6 × 6 matrix is reduced to the 2 × 2 transfer matrix using the relevant boundary conditions. This is combined with transfer matrices of other elements (upstream and downstream of this perforated element) to predict muffler performance like noise reduction, which is also measured. The correlation between experimental and theoretical values of noise reduction is shown to be satisfactory.

Switching Transients in Single-Phase Induction Motors

The subject of transients in polyphase induction motors and synchronous machines has been studied in very great detail by several investigators, but no published literature exists dealing exclusively with the analysis of the problem of transients in single-phase induction motors. This particular problem has been studied in this paper by applying the Laplace transform. The results of actual computation of the currents and developed electrical torque are compared with the data obtained by setting up the integro-differential equations of the machine on an electronic differential analyzer. It is shown that if the motor is switched on to the supply when the potential passes through its zero value, there is a pulsating fundamental frequency torque superimposed on the average steady-state unidirectional torque. If, on the other hand, the switch is closed when the applied potential passes through its maximum value, the developed electrical torque settles down to its final steady-state value during the first cycle of the supply voltage.

The x¿n-x Plane for Analysis of Certain Second-Order Nonlinear Systems

Analysis of certain second-order nonlinear systems, not easily amenable to the phase-plane methods, and described by either of the following differential equations xÿn-2ÿ+ f(x)xÿ2n+g(x)xÿn+h(x)=0 ÿ+f(x)xÿn+h(x)=0 n≫0 can be effected easily by drawing the entire portrait of trajectories on a new plane; that is, on one of the xÿnÿx planes. Simple equations are given to evaluate time from a trajectory on any of these n planes. Poincaré's fundamental phase plane xÿÿx is conceived of as the simplest case of the general xÿnÿx plane.

Techniques for Analysis of Certain Nonlinear Systems

This paper suggests the use of simple transformations like ÿ=kx, kx2 for second-order nonlinear differential equations to effect rapid plotting of the phase-plane trajectories. The method is particularly helpful in determining quickly the trajectory slopes along simple curves in any desired region of the phase plane. New planes such as the tÿ-x, tÿ2-x are considered for the study of some groups of nonlinear time-varying systems. Suggestions for solving certain higher-order nonlinear systems are also made.

Variation of the gas and radiation content in the sub-Keplerian accretion disk around black holes and its impact to the solutions

We investigate the variation of the gas and the radiation pressure in accretion disks during the infall of matter to the black hole and its effect to the flow. While the flow far away from the black hole might be non-relativistic, in the vicinity of the black hole it is expected to be relativistic behaving more like radiation. Therefore, the ratio of gas pressure to total pressure (beta) and the underlying polytropic index (gamma) should not be constant throughout the flow. We obtain that accretion flows exhibit significant variation of beta and then gamma, which affects solutions described in the standard literature based on constant beta. Certain solutions for a particular set of initial parameters with a constant beta do not exist when the variation of beta is incorporated appropriately. We model the viscous sub-Keplerian accretion disk with a nonzero component of advection and pressure gradient around black holes by preserving the conservations of mass, momentum, energy, supplemented by the evolution of beta. By solving the set of five coupled differential equations, we obtain the thermo-hydrodynamical properties of the flow. We show that during infall, beta of the flow could vary up to similar to 300%, while gamma up to similar to 20%. This might have a significant impact to the disk solutions in explaining observed data, e.g. super-luminal jets from disks, luminosity, and then extracting fundamental properties from them. Hence any conclusion based on constant gamma and beta should be taken with caution and corrected. (C) 2011 Elsevier B.V. All rights reserved.

Prediction of Breakout Noise from a Rectangular Duct with Compliant Walls

Breakout noise from HVAC ducts is important at low frequencies, and the coupling between the acoustic waves and the structural waves plays a critical role in the prediction of the transverse transmission loss. This paper describes the analytical calculation of breakout noise by incorporating three-dimensional effects along with the acoustical and structural wave coupling phenomena. The first step in the breakout noise prediction is to calculate the inside duct pressure field and the normal duct wall vibration by using the solution of the governing differential equations in terms of Green's function. The resultant equations are rearranged in terms of impedance and mobility, which results in a compact matrix formulation. The Green's function selected for the current problem is the cavity Green's function with modification of wave number in the longitudinal direction in order to incorporate the terminal impedance. The second step is to calculate the radiated sound power from the compliant duct walls by means of an ``equivalent unfolded plate'' model. The transverse transmission loss from the duct walls is calculated using the ratio of the incident power due to surface source inside the duct to the acoustic power radiated from the compliant duct walls. Analytical results are validated with the FE-BE numerical models.

Five-stage Milstein methods for SDEs

In this paper, we consider the problem of computing numerical solutions for Ito stochastic differential equations (SDEs). The five-stage Milstein (FSM) methods are constructed for solving SDEs driven by an m-dimensional Wiener process. The FSM methods are fully explicit methods. It is proved that the FSM methods are convergent with strong order 1 for SDEs driven by an m-dimensional Wiener process. The analysis of stability (with multidimensional Wiener process) shows that the mean-square stable regions of the FSM methods are unbounded. The analysis of stability shows that the mean-square stable regions of the methods proposed in this paper are larger than the Milstein method and three-stage Milstein methods.

Longitudinal Magnetic Field Effect on Nonlocal Ultrasonic Vibration Analysis of Single-Walled Carbon Nanotubes Based on Wave Propagation Approach

This paper studies the effect of longitudinal magnetic field on ultrasonic vibration in single walled carbon nanotubes (CNTs) based on nonlocal continuum medium theory. Governing partial differential equations of CNTs are derived by considering the Lorentz magnetic forces applied on CNTs induced by a longitudinal magnetic field through Maxwell equations. The vibration characteristics of CNTs under a longitudinal magnetic field are obtained by solving the governing equations via wave propagation approach. The effects of longitudinal magnetic field on vibration of CNTs are discussed through numerical experiments. The present analysis show that vibration frequencies of CNTs drops dramatically in the presence of the magnetic field for various circumferential wavenumbers. Such effect is also observed for various boundary conditions of the CNT. New features for the effect of longitudinal magnetic field on ultrasonic vibration of CNTs, presented in this paper are useful in the design of nano-drive device, nano-oscillator and actuators and nano-electron technology, where carbon nanotubes act as basic elements.

Determination of alpha(s)(M-tau(2)) from improved fixed order perturbation theory

We revisit the extraction of alpha(s)(M-tau(2)) from the QCD perturbative corrections to the hadronic tau branching ratio, using an improved fixed-order perturbation theory based on the explicit summation of all renormalization-group accessible logarithms, proposed some time ago in the literature. In this approach, the powers of the coupling in the expansion of the QCD Adler function are multiplied by a set of functions D-n, which depend themselves on the coupling and can be written in a closed form by iteratively solving a sequence of differential equations. We find that the new expansion has an improved behavior in the complex energy plane compared to that of the standard fixed-order perturbation theory (FOPT), and is similar but not identical to the contour-improved perturbation theory (CIPT). With five terms in the perturbative expansion we obtain in the (MS) over bar scheme alpha(s)(M-tau(2)) = 0.338 +/- 0.010, using as input a precise value for the perturbative contribution to the hadronic width of the tau lepton reported recently in the literature.

SCALE EFFECTS ON ULTRASONIC WAVE DISPERSION CHARACTERISTICS OF MONOLAYER GRAPHENE EMBEDDED IN AN ELASTIC MEDIUM

Ultrasonic wave propagation in a graphene sheet, which is embedded in an elastic medium, is studied using nonlocal elasticity theory incorporating small-scale effects. The graphene sheet is modeled as an one-atom thick isotropic plate and the elastic medium/substrate is modeled as distributed springs. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. After that, an ultrasonic type of wave propagation model is also derived. The explicit expressions for the cut-off frequencies are also obtained as functions of the nonlocal scaling parameter and the y-directional wavenumber. Local elasticity shows that the wave will propagate even at higher frequencies. But nonlocal elasticity predicts that the waves can propagate only up to certain frequencies (called escape frequencies), after which the wave velocity becomes zero. The results also show that the escape frequencies are purely a function of the nonlocal scaling parameter. The effect of the elastic medium is captured in the wave dispersion analysis and this analysis is explained with respect to both local and nonlocal elasticity. The simulations show that the elastic medium affects only the flexural wave mode in the graphene sheet. The presence of the elastic matrix increases the band gap of the flexural mode. The present results can provide useful guidance for the design of next-generation nanodevices in which graphene-based composites act as a major element.