The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part o...The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part of the operator and an implicit Euler discretization of the linear part. Finite difference schemes are used for the spatial part. This finally leads to the numerical solution of a sparse linear system that can be solved efficiently.展开更多
The calculation sequence of collision, propagation and macroscopic variables is not very clear in lattice Boltzmann method (LBM) code implementation. According to the definition, three steps should be computed on all ...The calculation sequence of collision, propagation and macroscopic variables is not very clear in lattice Boltzmann method (LBM) code implementation. According to the definition, three steps should be computed on all nodes respectively, which mean three loops are needed. While the “pull” scheme makes the only one loop possible for coding, this is called semi-implicit scheme in this study. The accuracy and efficiency of semi-implicit scheme are discussed in detail through the simulation of lid-driven cavity flow. Non-equilibrium extrapolation scheme is adopted on the boundary of simulation area. The results are compared with two classic articles, which show that semi-implicit scheme has good agreement with the classic scheme. When Re is less than 3000, the iterations steps of semi-scheme can be decreased by about 30% though comparing the semi-implicit scheme with standard scheme containing three loops. As the Re increases into more than 3400, the standard scheme is not converged. On the contrary, the iterations of semi-implicit scheme are approximately linear to Re.展开更多
The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearitie...The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearities.To discretize the continuous problem in space,we utilize a spectral Galerkin method.Subsequently,we introduce a nonlinear-tamed exponential integrator scheme,resulting in a fully discrete scheme.Within the framework of semigroup theory,this study provides precise estimations of the Sobolev regularity,L^(∞) regularity in space,and Hölder continuity in time for the mild solution,as well as for its semi-discrete and full-discrete approximations.Building upon these results,we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions.A numerical example is presented to validate the theoretical findings.展开更多
As blockchain technology rapidly evolves,smart contracts have seen widespread adoption in financial transactions and beyond.However,the growing prevalence of malicious Ponzi scheme contracts presents serious security ...As blockchain technology rapidly evolves,smart contracts have seen widespread adoption in financial transactions and beyond.However,the growing prevalence of malicious Ponzi scheme contracts presents serious security threats to blockchain ecosystems.Although numerous detection techniques have been proposed,existing methods suffer from significant limitations,such as class imbalance and insufficient modeling of transaction-related semantic features.To address these challenges,this paper proposes an oversampling-based detection framework for Ponzi smart contracts.We enhance the Adaptive Synthetic Sampling(ADASYN)algorithm by incorporating sample proximity to decision boundaries and ensuring realistic sample distributions.This enhancement facilitates the generation of high-quality minority class samples and effectively mitigates class imbalance.In addition,we design a Contract Transaction Graph(CTG)construction algorithm to preserve key transactional semantics through feature extraction from contract code.A graph neural network(GNN)is then applied for classification.This study employs a publicly available dataset from the XBlock platform,consisting of 318 verified Ponzi contracts and 6498 benign contracts.Sourced from real Ethereum deployments,the dataset reflects diverse application scenarios and captures the varied characteristics of Ponzi schemes.Experimental results demonstrate that our approach achieves an accuracy of 96%,a recall of 92%,and an F1-score of 94%in detecting Ponzi contracts,outperforming state-of-the-art methods.展开更多
Clouds play an important role in global atmospheric energy and water vapor budgets, and the low cloud simulations suffer from large biases in many atmospheric general circulation models. In this study, cloud microphys...Clouds play an important role in global atmospheric energy and water vapor budgets, and the low cloud simulations suffer from large biases in many atmospheric general circulation models. In this study, cloud microphysical processes such as raindrop evaporation and cloud water accretion in a double-moment six-class cloud microphysics scheme were revised to enhance the simulation of low clouds using the Global-Regional Integrated Forecast System(GRIST)model. The validation of the revised scheme using a single-column version of the GRIST demonstrated a reasonable reduction in liquid water biases. The revised parameterization simulated medium-and low-level cloud fractions that were in better agreement with the observations than the original scheme. Long-term global simulations indicate the mitigation of the originally overestimated low-level cloud fraction and cloud-water mixing ratio in mid-to high-latitude regions,primarily owing to enhanced accretion processes and weakened raindrop evaporation. The reduced low clouds with the revised scheme showed better consistency with satellite observations, particularly at mid-and high-latitudes. Further improvements can be observed in the simulated cloud shortwave radiative forcing and vertical distribution of total cloud cover. Annual precipitation in mid-latitude regions has also improved, particularly over the oceans, with significantly increased large-scale and decreased convective precipitation.展开更多
In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is de...In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.展开更多
This study proposes a class of augmented subspace schemes for the weak Galerkin(WG)finite element method used to solve eigenvalue problems.The augmented subspace is built with the conforming linear finite element spac...This study proposes a class of augmented subspace schemes for the weak Galerkin(WG)finite element method used to solve eigenvalue problems.The augmented subspace is built with the conforming linear finite element space defined on the coarse mesh and the eigen-function approximations in the WG finite element space defined on the fine mesh.Based on this augmented subspace,solving the eigenvalue problem in the fine WG finite element space can be reduced to the solution of the linear boundary value problem in the same WG finite element space and a low dimensional eigenvalue problem in the augmented sub-space.The proposed augmented subspace techniques have the second order convergence rate with respect to the coarse mesh size,as demonstrated by the accompanying error esti-mates.Finally,a few numerical examples are provided to validate the proposed numerical techniques.展开更多
In this paper,we propose and analyze two second-order accurate finite difference schemes for the one-dimensional heat equation with concentrated capacity on a computa-tional domain=[a,b].We first transform the target ...In this paper,we propose and analyze two second-order accurate finite difference schemes for the one-dimensional heat equation with concentrated capacity on a computa-tional domain=[a,b].We first transform the target equation into the standard heat equation on the domain excluding the singular point equipped with an inner interface matching(IIM)condition on the singular point x=ξ∈(a,b),then adopt Taylor’s ex-pansion to approximate the IIM condition at the singular point and apply second-order finite difference method to approximate the standard heat equation at the nonsingular points.This discrete procedure allows us to choose different grid sizes to partition the two sub-domains[a,ξ]and[ξ,b],which ensures that x=ξ is a grid point,and hence the pro-posed schemes can be generalized to the heat equation with more than one concentrated capacities.We prove that the two proposed schemes are uniquely solvable.And through in-depth analysis of the local truncation errors,we rigorously prove that the two schemes are second-order accurate both in temporal and spatial directions in the maximum norm without any constraint on the grid ratio.Numerical experiments are carried out to verify our theoretical conclusions.展开更多
This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called...This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.展开更多
It is a very common practice to use semi-implicit schemes in various computations,which treat selected linear terms implicitly and the nonlinear terms explicitly.For phase-field equations,the principal elliptic operat...It is a very common practice to use semi-implicit schemes in various computations,which treat selected linear terms implicitly and the nonlinear terms explicitly.For phase-field equations,the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step.However,very few recent numerical analysis is relevant to semi-implicit schemes,while”stabilized”schemes have become very popular.In this work,we will consider semiimplicit schemes for the Allen-Cahn equation with general potential function.It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions.This paper extends the result of Tang&Yang(J.Comput.Math.,34(5)(2016),pp.471–481),which studies the semi-implicit scheme for the Allen-Cahn equation with polynomial potentials.展开更多
We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equati...We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation （image selective smoothing model） given by Alvarez et al. （Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29： 845-866）. Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1＇2 （W1＇1） sense in isotropic （anisotropic） diffu- sion domain.展开更多
The regenerative braking energy utilization system(RBEUS)stands as a promising technique for improving the efficiency and power quality of electrified railways.Beyond the vital aspects of energy management and control...The regenerative braking energy utilization system(RBEUS)stands as a promising technique for improving the efficiency and power quality of electrified railways.Beyond the vital aspects of energy management and control strategies,ensuring fault protection is paramount for the secure and steady operation of the traction power supply system(TPSS)integrated with RBEUS.This paper introduces an innovative protection scheme tailored to diverse RBEUS application scenarios.Firstly,fault categories are streamlined into three levels:system,equipment,and warning.Subsequently,a novel multi-port active power differential protection method,aligned with RBEUS operational principles,is crafted to serve as a comprehensive and sensitive main protection.Building upon this foundation,a hierarchical protection structure for RBEUS is established,addressing the intricacies and variations in fault types while boosting anti-disturbance capabilities under faulty conditions.Embracing the principle of railway-oriented safety,a collaborative RBEUS-TPSS protection scheme is put forth.Finally,through simulated scenarios encompassing various fault conditions,the proposed scheme’s feasibility and effectiveness are convincingly validated.展开更多
BACKGROUND Elderly patients undergoing laparoscopic colorectal cancer surgery are at high risk for hypothermia-related complications.This study explores the efficacy of periop-erative composite insulation intervention...BACKGROUND Elderly patients undergoing laparoscopic colorectal cancer surgery are at high risk for hypothermia-related complications.This study explores the efficacy of periop-erative composite insulation interventions in maintaining normothermia and reducing postoperative risks in this vulnerable group.AIM To evaluate the efficacy of perioperative composite insulation in older patients undergoing colorectal cancer surgery.METHODS We selected 100 older patients who underwent laparoscopic surgery for colorectal cancer at Huzhou Central Hospital from September 2023 to April 2024.Using a random number table,patients were divided into a control group and inter-vention group of 50 patients each.After returning to the regular ward,the con-ventional group received traditional insulation intervention measures,while the intervention group received composite insulation nursing intervention.We ob-served and recorded postoperative blood pressure and heart rate changes,as well as postoperative anesthesia recovery time and incidence of complications.RESULTS The statistical results showed significant differences(P<0.05)in heart rate changes and systolic blood pressure between the two groups.There was a sig-nificant change in heart rate between the groups immediately after surgery and at 15 and 30 minutes after surgery(P<0.05).The heart rate and systolic blood pressure of the intervention group were significantly lower than those of the control group at 15 and 30 minutes after surgery(P<0.05).The rewarming time of the intervention group was shorter than that of the control group,and the overall incidence of postoperative complications was significantly lower than that of the control group(P<0.05).CONCLUSION For elderly patients undergoing laparoscopic colorectal cancer surgery,a composite insulation intervention during the perioperative period can maintain body temperature,reduce postoperative stress,and significantly reduce the incidence of hypothermia and related complications.展开更多
Multiple quantum well(MQW) Ⅲ-nitride diodes that can simultaneously emit and detect light feature an overlapping region between their electroluminescence and responsivity spectra, which allows them to be simultaneous...Multiple quantum well(MQW) Ⅲ-nitride diodes that can simultaneously emit and detect light feature an overlapping region between their electroluminescence and responsivity spectra, which allows them to be simultaneously used as both a transmitter and a receiver in a wireless light communication system. Here, we demonstrate a mobile light communication system using a time-division multiplexing(TDM) scheme to achieve bidirectional data transmission via the same optical channel.Two identical blue MQW diodes are defined by software as a transmitter or a receiver. To address the light alignment issue, an image identification module integrated with a gimbal stabilizer is used to automatically detect the locations of moving targets;thus, underwater audio communication is realized via a mobile blue-light TDM communication mode. This approach not only uses a single link but also integrates mobile nodes in a practical network.展开更多
This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunzia...This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.展开更多
ADER-WAF methods were first introduced by researchers E.F. Toro and V.A. Titarev. The linear stability criterion for the model equation for the ADER-WAF schemes is CCFL≤1, where CCFLdenotes the Courant-Friedrichs-Lew...ADER-WAF methods were first introduced by researchers E.F. Toro and V.A. Titarev. The linear stability criterion for the model equation for the ADER-WAF schemes is CCFL≤1, where CCFLdenotes the Courant-Friedrichs-Lewy (CFL) coefficient. Toro and Titarev employed CCFL=0.95for their experiments. Nonetheless, we noted that the experiments conducted in this study with CCFL=0.95produced solutions exhibiting spurious oscillations, particularly in the high-order ADER-WAF schemes. The homogeneous one-dimensional (1D) non-linear Shallow Water Equations (SWEs) are the subject of these experiments, specifically the solution of the Riemann Problem (RP) associated with the SWEs. The investigation was conducted on four test problems to evaluate the ADER-WAF schemes of second, third, fourth, and fifth order of accuracy. Each test problem constitutes a RP characterized by different wave patterns in its solution. This research has two primary objectives. We begin by illustrating the procedure for implementing the ADER-WAF schemes for the SWEs, providing the required relations. Afterward, following comprehensive testing, we present the range for the CFL coefficient for each test that yields solutions with diminished or eliminated spurious oscillations.展开更多
The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on mo...The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.展开更多
In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interfac...In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.展开更多
文摘The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part of the operator and an implicit Euler discretization of the linear part. Finite difference schemes are used for the spatial part. This finally leads to the numerical solution of a sparse linear system that can be solved efficiently.
文摘The calculation sequence of collision, propagation and macroscopic variables is not very clear in lattice Boltzmann method (LBM) code implementation. According to the definition, three steps should be computed on all nodes respectively, which mean three loops are needed. While the “pull” scheme makes the only one loop possible for coding, this is called semi-implicit scheme in this study. The accuracy and efficiency of semi-implicit scheme are discussed in detail through the simulation of lid-driven cavity flow. Non-equilibrium extrapolation scheme is adopted on the boundary of simulation area. The results are compared with two classic articles, which show that semi-implicit scheme has good agreement with the classic scheme. When Re is less than 3000, the iterations steps of semi-scheme can be decreased by about 30% though comparing the semi-implicit scheme with standard scheme containing three loops. As the Re increases into more than 3400, the standard scheme is not converged. On the contrary, the iterations of semi-implicit scheme are approximately linear to Re.
基金partially supported by the National Natural Science Foundation of China(Grant No.12071073)financial support by the Jiangsu Provincial Scientific Research Center of Applied Mathematics(Grant No.BK20233002).
文摘The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearities.To discretize the continuous problem in space,we utilize a spectral Galerkin method.Subsequently,we introduce a nonlinear-tamed exponential integrator scheme,resulting in a fully discrete scheme.Within the framework of semigroup theory,this study provides precise estimations of the Sobolev regularity,L^(∞) regularity in space,and Hölder continuity in time for the mild solution,as well as for its semi-discrete and full-discrete approximations.Building upon these results,we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions.A numerical example is presented to validate the theoretical findings.
基金supported by the Key Project of Joint Fund of the National Natural Science Foundation of China“Research on Key Technologies and Demonstration Applications for Trusted and Secure Data Circulation and Trading”(U24A20241)the National Natural Science Foundation of China“Research on Trusted Theories and Key Technologies of Data Security Trading Based on Blockchain”(62202118)+4 种基金the Major Scientific and Technological Special Project of Guizhou Province([2024]014)Scientific and Technological Research Projects from the Guizhou Education Department(Qian jiao ji[2023]003)the Hundred-Level Innovative Talent Project of the Guizhou Provincial Science and Technology Department(Qiankehe Platform Talent-GCC[2023]018)the Major Project of Guizhou Province“Research and Application of Key Technologies for Trusted Large Models Oriented to Public Big Data”(Qiankehe Major Project[2024]003)the Guizhou Province Computational Power Network Security Protection Science and Technology Innovation Talent Team(Qiankehe Talent CXTD[2025]029).
文摘As blockchain technology rapidly evolves,smart contracts have seen widespread adoption in financial transactions and beyond.However,the growing prevalence of malicious Ponzi scheme contracts presents serious security threats to blockchain ecosystems.Although numerous detection techniques have been proposed,existing methods suffer from significant limitations,such as class imbalance and insufficient modeling of transaction-related semantic features.To address these challenges,this paper proposes an oversampling-based detection framework for Ponzi smart contracts.We enhance the Adaptive Synthetic Sampling(ADASYN)algorithm by incorporating sample proximity to decision boundaries and ensuring realistic sample distributions.This enhancement facilitates the generation of high-quality minority class samples and effectively mitigates class imbalance.In addition,we design a Contract Transaction Graph(CTG)construction algorithm to preserve key transactional semantics through feature extraction from contract code.A graph neural network(GNN)is then applied for classification.This study employs a publicly available dataset from the XBlock platform,consisting of 318 verified Ponzi contracts and 6498 benign contracts.Sourced from real Ethereum deployments,the dataset reflects diverse application scenarios and captures the varied characteristics of Ponzi schemes.Experimental results demonstrate that our approach achieves an accuracy of 96%,a recall of 92%,and an F1-score of 94%in detecting Ponzi contracts,outperforming state-of-the-art methods.
基金National Natural Science Foundation of China(42375153,42105153,42205157)Development of Science and Technology at Chinese Academy of Meteorological Sciences(2023KJ038)。
文摘Clouds play an important role in global atmospheric energy and water vapor budgets, and the low cloud simulations suffer from large biases in many atmospheric general circulation models. In this study, cloud microphysical processes such as raindrop evaporation and cloud water accretion in a double-moment six-class cloud microphysics scheme were revised to enhance the simulation of low clouds using the Global-Regional Integrated Forecast System(GRIST)model. The validation of the revised scheme using a single-column version of the GRIST demonstrated a reasonable reduction in liquid water biases. The revised parameterization simulated medium-and low-level cloud fractions that were in better agreement with the observations than the original scheme. Long-term global simulations indicate the mitigation of the originally overestimated low-level cloud fraction and cloud-water mixing ratio in mid-to high-latitude regions,primarily owing to enhanced accretion processes and weakened raindrop evaporation. The reduced low clouds with the revised scheme showed better consistency with satellite observations, particularly at mid-and high-latitudes. Further improvements can be observed in the simulated cloud shortwave radiative forcing and vertical distribution of total cloud cover. Annual precipitation in mid-latitude regions has also improved, particularly over the oceans, with significantly increased large-scale and decreased convective precipitation.
基金The Natural Science Foundation of Xinjiang Uygur Autonomous Region of China“RBF-Hermite difference scheme for the time-fractional kdv-Burgers equation”(2024D01C43)。
文摘In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.
基金partly supported by the Beijing Natural Science Foundation(Grant No.Z200003)by the National Natural Science Foundation of China(Grant Nos.12331015,12301475,12301465)+1 种基金by the National Center for Mathematics and Interdisciplinary Science,Chinese Academy of Sciencesby the Research Foundation for the Beijing University of Technology New Faculty(Grant No.006000514122516).
文摘This study proposes a class of augmented subspace schemes for the weak Galerkin(WG)finite element method used to solve eigenvalue problems.The augmented subspace is built with the conforming linear finite element space defined on the coarse mesh and the eigen-function approximations in the WG finite element space defined on the fine mesh.Based on this augmented subspace,solving the eigenvalue problem in the fine WG finite element space can be reduced to the solution of the linear boundary value problem in the same WG finite element space and a low dimensional eigenvalue problem in the augmented sub-space.The proposed augmented subspace techniques have the second order convergence rate with respect to the coarse mesh size,as demonstrated by the accompanying error esti-mates.Finally,a few numerical examples are provided to validate the proposed numerical techniques.
基金supported by the National Natural Science Foundation of China(Grant No.11571181)by the Natural Science Foundation of Jiangsu Province(Grant No.BK20171454).
文摘In this paper,we propose and analyze two second-order accurate finite difference schemes for the one-dimensional heat equation with concentrated capacity on a computa-tional domain=[a,b].We first transform the target equation into the standard heat equation on the domain excluding the singular point equipped with an inner interface matching(IIM)condition on the singular point x=ξ∈(a,b),then adopt Taylor’s ex-pansion to approximate the IIM condition at the singular point and apply second-order finite difference method to approximate the standard heat equation at the nonsingular points.This discrete procedure allows us to choose different grid sizes to partition the two sub-domains[a,ξ]and[ξ,b],which ensures that x=ξ is a grid point,and hence the pro-posed schemes can be generalized to the heat equation with more than one concentrated capacities.We prove that the two proposed schemes are uniquely solvable.And through in-depth analysis of the local truncation errors,we rigorously prove that the two schemes are second-order accurate both in temporal and spatial directions in the maximum norm without any constraint on the grid ratio.Numerical experiments are carried out to verify our theoretical conclusions.
文摘This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.
文摘It is a very common practice to use semi-implicit schemes in various computations,which treat selected linear terms implicitly and the nonlinear terms explicitly.For phase-field equations,the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step.However,very few recent numerical analysis is relevant to semi-implicit schemes,while”stabilized”schemes have become very popular.In this work,we will consider semiimplicit schemes for the Allen-Cahn equation with general potential function.It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions.This paper extends the result of Tang&Yang(J.Comput.Math.,34(5)(2016),pp.471–481),which studies the semi-implicit scheme for the Allen-Cahn equation with polynomial potentials.
文摘We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation （image selective smoothing model） given by Alvarez et al. （Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29： 845-866）. Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1＇2 （W1＇1） sense in isotropic （anisotropic） diffu- sion domain.
基金supported by the National Natural Science Foundation of China(Nos.52107126 and52077179)the Key Regional Innovation and Development Joint Fund Project(No.2023YFB2303901)the funding of Chengdu Guojia Electrical Engineering Co.,Ltd.(No.NEEC-2022-B11).
文摘The regenerative braking energy utilization system(RBEUS)stands as a promising technique for improving the efficiency and power quality of electrified railways.Beyond the vital aspects of energy management and control strategies,ensuring fault protection is paramount for the secure and steady operation of the traction power supply system(TPSS)integrated with RBEUS.This paper introduces an innovative protection scheme tailored to diverse RBEUS application scenarios.Firstly,fault categories are streamlined into three levels:system,equipment,and warning.Subsequently,a novel multi-port active power differential protection method,aligned with RBEUS operational principles,is crafted to serve as a comprehensive and sensitive main protection.Building upon this foundation,a hierarchical protection structure for RBEUS is established,addressing the intricacies and variations in fault types while boosting anti-disturbance capabilities under faulty conditions.Embracing the principle of railway-oriented safety,a collaborative RBEUS-TPSS protection scheme is put forth.Finally,through simulated scenarios encompassing various fault conditions,the proposed scheme’s feasibility and effectiveness are convincingly validated.
文摘BACKGROUND Elderly patients undergoing laparoscopic colorectal cancer surgery are at high risk for hypothermia-related complications.This study explores the efficacy of periop-erative composite insulation interventions in maintaining normothermia and reducing postoperative risks in this vulnerable group.AIM To evaluate the efficacy of perioperative composite insulation in older patients undergoing colorectal cancer surgery.METHODS We selected 100 older patients who underwent laparoscopic surgery for colorectal cancer at Huzhou Central Hospital from September 2023 to April 2024.Using a random number table,patients were divided into a control group and inter-vention group of 50 patients each.After returning to the regular ward,the con-ventional group received traditional insulation intervention measures,while the intervention group received composite insulation nursing intervention.We ob-served and recorded postoperative blood pressure and heart rate changes,as well as postoperative anesthesia recovery time and incidence of complications.RESULTS The statistical results showed significant differences(P<0.05)in heart rate changes and systolic blood pressure between the two groups.There was a sig-nificant change in heart rate between the groups immediately after surgery and at 15 and 30 minutes after surgery(P<0.05).The heart rate and systolic blood pressure of the intervention group were significantly lower than those of the control group at 15 and 30 minutes after surgery(P<0.05).The rewarming time of the intervention group was shorter than that of the control group,and the overall incidence of postoperative complications was significantly lower than that of the control group(P<0.05).CONCLUSION For elderly patients undergoing laparoscopic colorectal cancer surgery,a composite insulation intervention during the perioperative period can maintain body temperature,reduce postoperative stress,and significantly reduce the incidence of hypothermia and related complications.
基金jointly supported by the National Natural Science Foundation of China (U21A20495)Natural Science Foundation of Jiangsu Province (BG2024023)+1 种基金National Key Research and Development Program of China (2022YFE0112000)111 Project (D17018)。
文摘Multiple quantum well(MQW) Ⅲ-nitride diodes that can simultaneously emit and detect light feature an overlapping region between their electroluminescence and responsivity spectra, which allows them to be simultaneously used as both a transmitter and a receiver in a wireless light communication system. Here, we demonstrate a mobile light communication system using a time-division multiplexing(TDM) scheme to achieve bidirectional data transmission via the same optical channel.Two identical blue MQW diodes are defined by software as a transmitter or a receiver. To address the light alignment issue, an image identification module integrated with a gimbal stabilizer is used to automatically detect the locations of moving targets;thus, underwater audio communication is realized via a mobile blue-light TDM communication mode. This approach not only uses a single link but also integrates mobile nodes in a practical network.
基金funded by the Italian Ministry of Education,University and Research(MIUR)in the frame of the Departments of Excellence Initiative 2018-2027 attributed to DICAM of the University of Trento(grant L.232/2016)in the frame of the PRIN 2017 project Innovative numerical methods for evolutionary partial differential equations and applications,the PRIN 2022 project High order structure-preserving semi-implicit schemes for hyperbolic equations.D.is member of INdAM GNCS and was also co-funded by the European Union NextGenerationEU(PNRR,Spoke 7 CN HPC).Views and opinions expressed are however those of the author(s)only and do not necessarily reflect those of the European Union or the European Research Council.Neither the European Union nor the granting authority can be held responsible for them.
文摘This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.
文摘ADER-WAF methods were first introduced by researchers E.F. Toro and V.A. Titarev. The linear stability criterion for the model equation for the ADER-WAF schemes is CCFL≤1, where CCFLdenotes the Courant-Friedrichs-Lewy (CFL) coefficient. Toro and Titarev employed CCFL=0.95for their experiments. Nonetheless, we noted that the experiments conducted in this study with CCFL=0.95produced solutions exhibiting spurious oscillations, particularly in the high-order ADER-WAF schemes. The homogeneous one-dimensional (1D) non-linear Shallow Water Equations (SWEs) are the subject of these experiments, specifically the solution of the Riemann Problem (RP) associated with the SWEs. The investigation was conducted on four test problems to evaluate the ADER-WAF schemes of second, third, fourth, and fifth order of accuracy. Each test problem constitutes a RP characterized by different wave patterns in its solution. This research has two primary objectives. We begin by illustrating the procedure for implementing the ADER-WAF schemes for the SWEs, providing the required relations. Afterward, following comprehensive testing, we present the range for the CFL coefficient for each test that yields solutions with diminished or eliminated spurious oscillations.
基金support of the National Natural Science Foundation of China(No.51675406)the Basic Research Project Group,China(No.514010106-205)。
文摘The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.
基金supported by the National Natural Science Foundation of China(Grant No.11925204).
文摘In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.
