Markov chains, R-trivial monoids and representation theory

We develop a general theory of Markov chains realizable as random walks on R-trivial monoids. It provides explicit and simple formulas for the eigenvalues of the transition matrix, for multiplicities of the eigenvalues via Mobius inversion along a lattice, a condition for diagonalizability of the transition matrix and some techniques for bounding the mixing time. In addition, we discuss several examples, such as Toom-Tsetlin models, an exchange walk for finite Coxeter groups, as well as examples previously studied by the authors, such as nonabelian sandpile models and the promotion Markov chain on posets. Many of these examples can be viewed as random walks on quotients of free tree monoids, a new class of monoids whose combinatorics we develop.

Computation of power released during corona treatment on polymeric insulators under ac and dc excitation

With increasing energy demand, it necessitates to generate and transmit the electrical power with minimal losses. High voltage power transmission is the most economical way of transmitting bulk power over long distances. Transmission insulator is one of the main components used as a mechanical support and to electrically isolate the conductor from the tower. Corona from the hardware and conductors can significantly affect the performance of the polymeric insulators. In the present investigation a methodology is presented to evaluate the corona performance of the polymeric shed material under different environment conditions for both ac and dc excitation. The results of the comprehensive analysis on various polymeric samples and the power released from the corona electrode for both the ac and dc excitation are presented. Some interesting results obtained from the chemical analysis confirmed the presence of nitric acid species on the treated sample which in long term will affect the strength of the insulator, also the morphological changes were found to be varying for different experimental conditions. (C) 2015 The Authors. Published by Elsevier Ltd.

Comparison of Matching Pursuit Algorithm with Other Signal Processing Techniques for Computation of the Time-Frequency Power Spectrum of Brain Signals

Signals recorded from the brain often show rhythmic patterns at different frequencies, which are tightly coupled to the external stimuli as well as the internal state of the subject. In addition, these signals have very transient structures related to spiking or sudden onset of a stimulus, which have durations not exceeding tens of milliseconds. Further, brain signals are highly nonstationary because both behavioral state and external stimuli can change on a short time scale. It is therefore essential to study brain signals using techniques that can represent both rhythmic and transient components of the signal, something not always possible using standard signal processing techniques such as short time fourier transform, multitaper method, wavelet transform, or Hilbert transform. In this review, we describe a multiscale decomposition technique based on an over-complete dictionary called matching pursuit (MP), and show that it is able to capture both a sharp stimulus-onset transient and a sustained gamma rhythm in local field potential recorded from the primary visual cortex. We compare the performance of MP with other techniques and discuss its advantages and limitations. Data and codes for generating all time-frequency power spectra are provided.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

A time domain computation method for dynamic behavior of mooring system

A quasi-steady time domain method is developed for the prediction of dynamic behavior of a mooring system under the environmental disturbances, such as regular or irregular waves, winds and currents. The mooring forces are obtained in a static sense at each instant. The dynamic feature of the mooring cables can be obtained by incorporating the extended 3-D lumped-mass method with the known ship motion history. Some nonlinear effects, such as the influence of the instantaneous change of the wetted hull surface on the hydrostatic restoring forces and Froude-Krylov forces, are included. The computational results show a satisfactory agreement with the experimental ones.

Filtering via approximate Bayesian computation

Approximate Bayesian computation (ABC) has become a popular technique to facilitate Bayesian inference from complex models. In this article we present an ABC approximation designed to perform biased filtering for a Hidden Markov Model when the likelihood function is intractable. We use a sequential Monte Carlo (SMC) algorithm to both fit and sample from our ABC approximation of the target probability density. This approach is shown to, empirically, be more accurate w.r.t.~the original filter than competing methods. The theoretical bias of our method is investigated; it is shown that the bias goes to zero at the expense of increased computational effort. Our approach is illustrated on a constrained sequential lasso for portfolio allocation to 15 constituents of the FTSE 100 share index.

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, K-I and K-II, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Sequential Monte Carlo computation of the score and observed information matrix in state-space models with application to parameter estimation

Sequential Monte Carlo (SMC) methods are popular computational tools for Bayesian inference in non-linear non-Gaussian state-space models. For this class of models, we propose SMC algorithms to compute the score vector and observed information matrix recursively in time. We propose two different SMC implementations, one with computational complexity $\mathcal{O}(N)$ and the other with complexity $\mathcal{O}(N^{2})$ where $N$ is the number of importance sampling draws. Although cheaper, the performance of the $\mathcal{O}(N)$ method degrades quickly in time as it inherently relies on the SMC approximation of a sequence of probability distributions whose dimension is increasing linearly with time. In particular, even under strong \textit{mixing} assumptions, the variance of the estimates computed with the $\mathcal{O}(N)$ method increases at least quadratically in time. The $\mathcal{O}(N^{2})$ is a non-standard SMC implementation that does not suffer from this rapid degrade. We then show how both methods can be used to perform batch and recursive parameter estimation.

用Symbolic Computation计算毛细重力波六阶解

本文首先运用Symbolic Computation在半物理平面(x,)上计算了毛细重力波的六阶解,得到了波形与色散关系,低阶解与 Hogan 结果一致。