Hamilton-Jacobi方程的小波Galerkin方法

本文选择Daubechies小波尺度函数空间作为Galerkin方法的测试函数空间,并将其应用于Hamilton-Jacobi方程,得到了求解Hamilton-Jacobi方程的小波Galerkin方法的数值格式．由于小波在时间和频率上的局部性,本方法适用于处理具有奇异解的问题,可以有效地防止数值振荡．数值试验显示,本方法是有效的．     

A wavelet-Galerkin method for high order numerical differentiation

Numerical differentiation is a classical ill-posed problem. In this paper, we propose a wavelet-Galerkin method for high order numerical differentiation. By an appropriate choice of the regularization parameter an order optimal stability estimate of Hölder type is obtained. Some numerical examples show that the method is effective and stable.     

Wavelet-Galerkin method for solving parabolic equations in finite domains

A novel wavelet-Galerkin method tailored to solve parabolic equations in finite domains is presented. The emphasis of the paper is on the development of the discretization formulations that are specific to finite domain parabolic equations with arbitrary boundary conditions based on weak form functionals. The proposed method also deals with the development of algorithms for computing the associated connection coefficients at arbitrary points. Here the Lagrange multiplier method is used to enforce the essential boundary conditions. The numerical results on a two-dimensional transient heat conducting problem are used to validate the proposed wavelet-Galerkin algorithm as an effective numerical method to solve finite domain parabolic equations.     

The weakly nonlocal limit of a one-population Wilson-Cowan model

We derive the weakly nonlocal limit of a one-population neuronal field model of the Wilson-Cowan type in one spatial dimension. By transforming this equation to an equation in the firing rate variable, it is shown that stationary periodic solutions exist by appealing to a pseudo-potential analysis. The solutions of the full nonlocal equation obey a uniform bound, and the stationary periodic solutions in the weakly nonlocal limit satisfying the same uniform bound are characterized by finite ranges of pseudo energy constants. The time dependent version of the model is reformulated as a Ginzburg-Landau-Khalatnikov type of equation in the firing rate variable where the maximum (minimum) points correspond stable (unstable) homogeneous solutions of the weakly nonlocal limit. Based on this formulation it is also conjectured that the stationary periodic solutions are unstable. We implement a numerical method for the weakly nonlocal limit of the Wilson-Cowan type of model based on the wavelet-Galerkin approach. We perform some numerical tests to illustrate the stability of homogeneous solutions and the evolution of the bumps.     

基于小波-Galerkin方法的结构随机动力响应功率谱的确定

在最近发展的周期广义谐和小波PGHW(Periodic Generalized Harmonic Wavelet)的基础上,通过小波-Galerkin方法推导得到了线性单自由度结构的随机动力响应功率谱密度。在此过程中,利用PGHW的解析形式及其在频域内的特殊性:(1)推导得出了PGHW的联系系数(Connection Coefficient)的解析形式;(2)基于PGHW及其联系系数,利用小波-Galerkin方法推导得到了线性单自由度系统在确定性激励下的响应;(3)得到了在具有演变功率谱的随机动力激励下单自由度线性振子的随机响应功率谱解答。数值算例表明,无论是确定性响应解答,还是随机动力响应的功率谱密度,小波-Galerkin法的计算结果均能较好地吻合数值解。     

Wavelet Solution of Plane Elasticity Problem in the Upper Half-plane

The plane elasticity problem includes plane strain problem and plane stress problem which are widely applied in mechanics and engineering. In this article, we first reduce the plane elasticity problem in the upper half-plane into natural boundary integral equation and then apply wavelet-Galerkin method to deal with the numerical solution of the natural boundary integral equation. The test and trial functions used are the scaling basis functions of Shannon wavelet. In our fast algorithm, the computational formulae of entries of the stiffness matrix yield simple close-form and only 3 K entries need to be computed for one 4 K ‐ 4 K stiffness matrix.     

Modeling vortex-shedding effects for the stochastic response of tall buildings in non-synoptic winds

This study derives a model for the vortex-induced vibration and the stochastic response of a tall building in strong non-synoptic wind regimes. The vortex-induced stochastic dynamics is obtained by combining turbulent-induced buffeting force, aeroelastic force and vortex-induced force. The governing equations of motion in non-synoptic winds account for the coupled motion with nonlinear aerodynamic damping and non-stationary wind loading. An engineering model, replicating the features of thunderstorm downbursts, is employed to simulate strong non-synoptic winds and non-stationary wind loading. This study also aims to examine the effectiveness of the wavelet-Galerkin (WG) approximation method to numerically solve the vortex-induced stochastic dynamics of a tall building with complex wind loading and coupled equations of motions. In the WG approximation method, the compactly supported Daubechies wavelets are used as orthonormal basis functions for the Galerkin projection, which transforms the time-dependent coupled, nonlinear, non-stationary stochastic dynamic equations into random algebraic equations in the wavelet space. An equivalent single-degree-of-freedom building model and a multi-degree-of-freedom model of the benchmark Commonwealth Advisory Aeronautical Research Council (CAARC) tall building are employed for the formulation and numerical analyses. Preliminary parametric investigations on the vortex-shedding effects and the stochastic dynamics of the two building models in non-synoptic downburst winds are discussed. The proposed WG approximation method proves to be very powerful and promising to approximately solve various cases of stochastic dynamics and the associated equations of motion accounting for vortex shedding effects, complex wind loads, coupling, nonlinearity and non-stationarity.     

一类生物模型方程的数值解及其分析

本文利用小波 -Galerkin方法 ,求解生物领域的一模型方程 :一类带有小位移的线性二阶微分 -差分方程 (DDE) ,并对解的边界层性质进行数值探讨。结果表明 ,当小位移增加但仍保持很小时 ,解的边界层结构发生改变 ,甚至遭破坏     

Wavelet based schemes for linear advection-dispersion equation

In this paper, two wavelet based adaptive solvers are developed for linear advection-dispersion equation. The localization properties and multilevel structure of the wavelets in the physical space are used for adaptive computational methods for solution of equation which exhibit both smooth and shock-like behaviour. The first framework is based on wavelet-Galerkin and the second is based on multiscale decomposition of finite element method. Coiflet wavelet filter is incorporated in both the methods. The main advantage of both the adaptive methods is the elimination of spurious oscillations at very high Peclet number.