Linearly implicit methods for nonlinear PDEs with linear dispersion and dissipation

The linear stability of IMEX (IMplicit–EXplicit) methods and exponential integrators for stiff systems of ODEs arising in the discrete solution of PDEs is examined for nonlinear PDEs with both linear dispersion and dissipation, and a clear method of visualization of the linear stability regions is proposed. Predictions are made based on these visualizations and are supported by a series of experiments on five PDEs including quasigeostrophic equations and stratified Boussinesq equations. The experiments, involving 24 IMEX and exponential methods of third and fourth order, confirm the predictions of the linear stability analysis, that the methods are typically limited by small eigenvalues of the linear term and by eigenvalues on or near the imaginary axis rather than by large eigenvalues near the negative real axis. The experiments also demonstrate that IMEX methods achieve comparable stability to exponential methods, and that exponential methods are significantly more accurate only when the problem is nearly linear. Novel IMEX predictor–corrector methods are also derived.     

Poincaré Map Based on Splitting Methods

Firstly, by using the Liouville formula, we prove that the Jacobian matrix determinants of splitting methods are equal to that of the exact flow. However, for the explicit Runge-Kutta methods, there is an error term of order p ＋ I for the Jacobian matrix determinants. Then, the volume evolution law of a given region in phase space is discussed for different algorithms. It is proved that splitting methods can exactly preserve the sum of Lyapunov exponents invariable. Finally, a Poincaré map and its energy distribution of the Duffing equation are computed by using the second-order splitting method and the Heun method （a second-order Runge-Kutta method）. Computation illustrates that the results by splitting methods can properly represent systems＇ chaotic phenomena.     

Numerical methods for the simulation of a corrosion model with moving oxide layer

In this paper, we design numerical methods for a PDE system arising in corrosion modeling. This system describes the evolution of a dense oxide layer. It is based on a drift–diffusion system and includes moving boundary equations. The choice of the numerical methods is justified by a stability analysis and by the study of their numerical performance. Finally, numerical experiments with real-life data shows the efficiency of the developed methods.     

二维三温能量方程组的高效代数解法

针对二维三温能量方程九点格式离散后形成的非线性方程组,研制了高效求解的代数解法器.主要思想是在部分Newton-Krylov(PNK)方法和Jacobi矩阵自由的Newton-Krylov(JFNK)方法的框架下,结合非精确Newton类方法和预条件Krylov子空间方法进行高效求解.数值结果显示,PNK方法比非线性块Gauss-Seidel方法快6倍以上,在PNK框架下比较了3种预条件子和4种Krylov子空间方法,得出不同组合的最佳方案.还比较了JFNK方法和PNK方法.     

Suppression of acoustic noise in the measurements of wall pressure fluctuations

It is shown that, in hydrodynamic noise measurements in the presence of acoustic noise acting upon the pressure fluctuation receiver, spatial filtering methods should provide the best results. Active methods are developed for suppressing the acoustic noise that affects a miniature receiver in the course of turbulent pressure fluctuation measurements. The methods are based on complicating the structure of the measuring transducer by introducing an extra compensating sensing element whose characteristics are identical with those of the main sensing element. The spatial filtering of small-scale turbulent pressure fluctuations by a finite-size electroacoustic transducer is used as the basis for the development of noise-compensated measuring systems, as well as methods of measuring the turbulent pressure fluctuations by receivers with noise compensation. A numerical study of the wave-number filtering of acoustic noise in wall pressure measurements by a noise-compensated receiver is performed.     

Direct computation of the noise radiated by a subsonic cavity flow and application of integral methods

The goal of this paper is to investigate the acoustic field generated by the flow over a cavity using two different and complementary numerical methods. First, a Direct Numerical Simulation of the 2-D compressible Navier-Stokes equations is performed to obtain directly the radiated noise. The results of the acoustic and aerodynamic fields are compared to the experimental data in the literature. Second, this reference solution is compared to solutions provided by hybrid methods using the flowfield computed inside the cavity combined with an integral formulation to evaluate the far-field noise. Numerical issues of three integral methods are studied: the Ffowcs Williams and Hawkings analogy that extends Lighthill's theory to account for solid boundaries and two Wave Extrapolation Methods from a control surface, the Kirchhoff and porous Ffowcs Williams and Hawkings methods. All methods show a good agreement with the Direct Numerical Simulation, but the first one is more expensive owing to an additional volume integral. However, the analogy can help in the analysis of wave patterns, by separating the direct waves from the reflected ones. The wave extrapolation methods from a surface are more efficient and provide a complementary tool to extend Computational Aeroacoustics near field to the very far field.     

Scattering of gaussian beam by a ferrite cylinder

The scattering of a Gaussian beam by a ferrite cylinder is analyzed using two methods. The spatial pictures in the amplitude of the total Gaussian beam near the cylinder resulting from the two methods are both clarified and compared. The scattering by a ferrite cylinder coating a dielectric sleeve is also presented.     

Synchronization in a unified fractional-order chaotic system

In this paper, the synchronization in a unified fractional-order chaotic system is investigated by two methods. One is the frequency-domain method that is analysed by using the Laplace transform theory. The other is the time-domain method that is analysed by using the Lyapunov stability theory. Finally, the numerical simulations are used to illustrate the effectiveness of the proposed synchronization methods.     

Remote sensing by open-path FTIR spectroscopy. Comparison of different analysis techniques applied to ozone and carbon monoxide detection

The analysis of FTIR spectra is usually carried out by commercial programs that use methods based on classical least-squares (CLS) procedures These procedures are normally appropriate to analyse O-P FTIR spectra, but some occasional discrepancies with standard extractive methods have been observed.In this work, a line-by-line method (SFIT) is additionally used. Our purpose is to explain the sporadic discrepancies observed between CLS and standard extractive method results and to study the capability of both methods (CLS and SFIT) to analyse open-path measurements. This capability has been studied for ozone and carbon monoxide.It has been demonstrated that in CLS-based methods discrepancies are induced by the experimental background reference spectrum. In these cases, SFIT results are in very good agreement with the standard extractive methods.     

光学材料折射率和色散的高精度测量技术研究

本文探讨了在一般精度的测角仪器上达到高精度测量折射率和色散的新原理和新方法,这些新原理和方法包括:光学角规补偿法测色散、直角照射法、自准封闭式最小偏向角法。通过实际大量的测量工作证明在1″测角仪(测角精度为3″)和测微望远镜上,上述原理方法都可以达到△n~±3×10~(-6),△(n_λ_1-n_λ_2)~±2×10~(-6)的精度。     

A novel knowledge-based potential for RNA 3D structure evaluation

Ribonucleic acids(RNAs) play a vital role in biology, and knowledge of their three-dimensional(3D) structure is required to understand their biological functions. Recently structural prediction methods have been developed to address this issue, but a series of RNA 3D structures are generally predicted by most existing methods. Therefore, the evaluation of the predicted structures is generally indispensable. Although several methods have been proposed to assess RNA 3D structures, the existing methods are not precise enough. In this work, a new all-atom knowledge-based potential is developed for more accurately evaluating RNA 3D structures. The potential not only includes local and nonlocal interactions but also fully considers the specificity of each RNA by introducing a retraining mechanism. Based on extensive test sets generated from independent methods, the proposed potential correctly distinguished the native state and ranked near-native conformations to effectively select the best. Furthermore, the proposed potential precisely captured RNA structural features such as base-stacking and base-pairing. Comparisons with existing potential methods show that the proposed potential is very reliable and accurate in RNA 3D structure evaluation.     

Analysis of speech-based Speech Transmission Index methods with implications for nonlinear operations

The Speech Transmission Index (STI) is a physical metric that is well correlated with the intelligibility of speech degraded by additive noise and reverberation. The traditional STI uses modulated noise as a probe signal and is valid for assessing degradations that result from linear operations on the speech signal. Researchers have attempted to extend the STI to predict the intelligibility of nonlinearly processed speech by proposing variations that use speech as a probe signal. This work considers four previously proposed speech-based STI methods and four novel methods, studied under conditions of additive noise, reverberation, and two nonlinear operations (envelope thresholding and spectral subtraction). Analyzing intermediate metrics in the STI calculation reveals why some methods fail for nonlinear operations. Results indicate that none of the previously proposed methods is adequate for all of the conditions considered, while four proposed methods produce qualitatively reasonable results and warrant further study. The discussion considers the relevance of this work to predicting the intelligibility of cochlear-implant processed speech.     

Crack Tip Opening Angle as a Fracture Resistance Parameter to Describe Ductile Crack Extension and Arrest in Steel Pipes under Service Pressure

The angle between two element sides representing the crack tip is defined as the crack tip opening angle (CTOA). Its critical value is used as a criterion of fracture resistance for characterizing stable tearing in thin metallic materials. Various methods are used for determination of the CTOA. Optical microscopy is one of the most common methods as well as fitting of experimental load-displacement diagrams by the finite element method (DIC). Additionally, analytical analysis using the experimental load-displacement curve method (SSM) derived from the plastic hinge model of deflection in three-point bending of a ductile specimen is applied. This approach assumes a constant rotation centre distance. Values of CTOA for API 5L X65 pipe steel found by three methods—DIC, CNM, and SSM—are given. Values of CTOA given by these three methods are similar and close to 20°. A discussion on the different parameters used to characterize the fracture resistance of running cracks in a pipe under service pressure is presented. The energy of fracture at impact determined by Charpy or drop-weight tear test (DWTT) tests and the critical J energy parameter are considered as well as the yield locus after damage, cohesive zone energy, and CTOA is another approach. One notes that CTOA is assumed to be constant during stable crack extension and decreases linearly with crack length during the instable and primary phase. A numerical technique to describe a ductile running crack using the node release technique and using CTOA as the fracture resistance criterion is presented. This method is compared with three different two-curve methods (TCMs): the Battelle, high strength line pipe (HLP), and HLP-Sumitomo methods. The Batelle TCM, as the oldest method, based on Charpy energy, gives a strongly conservative prediction. Predictions by the CTOA method are close to those obtained by the HLP-Sumitomo one.     

On high-order shock-fitting and front-tracking schemes for numerical simulation of shock–disturbance interactions

High-order methods that can resolve interactions of flow-disturbances with shock waves are critical for reliable numerical simulation of shock wave and turbulence interaction. Such problems are not well understood due to the limitations of numerical methods. Most of the popular shock-capturing methods are only first-order accurate at the shock and may incur spurious numerical oscillations near the shock. Shock-fitting algorithms have been proposed as an alternative which can achieve uniform high-order accuracy and can avoid possible spurious oscillations incurred in shock-capturing methods by treating shocks as sharp interfaces. We explore two ways for shock-fitting: conventional moving grid set-up and a new fixed grid set-up with front tracking. In the conventional shock-fitting method, a moving grid is fitted to the shock whereas in the newly developed fixed grid set-up the shock front is tracked using Lagrangian points and is free to move across the underlying fixed grid. Different implementations of shock-fitting methods have been published in the literature. However, uniform high-order accuracy of various shock-fitting methods has not been systematically established. In this paper, we carry out a rigorous grid-convergence analysis on different variations of shock-fitting methods with both moving and fixed grids. These shock-fitting methods consist of different combinations of numerical methods for computing flow away from the shock and those for computing the shock movement. Specifically, we consider fifth-order upwind finite-difference scheme and shock-capturing WENO schemes with conventional shock-fitting and show that a fifth-order convergence is indeed achieved for a canonical one-dimensional shock-entropy wave interaction problem. We also show that the method of finding shock velocity from one characteristic relation and Rankine–Hugoniot jump condition performs better than the other methods of computing shock velocities. A high-order front-tracking implementation of shock-fitting is also presented in this paper and nominal rate of convergence is shown. The front-tracking results are validated by comparing to results from the conventional shock-fitting method and a linear-interaction analysis for a two-dimensional shock disturbance interaction problem.     

光学形式量子化算符法与矩阵法的对应关系

本文从光学位置坐标算符和正则动量算符的定义以及线性光学系统光学算符的幺正性假设出发,由线性系统的普遍性质建立光学系统算符表示和矩阵表示的关联方程,进而利用该方程实现了光学算符和光学矩阵的互为导出,因此证明了光学算符法和矩阵法的等效对应关系。最后,本文说明了光学算符的物理意义。 关键词：     

A comparative study of the effective dielectric constant of a medium containing randomly distributed dielectric spheres embedded in a homogeneous background medium

The effective dielectric constant ∈_eff of a medium containing randomly distributed dielectric particles has been analysed by conventional methods: Foldy's approximation, the quasi-crystalline approximation (QCA) and the QCA with coherent potential. These conventional methods, however, have been indicated to become invalid for particles with a high dielectric constant; we have thus presented a new method that is valid for them. This paper compares ∈_eff of our method with those of the conventional methods by changing the volume fraction and the dielectric constant of spheres. As a result, our method is shown to be more powerful for the analysis of ∈_eff than the conventional methods.     

上升管中强烈段塞流消除方法试验研究

为了消除上升管中强烈段塞流,采用节流和扰动两种方法进行了试验研究。两种方法都能完全消除强烈段塞流, 能使强烈段塞流的发生范围显著减小。节流程度较小时将导致强烈段塞流压力波动更大、周期更长;但扰动法能使压力波动的幅度变小、周期变短。两种方法都使上游管道的压力显著增大,但扰动法产生的背压比节流法产生的背压小。     

Comparative analysis of side-illumination methods of measurement of emission decay of fluorescent polymer fibers

The fluorescence decay constant of a fluorescent polymer optical fiber is measured by two different side illumination methods. Comparative analysis of these methods is performed. It is shown that the difference in the measured parameters is caused by different conditions of energy exchange between the core, cladding, and radiation modes of the fiber waveguide.     

Krylov subspace iterative methods for boundary element method based near-field acoustic holography

The reconstruction of the acoustic field for general surfaces is obtained from the solution of a matrix system that results from a boundary integral equation discretized using boundary element methods. The solution to the resultant matrix system is obtained using iterative regularization methods that counteract the effect of noise on the measurements. These methods will not require the calculation of the singular value decomposition, which can be expensive when the matrix system is considerably large. Krylov subspace methods are iterative methods that have the phenomena known as "semi-convergence," i.e., the optimal regularization solution is obtained after a few iterations. If the iteration is not stopped, the method converges to a solution that generally is totally corrupted by errors on the measurements. For these methods the number of iterations play the role of the regularization parameter. We will focus our attention to the study of the regularizing properties from the Krylov subspace methods like conjugate gradients, least squares QR and the recently proposed Hybrid method. A discussion and comparison of the available stopping rules will be included. A vibrating plate is considered as an example to validate our results.     

A class of hybrid DG/FV methods for conservation laws I: Basic formulation and one-dimensional systems

By comparing the discontinuous Galerkin (DG) and the finite volume (FV) methods, a concept of ‘static reconstruction’ and ‘dynamic reconstruction’ is introduced for high-order numerical methods. Based on the new concept, a class of hybrid DG/FV schemes is presented for one-dimensional conservation law using a ‘hybrid reconstruction’ approach. In the hybrid DG/FV schemes, the lower-order derivatives of a piecewise polynomial solution are computed locally in a cell by the DG method based on Taylor basis functions (called as ‘dynamic reconstruction’), while the higher-order derivatives are re-constructed by the ‘static reconstruction’ of the FV method, using the known lower-order derivatives in the cell itself and its adjacent neighboring cells. The hybrid DG/FV methods can greatly reduce CPU time and memory required by the traditional DG methods with the same order of accuracy on the same mesh, and they can be extended directly to unstructured and hybrid grids in two and three dimensions similar to the DG and/or FV methods. The hybrid DG/FV methods are applied to one-dimensional conservation law, including linear and non-linear scalar equation and Euler equations. In order to capture the strong shock waves without spurious oscillations, a simple shock detection approach is developed to mark ‘trouble cells’, and a moment limiter is adopted for higher-order schemes. The numerical results demonstrate the accuracy, and the super-convergence property is shown for the third-order hybrid DG/FV schemes. In addition, by analyzing the eigenvalues of the semi-discretized system in one dimension, we discuss the spectral properties of the hybrid DG/FV schemes to explain the super-convergence phenomenon.