Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the ...In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
Recently,inspired by a modified generalized shift-splitting iteration method for complex symmetric linear systems,we propose two variants of the modified generalized shift-splitting iteration(MGSS)methods for solving ...Recently,inspired by a modified generalized shift-splitting iteration method for complex symmetric linear systems,we propose two variants of the modified generalized shift-splitting iteration(MGSS)methods for solving com-plex symmetric linear systems.One is a parameterized MGSS iteration method and the other is a modified parameterized MGSS iteration method.We prove that the proposed methods are convergent under appropriate constraints on the parameters.In addition,we also give the eigenvalue distributions of differ-ent preconditioned matrices to verify the effectiveness of the preconditioners proposed in this paper.展开更多
A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a n...A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a new mitigation method is urgently expected to be explored.This study proposed a novel asymptotic linear method(ALM)for micro pressure wave(MPW)mitigation to achieve a constant gradient of initial c ompression waves(ICWs),via a study with various open ratios on hoods.The properties of ICWs and MPWs under various open ratios of hoods were analyzed.The results show that as the open ratio increases,the MPW amplitude at the tunnel exit initially decreases before rising.At the open ratio of 2.28%,the slope of the ICW curve is linearly coincident with a supposed straight line in the ALM,which further reduces the MPW amplitude by 26.9%at 20 m and 20.0%at 50 m from the exit,as compared to the unvented hood.Therefore,the proposed method effectively mitigates MPW and quickly determines the upper limit of alleviation for the MPW amplitude at a fixed train-tunnel operation condition.All achievements provide a ne w potential measure for the adaptive design of tunnel hoods.展开更多
In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.1...In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.134:108344,2022).The convergence analysis of the proposed method is established,which can improve the existing results.Numerical examples show that the proposed method is efficient with the two-step technique.展开更多
In this paper, two kinds of chaotic systems are controlled respectively with and without time-delay to eliminate their chaotic behaviors. First of all, according to the first-order approximation method and the stabili...In this paper, two kinds of chaotic systems are controlled respectively with and without time-delay to eliminate their chaotic behaviors. First of all, according to the first-order approximation method and the stabilization condition of the linear system, one linear feedback controller is structured to control the chaotic system without time-delay, its chaotic behavior is eliminated and stabilized to its equilibrium. After that, based on the first-order approximation method, the Lyapunov stability theorem, and the matrix inequality theory, the other linear feedback controller is structured to control the chaotic system with time-delay and make it stabilized at its equilibrium. Finally, two numerical examples are given to illustrate the correctness and effectiveness of the two linear feedback controllers.展开更多
This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)metho...This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.展开更多
This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either...This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.展开更多
Piecewise linear systems are prevalent in engineering practice,and can be categorized into symmetric and asymmetric piecewise linear systems.Considering that symmetry is a special case of asymmetry,it is essential to ...Piecewise linear systems are prevalent in engineering practice,and can be categorized into symmetric and asymmetric piecewise linear systems.Considering that symmetry is a special case of asymmetry,it is essential to investigate the broader model,namely the asymmetric piecewise linear system.The traditional averaging method is frequently used for studying nonlinear systems,particularly symmetric piecewise linear systems,with the harmonic response of the oscillator serving as a key prerequisite for calculating steady-state solutions.However,asymmetric systems inherently exhibit nonharmonic,asymmetric responses,rendering the traditional averaging method inapplicable.To overcome this limitation,this paper introduces an improved averaging method tailored for an oscillator characterized by asymmetric gaps and springs.Unlike the traditional method,which assumes a purely harmonic response,the improved averaging method redefines the system response as a superposition of a direct current(DC)term and a first harmonic component.Herein,the DC term can be regarded as the offset induced by model asymmetry.Furthermore,the DC term is treated as a slow variable function of time,with its time derivative assumed to be zero when calculating the steady-state solution,akin to the amplitude and phase in the traditional averaging method.Numerical validation demonstrates that the responses computed in both time and frequency domains with the improved averaging method exhibit greater accuracy compared with those derived from the traditional method.Leveraging these improved results,the study also examines the parameter effect,stability,and bifurcation phenomena within the amplitude-frequency responses.展开更多
The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the soluti...The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the solution,characterized by||u_(n)||L^(2)(Ω)=O(t^(−αs)) as t→∞,is determined by the minimum order α_(s) of the time-fractional derivatives.Building on this foundational result,this article pursues two primary objectives.First,we introduce a strongly A-stable fractional linear multistep method and derive the numerical stability region for the governing equation.Second,we rigorously prove the long-term decay rate of the numerical solution through a detailed singularity analysis of its generating function.Notably,the numerical decay rate||u_(n)||L^(2)(Ω)=O(t_(n)^(−α_(s)) as t_(n)→∞aligns precisely with the continuous case.Theoretical findings are further validated through comprehensive numerical simulations,underscoring the robustness of our proposed method.展开更多
With the development of science and technology,the design and optimization of control systems are widely applied.This paper focuses on the application of matrix equations in linear time-invariant systems.Taking the in...With the development of science and technology,the design and optimization of control systems are widely applied.This paper focuses on the application of matrix equations in linear time-invariant systems.Taking the inverted pendulum model as an example,the algebraic Riccati equation is used to solve the optimal control problem,and the system performance and stability are achieved by selecting the closed-loop pole and designing the gain matrix.Then,the numerical methods for solving the stochastic algebraic Riccati equations are applied to practical problems,with Newton’s iteration method as the outer iteration and the solution of the mixed-type Lyapunov equations as the inner iteration.Two methods for solving the Lyapunov equations are introduced,providing references for related research.展开更多
Wave shoaling,which involves an increase in wave amplitude due to changes in water depth,can damage shore-lines.To mitigate this damage,we propose using porous structures such as mangrove forests.In this study,we use ...Wave shoaling,which involves an increase in wave amplitude due to changes in water depth,can damage shore-lines.To mitigate this damage,we propose using porous structures such as mangrove forests.In this study,we use a mathematical model to examine how mangroves,acting as porous breakwater,can reduce wave shoaling amplitude.The shallow water equations are used as the governing equations and are modified to account for the presence of porous media.To measure the wave reduction generated by the porous media,the wave transmis-sion coefficient is estimated using analytical and numerical approaches.The separation of variables method and the staggered finite volume method are utilized for each approach,respectively.The numerical results are then validated against the previously obtained analytical solutions.We then vary the friction and porosity parame-ters-determined by the presence and extent of porous media,to evaluate their effectiveness in reducing wave shoaling.展开更多
We study the thermodynamic properties of the classical one-dimensional generalized nonlinear Klein-Gordon lattice model(n≥2)by using the cluster variation method with linear response theory.The results of this method...We study the thermodynamic properties of the classical one-dimensional generalized nonlinear Klein-Gordon lattice model(n≥2)by using the cluster variation method with linear response theory.The results of this method are exact in the thermodynamic limit.We present the single-site reduced densityρ^((1))(z),averages such as(z^(2)),<|z^(n)|>,and<(z_(1)-z_(2))^(2)>,the specific heat C_(v),and the static correlation functions.We analyze the scaling behavior of these quantities and obtain the exact scaling powers at the low and high temperatures.Using these results,we gauge the accuracy of the projective truncation approximation for theφ^(4)lattice model.展开更多
Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate ...Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.展开更多
In this study,we design and numerically investigate a novel all optical D flip-flop(AODFF)based on linear photonic crystal(LPhC)structure that is composed of optical waveguides using the finite difference time domain(...In this study,we design and numerically investigate a novel all optical D flip-flop(AODFF)based on linear photonic crystal(LPhC)structure that is composed of optical waveguides using the finite difference time domain(FDTD)method.The proposed structure has the hexagonal close packed of 16×20 circular rods that are suspended in the air substrate with a lattice constant of 606 nm.The plane wave expansion(PWE)method is used to obtain the band diagram for AODFF at an operating wavelength of 1550 nm.The proposed optical flip-flop achieves a low delay time of 0.2 ps and a high contrast ratio(CR)of 10.33 dB.The main advantage of this design is that the input power as low as 1 mW/μm^(2) is sufficient for its operation,since no nonlinear rods are included.In addition,the footprint of the proposed AODFF is 100μm^(2),which is smaller compared to the structures reported in the literature,and it has a fast switching frequency of 5 Tbit/s.展开更多
The internal residual stress within a TC 17 titanium alloy joint welded by linear friction welding （LFW） was measured by the contour method, which is a relatively new and destructive technique to obtain a full map o...The internal residual stress within a TC 17 titanium alloy joint welded by linear friction welding （LFW） was measured by the contour method, which is a relatively new and destructive technique to obtain a full map of internal residual stress. The specimen was first cut into two parts; the out-of-plane displacement contour formed by the release of the residual stress was then measured; finally, taking the measured contour of the cut plane as the boundary conditions, a linear elastic finite element analysis was carried out to calculate the corresponding distribution of residual stress normal to the cut plane. The internal stress distribution of the TC 17 titanium alloy LFWjoint was also analyzed. The results show that the tensile residual stress in the TC17 LFW weld is mainly present within a region about 12 mm from the weld centerline; the peak tensile residual stress occurs at the weld centerline and reaches 360 MPa （about one third of the yield strength of TC17 alloy）; within the weld zone of the TC17 LFW weld, the through-thickness stress is not uniform, and the internal stress is larger than that near the top or bottom surface.展开更多
With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） meth...With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.展开更多
The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect ...The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.展开更多
The degree of numerical linear independence is proposed and discussed. Based on this linear independence theory, a modified limited memory BFGS method is deve loped. Similar to the standard limited memory method, thi...The degree of numerical linear independence is proposed and discussed. Based on this linear independence theory, a modified limited memory BFGS method is deve loped. Similar to the standard limited memory method, this new method determines the new update by applying the updating formula m times to an initial positive diagonal matrix using the m previous pairs of the change in iteration and gradient. Besides the most recent pair of the change, which guarantees the quadratic termination, the choice of the other ( m -1) pairs of the change in the new method is dependent on the degree of numerical linear independence of previous search directions. In addition, the numerical linear independence theory is further discussed and the computation of the degree of linear independence is simplified. Theoretical and numerical results show that this new modified method improves efficiently the standard limited memory method.展开更多
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金Supported by the Beijing Municipal Science&Technology Commission(Z211100004421012),the Key Reaserch and Development Pro⁃gram of China(2022YFF0605902)。
文摘In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金supported by the National Natural Science Foundation of China(Grant No.12371378)by the Natural Science Foundation of Fujian Province(Grant Nos.2024J01980,2024J08242).
文摘Recently,inspired by a modified generalized shift-splitting iteration method for complex symmetric linear systems,we propose two variants of the modified generalized shift-splitting iteration(MGSS)methods for solving com-plex symmetric linear systems.One is a parameterized MGSS iteration method and the other is a modified parameterized MGSS iteration method.We prove that the proposed methods are convergent under appropriate constraints on the parameters.In addition,we also give the eigenvalue distributions of differ-ent preconditioned matrices to verify the effectiveness of the preconditioners proposed in this paper.
基金Project(24A0006)supported by the Key Project of Scientific Research Fund of Hunan Provincial Department of Education,ChinaProject(2024JJ5430)supported by the Natural Science Foundation of Hunan Province,ChinaProjects(2024JK2045,2023RC3061)supported by the Science and Technology Innovation Program of Hunan Province,China。
文摘A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a new mitigation method is urgently expected to be explored.This study proposed a novel asymptotic linear method(ALM)for micro pressure wave(MPW)mitigation to achieve a constant gradient of initial c ompression waves(ICWs),via a study with various open ratios on hoods.The properties of ICWs and MPWs under various open ratios of hoods were analyzed.The results show that as the open ratio increases,the MPW amplitude at the tunnel exit initially decreases before rising.At the open ratio of 2.28%,the slope of the ICW curve is linearly coincident with a supposed straight line in the ALM,which further reduces the MPW amplitude by 26.9%at 20 m and 20.0%at 50 m from the exit,as compared to the unvented hood.Therefore,the proposed method effectively mitigates MPW and quickly determines the upper limit of alleviation for the MPW amplitude at a fixed train-tunnel operation condition.All achievements provide a ne w potential measure for the adaptive design of tunnel hoods.
基金supported by the Scientific Computing Research Innovation Team of Guangdong Province(no.2021KCXTD052)the Science and Technology Development Fund,Macao SAR(no.0096/2022/A,0151/2022/A)+3 种基金University of Macao(no.MYRG2020-00035-FST,MYRG2022-00076-FST)the Guangdong Key Construction Discipline Research Capacity Enhancement Project(no.2022ZDJS049)Technology Planning Project of Shaoguan(no.210716094530390)the ScienceFoundation of Shaoguan University(no.SZ2020KJ01).
文摘In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.134:108344,2022).The convergence analysis of the proposed method is established,which can improve the existing results.Numerical examples show that the proposed method is efficient with the two-step technique.
基金Supported by the National Natural Science Foundation of China (61863022)the Natural Science Foundation of Gansu Province(20JR10RA329)Scientific Research and Innovation Fund Project of Gansu University of Chinese Medicine in 2019 (2019KCYB-10)。
文摘In this paper, two kinds of chaotic systems are controlled respectively with and without time-delay to eliminate their chaotic behaviors. First of all, according to the first-order approximation method and the stabilization condition of the linear system, one linear feedback controller is structured to control the chaotic system without time-delay, its chaotic behavior is eliminated and stabilized to its equilibrium. After that, based on the first-order approximation method, the Lyapunov stability theorem, and the matrix inequality theory, the other linear feedback controller is structured to control the chaotic system with time-delay and make it stabilized at its equilibrium. Finally, two numerical examples are given to illustrate the correctness and effectiveness of the two linear feedback controllers.
基金the National Natural Science Foundation of China(Grant Nos.71961022,11902163,12265020,and 12262024)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(Grant Nos.2019BS01011 and 2022MS01003)+5 种基金2022 Inner Mongolia Autonomous Region Grassland Talents Project-Young Innovative and Entrepreneurial Talents(Mingjing Du)2022 Talent Development Foundation of Inner Mongolia Autonomous Region of China(Ming-Jing Du)the Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region Program(Grant No.NJYT-20-B18)the Key Project of High-quality Economic Development Research Base of Yellow River Basin in 2022(Grant No.21HZD03)2022 Inner Mongolia Autonomous Region International Science and Technology Cooperation High-end Foreign Experts Introduction Project(Ge Kai)MOE(Ministry of Education in China)Humanities and Social Sciences Foundation(Grants No.20YJC860005).
文摘This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.
基金Project supported by the National Natural Science Foundation of China(Nos.12272242 and U1934201)。
文摘Piecewise linear systems are prevalent in engineering practice,and can be categorized into symmetric and asymmetric piecewise linear systems.Considering that symmetry is a special case of asymmetry,it is essential to investigate the broader model,namely the asymmetric piecewise linear system.The traditional averaging method is frequently used for studying nonlinear systems,particularly symmetric piecewise linear systems,with the harmonic response of the oscillator serving as a key prerequisite for calculating steady-state solutions.However,asymmetric systems inherently exhibit nonharmonic,asymmetric responses,rendering the traditional averaging method inapplicable.To overcome this limitation,this paper introduces an improved averaging method tailored for an oscillator characterized by asymmetric gaps and springs.Unlike the traditional method,which assumes a purely harmonic response,the improved averaging method redefines the system response as a superposition of a direct current(DC)term and a first harmonic component.Herein,the DC term can be regarded as the offset induced by model asymmetry.Furthermore,the DC term is treated as a slow variable function of time,with its time derivative assumed to be zero when calculating the steady-state solution,akin to the amplitude and phase in the traditional averaging method.Numerical validation demonstrates that the responses computed in both time and frequency domains with the improved averaging method exhibit greater accuracy compared with those derived from the traditional method.Leveraging these improved results,the study also examines the parameter effect,stability,and bifurcation phenomena within the amplitude-frequency responses.
文摘The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the solution,characterized by||u_(n)||L^(2)(Ω)=O(t^(−αs)) as t→∞,is determined by the minimum order α_(s) of the time-fractional derivatives.Building on this foundational result,this article pursues two primary objectives.First,we introduce a strongly A-stable fractional linear multistep method and derive the numerical stability region for the governing equation.Second,we rigorously prove the long-term decay rate of the numerical solution through a detailed singularity analysis of its generating function.Notably,the numerical decay rate||u_(n)||L^(2)(Ω)=O(t_(n)^(−α_(s)) as t_(n)→∞aligns precisely with the continuous case.Theoretical findings are further validated through comprehensive numerical simulations,underscoring the robustness of our proposed method.
基金Supported by National Natural Science Foundation of China(Grant No.12571388)the Visiting Scholar Program of National Natural Science Foundation of China(Grant No.12426616)Natural Science Research Start-up Foundation of Recruiting Talents of Nanjing University of Posts and Telecommunications(Grant No.NY223127).
文摘With the development of science and technology,the design and optimization of control systems are widely applied.This paper focuses on the application of matrix equations in linear time-invariant systems.Taking the inverted pendulum model as an example,the algebraic Riccati equation is used to solve the optimal control problem,and the system performance and stability are achieved by selecting the closed-loop pole and designing the gain matrix.Then,the numerical methods for solving the stochastic algebraic Riccati equations are applied to practical problems,with Newton’s iteration method as the outer iteration and the solution of the mixed-type Lyapunov equations as the inner iteration.Two methods for solving the Lyapunov equations are introduced,providing references for related research.
基金support from Program Riset Kolaborasi Indonesia(RKI)2024(Grant No.1841/IT1.B07.1/TA.00/2024).
文摘Wave shoaling,which involves an increase in wave amplitude due to changes in water depth,can damage shore-lines.To mitigate this damage,we propose using porous structures such as mangrove forests.In this study,we use a mathematical model to examine how mangroves,acting as porous breakwater,can reduce wave shoaling amplitude.The shallow water equations are used as the governing equations and are modified to account for the presence of porous media.To measure the wave reduction generated by the porous media,the wave transmis-sion coefficient is estimated using analytical and numerical approaches.The separation of variables method and the staggered finite volume method are utilized for each approach,respectively.The numerical results are then validated against the previously obtained analytical solutions.We then vary the friction and porosity parame-ters-determined by the presence and extent of porous media,to evaluate their effectiveness in reducing wave shoaling.
基金supported by the National Natural Science Foundation of China(Grant No.11974420).
文摘We study the thermodynamic properties of the classical one-dimensional generalized nonlinear Klein-Gordon lattice model(n≥2)by using the cluster variation method with linear response theory.The results of this method are exact in the thermodynamic limit.We present the single-site reduced densityρ^((1))(z),averages such as(z^(2)),<|z^(n)|>,and<(z_(1)-z_(2))^(2)>,the specific heat C_(v),and the static correlation functions.We analyze the scaling behavior of these quantities and obtain the exact scaling powers at the low and high temperatures.Using these results,we gauge the accuracy of the projective truncation approximation for theφ^(4)lattice model.
基金National Natural Science Foundation of China under Grant No.U2139208。
文摘Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.
文摘In this study,we design and numerically investigate a novel all optical D flip-flop(AODFF)based on linear photonic crystal(LPhC)structure that is composed of optical waveguides using the finite difference time domain(FDTD)method.The proposed structure has the hexagonal close packed of 16×20 circular rods that are suspended in the air substrate with a lattice constant of 606 nm.The plane wave expansion(PWE)method is used to obtain the band diagram for AODFF at an operating wavelength of 1550 nm.The proposed optical flip-flop achieves a low delay time of 0.2 ps and a high contrast ratio(CR)of 10.33 dB.The main advantage of this design is that the input power as low as 1 mW/μm^(2) is sufficient for its operation,since no nonlinear rods are included.In addition,the footprint of the proposed AODFF is 100μm^(2),which is smaller compared to the structures reported in the literature,and it has a fast switching frequency of 5 Tbit/s.
基金Project(35061107)supported by the Doctoral Initiation Project of Jiangsu University of Science and Technology,China
文摘The internal residual stress within a TC 17 titanium alloy joint welded by linear friction welding （LFW） was measured by the contour method, which is a relatively new and destructive technique to obtain a full map of internal residual stress. The specimen was first cut into two parts; the out-of-plane displacement contour formed by the release of the residual stress was then measured; finally, taking the measured contour of the cut plane as the boundary conditions, a linear elastic finite element analysis was carried out to calculate the corresponding distribution of residual stress normal to the cut plane. The internal stress distribution of the TC 17 titanium alloy LFWjoint was also analyzed. The results show that the tensile residual stress in the TC17 LFW weld is mainly present within a region about 12 mm from the weld centerline; the peak tensile residual stress occurs at the weld centerline and reaches 360 MPa （about one third of the yield strength of TC17 alloy）; within the weld zone of the TC17 LFW weld, the through-thickness stress is not uniform, and the internal stress is larger than that near the top or bottom surface.
基金The National Natural Science Foundation of China(No.60702027)the Free Research Fund of the National Mobile Communications Research Laboratory of Southeast University (No.2008B07)the National Basic Research Program of China(973 Program)(No.2007CB310603)
文摘With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.
文摘The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.
文摘The degree of numerical linear independence is proposed and discussed. Based on this linear independence theory, a modified limited memory BFGS method is deve loped. Similar to the standard limited memory method, this new method determines the new update by applying the updating formula m times to an initial positive diagonal matrix using the m previous pairs of the change in iteration and gradient. Besides the most recent pair of the change, which guarantees the quadratic termination, the choice of the other ( m -1) pairs of the change in the new method is dependent on the degree of numerical linear independence of previous search directions. In addition, the numerical linear independence theory is further discussed and the computation of the degree of linear independence is simplified. Theoretical and numerical results show that this new modified method improves efficiently the standard limited memory method.
