Data-driven optimization study of the multi-relaxation-time lattice Boltzmann method for solid-liquid phase change

Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes.The multi-relaxation-time lattice Boltzmann method(MRT-LBM)shows great numerical performance during simulation;however,the value method of the relaxation parameters needs to be specified.Therefore,in this study,a random forest(RF)model is used to discriminate the importance of different relaxation parameters to the convergence,and a support vector machine(SVM)is used to explore the decision boundary of the convergent samples in each dimensional model.The results show that the convergence of the samples is consistent with the sign of the decision number,and two types of the numerical deviations appear,i.e.,the phase mushy zone and the non-physical heat transfer.The relaxation parameters chosen on the decision boundary can further suppress the numerical bias and improve numerical accuracy.

ASYNCHRONOUS RELAXED ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS OF EQUATIONS

In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.

ANALYSIS OF WAVEGUIDE PROBLEMS USING A RELAXED ITERATIVE DOMAIN DECOMPOSITION METHOD COMBINED WITH MULTIFRONTAL ALGORITHM

In this paper, an absorbing Fictitious Boundary Condition (FBC) is presented to generate an iterative Domain Decomposition Method (DDM) for analyzing waveguide problems.The relaxed algorithm is introduced to improve the iterative convergence. And the matrix equations are solved using the multifrontal algorithm. The resulting CPU time is greatly reduced.Finally, a number of numerical examples are given to illustrate its accuracy and efficiency.

ALGEBRAIC MULTILEVEL RELAXED PRECONDITIONING METHODS

Some ways of multilevel relaxed preconditioning matrices for the stiffness matrix in the discretization of selfad joint second order elliptic boundary value problems are proposed. For reason-able assumptions of the relaxed factor ω, smaller relative condition numbers are given. The optimal relaxed factor ω is derived, too.

Three-step relaxed hybrid steepest-descent methods for variational inequalities

The classical variational inequality problem with a Lipschitzian and strongly monotone operator on a nonempty closed convex subset in a real Hilbert space is studied. A new three-step relaxed hybrid steepest-descent method for this class of variational inequalities is introduced. Strong convergence of this method is established under suitable assumptions imposed on the algorithm parameters.

AN EXTRAGRADIENT METHOD FOR RELAXED COCOERCIVE VARIATIONAL INEQUALITY AND EQUILIBRIUM PROBLEMS

The purpose of this paper is to investigate the problem of finding the common element of the set of common fixed points of a countable family of nonexpansive mappings, the set of an equilibrium problem and the set of solutions of the variational inequality prob- lem for a relaxed cocoercive and Lipschitz continuous mapping in Hilbert spaces. Then, we show that the sequence converges strongly to a common element of the above three sets under some parameter controlling conditions, which are connected with Yao, Liou, Yao[17], Takahashi[12] and many others.

A General Projection Method for a System of Relaxed Coercive Variational Inequalities in Hilbert Spaces

In this paper,we consider a new algorithm for a generalized system for relaxed coercive nonlinear inequalities involving three different operators in Hilbert spaces by the convergence of projection methods.Our results include the previous results as special cases extend and improve the main results obtained by many others.

Modulus-Based Multisplitting Iteration Method for a Class of Weakly Nonlinear Complementarity Problem

In this paper,we present a modulus-based multisplitting iteration method based on multisplitting of the system matrix for a class of weakly nonlinear complementarity problem.And we prove the convergence of the method when the system matrix is an H_(+)-matrix.Finally,we give two numerical examples.

On Solving a System of Volterra Integral Equations with Relaxed Monte Carlo Method

A random simulation method was used for treatment of systems of Volterra integral equations of the second kind. Firstly, a linear algebra system was obtained by discretization using quadrature formula. Secondly, this algebra system was solved by using relaxed Monte Carlo method with importance sampling and numerical approximation solutions of the integral equations system were achieved. It is theoretically proved that the validity of relaxed Monte Carlo method is based on importance sampling to solve the integral equations system. Finally, some numerical examples from literatures are given to show the efficiency of the method.

A Class of Smoothing Modulus-Based Iterative Method for Solving Implicit Complementarity Problems

In this paper, a class of smoothing modulus-based iterative method was presented for solving implicit complementarity problems. The main idea was to transform the implicit complementarity problem into an equivalent implicit fixed-point equation, then introduces a smoothing function to obtain its approximation solutions. The convergence analysis of the algorithm was given, and the efficiency of the algorithms was verified by numerical experiments.

GENERALIZED MATRIX MULTISPLITTING RELAXATION METHODS AND THEIR CONVERGENCE

In this paper, we set up a general framework of parallel matrix mullisplitting relaxation methods for solving large scale system of linear equations. We investigate the convergence properties of this framework and give several sufficient conditions ensuring it to converge as well as diverge. At last, we conclude a necessary and sufficient condition for the convergence of this framework when the coefficient matrix is an L-matrix.

A time domain induced polarization relaxation time spectrum inversion method based on a damping factor and residual correction

Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.

Geometric interpretation of several classical iterative methods for linear system of equations and diverse relaxation parameter of the SOR method

Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.

Combined immersed boundary method and multiple-relaxation-time lattice Boltzmann flux solver for numerical simulations of incompressible flows

A method combining the immersed boundary technique and a multi- relaxation-time （MRT） lattice Boltzmann flux solver （LBFS） is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations （LBEs）. The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.

A multiple-relaxation-time lattice Boltzmann method for high-speed compressible flows

This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time （SRT） lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.

Influence of nanosized magnesia on the hydration of borehole-sealing cements prepared using different methods

Gas drainage is an efective technology for gas control in coal mines.A high borehole-sealing quality is the fundamental precondition for efcient gas drainage.The expansibilities of cement pastes used in borehole-sealing processes are critical for the borehole-sealing efect.Nanosized magnesia expansive agents are used to improve the expansibilities of cement pastes and improve the borehole-sealing efect.Nuclear magnetic resonance spectrometry and scanning electron microscopy were adopted to study the efects of nanosized magnesia on the hydration of borehole-sealing cements used with diferent preparation methods.The results showed that an increase in the mass fraction of the nanosized magnesia promoted cement hydration,and the mass fraction was positively correlated with the promotion efect.The use of diferent preparation methods did not change the water-phase distribution in the cement.When using the wet-mixing preparation method,nanosized magnesia promoted the induction,acceleration,and deceleration periods of hydration;when using the dry-mixing preparation method,the nanosized magnesia promoted the induction period of cement hydration,and the promotion efect was less obvious than that seen when using the wet-mixing method.When using the wet-mixing preparation method,the nanosized magnesia was uniformly dispersed,thus enlarging the surface area of the reaction,which provided more nucleation sites for the hydration products of the cement and therefore accelerated the hydration reaction.When using the dry-mixing preparation method,the nanosized magnesia powders were dispersed nonuniformly and aggregated.Under these conditions,only a few nanosized magnesia particles on the surfaces of the aggregated clusters took part in hydration,so only a small number of nucleation sites were provided for the hydration products of cement.This led to inconsistent hydration of cement pastes prepared using the dry-mixing method.The surface porosity of the cement prepared with the wet-mixing preparation method frst decreased and then increased with increases in the mass fraction of the nanosized magnesia.The cement surface exhibited compact hydration products and few pores,and the surface was relatively smooth.In comparison,the surface porosity of the cement prepared using the dry-mixing method fuctuated with increasing mass fraction of the nanosized magnesia,resulting in a rough cement surface and microfractures on some surfaces.The two preparation methods both reduced the surface porosity of the cement.The wet-mixing preparation was more efective and consistent in improving the compactness of the cement than the dry-mixing preparation.These results provide important guidance on the addition of nanosized magnesia in borehole-sealing engineering and the selection of cement preparation methods,and they also lay a solid foundation for realizing safe and efcient gas drainage.

PARALLEL NONLINEAR MULTISPLITTING RELAXATION METHODS

By further generalizing Frommer's results in the sense of nonlinear multisplitting, we build a class of nonlinear multisplitting AOR-type methods, which covers many rather practical nonlinear multisplitting relaxation methods such as multisplitting AOR-Newton method, multisplitting AOR-chord method and multisplitting AOR-Steffensen method, etc.. Furthermore,a general convergence theorem for the nonlinear multisplitting AOR-type methods and the local convergence for the multisplitting AOR-Newton method are discussed in detail.A lot of numerical tests show that our new methods are feasible and satisfactory.

A new system of nonlinear variational inequalities with relaxed cocoercive mappings in reflexive Banach spaces

A new system for relaxed cocoercive non-linear variational inequalities in uniformly smooth Banach spaces is introduced and studied using the convergence of projection methods.Our results generalize and improve the corresponding results of recent works.

Hydrodynamic Lubrication of Elastic Foil Gas Bearing Using Over Relaxation Iteration Method and Non-Dimensional Equation

The purpose is to accurately predict the performance of foil bearing and achieve accurate results in the design of foil bearing structure.A new type of foil bearing with surface microstructure is used as experimental material.First,the lubrication mechanism of elastic foil gas bearing is analyzed.Then,the numerical solution process of the static bearing capacity and friction torque is analyzed,including the discretization of the governing equation of rarefied gas pressure based on the non-dimensional modified Reynolds equation and the over relaxation iteration method,the grid planning within the calculation range,the static solution of boundary parameters and static solution of the numerical process.Finally,the solution program is analyzed.The experimental data in National Aeronautics and Space Administration(NASA)public literature are compared with the simulation results of this exploration,so as to judge the accuracy of the calculation process.The results show that under the same static load,the difference between the minimum film thickness calculated and the test results is not obvious;when the rotor speed of the bearing is 60000 r/min,the influence of the boundary slip effect increases with the increase of the micro groove depth on the flat foil surface;when the eccentricity or the micro groove depth of the bearing increases,the bearing capacity will be strengthened.When the eccentricity is 6µm and 14µm,the viscous friction torque of the new foil bearing increases significantly with the increase of the depth of the foil micro groove,but when the eccentricity is 22µm,the viscous friction torque does not change with the change of the depth of the foil micro groove.It shows that the bearing capacity and performance of foil bearing are improved.