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.

A Low Mach Number IMEX Flux Splitting for the Level Set Ghost Fluid Method

Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.

Nonlinear Algebraic Equations Solved by an Optimal Splitting-Linearizing Iterative Method

How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).

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.

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.

Construction of Conservative Numerical Fluxes for the Entropy Split Method

The entropy split method is based on the physical entropies of the thermally perfect gas Euler equations.The Euler flux derivatives are approximated as a sum of a conservative portion and a non-conservative portion in conjunction with summation-by-parts(SBP)difference boundary closure of(Gerritsen and Olsson in J Comput Phys 129:245-262,1996;Olsson and Oliger in RIACS Tech Rep 94.01,1994;Yee et al.in J Comp Phys 162:33-81,2000).Sj?green and Yee(J Sci Comput)recently proved that the entropy split method is entropy conservative and stable.Stand-ard high-order spatial central differencing as well as high order central spatial dispersion relation preserving(DRP)spatial differencing is part of the entropy stable split methodol-ogy framework.The current work is our first attempt to derive a high order conservative numerical flux for the non-conservative portion of the entropy splitting of the Euler flux derivatives.Due to the construction,this conservative numerical flux requires higher oper-ations count and is less stable than the original semi-conservative split method.However,the Tadmor entropy conservative(EC)method(Tadmor in Acta Numerica 12:451-512,2003)of the same order requires more operations count than the new construction.Since the entropy split method is a semi-conservative skew-symmetric splitting of the Euler flux derivative,a modified nonlinear filter approach of(Yee et al.in J Comput Phys 150:199-238,1999,J Comp Phys 162:3381,2000;Yee and Sj?green in J Comput Phys 225:910934,2007,High Order Filter Methods for Wide Range of Compressible flow Speeds.Proceedings of the ICOSAHOM09,June 22-26,Trondheim,Norway,2009)is proposed in conjunction with the entropy split method as the base method for problems containing shock waves.Long-time integration of 2D and 3D test cases is included to show the com-parison of these new approaches.

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.

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.

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.

Systematic Analysis Method of Shear-Wave Splitting:SAM Software System

In order to make a more effective use of the data from regional digital seismograph networks and to promote the study on shear wave splitting and its application to earthquake stress-forecasting, SAM software system, i.e., the software on systematic analysis method of shear wave splitting has been developed. This paper introduces the design aims, system structure, function and characteristics about the SAM software system and shows some graphical interfaces of data input and result output. Lastly, it discusses preliminarily the study of shear wave splitting and its application to earthquake forecasting.

Tensile strength of sea ice using splitting tests based on the digital image correlation method

The splitting test is a competitive alternative method to study the tensile strength of sea ice owing to its suitability for sampling.However,the approach was questioned to the neglect of local plastic deformation during the tests.In this study,splitting tests were performed on sea ice,with 32 samples subjected to the regular procedure and 8 samples subjected to the digital image correlation method.The salinity,density,and temperature were measured to determine the total porosity.With the advantage of the digital image correlation method,the full-field deformation of the ice samples could be determined.In the loading direction,the samples mainly deformed at the ice-platen contact area.In the direction vertical to the loading,deformation appears along the central line where the splitting crack occurs.Based on the distribution of the sample deformation,a modified solution was derived to calculate the tensile strength with the maximum load.Based on the modified solution,the tensile strength was further calculated together with the splitting test results.The results show that the tensile strength has a negative correlation with the total porosity,which agrees with previous studies based on uniaxial tension tests.

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.

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.

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.

Splitting Rolling Simulated by Reproducing Kernel Particle Method

During splitting rolling simulation, re-meshing is necessary to prevent the effect of severe mesh distortion when the conventional finite element method is used. However, extreme deformation cannot be solved by the finite element method in splitting rolling. The reproducing kernel particle method can solve this problem because the continuum body is discretized by a set of nodes, and a finite element mesh is unnecessary, and there is no explicit limitation of mesh when the metal is split. To ensure stability in the large deformation elastoplastic analysis, the Lagrange material shape function was introduced. The transformation method was utilized to impose the essential boundary conditions. The splitting rolling method was simulated and the simulation results were in accordance with the experimental ones in the 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.

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.

Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.

Operator Splitting Method for Coupled Problems:Transport and Maxwell Equations

In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.