In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependen...In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependence. Test problems including the single soliton wave motion, interaction of two solitary waves and interaction of three solitary waves will use to validate the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Finally, a Maxwellian initial condition pulse is then studied. The L<sub>2</sub> and L<sub>∞</sub> error norms are computed to study the accuracy and the simplicity of the presented method.展开更多
This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and ...This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and transition of 43×43 computational grids,a coordinate transformation is put up from physical panel to computational panel. Several zero turbulent models are computed comparatively. The results are credible when comparing with the previous methods.展开更多
Solar-powered aircraft have attracted great attention owing to their potential for longendurance flight and wide application prospects.Due to the particularity of energy system,flight strategy optimization is a signif...Solar-powered aircraft have attracted great attention owing to their potential for longendurance flight and wide application prospects.Due to the particularity of energy system,flight strategy optimization is a significant way to enhance the flight performance for solar-powered aircraft.In this study,a flight strategy optimization model for high-altitude long-endurance solar-powered aircraft was proposed.This model consists of three-dimensional kinematic model,aerodynamic model,energy collection model,energy store model and energy loss model.To solve the nonlinear optimal control problem with process constraints and terminal constraints,Gauss pseudo-spectral method was employed to discretize the state equations and constraint equations.Then a typical mission flying from given initial point to given final point within a time interval was considered.Results indicate that proper changes of the attitude angle contribute to increasing the energy gained by photovoltaic cells.Utilization of gravitational potential energy can partly take the role of battery pack.Integrating these two measures,the optimized flight strategy can improve the final state of charge compared with current constant-altitude constant-velocity strategy.The optimized strategy brings more profits on condition of lower sunlight intensity and shorter daytime.展开更多
Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. ...Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.展开更多
With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important me...With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important methods of wave-field simulation. Results of previous studies show that both methods have distinct advantages and disadvantages: Finite difference method has high precision but its dispersion is serious; pseudospectral method considers both computational efficiency and precision but has less precision than finite-difference. The authors consider the complex structural characteristics of the metal ore,furthermore add random media in order to simulate the complex effects produced by metal ore for wave field. First,the study introduced the theories of random media and two forward modelling methods. Second,it compared the simulation results of two methods on fault model. Then the authors established a complex metal ore model,added random media and compared computational efficiency and precision. As a result,it is found that finite difference method is better than pseudo-spectral method in precision and boundary treatment,but the computational efficiency of pseudospectral method is slightly higher than the finite difference method.展开更多
This paper presents a numerical scheme for space fractional diffusion equations (SFDEs) based on pseudo-spectral method. In this approach, using the Guass-Lobatto nodes, the unknown function is approximated by orthogo...This paper presents a numerical scheme for space fractional diffusion equations (SFDEs) based on pseudo-spectral method. In this approach, using the Guass-Lobatto nodes, the unknown function is approximated by orthogonal polynomials or interpolation polynomials. Then, by using pseudo-spectral method, the SFDE is reduced to a system of ordinary differential equations for time variable t. The high order Runge-Kutta scheme can be used to solve the system. So, a high order numerical scheme is derived. Numerical examples illustrate that the results obtained by this method agree well with the analytical solutions.展开更多
In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caput...In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caputo sense. The study is conducted through illustrative example to demonstrate the validity and applicability of the presented method. The results reveal that the proposed method is very effective and simple. Moreover, only a small number of shifted Legendre polynomials are needed to obtain a satisfactory result.展开更多
In this paper, an efficient numerical method is considered for solving the fractional wave equation (FWE). The fractional derivative is described in the Caputo sense. The method is based on Laguerre approximations. Th...In this paper, an efficient numerical method is considered for solving the fractional wave equation (FWE). The fractional derivative is described in the Caputo sense. The method is based on Laguerre approximations. The properties of Laguerre polynomials are utilized to reduce FWE to a system of ordinary differential equations, which is solved by the finite difference method. An approximate formula of the fractional derivative is given. Special attention is given to study the convergence analysis and estimate an error upper bound of the presented formula. Numerical solutions of FWE are given and the results are compared with the exact solution.展开更多
Abstract Recently,the numerical methods for long-time dynamics of PDEs with weak nonlinearity have received more and more attention.For the nonlinear Schrödinger equation(NLS)with wave operator(NLSW)and weak nonl...Abstract Recently,the numerical methods for long-time dynamics of PDEs with weak nonlinearity have received more and more attention.For the nonlinear Schrödinger equation(NLS)with wave operator(NLSW)and weak nonlinearity controlled by a small valueε∈(0,1],an exponential wave integrator Fourier pseudo-spectral(EWIFP)discretization has been developed(Guo et al.,2021)and proved to be uniformly accurate aboutεup to the time atΟ(1/ε^(2))However,the EWIFP method is not time symmetric and can not preserve the discrete energy.As we know,the time symmetry and energy-preservation are the important structural features of the true solution and we hope that this structure can be inherited along the numerical solution.In this work,we propose a time symmetric and energy-preserving exponential wave integrator Fourier pseudo-spectral(SEPEWIFP)method for the NLSW with periodic boundary conditions.Through rigorous error analysis,we establish uniform error bounds of the numerical solution atΟ(h^(mo)+ε^(2-βτ2))up to the time atΟ(1/ε^(β))forβ∈[0,2]where h andτare the mesh size and time step,respectively,and m0 depends on the regularity conditions.The tools for error analysis mainly include cut-off technique and the standard energy method.We also extend the results on error bounds,energy-preservation and time symmetry to the oscillatory NLSW with wavelength atΟ(1/ε^(2))in time which is equivalent to the NLSW with weak nonlinearity.Numerical experiments confirm that the theoretical results in this paper are correct.Our method is novel because that to the best of our knowledge there has not been any energy-preserving exponential wave integrator method for the NLSW.展开更多
This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the n...This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the nonlinear term is mere cubic.The new framework of convergence analysis consists of two steps.In the first step,by truncating the nonlinear term into a global Lipschitz function,an alternative numerical method is proposed and proved in a rigorous way to be convergent in the discrete L2 norm;followed in the second step,the maximum bound of the numerical solution of the alternative numerical method is obtained by using a lifting technique,as implies that the two numerical methods are the same one.Under our framework of convergence analysis,with neither any restriction on the grid ratio nor any requirement of the small initial value,we establish the error estimate of the proposed conservative Fourier pseudo-spectral method,while previous work requires the certain restriction for the focusing case.The error bound is proved to be of O(h^(r)+t^(2))with grid size h and time step t.In fact,the framework can be used to prove the unconditional convergence of many other Fourier pseudo-spectral methods for solving the nonlinear Schr¨odinger-type equations.Numerical results are conducted to indicate the accuracy and efficiency of the proposed method,and investigate the effect of the nonlinear term and initial data on the blow-up solution.展开更多
In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
In this paper,we propose a novel noncausal control framework to address the energy maximization problem of wave energy converters(WECs)subject to constraints.The energy maximization problem of WECs is a constrained op...In this paper,we propose a novel noncausal control framework to address the energy maximization problem of wave energy converters(WECs)subject to constraints.The energy maximization problem of WECs is a constrained optimal control problem.The proposed control framework converts this problem into a reference trajectory tracking problem through the Fourier pseudo-spectral method(FPSM)and utilizes the online tracking adaptive dynamic programming(OTADP)algorithm to realize real-time trajectory tracking for practical use in the ocean environment.Using the wave prediction technique,the optimal trajectory is generated online through a receding horizon(RH)implementation.A critic neural network(NN)is applied to approximate the optimal cost value function and calculate the error-tracking control by solving the associated Hamilton-Jacobi-Bellman(HJB)equation.The proposed WEC control framework improves computational efficiency and makes the online control feasible in practice.Simulation results show the effects of the receding horizon implementation of FPSM with different window lengths and window functions,while verifying the performances of tracking control and energy absorption of WECs in two different sea conditions.展开更多
In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical m...In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical model,and its validity was verified using a simple impact test.A crushable discrete element method(DEM)framework is built based on the previously established theoretical model.The tensile strength,which considers the fractal theory,size effect,and Weibull variation,was assigned to each generated particle.The assigned strength is then used for crush detection by comparing it with its maximum tensile stress.Mass conservation is ensured by inserting a series of sub-particles whose total mass was equal to the quality loss.Based on the crushable DEM framework,a numerical simulation of the crushing behavior of a pebble bed with hollow cylindrical geometry under a uniaxial compression test was performed.The results of this investigation showed that the particle withstands the external load by contact and sliding at the beginning of the compression process,and the results confirmed that crushing can be considered an important method of resisting the increasing external load.A relatively regular particle arrangement aids in resisting the load and reduces the occurrence of particle crushing.However,a limit exists to the promotion of resistance.When the strain increases beyond this limit,the distribution of the crushing position tends to be isotropic over the entire pebble bed.The theoretical model and crushable DEM framework provide a new method for exploring the pebble bed in a fusion reactor,considering particle crushing.展开更多
Effective partitioning is crucial for enabling parallel restoration of power systems after blackouts.This paper proposes a novel partitioning method based on deep reinforcement learning.First,the partitioning decision...Effective partitioning is crucial for enabling parallel restoration of power systems after blackouts.This paper proposes a novel partitioning method based on deep reinforcement learning.First,the partitioning decision process is formulated as a Markov decision process(MDP)model to maximize the modularity.Corresponding key partitioning constraints on parallel restoration are considered.Second,based on the partitioning objective and constraints,the reward function of the partitioning MDP model is set by adopting a relative deviation normalization scheme to reduce mutual interference between the reward and penalty in the reward function.The soft bonus scaling mechanism is introduced to mitigate overestimation caused by abrupt jumps in the reward.Then,the deep Q network method is applied to solve the partitioning MDP model and generate partitioning schemes.Two experience replay buffers are employed to speed up the training process of the method.Finally,case studies on the IEEE 39-bus test system demonstrate that the proposed method can generate a high-modularity partitioning result that meets all key partitioning constraints,thereby improving the parallelism and reliability of the restoration process.Moreover,simulation results demonstrate that an appropriate discount factor is crucial for ensuring both the convergence speed and the stability of the partitioning training.展开更多
The application of nitrogen fertilizers in agricultural fields can lead to the release of nitrogen-containing gases(NCGs),such as NO_(x),NH_(3) and N_(2)O,which can significantly impact regional atmospheric environmen...The application of nitrogen fertilizers in agricultural fields can lead to the release of nitrogen-containing gases(NCGs),such as NO_(x),NH_(3) and N_(2)O,which can significantly impact regional atmospheric environment and con-tribute to global climate change.However,there remain considerable research gaps in the accurate measurement of NCGs emissions from agricultural fields,hindering the development of effective emission reduction strategies.We improved an open-top dynamic chambers(OTDCs)system and evaluated the performance by comparing the measured and given fluxes of the NCGs.The results showed that the measured fluxes of NO,N_(2)O and NH_(3)were 1%,2%and 7%lower than the given fluxes,respectively.For the determination of NH_(3) concentration,we employed a stripping coil-ion chromatograph(SC-IC)analytical technique,which demonstrated an absorption efficiency for atmospheric NH_(3) exceeding 96.1%across sampling durations of 6 to 60 min.In the summer maize season,we utilized the OTDCs system to measure the exchange fluxes of NO,NH_(3),and N_(2)O from the soil in the North China Plain.Substantial emissions of NO,NH_(3) and N_(2)O were recorded following fertilization,with peaks of 107,309,1239 ng N/(m^(2)·s),respectively.Notably,significant NCGs emissions were observed following sus-tained heavy rainfall one month after fertilization,particularly with NH_(3) peak being 4.5 times higher than that observed immediately after fertilization.Our results demonstrate that the OTDCs system accurately reflects the emission characteristics of soil NCGs and meets the requirements for long-term and continuous flux observation.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
At present,there is currently a lack of unified standard methods for the determination of antimony content in groundwater in China.The precision and trueness of related detection technologies have not yet been systema...At present,there is currently a lack of unified standard methods for the determination of antimony content in groundwater in China.The precision and trueness of related detection technologies have not yet been systematically and quantitatively evaluated,which limits the effective implementation of environmental monitoring.In response to this key technical gap,this study aimed to establish a standardized method for determining antimony in groundwater using Hydride Generation–Atomic Fluorescence Spectrometry(HG-AFS).Ten laboratories participated in inter-laboratory collaborative tests,and the statistical analysis of the test data was carried out in strict accordance with the technical specifications of GB/T 6379.2—2004 and GB/T 6379.4—2006.The consistency and outliers of the data were tested by Mandel's h and k statistics,the Grubbs test and the Cochran test,and the outliers were removed to optimize the data,thereby significantly improving the reliability and accuracy.Based on the optimized data,parameters such as the repeatability limit(r),reproducibility limit(R),and method bias value(δ)were determined,and the trueness of the method was statistically evaluated.At the same time,precision-function relationships were established,and all results met the requirements.The results show that the lower the antimony content,the lower the repeatability limit(r)and reproducibility limit(R),indicating that the measurement error mainly originates from the detection limit of the method and instrument sensitivity.Therefore,improving the instrument sensitivity and reducing the detection limit are the keys to controlling the analytical error and improving precision.This study provides reliable data support and a solid technical foundation for the establishment and evaluation of standardized methods for the determination of antimony content in groundwater.展开更多
Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method ge...Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method generated a specific trajectory for the UAV to effectively induce the proportional navigation missile to successfully intercept the obstacle,thereby accomplishing the evasive maneuver.The evasive maneuver was divided into two distinct stages,namely the collision-inducing phase and the fast departure phase.The obstacle potential field-based target selection algorithm was employed to identify the most appropriate target obstacle,while the induced trajectory was determined through a combination of receding horizon optimization and the hp-adaptive pseudo-spectral method.Simulation experiments were carried out under three different types of obstacle environments and one multiobstacle environment,and the simulation results show that the method proposed in this paper greatly improves the success rate of UAV evasive maneuvers,proving the effectiveness of this method.展开更多
To improve the modeling accuracy of radiative transfer,the scattering properties of aerosol particles with irregular shapes and inhomogeneous compositions should be simulated accurately.To this end,a light-scattering ...To improve the modeling accuracy of radiative transfer,the scattering properties of aerosol particles with irregular shapes and inhomogeneous compositions should be simulated accurately.To this end,a light-scattering model for nonspherical particles is established based on the pseudo-spectral time domain(PSTD)technique.In this model,the perfectly matched layer with auxiliary differential equation(ADE-PML),an excellent absorption boundary condition(ABC)in the finite difference time domain generalized for the PSTD,and the weighted total field/scattered field(TF/SF)technique is employed to introduce the incident light into 3 D computational domain.To improve computational efficiency,the model is further parallelized using the Open MP technique.The modeling accuracy of the PSTD scheme is validated against Lorenz–Mie,Aden–Kerker,T-matrix theory and DDA for spheres,inhomogeneous particles and nonspherical particles,and the influence of the spatial resolution and thickness of ADE-PML on the modeling accuracy is discussed as well.Finally,the parallel computational efficiency of the model is also analyzed.The results show that an excellent agreement is achieved between the results of PSTD and well-tested scattering models,where the simulation errors of extinction efficiencies are generally smaller than 1%,indicating the high accuracy of our model.Despite its low spatial resolution,reliable modeling precision can still be achieved by using the PSTD technique,especially for large particles.To suppress the electromagnetic wave reflected by the absorption layers,a six-layer ADE-PML should be set in the computational domain at least.展开更多
文摘In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependence. Test problems including the single soliton wave motion, interaction of two solitary waves and interaction of three solitary waves will use to validate the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Finally, a Maxwellian initial condition pulse is then studied. The L<sub>2</sub> and L<sub>∞</sub> error norms are computed to study the accuracy and the simplicity of the presented method.
文摘This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and transition of 43×43 computational grids,a coordinate transformation is put up from physical panel to computational panel. Several zero turbulent models are computed comparatively. The results are credible when comparing with the previous methods.
文摘Solar-powered aircraft have attracted great attention owing to their potential for longendurance flight and wide application prospects.Due to the particularity of energy system,flight strategy optimization is a significant way to enhance the flight performance for solar-powered aircraft.In this study,a flight strategy optimization model for high-altitude long-endurance solar-powered aircraft was proposed.This model consists of three-dimensional kinematic model,aerodynamic model,energy collection model,energy store model and energy loss model.To solve the nonlinear optimal control problem with process constraints and terminal constraints,Gauss pseudo-spectral method was employed to discretize the state equations and constraint equations.Then a typical mission flying from given initial point to given final point within a time interval was considered.Results indicate that proper changes of the attitude angle contribute to increasing the energy gained by photovoltaic cells.Utilization of gravitational potential energy can partly take the role of battery pack.Integrating these two measures,the optimized flight strategy can improve the final state of charge compared with current constant-altitude constant-velocity strategy.The optimized strategy brings more profits on condition of lower sunlight intensity and shorter daytime.
文摘Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.
基金Supported by the National"863"Project(No.2014AA06A605)
文摘With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important methods of wave-field simulation. Results of previous studies show that both methods have distinct advantages and disadvantages: Finite difference method has high precision but its dispersion is serious; pseudospectral method considers both computational efficiency and precision but has less precision than finite-difference. The authors consider the complex structural characteristics of the metal ore,furthermore add random media in order to simulate the complex effects produced by metal ore for wave field. First,the study introduced the theories of random media and two forward modelling methods. Second,it compared the simulation results of two methods on fault model. Then the authors established a complex metal ore model,added random media and compared computational efficiency and precision. As a result,it is found that finite difference method is better than pseudo-spectral method in precision and boundary treatment,but the computational efficiency of pseudospectral method is slightly higher than the finite difference method.
文摘This paper presents a numerical scheme for space fractional diffusion equations (SFDEs) based on pseudo-spectral method. In this approach, using the Guass-Lobatto nodes, the unknown function is approximated by orthogonal polynomials or interpolation polynomials. Then, by using pseudo-spectral method, the SFDE is reduced to a system of ordinary differential equations for time variable t. The high order Runge-Kutta scheme can be used to solve the system. So, a high order numerical scheme is derived. Numerical examples illustrate that the results obtained by this method agree well with the analytical solutions.
文摘In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caputo sense. The study is conducted through illustrative example to demonstrate the validity and applicability of the presented method. The results reveal that the proposed method is very effective and simple. Moreover, only a small number of shifted Legendre polynomials are needed to obtain a satisfactory result.
文摘In this paper, an efficient numerical method is considered for solving the fractional wave equation (FWE). The fractional derivative is described in the Caputo sense. The method is based on Laguerre approximations. The properties of Laguerre polynomials are utilized to reduce FWE to a system of ordinary differential equations, which is solved by the finite difference method. An approximate formula of the fractional derivative is given. Special attention is given to study the convergence analysis and estimate an error upper bound of the presented formula. Numerical solutions of FWE are given and the results are compared with the exact solution.
基金supported in part by the Natural Science Foundation of Hebei Province(Grant No.A2021205036).
文摘Abstract Recently,the numerical methods for long-time dynamics of PDEs with weak nonlinearity have received more and more attention.For the nonlinear Schrödinger equation(NLS)with wave operator(NLSW)and weak nonlinearity controlled by a small valueε∈(0,1],an exponential wave integrator Fourier pseudo-spectral(EWIFP)discretization has been developed(Guo et al.,2021)and proved to be uniformly accurate aboutεup to the time atΟ(1/ε^(2))However,the EWIFP method is not time symmetric and can not preserve the discrete energy.As we know,the time symmetry and energy-preservation are the important structural features of the true solution and we hope that this structure can be inherited along the numerical solution.In this work,we propose a time symmetric and energy-preserving exponential wave integrator Fourier pseudo-spectral(SEPEWIFP)method for the NLSW with periodic boundary conditions.Through rigorous error analysis,we establish uniform error bounds of the numerical solution atΟ(h^(mo)+ε^(2-βτ2))up to the time atΟ(1/ε^(β))forβ∈[0,2]where h andτare the mesh size and time step,respectively,and m0 depends on the regularity conditions.The tools for error analysis mainly include cut-off technique and the standard energy method.We also extend the results on error bounds,energy-preservation and time symmetry to the oscillatory NLSW with wavelength atΟ(1/ε^(2))in time which is equivalent to the NLSW with weak nonlinearity.Numerical experiments confirm that the theoretical results in this paper are correct.Our method is novel because that to the best of our knowledge there has not been any energy-preserving exponential wave integrator method for the NLSW.
基金Jialing Wang’s work is supported by the National Natural Science Foundation of China(Grant No.11801277)Tingchun Wang’s work is supported by the National Natural Science Foundation of China(Grant No.11571181)+1 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20171454)Qing Lan Project.Yushun Wang’s work is supported by the National Natural Science Foundation of China(Grant Nos.11771213 and 12171245).
文摘This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the nonlinear term is mere cubic.The new framework of convergence analysis consists of two steps.In the first step,by truncating the nonlinear term into a global Lipschitz function,an alternative numerical method is proposed and proved in a rigorous way to be convergent in the discrete L2 norm;followed in the second step,the maximum bound of the numerical solution of the alternative numerical method is obtained by using a lifting technique,as implies that the two numerical methods are the same one.Under our framework of convergence analysis,with neither any restriction on the grid ratio nor any requirement of the small initial value,we establish the error estimate of the proposed conservative Fourier pseudo-spectral method,while previous work requires the certain restriction for the focusing case.The error bound is proved to be of O(h^(r)+t^(2))with grid size h and time step t.In fact,the framework can be used to prove the unconditional convergence of many other Fourier pseudo-spectral methods for solving the nonlinear Schr¨odinger-type equations.Numerical results are conducted to indicate the accuracy and efficiency of the proposed method,and investigate the effect of the nonlinear term and initial data on the blow-up solution.
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金supported by the Key R&D Program of Shandong Province,China(No.2021ZLGX04)the Taishan Industrial Experts Programme(No.tsls20231203)。
文摘In this paper,we propose a novel noncausal control framework to address the energy maximization problem of wave energy converters(WECs)subject to constraints.The energy maximization problem of WECs is a constrained optimal control problem.The proposed control framework converts this problem into a reference trajectory tracking problem through the Fourier pseudo-spectral method(FPSM)and utilizes the online tracking adaptive dynamic programming(OTADP)algorithm to realize real-time trajectory tracking for practical use in the ocean environment.Using the wave prediction technique,the optimal trajectory is generated online through a receding horizon(RH)implementation.A critic neural network(NN)is applied to approximate the optimal cost value function and calculate the error-tracking control by solving the associated Hamilton-Jacobi-Bellman(HJB)equation.The proposed WEC control framework improves computational efficiency and makes the online control feasible in practice.Simulation results show the effects of the receding horizon implementation of FPSM with different window lengths and window functions,while verifying the performances of tracking control and energy absorption of WECs in two different sea conditions.
基金supported by Anhui Provincial Natural Science Foundation(2408085QA030)Natural Science Research Project of Anhui Educational Committee,China(2022AH050825)+3 种基金Medical Special Cultivation Project of Anhui University of Science and Technology(YZ2023H2C008)the Excellent Research and Innovation Team of Anhui Province,China(2022AH010052)the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology,China(2021yjrc51)Collaborative Innovation Program of Hefei Science Center,CAS,China(2019HSC-CIP006).
文摘In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical model,and its validity was verified using a simple impact test.A crushable discrete element method(DEM)framework is built based on the previously established theoretical model.The tensile strength,which considers the fractal theory,size effect,and Weibull variation,was assigned to each generated particle.The assigned strength is then used for crush detection by comparing it with its maximum tensile stress.Mass conservation is ensured by inserting a series of sub-particles whose total mass was equal to the quality loss.Based on the crushable DEM framework,a numerical simulation of the crushing behavior of a pebble bed with hollow cylindrical geometry under a uniaxial compression test was performed.The results of this investigation showed that the particle withstands the external load by contact and sliding at the beginning of the compression process,and the results confirmed that crushing can be considered an important method of resisting the increasing external load.A relatively regular particle arrangement aids in resisting the load and reduces the occurrence of particle crushing.However,a limit exists to the promotion of resistance.When the strain increases beyond this limit,the distribution of the crushing position tends to be isotropic over the entire pebble bed.The theoretical model and crushable DEM framework provide a new method for exploring the pebble bed in a fusion reactor,considering particle crushing.
基金funded by the Beijing Engineering Research Center of Electric Rail Transportation.
文摘Effective partitioning is crucial for enabling parallel restoration of power systems after blackouts.This paper proposes a novel partitioning method based on deep reinforcement learning.First,the partitioning decision process is formulated as a Markov decision process(MDP)model to maximize the modularity.Corresponding key partitioning constraints on parallel restoration are considered.Second,based on the partitioning objective and constraints,the reward function of the partitioning MDP model is set by adopting a relative deviation normalization scheme to reduce mutual interference between the reward and penalty in the reward function.The soft bonus scaling mechanism is introduced to mitigate overestimation caused by abrupt jumps in the reward.Then,the deep Q network method is applied to solve the partitioning MDP model and generate partitioning schemes.Two experience replay buffers are employed to speed up the training process of the method.Finally,case studies on the IEEE 39-bus test system demonstrate that the proposed method can generate a high-modularity partitioning result that meets all key partitioning constraints,thereby improving the parallelism and reliability of the restoration process.Moreover,simulation results demonstrate that an appropriate discount factor is crucial for ensuring both the convergence speed and the stability of the partitioning training.
基金supported by the National Key Research and Develop-ment Program(No.2022YFC3701103)the National Natural Science Foundation of China(Nos.42130714 and 41931287).
文摘The application of nitrogen fertilizers in agricultural fields can lead to the release of nitrogen-containing gases(NCGs),such as NO_(x),NH_(3) and N_(2)O,which can significantly impact regional atmospheric environment and con-tribute to global climate change.However,there remain considerable research gaps in the accurate measurement of NCGs emissions from agricultural fields,hindering the development of effective emission reduction strategies.We improved an open-top dynamic chambers(OTDCs)system and evaluated the performance by comparing the measured and given fluxes of the NCGs.The results showed that the measured fluxes of NO,N_(2)O and NH_(3)were 1%,2%and 7%lower than the given fluxes,respectively.For the determination of NH_(3) concentration,we employed a stripping coil-ion chromatograph(SC-IC)analytical technique,which demonstrated an absorption efficiency for atmospheric NH_(3) exceeding 96.1%across sampling durations of 6 to 60 min.In the summer maize season,we utilized the OTDCs system to measure the exchange fluxes of NO,NH_(3),and N_(2)O from the soil in the North China Plain.Substantial emissions of NO,NH_(3) and N_(2)O were recorded following fertilization,with peaks of 107,309,1239 ng N/(m^(2)·s),respectively.Notably,significant NCGs emissions were observed following sus-tained heavy rainfall one month after fertilization,particularly with NH_(3) peak being 4.5 times higher than that observed immediately after fertilization.Our results demonstrate that the OTDCs system accurately reflects the emission characteristics of soil NCGs and meets the requirements for long-term and continuous flux observation.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金supported by the National Natural Science Foundation of China(Project No.42307555).
文摘At present,there is currently a lack of unified standard methods for the determination of antimony content in groundwater in China.The precision and trueness of related detection technologies have not yet been systematically and quantitatively evaluated,which limits the effective implementation of environmental monitoring.In response to this key technical gap,this study aimed to establish a standardized method for determining antimony in groundwater using Hydride Generation–Atomic Fluorescence Spectrometry(HG-AFS).Ten laboratories participated in inter-laboratory collaborative tests,and the statistical analysis of the test data was carried out in strict accordance with the technical specifications of GB/T 6379.2—2004 and GB/T 6379.4—2006.The consistency and outliers of the data were tested by Mandel's h and k statistics,the Grubbs test and the Cochran test,and the outliers were removed to optimize the data,thereby significantly improving the reliability and accuracy.Based on the optimized data,parameters such as the repeatability limit(r),reproducibility limit(R),and method bias value(δ)were determined,and the trueness of the method was statistically evaluated.At the same time,precision-function relationships were established,and all results met the requirements.The results show that the lower the antimony content,the lower the repeatability limit(r)and reproducibility limit(R),indicating that the measurement error mainly originates from the detection limit of the method and instrument sensitivity.Therefore,improving the instrument sensitivity and reducing the detection limit are the keys to controlling the analytical error and improving precision.This study provides reliable data support and a solid technical foundation for the establishment and evaluation of standardized methods for the determination of antimony content in groundwater.
基金Natural Science Foundation of Heilongjiang Province of China(Grant No.YQ2022F012)the Fundamental Research Funds for the Central Universities(Grant No.HIT.OCEF.2023010)to provide fund for conducting experiments.
文摘Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method generated a specific trajectory for the UAV to effectively induce the proportional navigation missile to successfully intercept the obstacle,thereby accomplishing the evasive maneuver.The evasive maneuver was divided into two distinct stages,namely the collision-inducing phase and the fast departure phase.The obstacle potential field-based target selection algorithm was employed to identify the most appropriate target obstacle,while the induced trajectory was determined through a combination of receding horizon optimization and the hp-adaptive pseudo-spectral method.Simulation experiments were carried out under three different types of obstacle environments and one multiobstacle environment,and the simulation results show that the method proposed in this paper greatly improves the success rate of UAV evasive maneuvers,proving the effectiveness of this method.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.41575025 and 41575024)
文摘To improve the modeling accuracy of radiative transfer,the scattering properties of aerosol particles with irregular shapes and inhomogeneous compositions should be simulated accurately.To this end,a light-scattering model for nonspherical particles is established based on the pseudo-spectral time domain(PSTD)technique.In this model,the perfectly matched layer with auxiliary differential equation(ADE-PML),an excellent absorption boundary condition(ABC)in the finite difference time domain generalized for the PSTD,and the weighted total field/scattered field(TF/SF)technique is employed to introduce the incident light into 3 D computational domain.To improve computational efficiency,the model is further parallelized using the Open MP technique.The modeling accuracy of the PSTD scheme is validated against Lorenz–Mie,Aden–Kerker,T-matrix theory and DDA for spheres,inhomogeneous particles and nonspherical particles,and the influence of the spatial resolution and thickness of ADE-PML on the modeling accuracy is discussed as well.Finally,the parallel computational efficiency of the model is also analyzed.The results show that an excellent agreement is achieved between the results of PSTD and well-tested scattering models,where the simulation errors of extinction efficiencies are generally smaller than 1%,indicating the high accuracy of our model.Despite its low spatial resolution,reliable modeling precision can still be achieved by using the PSTD technique,especially for large particles.To suppress the electromagnetic wave reflected by the absorption layers,a six-layer ADE-PML should be set in the computational domain at least.
