An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into ...An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into an equivalent system,and the k-order backward differentiation formula(BDF k)and central difference formula are used to discretize the temporal and spatial derivatives,respectively.Different from the traditional discrete method that adopts full implicit or full explicit for the nonlinear integral terms,the proposed scheme is based on the SAV idea and can be treated semi-implicitly,taking into account both accuracy and effectiveness.Numerical results are presented to demonstrate the high-order convergence(up to fourth-order)of the developed schemes and it is computationally efficient in long-time computations.展开更多
This paper develops a generalized scalar auxiliary variable(SAV)method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent ...This paper develops a generalized scalar auxiliary variable(SAV)method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent Ginzburg-Landau equations.In this method,the system is decoupled and linearized to avoid solving the non-linear equation at each step.The theoretical analysis proves that the generalized SAV method can preserve the maximum bound principle and energy stability,and this is confirmed by the numerical result,and also shows that the numerical algorithm is stable.展开更多
This manuscript deals with new class of almost unbiased ratio cum product estimators for the estimation of population mean of the study variable by using the known values auxiliary variable. The bias and mean squared ...This manuscript deals with new class of almost unbiased ratio cum product estimators for the estimation of population mean of the study variable by using the known values auxiliary variable. The bias and mean squared error of proposed estimators are obtained. An empirical study is carried out to assess the efficiency of proposed estimators over the existing estimators with the help of some known natural populations and it shows that the proposed estimators are almost unbiased and it perform better than the existing estimators.展开更多
In this paper, we have proposed estimators of finite population mean using generalized Ratio- cum-product estimator for two-Phase sampling using multi-auxiliary variables under full, partial and no information cases a...In this paper, we have proposed estimators of finite population mean using generalized Ratio- cum-product estimator for two-Phase sampling using multi-auxiliary variables under full, partial and no information cases and investigated their finite sample properties. An empirical study is given to compare the performance of the proposed estimators with the existing estimators that utilize auxiliary variable(s) for finite population mean. It has been found that the generalized Ra-tio-cum-product estimator in full information case using multiple auxiliary variables is more efficient than mean per unit, ratio and product estimator using one auxiliary variable, ratio and product estimator using multiple auxiliary variable and ratio-cum-product estimators in both partial and no information case in two phase sampling. A generalized Ratio-cum-product estimator in partial information case is more efficient than Generalized Ratio-cum-product estimator in No information case.展开更多
In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is int...In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is introduced for the physical system;the mass conservation and energy dissipation is proved for the CHHS system.Subsequently,a fully discrete SAV finite element scheme is proposed,with the mass conservation and energy dissipation laws established at a theoretical level.In addition,the convergence analysis and error estimate is provided for the proposed SAV numerical scheme.展开更多
A second order accurate(in time)numerical scheme is analyzed for the slope-selection(SS)equation of the epitaxial thin film growth model,with Fourier pseudo-spectral discretization in space.To make the numerical schem...A second order accurate(in time)numerical scheme is analyzed for the slope-selection(SS)equation of the epitaxial thin film growth model,with Fourier pseudo-spectral discretization in space.To make the numerical scheme linear while preserving the nonlinear energy stability,we make use of the scalar auxiliary variable(SAV)approach,in which a modified Crank-Nicolson is applied for the surface diffusion part.The energy stability could be derived a modified form,in comparison with the standard Crank-Nicolson approximation to the surface diffusion term.Such an energy stability leads to an H2 bound for the numerical solution.In addition,this H2 bound is not sufficient for the optimal rate convergence analysis,and we establish a uniform-in-time H3 bound for the numerical solution,based on the higher order Sobolev norm estimate,combined with repeated applications of discrete H¨older inequality and nonlinear embeddings in the Fourier pseudo-spectral space.This discrete H3 bound for the numerical solution enables us to derive the optimal rate error estimate for this alternate SAV method.A few numerical experiments are also presented,which confirm the efficiency and accuracy of the proposed scheme.展开更多
In this paper,we present a quadratic auxiliary variable(QAV)technique to develop a novel class of arbitrarily high-order energy-preserving algorithms for the Camassa-Holm equation.The QAV approach is first utilized to...In this paper,we present a quadratic auxiliary variable(QAV)technique to develop a novel class of arbitrarily high-order energy-preserving algorithms for the Camassa-Holm equation.The QAV approach is first utilized to transform the original equation into a reformulated QAV system with a consistent initial condition.Then the reformulated QAV system is discretized by applying the Fourier pseudo-spectral method in space and the symplectic Runge-Kutta methods in time,which arrives at a class of fully discrete schemes.Under the consistent initial condition,they can be rewritten as a new fully discrete system by eliminating the introduced auxiliary variable,which is rigorously proved to be energy-preserving and symmetric.Ample numerical experiments are conducted to confirm the expected order of accuracy,conservative property and efficiency of the proposed methods.The presented numerical strategy makes it possible to directly apply a special class of Runge-Kutta methods to develop energy-preserving algorithms for a general conservative system with any polynomial energy.展开更多
In the present paper,we propose an efficient scrambled estimator of population mean of quantitative sensitive study variable,using general linear transformation of nonsensitive auxiliary variable.Efficiency comparison...In the present paper,we propose an efficient scrambled estimator of population mean of quantitative sensitive study variable,using general linear transformation of nonsensitive auxiliary variable.Efficiency comparisons with the existing estimators have been carried out both theoretically and numerically.It has been found that our optimal scrambled estimator is always more efficient than most of the existing scrambled estimators and also it is more efficient than few other scrambled estimators under some conditions.展开更多
The second-order backward differential formula(BDF2)and the scalar auxiliary variable(SAV)approach are applied to con‐struct the linearly energy stable numerical scheme with the variable time steps for the epitaxial ...The second-order backward differential formula(BDF2)and the scalar auxiliary variable(SAV)approach are applied to con‐struct the linearly energy stable numerical scheme with the variable time steps for the epitaxial thin film growth models.Under the stepratio condition 0<τ_(n)/τ_(n-1)<4.864,the modified energy dissipation law is proven at the discrete levels with regardless of time step size.Nu‐merical experiments are presented to demonstrate the accuracy and efficiency of the proposed numerical scheme.展开更多
In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,t...In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.展开更多
In this work,we construct two efficient fully decoupled,linear,unconditionally stable numerical algorithms for the thermally coupled incompressible magnetohydrodynamic equations.Firstly,in order to obtain the desired ...In this work,we construct two efficient fully decoupled,linear,unconditionally stable numerical algorithms for the thermally coupled incompressible magnetohydrodynamic equations.Firstly,in order to obtain the desired algorithm,we introduce a scalar auxiliary variable(SAV)to get a new equivalent system.Secondly,by combining the pressure-correction method and the explicit-implicit method,we perform semi-discrete numerical algorithms of first and second order,respectively.Then,we prove that the obtained algorithms follow an unconditionally stable law in energy,and we provide a detailed implementation process,which we only need to solve a series of linear differential equations with constant coefficients at each time step.More importantly,with some powerful analysis,we give the order of convergence of the errors.Finally,to illustrate theoretical results,some numerical experiments are given.展开更多
In this paper,we construct a new class of efficient and high-order schemes for the Cahn-Hilliard-Navier-Stokes equations with periodic boundary conditions.These schemes are based on two types of scalar auxiliary varia...In this paper,we construct a new class of efficient and high-order schemes for the Cahn-Hilliard-Navier-Stokes equations with periodic boundary conditions.These schemes are based on two types of scalar auxiliary variable approaches.By using a new pressure correction method,the accuracy of the pressure has been greatly improved.Furthermore,one only needs to solve a series of fully decoupled linear equations with constant coefficients at each time step.In addition,we prove the unconditional energy stability of the schemes,rigorously.Finally,plenty of numerical simulations are carried out to verify the convergence rates,stability,and effectiveness of the proposed schemes numerically.展开更多
The spatial interpolation for soil texture does not necessarily satisfy the constant sum and nonnegativity constraints. Meanwhile, although numeric and categorical variables have been used as auxiliary variables to im...The spatial interpolation for soil texture does not necessarily satisfy the constant sum and nonnegativity constraints. Meanwhile, although numeric and categorical variables have been used as auxiliary variables to improve prediction accuracy of soil attributes such as soil organic matter, they （especially the categorical variables） are rarely used in spatial prediction of soil texture. The objective of our study was to comparing the performance of the methods for spatial prediction of soil texture with consideration of the characteristics of compositional data and auxiliary variables. These methods include the ordinary kriging with the symmetry logratio transform, regression kriging with the symmetry logratio transform, and compositional kriging （CK） approaches. The root mean squared error （RMSE）, the relative improvement value of RMSE and Aitchison＇s distance （DA） were all utilized to assess the accuracy of prediction and the mean squared deviation ratio was used to evaluate the goodness of fit of the theoretical estimate of error. The results showed that the prediction methods utilized in this paper could enable interpolation results of soil texture to satisfy the constant sum and nonnegativity constraints. Prediction accuracy and model fitting effect of the CK approach were better, suggesting that the CK method was more appropriate for predicting soil texture. The CK method is directly interpolated on soil texture, which ensures that it is optimal unbiased estimator. If the environment variables are appropriately selected as auxiliary variables, spatial variability of soil texture can be predicted reasonably and accordingly the predicted results will be satisfied.展开更多
In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-depe...In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.展开更多
In this paper, we propose an effective full array and sparse array adaptive beamforming scheme that can be applied for multiple desired signals based on the branch-and-bound algorithm. Adaptive beamforming for the mul...In this paper, we propose an effective full array and sparse array adaptive beamforming scheme that can be applied for multiple desired signals based on the branch-and-bound algorithm. Adaptive beamforming for the multiple desired signals is realized by the improved Capon method. At the same time,the sidelobe constraint is added to reduce the sidelobe level. To reduce the pointing errors of multiple desired signals, the array response phase of the desired signal is firstly optimized by using auxilary variables while keeping the response amplitude unchanged. The whole design is formulated as a convex optimization problem solved by the branch-and-bound algorithm. In addition,the beamformer weight vector is penalized with the modified reweighted l_(1)-norm to achieve sparsity. Theoretical analysis and simulation results show that the proposed algorithm has lower sidelobe level, higher SINR, and less pointing error than the stateof-the-art methods in the case of a single expected signal and multiple desired signals.展开更多
In this paper, we propose a class of estimators for estimating the finite population mean of the study variable under Ranked Set Sampling (RSS) when population mean of the auxiliary variable is known. The bias and Mea...In this paper, we propose a class of estimators for estimating the finite population mean of the study variable under Ranked Set Sampling (RSS) when population mean of the auxiliary variable is known. The bias and Mean Squared Error (MSE) of the proposed class of estimators are obtained to first degree of approximation. It is identified that the proposed class of estimators is more efficient as compared to [1] estimator and several other estimators. A simulation study is carried out to judge the performances of the estimators.展开更多
We present a decoupled,linearly implicit numerical scheme with energy stability and mass conservation for solving the coupled Cahn-Hilliard system.The time-discretization is done by leap-frog method with the scalar au...We present a decoupled,linearly implicit numerical scheme with energy stability and mass conservation for solving the coupled Cahn-Hilliard system.The time-discretization is done by leap-frog method with the scalar auxiliary variable(SAV)approach.It only needs to solve three linear equations at each time step,where each unknown variable can be solved independently.It is shown that the semi-discrete scheme has second-order accuracy in the temporal direction.Such convergence results are proved by a rigorous analysis of the boundedness of the numerical solution and the error estimates at different time-level.Numerical examples are presented to further confirm the validity of the methods.展开更多
In this paper,a linearized energy-stable scalar auxiliary variable(SAV)Galerkin scheme is investigated for a two-dimensional nonlinear wave equation and the unconditional superconvergence error estimates are obtained ...In this paper,a linearized energy-stable scalar auxiliary variable(SAV)Galerkin scheme is investigated for a two-dimensional nonlinear wave equation and the unconditional superconvergence error estimates are obtained without any certain time-step restrictions.The key to the analysis is to derive the boundedness of the numerical solution in the H^(1)-norm,which is different from the temporal-spatial error splitting approach used in the previous literature.Meanwhile,numerical results are provided to confirm the theoretical findings.展开更多
In this work, an unconditionally stable, decoupled, variable time step scheme is presentedfor the incompressible Navier-Stokes equations. Based on a scalar auxiliary variablein exponential function, this fully discret...In this work, an unconditionally stable, decoupled, variable time step scheme is presentedfor the incompressible Navier-Stokes equations. Based on a scalar auxiliary variablein exponential function, this fully discrete scheme combines the backward Euler schemefor temporal discretization with variable time step and a mixed finite element method forspatial discretization, where the nonlinear term is treated explicitly. Moreover, withoutany restriction on the time step, stability of the proposed scheme is discussed. Besides,error estimate is provided. Finally, some numerical results are presented to illustrate theperformances of the considered numerical scheme.展开更多
In this paper,we construct efficient schemes based on the scalar auxiliary variable block-centered finite difference method for the modified phase field crystal equation,which is a sixth-order nonlinear damped wave eq...In this paper,we construct efficient schemes based on the scalar auxiliary variable block-centered finite difference method for the modified phase field crystal equation,which is a sixth-order nonlinear damped wave equation.The schemes are linear,conserve mass and unconditionally dissipate a pseudo energy.We prove rigorously second-order error estimates in both time and space for the phase field variable in discrete norms.We also present some numerical experiments to verify our theoretical results and demonstrate the robustness and accuracy.展开更多
基金Supported by the National Natural Science Foundation of China(Grant Nos.12001210 and 12261103)the Natural Science Foundation of Henan(Grant No.252300420308)the Yunnan Fundamental Research Projects(Grant No.202301AT070117).
文摘An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into an equivalent system,and the k-order backward differentiation formula(BDF k)and central difference formula are used to discretize the temporal and spatial derivatives,respectively.Different from the traditional discrete method that adopts full implicit or full explicit for the nonlinear integral terms,the proposed scheme is based on the SAV idea and can be treated semi-implicitly,taking into account both accuracy and effectiveness.Numerical results are presented to demonstrate the high-order convergence(up to fourth-order)of the developed schemes and it is computationally efficient in long-time computations.
基金supported by the National Natural Science Foundation of China(12126318,12126302).
文摘This paper develops a generalized scalar auxiliary variable(SAV)method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent Ginzburg-Landau equations.In this method,the system is decoupled and linearized to avoid solving the non-linear equation at each step.The theoretical analysis proves that the generalized SAV method can preserve the maximum bound principle and energy stability,and this is confirmed by the numerical result,and also shows that the numerical algorithm is stable.
文摘This manuscript deals with new class of almost unbiased ratio cum product estimators for the estimation of population mean of the study variable by using the known values auxiliary variable. The bias and mean squared error of proposed estimators are obtained. An empirical study is carried out to assess the efficiency of proposed estimators over the existing estimators with the help of some known natural populations and it shows that the proposed estimators are almost unbiased and it perform better than the existing estimators.
文摘In this paper, we have proposed estimators of finite population mean using generalized Ratio- cum-product estimator for two-Phase sampling using multi-auxiliary variables under full, partial and no information cases and investigated their finite sample properties. An empirical study is given to compare the performance of the proposed estimators with the existing estimators that utilize auxiliary variable(s) for finite population mean. It has been found that the generalized Ra-tio-cum-product estimator in full information case using multiple auxiliary variables is more efficient than mean per unit, ratio and product estimator using one auxiliary variable, ratio and product estimator using multiple auxiliary variable and ratio-cum-product estimators in both partial and no information case in two phase sampling. A generalized Ratio-cum-product estimator in partial information case is more efficient than Generalized Ratio-cum-product estimator in No information case.
基金supported by NSFC(Grant No.11871441)supported by NSF(Grant No.DMS-2012669).
文摘In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is introduced for the physical system;the mass conservation and energy dissipation is proved for the CHHS system.Subsequently,a fully discrete SAV finite element scheme is proposed,with the mass conservation and energy dissipation laws established at a theoretical level.In addition,the convergence analysis and error estimate is provided for the proposed SAV numerical scheme.
文摘A second order accurate(in time)numerical scheme is analyzed for the slope-selection(SS)equation of the epitaxial thin film growth model,with Fourier pseudo-spectral discretization in space.To make the numerical scheme linear while preserving the nonlinear energy stability,we make use of the scalar auxiliary variable(SAV)approach,in which a modified Crank-Nicolson is applied for the surface diffusion part.The energy stability could be derived a modified form,in comparison with the standard Crank-Nicolson approximation to the surface diffusion term.Such an energy stability leads to an H2 bound for the numerical solution.In addition,this H2 bound is not sufficient for the optimal rate convergence analysis,and we establish a uniform-in-time H3 bound for the numerical solution,based on the higher order Sobolev norm estimate,combined with repeated applications of discrete H¨older inequality and nonlinear embeddings in the Fourier pseudo-spectral space.This discrete H3 bound for the numerical solution enables us to derive the optimal rate error estimate for this alternate SAV method.A few numerical experiments are also presented,which confirm the efficiency and accuracy of the proposed scheme.
基金Yuezheng Gong’s work is partially supported by the Foundation of Jiangsu Key Laboratory for Numerical Simulation of Large Scale Complex Systems(Grant No.202002)the Fundamental Research Funds for the Central Universities(Grant No.NS2022070)+7 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20220131)the National Natural Science Foundation of China(Grants Nos.12271252 and 12071216)Qi Hong’s work is partially supported by the National Natural Science Foundation of China(Grants No.12201297)the Foundation of Jiangsu Key Laboratory for Numerical Simulation of Large Scale Complex Systems(Grant No.202001)Chunwu Wang’s work is partially supported by Science Challenge Project(Grant No.TZ2018002)National Science and Technology Major Project(J2019-II-0007-0027)Yushun Wang’s work is partially supported by the National Key Research and Development Program of China(Grant No.2018YFC1504205)the National Natural Science Foundation of China(Grants No.12171245).
文摘In this paper,we present a quadratic auxiliary variable(QAV)technique to develop a novel class of arbitrarily high-order energy-preserving algorithms for the Camassa-Holm equation.The QAV approach is first utilized to transform the original equation into a reformulated QAV system with a consistent initial condition.Then the reformulated QAV system is discretized by applying the Fourier pseudo-spectral method in space and the symplectic Runge-Kutta methods in time,which arrives at a class of fully discrete schemes.Under the consistent initial condition,they can be rewritten as a new fully discrete system by eliminating the introduced auxiliary variable,which is rigorously proved to be energy-preserving and symmetric.Ample numerical experiments are conducted to confirm the expected order of accuracy,conservative property and efficiency of the proposed methods.The presented numerical strategy makes it possible to directly apply a special class of Runge-Kutta methods to develop energy-preserving algorithms for a general conservative system with any polynomial energy.
文摘In the present paper,we propose an efficient scrambled estimator of population mean of quantitative sensitive study variable,using general linear transformation of nonsensitive auxiliary variable.Efficiency comparisons with the existing estimators have been carried out both theoretically and numerically.It has been found that our optimal scrambled estimator is always more efficient than most of the existing scrambled estimators and also it is more efficient than few other scrambled estimators under some conditions.
文摘The second-order backward differential formula(BDF2)and the scalar auxiliary variable(SAV)approach are applied to con‐struct the linearly energy stable numerical scheme with the variable time steps for the epitaxial thin film growth models.Under the stepratio condition 0<τ_(n)/τ_(n-1)<4.864,the modified energy dissipation law is proven at the discrete levels with regardless of time step size.Nu‐merical experiments are presented to demonstrate the accuracy and efficiency of the proposed numerical scheme.
基金Supported by Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)+2 种基金Basic Research Plan of Shanxi Province(202203021211129)Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)。
文摘In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.
基金Supported by Research Project Supported by Shanxi Scholarship Council of China(2021-029)Shanxi Provincial International Cooperation Base and Platform Project(202104041101019)Shanxi Province Natural Science Foundation(202203021211129)。
文摘In this work,we construct two efficient fully decoupled,linear,unconditionally stable numerical algorithms for the thermally coupled incompressible magnetohydrodynamic equations.Firstly,in order to obtain the desired algorithm,we introduce a scalar auxiliary variable(SAV)to get a new equivalent system.Secondly,by combining the pressure-correction method and the explicit-implicit method,we perform semi-discrete numerical algorithms of first and second order,respectively.Then,we prove that the obtained algorithms follow an unconditionally stable law in energy,and we provide a detailed implementation process,which we only need to solve a series of linear differential equations with constant coefficients at each time step.More importantly,with some powerful analysis,we give the order of convergence of the errors.Finally,to illustrate theoretical results,some numerical experiments are given.
基金Supported by the Research Project Supported of Shanxi Scholarship Council of China(No.2021-029)Shanxi Provincial International Cooperation Base and Platform Project(202104041101019)Shanxi Province Natural Science Research(202203021211129)。
文摘In this paper,we construct a new class of efficient and high-order schemes for the Cahn-Hilliard-Navier-Stokes equations with periodic boundary conditions.These schemes are based on two types of scalar auxiliary variable approaches.By using a new pressure correction method,the accuracy of the pressure has been greatly improved.Furthermore,one only needs to solve a series of fully decoupled linear equations with constant coefficients at each time step.In addition,we prove the unconditional energy stability of the schemes,rigorously.Finally,plenty of numerical simulations are carried out to verify the convergence rates,stability,and effectiveness of the proposed schemes numerically.
基金supported by the National Natural Science Foundation of China (41071152)the Special Fund for Land and Resources Scientific Research in the Public Interest,China (201011006-3)the Special Fund for Agro-Scientific Research in the Public Interest,China (201103005-01-01)
文摘The spatial interpolation for soil texture does not necessarily satisfy the constant sum and nonnegativity constraints. Meanwhile, although numeric and categorical variables have been used as auxiliary variables to improve prediction accuracy of soil attributes such as soil organic matter, they （especially the categorical variables） are rarely used in spatial prediction of soil texture. The objective of our study was to comparing the performance of the methods for spatial prediction of soil texture with consideration of the characteristics of compositional data and auxiliary variables. These methods include the ordinary kriging with the symmetry logratio transform, regression kriging with the symmetry logratio transform, and compositional kriging （CK） approaches. The root mean squared error （RMSE）, the relative improvement value of RMSE and Aitchison＇s distance （DA） were all utilized to assess the accuracy of prediction and the mean squared deviation ratio was used to evaluate the goodness of fit of the theoretical estimate of error. The results showed that the prediction methods utilized in this paper could enable interpolation results of soil texture to satisfy the constant sum and nonnegativity constraints. Prediction accuracy and model fitting effect of the CK approach were better, suggesting that the CK method was more appropriate for predicting soil texture. The CK method is directly interpolated on soil texture, which ensures that it is optimal unbiased estimator. If the environment variables are appropriately selected as auxiliary variables, spatial variability of soil texture can be predicted reasonably and accordingly the predicted results will be satisfied.
文摘In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.
基金the National Key Research and Development Program(2021YFB3502500).
文摘In this paper, we propose an effective full array and sparse array adaptive beamforming scheme that can be applied for multiple desired signals based on the branch-and-bound algorithm. Adaptive beamforming for the multiple desired signals is realized by the improved Capon method. At the same time,the sidelobe constraint is added to reduce the sidelobe level. To reduce the pointing errors of multiple desired signals, the array response phase of the desired signal is firstly optimized by using auxilary variables while keeping the response amplitude unchanged. The whole design is formulated as a convex optimization problem solved by the branch-and-bound algorithm. In addition,the beamformer weight vector is penalized with the modified reweighted l_(1)-norm to achieve sparsity. Theoretical analysis and simulation results show that the proposed algorithm has lower sidelobe level, higher SINR, and less pointing error than the stateof-the-art methods in the case of a single expected signal and multiple desired signals.
文摘In this paper, we propose a class of estimators for estimating the finite population mean of the study variable under Ranked Set Sampling (RSS) when population mean of the auxiliary variable is known. The bias and Mean Squared Error (MSE) of the proposed class of estimators are obtained to first degree of approximation. It is identified that the proposed class of estimators is more efficient as compared to [1] estimator and several other estimators. A simulation study is carried out to judge the performances of the estimators.
基金supported by the National Natural Science Foundation of China(Grant Nos.12171442,12231003,12271215,12326378,11871248).
文摘We present a decoupled,linearly implicit numerical scheme with energy stability and mass conservation for solving the coupled Cahn-Hilliard system.The time-discretization is done by leap-frog method with the scalar auxiliary variable(SAV)approach.It only needs to solve three linear equations at each time step,where each unknown variable can be solved independently.It is shown that the semi-discrete scheme has second-order accuracy in the temporal direction.Such convergence results are proved by a rigorous analysis of the boundedness of the numerical solution and the error estimates at different time-level.Numerical examples are presented to further confirm the validity of the methods.
基金supported by the National Natural Science Foundation of China(No.12101568).
文摘In this paper,a linearized energy-stable scalar auxiliary variable(SAV)Galerkin scheme is investigated for a two-dimensional nonlinear wave equation and the unconditional superconvergence error estimates are obtained without any certain time-step restrictions.The key to the analysis is to derive the boundedness of the numerical solution in the H^(1)-norm,which is different from the temporal-spatial error splitting approach used in the previous literature.Meanwhile,numerical results are provided to confirm the theoretical findings.
基金supported by the Natural Science Foundation of China(Grant No.12361077)by the Tianshan Talent Training Program of Xinjiang Uygur Autonomous Region(Grant No.2023TSYCCX0103)by the Natural Science Foundation of Xinjiang Uygur Autonomous Region(Grant No.2023D14014).
文摘In this work, an unconditionally stable, decoupled, variable time step scheme is presentedfor the incompressible Navier-Stokes equations. Based on a scalar auxiliary variablein exponential function, this fully discrete scheme combines the backward Euler schemefor temporal discretization with variable time step and a mixed finite element method forspatial discretization, where the nonlinear term is treated explicitly. Moreover, withoutany restriction on the time step, stability of the proposed scheme is discussed. Besides,error estimate is provided. Finally, some numerical results are presented to illustrate theperformances of the considered numerical scheme.
基金supported by National Natural Science Foundation of China(Grant Nos.11901489 and 11971407)supported by National Science Foundation of USA(Grant No.DMS-1720442)。
文摘In this paper,we construct efficient schemes based on the scalar auxiliary variable block-centered finite difference method for the modified phase field crystal equation,which is a sixth-order nonlinear damped wave equation.The schemes are linear,conserve mass and unconditionally dissipate a pseudo energy.We prove rigorously second-order error estimates in both time and space for the phase field variable in discrete norms.We also present some numerical experiments to verify our theoretical results and demonstrate the robustness and accuracy.
