A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue prob...A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).展开更多
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.展开更多
This paper deals with Mckean-Vlasov backward stochastic differential equations with weak monotonicity coefficients.We first establish the existence and uniqueness of solutions to Mckean-Vlasov backward stochastic diff...This paper deals with Mckean-Vlasov backward stochastic differential equations with weak monotonicity coefficients.We first establish the existence and uniqueness of solutions to Mckean-Vlasov backward stochastic differential equations.Then we obtain a comparison theorem in one-dimensional situation.展开更多
The neutron diffusion equation plays a pivotal role in nuclear reactor analysis.Nevertheless,employing the physics-informed neural network(PINN)method for its solution entails certain limitations.Conventional PINN app...The neutron diffusion equation plays a pivotal role in nuclear reactor analysis.Nevertheless,employing the physics-informed neural network(PINN)method for its solution entails certain limitations.Conventional PINN approaches generally utilize a fully connected network(FCN)architecture that is susceptible to overfitting,training instability,and gradient vanishing as the network depth increases.These challenges result in accuracy bottlenecks in the solution.In response to these issues,the residual-based resample physics-informed neural network(R2-PINN)is proposed.It is an improved PINN architecture that replaces the FCN with a convolutional neural network with a shortcut(S-CNN).It incorporates skip connections to facilitate gradient propagation between network layers.Additionally,the incorporation of the residual adaptive resampling(RAR)mechanism dynamically increases the number of sampling points.This,in turn,enhances the spatial representation capabilities and overall predictive accuracy of the model.The experimental results illustrate that our approach significantly improves the convergence capability of the model and achieves high-precision predictions of the physical fields.Compared with conventional FCN-based PINN methods,R 2-PINN effectively overcomes the limitations inherent in current methods.Thus,it provides more accurate and robust solutions for neutron diffusion equations.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,t...In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.展开更多
Mesophase pitch carbon fibers have an ultra-high modulus and thermal conductivity that are unmatched by other carbon fibers,making it irreplaceable in many fields.However,due to the high temperature dependence of the ...Mesophase pitch carbon fibers have an ultra-high modulus and thermal conductivity that are unmatched by other carbon fibers,making it irreplaceable in many fields.However,due to the high temperature dependence of the viscosity of the melted pitch and the poor mechanical properties of pitch fibers,it is difficult to reduce the fiber diameter when using continuous spinning.We used the Mathworks Matlab software to optimize the mesophase pitch melt spinning model and to simulate the effects of spinning temperature,mass flow rate,winder speed,and quenching air temperature near the spinneret on the maximum shear rate during drawing.Simulation results demonstrate that applying gradient cooling to the melt upon exiting the spinneret significantly reduces the maximum shear rate and extends the drawing zone,thereby promoting the spinning stability and helping reduce the fiber diameter.In the experiment,instead of quenching in air,we applied gradient cooling to the melt,whose temperature decreased according to the equation Ta=298+278exp(−11.4z),where Ta is the final air temperature in Kelvin,and z is the distance from the spinneret in meters.It was found the gradient cooling greatly improved the draw-down ratio,reducing the average diameter of the pitch fibers from 20.8 to 13.1μm,along with improved process stability.The experimental results are in excellent agreement with the predictions.At the same time,the tensile strength of the 1150℃ carbonized fibers increased from 0.6 to 1.1 GPa.Although the degree of orientation of the fibers decreased slightly,the tight bonding between microcrystals,the suppression of splitting,and the smaller diameter improved the mechanical properties of carbon fibers.This study provides an effective method for reducing the fiber diameter while improving continuity.展开更多
In this work,we investigate the existence and asymptotic stability in mean square of mild solutions for non-linear impulsive neutral stochastic evolution equations with infinite delays in distribution in a real separa...In this work,we investigate the existence and asymptotic stability in mean square of mild solutions for non-linear impulsive neutral stochastic evolution equations with infinite delays in distribution in a real separable Hilbert space.By using the Banach fixed point principle,some sufficient conditions are derived to ensure the asymptotic stability of mild solutions.Moreover,we investigate the Hyers-Ulam stability for such stochastic system.Finally,an illustrative example is given to demonstrate the effectiveness of the obtained results.展开更多
This paper studies the Smoluchowski–Kramers approximation for a discrete-time dynamical system modeled as the motion of a particle in a force field.We show that the approximation holds for the drift-implicit Euler–M...This paper studies the Smoluchowski–Kramers approximation for a discrete-time dynamical system modeled as the motion of a particle in a force field.We show that the approximation holds for the drift-implicit Euler–Maruyama discretization and derive its convergence rate.In particular,the solution of the discretized system converges to the solution of the first-order limit equation in the mean-square sense,and this convergence is independent of the order in which the mass parameterμand the step size h tend to zero.展开更多
Fractional differential equations(FDEs)provide a powerful tool for modeling systems with memory and non-local effects,but understanding their underlying structure remains a significant challenge.While numerous numeric...Fractional differential equations(FDEs)provide a powerful tool for modeling systems with memory and non-local effects,but understanding their underlying structure remains a significant challenge.While numerous numerical and semi-analytical methods exist to find solutions,new approaches are needed to analyze the intrinsic properties of the FDEs themselves.This paper introduces a novel computational framework for the structural analysis of FDEs involving iterated Caputo derivatives.The methodology is based on a transformation that recasts the original FDE into an equivalent higher-order form,represented as the sum of a closed-form,integer-order component G(y)and a residual fractional power seriesΨ(x).This transformed FDE is subsequently reduced to a first-order ordinary differential equation(ODE).The primary novelty of the proposed methodology lies in treating the structure of the integer-order component G(y)not as fixed,but as a parameterizable polynomial whose coefficients can be determined via global optimization.Using particle swarm optimization,the framework identifies an optimal ODE architecture by minimizing a dual objective that balances solution accuracy against a high-fidelity reference and the magnitude of the truncated residual series.The effectiveness of the approach is demonstrated on both a linear FDE and a nonlinear fractional Riccati equation.Results demonstrate that the framework successfully identifies an optimal,low-degree polynomial ODE architecture that is not necessarily identical to the forcing function of the original FDE.This work provides a new tool for analyzing the underlying structure of FDEs and gaining deeper insights into the interplay between local and non-local dynamics in fractional systems.展开更多
The main purpose of this paper is to try to find all entire solutions of the Fermat type difference-differential equation[p1(z)f(z+c)]^(2)+[p2(z)f(z)+p3(z)f′(z)]^(2)=p(z);or[p1(z)f(z)]^(2)+[p2(z)f′(z)+p3(z)f(z+c)]^(...The main purpose of this paper is to try to find all entire solutions of the Fermat type difference-differential equation[p1(z)f(z+c)]^(2)+[p2(z)f(z)+p3(z)f′(z)]^(2)=p(z);or[p1(z)f(z)]^(2)+[p2(z)f′(z)+p3(z)f(z+c)]^(2)=p(z)or[p1(z)f′(z)]^(2)+[p2(z)f(z+c)+p3(z)f(z)]^(2)=p(z);where c is a nonzero complex number,p1;p2 and p3 are polynomials in C satisfying p1p3■0;and p is a nonzero irreducible polynomial in C.展开更多
Biomass is among the most important state variables used to characterize ecosystems. Estimation of tree biomass involves the development of species-specific “allometric equations” that describe the relationship betw...Biomass is among the most important state variables used to characterize ecosystems. Estimation of tree biomass involves the development of species-specific “allometric equations” that describe the relationship between tree biomass and tree diameter and/or height. While many allometric equations were developed for northern hemisphere and tropical species, rarely have they been developed for trees in arid ecosystems, limiting, amongst other things, our ability to estimate carbon stocks in arid regions. Acacia raddiana and A. tortilis are major components of savannas and arid regions in the Middle East and Africa, where they are considered keystone species. Using the opportunity that trees were being uprooted for land development, we measured height (H), north-south (C1) and east-west (C2) canopy diameters, stem diameter at 1.3 meters of the largest stem (D1.3 or DBH), and aboveground fresh and dry weight (FW and DW, respectively) of nine trees (n = 9) from each species. For A. tortilis only, we recorded the number of trunks, and measured the diameter of the largest trunk at ground level (D0). While the average crown (canopy) size (C1 + C2) was very similar among the two species, Acacia raddiana trees were found to be significantly taller than their Acacia tortilis counterparts. Results show that in the arid Arava (southern Israel), an average adult acacia tree has ~200 kg of aboveground dry biomass and that a typical healthy acacia ecosystem in this region, may include ~41 tons of tree biomass per km2. The coefficients of DBH (tree diameter at breast height) to biomass and wood volume, could be used by researchers studying acacia trees throughout the Middle East and Africa, enabling them to estimate biomass of acacia trees and to evaluate their importance for carbon stocks in their arid regions. Highlights: 1) Estimations of tree biomass in arid regions are rare. 2) Biomass allometric equations were developed for A. raddiana and A. tortilis trees. 3) Equations contribute to the estimation of carbon stocks in arid regions.展开更多
In this paper,we study the weighted higher order semilinear equation in an exterior domain(-△)^(m)u=|x|^(α)g(u)in R^(N)\B_(R_(0)),where N≥1,m≥2 are integers,α>-2m,g is a continuous and nondecreasing function i...In this paper,we study the weighted higher order semilinear equation in an exterior domain(-△)^(m)u=|x|^(α)g(u)in R^(N)\B_(R_(0)),where N≥1,m≥2 are integers,α>-2m,g is a continuous and nondecreasing function in(0,+∞)and positive in(0,+∞),B_(R_(0))is the ball of the radius R0 centered at the origin.We prove that a positive supersolution of the problem which verifies(-△)_(i)u>0 in R^(N)\B_(R_(0))(i=0,…,m-1)exists if and only if N>2m and∫_(0)^(δ)g(t)/t^(2(N-m)+α/N-2m)dt<∞,,for someδ>0.We further obtain some existence and nonexistence results for the positive solution to the Dirichlet problem when g(u)=u^(p)with p>1,by using the Pohozaev identity and an embedding lemma of radial Sobolev spaces.展开更多
In this paper,we consider the inhomogeneous pressureless Euler equations.First,we present a class of self-similar analytical solutions to the 1D Cauchy problem and investigate the large-time behavior of the solutions,...In this paper,we consider the inhomogeneous pressureless Euler equations.First,we present a class of self-similar analytical solutions to the 1D Cauchy problem and investigate the large-time behavior of the solutions,and particularly,we obtain slant kink-wave solutions for the inhomogeneous Burgers(InhB)type equation.Next,we prove the integrability of the InhB equation in the sense of Lax pair.Furthermore,we study the spreading rate of the moving domain occupied by mass for the 1D Cauchy problem with compact support initial density.We find that the expanding domain grows exponentially in time,provided that the solutions exist and smooth at all time.Finally,we extend the corresponding results of the inhomogeneous pressureless Euler equations to the radially symmetric multi-dimensional case.展开更多
The Andean montane forests provide a wide range of ecosystem services like water supply, carbon sequestration, and biodiversity preservation. Restoration of these forests is critical due to their degraded state and th...The Andean montane forests provide a wide range of ecosystem services like water supply, carbon sequestration, and biodiversity preservation. Restoration of these forests is critical due to their degraded state and the need to recover, maintain and enhance the ecosystem services they provide. However, we lack understanding of aboveground biomass (AGB) accumulation in restored Andean montane forests. AGB is a key indicator of ecosystem productivity and provides essential data on vegetation carbon stocks, permitting the assess successfulness of restoration efforts. In 2010 the initiative Más Bosques para Medellín was formulated in Medellín City, tropical Andes, Colombia, aiming to restore the forests located in the surrounding rural areas of the city, with interest in preserving the ecosystems services like water supply. The project established 548 ha of mixed plantations with native species. After 13 years, our study aims to developed in situ allometric equations and to evaluate AGB accumulation to assess restoration performance. We measured, harvested, and weighted 144 individuals from these arrangements to fit a general equation for the project and six specific equations for each one of the six most frequent species. The AGB had a positive correlation with diameter at breast height (D), total height (H) and specific wood density (WD). The best general equation uses D and WD as predictors (R^(2) = 0.928). The specific species equations certainly responded to the functional traits of each species. Using the latest inventory of permanent plots of the project we estimated a mean AGB accumulation of 41.91 ± 30.34 Mg ha^(–1) and a total accumulation of 22,996.05 Mg of AGB for the 548 ha. We compared these results with studies developed for natural forest in the region and other land covers. We found contrast behaviors in the AGB accumulation across our study zones. The developed equations have broad applicability across the Andes region, offering valuable insights for similar restoration initiatives. Furthermore, will facilitate the assessment of current restoration efforts and inform scientifically based decisions for future mixed plantation practices.展开更多
The stability of solutions of Herglotz-type equations for non-autonomous non-conservative systems is studied by means of generalized gradient method.Firstly,Herglotz-type equations for non-conservative systems are giv...The stability of solutions of Herglotz-type equations for non-autonomous non-conservative systems is studied by means of generalized gradient method.Firstly,Herglotz-type equations for non-conservative systems are given and expressed as contravariant algebraic form.Secondly,two classes of generalized gradient systems are introduced.Thirdly,the conditions for the transformation of Herglotz-type equations into generalized gradient systems are given,and the solutions of Herglotz-type equations and the stability of the solutions are analyzed.Finally,for each case discussed in this paper,the calculation process is demonstrated in detail to show that the method is effective.展开更多
In recent years,variable-order fractional partial differential equations have attracted growing interest due to their enhanced ability tomodel complex physical phenomena withmemory and spatial heterogeneity.However,ex...In recent years,variable-order fractional partial differential equations have attracted growing interest due to their enhanced ability tomodel complex physical phenomena withmemory and spatial heterogeneity.However,existing numerical methods often struggle with the computational challenges posed by such equations,especially in nonlinear,multi-term formulations.This study introduces two hybrid numerical methods—the Linear-Sine and Cosine(L1-CAS)and fast-CAS schemes—for solving linear and nonlinear multi-term Caputo variable-order(CVO)fractional partial differential equations.These methods combine CAS wavelet-based spatial discretization with L1 and fast algorithms in the time domain.A key feature of the approach is its ability to efficiently handle fully coupled spacetime variable-order derivatives and nonlinearities through a second-order interpolation technique.In addition,we derive CAS wavelet operational matrices for variable-order integration and for boundary value problems,forming the foundation of the spatial discretization.Numerical experiments confirm the accuracy,stability,and computational efficiency of the proposed methods.展开更多
Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremend...Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremendous time due to the extremely large size encountered in most real-world engineering applications.So,practical solvers for systems of linear and nonlinear equations based on multi graphic process units（GPUs）are proposed in order to accelerate the solving process.In the linear and nonlinear solvers,the preconditioned bi-conjugate gradient stable（PBi-CGstab）method and the Inexact Newton method are used to achieve the fast and stable convergence behavior.Multi-GPUs are utilized to obtain more data storage that large size problems need.展开更多
In this paper,the physics informed neural network(PINN)deep learning method is applied to solve two-dimensional nonlocal equations,including the partial reverse space y-nonlocal Mel'nikov equation,the partial reve...In this paper,the physics informed neural network(PINN)deep learning method is applied to solve two-dimensional nonlocal equations,including the partial reverse space y-nonlocal Mel'nikov equation,the partial reverse space-time nonlocal Mel'nikov equation and the nonlocal twodimensional nonlinear Schr?dinger(NLS)equation.By the PINN method,we successfully derive a data-driven two soliton solution,lump solution and rogue wave solution.Numerical simulation results indicate that the error range between the data-driven solution and the exact solution is relatively small,which verifies the effectiveness of the PINN deep learning method for solving high dimensional nonlocal equations.Moreover,the parameter discovery of the partial reverse space-time nonlocal Mel'nikov equation is analysed in terms of its soliton solution for the first time.展开更多
基金supported by the National Natural Science Foundation of China(12401482)the second author was supported by the National Natural Science Foundation of China(12371371,12261160361,11971366)supported by the Open Research Fund of Hubei Key Laboratory of Computational Science,Wuhan University.
文摘A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).
文摘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.
基金Supported by the National Natural Science Foundation of China(12001074)the Research Innovation Program of Graduate Students in Hunan Province(CX20220258)+1 种基金the Research Innovation Program of Graduate Students of Central South University(1053320214147)the Key Scientific Research Project of Higher Education Institutions in Henan Province(25B110025)。
文摘This paper deals with Mckean-Vlasov backward stochastic differential equations with weak monotonicity coefficients.We first establish the existence and uniqueness of solutions to Mckean-Vlasov backward stochastic differential equations.Then we obtain a comparison theorem in one-dimensional situation.
基金supported by the Science and Technology on Reactor System Design Technology Laboratory(No.LRSDT12023108)supported in part by the Chongqing Postdoctoral Science Foundation(No.cstc2021jcyj-bsh0252)+2 种基金the National Natural Science Foundation of China(No.12005030)Sichuan Province to unveil the list of marshal industry common technology research projects(No.23jBGOV0001)Special Program for Stabilizing Support to Basic Research of National Basic Research Institutes(No.WDZC-2023-05-03-05).
文摘The neutron diffusion equation plays a pivotal role in nuclear reactor analysis.Nevertheless,employing the physics-informed neural network(PINN)method for its solution entails certain limitations.Conventional PINN approaches generally utilize a fully connected network(FCN)architecture that is susceptible to overfitting,training instability,and gradient vanishing as the network depth increases.These challenges result in accuracy bottlenecks in the solution.In response to these issues,the residual-based resample physics-informed neural network(R2-PINN)is proposed.It is an improved PINN architecture that replaces the FCN with a convolutional neural network with a shortcut(S-CNN).It incorporates skip connections to facilitate gradient propagation between network layers.Additionally,the incorporation of the residual adaptive resampling(RAR)mechanism dynamically increases the number of sampling points.This,in turn,enhances the spatial representation capabilities and overall predictive accuracy of the model.The experimental results illustrate that our approach significantly improves the convergence capability of the model and achieves high-precision predictions of the physical fields.Compared with conventional FCN-based PINN methods,R 2-PINN effectively overcomes the limitations inherent in current methods.Thus,it provides more accurate and robust solutions for neutron diffusion equations.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金Supported by Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)+2 种基金Basic Research Plan of Shanxi Province(202203021211129)Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)。
文摘In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.
文摘Mesophase pitch carbon fibers have an ultra-high modulus and thermal conductivity that are unmatched by other carbon fibers,making it irreplaceable in many fields.However,due to the high temperature dependence of the viscosity of the melted pitch and the poor mechanical properties of pitch fibers,it is difficult to reduce the fiber diameter when using continuous spinning.We used the Mathworks Matlab software to optimize the mesophase pitch melt spinning model and to simulate the effects of spinning temperature,mass flow rate,winder speed,and quenching air temperature near the spinneret on the maximum shear rate during drawing.Simulation results demonstrate that applying gradient cooling to the melt upon exiting the spinneret significantly reduces the maximum shear rate and extends the drawing zone,thereby promoting the spinning stability and helping reduce the fiber diameter.In the experiment,instead of quenching in air,we applied gradient cooling to the melt,whose temperature decreased according to the equation Ta=298+278exp(−11.4z),where Ta is the final air temperature in Kelvin,and z is the distance from the spinneret in meters.It was found the gradient cooling greatly improved the draw-down ratio,reducing the average diameter of the pitch fibers from 20.8 to 13.1μm,along with improved process stability.The experimental results are in excellent agreement with the predictions.At the same time,the tensile strength of the 1150℃ carbonized fibers increased from 0.6 to 1.1 GPa.Although the degree of orientation of the fibers decreased slightly,the tight bonding between microcrystals,the suppression of splitting,and the smaller diameter improved the mechanical properties of carbon fibers.This study provides an effective method for reducing the fiber diameter while improving continuity.
基金supported by the Excellent Talents Support Program for Colleges and Universities in Anhui Province(Grant No.gxyqZD2020003)。
文摘In this work,we investigate the existence and asymptotic stability in mean square of mild solutions for non-linear impulsive neutral stochastic evolution equations with infinite delays in distribution in a real separable Hilbert space.By using the Banach fixed point principle,some sufficient conditions are derived to ensure the asymptotic stability of mild solutions.Moreover,we investigate the Hyers-Ulam stability for such stochastic system.Finally,an illustrative example is given to demonstrate the effectiveness of the obtained results.
基金supported by the PhD Research Startup Foundation of Hubei University of Economics(Grand No.XJ23BS42).
文摘This paper studies the Smoluchowski–Kramers approximation for a discrete-time dynamical system modeled as the motion of a particle in a force field.We show that the approximation holds for the drift-implicit Euler–Maruyama discretization and derive its convergence rate.In particular,the solution of the discretized system converges to the solution of the first-order limit equation in the mean-square sense,and this convergence is independent of the order in which the mass parameterμand the step size h tend to zero.
基金Research Council of Lithuania(LMTLT),agreement No.S-PD-24-120Research Council of Lithuania(LMTLT),agreement No.S-PD-24-120funded by the Research Council of Lithuania.
文摘Fractional differential equations(FDEs)provide a powerful tool for modeling systems with memory and non-local effects,but understanding their underlying structure remains a significant challenge.While numerous numerical and semi-analytical methods exist to find solutions,new approaches are needed to analyze the intrinsic properties of the FDEs themselves.This paper introduces a novel computational framework for the structural analysis of FDEs involving iterated Caputo derivatives.The methodology is based on a transformation that recasts the original FDE into an equivalent higher-order form,represented as the sum of a closed-form,integer-order component G(y)and a residual fractional power seriesΨ(x).This transformed FDE is subsequently reduced to a first-order ordinary differential equation(ODE).The primary novelty of the proposed methodology lies in treating the structure of the integer-order component G(y)not as fixed,but as a parameterizable polynomial whose coefficients can be determined via global optimization.Using particle swarm optimization,the framework identifies an optimal ODE architecture by minimizing a dual objective that balances solution accuracy against a high-fidelity reference and the magnitude of the truncated residual series.The effectiveness of the approach is demonstrated on both a linear FDE and a nonlinear fractional Riccati equation.Results demonstrate that the framework successfully identifies an optimal,low-degree polynomial ODE architecture that is not necessarily identical to the forcing function of the original FDE.This work provides a new tool for analyzing the underlying structure of FDEs and gaining deeper insights into the interplay between local and non-local dynamics in fractional systems.
基金Supported by the National Natural Science Foundation of China(11871260,11761050)the Jiangxi Natural Science Foundation(#20232ACB201005)+1 种基金the Shandong Natural Science Foundation(#ZR2024MA024)Doctoral Startup Fund of Jiangxi Science and Technology Normal University(#2021BSQD30).
文摘The main purpose of this paper is to try to find all entire solutions of the Fermat type difference-differential equation[p1(z)f(z+c)]^(2)+[p2(z)f(z)+p3(z)f′(z)]^(2)=p(z);or[p1(z)f(z)]^(2)+[p2(z)f′(z)+p3(z)f(z+c)]^(2)=p(z)or[p1(z)f′(z)]^(2)+[p2(z)f(z+c)+p3(z)f(z)]^(2)=p(z);where c is a nonzero complex number,p1;p2 and p3 are polynomials in C satisfying p1p3■0;and p is a nonzero irreducible polynomial in C.
文摘Biomass is among the most important state variables used to characterize ecosystems. Estimation of tree biomass involves the development of species-specific “allometric equations” that describe the relationship between tree biomass and tree diameter and/or height. While many allometric equations were developed for northern hemisphere and tropical species, rarely have they been developed for trees in arid ecosystems, limiting, amongst other things, our ability to estimate carbon stocks in arid regions. Acacia raddiana and A. tortilis are major components of savannas and arid regions in the Middle East and Africa, where they are considered keystone species. Using the opportunity that trees were being uprooted for land development, we measured height (H), north-south (C1) and east-west (C2) canopy diameters, stem diameter at 1.3 meters of the largest stem (D1.3 or DBH), and aboveground fresh and dry weight (FW and DW, respectively) of nine trees (n = 9) from each species. For A. tortilis only, we recorded the number of trunks, and measured the diameter of the largest trunk at ground level (D0). While the average crown (canopy) size (C1 + C2) was very similar among the two species, Acacia raddiana trees were found to be significantly taller than their Acacia tortilis counterparts. Results show that in the arid Arava (southern Israel), an average adult acacia tree has ~200 kg of aboveground dry biomass and that a typical healthy acacia ecosystem in this region, may include ~41 tons of tree biomass per km2. The coefficients of DBH (tree diameter at breast height) to biomass and wood volume, could be used by researchers studying acacia trees throughout the Middle East and Africa, enabling them to estimate biomass of acacia trees and to evaluate their importance for carbon stocks in their arid regions. Highlights: 1) Estimations of tree biomass in arid regions are rare. 2) Biomass allometric equations were developed for A. raddiana and A. tortilis trees. 3) Equations contribute to the estimation of carbon stocks in arid regions.
文摘In this paper,we study the weighted higher order semilinear equation in an exterior domain(-△)^(m)u=|x|^(α)g(u)in R^(N)\B_(R_(0)),where N≥1,m≥2 are integers,α>-2m,g is a continuous and nondecreasing function in(0,+∞)and positive in(0,+∞),B_(R_(0))is the ball of the radius R0 centered at the origin.We prove that a positive supersolution of the problem which verifies(-△)_(i)u>0 in R^(N)\B_(R_(0))(i=0,…,m-1)exists if and only if N>2m and∫_(0)^(δ)g(t)/t^(2(N-m)+α/N-2m)dt<∞,,for someδ>0.We further obtain some existence and nonexistence results for the positive solution to the Dirichlet problem when g(u)=u^(p)with p>1,by using the Pohozaev identity and an embedding lemma of radial Sobolev spaces.
基金Supported by the Henan Natural Science Foundation(242300421397)the Basic Research Projects of Key Scientific Research Projects Plan in Henan Higher Education Institutions(25ZX013)+2 种基金the Scientific Research Team Plan of Zhengzhou University of Aeronautics(23ZHTD01003)the National Natural Science Foundation of China(11971475)the FLASS Internationalization and Exchange Scheme(FLASS/IE−D09/19-20−FLASS).
文摘In this paper,we consider the inhomogeneous pressureless Euler equations.First,we present a class of self-similar analytical solutions to the 1D Cauchy problem and investigate the large-time behavior of the solutions,and particularly,we obtain slant kink-wave solutions for the inhomogeneous Burgers(InhB)type equation.Next,we prove the integrability of the InhB equation in the sense of Lax pair.Furthermore,we study the spreading rate of the moving domain occupied by mass for the 1D Cauchy problem with compact support initial density.We find that the expanding domain grows exponentially in time,provided that the solutions exist and smooth at all time.Finally,we extend the corresponding results of the inhomogeneous pressureless Euler equations to the radially symmetric multi-dimensional case.
基金supported by resources of the Secretary of Environment of the Medellín district,through the agreement N°001-2023 between this entity,the Faculty of Agricultural Sciences of National University of Colombia at Medellín
文摘The Andean montane forests provide a wide range of ecosystem services like water supply, carbon sequestration, and biodiversity preservation. Restoration of these forests is critical due to their degraded state and the need to recover, maintain and enhance the ecosystem services they provide. However, we lack understanding of aboveground biomass (AGB) accumulation in restored Andean montane forests. AGB is a key indicator of ecosystem productivity and provides essential data on vegetation carbon stocks, permitting the assess successfulness of restoration efforts. In 2010 the initiative Más Bosques para Medellín was formulated in Medellín City, tropical Andes, Colombia, aiming to restore the forests located in the surrounding rural areas of the city, with interest in preserving the ecosystems services like water supply. The project established 548 ha of mixed plantations with native species. After 13 years, our study aims to developed in situ allometric equations and to evaluate AGB accumulation to assess restoration performance. We measured, harvested, and weighted 144 individuals from these arrangements to fit a general equation for the project and six specific equations for each one of the six most frequent species. The AGB had a positive correlation with diameter at breast height (D), total height (H) and specific wood density (WD). The best general equation uses D and WD as predictors (R^(2) = 0.928). The specific species equations certainly responded to the functional traits of each species. Using the latest inventory of permanent plots of the project we estimated a mean AGB accumulation of 41.91 ± 30.34 Mg ha^(–1) and a total accumulation of 22,996.05 Mg of AGB for the 548 ha. We compared these results with studies developed for natural forest in the region and other land covers. We found contrast behaviors in the AGB accumulation across our study zones. The developed equations have broad applicability across the Andes region, offering valuable insights for similar restoration initiatives. Furthermore, will facilitate the assessment of current restoration efforts and inform scientifically based decisions for future mixed plantation practices.
基金Supported by the National Natural Science Foundation of China(12272248)。
文摘The stability of solutions of Herglotz-type equations for non-autonomous non-conservative systems is studied by means of generalized gradient method.Firstly,Herglotz-type equations for non-conservative systems are given and expressed as contravariant algebraic form.Secondly,two classes of generalized gradient systems are introduced.Thirdly,the conditions for the transformation of Herglotz-type equations into generalized gradient systems are given,and the solutions of Herglotz-type equations and the stability of the solutions are analyzed.Finally,for each case discussed in this paper,the calculation process is demonstrated in detail to show that the method is effective.
基金supported by the National Research Foundation of Korea(NRF)grant funded by the Korean government(MSIT)(NRF-2021R1A2C1011817)the BK21 Program(Next Generation Education Program for Mathematical Sciences,4299990414089)funded by the Ministry of Education(MOE,Republic of Korea).
文摘In recent years,variable-order fractional partial differential equations have attracted growing interest due to their enhanced ability tomodel complex physical phenomena withmemory and spatial heterogeneity.However,existing numerical methods often struggle with the computational challenges posed by such equations,especially in nonlinear,multi-term formulations.This study introduces two hybrid numerical methods—the Linear-Sine and Cosine(L1-CAS)and fast-CAS schemes—for solving linear and nonlinear multi-term Caputo variable-order(CVO)fractional partial differential equations.These methods combine CAS wavelet-based spatial discretization with L1 and fast algorithms in the time domain.A key feature of the approach is its ability to efficiently handle fully coupled spacetime variable-order derivatives and nonlinearities through a second-order interpolation technique.In addition,we derive CAS wavelet operational matrices for variable-order integration and for boundary value problems,forming the foundation of the spatial discretization.Numerical experiments confirm the accuracy,stability,and computational efficiency of the proposed methods.
文摘Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremendous time due to the extremely large size encountered in most real-world engineering applications.So,practical solvers for systems of linear and nonlinear equations based on multi graphic process units（GPUs）are proposed in order to accelerate the solving process.In the linear and nonlinear solvers,the preconditioned bi-conjugate gradient stable（PBi-CGstab）method and the Inexact Newton method are used to achieve the fast and stable convergence behavior.Multi-GPUs are utilized to obtain more data storage that large size problems need.
文摘In this paper,the physics informed neural network(PINN)deep learning method is applied to solve two-dimensional nonlocal equations,including the partial reverse space y-nonlocal Mel'nikov equation,the partial reverse space-time nonlocal Mel'nikov equation and the nonlocal twodimensional nonlinear Schr?dinger(NLS)equation.By the PINN method,we successfully derive a data-driven two soliton solution,lump solution and rogue wave solution.Numerical simulation results indicate that the error range between the data-driven solution and the exact solution is relatively small,which verifies the effectiveness of the PINN deep learning method for solving high dimensional nonlocal equations.Moreover,the parameter discovery of the partial reverse space-time nonlocal Mel'nikov equation is analysed in terms of its soliton solution for the first time.
