We developed a modified stochastic finite-fault method for estimating strong ground motions.An adjustment to the dynamic corner frequency was introduced,which accounted for the effect of the location of the subfault r...We developed a modified stochastic finite-fault method for estimating strong ground motions.An adjustment to the dynamic corner frequency was introduced,which accounted for the effect of the location of the subfault relative to the hypocenter and rupture propagation direction,to account for the influence of the rupture propagation direction on the subfault dynamic corner frequency.By comparing the peak ground acceleration(PGA),pseudo-absolute response spectra acceleration(PSA,damping ratio of 5%),and duration,the results of the modified and existing methods were compared,demonstrating that our proposed adjustment to the dynamic corner frequency can accurately reflect the rupture directivity effect.We applied our modified method to simulate near-field strong motions within 150 km of the 2008 MW7.9 Wenchuan earthquake rupture.Our modified method performed well over a broad period range,particularly at 0.04-4 s.The total deviations of the stochastic finite-fault method(EXSIM)and the modified EXSIM were 0.1676 and 0.1494,respectively.The modified method can effectively account for the influence of the rupture propagation direction and provide more realistic ground motion estimations for earthquake disaster mitigation.展开更多
Proximal gradient descent and its accelerated version are resultful methods for solving the sum of smooth and non-smooth problems. When the smooth function can be represented as a sum of multiple functions, the stocha...Proximal gradient descent and its accelerated version are resultful methods for solving the sum of smooth and non-smooth problems. When the smooth function can be represented as a sum of multiple functions, the stochastic proximal gradient method performs well. However, research on its accelerated version remains unclear. This paper proposes a proximal stochastic accelerated gradient (PSAG) method to address problems involving a combination of smooth and non-smooth components, where the smooth part corresponds to the average of multiple block sums. Simultaneously, most of convergence analyses hold in expectation. To this end, under some mind conditions, we present an almost sure convergence of unbiased gradient estimation in the non-smooth setting. Moreover, we establish that the minimum of the squared gradient mapping norm arbitrarily converges to zero with probability one.展开更多
Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate b...Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.展开更多
In this paper, we investigate an SIS model with treatment and immigration. Firstly, the two-dimensional model is simplified by using the stochastic averaging method. Then, we derive the local stability of the stochast...In this paper, we investigate an SIS model with treatment and immigration. Firstly, the two-dimensional model is simplified by using the stochastic averaging method. Then, we derive the local stability of the stochastic system by computing the Lyapunov exponent of the linearized system. Further, the global stability of the stochastic model is analyzed based on the singular boundary theory. Moreover, we prove that the model undergoes a Hopf bifurcation and a pitchfork bifurcation. Finally, several numerical examples are provided to illustrate the theoretical results. .展开更多
To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltai...To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.展开更多
A stochastic wheelset model with a nonlinear wheel-rail contact relationship is established to investigate the stochastic stability and stochastic bifurcation of the wheelset system with the consideration of the stoch...A stochastic wheelset model with a nonlinear wheel-rail contact relationship is established to investigate the stochastic stability and stochastic bifurcation of the wheelset system with the consideration of the stochastic parametric excitations of equivalent conicity and suspension stiffness.The wheelset is systematized into a onedimensional(1D)diffusion process by using the stochastic average method,the behavior of the singular boundary is analyzed to determine the hunting stability condition of the wheelset system,and the critical speed of stochastic bifurcation is obtained.The stationary probability density and joint probability density are derived theoretically.Based on the topological structure change of the probability density function,the stochastic Hopf bifurcation form and bifurcation condition of the wheelset system are determined.The effects of stochastic factors on the hunting stability and bifurcation characteristics are analyzed,and the simulation results verify the correctness of the theoretical analysis.The results reveal that the boundary behavior of the diffusion process determines the hunting stability of the stochastic wheelset system,and the left boundary characteristic value cL=1 is the critical state of hunting stability.Besides,stochastic D-bifurcation and P-bifurcation will appear in the wheelset system,and the critical speeds of the two kinds of stochastic bifurcation decrease with the increase in the stochastic parametric excitation intensity.展开更多
A new fuzzy stochastic finite element method based on the fuzzy factor method and random factor method is given and the analysis of structural dynamic characteristic for fuzzy stochastic truss structures is presented....A new fuzzy stochastic finite element method based on the fuzzy factor method and random factor method is given and the analysis of structural dynamic characteristic for fuzzy stochastic truss structures is presented. Considering the fuzzy randomness of the structural physical parameters and geometric dimensions simultaneously, the structural stiffness and mass matrices axe constructed based on the fuzzy factor method and random factor method; from the Rayleigh's quotient of structural vibration, the structural fuzzy random dynamic characteristic is obtained by means of the interval arithmetic; the fuzzy numeric characteristics of dynamic characteristic axe then derived by using the random variable's moment function method and algebra synthesis method. Two examples axe used to illustrate the validity and rationality of the method given. The advantage of this method is that the effect of the fuzzy randomness of one of the structural parameters on the fuzzy randomness of the dynamic characteristic can be reflected expediently and objectively.展开更多
To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random ...To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory (MATLAB) sottware, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by compar- ing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatial- ly random fields.展开更多
Nonlinear and stochastic dynamics is essential and practical in almost all branches of natural science and engineering.It has been a central subject to understand various complex dynamics,such as random vibration,stoc...Nonlinear and stochastic dynamics is essential and practical in almost all branches of natural science and engineering.It has been a central subject to understand various complex dynamics,such as random vibration,stochastic transition,synchronization,et al.,in the areas of mechanical and aerospace engineering,physics and chemistry.For example,the stochastic resonance has been utilized effectively in mechanical fault diagnosis and weak signal detection.Naturally,many methods have been developed to get the stochastic responses,including stochastic averaging method.展开更多
The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed....The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed.Random sea surface is treated as a narrow-band stochastic process,and the stochastic parametric excitation is studied based on the effective wave theory.The nonlinear restored arm function obtained from the numerical simulation is expressed as the approximate analytic function.By using the stochastic averaging method,the differential equation of motion is transformed into Ito’s stochastic differential equation.The steady-state probability density function of roll motion is obtained,and the results are validated with the numerical simulation and model test.展开更多
The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansi...The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displace- ments are derived, respectively. It is found that the randomness of material param- eters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods.展开更多
Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperat...Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperature fields and thermal stresses.Therefore,it is reasonable to consider the variability of thermal environments while conducting thermal analysis.However,for ambient thermal excitations,only stationary random processes have been investigated thus far.In this study,the highly efficient explicit time-domain method(ETDM)is proposed for the analysis of non-stationary stochastic transient heat conduction and thermal stress problems.The explicit time-domain expressions of thermal responses are first constructed for a thermoelastic body.Then the statistical moments of thermal displacements and stresses can be directly obtained based on the explicit expressions of thermal responses.A numerical example involving non-stationary stochastic internal heat generation rate is investigated.The accuracy and efficiency of the proposed method are validated by comparison with the Monte-Carlo simulation.展开更多
This paper considers distributed stochastic optimization,in which a number of agents cooperate to optimize a global objective function through local computations and information exchanges with neighbors over a network...This paper considers distributed stochastic optimization,in which a number of agents cooperate to optimize a global objective function through local computations and information exchanges with neighbors over a network.Stochastic optimization problems are usually tackled by variants of projected stochastic gradient descent.However,projecting a point onto a feasible set is often expensive.The Frank-Wolfe(FW)method has well-documented merits in handling convex constraints,but existing stochastic FW algorithms are basically developed for centralized settings.In this context,the present work puts forth a distributed stochastic Frank-Wolfe solver,by judiciously combining Nesterov's momentum and gradient tracking techniques for stochastic convex and nonconvex optimization over networks.It is shown that the convergence rate of the proposed algorithm is O(k^(-1/2))for convex optimization,and O(1/log_(2)(k))for nonconvex optimization.The efficacy of the algorithm is demonstrated by numerical simulations against a number of competing alternatives.展开更多
Nonlinear stochastic modeling plays a significant role in disciplines such as psychology,finance,physical sciences,engineering,econometrics,and biological sciences.Dynamical consistency,positivity,and boundedness are ...Nonlinear stochastic modeling plays a significant role in disciplines such as psychology,finance,physical sciences,engineering,econometrics,and biological sciences.Dynamical consistency,positivity,and boundedness are fundamental properties of stochastic modeling.A stochastic coronavirus model is studied with techniques of transition probabilities and parametric perturbation.Well-known explicit methods such as Euler Maruyama,stochastic Euler,and stochastic Runge–Kutta are investigated for the stochastic model.Regrettably,the above essential properties are not restored by existing methods.Hence,there is a need to construct essential properties preserving the computational method.The non-standard approach of finite difference is examined to maintain the above basic features of the stochastic model.The comparison of the results of deterministic and stochastic models is also presented.Our proposed efficient computational method well preserves the essential properties of the model.Comparison and convergence analyses of the method are presented.展开更多
In this paper, the cubic and quintic diffusion equation under stochastic non homogeneity is solved using Wiener- Hermite expansion and perturbation (WHEP) technique, Homotopy perturbation method (HPM) and Pickard appr...In this paper, the cubic and quintic diffusion equation under stochastic non homogeneity is solved using Wiener- Hermite expansion and perturbation (WHEP) technique, Homotopy perturbation method (HPM) and Pickard approximation technique. The analytic solution of the linear case is obtained using Eigenfunction expansion .The Picard approximation method is used to introduce the first and second order approximate solution for the non linear case. The WHEP technique is also used to obtain approximate solution under different orders and different corrections. The Homotopy perturbation method (HPM) is also used to obtain some approximation orders for mean and variance. Using mathematica-5, the methods of solution are illustrated through figures, comparisons among different methods and some parametric studies.展开更多
In this paper,we establish a large deviation principle for the stochastic generalized Ginzburg-Landau equation driven by jump noise.The main difficulties come from the highly non-linear coefficient and the jump noise....In this paper,we establish a large deviation principle for the stochastic generalized Ginzburg-Landau equation driven by jump noise.The main difficulties come from the highly non-linear coefficient and the jump noise.Here,we adopt a new sufficient condition for the weak convergence criterion of the large deviation principle,which was initially proposed by Matoussi,Sabbagh and Zhang(2021).展开更多
We propose a novel stochastic modeling framework for coal production and logistics using option pricing theory.The problem of valuing the inherent real optionality a coal producer has when mining and processing therma...We propose a novel stochastic modeling framework for coal production and logistics using option pricing theory.The problem of valuing the inherent real optionality a coal producer has when mining and processing thermal coal is modelled as pricing spread options of three assets under the stochastic volatility model.We derive a three-dimensional Fast Fourier Transform(“FFT”)lower bound approximation to value the inherent real optionality and for robustness check,we compare the semi-analytical pricing accuracy with the Monte Carlo simulation.Model parameters are estimated from the historical monthly data,and stochastic volatility parameters are obtained by matching the Kurtosis of the low-ash diff data to the Kurtosis of the stochastic volatility process which is assumed to follow Cox–Ingersoll–Ross(“CIR”)model.展开更多
The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system ...The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system by orthogonal polynomial approximation. Then, the linear stability theory and Routh-Hurwitz criterion for nonlinear deterministic systems are applied to the equivalent one. At last, at the aid of Lyapunov second method, we obtain that as the random intensity or statistical parameter of random variable is changed, the stability about stochastic Lotka-Volterra model is different from the deterministic system.展开更多
Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong conve...Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong convergence can preserve the mean square stability with the sufficiently small stepsize. Weak variants and their mean square stability are also considered. Several numerical experiments are given for illustration and show that the fully implicit methods are superior to those of the explicit methods in terms of mean-square stabilities for relatively large stepsizes especially.展开更多
Utilizing the well-known aggregation technique, we propose a smoothing sample average approximation (SAA) method for a stochastic linear complementarity problem, where the underlying functions are represented by exp...Utilizing the well-known aggregation technique, we propose a smoothing sample average approximation (SAA) method for a stochastic linear complementarity problem, where the underlying functions are represented by expectations of stochastic functions. The method is proved to be convergent and the preliminary numerical results are reported.展开更多
文摘We developed a modified stochastic finite-fault method for estimating strong ground motions.An adjustment to the dynamic corner frequency was introduced,which accounted for the effect of the location of the subfault relative to the hypocenter and rupture propagation direction,to account for the influence of the rupture propagation direction on the subfault dynamic corner frequency.By comparing the peak ground acceleration(PGA),pseudo-absolute response spectra acceleration(PSA,damping ratio of 5%),and duration,the results of the modified and existing methods were compared,demonstrating that our proposed adjustment to the dynamic corner frequency can accurately reflect the rupture directivity effect.We applied our modified method to simulate near-field strong motions within 150 km of the 2008 MW7.9 Wenchuan earthquake rupture.Our modified method performed well over a broad period range,particularly at 0.04-4 s.The total deviations of the stochastic finite-fault method(EXSIM)and the modified EXSIM were 0.1676 and 0.1494,respectively.The modified method can effectively account for the influence of the rupture propagation direction and provide more realistic ground motion estimations for earthquake disaster mitigation.
文摘Proximal gradient descent and its accelerated version are resultful methods for solving the sum of smooth and non-smooth problems. When the smooth function can be represented as a sum of multiple functions, the stochastic proximal gradient method performs well. However, research on its accelerated version remains unclear. This paper proposes a proximal stochastic accelerated gradient (PSAG) method to address problems involving a combination of smooth and non-smooth components, where the smooth part corresponds to the average of multiple block sums. Simultaneously, most of convergence analyses hold in expectation. To this end, under some mind conditions, we present an almost sure convergence of unbiased gradient estimation in the non-smooth setting. Moreover, we establish that the minimum of the squared gradient mapping norm arbitrarily converges to zero with probability one.
基金supported by the National Natural Science Foundation of China(Grant No.52109010)the Postdoctoral Science Foundation of China(Grant No.2021M701047)the China National Postdoctoral Program for Innovative Talents(Grant No.BX20200113).
文摘Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.
文摘In this paper, we investigate an SIS model with treatment and immigration. Firstly, the two-dimensional model is simplified by using the stochastic averaging method. Then, we derive the local stability of the stochastic system by computing the Lyapunov exponent of the linearized system. Further, the global stability of the stochastic model is analyzed based on the singular boundary theory. Moreover, we prove that the model undergoes a Hopf bifurcation and a pitchfork bifurcation. Finally, several numerical examples are provided to illustrate the theoretical results. .
基金This research is supported by the Deputyship forResearch&Innovation,Ministry of Education in Saudi Arabia under Project Number(IFP-2022-35).
文摘To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.
基金Project supported by the National Natural Science Foundation of China(Nos.11790282,12172235,12072208,and 52072249)the Opening Foundation of State Key Laboratory of Shijiazhuang Tiedao University of China(No.ZZ2021-13)。
文摘A stochastic wheelset model with a nonlinear wheel-rail contact relationship is established to investigate the stochastic stability and stochastic bifurcation of the wheelset system with the consideration of the stochastic parametric excitations of equivalent conicity and suspension stiffness.The wheelset is systematized into a onedimensional(1D)diffusion process by using the stochastic average method,the behavior of the singular boundary is analyzed to determine the hunting stability condition of the wheelset system,and the critical speed of stochastic bifurcation is obtained.The stationary probability density and joint probability density are derived theoretically.Based on the topological structure change of the probability density function,the stochastic Hopf bifurcation form and bifurcation condition of the wheelset system are determined.The effects of stochastic factors on the hunting stability and bifurcation characteristics are analyzed,and the simulation results verify the correctness of the theoretical analysis.The results reveal that the boundary behavior of the diffusion process determines the hunting stability of the stochastic wheelset system,and the left boundary characteristic value cL=1 is the critical state of hunting stability.Besides,stochastic D-bifurcation and P-bifurcation will appear in the wheelset system,and the critical speeds of the two kinds of stochastic bifurcation decrease with the increase in the stochastic parametric excitation intensity.
基金Project supported by the Natural Science Foundation of Shaanxi Province of China (No,A200214)
文摘A new fuzzy stochastic finite element method based on the fuzzy factor method and random factor method is given and the analysis of structural dynamic characteristic for fuzzy stochastic truss structures is presented. Considering the fuzzy randomness of the structural physical parameters and geometric dimensions simultaneously, the structural stiffness and mass matrices axe constructed based on the fuzzy factor method and random factor method; from the Rayleigh's quotient of structural vibration, the structural fuzzy random dynamic characteristic is obtained by means of the interval arithmetic; the fuzzy numeric characteristics of dynamic characteristic axe then derived by using the random variable's moment function method and algebra synthesis method. Two examples axe used to illustrate the validity and rationality of the method given. The advantage of this method is that the effect of the fuzzy randomness of one of the structural parameters on the fuzzy randomness of the dynamic characteristic can be reflected expediently and objectively.
基金funded by the National Basic Research Program of China (No. 2012CB026103)the National High Technology Research and Development Program of China (No. 2012AA06A401)the National Natural Science Foundation of China (No. 41271096)
文摘To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory (MATLAB) sottware, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by compar- ing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatial- ly random fields.
文摘Nonlinear and stochastic dynamics is essential and practical in almost all branches of natural science and engineering.It has been a central subject to understand various complex dynamics,such as random vibration,stochastic transition,synchronization,et al.,in the areas of mechanical and aerospace engineering,physics and chemistry.For example,the stochastic resonance has been utilized effectively in mechanical fault diagnosis and weak signal detection.Naturally,many methods have been developed to get the stochastic responses,including stochastic averaging method.
基金the State Administration of Science,Technology and Industry for National Defense of China(Grant No.B2420132001).
文摘The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed.Random sea surface is treated as a narrow-band stochastic process,and the stochastic parametric excitation is studied based on the effective wave theory.The nonlinear restored arm function obtained from the numerical simulation is expressed as the approximate analytic function.By using the stochastic averaging method,the differential equation of motion is transformed into Ito’s stochastic differential equation.The steady-state probability density function of roll motion is obtained,and the results are validated with the numerical simulation and model test.
文摘The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displace- ments are derived, respectively. It is found that the randomness of material param- eters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods.
基金funded by the National Natural Science Foundation of China (51678252)the Guangzhou Science and Technology Project (201804020069)
文摘Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperature fields and thermal stresses.Therefore,it is reasonable to consider the variability of thermal environments while conducting thermal analysis.However,for ambient thermal excitations,only stationary random processes have been investigated thus far.In this study,the highly efficient explicit time-domain method(ETDM)is proposed for the analysis of non-stationary stochastic transient heat conduction and thermal stress problems.The explicit time-domain expressions of thermal responses are first constructed for a thermoelastic body.Then the statistical moments of thermal displacements and stresses can be directly obtained based on the explicit expressions of thermal responses.A numerical example involving non-stationary stochastic internal heat generation rate is investigated.The accuracy and efficiency of the proposed method are validated by comparison with the Monte-Carlo simulation.
基金supported in part by the National Key R&D Program of China(2021YFB1714800)the National Natural Science Foundation of China(62222303,62073035,62173034,61925303,62088101,61873033)+1 种基金the CAAI-Huawei MindSpore Open Fundthe Chongqing Natural Science Foundation(2021ZX4100027)。
文摘This paper considers distributed stochastic optimization,in which a number of agents cooperate to optimize a global objective function through local computations and information exchanges with neighbors over a network.Stochastic optimization problems are usually tackled by variants of projected stochastic gradient descent.However,projecting a point onto a feasible set is often expensive.The Frank-Wolfe(FW)method has well-documented merits in handling convex constraints,but existing stochastic FW algorithms are basically developed for centralized settings.In this context,the present work puts forth a distributed stochastic Frank-Wolfe solver,by judiciously combining Nesterov's momentum and gradient tracking techniques for stochastic convex and nonconvex optimization over networks.It is shown that the convergence rate of the proposed algorithm is O(k^(-1/2))for convex optimization,and O(1/log_(2)(k))for nonconvex optimization.The efficacy of the algorithm is demonstrated by numerical simulations against a number of competing alternatives.
基金the Research and initiative center COVID-19-DES-2020-65,Prince Sultan University.
文摘Nonlinear stochastic modeling plays a significant role in disciplines such as psychology,finance,physical sciences,engineering,econometrics,and biological sciences.Dynamical consistency,positivity,and boundedness are fundamental properties of stochastic modeling.A stochastic coronavirus model is studied with techniques of transition probabilities and parametric perturbation.Well-known explicit methods such as Euler Maruyama,stochastic Euler,and stochastic Runge–Kutta are investigated for the stochastic model.Regrettably,the above essential properties are not restored by existing methods.Hence,there is a need to construct essential properties preserving the computational method.The non-standard approach of finite difference is examined to maintain the above basic features of the stochastic model.The comparison of the results of deterministic and stochastic models is also presented.Our proposed efficient computational method well preserves the essential properties of the model.Comparison and convergence analyses of the method are presented.
文摘In this paper, the cubic and quintic diffusion equation under stochastic non homogeneity is solved using Wiener- Hermite expansion and perturbation (WHEP) technique, Homotopy perturbation method (HPM) and Pickard approximation technique. The analytic solution of the linear case is obtained using Eigenfunction expansion .The Picard approximation method is used to introduce the first and second order approximate solution for the non linear case. The WHEP technique is also used to obtain approximate solution under different orders and different corrections. The Homotopy perturbation method (HPM) is also used to obtain some approximation orders for mean and variance. Using mathematica-5, the methods of solution are illustrated through figures, comparisons among different methods and some parametric studies.
基金partially supported by the National Natural Science Foundation of China(11871382,12071361)partially supported by the National Natural Science Foundation of China(11971361,11731012)。
文摘In this paper,we establish a large deviation principle for the stochastic generalized Ginzburg-Landau equation driven by jump noise.The main difficulties come from the highly non-linear coefficient and the jump noise.Here,we adopt a new sufficient condition for the weak convergence criterion of the large deviation principle,which was initially proposed by Matoussi,Sabbagh and Zhang(2021).
文摘We propose a novel stochastic modeling framework for coal production and logistics using option pricing theory.The problem of valuing the inherent real optionality a coal producer has when mining and processing thermal coal is modelled as pricing spread options of three assets under the stochastic volatility model.We derive a three-dimensional Fast Fourier Transform(“FFT”)lower bound approximation to value the inherent real optionality and for robustness check,we compare the semi-analytical pricing accuracy with the Monte Carlo simulation.Model parameters are estimated from the historical monthly data,and stochastic volatility parameters are obtained by matching the Kurtosis of the low-ash diff data to the Kurtosis of the stochastic volatility process which is assumed to follow Cox–Ingersoll–Ross(“CIR”)model.
文摘The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system by orthogonal polynomial approximation. Then, the linear stability theory and Routh-Hurwitz criterion for nonlinear deterministic systems are applied to the equivalent one. At last, at the aid of Lyapunov second method, we obtain that as the random intensity or statistical parameter of random variable is changed, the stability about stochastic Lotka-Volterra model is different from the deterministic system.
基金National Natural Science Foundations of China(Nos.11561028,11101101,11461032,11401267)Natural Science Foundations of Jiangxi Province,China(Nos.20151BAB201013,20151BAB201010,20151BAB201015)
文摘Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong convergence can preserve the mean square stability with the sufficiently small stepsize. Weak variants and their mean square stability are also considered. Several numerical experiments are given for illustration and show that the fully implicit methods are superior to those of the explicit methods in terms of mean-square stabilities for relatively large stepsizes especially.
文摘Utilizing the well-known aggregation technique, we propose a smoothing sample average approximation (SAA) method for a stochastic linear complementarity problem, where the underlying functions are represented by expectations of stochastic functions. The method is proved to be convergent and the preliminary numerical results are reported.