Forced response statistics of mistuned bladed disks: a stochastic reduced basis approach

This paper presents a stochastic reduced basis approach for predicting the forced response statistics of mistuned bladed-disk assemblies. In this approach, the system response in the frequency domain is represented using a linear combination of complex stochastic basis vectors with undermined coefficients. The terms of the preconditioned stochastic Krylov subspace are used here as basis vectors. Two variants of the stochastic Bubnov-Galerkin scheme are employed for computing the undetermined terms in the reduced basis representation, which arise from how the condition for orthogonality between two random vectors is interpreted. Explicit expressions for the response quantities can then be derived in terms of the random system parameters, which allow for the possibility of efficiently computing the response statistics in the post-processing stage. Numerical studies are presented for mistuned cyclic assemblies of mono-coupled single-mode components. It is demonstrated that the accuracy of the response statistical moments computed using stochastic reduced basis methods can be orders of magnitude better than classical perturbation methods.     

Approximate vacuums in (:φ2m :g(x))2 field theories and perturbation series

A class of perturbation problems is considered, in which the Rayleigh-Schrödinger perturbation series for the ground state eigenvalue and eigenvector are presumed to diverge. This class includes the (:^2m :g(x))₂, (m=2, 3) quantum field theory models and the quantum mechanical anharmonic oscillator. It is shown that, using matrix elements and vectors which occur in the series coefficients, one may construct convergent approximants to the eigenvalue and eigenvector. Results of a calculation of the ground state energy of thex ⁴ anharmonic oscillator are given.Supported in part by the National Research Council of Canada.     

Application of h-adaptive,high order finite element method to solve radial Schrödinger equation

The h-adaptive, high order finite element method is applied to solve a second order one dimension eigenvalue problem. The finite element formulation for the Lobatto basis is given, for which basis functions of arbitrary order can be constructed. The adaptive algorithm is simple, yet very efficient and straightforward to implement. The algorithm is based on the observation that the expansion coefficients of Lobatto basis functions decay rapidly. It allows evaluating the smallest eigenvalues simultaneously with the comparable accuracy for all eigenvalues. The presented algorithm is applied to solve the radial Schrödinger equation with the Coulomb and the Woods–Saxon potentials. For both potentials the convergence rate is presented. After seven adaptive iterations nine-digit accuracy was obtained.     

Depolarization by aerosols: Entropy of the Amsterdam light scattering database

In this paper we apply an entropy analysis to measured scattering matrices from the Amsterdam light scattering database. We select examples of mineral aerosols from the database and use them to demonstrate differences in polarization behavior between the particle clouds using a new coherency matrix formulation. These differences are further investigated by analyzing the polarized component of the matrices via two new eigenvector parameters, which can be mapped conveniently onto the surface of a sphere, analogous to the Poincaré sphere used for wave states. We conclude by considering the potential for discriminating different aerosols on the basis of their entropy/eigenvector signatures by solving the contrast optimization problem for clouds with different scattering matrices by using a novel generalized coherency eigenvalue formulation.     

基于随机矩阵理论的城市人群移动行为分析

移动通信应用为人类移动规律的研究提供了独特的数据来源. 本文通过城市手机用户的分布数据,研究城市移动人群的整体动力学行为. 借助随机矩阵理论的方法,通过比较移动人群数据与随机数据在互相关矩阵谱分布上的差异,发现移动人群数据互相关矩阵的相关系数均值、最大本征值及其对应的本征向量明显偏离于随机互相关矩阵的分布,指出这种差异体现了城市移动人群的整体行为特性,且这种差异在不同时间段也会有所不同. 研究结果体现出相关系数的均值和最大本征值的波动趋势,并指出本征向量成员权重的时空模式与城市移动人群整体行为特征的波动过 关键词： 随机矩阵理论 移动人群 宏观行为     

A method for solving stochastic equations by reduced order models and local approximations

A method is proposed for solving equations with random entries, referred to as stochastic equations (SEs). The method is based on two recent developments. The first approximates the response surface giving the solution of a stochastic equation as a function of its random parameters by a finite set of hyperplanes tangent to it at expansion points selected by geometrical arguments. The second approximates the vector of random parameters in the definition of a stochastic equation by a simple random vector, referred to as stochastic reduced order model (SROM), and uses it to construct a SROM for the solution of this equation.The proposed method is a direct extension of these two methods. It uses SROMs to select expansion points, rather than selecting these points by geometrical considerations, and represents the solution by linear and/or higher order local approximations. The implementation and the performance of the method are illustrated by numerical examples involving random eigenvalue problems and stochastic algebraic/differential equations. The method is conceptually simple, non-intrusive, efficient relative to classical Monte Carlo simulation, accurate, and guaranteed to converge to the exact solution.     

Modified perturbation method for eigenvalues of structure with interval parameters

In overcoming the drawbacks of traditional interval perturbation method due to the unpredictable effect of ignoring higher order terms,a modified parameter perturbation method is presented to predict the eigenvalue intervals of the uncertain structures with interval parameters.In the proposed method,interval variables are used to quantitatively describe all the uncertain parameters.Different order perturbations in both eigenvalues and eigenvectors are fully considered.By retaining higher order terms,the original dynamic eigenvalue equations are transformed into interval linear equations based on the orthogonality and regularization conditions of eigenvectors.The eigenvalue ranges and corresponding eigenvectors can be approximately predicted by the parameter combinatorial approach.Compared with the Monte Carlo method,two numerical examples are given to demonstrate the accuracy and efficiency of the proposed algorithm to solve both the real eigenvalue problem and complex eigenvalue problem.     

Radiative temperature dissipation in a finite atmosphere—I. The homogeneous case

We have developed a new approach toward solving problems of linear radiative relaxation of LTE temperature perturbations in a plane-parallel atmosphere of finite extent. We show that the mathematical problem is one of solving an integral eigenvalue equation, for which non-trivial solutions exist only for discrete values of the radiative relaxation time. The solutions for the spatial part of the perturbation constitute a complete and orthogonal set of basis functions, making it possible to solve more general problems of temperature relaxation. In applying this method to radiative relaxation in the middle atmosphere of earth, we show how the additional influences of photochemical coupling, advection by winds, and eddy diffusion by small-scale turbulence may be easily included using matrix perturbation techniques. We have solved the homogeneous integral equation for a wide variety of vertical thicknesses in an idealized homogeneous slab medium. Adopting a number of different analytic line profiles (rectangular, Doupler, Voigt, and Lorentz) we have obtained numerical solutions using an exponential-kernel method for solving the integral equation. The discrete eigenvalue “spectrum” is presented for vertical optical depths (0–10³) at line-center, and is used in solving several initial-value problems for a decaying temperature perturbation. We find that the eigenvalue spectrum is bounded from above by the lowest-order eigenvalue, and bounded from below by the familiar transparent approximation. The dependence of the lowest even eigenvalue on optical depth and the relative separation of the higher eigenvalues are found to depend sensitively on the line profile.     

Eigenvalues of structures with uncertain elastic boundary restraints

The eigenvalue problems of structures with random elastic boundary supports are studied in this paper. Using the perturbation method, the differential equations including stochastic distributed parameters and random boundary conditions that govern the eigenproblems are transformed to a series of deterministic differential equations and boundary conditions. Then the differential equations and boundary conditions are discretized utilizing the finite element method (FEM). This process is easy to be implemented since the resulting perturbation equations with different orders own the same FEM meshes. The first-order and second-order sensitivities of eigenvalues are derived through solving the deterministic algebraic equations. Upon determining these sensitivities of eigenvalues, the approximate statistic expressions of random eigenvalues considering uncertain elastic boundary supports are established. At the end, several numerical examples are given to illustrate the application and effectiveness of the present method. Comparison of the present numerical results with those from the Monte-Carlo simulation method verifies the accuracy of the developed method.     

Eigenvector approximation leading to exponential speedup of quantum eigenvalue calculation

We present an efficient method for preparing the initial state required by the eigenvalue approximation quantum algorithm of Abrams and Lloyd. Our method can be applied when solving continuous Hermitian eigenproblems, e.g., the Schr?dinger equation, on a discrete grid. We start with a classically obtained eigenvector for a problem discretized on a coarse grid, and we efficiently construct, quantum mechanically, an approximation of the same eigenvector on a fine grid. We use this approximation as the initial state for the eigenvalue estimation algorithm, and show the relationship between its success probability and the size of the coarse grid.     

一类新型激波捕捉格式的耗散性与稳定性分析

高超声速流动是高复杂性的可压缩黏性流动, 其中存在激波、剪切层、激波/激波干扰、激波/边界层干扰、旋涡与分离流动等复杂流场结构. 对其进行准确模拟需要使用低耗散、强鲁棒性的激波捕捉方法. 本文基于一类新型的通量项分裂方法, 提出了一种耗散低且鲁棒性好的激波捕捉格式K-CUSP-X. 对该格式的耗散性和激波稳定性进行了详细的理论分析, 得到了格式激波稳定的数值条件. 推论认为, 迎风格式激波稳定的充分条件为速度扰动量具有衰减性, 数值实验验证了该推论. 研究表明, 该格式与Toro提出的通量分裂格式K-CUSP-T相比, 在保证精确捕捉接触间断的同时, 又具有更好的稳定性, 在激波处不会产生“红玉”现象.     

Finite element reduced order model for noise and vibration reduction of double sandwich panels using shunted piezoelectric patches

This work concerns the control of sound transmission through double laminated panels with viscoelastic core using semi-passive piezoelectric shunt technique. More specifically, the system consists of two laminated walls, each one composed of three layers and called sandwich panel with an air cavity in between. The external sandwich panel has a surface-mounted piezoelectric patches. The piezoelectric elements, connected with resonant shunt circuits, are used for the vibration damping of some specific resonance frequencies of the coupled system. Firstly, a finite element formulation of the fully coupled visco-electro-mechanical-acoustic system is presented. This formulation takes into account the frequency dependence of the viscoelastic material. A modal reduction approach is then proposed to solve the problem at a lower cost. In the proposed technique, the coupled system is solved by projecting the mechanical displacement unknown on a truncated basis composed by the first real short-circuit structural normal modes and the pressure unknown on a truncated basis composed by the first acoustic modes with rigid boundaries conditions. The few initial electrical unknowns are kept in the reduced system. A static correction is also introduced in order to take into account the effect of higher modes. Various results are presented in order to validate and illustrate the efficiency of the proposed finite element reduced order formulation.     

Photonic crystal device sensitivity analysis with Wannier basis gradients

We present a powerful sensitivity analysis method for devices in a photonic crystal. The method is based on a Wannier basis field expansion and efficient matrix analysis techniques for finding eigenvalue and transmission gradients with respect to the perturbation. The method permits fast analysis of a large number of dielectric perturbation situations for multiple devices in a photonic crystal. We verify the method with finite-difference time-domain and plane-wave expansion calculations.     

Second-order small perturbation method for transmission from dielectric rough surfaces

Although the small perturbation method (SPM) for rough surface scattering has been studied extensively in problems in optics, remote sensing and propagation, there are fewer studies on rough surface transmission by the SPM. In this paper, from Huygens’ principle and the extinction theorem, the SPM is used to derive the transmitted field to the second order, and expressions for the bidirectional transmission coefficient and the total surface transmittance to the second order are developed for the random rough surface. The refined expressions can be applied to the situations where the transmission characteristics of a random rough surface need to be more accurately calculated. For example, to calculate the brightness temperature of stratified rough media by the incoherent method, we have to know the bidirectional transmission coefficient or transmittance of random rough surface as accurately as possible. The accuracy of the presented expressions is verified through the conservation of energy. It is shown that the transmission characteristics calculated by SPM to the first order violate conservation of energy, whereas solutions to the second order conform to energy conservation much better. This is particularly important for the calculation of transmittance or emissivity.     

A new formulation of the linear dynamic response to random excitations

With the damped harmonic oscillator as a model, a new formulation is presented of the stochastic response of linear dynamic systems under random excitation. The starting point of the method is the Fourier-Stieltjes representation for a general random process first introduced by Priestley. A unique feature of the present formulation is analytical compactness. This becomes a significant advantage over existing methods when one attempts to compute the covariance matrix elements which are often necessary for crossing problems.     

Iterative approach for the eigenvalue problems

An approximation method based on the iterative technique is developed within the framework of linear delta expansion (LDE) technique for the eigenvalues and eigenfunctions of the one-dimensional and three-dimensional realistic physical problems. This technique allows us to obtain the coefficient in the perturbation series for the eigenfunctions and the eigenvalues directly by knowing the eigenfunctions and the eigenvalues of the unperturbed problems in quantum mechanics. Examples are presented to support this. Hence, the LDE technique can be used for nonperturbative as well as perturbative systems to find approximate solutions of eigenvalue problems.     

NON-LINEAR FORCED VIBRATIONS OF PLATES BY AN ASYMPTOTIC-NUMERICAL METHOD

Non-linear forced vibrations of thin elastic plates have been investigated by an asymptotic-numerical method (ANM). Various types of harmonic excitation forces such as distributed and concentrated are considered. Using the harmonic balance method and Hamilton's principle, the equation of motion is converted into an operational formulation. Based on the finite element method a starting point corresponding to a non-linear solution associated to a given frequency and amplitude of excitation is computed. Applying perturbation techniques in the vicinity of this solution, the non-linear governing equation obtained is transformed into a sequence of linear problems having the same stiffness matrix. Employing one matrix inversion, a large number of terms of the perturbation series of the displacement and frequency can be easily computed with a small computation time. Iterations of this method lead to a powerful path-following technique. Comprehensive numerical tests for forced vibrations of plates subjected to time-harmonic lateral excitations are reported.     

Stochastic interval analysis of natural frequency and mode shape of structures with uncertainties

In this paper, natural frequencies and mode shapes of structures with mixed random and interval parameters are investigated by using a hybrid stochastic and interval approach. Expressions for the mean value and variance of natural frequencies and mode shapes are derived by using perturbation method and random interval moment method. The bounds of these probabilistic characteristics are then determined by interval arithmetic. Two examples are given first to illustrate the feasibility of the presented method and the results are verified by Monte Carlo Simulations. The presented approach is also applicable to solve pure random and pure interval problems. This capability is demonstrated in the third and fourth examples through the comparisons with the peer research outcomes.     

Derivatives with respect to parameters of the eigenvectors of a Hamiltonian and their application

A method is described for finding the first- and higher-order derivatives of the eigenvectors of a Hamiltonian with respect to its parameters. This method is useful even when the explicit dependence of the eigenvectors on the parameters is not known. The method is based on a transfer of the differentiation from the eigenvector to the Hamiltonian and on a separate analysis of the derivatives of the projections of the eigenvector onto the corresponding subspace and onto the orthogonal complement of this subspace. Conditions governing the position of the eigenvector being differentiated in its degenerate subspace are analyzed. This method can be used in certain fundamental problems, and it can be related to steady-state Rayleigh-Schrö-dinger perturbation theory.Translated from Izvestiya Vysshikh Uchebnykh Zave-denii Fizika, No. 3, pp. 74–78, March, 1971.