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.

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.

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.

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.

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.

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).

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 symplecticity exactly.Numerical experiments for the plane wave solution and singular solution of the 2D-NLSE show the accuracy and effectiveness of the proposed method.

Engineering of Self-Supported Electrocatalysts on a Three-Dimensional Nickel Foam Platform for Efficient Water Electrolysis

Economical water electrolysis requires highly active non-noble electrocatalysts to overcome the sluggish kinetics of the two half-cell reactions,oxygen evolution reaction,and hydrogen evolution reaction.Although intensive efforts have been committed to achieve a hydrogen economy,the expensive noble metal-based catalysts remain under consideration.Therefore,the engineering of self-supported electrocatalysts prepared using a direct growth strategy on three-dimensional(3D)nickel foam(NF)as a conductive substrate has garnered significant interest.This is due to the large active surface area and 3D porous network offered by these electrocatalysts,which can enhance the synergistic eff ect between the catalyst and the substrate,as well as improve electrocatalytic performance.Hydrothermal-assisted growth,microwave heating,electrodeposition,and other physical methods(i.e.,chemical vapor deposition and plasma treatment)have been applied to NF to fabricate competitive electrocatalysts with low overpotential and high stability.In this review,recent advancements in the development of self-supported electrocatalysts on 3D NF are described.Finally,we provide future perspectives of self-supported electrode platforms in electrochemical water splitting.

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.

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.

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.

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.

A Successive Shift Box-Counting Method for Calculating Fractal Dimensions and Its Application in Identification of Faults

Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This paper proposes a successive shift box-counting method,in which the studied object is divided into small sub-objects that are composed of a series of gridsaccording to its characteristic scaling. The terrain fractal dimensions in the grids are calculatedwith the successive shift box-counting method and the scattered points with values of fractaldimensions are obtained. The present research shows that the planar variation of fractal dimensionsis well consistent with fault traces and geological boundaries.

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.

Three-dimensional Extension of the Unit-Feature Spatial Classification Method for Cloud Type

We describe how the Unit-Feature Spatial Classification Method(UFSCM) can be used operationally to classify cloud types in satellite imagery efficiently and conveniently.By using a combination of Interactive Data Language(IDL) and Visual C++(VC) code in combination to extend the technique in three dimensions(3-D),this paper provides an efficient method to implement interactive computer visualization of the 3-D discrimination matrix modification,so as to deal with the bi-spectral limitations of traditional two dimensional(2-D) UFSCM.The case study of cloud-type classification based on FY-2C satellite data (0600 UTC 18 and 0000 UTC 10 September 2007) is conducted by comparison with ground station data, and indicates that 3-D UFSCM makes more use of the pattern recognition information in multi-spectral imagery,resulting in more reasonable results and an improvement over the 2-D 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.

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.