This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either...This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.展开更多
Boussinesq type equations have been widely studied to model the surface water wave.In this paper,we consider the abcd Boussinesq system which is a family of Boussinesq type equations including many well-known models s...Boussinesq type equations have been widely studied to model the surface water wave.In this paper,we consider the abcd Boussinesq system which is a family of Boussinesq type equations including many well-known models such as the classical Boussinesq system,the BBM-BBM system,the Bona-Smith system,etc.We propose local discontinuous Galerkin(LDG)methods,with carefully chosen numerical fluxes,to numerically solve this abcd Boussinesq system.The main focus of this paper is to rigorously establish a priori error estimate of the proposed LDG methods for a wide range of the parameters a,b,c,d.Numerical experiments are shown to test the convergence rates,and to demonstrate that the proposed methods can simulate the head-on collision of traveling wave and finite time blow-up behavior well.展开更多
In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materi...In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen- Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to dem-onstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use rela-tively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.展开更多
In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We the...In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.展开更多
The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes ...The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes fan be assumed equal to arbitrary points,where the integrand function f is known; iii) the number of the requested evaluations of f at the nodes is low,iv) a satisfactory convergence theory can be proved.展开更多
In this paper,we develop novel local discontinuous Galerkin(LDG)methods for fractional diffusion equations with non-smooth solutions.We consider such problems,for which the solutions are not smooth at boundary,and the...In this paper,we develop novel local discontinuous Galerkin(LDG)methods for fractional diffusion equations with non-smooth solutions.We consider such problems,for which the solutions are not smooth at boundary,and therefore the traditional LDG methods with piecewise polynomial solutions suffer accuracy degeneracy.The novel LDG methods utilize a solution information enriched basis,simulate the problem on a paired special mesh,and achieve optimal order of accuracy.We analyze the L2 stability and optimal error estimate in L2-norm.Finally,numerical examples are presented for validating the theoretical conclusions.展开更多
The localization characteristics of the in-plane elastic waves in locally resonant aperiodic phononic crystals are examined in this study. In particular, the phononic crystals generated according to the Thue-Morse, Ru...The localization characteristics of the in-plane elastic waves in locally resonant aperiodic phononic crystals are examined in this study. In particular, the phononic crystals generated according to the Thue-Morse, Rudin-Shapiro and Period-Doubling sequences are theoretically investigated by using the extended transfer matrix method. For comparison, the binary and ternary locally resonant systems are considered, and their band structures are characterized by using the localization factors. Moreover, the influences of structural arrangement, material combination, incidence angle, number of components, length ratio, and random disorder on the band structures are also discussed. Some novel and interesting phenomena are observed and discussed.展开更多
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 study examines the multicriteria scheduling problem on a single machine to minimize three criteria: the maximum cost function, denoted by maximum late work (V<sub>max</sub>), maximum tardy job, denote...This study examines the multicriteria scheduling problem on a single machine to minimize three criteria: the maximum cost function, denoted by maximum late work (V<sub>max</sub>), maximum tardy job, denoted by (T<sub>max</sub>), and maximum earliness (E<sub>max</sub>). We propose several algorithms based on types of objectives function to be optimized when dealing with simultaneous minimization problems with and without weight and hierarchical minimization problems. The proposed Algorithm (3) is to find the set of efficient solutions for 1//F (V<sub>max</sub>, T<sub>max</sub>, E<sub>max</sub>) and 1//(V<sub>max</sub> + T<sub>max</sub> + E<sub>max</sub>). The Local Search Heuristic Methods (Descent Method (DM), Simulated Annealing (SA), Genetic Algorithm (GA), and the Tree Type Heuristics Method (TTHM) are applied to solve all suggested problems. Finally, the experimental results of Algorithm (3) are compared with the results of the Branch and Bound (BAB) method for optimal and Pareto optimal solutions for smaller instance sizes and compared to the Local Search Heuristic Methods for large instance sizes. These results ensure the efficiency of Algorithm (3) in a reasonable time.展开更多
This paper presents an improved approach based on the equivalent-weights particle filter(EWPF)that uses the proposal density to effectively improve the traditional particle filter.The proposed approach uses historical...This paper presents an improved approach based on the equivalent-weights particle filter(EWPF)that uses the proposal density to effectively improve the traditional particle filter.The proposed approach uses historical data to calculate statistical observations instead of the future observations used in the EWPF’s proposal density and draws on the localization scheme used in the localized PF(LPF)to construct the localized EWPF.The new approach is called the statistical observation localized EWPF(LEWPF-Sobs);it uses statistical observations that are better adapted to the requirements of real-time assimilation and the localization function is used to calculate weights to reduce the effect of missing observations on the weights.This approach not only retains the advantages of the EWPF,but also improves the assimilation quality when using sparse observations.Numerical experiments performed with the Lorenz 96 model show that the statistical observation EWPF is better than the EWPF and EAKF when the model uses standard distribution observations.Comparisons of the statistical observation localized EWPF and LPF reveal the advantages of the new method,with fewer particles giving better results.In particular,the new improved filter performs better than the traditional algorithms when the observation network contains densely spaced measurements associated with model state nonlinearities.展开更多
To study the association of oxytocin (OT)'s distribution in hypothalamatic,pituitary and ovary,and understand how the OT secrete releasing in hypothalamus,pituitary and ovaries,the paraffin section immunohistochem...To study the association of oxytocin (OT)'s distribution in hypothalamatic,pituitary and ovary,and understand how the OT secrete releasing in hypothalamus,pituitary and ovaries,the paraffin section immunohistochemistry SuperPicTureTM two step method was used to detect the distribution of OT in hypothalamatic-pituitary-ovary axis of five femal Guangxi local buffalo. The test results could provide morphology according to study the OT's synthesis and mechanism of action,and could play reference and directions part in breeding Guangxi local buffalo. The test results display:oxytocin immuno reactive (OT-IR) neuronsw eremainly distributed arcuate nucleus,supraoptic nucleus and paraventricular nucleus,and OT-IR neurons was also found in ventromedial nucleus,ventrolateralis nucleus,suprachiasmaticus nucleus,dorsomedial nucleus,mamillary body,anterior hypothalamic nucleus and so on. The OT immunoactive production was found in pituitary and few OT-IR nerve fibers extended to post pituitary from hypophyseal stalk and medium eminence. In ovaries,OT immunoactive productions were only distributed in germinal epithelium cells,granulosa cells and lutein cells. The OT was first discovered in singulorum link of hypothalamatic-pituitary-ovary axis of Guangxi local buffalo. The OT immunoactive neurons were first discovered in every main nucleus of Guangxi local buffalo hypothalamus,especially distributed in arcuate nucleus,supraoptic nucleus and paraventricular nucleus.展开更多
Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficul...Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.展开更多
A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are ...A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.展开更多
In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. ...In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O（ε1/2 ＋ h1/2）hk） convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.展开更多
In this article,a procedure for estimating the coefficient functions on the functional-coefficient regression models with different smoothing variables in different coefficient functions is defined.First step,by the l...In this article,a procedure for estimating the coefficient functions on the functional-coefficient regression models with different smoothing variables in different coefficient functions is defined.First step,by the local linear technique and the averaged method,the initial estimates of the coefficient functions are given.Second step,based on the initial estimates,the efficient estimates of the coefficient functions are proposed by a one-step back-fitting procedure.The efficient estimators share the same asymptotic normalities as the local linear estimators for the functional-coefficient models with a single smoothing variable in different functions.Two simulated examples show that the procedure is effective.展开更多
In this study, a multivariate local quadratic polynomial regression(MLQPR) method is proposed to design a model for the sludge volume index(SVI). In MLQPR, a quadratic polynomial regression function is established to ...In this study, a multivariate local quadratic polynomial regression(MLQPR) method is proposed to design a model for the sludge volume index(SVI). In MLQPR, a quadratic polynomial regression function is established to describe the relationship between SVI and the relative variables, and the important terms of the quadratic polynomial regression function are determined by the significant test of the corresponding coefficients. Moreover, a local estimation method is introduced to adjust the weights of the quadratic polynomial regression function to improve the model accuracy. Finally, the proposed method is applied to predict the SVI values in a real wastewater treatment process(WWTP). The experimental results demonstrate that the proposed MLQPR method has faster testing speed and more accurate results than some existing methods.展开更多
Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).How...Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).However,the deformation field obtained by GSRM could not reflect the real deformation of a slope when the slope became unstable.For most slopes,failure occurs once the strength of some regional soil is sufficiently weakened; thus,the local strength reduction method(LSRM)was proposed to analyze slope stability.In contrast with GSRM,LSRM only reduces the strength of local soil,while the strength of other soil remains unchanged.Therefore,deformation by LSRM is more reasonable than that by GSRM.In addition,the accuracy of the slope's deformation depends on the constitutive model to a large degree,and the variable-modulus elasto-plastic model was thus adopted.This constitutive model was an improvement of the Duncan–Chang model,which modified soil's deformation modulus according to stress level,and it thus better reflected the plastic feature of soil.Most importantly,the parameters of the variable-modulus elasto-plastic model could be determined through in-situ tests,and parameters determination by plate loading test and pressuremeter test were introduced.Therefore,it is easy to put this model into practice.Finally,LSRM and the variable-modulus elasto-plastic model were used to analyze Egongdai ancient landslide.Safety factor,deformation field,and optimal reinforcement measures for Egongdai ancient landslide were obtained based on the proposed method.展开更多
Fatigue fracture is one of the most common failure modes of engineering compo-nents,and the combined action of geometrie discontinuity and multiaxial loading is more likely to cause severe fatigue damage of components...Fatigue fracture is one of the most common failure modes of engineering compo-nents,and the combined action of geometrie discontinuity and multiaxial loading is more likely to cause severe fatigue damage of components.This work focuses on the fatigue behavior of U-notched Q345 steel specimens with differen t notch sizes under proportional cyclic tension-torsion.Firstly,based on the concept of strain energy,the calculation method of critical plane is given and the equivalent stress of the specified path on the critical plane is extracted to char-acterize the equivalent stress distribution state and the stress gradient effect.Then,based on the high stress volume method and theory of critical distance,a simple method for determining the critical distance is given considering the contribution of stress at the dangerous point and the critical point.In addition,based on the idea of stress-distance normalization,a new stress gradient impact factor is defined and a new method for predicting the multiaxial fatigue life of notched specimens is given.The prediction results of the proposed model,the local stress-strain method and the point method of theory of critical distance are compared with the experimental results.The comparisons show that the prediction results of the proposed model are closer to experimentai life,and the calculation accuracy is higher.展开更多
The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many fact...The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.展开更多
We study a nonlinear periodic problem driven by the p（t）-Laplacian and having a nonsmooth potential （hemivariational inequalities）. Using a variational method based on nonsmooth critical point theory for locally L...We study a nonlinear periodic problem driven by the p（t）-Laplacian and having a nonsmooth potential （hemivariational inequalities）. Using a variational method based on nonsmooth critical point theory for locally Lipschitz functions, we first prove the existence of at least two nontrivial solutions under the generalized subquadratic and then establish the existence of at least one nontrivial solution under the generalized superquadratic.展开更多
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.
基金The work of J.Sun and Y.Xing is partially sponsored by NSF grant DMS-1753581.
文摘Boussinesq type equations have been widely studied to model the surface water wave.In this paper,we consider the abcd Boussinesq system which is a family of Boussinesq type equations including many well-known models such as the classical Boussinesq system,the BBM-BBM system,the Bona-Smith system,etc.We propose local discontinuous Galerkin(LDG)methods,with carefully chosen numerical fluxes,to numerically solve this abcd Boussinesq system.The main focus of this paper is to rigorously establish a priori error estimate of the proposed LDG methods for a wide range of the parameters a,b,c,d.Numerical experiments are shown to test the convergence rates,and to demonstrate that the proposed methods can simulate the head-on collision of traveling wave and finite time blow-up behavior well.
基金Jie Du is supported by the National Natural Science Foundation of China under Grant Number NSFC 11801302Tsinghua University Initiative Scientific Research Program+1 种基金Eric Chung is supported by Hong Kong RGC General Research Fund(Projects 14304217 and 14302018)The third author is supported by the NSF grant DMS-1818467.
文摘In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen- Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to dem-onstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use rela-tively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.
基金The project supported by the China NKBRSF(2001CB409604)
文摘In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.
基金Work sponsored by"Ministero dell' University"CNR of Italy
文摘The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes fan be assumed equal to arbitrary points,where the integrand function f is known; iii) the number of the requested evaluations of f at the nodes is low,iv) a satisfactory convergence theory can be proved.
文摘In this paper,we develop novel local discontinuous Galerkin(LDG)methods for fractional diffusion equations with non-smooth solutions.We consider such problems,for which the solutions are not smooth at boundary,and therefore the traditional LDG methods with piecewise polynomial solutions suffer accuracy degeneracy.The novel LDG methods utilize a solution information enriched basis,simulate the problem on a paired special mesh,and achieve optimal order of accuracy.We analyze the L2 stability and optimal error estimate in L2-norm.Finally,numerical examples are presented for validating the theoretical conclusions.
基金the financial support from the National Natural Science Foundation of China (No. 11002026, 11372039)Beijing Natural Science Foundation (No. 3133039)the Scientific Research Foundation for the Returned Overseas Chinese Scholars (No. 20121832001)
文摘The localization characteristics of the in-plane elastic waves in locally resonant aperiodic phononic crystals are examined in this study. In particular, the phononic crystals generated according to the Thue-Morse, Rudin-Shapiro and Period-Doubling sequences are theoretically investigated by using the extended transfer matrix method. For comparison, the binary and ternary locally resonant systems are considered, and their band structures are characterized by using the localization factors. Moreover, the influences of structural arrangement, material combination, incidence angle, number of components, length ratio, and random disorder on the band structures are also discussed. Some novel and interesting phenomena are observed and discussed.
基金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.
文摘This study examines the multicriteria scheduling problem on a single machine to minimize three criteria: the maximum cost function, denoted by maximum late work (V<sub>max</sub>), maximum tardy job, denoted by (T<sub>max</sub>), and maximum earliness (E<sub>max</sub>). We propose several algorithms based on types of objectives function to be optimized when dealing with simultaneous minimization problems with and without weight and hierarchical minimization problems. The proposed Algorithm (3) is to find the set of efficient solutions for 1//F (V<sub>max</sub>, T<sub>max</sub>, E<sub>max</sub>) and 1//(V<sub>max</sub> + T<sub>max</sub> + E<sub>max</sub>). The Local Search Heuristic Methods (Descent Method (DM), Simulated Annealing (SA), Genetic Algorithm (GA), and the Tree Type Heuristics Method (TTHM) are applied to solve all suggested problems. Finally, the experimental results of Algorithm (3) are compared with the results of the Branch and Bound (BAB) method for optimal and Pareto optimal solutions for smaller instance sizes and compared to the Local Search Heuristic Methods for large instance sizes. These results ensure the efficiency of Algorithm (3) in a reasonable time.
基金The National Basic Research Program of China under contract Nos 2017YFC1404100,2017YFC1404103 and 2017YFC1404104the National Natural Science Foundation of China under contract No.41676088。
文摘This paper presents an improved approach based on the equivalent-weights particle filter(EWPF)that uses the proposal density to effectively improve the traditional particle filter.The proposed approach uses historical data to calculate statistical observations instead of the future observations used in the EWPF’s proposal density and draws on the localization scheme used in the localized PF(LPF)to construct the localized EWPF.The new approach is called the statistical observation localized EWPF(LEWPF-Sobs);it uses statistical observations that are better adapted to the requirements of real-time assimilation and the localization function is used to calculate weights to reduce the effect of missing observations on the weights.This approach not only retains the advantages of the EWPF,but also improves the assimilation quality when using sparse observations.Numerical experiments performed with the Lorenz 96 model show that the statistical observation EWPF is better than the EWPF and EAKF when the model uses standard distribution observations.Comparisons of the statistical observation localized EWPF and LPF reveal the advantages of the new method,with fewer particles giving better results.In particular,the new improved filter performs better than the traditional algorithms when the observation network contains densely spaced measurements associated with model state nonlinearities.
基金Supported by Guangxi Scientific Fund Project (Guikezi0991042, Guikezi 0640015 and Guikezi 0832043)Guangxi Area Education Department Educational and Scientific Layout Project (C, 2006C3)+1 种基金Guangxi Education Department Scientific Research Fund (200709LX075)Guangxi Large Apparatus Collaborated Sharing Net~~
文摘To study the association of oxytocin (OT)'s distribution in hypothalamatic,pituitary and ovary,and understand how the OT secrete releasing in hypothalamus,pituitary and ovaries,the paraffin section immunohistochemistry SuperPicTureTM two step method was used to detect the distribution of OT in hypothalamatic-pituitary-ovary axis of five femal Guangxi local buffalo. The test results could provide morphology according to study the OT's synthesis and mechanism of action,and could play reference and directions part in breeding Guangxi local buffalo. The test results display:oxytocin immuno reactive (OT-IR) neuronsw eremainly distributed arcuate nucleus,supraoptic nucleus and paraventricular nucleus,and OT-IR neurons was also found in ventromedial nucleus,ventrolateralis nucleus,suprachiasmaticus nucleus,dorsomedial nucleus,mamillary body,anterior hypothalamic nucleus and so on. The OT immunoactive production was found in pituitary and few OT-IR nerve fibers extended to post pituitary from hypophyseal stalk and medium eminence. In ovaries,OT immunoactive productions were only distributed in germinal epithelium cells,granulosa cells and lutein cells. The OT was first discovered in singulorum link of hypothalamatic-pituitary-ovary axis of Guangxi local buffalo. The OT immunoactive neurons were first discovered in every main nucleus of Guangxi local buffalo hypothalamus,especially distributed in arcuate nucleus,supraoptic nucleus and paraventricular nucleus.
基金Project supported by the National 973 Program (No.2004CB719402), the National Natural Science Foundation of China (No. 10372030)the Open Research Projects supported by the Project Fund of the Hubei Province Key Lab of Mechanical Transmission & Manufacturing Engineering Wuhan University of Science & Technology (No.2003A16).
文摘Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.
基金supported by the National Natural Science Foundation of China(11390363 and 11172041)Beijing Higher Education Young Elite Teacher Project(YETP1190)
文摘A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.
基金Supported by National Natural Science Foundation of China (10571046, 10571053, and 10871066)Program for New Century Excellent Talents in University (NCET-06-0712)+2 种基金Key Laboratory of Computational and Stochastic Mathematics and Its Applications, Universities of Hunan Province, Hunan Normal Universitythe Project of Scientific Research Fund of Hunan Provincial Education Department (09K025)the Key Scientific Research Topic of Jiaxing University (70110X05BL)
文摘In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O（ε1/2 ＋ h1/2）hk） convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.
文摘In this article,a procedure for estimating the coefficient functions on the functional-coefficient regression models with different smoothing variables in different coefficient functions is defined.First step,by the local linear technique and the averaged method,the initial estimates of the coefficient functions are given.Second step,based on the initial estimates,the efficient estimates of the coefficient functions are proposed by a one-step back-fitting procedure.The efficient estimators share the same asymptotic normalities as the local linear estimators for the functional-coefficient models with a single smoothing variable in different functions.Two simulated examples show that the procedure is effective.
文摘In this study, a multivariate local quadratic polynomial regression(MLQPR) method is proposed to design a model for the sludge volume index(SVI). In MLQPR, a quadratic polynomial regression function is established to describe the relationship between SVI and the relative variables, and the important terms of the quadratic polynomial regression function are determined by the significant test of the corresponding coefficients. Moreover, a local estimation method is introduced to adjust the weights of the quadratic polynomial regression function to improve the model accuracy. Finally, the proposed method is applied to predict the SVI values in a real wastewater treatment process(WWTP). The experimental results demonstrate that the proposed MLQPR method has faster testing speed and more accurate results than some existing methods.
基金Project([2005]205)supported by the Science and Technology Planning Project of Water Resources Department of Guangdong Province,ChinaProject(2012-7)supported by Guangdong Bureau of Highway Administration,ChinaProject(2012210020203)supported by the Fundamental Research Funds for the Central Universities,China
文摘Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).However,the deformation field obtained by GSRM could not reflect the real deformation of a slope when the slope became unstable.For most slopes,failure occurs once the strength of some regional soil is sufficiently weakened; thus,the local strength reduction method(LSRM)was proposed to analyze slope stability.In contrast with GSRM,LSRM only reduces the strength of local soil,while the strength of other soil remains unchanged.Therefore,deformation by LSRM is more reasonable than that by GSRM.In addition,the accuracy of the slope's deformation depends on the constitutive model to a large degree,and the variable-modulus elasto-plastic model was thus adopted.This constitutive model was an improvement of the Duncan–Chang model,which modified soil's deformation modulus according to stress level,and it thus better reflected the plastic feature of soil.Most importantly,the parameters of the variable-modulus elasto-plastic model could be determined through in-situ tests,and parameters determination by plate loading test and pressuremeter test were introduced.Therefore,it is easy to put this model into practice.Finally,LSRM and the variable-modulus elasto-plastic model were used to analyze Egongdai ancient landslide.Safety factor,deformation field,and optimal reinforcement measures for Egongdai ancient landslide were obtained based on the proposed method.
基金This research was supported by the National Natural Science Foundation of China(Grant No.51605212)the Natural Science Foundation of Gansu Province(Grant No.20JR10RA161)the Project of Hongliu Excellent Youth Program of Lanzhou University of Technology(Grant No.2020062001).
文摘Fatigue fracture is one of the most common failure modes of engineering compo-nents,and the combined action of geometrie discontinuity and multiaxial loading is more likely to cause severe fatigue damage of components.This work focuses on the fatigue behavior of U-notched Q345 steel specimens with differen t notch sizes under proportional cyclic tension-torsion.Firstly,based on the concept of strain energy,the calculation method of critical plane is given and the equivalent stress of the specified path on the critical plane is extracted to char-acterize the equivalent stress distribution state and the stress gradient effect.Then,based on the high stress volume method and theory of critical distance,a simple method for determining the critical distance is given considering the contribution of stress at the dangerous point and the critical point.In addition,based on the idea of stress-distance normalization,a new stress gradient impact factor is defined and a new method for predicting the multiaxial fatigue life of notched specimens is given.The prediction results of the proposed model,the local stress-strain method and the point method of theory of critical distance are compared with the experimental results.The comparisons show that the prediction results of the proposed model are closer to experimentai life,and the calculation accuracy is higher.
基金the Scientific Foundation of National Outstanding Youth of China(No.50225520)the Science Foundation of Shandong University of Technology of China(No.2006KJM33).
文摘The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.
基金supported by the National Science Foundation of China (11001063, 10971043)the Fundamental Research Funds for the Central Universities (HEUCF 20111134)+2 种基金China Postdoctoral Science Foundation Funded Project (20110491032)Heilongjiang Provincial Science Foundation for Distinguished Young Scholars (JC200810)Program of Excellent Team in Harbin Institute of Technology and the Natural Science Foundation of Heilongjiang Province (A200803)
文摘We study a nonlinear periodic problem driven by the p（t）-Laplacian and having a nonsmooth potential （hemivariational inequalities）. Using a variational method based on nonsmooth critical point theory for locally Lipschitz functions, we first prove the existence of at least two nontrivial solutions under the generalized subquadratic and then establish the existence of at least one nontrivial solution under the generalized superquadratic.
