Average vector field methods for the coupled Schrdinger KdV equations

The energy preserving average vector field(AVF) method is applied to the coupled Schro¨dinger–KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.     

Average vector field methods for the coupled Schr/idinger-KdV equations

The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.     

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation is proposed in this paper to improve the efficiency and accuracy of electromagnetic calculation. An image boundary extraction algorithm is developed to map the image on the geometric domain. Identification algorithm for the location of nodes in polygon area is proposed to determine the state of the node. To promote the average quality of the mesh and the efficiency of mesh generation, a novel force-based mesh smoothing algorithm is proposed. One test case and a typical electromagnetic calculation are used to testify the effectiveness and efficiency of the proposed algorithm. The results demonstrate that the proposed algorithm can produce a high-quality mesh with less iteration.     

A Simple Framework of Conservative Algorithms for the Coupled Nonlinear Schrdinger Equations with Multiply Components

Considering the coupled nonlinear Schr¨odinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite difference method, Fourier pseudospectral method and wavelet collocation method for spatial discretizations, a series of high accurate conservative algorithms are presented. The proposed algorithms can preserve the corresponding discrete charge and energy conservation laws exactly, which would guarantee their numerical stabilities during long time computations.Furthermore, several analogous multi-symplectic algorithms are constructed as comparison. Numerical experiments for the unstable plane waves will show the advantages of the proposed algorithms over long time and verify the theoretical analysis.     

A high order energy preserving scheme for the strongly coupled nonlinear Schr¨odinger system

A high order energy preserving scheme for a strongly coupled nonlinear Schr¨odinger system is proposed by using the average vector field method. The high order energy preserving scheme is applied to simulate the soliton evolution of the strongly coupled Schr¨odinger system. Numerical results show that the high order energy preserving scheme can well simulate the soliton evolution, moreover, it preserves the discrete energy of the strongly coupled nonlinear Schr¨odinger system exactly.     

Time reversal focusing with variable depth and range based on mode extraction

A method for time reversal focusing with variable depth and range based on mode extraction was proposed.First,the normal modes of acoustic propagation in the shallow water are extracted by modal decomposition from the probe signals received by a source receiver array.Furthermore,a diagonal matrix and a vector determined separately by the depth and the range of the probe source are extracted from the received acoustic field data.And time reversal focusing at different depths and ranges can be achieved by modulating the depth-dependent diagonal matrix and the range-dependent vector properly.Then the diagonal matrix and the vector are modulated separately according to the depth and the range of the expected focal location to construct a new acoustic field vector.When this new acoustic field vector is retransmitted by the source receiver array in time reversal order(or phase conjugation in frequency domain),focusing of the resulting acoustic field at the expected location rather than the origin of the probe source can be obtained.Numerical simulations in typical shallow water environment demonstrate the effectiveness of the proposed method.     

Chaotic time series prediction using fuzzy sigmoid kernel-based support vector machines

Support vector machines (SVM) have been widely used in chaotic time series predictions in recent years. In order to enhance the prediction efficiency of this method and implement it in hardware, the sigmoid kernel in SVM is drawn in a more natural way by using the fuzzy logic method proposed in this paper. This method provides easy hardware implementation and straightforward interpretability. Experiments on two typical chaotic time series predictions have been carried out and the obtained results show that the average CPU time can be reduced significantly at the cost of a small decrease in prediction accuracy, which is favourable for the hardware implementation for chaotic time series prediction.     

Echo-to-reverberation enhancement based on decomposition of time reversal operator with forward and backward reverberation nulling constrains

Combined the decomposition of time reversal operator and the time reversal reverberation nulling, a new time reversal processing approach for echo-to-reverberation ratio enhancement is proposed. In this method, a 2-dimensional signal subspace for the range of the target and two bottom focusing weight vectors for the ranges near the target are obtained by the decomposition of time reversal operator. From the signal subspace and focusing weight vectors, a constrained optimal excitation weight vector of source receiver array can be deduced to null the acoustic energy on the corresponding bottom and maximize the energy at the tar- get. This method remedies the shortages of conventional time reversal processing, time reversal reverberation nulling and time reversal selective focusing method. It focuses sound energy at the target and nulls the energy at the bottom near the target range simultaneously, therefore enhancing the echo-to-reverberation ratio without probe source and prior-knowledge of the relative scattering intensity of target and bottom. Numerical simulations in typical shallow water environments showed the effectiveness of the proposed method and its improved performance for echo-reverberation enhancement than conventional time reversal processing.     

A Simple Framework of Conservative Algorithms for the Coupled Nonlinear Schrodinger Equations with Multiply Components

Considering the coupled nonlinear Schrodinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite difference method, Fourier pseudospectral method and wavelet collocation method for spatial discretizations, a series of high accurate conservative algorithms are presented. The proposed algorithms can preserve the corresponding discrete charge and energy conservation laws exactly, which would guarantee their numerical stabilities during long time computations. Furthermore, several analogous multi-symplectic algorithms are constructed as comparison. Numerical experiments for the unstable plane waves will show the advantages of the proposed algorithms over long time and verify the theoretical analysis.     

An improved calculation method for fiber Raman amplifier equations with multi-wavelength pumping

A novel numerical method for fiber Raman amplifier (FRA) from standard propagation equations is presented and derived based on the one-step method for ordinary differential equation (ODE). The proposed algorithm is effective in solving FRA equations including all the interactions among pumps, signals, and noises. Applications of the numerical analysis to practical FRA-based systems show a great reduction in computation time in comparison with the average power method and the fourth-order Runge-Kutta (RK) method, under the same condition. Also the proposed method can decrease the computing time over three orders of magnitude with excellent accuracy promises in comparison with the direct integration method.     

Radiation diffusion for multi-fluid Eulerian hydrodynamics with adaptive mesh refinement

Block-structured meshes provide the ability to concentrate grid points and computational effort in interesting regions of a flow field, without sacrificing the efficiency and low memory requirements of a regular grid. We describe an algorithm for simulating radiation diffusion on such a mesh, coupled to multi-fluid gasdynamics. Conservation laws are enforced by using locally conservative difference schemes along with explicit synchronization operations between different levels of refinement. In unsteady calculations each refinement level is advanced at its own optimal timestep. Particular attention is given to the appropriate coupling between the fluid energy and the radiation field, the behavior of the discretization at sharp interfaces, and the form of synchronization between levels required for energy conservation in the diffusion process. Two- and three-dimensional examples are presented, including parallel calculations performed on an IBM SP-2.     

Minimization of the Truncation Error by Grid Adaptation

A new grid adaptation strategy, which minimizes the truncation error of a pth-order finite difference approximation, is proposed. The main idea of the method is based on the observation that the global truncation error associated with discretization on nonuniform meshes can be minimized if the interior grid points are redistributed in an optimal sequence. The method does not explicitly require the truncation error estimate, and at the same time, it allows one to increase the design order of approximation globally by one, so that the same finite difference operator reveals superconvergence properties on the optimal grid. Another very important characteristic of the method is that if the differential operator and the metric coefficients are evaluated identically by some hybrid approximation, then the single optimal grid generator can be employed in the entire computational domain independently of points where the hybrid discretization switches from one approximation to another. Generalization of the present method to multiple dimensions is presented. Numerical calculations of several one-dimensional and one two-dimensional test examples demonstrate the performance of the method and corroborate the theoretical results.     

Immersed-boundary methods for general finite-difference and finite-volume Navier–Stokes solvers

We present an immersed-boundary algorithm for incompressible flows with complex boundaries, suitable for Cartesian or curvilinear grid system. The key stages of any immersed-boundary technique are the interpolation of a velocity field given on a mesh onto a general boundary (a line in 2D, a surface in 3D), and the spreading of a force field from the immersed boundary to the neighboring mesh points, to enforce the desired boundary conditions on the immersed-boundary points. We propose a technique that uses the Reproducing Kernel Particle Method [W.K. Liu, S. Jun, Y.F. Zhang, Reproducing kernel particle methods, Int. J. Numer. Methods Fluids 20(8) (1995) 1081–1106] for the interpolation and spreading. Unlike other methods presented in the literature, the one proposed here has the property that the integrals of the force field and of its moment on the grid are conserved, independent of the grid topology (uniform or non-uniform, Cartesian or curvilinear). The technique is easy to implement, and is able to maintain the order of the original underlying spatial discretization. Applications to two- and three-dimensional flows in Cartesian and non-Cartesian grid system, with uniform and non-uniform meshes are presented.     

A low numerical dissipation immersed interface method for the compressible Navier–Stokes equations

A numerical method to solve the compressible Navier–Stokes equations around objects of arbitrary shape using Cartesian grids is described. The approach considered here uses an embedded geometry representation of the objects and approximate the governing equations with a low numerical dissipation centered finite-difference discretization. The method is suitable for compressible flows without shocks and can be classified as an immersed interface method. The objects are sharply captured by the Cartesian mesh by appropriately adapting the discretization stencils around the irregular grid nodes, located around the boundary. In contrast with available methods, no jump conditions are used or explicitly derived from the boundary conditions, although a number of elements are adopted from previous immersed interface approaches. A new element in the present approach is the use of the summation-by-parts formalism to develop stable non-stiff first-order derivative approximations at the irregular grid points. Second-order derivative approximations, as those appearing in the transport terms, can be stiff when irregular grid points are located too close to the boundary. This is addressed using a semi-implicit time integration method. Moreover, it is shown that the resulting implicit equations can be solved explicitly in the case of constant transport properties. Convergence studies are performed for a rotating cylinder and vortex shedding behind objects of varying shapes at different Mach and Reynolds numbers.     

Electromagnetics via the Taylor-Galerkin Finite Element Method on Unstructured Grids

Traditional techniques for computing electromagnetic solutions in the time domain rely on finite differences. These so-called FDTD (finite-difference time-domain) methods are usually defined only on regular lattices of points and can be too restrictive for geometrically demanding problems. Great geometric flexibility can be achieved by abandoning the regular latticework of sample points and adopting an unstructured grid. An unstructured grid allows one to place the grid points anywhere one chooses, so that curved boundaries can be fit with ease and local regions in which the field gradients are steep can be selectively resolved with a fine mesh. In this paper we present a technique for solving Maxwell's equations on an unstructured grid based on the Taylor-Galerkin finite-element method. We present several numerical examples which reveal the fundamental accuracy and adaptability of the method. Although our examples are in two dimensions, the techniques and results generalize readily to 3D.     

A grid based particle method for solving partial differential equations on evolving surfaces and modeling high order geometrical motion

We develop numerical methods for solving partial differential equations (PDE) defined on an evolving interface represented by the grid based particle method (GBPM) recently proposed in [S. Leung, H.K. Zhao, A grid based particle method for moving interface problems, J. Comput. Phys. 228 (2009) 7706–7728]. In particular, we develop implicit time discretization methods for the advection–diffusion equation where the time step is restricted solely by the advection part of the equation. We also generalize the GBPM to solve high order geometrical flows including surface diffusion and Willmore-type flows. The resulting algorithm can be easily implemented since the method is based on meshless particles quasi-uniformly sampled on the interface. Furthermore, without any computational mesh or triangulation defined on the interface, we do not require remeshing or reparametrization in the case of highly distorted motion or when there are topological changes. As an interesting application, we study locally inextensible flows governed by energy minimization. We introduce tension force via a Lagrange multiplier determined by the solution to a Helmholtz equation defined on the evolving interface. Extensive numerical examples are also given to demonstrate the efficiency of the proposed approach.     

Grid-based photogrammetry system for large scale sheet metal strain measurement

This paper presents a new approach to fast strain measurement with high accuracy for large scale sheet metal based on the surface circular grid and digital close range photogrammetry. A multi-block measuring method of discretization is implemented to archive large scale measurement. The sheet metal is separated into several blocks for respective calculating and joined together by common reference points. A surface circular searching method is presented for fast and robust 3D grid generation. A flexible bundle adjustment method is proposed for large amount 3D grid nodes reconstruction, which employs the conception of sampling points and is proved to be efficient. Furthermore, a multi-stage grid registration method is introduced to improve the accuracy of strain field by correcting the true deformation gradient tensor. A novel system is developed and performances well in actual large scale sheet metal strain measurement. Two accuracy tests confirm that the system strain measurement error is less than 0.2%.     

Conformal structure-preserving method for damped nonlinear Schrdinger equation

In this paper,we propose a conformal momentum-preserving method to solve a damped nonlinear Schrodinger(DNLS) equation.Based on its damped multi-symplectic formulation,the DNLS system can be split into a Hamiltonian part and a dissipative part.For the Hamiltonian part,the average vector field(AVF) method and implicit midpoint method are employed in spatial and temporal discretizations,respectively.For the dissipative part,we can solve it exactly.The proposed method conserves the conformal momentum conservation law in any local time-space region.With periodic boundary conditions,this method also preserves the total conformal momentum and the dissipation rate of momentum exactly.Numerical experiments are presented to demonstrate the conservative properties of the proposed method.     

NS方程的各向异性笛卡尔网格法研究

将近年发展起来的用于Euler方程求解的具有局部均匀网格总体非结构特性的笛卡尔网格法推广到NS方程的求解。为了与流场的各向异性相适应、减少网格点数量,提出了一种各向异性网格加密法。另外还研究了分级笛卡尔网格对内点格式稳定性的影响和插值固体边界条件的稳定性。数值结果表明各向异性笛卡尔网格法相对于传统的各向同性网格方法能大量节省网格点数量而且与后者具有同样的精度。     

Accuracy and Conservation Properties of a Three-Dimensional Unstructured Staggered Mesh Scheme for Fluid Dynamics

A three-dimensional unstructured mesh discretization of the rotational from of the incompressible Navier–Stokes is presented. The method uses novel and highly efficient algorithms for interpolating the velocity vector and constructing the convention term. The resulting discretization is shown to conserve mass, kinetic energy, and vorticity to machine precision both locally and globally. The spatial accuracy of the method is analyzed and found to be second order on regular meshes and first order on irregular meshes. The numerical efficiency, accuracy, and conservation properties of the method are tested on three-dimensional meshes and found to be in agreement with theory.