ON SMOOTH LU DECOMPOSITIONS WITH APPLICATIONS TO SOLUTIONS OF NONLINEAR EIGENVALUE PROBLEMS

We study the smooth LU decomposition of a given analytic functional A-matrix A（A） and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A（A）, and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.     

Bootstrap approximation of wavelet estimates in a semiparametric regression model

The inference for the parameters in a semiparametric regression model is studied by using the wavelet and the bootstrap methods. The bootstrap statistics are constructed by using Efron＇s resampling technique, and the strong uniform convergence of the bootstrap approximation is proved. Our results can be used to construct the large sample confidence intervals for the parameters of interest. A simulation study is conducted to evaluate the finite-sample performance of the bootstrap method and to compare it with the normal approximation-based method.     

Cross approximation in tensor electron density computations

We propose new tensor approximation algorithms for certain discrete functions related with Hartree–Fock/Kohn–Sham equations. Given a canonical tensor representation for the electron density function (for example, produced by quantum chemistry packages such as MOLPRO), we obtain its Tucker approximation with much fewer parameters than the input data and the Tucker approximation for the cubic root of this function, which is part of the Kohn–Sham exchange operator. The key idea is in the fast and accurate prefiltering of possibly large‐scale factors of the canonical tensor input. The new algorithms are based on the incomplete cross approximation method applied to matrices and tensors of order 3 and outperform other tools for the same purpose. First, we show that the cross approximation method is robust and much faster than the singular value decomposition‐based approach. As a consequence, it becomes possible to increase the resolution of grid and the complexity of molecules that can be handled by the Hartree–Fock chemical models. Second, we propose a new fast approximation method for f^1/3(x, y, z), based on the factor prefiltering method for f(x, y, z) and certain mimic approximation hypothesis. Third, we conclude that the Tucker format has advantages in the storage and computation time compared with the ubiquitous canonical format. Copyright © 2009 John Wiley & Sons, Ltd.     

A novel approach to solve linear systems arising from BEM using fast wavelet transforms

This paper uses Daubechies orthogonal wavelets to change dense and fully populated matrices of boundary element method (BEM) systems into sparse and semi‐banded matrices. Then a novel algorithm based on hierarchical nature of multiresolution analysis is introduced to solving resultant sparse linear systems. This algorithm decomposes NS‐form of transformed parent matrix into descendant systems with reduced sizes and solves them iteratively using GMRES algorithm. Both parts, changing dense matrices to sparse systems and the novel solver, can be added as a black box to the existing BEM codes. Transforming matrices into wavelet space needs less time than saved by solving sparse large systems. Numerical results with a precise study on sensitivity of solution for physical variables to the thresholding parameter, and savings in computer time and memory are presented. Also, the suitable value for thresholding parameter is recommended for elasticity problems. The results indicate that the proposed method is efficient for large problems. Copyright © 2009 John Wiley & Sons, Ltd.     

Global well‐posedness of the Cauchy problem for certain magnetohydrodynamic‐α models

This paper is devoted to study the Cauchy problem for certain incompressible magnetohydrodynamics‐α model. In the Sobolev space with fractional index s>1, we proved the local solutions for any initial data, and global solutions for small initial data. Furthermore, we also prove that as α→0, the MHD‐α model reduces to the MHD equations, and the solutions of the MHD‐α model converge to a pair of solutions for the MHD equations. Copyright © 2010 John Wiley & Sons, Ltd.     

Solving evolution equations using a new iterative method

In this article, we apply the new iterative method proposed by Daftardar‐Gejji and Jafari (J Math Anal Appl 316, (2006), 753–763) for solving various linear and nonlinear evolution equations. The results obtained are compared with the results by existing methods. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

基于Adaboost-SVR模型的我国碳排放强度分析与预测

首先对我国1960-2017年的碳排放趋势分5个阶段分析,发现虽然在不同时期存在波动,但长期来看,我国碳排放强度呈逐步下降趋势.然后对差分平稳后的序列数据建立Adaboost-SVR预测模型,采用RMSE、MAPE、MAE、MSE四个评价指标比较Adaboost-SVR模型与Adaboost-DT、SVR、BP神经网络对碳排放强度的预测精度.结果表明,组合模型明显优于其他3种模型,对于碳排放强度预测具有很高的可靠性.另外,通过使用Adaboost-SVR模型进行后续年份预测,发现我国未来碳排放强度总体将继续缓慢下降.最后,基于二氧化碳排放量的LMID分解结果,提出调整能源产业结构, 促进可再生能源利用等节能减排建议.     

梁方程时间依赖全局吸引子的存在性

研究了梁方程时间依赖吸引子的存在性,在非线性项f满足临界增长条件时,基于时间依赖全局吸引子的存在性定理,应用先验估计和算子分解方法验证了系数参数与时间t有关时,梁方程对应的过程族{U(t,τ)}的渐近紧性,从而得到梁方程时间依赖全局吸引子的存在性及正则性.     

Numerical Investigation of Sound Generation due to Laminar Flow Past Elliptic Cylinders

Numerical investigation of sound generation due to unsteady laminar flow past elliptic cylinders has been carried out using direct numerical simulation $(DNS)$ approach at a free-stream Mach number of $0.2$. Effects of aspect ratio $(0.6\le AR\le 1.0)$ and Reynolds number $(100\le Re \le 160)$ on the characteristics of radiated sound fields are analyzed. Two-dimensional compressible fluid flow equations are solved on a refined grid using high resolution dispersion relation preserving $(DRP)$ schemes. Using present $DNS$ data, equivalent noise sources as given by various acoustic analogies are evaluated. Amplitudes and frequencies associated with these noise sources are further related to characteristics of disturbance pressure fields. Disturbance pressure fields are intensified with increase in Reynolds number and aspect ratio. Thus, radiated sound power increases with increase in Reynolds number and aspect ratio. Among various cases studied here, minimum and maximum values of radiated sound power are found at $Re=120$ & $AR=0.6$ and $Re=160$ & $AR=1.0$, respectively. Directivity patterns show that the generated sound fields are dominated by the lift dipole for all cases. Next, proper orthogonal decomposition $(POD)$ technique has been implemented for decomposing disturbance pressure fields. The $POD$ modes associated with the lift and the drag dipoles have been identified. $POD$ analyses also clearly display that the radiated sound fields are dominated by the lift dipole only. Further, acoustic and hydrodynamic modes obtained using Doak's decomposition method have confirmed the patterns of radiated sound field intensities.