Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple the...Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.展开更多
This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)metho...This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.展开更多
In this paper,we present a semi-Lagrangian(SL)method based on a non-polynomial function space for solving the Vlasov equation.We fnd that a non-polynomial function based scheme is suitable to the specifcs of the targe...In this paper,we present a semi-Lagrangian(SL)method based on a non-polynomial function space for solving the Vlasov equation.We fnd that a non-polynomial function based scheme is suitable to the specifcs of the target problems.To address issues that arise in phase space models of plasma problems,we develop a weighted essentially non-oscillatory(WENO)scheme using trigonometric polynomials.In particular,the non-polynomial WENO method is able to achieve improved accuracy near sharp gradients or discontinuities.Moreover,to obtain a high-order of accuracy in not only space but also time,it is proposed to apply a high-order splitting scheme in time.We aim to introduce the entire SL algorithm with high-order splitting in time and high-order WENO reconstruction in space to solve the Vlasov-Poisson system.Some numerical experiments are presented to demonstrate robustness of the proposed method in having a high-order of convergence and in capturing non-smooth solutions.A key observation is that the method can capture phase structure that require twice the resolution with a polynomial based method.In 6D,this would represent a signifcant savings.展开更多
A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a G...A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a Galerkin procedure is developed for the spatial discretization of the generalized nonlinear Schr6dinger (NLS) equa- tions, and a system of ordinary differential equations for the time dependent unknowns is obtained. Then, the classical fourth-order explicit Runge-Kutta method is used to solve this semi-discretization system. To justify the present method, several widely considered problems are solved as the test examples, and the results demonstrate that the proposed wavelet algorithm has much better accuracy and a faster convergence rate in space than many existing numerical methods.展开更多
This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink (soliton)...This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink (soliton) solutions and multiple-singular-kink (soliton) solutions are derived for the two equations.展开更多
This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been...This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chemically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic properties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicomponent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.展开更多
This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept...This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept of equivalent shear stiffness which can meet the requirement of the HSM algorithm.A study is done to theoretically validate the technique by the numerical analysis of two-storey shear building structure,in comparison of the proposed substructure pseudo-dynamic testing algorithm with the central difference method(CDM).Then,a full-scale SPDT model,the three-storey frame-supported reinforced concrete short-limb masonry shear wall structure,is designed and tested to simulate the seismic response of the corresponding six-storey structure and verify the proposed force control HSM technique.Meanwhile,the techniques of both stiffness correction and force control are suggested to control algorithmic error,control error and measurement error.The results indicate that the force control HSM can be used in the full-scale multi-degree-of-freedom(MDOF)substructure pseudo-dynamic testing before descent segment of structure restoring force properties.展开更多
The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass cons...The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.展开更多
An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. I...An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.展开更多
A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order pro...A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.展开更多
An effective numerical model for wave propagation over three-dimensional(3D)bathymetry was developed based on the High-Order Spectral(HOS)method and combined with a moving bottom boundary.Based on this model,tsunami w...An effective numerical model for wave propagation over three-dimensional(3D)bathymetry was developed based on the High-Order Spectral(HOS)method and combined with a moving bottom boundary.Based on this model,tsunami waves caused by various mechanisms were simulated and analyzed.Two-dimensional bed upthrust and the effect of the uplift velocity of the bathymetry on the wave profiles of tsunami waves were studied.Next,tsunami waves caused by 3D submarine slides were generated and the effects of the slide velocity,slide dimension and water depth on the tsunami waves were analyzed.Based on wavelet analysis,the properties of the tsunami wave propagation were investigated.The results show that the bottom movement can significantly affect the generation and propagation of tsunami waves and the studies could help understand the mechanisms of tsunamis caused by a moving bottom boundary.展开更多
The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi...The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi_harmonic operator and Laplace operator, applying the multiple reciprocity method(MRM), the MRM boundary integral equation for plate bending problem was constructed. It proves that the boundary integral equation derived from MRM is essentially identical to the conventional boundary integral equation. Hence the convergence analysis of MRM for plate bending problem can be obtained by the error estimation for the conventional boundary integral equation. In addition, this method can extend to the case of more series of the high_order fundamental solution sequences.展开更多
In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the...In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the time and space derivatives by compact finite-difference approximations. The system of resulting nonlinear finite difference equations are solved by preconditioned Krylov subspace methods. Numerical results are given to verify the behavior of high-order compact approximations in combination preconditioned methods for stability, convergence. Also, the accuracy and efficiency of the proposed scheme are considered.展开更多
An efficient high-order immersed interface method (IIM) is proposed to solve two-dimensional (2D) heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in ...An efficient high-order immersed interface method (IIM) is proposed to solve two-dimensional (2D) heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact (HOC) difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit (ADI) scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.展开更多
This paper considers practical, high-order methods for the iterative location of the roots of nonlinear equations, one at a time. Special attention is being paid to algorithms also applicable to multiple roots of init...This paper considers practical, high-order methods for the iterative location of the roots of nonlinear equations, one at a time. Special attention is being paid to algorithms also applicable to multiple roots of initially known and unknown multiplicity. Efficient methods are presented in this note for the evaluation of the multiplicity index of the root being sought. Also reviewed here are super-linear and super-cubic methods that converge contrarily or alternatingly, enabling us, not only to approach the root briskly and confidently but also to actually bound and bracket it as we progress.展开更多
We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-or...We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.展开更多
The Gerdjikov-Ivanov(GI)hierarchy is derived via recursion operator,in this article,we mainly investigate the third-order flow GI equation.In the framework of the Riemann-Hilbert method,the soliton matrices of the thi...The Gerdjikov-Ivanov(GI)hierarchy is derived via recursion operator,in this article,we mainly investigate the third-order flow GI equation.In the framework of the Riemann-Hilbert method,the soliton matrices of the third-order flow GI equation with simple zeros and elementary high-order zeros of Riemann-Hilbert problem are constructed through the standard dressing process.Taking advantage of this result,some properties and asymptotic analysis of single soliton solution and two soliton solution are discussed,and the simple elastic interaction of two soliton are proved.Compared with soliton solution of the classical second-order flow,we find that the higher-order dispersion term affects the propagation velocity,propagation direction and amplitude of the soliton.Finally,by means of a certain limit technique,the high-order soliton solution matrix for the third-order flow GI equation is derived.展开更多
A transmission grating coupled with an X-ray charge coupled device (CCD) is used to quantitatively measure the proportion of high-order harmonics of the soft-X-ray source of beam line 4B7B. The results show that the...A transmission grating coupled with an X-ray charge coupled device (CCD) is used to quantitatively measure the proportion of high-order harmonics of the soft-X-ray source of beam line 4B7B. The results show that the monochromatic X-ray has third-order and second-order harmonics. The proportion of second-order harmonic of 4B7B is less than 9.0% and the third- order harmonic is below 0.7% when no suppressing method is applied. When suppression methods are used, the proportion of second-order harmonic is less than 1.7% and the third-order harmonic is ignorable.展开更多
In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solution...In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solutions are only valid for small values of the independent variable. The DTM solutions diverge for some differential equations that extremely have nonlinear behaviors or have boundary-conditions at infinity. For this reason the governing boundary-layer equations are solved by the Multi-step Differential Transform Method (MDTM). The main advantage of this method is that it can be applied directly to nonlinear differential equations without requiring linearization, discretization, or perturbation. It is a semi analytical-numerical technique that formulizes Taylor series in a very different manner. By applying the MDTM the interval of convergence for the series solution is increased. The MDTM is treated as an algorithm in a sequence of intervals for finding accurate approximate solutions for systems of differential equations. It is predicted that the MDTM can be applied to a wide range of engineering applications.展开更多
This paper mainly studies the data characteristics of high order cumulants using digitally modulated signals, and constructs the identification feature parameters that can distinguish the signal modulation type by the...This paper mainly studies the data characteristics of high order cumulants using digitally modulated signals, and constructs the identification feature parameters that can distinguish the signal modulation type by the high-order cumulants data of the digital modulation signal. Set the identification signal modulation type determination threshold based on the value of the identification feature parameter. The identification feature parameter value of the signal modulation type is compared with the set determination threshold, to realize the recognition of digital modulation signal. This identification method is implemented based on MATLAB design, with a 2ASK (2-ary Amplitude Shift Keying) signal, 4ASK (4-ary Amplitude Shift Keying) signal, 2PSK (2-ary Phase Shift Keying) signal, 4PSK (4-ary Phase Shift Keying) signal, 2FSK (2-ary Frequency Shift Keying) signal, 4FSK (4-ary Frequency Shift Keying) signal. The second, fourth and sixth order cumulants of the six signals were analyzed. Calculate the selected identification feature parameter value and the determination threshold to identify the six signals. The six signals have made MATLAB identification simulation. Simulation results show that this method is feasible and has high recognition rate. Simulation results verify that such recognition methods maintain a high recognition rate under conditions with low signal-to-noise ratio. This identification method can be extended to more MASK (M-ary Amplitude Shift Keying), MPSK (M-ary Phase Shift Keying), MFSK (M-ary Frequency Shift Keying), MQAM (M-ary Quadrature Amplitude Modulation) signal identification.展开更多
基金supported by the NSFC Grant no.12271492the Natural Science Foundation of Henan Province of China Grant no.222300420550+1 种基金supported by the NSFC Grant no.12271498the National Key R&D Program of China Grant no.2022YFA1005202/2022YFA1005200.
文摘Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.
基金the National Natural Science Foundation of China(Grant Nos.71961022,11902163,12265020,and 12262024)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(Grant Nos.2019BS01011 and 2022MS01003)+5 种基金2022 Inner Mongolia Autonomous Region Grassland Talents Project-Young Innovative and Entrepreneurial Talents(Mingjing Du)2022 Talent Development Foundation of Inner Mongolia Autonomous Region of China(Ming-Jing Du)the Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region Program(Grant No.NJYT-20-B18)the Key Project of High-quality Economic Development Research Base of Yellow River Basin in 2022(Grant No.21HZD03)2022 Inner Mongolia Autonomous Region International Science and Technology Cooperation High-end Foreign Experts Introduction Project(Ge Kai)MOE(Ministry of Education in China)Humanities and Social Sciences Foundation(Grants No.20YJC860005).
文摘This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.
基金AFOSR and NSF for their support of this work under grants FA9550-19-1-0281 and FA9550-17-1-0394 and NSF grant DMS 191218。
文摘In this paper,we present a semi-Lagrangian(SL)method based on a non-polynomial function space for solving the Vlasov equation.We fnd that a non-polynomial function based scheme is suitable to the specifcs of the target problems.To address issues that arise in phase space models of plasma problems,we develop a weighted essentially non-oscillatory(WENO)scheme using trigonometric polynomials.In particular,the non-polynomial WENO method is able to achieve improved accuracy near sharp gradients or discontinuities.Moreover,to obtain a high-order of accuracy in not only space but also time,it is proposed to apply a high-order splitting scheme in time.We aim to introduce the entire SL algorithm with high-order splitting in time and high-order WENO reconstruction in space to solve the Vlasov-Poisson system.Some numerical experiments are presented to demonstrate robustness of the proposed method in having a high-order of convergence and in capturing non-smooth solutions.A key observation is that the method can capture phase structure that require twice the resolution with a polynomial based method.In 6D,this would represent a signifcant savings.
基金supported by the National Natural Science Foundation of China(Nos.11502103 and11421062)the Open Fund of State Key Laboratory of Structural Analysis for Industrial Equipment of China(No.GZ15115)
文摘A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a Galerkin procedure is developed for the spatial discretization of the generalized nonlinear Schr6dinger (NLS) equa- tions, and a system of ordinary differential equations for the time dependent unknowns is obtained. Then, the classical fourth-order explicit Runge-Kutta method is used to solve this semi-discretization system. To justify the present method, several widely considered problems are solved as the test examples, and the results demonstrate that the proposed wavelet algorithm has much better accuracy and a faster convergence rate in space than many existing numerical methods.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10871117 and 10571110)
文摘This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink (soliton) solutions and multiple-singular-kink (soliton) solutions are derived for the two equations.
基金supported by an Early Career Faculty grant from NASA's Space Technology Research Grants Programprovided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center
文摘This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chemically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic properties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicomponent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.
基金Sponsored by the National Natural Science Foundation of China(Grant No.50508012)
文摘This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept of equivalent shear stiffness which can meet the requirement of the HSM algorithm.A study is done to theoretically validate the technique by the numerical analysis of two-storey shear building structure,in comparison of the proposed substructure pseudo-dynamic testing algorithm with the central difference method(CDM).Then,a full-scale SPDT model,the three-storey frame-supported reinforced concrete short-limb masonry shear wall structure,is designed and tested to simulate the seismic response of the corresponding six-storey structure and verify the proposed force control HSM technique.Meanwhile,the techniques of both stiffness correction and force control are suggested to control algorithmic error,control error and measurement error.The results indicate that the force control HSM can be used in the full-scale multi-degree-of-freedom(MDOF)substructure pseudo-dynamic testing before descent segment of structure restoring force properties.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11862003 and 81860635)the Key Project of the Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2017GXNSFDA198038)+3 种基金the Project of Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2018GXNSFAA281302)the Project for Promotion of Young and Middle-aged Teachers’Basic Scientific Research Ability in Guangxi Universities,China(Grant No.2019KY0084)the“Bagui Scholar”Teams for Innovation and Research Project of Guangxi Zhuang Autonomous Region,Chinathe Graduate Innovation Program of Guangxi Normal University,China(Grant No.JXYJSKT-2019-007)。
文摘The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.
基金co-supported by the National High Technology Research and Development Program of China(No.2015AA015303)the National Natural Science Foundation of China(No.11272152)+1 种基金the Aeronautical Science Foundation of China(No.20152752033)the Open Project of Key Laboratory of Aerodynamic Noise Control
文摘An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.
文摘A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.
基金This study was financially supported by the National Natural Science Foundation of China(Grant Nos.51739010 and 51879037).
文摘An effective numerical model for wave propagation over three-dimensional(3D)bathymetry was developed based on the High-Order Spectral(HOS)method and combined with a moving bottom boundary.Based on this model,tsunami waves caused by various mechanisms were simulated and analyzed.Two-dimensional bed upthrust and the effect of the uplift velocity of the bathymetry on the wave profiles of tsunami waves were studied.Next,tsunami waves caused by 3D submarine slides were generated and the effects of the slide velocity,slide dimension and water depth on the tsunami waves were analyzed.Based on wavelet analysis,the properties of the tsunami wave propagation were investigated.The results show that the bottom movement can significantly affect the generation and propagation of tsunami waves and the studies could help understand the mechanisms of tsunamis caused by a moving bottom boundary.
文摘The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi_harmonic operator and Laplace operator, applying the multiple reciprocity method(MRM), the MRM boundary integral equation for plate bending problem was constructed. It proves that the boundary integral equation derived from MRM is essentially identical to the conventional boundary integral equation. Hence the convergence analysis of MRM for plate bending problem can be obtained by the error estimation for the conventional boundary integral equation. In addition, this method can extend to the case of more series of the high_order fundamental solution sequences.
文摘In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the time and space derivatives by compact finite-difference approximations. The system of resulting nonlinear finite difference equations are solved by preconditioned Krylov subspace methods. Numerical results are given to verify the behavior of high-order compact approximations in combination preconditioned methods for stability, convergence. Also, the accuracy and efficiency of the proposed scheme are considered.
基金supported by the National Natural Science Foundation of China(No.51174236)the National Basic Research Program of China(973 Program)(No.2011CB606306)the Opening Project of State Key Laboratory of Porous Metal Materials(No.PMM-SKL-4-2012)
文摘An efficient high-order immersed interface method (IIM) is proposed to solve two-dimensional (2D) heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact (HOC) difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit (ADI) scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.
文摘This paper considers practical, high-order methods for the iterative location of the roots of nonlinear equations, one at a time. Special attention is being paid to algorithms also applicable to multiple roots of initially known and unknown multiplicity. Efficient methods are presented in this note for the evaluation of the multiplicity index of the root being sought. Also reviewed here are super-linear and super-cubic methods that converge contrarily or alternatingly, enabling us, not only to approach the root briskly and confidently but also to actually bound and bracket it as we progress.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11571366 and 11501570the Open Foundation of State Key Laboratory of High Performance Computing of China+1 种基金the Research Fund of National University of Defense Technology under Grant No JC15-02-02the Fund from HPCL
文摘We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.
基金supported by the National Natural Science Foundation of China(No.12175069 and No.12235007)Science and Technology Commission of Shanghai Municipality(No.21JC1402500 and No.22DZ2229014)Natural Science Foundation of Shanghai,China(No.23ZR1418100).
文摘The Gerdjikov-Ivanov(GI)hierarchy is derived via recursion operator,in this article,we mainly investigate the third-order flow GI equation.In the framework of the Riemann-Hilbert method,the soliton matrices of the third-order flow GI equation with simple zeros and elementary high-order zeros of Riemann-Hilbert problem are constructed through the standard dressing process.Taking advantage of this result,some properties and asymptotic analysis of single soliton solution and two soliton solution are discussed,and the simple elastic interaction of two soliton are proved.Compared with soliton solution of the classical second-order flow,we find that the higher-order dispersion term affects the propagation velocity,propagation direction and amplitude of the soliton.Finally,by means of a certain limit technique,the high-order soliton solution matrix for the third-order flow GI equation is derived.
文摘A transmission grating coupled with an X-ray charge coupled device (CCD) is used to quantitatively measure the proportion of high-order harmonics of the soft-X-ray source of beam line 4B7B. The results show that the monochromatic X-ray has third-order and second-order harmonics. The proportion of second-order harmonic of 4B7B is less than 9.0% and the third- order harmonic is below 0.7% when no suppressing method is applied. When suppression methods are used, the proportion of second-order harmonic is less than 1.7% and the third-order harmonic is ignorable.
文摘In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solutions are only valid for small values of the independent variable. The DTM solutions diverge for some differential equations that extremely have nonlinear behaviors or have boundary-conditions at infinity. For this reason the governing boundary-layer equations are solved by the Multi-step Differential Transform Method (MDTM). The main advantage of this method is that it can be applied directly to nonlinear differential equations without requiring linearization, discretization, or perturbation. It is a semi analytical-numerical technique that formulizes Taylor series in a very different manner. By applying the MDTM the interval of convergence for the series solution is increased. The MDTM is treated as an algorithm in a sequence of intervals for finding accurate approximate solutions for systems of differential equations. It is predicted that the MDTM can be applied to a wide range of engineering applications.
文摘This paper mainly studies the data characteristics of high order cumulants using digitally modulated signals, and constructs the identification feature parameters that can distinguish the signal modulation type by the high-order cumulants data of the digital modulation signal. Set the identification signal modulation type determination threshold based on the value of the identification feature parameter. The identification feature parameter value of the signal modulation type is compared with the set determination threshold, to realize the recognition of digital modulation signal. This identification method is implemented based on MATLAB design, with a 2ASK (2-ary Amplitude Shift Keying) signal, 4ASK (4-ary Amplitude Shift Keying) signal, 2PSK (2-ary Phase Shift Keying) signal, 4PSK (4-ary Phase Shift Keying) signal, 2FSK (2-ary Frequency Shift Keying) signal, 4FSK (4-ary Frequency Shift Keying) signal. The second, fourth and sixth order cumulants of the six signals were analyzed. Calculate the selected identification feature parameter value and the determination threshold to identify the six signals. The six signals have made MATLAB identification simulation. Simulation results show that this method is feasible and has high recognition rate. Simulation results verify that such recognition methods maintain a high recognition rate under conditions with low signal-to-noise ratio. This identification method can be extended to more MASK (M-ary Amplitude Shift Keying), MPSK (M-ary Phase Shift Keying), MFSK (M-ary Frequency Shift Keying), MQAM (M-ary Quadrature Amplitude Modulation) signal identification.