Application of Fourier Series Expansion Method with PMLs to the Microcavities on Two-Dimensional Photonic Crystals

By using a Fourier series expansion method combined with Chew＇s perfectly matched layers （PMLs）, we analyze the frequency and quality factor of a micro-cavity on a two-dimensional photonic crystal is analyzed. Compared with the results by the method without PML and finite-difference time-domain （FDTD） based on supercell approximation, it can be shown that by the present method with PMLs, the resonant frequency and the quality factor values can be calculated satisfyingly and the characteristics of the micro-cavity can be obtained by changing the size and permittivity of the point defect in the micro-cavity.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

Multi-Symplectic Splitting Method for Two-Dimensional Nonlinear Schrodinger Equation

Using the idea of splitting numerical methods and the multi-symplectic methods, we propose a multisymplectic splitting （MSS） method to solve the two-dimensional nonlinear Schrodinger equation （2D-NLSE） in this paper. It is further shown that the method constructed in this way preserve the global symplectieity exactly. Numerical experiments for the plane wave solution and singular solution of the 2D-NLSE show the accuracy and effectiveness of the proposed method.

Experimental Consideration of Two-Dimensional Fourier Transform Spectroscopy

Two-dimensional Fourier transform(2D FT) spectroscopy is an important technology that developed in recent decades and has many advantages over other ultrafast spectroscopy methods. Although 2D FT spectroscopy provides great opportunities for studying various complex systems, the experimental implementation and theoretical description of 2D FT spectroscopy measurement still face many challenges, which limits their wide application.Recently, the 2D FT spectroscopy reaches maturity due to many new developments which greatly reduces the technical barrier in the experimental implementation of the 2D FT spectrometer. There have been several different approaches developed for the optical design of the 2D FT spectrometer, each with its own advantages and limitations. Thus, a procedure to help an experimentalist to build a 2D FT spectroscopy experimental apparatus is needed.This tutorial review is intending to provide an accessible introduction for a beginner to build a 2D FT spectrometer.

Comparison of different inversion methods of D-T_(2)two-dimensional nuclear magnetic resonance logging and applicability analysis

D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization methods)have been proposed successively,but most are limited to numerical simulations.This study focused on the applicability of different inversion methods for NMR logging data of various acquisition sequences,from which the optimal inversion method was selected based on the comparative analysis.First,the two-dimensional NMR logging principle was studied.Then,these inversion methods were studied in detail,and the precision and computational efficiency of CPMG and diffusion editing(DE)sequences obtained from oil-water and gas-water models were compared,respectively.The inversion results and calculation time of truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization were compared and analyzed through numerical simulations.The inversion method was optimized to process SP mode logging data from the MR Scanner instrument.The results showed that the TIST-regularization and LM-norm smoothing methods were more accurate for the CPMG and DE sequence echo trains of the oil-water and gas-water models.However,the LM-norm smoothing method was less time-consuming,making it more suitable for logging data processing.A case study in well A25 showed that the processing results by the LM-norm smoothing method were consistent with GEOLOG software.This demonstrates that the LM-norm smoothing method is applicable in practical NMR logging processing.

Eigenfunction expansion method and its application to two-dimensional elasticity problems based on stress formulation

This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.

A simplified two-dimensional boundary element method with arbitrary uniform mean flow

To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.

THE CHEBYSHEV PSEUDOSPECTRAL DOMAIN DECOMPO SITION METHOD FOR SOLVING TWO-DIMENSIONAL ELLIPTIC EQUATION

This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stability and convergence of this new scheme.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

Solution of two-dimensional scattering problem in piezoelectric/piezomagnetic media using a polarization method

Using a polarization method, the scattering problem for a two-dimensional inclusion embedded in infinite piezoelectric/piezomagnetic matrices is investigated. To achieve the purpose, the polarization method for a two-dimensional piezoelectric/piezomagnetic ＂comparison body＂ is formulated. For simple harmonic motion, kernel of the polarization method reduces to a 2-D time-harmonic Green＇s function, which is obtained using the Radon transform. The expression is further simplified under conditions of low frequency of the incident wave and small diameter of the inclusion. Some analytical expressions are obtained. The analytical solutions for generalized piezoelectric/piezomagnetic anisotropic composites are given followed by simplified results for piezoelectric composites. Based on the latter results, two numerical results are provided for an elliptical cylindrical inclusion in a PZT-5H-matrix, showing the effect of different factors including size, shape, material properties, and piezoelectricity on the scattering cross-section.

Eigenfunction expansion method of upper triangular operator matrixand application to two-dimensional elasticity problems based onstress formulation

This paper studies the eigenfunction expansion method to solve the two dimensional （2D） elasticity problems based on the stress formulation. The fundamental system of partial differential equations of the 2D problems is rewritten as an upper tri angular differential system based on the known results, and then the associated upper triangular operator matrix matrix is obtained. By further research, the two simpler com plete orthogonal systems of eigenfunctions in some space are obtained, which belong to the two block operators arising in the operator matrix. Then, a more simple and conve nient general solution to the 2D problem is given by the eigenfunction expansion method. Furthermore, the boundary conditions for the 2D problem, which can be solved by this method, are indicated. Finally, the validity of the obtained results is verified by a specific example.

A Local Discontinuous Galerkin Method for Two-Dimensional Time Fractional Diffusion Equations

For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.

COMPACT FINITE DIFFERENCE-FOURIER SPECTRAL METHOD FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Adaptive split-step Fourier method for simulating ultrashort laser pulse propagation in photonic crystal fibres

In this paper, the generalized nonlinear Schrodinger equation （GNLSE） is solved by an adaptive split-step Fourier method （ASSFM）. It is found that ASSFM must be used to solve GNLSE to ensure precision when the soliton selffrequency shift is remarkable and the photonic crystal fibre （PCF） parameters vary with the frequency considerably. The precision of numerical simulation by using ASSFM is higher than that by using split-step Fourier method in the process of laser pulse propagation in PCFs due to the fact that the variation of fibre parameters with the peak frequency in the pulse spectrum can be taken into account fully.

On the Fourier approximation method for steady water waves

A computational method for steady water waves is presented on the basis of potential theory in the physical plane with spatial variables as independent quantities. The finite Fourier series are applied to approximating the free surface and potential function. A set of nonlinear algebraic equations for the Fourier coefficients are derived from the free surface kinetic and dynamic boundary conditions. These algebraic equations are numerically solved through Newton＇s iterative method, and the iterative stability is further improved by a relaxation technology. The integral properties of steady water waves are numerically analyzed, showing that （1） the set-up and the set-down are both non-monotonic quantities with the wave steepness, and （2） the Fourier spectrum of the free surface is broader than that of the potential function. The latter further leads us to explore a modification for the present method by approximating the free surface and potential function through different Fourier series, with the truncation of the former higher than that of the latter. Numerical tests show that this modification is effective, and can notably reduce the errors of the free surface boundary conditions.

Rational and Periodic Wave Solutions of Two-Dimensional Boussinesq Equation

Two new exact, rational and periodic wave solutions are derived for the two-dimensional Boussinesq equation. For the first solution it is obtained by performing an appropriate limiting procedure on the soliton solutions obtained by Hirota bilinear method. The second one in terms of Riemann theta function is explicitly presented by virtue of Hirota bilinear method and its asymptotic property is also analyzed in detail. Moreover, it is of interest to note that classical soliton solutions can be reduced from the periodic wave solutions.

Ultralow-weight loading Ni catalyst supported on two-dimensional vermiculite for carbon monoxide methanation

Nickel-based catalysts represent the most commonly used systems for CO methanation.We have successfully prepared a Ni catalyst system supported on two-dimensional plasma-treated vermiculite（2D-PVMT）with a very low Ni loading（0.5 wt%）.The catalyst precursor was subjected to heat treatment via either conventional heat treatment（CHT）or the plasma irradiation method（PIM）.The as-obtained CHT-Ni/PVMT and PIM-Ni/PVMT catalysts were characterized with scanning electron microscopy（SEM）,energy dispersive X-ray（EDX）,X-ray diffraction（XRD）,X-ray photoelectron spectroscopy（XPS）,inductively coupled plasma-atomic emission spectroscopy（ICP-AES）and high-angle annular dark field scanning transmission electron microscopy（HAADF-STEM）.Additionally,CHT-NiO/PVMT and PIM-NiO/PVMT catalysts were characterized with hydrogen temperature programmed reduction（H2-TPR）.Compared with CHT-Ni/PVMT,PIM-Ni/PVMT exhibited superior catalytic performance.The plasma treated catalyst PIM-Ni/PVMT achieved a CO conversion of93.5%and a turnover frequency（TOF）of 0.8537 s^-1,at a temperature of 450℃,a gas hourly space velocity of 6000 ml·g^-1·h^-1,a synthesis gas flow rate of 65 ml·min^-1,and a pressure of 1.5 MPa.Plasma irradiation may provide a successful strategy for the preparation of catalysts with very low metal loadings which exhibit excellent properties.

Inversion of two-dimensional tidal open boundary conditions of M_2 constituent in the Bohai and Yellow Seas

Two-dimensional tidal open boundary conditions of the M2 constituent in the Bohai and Yellow Seas(BYS) have been estimated by assimilating T/P altimeter data.During inversion,independent point(IP) strategy was used,in which several IPs on the open boundary is assumed,values at these IPs can be optimized with an adjoint method,and those at other grid points are determined by linearly interpolating the values at IPs.The reasonability and feasibility of the model are tested by ideal twin experiments.In the practical experiment(PE) after assimilation,the cost function may reach 1% or less of its initial value.Mean absolute errors in amplitude and phase can be less than 5 cm and 5°,respectively,and the obtained co-chart can show the character of the M2 constituent in the BYS.The results of the PE indicate that using only two IPs on the open boundary can yield better simulated results.

Numerical Solution of the Coupled Viscous Burgers’ Equation Using Differential Quadrature Method Based on Fourier Expansion Basis

The differential quadrature method based on Fourier expansion basis is applied in this work to solve coupled viscous Burgers’ equation with appropriate initial and boundary conditions. In the first step for the given problem we have discretized the interval and replaced the differential equation by the Differential quadrature method based on Fourier expansion basis to obtain a system of ordinary differential equation (ODE) then we implement the numerical scheme by computer programing and perform numerical solution. Finally the validation of the present scheme is demonstrated by numerical example and compared with some existing numerical methods in literature. The method is analyzed for stability and convergence. It is found that the proposed numerical scheme produces a good result as compared to other researcher’s result and even generates a value at the nodes or mesh points that the results have not seen yet.

Interfacial fracture analysis for a two-dimensional decagonal quasi-crystal coating layer structure

The interface crack problems in the two-dimensional(2D)decagonal quasicrystal(QC)coating are theoretically and numerically investigated with a displacement discontinuity method.The 2D general solution is obtained based on the potential theory.An analogy method is proposed based on the relationship between the general solutions for 2D decagonal and one-dimensional(1D)hexagonal QCs.According to the analogy method,the fundamental solutions of concentrated point phonon displacement discontinuities are obtained on the interface.By using the superposition principle,the hypersingular boundary integral-differential equations in terms of displacement discontinuities are determined for a line interface crack.Further,Green’s functions are found for uniform displacement discontinuities on a line element.The oscillatory singularity near a crack tip is eliminated by adopting the Gaussian distribution to approximate the delta function.The stress intensity factors(SIFs)with ordinary singularity and the energy release rate(ERR)are derived.Finally,a boundary element method is put forward to investigate the effects of different factors on the fracture.