Explicit multi-symplectic method for the Zakharov-Kuznetsov equation

We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.

Multi-symplectic method for generalized Boussinesq equation

The generalized Boussinesq equation that represents a group of important nonlinear equations possesses many interesting properties. Multi-symplectic formulations of the generalized Boussinesq equation in the Hamilton space are introduced in this paper. And then an implicit multi-symplectic scheme equivalent to the multi-symplectic Box scheme is constructed to solve the partial differential equations （PDEs） derived from the generalized Boussinesq equation. Finally, the numerical experiments on the soliton solutions of the generalized Boussinesq equation are reported. The results show that the multi-symplectic method is an efficient algorithm with excellent long-time numerical behaviors for nonlinear partial differential equations.

Multi-symplectic methods for membrane free vibration equation

In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws-a multi-symplectic conservation law （CLS）, a local energy conservation law （ECL） as well as a local momentum conservation law （MCL） --is constructed to discrete the PDEs that are derived from the membrane free vibration equation. The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior,

Multi-symplectic method for generalized fifth-order KdV equation

This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete multi-symplectic conservation law to solve the partial differential equations which are derived from the generalized fifth-order KdV equation numerically. The results of the numerical experiments show that this multi-symplectic algorithm is good in accuracy and its long-time numerical behaviour is also perfect.

Second order conformal multi-symplectic method for the damped Korteweg–de Vries equation

A conformal multi-symplectic method has been proposed for the damped Korteweg–de Vries(DKdV) equation, which is based on the conformal multi-symplectic structure. By using the Strang-splitting method and the Preissmann box scheme,we obtain a conformal multi-symplectic scheme for multi-symplectic partial differential equations(PDEs) with added dissipation. Applying it to the DKdV equation, we construct a conformal multi-symplectic algorithm for it, which is of second order accuracy in time. Numerical experiments demonstrate that the proposed method not only preserves the dissipation rate of mass exactly with periodic boundary conditions, but also has excellent long-time numerical behavior.

Multi-symplectic method for the coupled Schrdinger–KdV equations

In this paper, we present a multi-symplectic Hamiltonian formulation of the coupled Schrtidinger-KdV equations （CS＇KE） with periodic boundary conditions. Then we develop a novel multi-symplectic Fourier pseudospectral （MSFP） scheme for the CSKE. In numerical experiments, we compare the MSFP method with the Crank-Nicholson （CN） method. Our results show high accuracy, effectiveness, and good ability of conserving the invariants of the MSFP method.

A Review on Stochastic Multi-symplectic Methods for Stochastic Maxwell Equations

Stochastic multi-symplectic methods are a class of numerical methods preserving the discrete stochastic multi-symplectic conservation law. These methods have the remarkable superiority to conventional numerical methods when applied to stochastic Hamiltonian partial differential equations (PDEs), such as long-time behavior, geometric structure preserving, and physical properties preserving. Stochastic Maxwell equations driven by either additive noise or multiplicative noise are a system of stochastic Hamiltonian PDEs intrinsically, which play an important role in fields such as stochastic electromagnetism and statistical radiophysics. Thereby, the construction and the analysis of various numerical methods for stochastic Maxwell equations which inherit the stochastic multi-symplecticity, the evolution laws of energy and divergence of the original system are an important and promising subject. The first stochastic multi-symplectic method is designed and analyzed to stochastic Maxwell equations by Hong et al.(A stochastic multi-symplectic scheme for stochastic Maxwell equations with additive noise. J. Comput. Phys. 268:255-268, 2014). Subsequently, there have been developed various stochastic multi-symplectic methods to solve stochastic Maxwell equations. In this paper, we make a review on these stochastic multi-symplectic methods for solving stochastic Maxwell equations driven by a stochastic process. Meanwhile, the theoretical results of well-posedness and conservation laws of the stochastic Maxwell equations are included.

Explicit Multi-symplectic Method for a High Order Wave Equation of KdV Type

In this paper, we consider multi-symplectic Fourier pseudospectral method for a high order integrable equation of KdV type, which describes many important physical phenomena. The multi-symplectic structure are constructed for the equation, and the conservation laws of the continuous equation are presented. The multisymplectic discretization of each formulation is exemplified by the multi-symplectic Fourier pseudospectral scheme. The numerical experiments are given, and the results verify the efficiency of the Fourier pseudospectral method.

Multi-symplectic method for the generalized(2+1)-dimensionalKdV-mKdV equation

In the present paper, a general solution involv- ing three arbitrary functions for the generalized （2＋1）- dimensional KdV-mKdV equation, which is derived from the generalized （1＋1）-dimensional KdV-mKdV equa- tion, is first introduced by means of the Wiess, Tabor, Carnevale （WTC） truncation method. And then multi- symplectic formulations with several conservation laws taken into account are presented for the generalized （2＋1）- dimensional KdV-mKdV equation based on the multi- symplectic theory of Bridges. Subsequently, in order to simulate the periodic wave solutions in terms of rational functions of the Jacobi elliptic functions derived from thegeneral solution, a semi-implicit multi-symplectic scheme is constructed that is equivalent 1：o the Preissmann scheme. From the results of the numerical experiments, we can con- clude that the multi-symplectic schemes can accurately sim- ulate the periodic wave solutions of the generalized （2＋1）- dimensional KdV-mKdV equation while preserve approxi- mately the conservation laws.

Multi-symplectic Runge-Kutta methods for Landau-Ginzburg-Higgs equation

Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.

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.

Multi-symplectic method to analyze the mixed state of Ⅱ-superconductors

The mixed state of two-band II-superconductor is analyzed by the multi-symplectic method. As to the Ginzburg-Landau equation depending on time that describes the mixed state of two-band II-superconductor, the multi-symplectic formulations with several conservation laws: a multi-symplectic conservation law, an energy con- servation law, as well as a momentum conservation law, are presented firstly; then an eighteen points scheme is constructed to simulate the multi-symplectic partial differential equations (PDEs) that are derived from the Ginzburg-Landau equation; finally, based on the simulation results, the volt-ampere characteristic curves are obtained, as well as the relationships between the temperature and resistivity of a suppositional two-band II-superconductor model under different magnetic intensi- ties. From the results of the numerical experiments, it is concluded that the notable property of the mixed state of the two-band II-superconductor is that: The trans- formation temperature decreases and the resistivity ρ increases rapidly with the increase of the magnetic intensity B. In addition, the simulation results show that the multi-symplectic method has two remarkable advantages: high accuracy and excellent long-time numerical behavior.

A Review on Sources,Extractions and Analysis Methods of a Sustainable Biomaterial:Tannins

Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

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.

Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

Sparse Modal Decomposition Method Addressing Underdetermined Vortex-Induced Vibration Reconstruction Problem for Marine Risers

When investigating the vortex-induced vibration(VIV)of marine risers,extrapolating the dynamic response on the entire length based on limited sensor measurements is a crucial step in both laboratory experiments and fatigue monitoring of real risers.The problem is conventionally solved using the modal decomposition method,based on the principle that the response can be approximated by a weighted sum of limited vibration modes.However,the method is not valid when the problem is underdetermined,i.e.,the number of unknown mode weights is more than the number of known measurements.This study proposed a sparse modal decomposition method based on the compressed sensing theory and the Compressive Sampling Matching Pursuit(Co Sa MP)algorithm,exploiting the sparsity of VIV in the modal space.In the validation study based on high-order VIV experiment data,the proposed method successfully reconstructed the response using only seven acceleration measurements when the conventional methods failed.A primary advantage of the proposed method is that it offers a completely data-driven approach for the underdetermined VIV reconstruction problem,which is more favorable than existing model-dependent solutions for many practical applications such as riser structural health monitoring.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.