A High-Accuracy Curve Boundary Recognition Method Based on the Lattice Boltzmann Method and Immersed Moving Boundary Method

Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.

A moving Kriging interpolation-based boundary node method for two-dimensional potential problems

In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.

A meshless method based on moving Kriging interpolation for a two-dimensional time-fractional diffusion equation

Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

High-precision solution to the moving load problem using an improved spectral element method

In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.

BOUNDARY ELEMENT METHOD FOR MOVING AND ROLLING CONTACT OF 2D ELASTIC BODIES WITH DEFECTS

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.

Assessing dynamic response of multispan viscoelastic thin beams under a moving mass via generalized moving least square method

Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatial domain using generalized moving least square method （GMLSM） and then, discrete equations of motion based on Lagrange＇s equation are obtained. Maximum deflection and bending moments are considered as the important design parameters. The design parameter spectra in terms of mass weight and velocity of the moving mass are presented for multispan viscoelastic beams as well as various values of relaxation rate and beam span number. A reasonable good agreement is achieved between the results of the proposed solution and those obtained by other researchers. The results indicate that, although the load inertia effects in beams with higher span number would be intensified for higher levels of moving mass velocity, the maximum values of design parameters would increase either. Moreover, the possibility of mass separation is shown to be more critical as the span number of the beam increases. This fact also violates the linear relation between the mass weight of the moving load and the associated design parameters, especially for high moving mass velocities. However, as the relaxation rate of the beam material increases, the load inertia effects as well as the possibility of moving mass separation reduces.

ELASTO-PLASTICITY ANALYSIS BASED ON COLLOCATION WITH THE MOVING LEAST SQUARE METHOD

A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis.

Determination of the natural frequencies of axially moving beams by the method of multiple scales

The natural frequencies of an axially moving beam were determined by using the method of multiple scales. The method of second-order multiple scales could be directly applied to the governing equation if the axial motion of the beam is assumed to be small. It can be concluded that the natural frequencies affected by the axial motion are proportional to the square of the velocity of the axially moving beam. The results obtained by the perturbation method were compared with those given with a numerical method and the comparison shows the correctness of the multiple-scale method if the velocity is rather small.

Large deformation simulations of geomaterials using moving particle semi-implicit method

Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.

A Size-Dependent Functionally Graded Higher Order Plate Analysis Based on Modified Couple Stress Theory and Moving Kriging Meshfree Method

A size-dependent computational approach for bending,free vibration and buckling analyses of isotropic and sandwich functionally graded(FG)microplates is in this study presented.We consider both shear deformation and small scale effects through the generalized higher order shear deformation theory and modified couple stress theory(MCST).The present model only retains a single material length scale parameter for capturing properly size effects.A rule of mixture is used to model material properties varying through the thickness of plates.The principle of virtual work is used to derive the discrete system equations which are approximated by moving Kriging interpolation(MKI)meshfree method.Numerical examples consider the inclusions of geometrical parameters,volume fraction,boundary conditions and material length scale parameter.Reliability and effectiveness of the present method are confirmed through numerical results.

Probability distribution of wind power volatility based on the moving average method and improved nonparametric kernel density estimation

In the process of large-scale,grid-connected wind power operations,it is important to establish an accurate probability distribution model for wind farm fluctuations.In this study,a wind power fluctuation modeling method is proposed based on the method of moving average and adaptive nonparametric kernel density estimation(NPKDE)method.Firstly,the method of moving average is used to reduce the fluctuation of the sampling wind power component,and the probability characteristics of the modeling are then determined based on the NPKDE.Secondly,the model is improved adaptively,and is then solved by using constraint-order optimization.The simulation results show that this method has a better accuracy and applicability compared with the modeling method based on traditional parameter estimation,and solves the local adaptation problem of traditional NPKDE.

A Comprehensive Method for Evaluating Precision of Transfer Alignment on a Moving Base

In this study, we propose the use of the Degree of Alignment(DOA) in engineering applications for evaluating the precision of and identifying the transfer alignment on a moving base. First, we derive the statistical formula on the basis of estimations. Next, we design a scheme for evaluating the transfer alignment on a moving base, for which the attitude error cannot be directly measured. Then, we build a mathematic estimation model and discuss Fixed Point Smoothing(FPS), Returns to Scale(RTS), Inverted Sequence Recursive Estimation(ISRE), and Kalman filter estimation methods, which can be used when evaluating alignment accuracy. Our theoretical calculations and simulated analyses show that the DOA reflects not only the alignment time and accuracy but also differences in the maneuver schemes, and is suitable for use as an integrated evaluation index. Furthermore, all four of these algorithms can be used to identify the transfer alignment and evaluate its accuracy. We recommend RTS in particular for engineering applications. Generalized DOAs should be calculated according to the tactical requirements.

Calculation of stratum surface principal curvature based on a moving least square method

With the east section of the Changji sag Zhunger Basin as a case study, both a principal curvature method and a moving least square method are elaborated. The moving least square method is introduced, for the first time, to fit a stratum surface. The results show that, using the same-degree base function, compared with a traditional least square method, the moving least square method can produce lower fitting errors, the fitting surface can describe the morphological characteristics of stratum surfaces more accurately and the principal curvature values vary within a wide range and may be more suitable for the prediction of the distribution of structural fractures. The moving least square method could be useful in curved surface fitting and stratum curvature analysis.

Moving-Particle Semi-Implicit Method for Vortex Patterns and Roll Damping of 2D Ship Sections

A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.

DYNAMIC RESPONSE OF PLATES DUE TO MOVING VEHICLES USING FINITE STRIP METHOD

Dynamic response of beam-like structures to moving vehicles has been extensively studied. However, the study on dynamic response of plates to moving vehicles has so far received but scant attention. A plate-vehicle strip for simulating the interaction between a rectangular plate and moving vehicles was described. For the portion of strips that are in direct contact with the moving vehicles, the plate-vehicle strips were employed. Conventional plate finite strips were used to model the portion of strips that are not directly under the action of moving vehicles. In the analysis, each moving vehicle is idealized as a one-foot dynamic system with tire unsprung mass and sprund mass interconnected by a spring and a dashpot. The numerical results obtained from the proposed method agree well with available results.

TRACKING AND RECONSTRUCTION METHODS FOR MOVING-INTERFACES

A simple,efficient and accurate high resolution method to tracking moving-interfaces-the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluation of six main efficient methods for interface reconstruction are made. Through strict numerical simulation,their characters,advantages and shortcomings are compared,analyzed and commended in particular.

A baroclinic typhoon model with a moving multi-nested grid and variational adjustment initialisationI. Numerical method

A baroclinic typhoon model with a moving multi--nested grid is applied in marine environmental forecasts. This paper describes the numerical methods of the model including governing equations, finite differencing, split scheme and time integration.

Chebyshev finite spectral method with extended moving grids

A Chebyshev finite spectral method on non-uniform meshes is proposed. An equidistribution scheme for two types of extended moving grids is used to generate grids. One type is designed to provide better resolution for the wave surface, and the other type is for highly variable gradients. The method has high-order accuracy because of the use of the Chebyshev polynomial as the basis function. The polynomial is used to interpolate the values between the two non-uniform meshes from a previous time step to the current time step. To attain high accuracy in the time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme is used. To avoid numerical oscillations caused by the dispersion term in the Korteweg-de Vries （KdV） equation, a numerical technique on non-uniform meshes is introduced. The proposed numerical scheme is validated by the applications to the Burgers equation （nonlinear convectiondiffusion problems） and the KdV equation （single solitary and 2-solitary wave problems）, where analytical solutions are available for comparisons. Numerical results agree very well with the corresponding analytical solutions in all cases.

THE UPWIND FINITE DIFFERENCE METHOD FOR MOVING BOUNDARY VALUE PROBLEM OF COUPLED SYSTEM

Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.