• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    On the Boundary Condition and Related Instability in the Smoothed Particle Hydrodynamics

    2019-11-07 02:58:20ChongYe葉翀PhilipeMotaJinLi李瑾KaiLin林愷andWeiLiangQian錢衛(wèi)良
    Communications in Theoretical Physics 2019年11期

    Chong Ye (葉翀), Philipe Mota, Jin Li (李瑾),Kai Lin (林愷), and Wei-Liang Qian (錢衛(wèi)良)

    1College of Physics, Chongqing University, Chongqing 401331, China

    2Centro Brasileiro de Pesquisas Físicas, 22290-180, Rio de Janeiro, RJ, Brazil

    3Institute of Geophysics and Geoinformatics, China University of Geosciences, Wuhan 430074, China

    4Escola de Engenharia de Lorena, Universidade de S?o Paulo, 12602-810, Lorena, SP, Brazil

    5Faculdade de Engenharia de Guaratinguetá, Universidade Estadual Paulista, 12516-410, Guaratinguetá, SP, Brazil

    6School of Physical Science and Technology, Yangzhou University, Yangzhou 225002, China

    Abstract In this work,we explore various relevant aspects of the Smoothed Particle Hydrodynamics regarding Burger’s equation.The stability, precision, and efficiency of the algorithm are investigated in terms of different implementations.In particular, we argue that the boundary condition plays an essential role in the stability of numerical implementation.Besides, the issue is shown to be closely associated with the initial particle distribution and the interpolation scheme.Among others, we introduce an interpolation scheme termed symmetrized finite particle method.The main advantage of the scheme is that its implementation does not involve any derivative of the kernel function.Concerning the equation of motion,the calculations are carried out using two distinct scenarios,where the particles are chosen to be either stationary or dynamically evolved.The obtained results are compared with those obtained by using the standard finite difference method for spatial derivatives.Our numerical results indicate subtle differences between different schemes regarding the choice of boundary condition.In particular, a novel type of instability is observed where the regular distribution is compromised as the particles start to traverse each other.Implications and further discussions of the present study are also addressed.

    Key words:hydrodynamic model, SPH algorithm, pair instability

    1 Introduction

    The smoothed particle hydrodynamics (SPH)[1?9]is a meshfree, Lagrangian method for the partial different equations.Distinct from traditional approaches such as the finite difference method (FDM) or finite element method,the SPH approximates a partial differential equation in terms of ordinary differential equations of freely moving interpolation points, known as SPH particles.Subsequently, the infinite number of degrees of freedom of the original system is represented by a finite number of SPH particles, aiming to describe any arbitrary configuration of the continuum medium.As one of the most important applications of the partial differential equation is to describe the dynamics of fluid, SPH method naturally finds extensive applications in the context of hydrodynamics, as inferred by the terminology.Indeed,due to the significance of hydrodynamics, as well as recent advance in computing power, the SPH method is widely employed in many distinct areas of science and engineering.Among traditional applications to fluid dynamics,the SPH method has been employed to investigate a great variety of topics.In astrophysics, the SPH is implemented for the study of galaxy formation, supernova explosion.[7,9?11]In solid state mechanics,[12]the method has been applied to problems such as metal forming,[13]fracture and fragmentation,[14]and shock wave propagation in solids.[15]Moreover, concerning the high energy nuclear physics, the hydrodynamic approach of heavy-ion nuclear collisions plays an essential role associated with RHIC and LHC.[16?18]The latter further inspires onging investigations of fluid/gravity duality,[19?22]one realization of the holographic principle.[23,24]Numerical simulations based on implementations of SPH algorithm[25?29]have shown to give reasonable description on various experimental observables.The algorithm is also employed in chemistry where a variety of chemical evolutions simulations has been implemented in terms of the SPH method.

    The SPH was initially suggested to investigate astrophysical problems by Gingold and Monaghan[2]and by Lucy.[1]In the standard implementation, each SPH particle is assigned a finite size,h, over the extension of which the relevant physical properties are smoothed in terms of a kernel function.By kernel function, the contribution of each particle is weighted regarding their respective distance to the coordinate position of interest.As a consequence,a physical quantity at a given position is evaluated by summing over all the contributions from all the relevant SPH particles in the neighborhood within reach of the kernel function.In terms of the above interpolation regarding SPH particles, a partial differential equation is transformed into a system of ordinary differential equations of the SPH degrees of freedom.

    The SPH algorithm is a meshfree approach, which implies several essential features.First,the method is suited to handle problems with sophisticated boundary conditions, inclusively, scenarios with free boundary.The algorithm can achieve exact advection mostly in a trivial manner, which is a rather difficult task for any Eulerian grid-based approaches.Meanwhile, the absence of the structured computational mesh, by and large, simplifies model implementation.In particular, the spatial derivatives are implemented in a flexible and straightforward fashion, independent of any grid-based configuration.Besides, one of the most remarkable features of the SPH method is that it naturally guarantees the conservation law.In other words, independent of any specific application, the algorithm rigorously conserves quantities such as mass, energy, momentum, and angular momentum by construction.Also, a certain degree of re-meshing takes place automatically as time evolves, as a consequence of the kernel function.It is distinct from other Lagrangian grid-based methods where the grid can sometimes become arbitrarily disordered in time.[30]Moreover, when appropriately configured, it is understood that the computational cost of numerical simulations per SPH particle is usually less than that per cell of grid-based simulations.This is because the resolution of the method follows naturally the mass flow rather than the volume as in the case of most grid-based techniques.It happens as each SPH particle carries a fixed amount of mass and dislocate itself according to the flow velocity, and subsequently, the density only goes up related to an increase in the local matter concentration.[30]To be more specific,it is the case especially when the metric of interest is associated with the relevant density, and thus an automatical refinement on density is achieved.Furthermore, the SPH method is highly flexible and can be combined with other algorithms straight forwardly.[10]

    However, the SPH formulation also faces a few challenges.Several notable problems potentially hamper further exploiation of the method, such as accuracy, tensile instability,[31?33]pairing instability,[34]the kernel as well as particle consistency,[8,35?38]and artificial viscosity.[39]From a physical viewpoint, tensile instability is triggered when SPH particles become mutually attractive due to,by and large,the presence of negative pressure.The eventual instability manifested regarding particle distribution,known as tensile instability,can be viewed as equivalent to the mesh distortion in a mesh-based scenario, is typically found in simulations of elastic material under tension.One possible approach to deal with such instability associated with irregular particle distribution is to sacrifice exact conservation of the algorithm by a slim amount.[40]Besides,other approaches were proposed to handle the instability by improving the estimation of spatial gradients, such as the corrected derivatives method.[41]As the number of neighbors of an individual SPH particle increases,the ratio of the particle spacing to the smoothing length decreases.However, for standard “Bell-shaped” or “Gaussian-like”kernels, their first derivative, and equivalently the force between SPH particles, vanishes at the origin.Consequently,particle pairing occurs.This problem is relatively benign, as it usually does not lead to any severe consequence.It can mostly be avoided by just using an adequate initial particle distribution with a smaller number of neighbors or adopting a short-range kernel function.Concerning kernel and particle consistency, much effort has been devoted to restore the consistency and subsequently improve accuracy.Monaghan derived a symmetrization formulation for the spatial gradient term[3]consistent with the variational approach,which is shown to provide better accuracy.Randles and Libersky proposed a normalization scheme for the density approximation.[35]The concept is further developed by Chenet al.[33,36,42]in terms of a corrective smoothed particle method (CSPM) as well as by Liuet al.[37]in terms of a finite particle method (FPM).For some physical phenomena such as shockwave, the entropy is increased as a physical consequence.In fact, unlike in a Eulerian scheme where numerical dissipation is an intrinsic feature related to the discretization process and mostly unavoidable, diffusion terms in SPH must be inserted explicitly, and as a result, their physical content is transparent.In this regard, it has been argued that it is essential to supply an adequate amount of artificial viscosity which works for not only scenarios where artificial viscosity is required but also for other cases where its presence ought to be minimized.[43,44]

    As a matter of fact, there does not exist a single numerical method which performs satisfactorily well at all different aspects.Therefore, in practice, the choice of the best numerical approach can often depend entirely on which feature of the problem is more relevant to the present study.Nonetheless, stability is usually the primary concern of a numerical scheme, while precision and efficiency are two essential performance metrics in most cases.Our present study thus aims to provide a quantitative analysis of these aspects of the SPH method.In particular, we explore the role of boundary condition, as well as its relationship with initial particle distribution and interpolation scheme.Regarding the specific problem we choose to investigate, the equation of motion is implemented using two distinct scenarios, where the SPH particles are selected to be either stationary or allowed to evolve dynamically in time.The first scenario is shown to be an excellent indicator of the efficiency and accuracy of the individual interpolation algorithm.The second scenario, on the other hand, resides on a more general basis and is subsequently found to be sensitively associated with the stability of different implementations.

    The paper is organized as follows.In the following section, we briefly discussed the main features of different implementations of the SPH algorithms involves in the present study, namely, the standard SPH, CSPM, FPM,and the symmetrized finite particle method (SFPM).The calculations are carried out in Secs.3 and 4 where Burger’s equation is chosen to be the object of the numerical simulations.The results are then compared against the analytic solutions.Two specific scenarios for the equation of motion are made use of with either stationary and dynamical SPH particles.In Sec.3 we show the results concerning the first scenario, where we mainly focus on the precision and efficiency of different interpolation schemes.The second scenario is studied in Sec.4, which is mostly related to the stability of the algorithm regarding the boundary conditions.We discuss how the boundary condition affects the stability of the numerical calculations for different schemes.Also, a type of novel instability is observed where SPH particles are found to traverse each other,which, in turn, undermine the regular particle distribution.The last section is devoted to further discussions and concluding remarks.

    2 The SPH Method

    In the present study, we mainly focus on four different SPH interpolation schemes, namely, standard SPH,CSPM, FPM, and SFPM.In what follows, what briefly review these schemes concerning their derivations.For the spatial distribution of a given physical quantity,f(x),all the above schemes can be viewed as specific implementations to interpolate the distribution by using a finite number of SPH particles as well as the kernel functionW.The latter is assumed to be an even function and its normalization reads

    The Taylor series expansion off(x) reads

    We proceed by multiplying both sides by the kernel functionW(x ?xa) and integrating overx, we have

    As one keeps the first term on the r.h.s.of the above equation and discard the higher order ones, one obtains

    By substituting the normalization condition Eq.(1),we obtain the standard SPH formula, namely,

    We note that the particle approximation is utilized,where a summation enumerating all SPH particle is used to approximate the integral.Hereνbis the “mass” of theb-th particle, usually it is a given value assign to the particle at the beginning of the algorithm.ρbis the density of the SPH particles, iffitself is the density, one has

    In practice, the ratiofb/ρbcan be determined by the physical properties of the system in question.For instance, iffrepresents pressure andρis the number density, their ratio might be determined by the equation of state.[17]Otherwise,fbandρbcan be that evaluated for theb-th particle at the last time step.The first and second order derivatives can be obtained by carrying out spatial derivative of Eq.(5), and are found to be

    where symmetric formulism[3]has been employed.

    One might be aware of the fact that in the particle approximation of Eq.(1),the kernel function is not precisely identical to the Diracδ-function, and instead employ the following modified particle approximation in the place of Eq.(4),

    One finds

    Equation (10) is known as CSPM in literature, which can be shown to preserve the zeroth-order kernel and particle consistency.[8]The first-order and second-order derivatives can be derived by retaining until the second order term on the r.h.s.of Eq.(2), multiplying both sides byW′andW′′, and making use of the symmetry of these functions.One finds

    In fact, one may proceed further by generalizing the above arguments to higher order in a slightly more consistent fashion.This is achieved by consistently expanding the r.h.s.of Eq.(2) to the second order, namely,

    The present case is different from that of CSPM scheme, since now Eqs.(12)–(13) furnishes a system of coupled equations, as follows

    By inverting the above matrix equation, one obtains for the desired expressions for the FPM scheme, which read

    with

    The FPM scheme is understood to possess the zeroth and first-order particle consistency.[8]

    It is worth mentioning that in deriving Eq.(13), one has the freedom to utilize any basis function.Another particular choice, which is meaningful by its own right, is to make use of(x?xa)W(x?xa)and(x?xa)2W(x?xa)in the place ofW′(x ?xa) andW′(x ?xa).The only difference is that the definitions in Eq.(15) are replaced by

    We note that the alternative basis functions maintain the symmetries (even or odd with respect to the origin) of the integral kernel.Meanwhile, the formalism has avoided using the derivatives of the kernel function,which turns out to be quite advantageous as shown below.This last scheme was initially proposed in Ref.[45]and will be referred to as the symmetrized finite particle method(SFPM) in the present work.

    To sum up, Eqs.(5)–(7), Eq.(10), Eqs.(15)–(16) and Eqs.(16)–(17) give the specific formalism for standard SPH, CSPM, FPM and SFPM.In the remainder of the paper, these differences will be put into practice in the following two sections.

    3 Burgers’ Equation with Stationary SPH Particles:on Precision and Efficiency of Different Schemes

    In this section, we focus ourselves on the static aspect of the algorithm regarding the precision as well as the efficiency of different schemes.By and large, the precision of a given scheme can be measured in terms of kernel consistency and particle consistency.The former is determined by then-th moment of the kernel functions, namely,

    forn=0,1,2,...The above conditions are rigorous when the interpolation is perfect.When it is satisfied for a given ordern, the corresponding SPH approximation is said to possessn-th order orCnconsistency.

    The particle consistency, on the other hand, is related to the fact that in practice, the integral in Eq.(18) is implemented by a summation.As a result, one requires

    As mentioned above, the present section is primarily dedicated to investigating to what degree the SPH algorithm is capable of reproducing a well-defined distribution of certain physical quantity in terms of the SPH particles.In this context, the precision and efficiency are affected by various factors, such as the specific form of the kernel function, the number of particles, and the interpolation scheme.In what follows, we carry out numerical simulations while tuning parameters regarding these aspects.We note that, for the time being, we will temporarily discard another vital factor, the stability of the algorithm.In other words, we only concentrate on the precision and efficiency of specific problems where the stability is guaranteed or a minor issue.In this regard, we study the following one-dimensional Burgers’ equation:

    whereγis the viscosity coefficient.The reason for our choice of Burgers’ equation is as follows.First, the equation possesses an analytic solution, corresponding to smooth initial conditions whereuvanishes on the boundary.Secondly, usually, the physical interpretation of the equation is the same as the Navier-Stokes equation.In other words,the Burgers’equation determines the temporal evolution of,u,the velocity of the fluid.Therefore,the physical system is governed by two equations for the fluid velocityu(x,t) and matter densityρ(x,t), respectively.The latter is nothing but the continuity equation of the fluid, which takes the following form in one-dimensional space,

    However, it is interesting to point out that one may simply viewuas an intensive physical quantity, and solve Eq.(20) alone independent of Eq.(21).This is possible since, for the specific case of the Burgers’ equation, the equation ofudoes not depend onρ.Therefore, when implementing the SPH scheme, we choose to assign the SPH particles to fixed spatial coordinates.Here,uaof thea-th particle is a function of time, determined by the r.h.s.of Eq.(20), while the coordinates of the SPH particles remain unchanged during the simulation.Lastly,the boundary condition of the present problem is that the distribution attains zero at both end points, namely,

    But since the SPH particle do not dislocate in time, for the present section, one does not need to explicitly implement Eq.(22) into the code.As shown in the next section, the latter is closely related to the issue of stability.For simplicity,throughout our calculations,we employ the fifth order spline kernel function.[46]We also note that the present application involves fixed boundary condition,while on the other hand, SPH algorithm is most advantageous for problems with generic free boundary condition,the resultant precision and efficiency of SPH is not expected to be superior to the outcome of traditional FDM.In what follows,the different schemes discussed in the previous section will be applied by using the above setup.

    In Fig.1, we show the calculated temporal evolutions of the distribution in comparison with those obtained by using the FDM as well as the analytic solutions.It is found that most methods give consistent and reasonable results for different viscosity coefficients.Among all the interpolations schemes, the standard SPH shows the most substantial deviations from the analytic results.To precisely identify the precision and efficiency of different schemes,we numerically evaluate the deviations from analytic results.In Table 1, the averaged relative precision and efficiency for different numbers of SPH particles or grid numbers (for FDM) are summarized.In our calculations, the particles are distributed uniformly, and the size of SPH particles,h, varies according to the total number of SPH particles.To be specific, we chooseh=?xwhere ?xis the distance between neighbor particles.In the calculations, the time interval is taken to be ?t=h/500, and we considerγ=1.0.The precision is estimated for the instantt=1.0 while regarding an average of the relative deviation

    whereNis the total number of selected spatial points where the quality of the interpolation is evaluated by comparing to analytic result.In our calculations,we have used seven evenly distributed spatial points for the evaluation.As the efficiency, one additionally takes into account the CPU time consumption, ?t.The efficiency of a given method is defined as follows

    which is measured with respect to that of the standard SPH method.

    Fig.1 (Color online) The numerical solutions of the Burgers’ equation in comparison with the analytic ones.The calculations are carried out by using stationary SPH particles for standard SPH, CSPM, FPM, FDM, and SFPM,respectively.The results are presented for different values of the viscosity coefficients γ.

    Table 1 The relative precision and efficiency of different methods.The first row indicates the number of SPH particles or grids used in the calculations.

    It is observed in Table 1,for all the methods,the calculatedPmonotonically decrease with increasing SPH particles.This is as expected, which indicates that the precision improves as the particle number increases.Also, the efficiency,which is defined as an overall evaluation of speed and accuracy, improves as the SPH particle increase.The modified SPH methods included in this study; namely,CSPM,FPM,and SFPM show significantly better performance than that of the standard SPH method.It is found that as the number of SPH particle further increases, the improvement of efficiency gradually slows down.This happens due to that the precision of the interpolation starts to be saturated while the computational time increases linearly.It is further expected for the higher dimensional case since the computation cost increases much faster with increasing particle number, the efficiency of the methods may even start to decrease as the particle number becomes significantly large.

    4 Burgers’ Equation with Dynamical SPH Particles:on the Stability of the Algorithm

    This section is devoted to the stability analysis of different schemes regarding the boundary condition and initial particle distribution.When not appropriately configured, a numerical algorithm might be plagued by instability.Here we investigate particularly how the boundary condition may affect the numerical instability.As discussed in the introduction, there are at least three different types of characteristics of the SPH algorithm closely related to stability issue, namely, self-regularization, pairing instability, and tensile instability.All of them are related to the particle distribution and the specific form of the kernel function.SPH can automatically maintain a reasonable particle arrangement[10,30]during the course of the temporal evolution.Owing to the form of the kernel function, mutual repelling forces are exerted between neighboring SPH particles,which automatically drives the system towards a regular particle distribution.Intuitively,it corresponds to “re-meshing” procedure that is utilized in other Lagrangian mesh methods.As the SPH particles come closer, paring instability takes place as soon as the repelling force start to decrease with decreasing neighbor distance.Furthermore, under certain conditions[31,32]tensile instability takes place and disrrupts the particles distribution.In what follows, we attempt to configure an initial SPH distribution where the distance of neighboring SPH particles is randomized but within the range where particles repelling each other and the strength decreases with increasing distance.However,our calculation does not reveal any automatical re-meshing.Moreover,we show that boundary condition may introduce small perturbation which is accumulated over time and turned into oscillations.Such instability is observed to be accompanied by particle crossing.In order to carry out the analysis, one allows the SPH particle to move according to the continuity equation, namely, Eq.(21).In other words, we will solve Eqs.(20) and (21) simultaneously, as in most SPH implementations.The numerical results are presented in Figs.2?5.

    Fig.2 (Color online) The numerical solutions of the Burgers’ equation in comparison with the analytic ones.The calculations are carried out by using different values of the viscosity coefficients for standard SPH, CSPM, FPM, FDM,and SFPM, respectively.The boundary condition is implemented by using fixed SPH particles.

    Fig.3 (Color online) The numerical solutions of the Burgers’ equation in comparison with the analytic ones.The calculations are carried out by using different values of the viscosity coefficients for standard SPH, CSPM, FPM, FDM,and SFPM, respectively.The boundary condition is implemented by adding extra SPH particles beyond the physical boundary.

    Fig.4 (Color online) The numerical solutions of the Burgers’ equation in comparison with the analytic ones.The calculations are carried out by using different values of the viscosity coefficients for standard SPH, CSPM, FPM, FDM,and SFPM, respectively.The boundary condition is implemented by adding mirrored imaginary SPH particles beyond the boundary.

    A critical aspect of the present study is to show that the boundary condition plays an essential role in the SPH algorithm.From its appearance, the boundary condition Eq.(22)of the problem looks“benign”.However,it turns out that the numerical stability sensitively depends on its implementation.We study three types of implementation for the boundary condition.First, we freeze the coordinates of the SPH particles on the boundary.Secondly, we extend the system beyond the boundary by adding more SPH particles and let those additional particles evolve freely.To be specific, we do not explicitly exert any additional condition for the most outward particles.Lastly,we implement mirroring boundary condition, by adding imaginary particles outside of the boundary to ensure the validity of Eq.(22).To be more precise, the imaginary particles are mirror images; namely, their positions are mirrored with respect to the location of the boundary,and their velocities are inverted regarding their counterparts.The obtained results are shown in Figs.2?4.

    The calculated results indicate that it is essential to appropriately implement the boundary condition, determined by Eq.(22).As shown in Fig.2, the calculations with fixed boundary SPH particles start to fail as the coefficient of viscosity becomes significantly small.Asγdecreases, it takes a longer time for the system to evolve,and subsequently,small deviations are accumulated for an extended period,and the results demonstrate more significant discrepancy.It is observed that most schemes show significant deviations forγ <0.1 and fort >0.75, where the standard SPH fails dramatically.While one adds additional SPH particles beyond the boundary of the physical problem and apply free boundary condition for those particles,the results show a very similar feature asγbecomes sufficiently small.This can be observed in Fig.3.In Fig.4,the results are shown by employing mirrored boundary condition.We did not draw the mirrored SPH particles beyond the boundary as they are symmetrically positioned regarding their counterparts.In general, the results are significantly improved.It is found that the implemented boundary condition gives a significant impact on the resulting temporal evolution.For CSPM, the results are more stable and agree reasonably with the analytic ones for all different viscosity coefficients.Although, in some region,the numerical results are still found to deviate from the analytic ones.For FPM and SFPM, small oscillations are observed when the time becomes more substantial as SPH particles start to pile up at the right boundaryx=1.For the standard SPH method,the results are less satisfactory compared to the other approaches.We understand that it is mostly related to the fact that the standard SPH method suffers a severe problem with kernel and particle consistency.

    In order to investigate the instability occurs near the rightmost boundary, we carry out further calculations focusing the region near the vicinity ofx=1.To be specific, the numerical calculations are carried out for the range ofx ∈(0,2) with the boundary beingx=0 andx=2.The imaginary SPH particles beyond the original boundary are now treated as real particles and evolve according to the Burgers’ equation.In this context, the particles on both sides are moving toward the center atx=1.The numerical results are shown in Fig.5.The two plots on the first row are the trajectories the positions and velocities of SPH particles, shown by curves in different colors.A total of 125 SPH particles are initially distributed uniformly, whereh=0.2 and the size of the time step dt=0.02.The plots on the second as well as third rows are the same as the first row, except for the following difference.For the second row,the SPH particle distribution is also uniform initially, but a total number of 250 particles,withh=0.1 and dt=0.01.For the third row,regarding the initial condition for the case of the second row, the positions and velocity of SPH particles are perturbed by a small amount at the initial time.While the calculations are carried out for the rangex ∈(0,2)with the boundary are defined by the pointsx=0 and 2,the results are only shown in the vicinity region ofx ~1.It is found that on the first row of Fig.5, particles simple travel across each other after they come close.This is a troublesome result, especially when comparing to the analytic results, the fluid on either side of the boundary should never traversex=1.Therefore,some fundamental problem regarding interpolation is implied.As will be discussed below, the stability mentioned above in Figs.2?4 is likely related to the present phenomenon.On the second and third rows, particles bound back from each other due to the repelling force between them, as intended.As discussed before, if particles maintain at an appropriate distance between each other, particle auto-remeshing is expected.In the plots on the second row of Fig.5, this is manifested, to a certain degree, in terms of the fact that all the SPH particles gradually decelerate and come to a full stop simultaneously att ~0.6.Further, in the plots on the third row,it is shown that the particles also bounce back from each other, owing to the auto-remeshing mechanism.However, our results indicate that auto-remeshing does not always work in the present context, in particular, where the fluid becomes piled up atx=1 and is forced to squeeze into a small region, as governed by the Burgers’ equation.In this case, the results become quite sensitive to the specific parameter used in the calculations,as demonstrated by the plots on the first and second rows.When the auto-remeshing fails, as shown in the first row,the SPH particles are found to traverse each other when they stay too close to each other.Judging from the form of the kernel function, when ?x < h, the repelling between two neighboring SPH particles decrease with decreasing distance.Therefore, as the neighboring particles become too close,as it inevitably happens in the case of the Burgers’ equation, the mechanism of particle auto-remeshing may fail.In fact, we observe that the resultant particle traversing occurs for a broad choice of parameters in the present scenario of the Burgers’ equation.We note that such a phenomenon is observed to accompany the instability discussed previously concerning Figs.2?4 closely.

    5 Concluding Remarks

    To summarize, in this paper, we carry out an analysis regarding the stability, precision, and efficiency of the SPH method.In particular, various aspects of the algorithm are explored such as the boundary condition,initial particle distribution,and interpolation scheme.Besides,a generalized FPM scheme, the symmetrized finite particle method,is employed in our study.The main advantage of the SFPM is that its implementation does not explicitly involve the derivatives of the kernel function.We studied two different scenarios for the implementation of the equation of motion, where the SPH particles are either stationary or dynamically evolved in time.In the first scenario, the obtained results are also compared with those obtained by using the FDM as well as analytic ones.All numerical schemes provide reasonable results in the case.Among others, it is indicated that CSPM and SFPM provides an overall better performance regarding precision as well as stability.

    For the second scenario, the positions and velocity of the SPH particles are determined by the equations of motion.Though it is more meaningful to apply this scheme to topics with dynamical boundary conditions,in the present study, it is adopted to the same problem of the Burgers’equation.It is subsequently found that the boundary condition implementation plays a crucial role in the success of the algorithm.Among the three implementations, the mirrored boundary condition leads to mostly satisfactory results.Still,instability in terms of oscillations is observed when particle crossing starts to take place.Even for the case of mirrored boundary condition,increasing oscillation is observed for FPM and SFPM methods in regions with more significant particle density.Since the SFPM does not explicitly involve any derivative of the kernel function, the observed instability is likely of a distinct nature.This is because the sufficient condition of tensile instability involves the second order of the kernel function when derived regarding small perturbations.[31]One crucial feature of observed instability is that the oscillations take place when SPH particles start to cross each other.We will try to address the issue of particle crossing and the related instability in future work.

    Acknowledgment

    We gratefully acknowledge the financial support from Brazilian funding agencies Funda?c?o de Amparo à Pesquisa do Estado de S?o Paulo (FAPESP), Conselho Nacional de Desenvolvimento Científico e Tecnológico(CNPq), Coordena?c?o de Aperfei?coamento de Pessoal de Nível Superior (CAPES), National Natural Science Foundation of China(NNSFC)under Grant Nos.11805166 and 11873001, and Nature Science Fund of Chongqing under Grant No.cstc2018jcyjAX0767.

    亚洲欧美精品综合一区二区三区| 露出奶头的视频| 无遮挡黄片免费观看| 亚洲成人免费电影在线观看| 久久久久久久久免费视频了| 黄色片一级片一级黄色片| 男人操女人黄网站| 国产精品久久电影中文字幕| 精品卡一卡二卡四卡免费| 免费在线观看视频国产中文字幕亚洲| 欧美激情久久久久久爽电影 | 午夜福利高清视频| 午夜福利欧美成人| 国产91精品成人一区二区三区| 精品国产美女av久久久久小说| 女生性感内裤真人,穿戴方法视频| 国内毛片毛片毛片毛片毛片| 久久性视频一级片| 操出白浆在线播放| 中文字幕最新亚洲高清| 黑人操中国人逼视频| 欧美性长视频在线观看| 男人舔女人的私密视频| 日韩 欧美 亚洲 中文字幕| 两个人免费观看高清视频| 国产一区二区三区综合在线观看| 满18在线观看网站| 19禁男女啪啪无遮挡网站| 两个人视频免费观看高清| 成人手机av| 亚洲一区高清亚洲精品| 精品第一国产精品| 在线十欧美十亚洲十日本专区| 妹子高潮喷水视频| 一a级毛片在线观看| 亚洲成av片中文字幕在线观看| 亚洲性夜色夜夜综合| 一个人观看的视频www高清免费观看 | 国内久久婷婷六月综合欲色啪| 丝袜美足系列| 亚洲熟女毛片儿| 日本五十路高清| 国产欧美日韩一区二区三| 露出奶头的视频| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美三级三区| 一个人观看的视频www高清免费观看 | 国产精品一区二区免费欧美| 亚洲专区字幕在线| 在线观看舔阴道视频| 黄色成人免费大全| 国产欧美日韩一区二区三| 亚洲九九香蕉| 日本一区二区免费在线视频| 久久午夜综合久久蜜桃| 天天一区二区日本电影三级 | 黑丝袜美女国产一区| 又大又爽又粗| 免费久久久久久久精品成人欧美视频| 午夜日韩欧美国产| 精品国产超薄肉色丝袜足j| 午夜福利18| 亚洲av美国av| 亚洲精华国产精华精| 国产欧美日韩综合在线一区二区| 色综合欧美亚洲国产小说| 亚洲色图综合在线观看| 亚洲电影在线观看av| 亚洲欧美日韩无卡精品| 琪琪午夜伦伦电影理论片6080| 亚洲中文av在线| 久久亚洲真实| 久久久久国产一级毛片高清牌| 欧美最黄视频在线播放免费| 国产亚洲av高清不卡| 午夜福利18| 老汉色av国产亚洲站长工具| 一边摸一边抽搐一进一小说| 欧美色视频一区免费| 国产精品,欧美在线| 国产精品野战在线观看| 中文字幕久久专区| 亚洲欧美激情综合另类| 国产一区二区三区在线臀色熟女| 亚洲久久久久久中文字幕| 国产精品久久视频播放| 欧美日韩中文字幕国产精品一区二区三区| 欧美一区二区国产精品久久精品| 一区二区三区激情视频| 亚洲最大成人手机在线| 少妇的逼好多水| 亚洲精品色激情综合| АⅤ资源中文在线天堂| 高清日韩中文字幕在线| 日韩一本色道免费dvd| 成人毛片a级毛片在线播放| 不卡一级毛片| 少妇高潮的动态图| 欧美激情国产日韩精品一区| 天堂av国产一区二区熟女人妻| 九色国产91popny在线| 亚洲久久久久久中文字幕| 男人和女人高潮做爰伦理| 日韩欧美精品免费久久| 一个人观看的视频www高清免费观看| 国产精品女同一区二区软件 | 岛国在线免费视频观看| 婷婷精品国产亚洲av在线| 日本-黄色视频高清免费观看| 国产精品乱码一区二三区的特点| 国产在线精品亚洲第一网站| 特大巨黑吊av在线直播| 午夜福利视频1000在线观看| 性色avwww在线观看| 亚洲无线在线观看| 国产精品三级大全| 国产高清激情床上av| 3wmmmm亚洲av在线观看| 精品人妻1区二区| 久久热精品热| 人妻夜夜爽99麻豆av| 精品一区二区三区人妻视频| 一级黄色大片毛片| 看免费成人av毛片| 国产亚洲精品av在线| 午夜精品一区二区三区免费看| 亚洲不卡免费看| 老司机福利观看| 简卡轻食公司| 男插女下体视频免费在线播放| 校园人妻丝袜中文字幕| 色哟哟·www| 免费电影在线观看免费观看| 国产 一区精品| 日本与韩国留学比较| 久久精品人妻少妇| 日韩国内少妇激情av| 少妇熟女aⅴ在线视频| 日韩,欧美,国产一区二区三区 | 在线a可以看的网站| 久久久午夜欧美精品| 搡女人真爽免费视频火全软件 | 久久热精品热| 国产一区二区三区av在线 | 色综合婷婷激情| 欧美成人免费av一区二区三区| 免费看日本二区| 97超级碰碰碰精品色视频在线观看| 中国美女看黄片| 99国产精品一区二区蜜桃av| 3wmmmm亚洲av在线观看| 亚洲国产精品sss在线观看| 男人舔奶头视频| 一进一出抽搐gif免费好疼| 色哟哟哟哟哟哟| 精品一区二区三区人妻视频| 三级国产精品欧美在线观看| 麻豆成人av在线观看| 亚洲专区国产一区二区| 午夜视频国产福利| 国产精品不卡视频一区二区| 窝窝影院91人妻| 男人狂女人下面高潮的视频| 久久久国产成人免费| 欧美不卡视频在线免费观看| 亚洲av免费高清在线观看| 午夜精品在线福利| 最后的刺客免费高清国语| 听说在线观看完整版免费高清| 91av网一区二区| 五月伊人婷婷丁香| 在线国产一区二区在线| 亚洲午夜理论影院| 国产久久久一区二区三区| 制服丝袜大香蕉在线| 国产私拍福利视频在线观看| 日韩欧美国产一区二区入口| 99久久精品一区二区三区| 成人精品一区二区免费| 国产亚洲精品av在线| 97超级碰碰碰精品色视频在线观看| 精品免费久久久久久久清纯| 一个人观看的视频www高清免费观看| 久久午夜福利片| 一区二区三区免费毛片| 精品久久久久久久末码| 久久国内精品自在自线图片| 可以在线观看的亚洲视频| 欧美激情在线99| 亚洲aⅴ乱码一区二区在线播放| 97热精品久久久久久| 成人av在线播放网站| 国产美女午夜福利| 韩国av在线不卡| 成人无遮挡网站| 夜夜爽天天搞| 麻豆国产av国片精品| 麻豆av噜噜一区二区三区| 国产久久久一区二区三区| 精品久久久噜噜| 看黄色毛片网站| 婷婷亚洲欧美| 少妇高潮的动态图| 国产老妇女一区| 国产高潮美女av| 22中文网久久字幕| 国产av在哪里看| 午夜福利成人在线免费观看| 久久久精品大字幕| 日本欧美国产在线视频| 国产美女午夜福利| 最近最新免费中文字幕在线| 97热精品久久久久久| 精品人妻熟女av久视频| 在线免费十八禁| 精品人妻熟女av久视频| 偷拍熟女少妇极品色| av在线天堂中文字幕| 免费av毛片视频| 两个人的视频大全免费| 成人美女网站在线观看视频| 又爽又黄无遮挡网站| 欧美丝袜亚洲另类 | 日韩强制内射视频| 欧美一区二区国产精品久久精品| 色播亚洲综合网| 久久久国产成人免费| 久久精品国产亚洲网站| 91麻豆av在线| 久久午夜福利片| 久久午夜亚洲精品久久| 日本免费a在线| 观看美女的网站| 女人十人毛片免费观看3o分钟| 一本一本综合久久| 免费看av在线观看网站| 中文在线观看免费www的网站| 午夜精品一区二区三区免费看| 99精品在免费线老司机午夜| 国产亚洲精品久久久久久毛片| 日本 av在线| 制服丝袜大香蕉在线| 国产单亲对白刺激| 十八禁国产超污无遮挡网站| 欧美不卡视频在线免费观看| 亚洲欧美日韩无卡精品| 日韩欧美三级三区| 久久九九热精品免费| 中文字幕人妻熟人妻熟丝袜美| 在线a可以看的网站| av国产免费在线观看| 欧美性感艳星| 99国产极品粉嫩在线观看| 国产不卡一卡二| 国产精品精品国产色婷婷| 亚洲aⅴ乱码一区二区在线播放| 俄罗斯特黄特色一大片| 在线观看美女被高潮喷水网站| 有码 亚洲区| 免费观看的影片在线观看| 99热这里只有是精品在线观看| 亚洲久久久久久中文字幕| 欧美最黄视频在线播放免费| 国产精品一区二区免费欧美| 日韩 亚洲 欧美在线| 性欧美人与动物交配| 国产aⅴ精品一区二区三区波| 国产在线精品亚洲第一网站| 国产又黄又爽又无遮挡在线| 99久久精品一区二区三区| 联通29元200g的流量卡| 国产在线男女| 一级黄色大片毛片| 国产精品福利在线免费观看| 一级a爱片免费观看的视频| 波多野结衣高清无吗| 日韩中字成人| 亚洲精品日韩av片在线观看| 少妇人妻一区二区三区视频| 精品一区二区免费观看| 美女免费视频网站| 久久中文看片网| 欧美性猛交黑人性爽| 乱人视频在线观看| 亚洲精品成人久久久久久| 久久精品人妻少妇| 午夜精品一区二区三区免费看| 国产成人一区二区在线| 成人亚洲精品av一区二区| 一边摸一边抽搐一进一小说| 免费观看人在逋| 久久午夜福利片| 蜜桃亚洲精品一区二区三区| 美女cb高潮喷水在线观看| 色精品久久人妻99蜜桃| 少妇猛男粗大的猛烈进出视频 | 此物有八面人人有两片| 99久久久亚洲精品蜜臀av| 久久精品国产鲁丝片午夜精品 | 国产精品福利在线免费观看| 欧美日韩乱码在线| 日本三级黄在线观看| 国产熟女欧美一区二区| 国产成人福利小说| 亚洲精华国产精华精| 国产成人a区在线观看| 成人高潮视频无遮挡免费网站| 国产乱人视频| 久久久久久国产a免费观看| 午夜福利视频1000在线观看| 老司机午夜福利在线观看视频| 久久人妻av系列| 天天一区二区日本电影三级| 老熟妇仑乱视频hdxx| 欧美bdsm另类| 国产精品免费一区二区三区在线| 免费看美女性在线毛片视频| 99久久精品国产国产毛片| 一卡2卡三卡四卡精品乱码亚洲| 日韩大尺度精品在线看网址| 特级一级黄色大片| 久久久久久九九精品二区国产| 精品乱码久久久久久99久播| 男插女下体视频免费在线播放| 成人二区视频| 啦啦啦观看免费观看视频高清| 一个人看视频在线观看www免费| 亚洲欧美日韩无卡精品| 亚洲av成人av| 夜夜爽天天搞| 精品无人区乱码1区二区| 精品国内亚洲2022精品成人| 黄色日韩在线| 成人av一区二区三区在线看| 九九热线精品视视频播放| 国产精华一区二区三区| 欧美极品一区二区三区四区| 欧美日韩国产亚洲二区| 国内精品美女久久久久久| 欧美日韩精品成人综合77777| 97人妻精品一区二区三区麻豆| 国产成人福利小说| 欧美性感艳星| xxxwww97欧美| 看片在线看免费视频| 国产老妇女一区| 日韩欧美一区二区三区在线观看| 性插视频无遮挡在线免费观看| 亚洲av电影不卡..在线观看| 亚洲精品一卡2卡三卡4卡5卡| 人人妻人人澡欧美一区二区| 色哟哟·www| 亚洲精品456在线播放app | 搡老岳熟女国产| 变态另类丝袜制服| 午夜精品在线福利| 中文亚洲av片在线观看爽| 国产极品精品免费视频能看的| 亚洲综合色惰| 可以在线观看的亚洲视频| 99热精品在线国产| 熟女电影av网| 国产成人福利小说| 久久久久九九精品影院| 中文字幕免费在线视频6| 国产免费男女视频| 国产高清视频在线播放一区| 亚洲熟妇中文字幕五十中出| 国产成人av教育| 亚洲天堂国产精品一区在线| 黄色配什么色好看| 国产熟女欧美一区二区| 亚洲精品日韩av片在线观看| 少妇裸体淫交视频免费看高清| 中文字幕免费在线视频6| 亚洲无线观看免费| 91久久精品电影网| 香蕉av资源在线| 精品久久国产蜜桃| 人妻少妇偷人精品九色| 乱码一卡2卡4卡精品| 色综合色国产| 午夜福利18| 俄罗斯特黄特色一大片| 一进一出抽搐动态| 免费一级毛片在线播放高清视频| 九九在线视频观看精品| ponron亚洲| 香蕉av资源在线| 97热精品久久久久久| 欧美另类亚洲清纯唯美| 又紧又爽又黄一区二区| 久久久久久久久久黄片| 国产美女午夜福利| 欧美日本视频| 少妇被粗大猛烈的视频| 日本免费a在线| 久久草成人影院| 亚洲avbb在线观看| 欧美区成人在线视频| 亚洲美女视频黄频| 亚洲欧美日韩卡通动漫| 欧美3d第一页| av专区在线播放| 成人av一区二区三区在线看| 搡女人真爽免费视频火全软件 | 色播亚洲综合网| 69人妻影院| 一夜夜www| 亚洲av熟女| 亚洲最大成人av| 99久久精品国产国产毛片| 欧美3d第一页| 日韩欧美国产在线观看| 亚洲不卡免费看| 啦啦啦韩国在线观看视频| 日本成人三级电影网站| 亚洲乱码一区二区免费版| 亚洲精华国产精华精| 日韩欧美国产在线观看| 午夜福利高清视频| 深爱激情五月婷婷| 亚洲av不卡在线观看| 国产毛片a区久久久久| 亚洲18禁久久av| 少妇裸体淫交视频免费看高清| aaaaa片日本免费| 精品久久久久久成人av| 久久久久久久久久成人| 国产真实伦视频高清在线观看 | 天堂动漫精品| 亚洲av熟女| 精品久久久久久久久久免费视频| 91久久精品电影网| 十八禁网站免费在线| 久久国内精品自在自线图片| 99精品在免费线老司机午夜| 大型黄色视频在线免费观看| a级毛片a级免费在线| 国产日本99.免费观看| 欧美+日韩+精品| 一a级毛片在线观看| 97热精品久久久久久| 成人美女网站在线观看视频| 日韩一本色道免费dvd| av在线天堂中文字幕| 中出人妻视频一区二区| 国产高清不卡午夜福利| 亚州av有码| 狂野欧美白嫩少妇大欣赏| 中文字幕av在线有码专区| 亚洲精品一区av在线观看| 噜噜噜噜噜久久久久久91| 亚洲一区高清亚洲精品| 日韩国内少妇激情av| 亚洲经典国产精华液单| 亚洲经典国产精华液单| 97人妻精品一区二区三区麻豆| 99九九线精品视频在线观看视频| x7x7x7水蜜桃| 亚洲欧美日韩东京热| 男女下面进入的视频免费午夜| 成人无遮挡网站| 九色国产91popny在线| 1000部很黄的大片| 嫩草影院精品99| 欧美中文日本在线观看视频| 黄色一级大片看看| 国产伦精品一区二区三区视频9| 免费搜索国产男女视频| 99热这里只有精品一区| 成人av一区二区三区在线看| 久久久久久久精品吃奶| 日本三级黄在线观看| 12—13女人毛片做爰片一| 亚洲狠狠婷婷综合久久图片| 精品欧美国产一区二区三| 国产av一区在线观看免费| 午夜福利在线观看免费完整高清在 | 免费在线观看影片大全网站| 久久久久久久午夜电影| 乱系列少妇在线播放| 99视频精品全部免费 在线| 亚洲熟妇熟女久久| 又爽又黄无遮挡网站| 在线观看免费视频日本深夜| 男女下面进入的视频免费午夜| 亚洲欧美激情综合另类| h日本视频在线播放| 九色成人免费人妻av| 亚洲成a人片在线一区二区| 日韩国内少妇激情av| 久久精品综合一区二区三区| 亚洲成人久久爱视频| 九九爱精品视频在线观看| 欧美3d第一页| 国产精品1区2区在线观看.| 久久久久久久久久黄片| 日日干狠狠操夜夜爽| 日本在线视频免费播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美日本亚洲视频在线播放| 国产成人影院久久av| 美女被艹到高潮喷水动态| 国产精品综合久久久久久久免费| 免费看美女性在线毛片视频| 99久久成人亚洲精品观看| 日本免费一区二区三区高清不卡| 最近中文字幕高清免费大全6 | 国内毛片毛片毛片毛片毛片| 亚洲国产精品久久男人天堂| 在线a可以看的网站| 亚洲av免费在线观看| 国产成人aa在线观看| 看黄色毛片网站| 日本色播在线视频| 国产黄色小视频在线观看| 18禁黄网站禁片午夜丰满| 少妇熟女aⅴ在线视频| 国内揄拍国产精品人妻在线| 高清毛片免费观看视频网站| a级一级毛片免费在线观看| 日本色播在线视频| 波多野结衣高清作品| 免费在线观看成人毛片| 欧美激情国产日韩精品一区| 久久热精品热| 欧美黑人巨大hd| 看片在线看免费视频| 国产真实伦视频高清在线观看 | 在线观看午夜福利视频| 天堂影院成人在线观看| 亚洲第一区二区三区不卡| 国国产精品蜜臀av免费| 日本三级黄在线观看| 亚洲国产精品合色在线| 97超视频在线观看视频| aaaaa片日本免费| 国产黄片美女视频| 国产男人的电影天堂91| 日韩 亚洲 欧美在线| 人人妻人人看人人澡| 最好的美女福利视频网| 欧美激情在线99| 中文在线观看免费www的网站| av中文乱码字幕在线| 性插视频无遮挡在线免费观看| 精品久久久噜噜| 成人三级黄色视频| 男人的好看免费观看在线视频| 亚洲性久久影院| 亚洲av美国av| 国产单亲对白刺激| 在线观看免费视频日本深夜| 欧美日韩中文字幕国产精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 国产乱人视频| 久久6这里有精品| 午夜精品一区二区三区免费看| 男人舔奶头视频| 欧美成人性av电影在线观看| 亚洲第一电影网av| 免费在线观看影片大全网站| 亚洲国产欧美人成| 久久99热6这里只有精品| 国产老妇女一区| 俄罗斯特黄特色一大片| 亚洲人成网站高清观看| 91久久精品电影网| 91久久精品国产一区二区成人| 99热精品在线国产| 久久亚洲精品不卡| 无人区码免费观看不卡| 伦理电影大哥的女人| 村上凉子中文字幕在线| 国产成人影院久久av| 国产成人一区二区在线| 麻豆久久精品国产亚洲av| 搡老妇女老女人老熟妇| av专区在线播放| 日韩欧美一区二区三区在线观看| 能在线免费观看的黄片| 日韩精品中文字幕看吧| a级毛片a级免费在线| 一本精品99久久精品77| 少妇人妻一区二区三区视频| av女优亚洲男人天堂| 又黄又爽又免费观看的视频| 琪琪午夜伦伦电影理论片6080| 国产大屁股一区二区在线视频| 日本与韩国留学比较| 麻豆国产av国片精品| 女同久久另类99精品国产91| 超碰av人人做人人爽久久| 免费在线观看日本一区| 国产蜜桃级精品一区二区三区| 啦啦啦啦在线视频资源| 亚洲第一区二区三区不卡| 我要搜黄色片| 啦啦啦观看免费观看视频高清| 少妇被粗大猛烈的视频| 亚洲最大成人手机在线| 亚洲性久久影院| 长腿黑丝高跟| 九九爱精品视频在线观看| 最好的美女福利视频网| 欧美一区二区国产精品久久精品| 男女那种视频在线观看| 国产一区二区在线观看日韩| 韩国av在线不卡| 日日摸夜夜添夜夜添av毛片 | 观看美女的网站| 国产大屁股一区二区在线视频| 欧美日韩精品成人综合77777| 久久久色成人|