Modulations of Deep Water Waves and Spectral Filtering

Modulations of deep water waves are studied by a new formalism of spectral filtering. For single-mode dynamics, spectral filtering results in computable equations, which are counterpart to the nonlinear Schrödinger (NLS) equations. An essential feature of new equations is that bandwidth limitation is decoupled from small-amplitude assumption. The filtered equations have a substantially broader range of validity than the NLS equations, and may be viewed as intermediate between the NLS and Zakharov equations. The new single-mode equations reproduce exactly the conditions for nonlinear four-wave resonance ("figure 8" of Phillips [ 1 ]) even for bandwidths greater than unity. Sideband instability for uniform Stokes waves is limited to finite bandwidths only, and agrees well with exact results of McLean [ 2 ].     

Application of spectral filtering to discontinuous Galerkin methods on triangulations

We adapt the spectral viscosity (SV) formulation implemented as a modal filter to a discontinuous Galerkin (DG) method solving hyperbolic conservation laws on triangular grids. The connection between SV and spectral filtering, which is undertaken for the first time in the context of DG methods on unstructured grids, allows to specify conditions on the filter strength regarding time step choice and mesh refinement. A crucial advantage of this novel damping strategy is its low computational cost. We furthermore obtain new error bounds for filtered Dubiner expansions of smooth functions. While high order accuracy with respect to the polynomial degree N is proven for the filtering procedure in this case, an adaptive application is proposed to retain the high spatial approximation order. Although spectral filtering stabilizes the scheme, it leaves weaker oscillations. Therefore, as a postprocessing step, we apply the image processing technique of digital total variation (DTV) filtering in the new context of DG solutions and prove conservativity in the limit for this filtering procedure. Numerical experiments for scalar conservation laws confirm the designed order of accuracy of the DG scheme with adaptive modal filtering for polynomial degrees up to 8 and the viability of spectral and DTV filtering in case of shocks. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

A Green’s function model for the analysis of laser heating of materials

An analytical model based on Green’s function method is developed to analyze the temperature distribution and heated regions in a material irradiated by a high-energy laser beam. The model is multi-dimensional, transient and incorporates different types of beam characteristics and boundary conditions. The multi-dimensional integration formulas in the Green’s function solution equation are evaluated using an adaptive numerical integration algorithm. A parametric study is conducted to show the effect of various laser beam parameters and material properties on the laser heating process.     

VCSEL激光器仿真模型在家庭固定网络的应用

探讨了VCSEL激光器工作环境的温度问题和激光器的带宽问题,提出VCSEL激光器的光功率与电流(L-I)模型和器件工作时的电压与电流(U-I)特性曲线模型,利用最小二乘法和高斯牛顿法,对模型初步求解并改进,最终得到最优的激光器工作环境的温度.提出一种基于速率方程的VCSEL的带宽模型(小信号响应模型),并考虑激光器的温度和偏置电流对带宽的影响,利用非线性最优化的方法确定带宽模型,最终得到相应的结论.     

The Schrödinger operator: Perturbation determinants,the spectral shift function,trace identities,and all that

We discuss applications of the M. G. Krein theory of the spectral shift function to the multidimensional Schrödinger operator. Specific properties of this function, for example, its high-energy asymptotics are studied. Trace identities are derived.     

A trace formula and application to Stark Hamiltonians with nonconstant magnetic fields

In this paper we generalize a trace formula due to D. Robert involving the spectral shift function. We then give applications to the study of the high-energy and the semiclassical asymptotics of the spectral shift function for a Stark Hamiltonian with nonconstant magnetic field.     

Spectral Theory for a Mathematical Model of the Weak Interaction: The Decay of the Intermediate Vector Bosons W ��, II

We do the spectral analysis of the Hamiltonian for the weak leptonic decay of the gauge bosons W ^±. Using Mourre theory, it is shown that the spectrum between the unique ground state and the first threshold is purely absolutely continuous. Neither sharp neutrino high-energy cutoff nor infrared regularization is assumed.     

Filtering in Legendre spectral methods

We discuss the impact of modal filtering in Legendre spectral methods, both on accuracy and stability. For the former, we derive sufficient conditions on the filter to recover high order accuracy away from points of discontinuity. Computational results confirm that less strict necessary conditions appear to be adequate. We proceed to discuss a instability mechanism in polynomial spectral methods and prove that filtering suffices to ensure stability. The results are illustrated by computational experiments.

     

An adaptive spectral least-squares scheme for the Burgers equation

A least-squares spectral collocation method for the one-dimensional inviscid Burgers equation is proposed. This model problem shows the stability and high accuracy of these schemes for nonlinear hyperbolic scalar equations. Here we make use of a least-squares spectral approach which was already used in an earlier paper for discontinuous and singular perturbation problems (Heinrichs, J. Comput. Appl. Math. 157:329–345, 2003). The domain is decomposed in subintervals where continuity is enforced at the interfaces. Equal order polynomials are used on all subdomains. For the spectral collocation scheme Chebyshev polynomials are employed which allow the efficient implementation with Fast Fourier Transforms (FFTs). The collocation conditions and the interface conditions lead to an overdetermined system which can be efficiently solved by least-squares. The solution technique will only involve symmetric positive definite linear systems. The scheme exhibits exponential convergence where the exact solution is smooth. In parts of the domain where the solution contains discontinuities (shocks) the spectral solution displays a Gibbs-like behavior. Here this is overcome by some suitable exponential filtering at each time level. Here we observe that by over-collocation the results remain stable also for increasing filter parameters and also without filtering. Furthermore by an adaptive grid refinement we were able to locate the precise position of the discontinuity. Numerical simulations confirm the high accuracy of our spectral least-squares scheme.      

非线性守恒方程的两种稳定化谱元方法

考虑非线性守恒方程的高阶数值解法,介绍了基于谱元法的两种稳定性方法,一种是谱粘性消去法(SVV),另一种是过滤法.在SVV方法中,我们推广并分析了传统的基于单区域的SVV算子的定义.在过滤法中,我们分析了SVV-H elm holtz过滤算子的性质.文中从分析和计算两方面对两种方法进行了比较,建立了两者之间的关系.最后通过一系列数值试验说明方法的有效性.     

A projective method for the construction of solutions of the problem of normal symmetric oscillations of viscous liquid

We propose a variational formulation of the spectral problem of normal symmetric oscillations of viscous liquid. On the basis of this formulation, we construct a projective method for the determination of real eigenvalues of the problem. We present the numerical realization of this method in the case of a spherical cavity.     

New GMRES(k)-Type Algorithms with Explicit Restarts and the Analysis of Their Convergence Properties Based on Matrix Relations in QR Form

The paper considers the problem of constructing a basic iterative scheme for solving systems of linear algebraic equations with unsymmetric and indefinite coefficient matrices. A new GMRES-type algorithm with explicit restarts is suggested. When restarting, this algorithm takes into account the spectral/singular data transferred using orthogonal matrix relations in the so-called QR form, which arise when performing inner iterations of Arnoldi type. The main idea of the algorithm developed is to organize inner iterations and the filtering of directions before restarting in such a way that, from one restart to another, matrix relations effectively accumulate information concerning the current approximate solution and, simultaneously, spectral/singular data, which allow the algorithm to converge with a rate comparable with that of the GMRES algorithm without restarts. Convergence theory is provided for the case of nonsingular, unsymmetric, and indefinite matrices. A bound for the rate of decrease of the residual in the course of inner Arnoldi-type iterations is obtained. This bound depends on the spectral/singular characterization of the subspace spanned by the directions retained upon filtering and is used in developing efficient filtering procedures. Numerical results are provided. Bibliography: 9 titles.     

Data-Adaptive Estimation of Time-Varying Spectral Densities

This article introduces a data-adaptive nonparametric approach for the estimation of time-varying spectral densities from nonstationary time series. Time-varying spectral densities are commonly estimated by local kernel smoothing. The performance of these nonparametric estimators, however, depends crucially on the smoothing bandwidths that need to be specified in both time and frequency direction. As an alternative and extension to traditional bandwidth selection methods, we propose an iterative algorithm for constructing localized smoothing kernels data-adaptively. The main idea, inspired by the concept of propagation-separation, is to determine for a point in the time-frequency plane the largest local vicinity over which smoothing is justified by the data. By shaping the smoothing kernels nonparametrically, our method not only avoids the problem of bandwidth selection in the strict sense but also becomes more flexible. It not only adapts to changing curvature in smoothly varying spectra but also adjusts for structural breaks in the time-varying spectrum. Supplementary materials, including the R package tvspecAdapt containing an implementation of the routine, are available online.     

Effective bandwidths: Call admission,traffic policing and filtering for ATM networks

In this paper we review and extend the effective bandwidth results of Kelly [28], and Kesidis, Walrand and Chang [29, 6]. These results provide a framework for call admission schemes which are sensitive to constraints on the mean delay or the tail distribution of the workload in buffered queues. We present results which are valid for a wide variety of traffic streams and discuss their applicability for traffic management in ATM networks. We discuss the impact of traffic policing schemes, such as thresholding and filtering, on the effective bandwidth of sources. Finally we discuss effective bandwidth results for Brownian traffic models for which explicit results reveal the interaction arising in finite buffers.     

Resolvent and spectral measure on non-trapping asymptotically hyperbolic manifolds I: Resolvent construction at high energy

This is the first in a series of papers in which we investigate the resolvent and spectral measure on non-trapping asymptotically hyperbolic manifolds with applications to the restriction theorem, spectral multiplier results and Strichartz estimates. In this first paper, we construct the high energy resolvent on general non-trapping asymptotically hyperbolic manifolds, using semiclassical Lagrangian distributions and semiclassical intersecting Lagrangian distributions, together with the 0-calculus of Mazzeo-Melrose.

Our results generalize recent work of Melrose, Sá Barreto and Vasy, which applies to metrics close to the exact hyperbolic metric. We note that there is an independent work by Y. Wang which also constructs the high-energy resolvent.     

Higher-order accurate polyspectral estimation with flat-top lag-windows

Improved performance in higher-order spectral density estimation is achieved using a general class of infinite-order kernels. These estimates are asymptotically less biased but with the same order of variance as compared to the classical estimators with second-order kernels. A simple, data-dependent algorithm for selecting the bandwidth is introduced and is shown to be consistent with estimating the optimal bandwidth. The combination of the specialized family of kernels with the new bandwidth selection algorithm yields a considerably improved polyspectral estimator surpassing the performances of existing estimators using second-order kernels. Bispectral simulations with several standard models are used to demonstrate the enhanced performance with the proposed methodology.     

Solving Systems of Linear Equations with Normal Coefficient Matrices and the Degree of the Minimal Polyanalytic Polynomial

The generalized Lanczos process applied to a normal matrix A builds up a condensed form of A, which can be described as a band matrix with slowly growing bandwidth. For certain classes of normal matrices, the bandwidth turns out to be constant. It is shown that, in such cases, the bandwidth is determined by the degree of the minimal polyanalytic polynomial of A. It was in relation to the generalized Lanczos process thatM.Huhtanen introduced the concept of the minimal polyanalytic polynomial of a normal matrix.     

A four-wing chaotic attractor generated from a new 3-D quadratic autonomous system

This paper introduces a new 3-D quadratic autonomous system, which can generate two coexisting single-wing chaotic attractors and a pair of diagonal double-wing chaotic attractors. More importantly, the system can generate a four-wing chaotic attractor with very complicated topological structures over a large range of parameters. Some basic dynamical behaviors and the compound structure of the new 3-D system are investigated. Detailed bifurcation analysis illustrates the evolution processes of the system among two coexisting sinks, two coexisting periodic orbits, two coexisting single-wing chaotic attractors, major and minor diagonal double-wing chaotic attractors, and a four-wing chaotic attractor. Poincaré-map analysis shows that the system has extremely rich dynamics. The physical existence of the four-wing chaotic attractor is verified by an electronic circuit. Finally, spectral analysis shows that the system has an extremely broad frequency bandwidth, which is very desirable for engineering applications such as secure communications.     

Asymptotic expansion for harmonic functions in the half‐space with a pressurized cavity

In this paper, we address a simplified version of a problem arising from volcanology. Specifically, as a reduced form of the boundary value problem for the Lamé system, we consider a Neumann problem for harmonic functions in the half‐space with a cavity C. Zero normal derivative is assumed at the boundary of the half‐space; differently, at ?C, the normal derivative of the function is required to be given by an external datum g, corresponding to a pressure term exerted on the medium at ?C. Under the assumption that the (pressurized) cavity is small with respect to the distance from the boundary of the half‐space, we establish an asymptotic formula for the solution of the problem. Main ingredients are integral equation formulations of the harmonic solution of the Neumann problem and a spectral analysis of the integral operators involved in the problem. In the special case of a datum g, which describes a constant pressure at ?C, we recover a simplified representation based on a polarization tensor. Copyright © 2016 John Wiley & Sons, Ltd.