The modified Korteweg-de Vries （mKdV） typed equations can be used to describe certain nonlinear phenomena in fluids, plasmas, and optics. In this paper, the discretized mKdV lattice equation is investigated. With th...The modified Korteweg-de Vries （mKdV） typed equations can be used to describe certain nonlinear phenomena in fluids, plasmas, and optics. In this paper, the discretized mKdV lattice equation is investigated. With the aid of symbolic computation, the discrete matrix spectral problem for that system is constructed. Darboux transformation for that system is established based on the resulting spectral problem. Explicit solutions are derived via the Darboux transformation. Structures of those solutions are shown graphically, which might be helpful to understand some physical processes in fluids, plasmas, and optics.展开更多
In order to achieve the complex dynamic analysis of the self-propelled seafloor pilot miner moving on the seafloor of extremely cohesive soft soil and further to make it possible to integrate the miner system with som...In order to achieve the complex dynamic analysis of the self-propelled seafloor pilot miner moving on the seafloor of extremely cohesive soft soil and further to make it possible to integrate the miner system with some subsystems to form the complete integrated deep ocean mining pilot system and perform dynamic analysis, a new method for the dynamic modeling and analysis of the miner is proposed and developed in this paper, resulting in a simplified 3D single-body vehicle model with three translational and three rotational degrees of freedom, while the track-terrain interaction model is built by partitioning the track-terrain interface into discrete elements with parameterized force dements built on the theory of terramechanics acting on each discrete dement. To evaluate and verify the correctness and effectiveness of this new modeling and analysis method, typical comparative studies with regard to computational efficiency and solution accuracy are carried out between the traditional modeling method of building the tracked vehicle as a multi-body model and the new modeling method. In full consideration of the particMar structure design of the pilot miner, the special characteristics of the seafioor soil and the hydrodynamic force of near-seafloor currnt, the dynamic simulation analysis of the miner is performed and discussed, which can provide useful guidance and reference for the practical miner system in design and operation. This new method can not only realize the rapid dynamic simulation analysis of the miner but also make possible the integration and rapid dynamic analysis of the complete integrated deep ocean mining pilot system in further researches.展开更多
In this letter, we study discretized mKdV lattice equation by using a new generalized ansatz. As a result,many explicit rational exact solutions, including some new solitary wave solutions, are obtained by symbolic co...In this letter, we study discretized mKdV lattice equation by using a new generalized ansatz. As a result,many explicit rational exact solutions, including some new solitary wave solutions, are obtained by symbolic computation code Maple.展开更多
The rock fragmentation involves the inter-block and the intra-block fracture.A simulation method for rock fragmentation is developed by coupling Voronoi diagram(VD)and discretized virtual internal bond(DVIB).The DVIB ...The rock fragmentation involves the inter-block and the intra-block fracture.A simulation method for rock fragmentation is developed by coupling Voronoi diagram(VD)and discretized virtual internal bond(DVIB).The DVIB is a lattice model that consists of bonds.The VD is used to generate the potential block structure in the DVIB mesh.Each potential block may contain any number of bond cells.To characterize the inter-block fracture,a hyperelastic bond potential is employed for the bond cells that are cut by the VD edges.While to characterize the intra-block fracture,an elastobrittle bond potential is adopted for the bonds in a block.By this method,both the inter-block and intra-block fracture can be well simulated.The simulation results suggest that this method is a simple and efficient approach to rock fragmentation simulation with block smash.展开更多
This paper focuses on the robust stability for time-delay systems of neutral type. A new complete Lyapunov-Krasovskii function (LKF) is developed. Based on this function and discretization, stability conditions in ter...This paper focuses on the robust stability for time-delay systems of neutral type. A new complete Lyapunov-Krasovskii function (LKF) is developed. Based on this function and discretization, stability conditions in terms of linear matrix inequalities are obtained. A class of time-varying uncertainty of system matrices can be studied by the method.展开更多
Discrete software reliability measurement has a proper characteristic for describing a software reliability growth process which depends on a unit of the software fault-detection period, such as the number of test run...Discrete software reliability measurement has a proper characteristic for describing a software reliability growth process which depends on a unit of the software fault-detection period, such as the number of test runs, the number of executed test cases. This paper discusses discrete software reliability measurement based on a discretized nonhomogeneous Poisson process (NHPP) model. Especially, we use a bootstrapping method in our discrete software reliability measurement for discussing the statistical inference on parameters and software reliability assessment measures of our model. Finally we show numerical examples of interval estimations based on our bootstrapping method for the several software reliability assessment measures by using actual data.展开更多
For nonlinear stability problems of discretized conservative systems with multiple parameter variables and multiple state variables, the activation method is put forward, by which activated potential functions and act...For nonlinear stability problems of discretized conservative systems with multiple parameter variables and multiple state variables, the activation method is put forward, by which activated potential functions and activated equilibrium equations are derived. The activation method is the improvement and enhancement of Liapunov-Schmidt method in elastic stability theory. It is more generalized and more normalized than conventional perturbation methods. The activated potential functions may be transformed into normalized catastrophe potential functions. The activated equilibrium equations may be treated as bifurcation equations. The researches in this paper will motivate the combination of elastic stability theory with catastrophe theory and bifurcation theory展开更多
The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
In this paper, a discretized SIR model with pulse vaccination and time delay is proposed. We introduce two thresholds R* and R<sub>*</sub>, and further prove that the disease-free periodic solution is glob...In this paper, a discretized SIR model with pulse vaccination and time delay is proposed. We introduce two thresholds R* and R<sub>*</sub>, and further prove that the disease-free periodic solution is globally attractive if R* is less than unit and the disease can invade if R<sub>*</sub> is larger than unit. The numerical simulations not only illustrate the validity of our main results, but also exhibit bifurcation phenomenon. Our result shows that decreasing infection rate can put off the disease outbreak and reduce the number of infected individuals.展开更多
We present a systematic study of 6Li elastic scattering and total reaction cross sections at incident energies around the Coulomb barrier within the continuum discretized coupled-channels(CDCC)framework,where 6Li is t...We present a systematic study of 6Li elastic scattering and total reaction cross sections at incident energies around the Coulomb barrier within the continuum discretized coupled-channels(CDCC)framework,where 6Li is treated in anα+d two-body model.Collisions with 27Al,64Zn,138Ba,and 208Pa are analyzed.The microscopic optical potentials(MOP)based on Skyrme nucleon-nucleon interaction forαand d are adopted in CDCC calculations and satisfactory agreement with the experimental data is obtained without any adjustment on MOPs.For comparison,αand d global phenomenological optical potentials(GOP)are also used in CDCC analysis and a reduction of no less than 50%on the surface imaginary part of deuteron GOP is required for describing the data.In all cases,the 6Li breakup effect is significant and provides repulsive correction to the folding model potential.The reduction on the surface imaginary part of GOP of deuteron reveals a strong suppression of the reaction probability of deuteron as a component of 6Li when compared with that of a free deuteron.Further investigation is performed by considering the d breakup process equivalently within the dynamic polarization potential approach,and the results show that d behaves in a manner similar to a tightly bound nucleus in 6Li induced reactions.展开更多
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
This paper presents an Euler discretized inertial delayed neuron model, and its bifurcation dynamical behaviors are discussed. By using the associated characteristic model, center manifold theorem and the normal form ...This paper presents an Euler discretized inertial delayed neuron model, and its bifurcation dynamical behaviors are discussed. By using the associated characteristic model, center manifold theorem and the normal form method, it is shown that the model not only undergoes codimension one(flip, Neimark-Sacker) bifurcation, but also undergoes cusp and resonance bifurcation(1:1 and 1:2) of codimension two. Further, it is found that the parity of delay has some effect on bifurcation behaviors. Finally, some numerical simulations are given to support the analytic results and explore complex dynamics, such as periodic orbits near homoclinic orbits, quasiperiodic orbits, and chaotic orbits.展开更多
Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hyd...Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.展开更多
Fiber-reinforced composites are an ideal material for the lightweight design of aerospace structures. Especially in recent years, with the rapid development of composite additive manufacturing technology, the design o...Fiber-reinforced composites are an ideal material for the lightweight design of aerospace structures. Especially in recent years, with the rapid development of composite additive manufacturing technology, the design optimization of variable stiffness of fiber-reinforced composite laminates has attracted widespread attention from scholars and industry. In these aerospace composite structures, numerous cutout panels and shells serve as access points for maintaining electrical, fuel, and hydraulic systems. The traditional fiber-reinforced composite laminate subtractive drilling manufacturing inevitably faces the problems of interlayer delamination, fiber fracture, and burr of the laminate. Continuous fiber additive manufacturing technology offers the potential for integrated design optimization and manufacturing with high structural performance. Considering the integration of design and manufacturability in continuous fiber additive manufacturing, the paper proposes linear and nonlinear filtering strategies based on the Normal Distribution Fiber Optimization (NDFO) material interpolation scheme to overcome the challenge of discrete fiber optimization results, which are difficult to apply directly to continuous fiber additive manufacturing. With minimizing structural compliance as the objective function, the proposed approach provides a strategy to achieve continuity of discrete fiber paths in the variable stiffness design optimization of composite laminates with regular and irregular holes. In the variable stiffness design optimization model, the number of candidate fiber laying angles in the NDFO material interpolation scheme is considered as design variable. The sensitivity information of structural compliance with respect to the number of candidate fiber laying angles is obtained using the analytical sensitivity analysis method. Based on the proposed variable stiffness design optimization method for complex perforated composite laminates, the numerical examples consider the variable stiffness design optimization of typical non-perforated and perforated composite laminates with circular, square, and irregular holes, and systematically discuss the number of candidate discrete fiber laying angles, discrete fiber continuous filtering strategies, and filter radius on structural compliance, continuity, and manufacturability. The optimized discrete fiber angles of variable stiffness laminates are converted into continuous fiber laying paths using a streamlined process for continuous fiber additive manufacturing. Meanwhile, the optimized non-perforated and perforated MBB beams after discrete fiber continuous treatment, are manufactured using continuous fiber co-extrusion additive manufacturing technology to verify the effectiveness of the variable stiffness fiber optimization framework proposed in this paper.展开更多
This paper focuses on the use of rotary-percussive drilling for hard rocks.In order to improve efficiency and reduce costs,it is essential to understand how operational parameters,bit wear,and drilling performance are...This paper focuses on the use of rotary-percussive drilling for hard rocks.In order to improve efficiency and reduce costs,it is essential to understand how operational parameters,bit wear,and drilling performance are related.A model is presented therein that combines multibody dynamics and discrete element method(DEM)to investigate the influences of operational parameters and bit wear on the rate of penetration and wear characteristics.The model accurately captures the motion of the bit and recreates rock using the cutting sieving result.Field experimental results validate the rod dynamic behavior,rock recreating model,and coupling model in the simulation.The findings indicate that hammer pressure significantly influences the rate of penetration and wear depth of the bit,and there is an optimal range for economical hammer pressure.The wear coefficient has a major effect on the rate of penetration,when wear coefficient is between 1/3 and 2/3.Increasing the wear coefficient can reduce drill bit button pressure and wear depth at the same drill distance.Gauge button loss increases the rate of penetration due to higher pressure on the remaining buttons,which also accelerates destruction of the bit.Furthermore,a more evenly distributed button on the bit enhances the rate of penetration(ROP)when the same number of buttons is lost.展开更多
An H-polygon is a polygon with corners in the set of vertices of a tiling of R~2 by regular hexagons of unit edge.It is known that any H-triangle with exactly k interior H-points can have at most 3k+7 boundary H-point...An H-polygon is a polygon with corners in the set of vertices of a tiling of R~2 by regular hexagons of unit edge.It is known that any H-triangle with exactly k interior H-points can have at most 3k+7 boundary H-points.In this note we prove that any convex H-quadrilateral with exactly k interior H-points can have at most 3k+7 boundary H-points.展开更多
As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness dis...As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness distribution functions of the Bakhtiary dam site and Oskarshamn/Forsmark mountain were fitted using statistical methods.The COMSOL Multiphysics finite element software was utilized to analyze the effects of fracture roughness distribution types and empirical formulas for fracture hydraulic aperture on the seepage field and temperature field of rock masses.The results show that:(1)The fracture roughness at the Bakhtiary dam site and Oskarshamn/Forsmark mountain follows lognormal and normal distributions,respectively;(2)For rock masses with the same expected value and standard deviation of fracture roughness,the outflow from rock masses with lognormal distribution of fracture roughness is significantly larger than that of rock masses with normal distribution of fracture roughness;(3)The fracture hydraulic aperture,outflow,and cold front distance of the Li and Jiang model are significantly larger than those of the Barton model;(4)The outflow,hydraulic pressure distribution,and temperature distribution of the Barton model are more sensitive to the fracture roughness distribution type than those of the Li and Jiang model.展开更多
Precast concrete pavements(PCPs)represent an innovative solution in the construction industry,addressing the need for rapid,intelligent,and low-carbon pavement technologies that significantly reduce construction time ...Precast concrete pavements(PCPs)represent an innovative solution in the construction industry,addressing the need for rapid,intelligent,and low-carbon pavement technologies that significantly reduce construction time and environmental impact.However,the integration of prefabricated technology in pavement surface and base layers lacks systematic classification and understanding.This paper aims to fill this gap by introducing a detailed analysis of discretization and assembly connection technology for cement concrete pavement(CCP)structures.Through a comprehensive review of domestic and international literature,the study classifies prefabricated pavement technology based on discrete assembly structural layers and presents specific conclusions(i)surface layer discrete units are categorized into bottom plates,top plates,plate-rod separated assemblies,and prestressed connections,with optimal material compositions identified to enhance mechanical properties;(ii)base layer discrete units include block-type,plate-type,and beam-type elements,highlighting their contributions to sustainability by incorporating recycled materials(iii)planar assembly connection types are assessed,ranking them by load transfer efficiency,with specific dimensions provided for optimal performance;and(iv)vertical assembly connections are defined by their leveling and sealing layers,suitable for both new constructions and repairs of existing roads.The insights gained from this review not only clarify the distinctions between various structural layers but also provide practical guidelines for enhancing the design and implementation of PCP.This work contributes to advancing sustainable and resilient road construction practices,making it a significant reference for researchers and practitioners in the field.展开更多
Ensuring digital media security through robust image watermarking is essential to prevent unauthorized distribution,tampering,and copyright infringement.This study introduces a novel hybrid watermarking framework that...Ensuring digital media security through robust image watermarking is essential to prevent unauthorized distribution,tampering,and copyright infringement.This study introduces a novel hybrid watermarking framework that integrates Discrete Wavelet Transform(DWT),Redundant Discrete Wavelet Transform(RDWT),and Möbius Transformations(MT),with optimization of transformation parameters achieved via a Genetic Algorithm(GA).By combining frequency and spatial domain techniques,the proposed method significantly enhances both the imper-ceptibility and robustness of watermark embedding.The approach leverages DWT and RDWT for multi-resolution decomposition,enabling watermark insertion in frequency subbands that balance visibility and resistance to attacks.RDWT,in particular,offers shift-invariance,which improves performance under geometric transformations.Möbius transformations are employed for spatial manipulation,providing conformal mapping and spatial dispersion that fortify watermark resilience against rotation,scaling,and translation.The GA dynamically optimizes the Möbius parameters,selecting configurations that maximize robustness metrics such as Peak Signal-to-Noise Ratio(PSNR),Structural Similarity Index Measure(SSIM),Bit Error Rate(BER),and Normalized Cross-Correlation(NCC).Extensive experiments conducted on medical and standard benchmark images demonstrate the efficacy of the proposed RDWT-MT scheme.Results show that PSNR exceeds 68 dB,SSIM approaches 1.0,and BER remains at 0.0000,indicating excellent imperceptibility and perfect watermark recovery.Moreover,the method exhibits exceptional resilience to a wide range of image processing attacks,including Gaussian noise,JPEG compression,histogram equalization,and cropping,achieving NCC values close to or equal to 1.0.Comparative evaluations with state-of-the-art watermarking techniques highlight the superiority of the proposed method in terms of robustness,fidelity,and computational efficiency.The hybrid framework ensures secure,adaptive watermark embedding,making it highly suitable for applications in digital rights management,content authentication,and medical image protection.The integration of spatial and frequency domain features with evolutionary optimization presents a promising direction for future watermarking technologies.展开更多
This paper proposes a universal impulse-function-based method for extending discrete chaotic maps,enabling flexible construction of multicavity chaotic attractors.The proposed method achieves one-directional(1D)/two-d...This paper proposes a universal impulse-function-based method for extending discrete chaotic maps,enabling flexible construction of multicavity chaotic attractors.The proposed method achieves one-directional(1D)/two-directional(2D)extensions without introducing additional nonlinear terms or altering system stability.Theoretically,the cavity quantity in arbitrary directions is controlled by adjusting impulse levels N,while the amplitude regulation is implemented through modifications to the proportionality parameter r.Theoretical analyses,including Lyapunov exponents(LEs)and bifurcation diagrams,are conducted,confirming that the extended maps retain the intrinsic dynamics of five rational map classes.The field-programmable gate array(FPGA)implementation results are consistent with the numerical simulation results,verifying the correctness of the theoretical analysis.The method enables the expansion of unipolar attractors and enhances entropy metrics,offering a robust framework for applications in secure communication,encryption,and chaos-based technologies.展开更多
基金Supported by the National Natural Science Foundation of China under Grant No. 60772023by the Open Fund of the State Key Laboratory of Software Development Environment under Grant No. BUAA-SKLSDE-09KF-04+2 种基金Beijing University of Aeronautics and Astronautics, by the National Basic Research Program of China (973 Program) under Grant No. 2005CB321901the Specialized Research Fund for the Doctoral Program of Higher Education under Grant Nos. 20060006024 and 200800130006Chinese Ministry of Education, and Scientific Research Common Program of Beijing Municipal Commission of Education under Grant No. KM201010772020
文摘The modified Korteweg-de Vries （mKdV） typed equations can be used to describe certain nonlinear phenomena in fluids, plasmas, and optics. In this paper, the discretized mKdV lattice equation is investigated. With the aid of symbolic computation, the discrete matrix spectral problem for that system is constructed. Darboux transformation for that system is established based on the resulting spectral problem. Explicit solutions are derived via the Darboux transformation. Structures of those solutions are shown graphically, which might be helpful to understand some physical processes in fluids, plasmas, and optics.
基金supported by the National High Technology Research and Development Program of China(863 Program, Grant No.2006AA09Z240)the National Deep-Sea Technology Project of Development and Re-search(Grant No.DYXM-115-04-02-01)
文摘In order to achieve the complex dynamic analysis of the self-propelled seafloor pilot miner moving on the seafloor of extremely cohesive soft soil and further to make it possible to integrate the miner system with some subsystems to form the complete integrated deep ocean mining pilot system and perform dynamic analysis, a new method for the dynamic modeling and analysis of the miner is proposed and developed in this paper, resulting in a simplified 3D single-body vehicle model with three translational and three rotational degrees of freedom, while the track-terrain interaction model is built by partitioning the track-terrain interface into discrete elements with parameterized force dements built on the theory of terramechanics acting on each discrete dement. To evaluate and verify the correctness and effectiveness of this new modeling and analysis method, typical comparative studies with regard to computational efficiency and solution accuracy are carried out between the traditional modeling method of building the tracked vehicle as a multi-body model and the new modeling method. In full consideration of the particMar structure design of the pilot miner, the special characteristics of the seafioor soil and the hydrodynamic force of near-seafloor currnt, the dynamic simulation analysis of the miner is performed and discussed, which can provide useful guidance and reference for the practical miner system in design and operation. This new method can not only realize the rapid dynamic simulation analysis of the miner but also make possible the integration and rapid dynamic analysis of the complete integrated deep ocean mining pilot system in further researches.
基金the National Key Basic Research Project of China under
文摘In this letter, we study discretized mKdV lattice equation by using a new generalized ansatz. As a result,many explicit rational exact solutions, including some new solitary wave solutions, are obtained by symbolic computation code Maple.
基金the National Natural ScienceFoundation of China(Grant 11772190),which is gratefully acknowledged.
文摘The rock fragmentation involves the inter-block and the intra-block fracture.A simulation method for rock fragmentation is developed by coupling Voronoi diagram(VD)and discretized virtual internal bond(DVIB).The DVIB is a lattice model that consists of bonds.The VD is used to generate the potential block structure in the DVIB mesh.Each potential block may contain any number of bond cells.To characterize the inter-block fracture,a hyperelastic bond potential is employed for the bond cells that are cut by the VD edges.While to characterize the intra-block fracture,an elastobrittle bond potential is adopted for the bonds in a block.By this method,both the inter-block and intra-block fracture can be well simulated.The simulation results suggest that this method is a simple and efficient approach to rock fragmentation simulation with block smash.
基金the National High Technology Research and Development Program (863) of China(No. 2006AA05Z148)
文摘This paper focuses on the robust stability for time-delay systems of neutral type. A new complete Lyapunov-Krasovskii function (LKF) is developed. Based on this function and discretization, stability conditions in terms of linear matrix inequalities are obtained. A class of time-varying uncertainty of system matrices can be studied by the method.
文摘Discrete software reliability measurement has a proper characteristic for describing a software reliability growth process which depends on a unit of the software fault-detection period, such as the number of test runs, the number of executed test cases. This paper discusses discrete software reliability measurement based on a discretized nonhomogeneous Poisson process (NHPP) model. Especially, we use a bootstrapping method in our discrete software reliability measurement for discussing the statistical inference on parameters and software reliability assessment measures of our model. Finally we show numerical examples of interval estimations based on our bootstrapping method for the several software reliability assessment measures by using actual data.
基金Project supported by the National Natural Science Foundation and of the Ministry of Construction of China
文摘For nonlinear stability problems of discretized conservative systems with multiple parameter variables and multiple state variables, the activation method is put forward, by which activated potential functions and activated equilibrium equations are derived. The activation method is the improvement and enhancement of Liapunov-Schmidt method in elastic stability theory. It is more generalized and more normalized than conventional perturbation methods. The activated potential functions may be transformed into normalized catastrophe potential functions. The activated equilibrium equations may be treated as bifurcation equations. The researches in this paper will motivate the combination of elastic stability theory with catastrophe theory and bifurcation theory
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
文摘In this paper, a discretized SIR model with pulse vaccination and time delay is proposed. We introduce two thresholds R* and R<sub>*</sub>, and further prove that the disease-free periodic solution is globally attractive if R* is less than unit and the disease can invade if R<sub>*</sub> is larger than unit. The numerical simulations not only illustrate the validity of our main results, but also exhibit bifurcation phenomenon. Our result shows that decreasing infection rate can put off the disease outbreak and reduce the number of infected individuals.
基金Supported by the National Natural Science Foundation of China(U2067205)。
文摘We present a systematic study of 6Li elastic scattering and total reaction cross sections at incident energies around the Coulomb barrier within the continuum discretized coupled-channels(CDCC)framework,where 6Li is treated in anα+d two-body model.Collisions with 27Al,64Zn,138Ba,and 208Pa are analyzed.The microscopic optical potentials(MOP)based on Skyrme nucleon-nucleon interaction forαand d are adopted in CDCC calculations and satisfactory agreement with the experimental data is obtained without any adjustment on MOPs.For comparison,αand d global phenomenological optical potentials(GOP)are also used in CDCC analysis and a reduction of no less than 50%on the surface imaginary part of deuteron GOP is required for describing the data.In all cases,the 6Li breakup effect is significant and provides repulsive correction to the folding model potential.The reduction on the surface imaginary part of GOP of deuteron reveals a strong suppression of the reaction probability of deuteron as a component of 6Li when compared with that of a free deuteron.Further investigation is performed by considering the d breakup process equivalently within the dynamic polarization potential approach,and the results show that d behaves in a manner similar to a tightly bound nucleus in 6Li induced reactions.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金supported by the National Priorities Research Program through the Qatar National Research Funda member of Qatar Foundation(Grant No.NPRP 4-1162-1-181)+2 种基金the Natural Science Foundation of China(Grant Nos.6140331361374078&61375102)the Natural Science Foundation Project of Chongqing CSTC(Grant No.cstc2014jcyj A40014)
文摘This paper presents an Euler discretized inertial delayed neuron model, and its bifurcation dynamical behaviors are discussed. By using the associated characteristic model, center manifold theorem and the normal form method, it is shown that the model not only undergoes codimension one(flip, Neimark-Sacker) bifurcation, but also undergoes cusp and resonance bifurcation(1:1 and 1:2) of codimension two. Further, it is found that the parity of delay has some effect on bifurcation behaviors. Finally, some numerical simulations are given to support the analytic results and explore complex dynamics, such as periodic orbits near homoclinic orbits, quasiperiodic orbits, and chaotic orbits.
基金Australian Research Council Linkage Program(LP200301404)for sponsoring this researchthe financial support provided by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology,SKLGP2021K002)National Natural Science Foundation of China(52374101,32111530138).
文摘Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.
基金supports for this research were provided by the National Natural Science Foundation of China(No.12272301,12002278,U1906233)the Guangdong Basic and Applied Basic Research Foundation,China(Nos.2023A1515011970,2024A1515010256)+1 种基金the Dalian City Supports Innovation and Entrepreneurship Projects for High-Level Talents,China(2021RD16)the Key R&D Project of CSCEC,China(No.CSCEC-2020-Z-4).
文摘Fiber-reinforced composites are an ideal material for the lightweight design of aerospace structures. Especially in recent years, with the rapid development of composite additive manufacturing technology, the design optimization of variable stiffness of fiber-reinforced composite laminates has attracted widespread attention from scholars and industry. In these aerospace composite structures, numerous cutout panels and shells serve as access points for maintaining electrical, fuel, and hydraulic systems. The traditional fiber-reinforced composite laminate subtractive drilling manufacturing inevitably faces the problems of interlayer delamination, fiber fracture, and burr of the laminate. Continuous fiber additive manufacturing technology offers the potential for integrated design optimization and manufacturing with high structural performance. Considering the integration of design and manufacturability in continuous fiber additive manufacturing, the paper proposes linear and nonlinear filtering strategies based on the Normal Distribution Fiber Optimization (NDFO) material interpolation scheme to overcome the challenge of discrete fiber optimization results, which are difficult to apply directly to continuous fiber additive manufacturing. With minimizing structural compliance as the objective function, the proposed approach provides a strategy to achieve continuity of discrete fiber paths in the variable stiffness design optimization of composite laminates with regular and irregular holes. In the variable stiffness design optimization model, the number of candidate fiber laying angles in the NDFO material interpolation scheme is considered as design variable. The sensitivity information of structural compliance with respect to the number of candidate fiber laying angles is obtained using the analytical sensitivity analysis method. Based on the proposed variable stiffness design optimization method for complex perforated composite laminates, the numerical examples consider the variable stiffness design optimization of typical non-perforated and perforated composite laminates with circular, square, and irregular holes, and systematically discuss the number of candidate discrete fiber laying angles, discrete fiber continuous filtering strategies, and filter radius on structural compliance, continuity, and manufacturability. The optimized discrete fiber angles of variable stiffness laminates are converted into continuous fiber laying paths using a streamlined process for continuous fiber additive manufacturing. Meanwhile, the optimized non-perforated and perforated MBB beams after discrete fiber continuous treatment, are manufactured using continuous fiber co-extrusion additive manufacturing technology to verify the effectiveness of the variable stiffness fiber optimization framework proposed in this paper.
基金supported by the National Natural Science Foundation of China Youth Science Foundation of China(Grant No.52308388)the Key Project of High-speed Rail Joint Fund of National Natural Science Foundation of China(Grant No.U1934210).
文摘This paper focuses on the use of rotary-percussive drilling for hard rocks.In order to improve efficiency and reduce costs,it is essential to understand how operational parameters,bit wear,and drilling performance are related.A model is presented therein that combines multibody dynamics and discrete element method(DEM)to investigate the influences of operational parameters and bit wear on the rate of penetration and wear characteristics.The model accurately captures the motion of the bit and recreates rock using the cutting sieving result.Field experimental results validate the rod dynamic behavior,rock recreating model,and coupling model in the simulation.The findings indicate that hammer pressure significantly influences the rate of penetration and wear depth of the bit,and there is an optimal range for economical hammer pressure.The wear coefficient has a major effect on the rate of penetration,when wear coefficient is between 1/3 and 2/3.Increasing the wear coefficient can reduce drill bit button pressure and wear depth at the same drill distance.Gauge button loss increases the rate of penetration due to higher pressure on the remaining buttons,which also accelerates destruction of the bit.Furthermore,a more evenly distributed button on the bit enhances the rate of penetration(ROP)when the same number of buttons is lost.
文摘An H-polygon is a polygon with corners in the set of vertices of a tiling of R~2 by regular hexagons of unit edge.It is known that any H-triangle with exactly k interior H-points can have at most 3k+7 boundary H-points.In this note we prove that any convex H-quadrilateral with exactly k interior H-points can have at most 3k+7 boundary H-points.
基金College Students Innovation and Entrepreneurship Project of Guangzhou Railway Polytechnic(2025CXCY015)。
文摘As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness distribution functions of the Bakhtiary dam site and Oskarshamn/Forsmark mountain were fitted using statistical methods.The COMSOL Multiphysics finite element software was utilized to analyze the effects of fracture roughness distribution types and empirical formulas for fracture hydraulic aperture on the seepage field and temperature field of rock masses.The results show that:(1)The fracture roughness at the Bakhtiary dam site and Oskarshamn/Forsmark mountain follows lognormal and normal distributions,respectively;(2)For rock masses with the same expected value and standard deviation of fracture roughness,the outflow from rock masses with lognormal distribution of fracture roughness is significantly larger than that of rock masses with normal distribution of fracture roughness;(3)The fracture hydraulic aperture,outflow,and cold front distance of the Li and Jiang model are significantly larger than those of the Barton model;(4)The outflow,hydraulic pressure distribution,and temperature distribution of the Barton model are more sensitive to the fracture roughness distribution type than those of the Li and Jiang model.
基金supported by the Research Program of Wuhan Building Energy Efficiency Office(grant number 202331).
文摘Precast concrete pavements(PCPs)represent an innovative solution in the construction industry,addressing the need for rapid,intelligent,and low-carbon pavement technologies that significantly reduce construction time and environmental impact.However,the integration of prefabricated technology in pavement surface and base layers lacks systematic classification and understanding.This paper aims to fill this gap by introducing a detailed analysis of discretization and assembly connection technology for cement concrete pavement(CCP)structures.Through a comprehensive review of domestic and international literature,the study classifies prefabricated pavement technology based on discrete assembly structural layers and presents specific conclusions(i)surface layer discrete units are categorized into bottom plates,top plates,plate-rod separated assemblies,and prestressed connections,with optimal material compositions identified to enhance mechanical properties;(ii)base layer discrete units include block-type,plate-type,and beam-type elements,highlighting their contributions to sustainability by incorporating recycled materials(iii)planar assembly connection types are assessed,ranking them by load transfer efficiency,with specific dimensions provided for optimal performance;and(iv)vertical assembly connections are defined by their leveling and sealing layers,suitable for both new constructions and repairs of existing roads.The insights gained from this review not only clarify the distinctions between various structural layers but also provide practical guidelines for enhancing the design and implementation of PCP.This work contributes to advancing sustainable and resilient road construction practices,making it a significant reference for researchers and practitioners in the field.
文摘Ensuring digital media security through robust image watermarking is essential to prevent unauthorized distribution,tampering,and copyright infringement.This study introduces a novel hybrid watermarking framework that integrates Discrete Wavelet Transform(DWT),Redundant Discrete Wavelet Transform(RDWT),and Möbius Transformations(MT),with optimization of transformation parameters achieved via a Genetic Algorithm(GA).By combining frequency and spatial domain techniques,the proposed method significantly enhances both the imper-ceptibility and robustness of watermark embedding.The approach leverages DWT and RDWT for multi-resolution decomposition,enabling watermark insertion in frequency subbands that balance visibility and resistance to attacks.RDWT,in particular,offers shift-invariance,which improves performance under geometric transformations.Möbius transformations are employed for spatial manipulation,providing conformal mapping and spatial dispersion that fortify watermark resilience against rotation,scaling,and translation.The GA dynamically optimizes the Möbius parameters,selecting configurations that maximize robustness metrics such as Peak Signal-to-Noise Ratio(PSNR),Structural Similarity Index Measure(SSIM),Bit Error Rate(BER),and Normalized Cross-Correlation(NCC).Extensive experiments conducted on medical and standard benchmark images demonstrate the efficacy of the proposed RDWT-MT scheme.Results show that PSNR exceeds 68 dB,SSIM approaches 1.0,and BER remains at 0.0000,indicating excellent imperceptibility and perfect watermark recovery.Moreover,the method exhibits exceptional resilience to a wide range of image processing attacks,including Gaussian noise,JPEG compression,histogram equalization,and cropping,achieving NCC values close to or equal to 1.0.Comparative evaluations with state-of-the-art watermarking techniques highlight the superiority of the proposed method in terms of robustness,fidelity,and computational efficiency.The hybrid framework ensures secure,adaptive watermark embedding,making it highly suitable for applications in digital rights management,content authentication,and medical image protection.The integration of spatial and frequency domain features with evolutionary optimization presents a promising direction for future watermarking technologies.
基金supported by the National Natural Science Foundation of China(Grant No.62001391)the Guangdong Basic and Applied Basic Research Foundation(Grant No.2024A1515010308)+1 种基金the Natural Science Basic Research Program of Shaanxi(Grant No.2024JC-YBQN-0464)the Scientific Research Program Funded by Education Department of Shaanxi Provincial Government(Grant No.24JK0559).
文摘This paper proposes a universal impulse-function-based method for extending discrete chaotic maps,enabling flexible construction of multicavity chaotic attractors.The proposed method achieves one-directional(1D)/two-directional(2D)extensions without introducing additional nonlinear terms or altering system stability.Theoretically,the cavity quantity in arbitrary directions is controlled by adjusting impulse levels N,while the amplitude regulation is implemented through modifications to the proportionality parameter r.Theoretical analyses,including Lyapunov exponents(LEs)and bifurcation diagrams,are conducted,confirming that the extended maps retain the intrinsic dynamics of five rational map classes.The field-programmable gate array(FPGA)implementation results are consistent with the numerical simulation results,verifying the correctness of the theoretical analysis.The method enables the expansion of unipolar attractors and enhances entropy metrics,offering a robust framework for applications in secure communication,encryption,and chaos-based technologies.
