THE NAVIER-STOKES EQUATIONS IN STREAM LAYER AND ON STREAM SURFACE AND A DIMENSION SPLIT METHODS

In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on the stream surface and 2-D Navier-Stokes equations on a two dimensional manifold. After introducing stream function on the stream surface,a nonlinear initial-boundary value problem satisfies by stream function is obtained,existence and uniqueness of its solution are proven.Based this theory we proposal a new method called"dimension split method"to solve 3D NSE.

The Boundary Layer Equations and a Dimensional Split Method for Navier-Stokes Equations in Exterior Domain of a Spheroid and Ellipsoid

In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of the obstacle. BLE are nonlinear partial differential equations on unknown normal viscous stress tensor and pressure on the obstacle and the existence of solution of BLE is proved. In addition a dimensional split method for dimensional three Navier-Stokes equations is established by applying several 2D-3C partial differential equations on two dimensional manifolds to approach 3D Navier-Stokes equations. The examples for the exterior flow around spheroid and ellipsoid are presents here.

A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.

A SUPERLINEARLY CONVERGENT SPLITTING FEASIBLE SEQUENTIAL QUADRATIC OPTIMIZATION METHOD FOR TWO-BLOCK LARGE-SCALE SMOOTH OPTIMIZATION

This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

Another SSOR Iteration Method

Kellogg gave a version of the Peaceman-Radford method. In this paper, we introduce a SSOR iteration method which uses Kellogg’s method. The new algorithm has some advantages over the traditional SSOR algorithm. A Cyclic Reduction algorithm is introduced via a decoupling in Kellogg’s method.

The dimension split element-free Galerkin method for three-dimensional potential problems

This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.

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.

An Improved Splitting Method

In this paper, an improved splitting method, based on the completely square-conservative explicit difference schemes, is established. Not only can the time-direction precision of this method be higher than that of the traditional splitting methods but also can the physical feature of mutual dependence of the fast and the slow stages that are calculated separately and splittingly be kept as well. Moreover, the method owns an universality, it can be generalized to other square-conservative difference schemes, such as the implicit and complete ones and the explicit and instantaneous ones. Good time benefits can be acquired when it is applied in the numerical simulations of the monthly mean currents of the South China Sea.

Some Splitting Methods for Equations of Geophysical Fluid Dynamics

In this paper, equations of atmospheric and oceanic dynamics are reduced to a kind of evolutionary equation in operator form, based on which a conclusion that the separability of motion stages is relative is made and an issue that the tractional splitting methods established on the physical separability of the fast stage and the slow stage neglect the interaction between the two stages to some extent is shown. Also, three splitting patterns are summed up from the splitting methods in common use so that a comparison between them is carried out. The comparison shows that only the improved splitting pattern (ISP) can be in second order and keep the interaction well. Finally, the applications of some splitting methods on numerical simulations of typhoon tracks made clear that ISP owns the best effect and can save more than 80% CPU time.

APPROXIMATION OF THE CONCENTRATION BY A VISCOSITY SPLITTING METHOD FOR TWO-PHASE MISCIBLE DISPLACEMENT PROBLEM

The miscible displacement of one incompressible fluid by another in a porous medium is considered in this paper. The concentration is split in a first-order hyberbolic equation and a homogeneous parabolic equation within each lime step. The pressure and Us velocity field is computed by a mixed finite element method. Optimal order estimates are derived for the no diffusion case and the diffusion case.

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.

Preconditioned Iterative Method for Regular Splitting

Several preconditioners are proposed for improving the convergence rate of the iterative method derived from splitting. In this paper, the comparison theorem of preconditioned iterative method for regular splitting is proved. And the convergence and comparison theorem for any preconditioner are indicated. This comparison theorem indicates the possibility of finding new preconditioner and splitting. The purpose of this paper is to show that the preconditioned iterative method yields a new splitting satisfying the regular or weak regular splitting. And new combination preconditioners are proposed. In order to denote the validity of the comparison theorem, some numerical examples are shown.

Application of Signal Processing Methods to Shear-Wave Splitting Study

This paper briefly discusses the new methods that the authors have put forward to distinguish splitting shear-waves.By combining these new methods with other methods,the authors have processed the recorded data of an earthquake.The study results are consistent with each other.

Multidimensional Numerical Simulation of Glow Discharge by Using the N-BEE-Time Splitting Method

In this work, a new numerical technique is proposed for the resolution of a fluid model based on three Boltzmann moments. The main purpose of this technique is to calculate electric and physical properties in the non-equilibrium electric discharge at low pressure. The transport and Poisson＇s equations form a self-consistent model. This equation system is written in cylindrical coordinates following the geometric shape of a plasma reactor. Our transport equation system is discretized using the finite volume approach and resolved by the N-BEE explicit scheme coupled to the time splitting method. This programming structure reduces computation time considerably. The 2D code is carried out and tested by comparing our results with those found in literature.

Multi-symplectic wavelet splitting method for the strongly coupled Schrodinger system

We propose a multi-symplectic wavelet splitting equations. Based on its mu]ti-symplectic formulation, method to solve the strongly coupled nonlinear SchrSdinger the strongly coupled nonlinear SchrSdinger equations can be split into one linear multi-symplectic subsystem and one nonlinear infinite-dimensional Hamiltonian subsystem. For the linear subsystem, the multi-symplectic wavelet collocation method and the symplectic Euler method are employed in spatial and temporal discretization, respectively. For the nonlinear subsystem, the mid-point symplectic scheme is used. Numerical simulations show the effectiveness of the proposed method during long-time numerical calculation.

THE SOLUTION FOR TRANSIENT TWO-PHASE FLOW BY SPLIT FLUX VECTOR METHOD

In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difference operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.

SPLITTING MODULUS FINITE ELEMENT METHOD FOR ORTHOGONAL ANISOTROPIC PLATE BENGING

Splitting modulus variational principle in linear theory of solid mechanics was introduced, the principle for thin plate was derived, and splitting modulus finite element method of thin plate was established too. The distinctive feature of the splitting model is that its functional contains one or more arbitrary additional parameters, called splitting factors, so stiffness of the model can be adjusted by properly selecting the splitting factors. Examples show that splitting modulus method has high precision and the ability to conquer some ill-conditioned problems in usual finite elements. The cause why the new method could transform the ill-conditioned problems into well-conditioned problem, is analyzed finally.

Tip-splitting instability in directional solidification based on bias field method

Tip splitting instability of cellular interface morphology in directional solidification is analyzed based on the bias field method proposed recently by Glicksman. The physical mechanism of tip instability is explained by analyzing the interface potential, the tangential energy flux, and the normal energy flux. A rigorous criterion for tip-splitting instability is established analytically, i.e., the ratio of the cellular tip radius to the cellular width α 〉3/2/π≈ 0.3899, which is in good agreement with simulation results. This study also reveals that the cellular tip splitting instability is attributable to weak Gibbs–Thomson energy acting on the interface.

OPERATOR-SPLITTING GALERKIN METHOD FOR ONE KIND OF OIN REACTION MODEL FOR THE POLLUTION IN GROUNDWATER

The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solution.