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 article,we consider the numerical prediction of the noise emission from a wheelset in laboratory conditions.We focus on the fluid-structure interaction leading to sound emission in the fluid domain by analyzin...In this article,we consider the numerical prediction of the noise emission from a wheelset in laboratory conditions.We focus on the fluid-structure interaction leading to sound emission in the fluid domain by analyzing three different methods to account for acoustic sources.These are a discretized baffled piston using the discrete calculation method(DCM),a closed cylindrical volume using the boundary element method(BEM)and radiating elastic disks in a cubic enclosure solved with the finite element method(FEM).We provide the validation of the baffled piston and the BEM using measurements of the noise emission of a railway wheel by considering ground reflections in the numerical models.Selected space-resolved waveforms are compared with experimental results as well as with a fluid-structure interaction finite element model.The computational advantage of a discretized disk mounted on a baffle and BEM compared to FEM is highlighted,and the baffled pistons limitations caused by a lack of edge radiation effects are investigated.展开更多
The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on th...The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.展开更多
Deep shale gas reservoirs have geological characteristics of high temperature,high pressure,high stress,and inferior ability to pass through fluids.The multi-stage fractured horizontal well is the key to exploiting th...Deep shale gas reservoirs have geological characteristics of high temperature,high pressure,high stress,and inferior ability to pass through fluids.The multi-stage fractured horizontal well is the key to exploiting the deep shale gas reservoir.However,during the production process,the effectiveness of the hydraulic fracture network decreases with the closure of fractures,which accelerates the decline of shale gas production.In this paper,we addressed the problems of unclear fracture closure mechanisms and low accuracy of shale gas production prediction during deep shale gas production.Then we established the fluid—solid—heat coupled model coupling the deformation and fluid flow among the fracture surface,proppant and the shale matrix.When the fluid—solid—heat coupled model was applied to the fracture network,it was well solved by our numerical method named discontinuous discrete fracture method.Compared with the conventional discrete fracture method,the discontinuous discrete fracture method can describe the three-dimensional morphology of the fracture while considering the effect of the change of fracture surface permeation coefficient on the coupled fracture—matrix flow and describing the displacement discontinuity across the fracture.Numerical simulations revealed that the degree of fracture closure increases as the production time proceeds,and the degree of closure of the secondary fractures is higher than that of the primary fractures.Shale creep and proppant embedment both increase the degree of fracture closure.The reduction in fracture surface permeability due to proppant embedment reduces the rate of fluid transfer between matrix and fracture,which has often been overlooked in the past.However,it significantly impacts shale gas production,with calculations showing a 24.7%cumulative three-year yield reduction.This study is helpful to understand the mechanism of hydraulic fracture closure.Therefore,it provides the theoretical guidance for maintaining the long-term effectiveness of hydraulic fractures.展开更多
The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scann...The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.展开更多
Research on discrete memristor-based neural networks has received much attention.However,current research mainly focuses on memristor–based discrete homogeneous neuron networks,while memristor-coupled discrete hetero...Research on discrete memristor-based neural networks has received much attention.However,current research mainly focuses on memristor–based discrete homogeneous neuron networks,while memristor-coupled discrete heterogeneous neuron networks are rarely reported.In this study,a new four-stable discrete locally active memristor is proposed and its nonvolatile and locally active properties are verified by its power-off plot and DC V–I diagram.Based on two-dimensional(2D)discrete Izhikevich neuron and 2D discrete Chialvo neuron,a heterogeneous discrete neuron network is constructed by using the proposed discrete memristor as a coupling synapse connecting the two heterogeneous neurons.Considering the coupling strength as the control parameter,chaotic firing,periodic firing,and hyperchaotic firing patterns are revealed.In particular,multiple coexisting firing patterns are observed,which are induced by different initial values of the memristor.Phase synchronization between the two heterogeneous neurons is discussed and it is found that they can achieve perfect synchronous at large coupling strength.Furthermore,the effect of Gaussian white noise on synchronization behaviors is also explored.We demonstrate that the presence of noise not only leads to the transition of firing patterns,but also achieves the phase synchronization between two heterogeneous neurons under low coupling strength.展开更多
Synaptic crosstalk is a prevalent phenomenon among neuronal synapses,playing a crucial role in the transmission of neural signals.Therefore,considering synaptic crosstalk behavior and investigating the dynamical behav...Synaptic crosstalk is a prevalent phenomenon among neuronal synapses,playing a crucial role in the transmission of neural signals.Therefore,considering synaptic crosstalk behavior and investigating the dynamical behavior of discrete neural networks are highly necessary.In this paper,we propose a heterogeneous discrete neural network(HDNN)consisting of a three-dimensional KTz discrete neuron and a Chialvo discrete neuron.These two neurons are coupled mutually by two discrete memristors and the synaptic crosstalk is considered.The impact of crosstalk strength on the firing behavior of the HDNN is explored through bifurcation diagrams and Lyapunov exponents.It is observed that the HDNN exhibits different coexisting attractors under varying crosstalk strengths.Furthermore,the influence of different crosstalk strengths on the synchronized firing of the HDNN is investigated,revealing a gradual attainment of phase synchronization between the two discrete neurons as the crosstalk strength decreases.展开更多
In the era of big data,the number of images transmitted over the public channel increases exponentially.As a result,it is crucial to devise the efficient and highly secure encryption method to safeguard the sensitive ...In the era of big data,the number of images transmitted over the public channel increases exponentially.As a result,it is crucial to devise the efficient and highly secure encryption method to safeguard the sensitive image.In this paper,an improved sine map(ISM)possessing a larger chaotic region,more complex chaotic behavior and greater unpredictability is proposed and extensively tested.Drawing upon the strengths of ISM,we introduce a lightweight symmetric image encryption cryptosystem in wavelet domain(WDLIC).The WDLIC employs selective encryption to strike a satisfactory balance between security and speed.Initially,only the low-frequency-low-frequency component is chosen to encrypt utilizing classic permutation and diffusion.Then leveraging the statistical properties in wavelet domain,Gaussianization operation which opens the minds of encrypting image information in wavelet domain is first proposed and employed to all sub-bands.Simulations and theoretical analysis demonstrate the high speed and the remarkable effectiveness of WDLIC.展开更多
Internal polyhedral structures of a granular system can be investigated using the Voronoi tessellations.This technique has gained increasing recognition in research of kinetic properties of granular flows.For systems ...Internal polyhedral structures of a granular system can be investigated using the Voronoi tessellations.This technique has gained increasing recognition in research of kinetic properties of granular flows.For systems with mono-sized spherical particles,Voronoi tessellations can be utilized,while radial Voronoi tessellations are necessary for analyzing systems with multi-sized spherical particles.However,research about polyhedral structures of non-spherical particle systems is limited.We utilize the discrete element method to simulate a system of ellipsoidal particles,defined by the equation(x/a)^(2)+(y/1)^(2)+(z/1/a)^(2)=1,where a ranges from 1.1 to 2.0.The system is then dissected by using tangent planes at the contact points,and the geometric quantities of the resulting polyhedra in different shaped systems,such as surface area,volume,number of vertices,number of edges,and number of faces,are calculated.Meanwhile,the longitudinal and transverse wave velocities within the system are calculated with the time-of-flight method.The results demonstrate a strong correlation between the sound velocity of the system and the geometry of the dissected polyhedra.The sound velocity of the system increases with the increase in a,peaking at a=1.3,and then decreases as a continues to increase.The average volume,surface area,number of vertices,number of edges,and number of faces of the polyhedra decrease with the increase in sound velocity.That is,these quantities initially decrease with the increase in a,reaching minima at a=1.3,and then increase with further increase of a.The relationship between sound velocity and the geometric quantities of the dissected polyhedra can serve as a reference for acoustic material design.展开更多
We propose an optical image watermarking scheme based on orbital angular momentum(OAM)holography.Multiple topological charges(TCs,l)of OAM,as multiple cryptographic sub-keys,are embedded into the host image along with...We propose an optical image watermarking scheme based on orbital angular momentum(OAM)holography.Multiple topological charges(TCs,l)of OAM,as multiple cryptographic sub-keys,are embedded into the host image along with the watermark information.Moreover,the Arnold transformation is employed to further enhance the security and the scrambling time(m)is also served as another cryptographic key.The watermark image is embedded into the host image by using the discrete wavelet transformation(DWT)and singular value decomposition(SVD)methods.Importantly,the interference image is utilized to further enhance security.The imperceptibility of our proposed method is analyzed by using the peak signal-to-noise ratio(PSNR)and the histogram of the watermarked host image.To demonstrate robustness,a series of attack tests,including Gaussian noise,Poisson noise,salt-and-pepper noise,JPEG compression,Gaussian lowpass filtering,cropping,and rotation,are conducted.The experimental results show that our proposed method has advanced security,imperceptibility,and robustness,making it a promising option for optical image watermarking applications.展开更多
We study the Nadaraya-Watson estimators for the drift function of two-sided reflected stochastic differential equations.The estimates,based on either the continuously observed process or the discretely observed proces...We study the Nadaraya-Watson estimators for the drift function of two-sided reflected stochastic differential equations.The estimates,based on either the continuously observed process or the discretely observed process,are considered.Under certain conditions,we prove the strong consistency and the asymptotic normality of the two estimators.Our method is also suitable for one-sided reflected stochastic differential equations.Simulation results demonstrate that the performance of our estimator is superior to that of the estimator proposed by Cholaquidis et al.(Stat Sin,2021,31:29-51).Several real data sets of the currency exchange rate are used to illustrate our proposed methodology.展开更多
This study introduces a novel method integrating CO_(2)flooding with radial borehole fracturing for enhanced oil recovery and CO_(2)underground storage,a solution to the limited vertical stimulation reservoir volume i...This study introduces a novel method integrating CO_(2)flooding with radial borehole fracturing for enhanced oil recovery and CO_(2)underground storage,a solution to the limited vertical stimulation reservoir volume in horizontal well fracturing.A numerical model is established to investigate the production rate,reservoir pressure field,and CO_(2)saturation distribution corresponding to changing time of CO_(2)flooding with radial borehole fracturing.A sensitivity analysis on the influence of CO_(2)injection location,layer spacing,pressure difference,borehole number,and hydraulic fractures on oil production and CO_(2)storage is conducted.The CO_(2)flooding process is divided into four stages.Reductions in layer spacing will significantly improve oil production rate and gas storage capacity.However,serious gas channeling can occur when the spacing is lower than 20 m.Increasing the pressure difference between the producer and injector,the borehole number,the hydraulic fracture height,and the fracture width can also increase the oil production rate and gas storage rate.Sensitivity analysis shows that layer spacing and fracture height greatly influence gas storage and oil production.Research outcomes are expected to provide a theoretical basis for the efficient development of shale oil reservoirs in the vertical direction.展开更多
Acoustic holograms can recover wavefront stored acoustic field information and produce high-fidelity complex acoustic fields. Benefiting from the huge spatial information that traditional acoustic elements cannot matc...Acoustic holograms can recover wavefront stored acoustic field information and produce high-fidelity complex acoustic fields. Benefiting from the huge spatial information that traditional acoustic elements cannot match, acoustic holograms pursue the realization of high-resolution complex acoustic fields and gradually tend to high-frequency ultrasound applications. However, conventional continuous phase holograms are limited by three-dimensional(3D) printing size, and the presence of unavoidable small printing errors makes it difficult to achieve acoustic field reconstruction at high frequency accuracy. Here, we present an optimized discrete multi-step phase hologram. It can ensure the reconstruction quality of image with high robustness, and properly lower the requirement for the 3D printing accuracy. Meanwhile, the concept of reconstruction similarity is proposed to refine a measure of acoustic field quality. In addition, the realized complex acoustic field at 20 MHz promotes the application of acoustic holograms at high frequencies and provides a new way to generate high-fidelity acoustic fields.展开更多
Discrete feedback control was designed to stabilize an unstable hybrid neutral stochastic differential delay system(HNSDDS) under a highly nonlinear constraint in the H_∞ and exponential forms.Nevertheless,the existi...Discrete feedback control was designed to stabilize an unstable hybrid neutral stochastic differential delay system(HNSDDS) under a highly nonlinear constraint in the H_∞ and exponential forms.Nevertheless,the existing work just adapted to autonomous cases,and the obtained results were mainly on exponential stabilization.In comparison with autonomous cases,non-autonomous systems are of great interest and represent an important challenge.Accordingly,discrete feedback control has here been adjusted with a time factor to stabilize an unstable non-autonomous HNSDDS,in which new Lyapunov-Krasovskii functionals and some novel technologies are adopted.It should be noted,in particular,that the stabilization can be achieved not only in the routine H_∞ and exponential forms,but also the polynomial form and even a general form.展开更多
In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. ...In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.展开更多
In order to avoid the complexity of Gaussian modulation and the problem that the traditional point-to-point communication DM-CVQKD protocol cannot meet the demand for multi-user key sharing at the same time, we propos...In order to avoid the complexity of Gaussian modulation and the problem that the traditional point-to-point communication DM-CVQKD protocol cannot meet the demand for multi-user key sharing at the same time, we propose a multi-ring discrete modulation continuous variable quantum key sharing scheme(MR-DM-CVQSS). In this paper, we primarily compare single-ring and multi-ring M-symbol amplitude and phase-shift keying modulations. We analyze their asymptotic key rates against collective attacks and consider the security key rates under finite-size effects. Leveraging the characteristics of discrete modulation, we improve the quantum secret sharing scheme. Non-dealer participants only require simple phase shifters to complete quantum secret sharing. We also provide the general design of the MR-DM-CVQSS protocol.We conduct a comprehensive analysis of the improved protocol's performance, confirming that the enhancement through multi-ring M-PSK allows for longer-distance quantum key distribution. Additionally, it reduces the deployment complexity of the system, thereby increasing the practical value.展开更多
Extensive numerical simulations and scaling analysis are performed to investigate competitive growth between the linear and nonlinear stochastic dynamic growth systems, which belong to the Edwards–Wilkinson(EW) and K...Extensive numerical simulations and scaling analysis are performed to investigate competitive growth between the linear and nonlinear stochastic dynamic growth systems, which belong to the Edwards–Wilkinson(EW) and Kardar–Parisi–Zhang(KPZ) universality classes, respectively. The linear growth systems include the EW equation and the model of random deposition with surface relaxation(RDSR), the nonlinear growth systems involve the KPZ equation and typical discrete models including ballistic deposition(BD), etching, and restricted solid on solid(RSOS). The scaling exponents are obtained in both the(1 + 1)-and(2 + 1)-dimensional competitive growth with the nonlinear growth probability p and the linear proportion 1-p. Our results show that, when p changes from 0 to 1, there exist non-trivial crossover effects from EW to KPZ universality classes based on different competitive growth rules. Furthermore, the growth rate and the porosity are also estimated within various linear and nonlinear growths of cooperation and competition.展开更多
The cross-level and twist irregularities are the most dangerous irregularity types that could cause wheel unloading with the risk of derailments and additional maintenance expenses.However,the mechanism of the irregul...The cross-level and twist irregularities are the most dangerous irregularity types that could cause wheel unloading with the risk of derailments and additional maintenance expenses.However,the mechanism of the irregularities initiation and development is unclear.The motivation of the present study was the previous experimental studies on the application of wide sleepers in the ballasted track.The long-term track geometry measurements with wide sleepers show an enormous reduction of the vertical longitudinal irregularities compared to the conventional track.However,wide sleepers had higher twist and cross-section level irregularities.The present paper aims to explain the phenomenon by discrete element method(DEM)modeling the development process of sleeper inhomogeneous support at cross-level depending on the sleeper form.The DEM simulations show that the maximal settlement intensity is up to 3.5 times lower for a wide sleeper in comparison with the conventional one.Nevertheless,the cross-level differential settlements are almost the same for both sleepers.The particle loading distribution after all loading cycles is concentrated on the smaller area,up to the half sleeper length,with fully unloaded zones under sleeper ends.Ballast flow limitation under the central part of the sleeper could improve the resilience of wide sleepers to the development of cross-level irregularities.The mechanism of initiation of the cross-level irregularity is proposed,which assumes the loss of sleeper support under sleeper ends.The further growth of inhomogeneous settlements along the sleeper is assumed as a result of the interaction of two processes:ballast flow due to dynamic impact during void closing and on the other side high pressure due to the concentration of the pressure under the middle part of the sleeper.The DEM simulation results support the assumption of the mechanism and agree with the experimental studies.展开更多
A novel color image encryption scheme is developed to enhance the security of encryption without increasing the complexity. Firstly, the plain color image is decomposed into three grayscale plain images, which are con...A novel color image encryption scheme is developed to enhance the security of encryption without increasing the complexity. Firstly, the plain color image is decomposed into three grayscale plain images, which are converted into the frequency domain coefficient matrices(FDCM) with discrete cosine transform(DCT) operation. After that, a twodimensional(2D) coupled chaotic system is developed and used to generate one group of embedded matrices and another group of encryption matrices, respectively. The embedded matrices are integrated with the FDCM to fulfill the frequency domain encryption, and then the inverse DCT processing is implemented to recover the spatial domain signal. Eventually,under the function of the encryption matrices and the proposed diagonal scrambling algorithm, the final color ciphertext is obtained. The experimental results show that the proposed method can not only ensure efficient encryption but also satisfy various sizes of image encryption. Besides, it has better performance than other similar techniques in statistical feature analysis, such as key space, key sensitivity, anti-differential attack, information entropy, noise attack, etc.展开更多
Laying the under-sleeper pad(USP)is one of the effective measures commonly used to delay ballast degradation and reduce maintenance workload.To explore the impact of application of the USP on the dynamic and static me...Laying the under-sleeper pad(USP)is one of the effective measures commonly used to delay ballast degradation and reduce maintenance workload.To explore the impact of application of the USP on the dynamic and static mechanical behavior of the ballast track in the heavy-haul railway system,numerical simulation models of the ballast bed with USP and without USP are presented in this paper by using the discrete element method(DEM)-multi-flexible body dynamic(MFBD)coupling analysis method.The ballast bed support stiffness test and dynamic displacement tests were carried out on the actual operation of a heavy-haul railway line to verify the validity of the models.The results show that using the USP results in a 43.01%reduction in the ballast bed support stiffness and achieves a more uniform distribution of track loads on the sleepers.It effectively reduces the load borne by the sleeper directly under the wheel load,with a 7.89%reduction in the pressure on the sleeper.Furthermore,the laying of the USP changes the lateral resistance sharing ratio of the ballast bed,significantly reducing the stress level of the ballast bed under train loads,with an average stress reduction of 42.19 kPa.It also reduces the plastic displacement of ballast particles and lowers the peak value of rotational angular velocity by about 50%to 70%,which is conducive to slowing down ballast bed settlement deformation and reducing maintenance costs.In summary,laying the USP has a potential value in enhancing the stability and extending the lifespan of the ballast bed in heavy-haul railway systems.展开更多
文摘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.
基金The project was commissioned and supported by the funding of the Federal Office of Environment(No.1337000438).
文摘In this article,we consider the numerical prediction of the noise emission from a wheelset in laboratory conditions.We focus on the fluid-structure interaction leading to sound emission in the fluid domain by analyzing three different methods to account for acoustic sources.These are a discretized baffled piston using the discrete calculation method(DCM),a closed cylindrical volume using the boundary element method(BEM)and radiating elastic disks in a cubic enclosure solved with the finite element method(FEM).We provide the validation of the baffled piston and the BEM using measurements of the noise emission of a railway wheel by considering ground reflections in the numerical models.Selected space-resolved waveforms are compared with experimental results as well as with a fluid-structure interaction finite element model.The computational advantage of a discretized disk mounted on a baffle and BEM compared to FEM is highlighted,and the baffled pistons limitations caused by a lack of edge radiation effects are investigated.
基金This work was supported by the National Key R&D Program‘Transportation Infrastructure’project(No.2022YFB2603400).
文摘The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.
基金the supports provided by China University of Petroleum,Beijing(Grand No.ZX20230042)the National Natural Science Foundation of China(Grand No.52334001and Grand No.51904314)。
文摘Deep shale gas reservoirs have geological characteristics of high temperature,high pressure,high stress,and inferior ability to pass through fluids.The multi-stage fractured horizontal well is the key to exploiting the deep shale gas reservoir.However,during the production process,the effectiveness of the hydraulic fracture network decreases with the closure of fractures,which accelerates the decline of shale gas production.In this paper,we addressed the problems of unclear fracture closure mechanisms and low accuracy of shale gas production prediction during deep shale gas production.Then we established the fluid—solid—heat coupled model coupling the deformation and fluid flow among the fracture surface,proppant and the shale matrix.When the fluid—solid—heat coupled model was applied to the fracture network,it was well solved by our numerical method named discontinuous discrete fracture method.Compared with the conventional discrete fracture method,the discontinuous discrete fracture method can describe the three-dimensional morphology of the fracture while considering the effect of the change of fracture surface permeation coefficient on the coupled fracture—matrix flow and describing the displacement discontinuity across the fracture.Numerical simulations revealed that the degree of fracture closure increases as the production time proceeds,and the degree of closure of the secondary fractures is higher than that of the primary fractures.Shale creep and proppant embedment both increase the degree of fracture closure.The reduction in fracture surface permeability due to proppant embedment reduces the rate of fluid transfer between matrix and fracture,which has often been overlooked in the past.However,it significantly impacts shale gas production,with calculations showing a 24.7%cumulative three-year yield reduction.This study is helpful to understand the mechanism of hydraulic fracture closure.Therefore,it provides the theoretical guidance for maintaining the long-term effectiveness of hydraulic fractures.
基金supported by the National Natural Science Foundation of China under Grants Nos.52165013 and 51565021.
文摘The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.
基金Project supported by the National Natural Science Foundations of China(Grant Nos.62171401 and 62071411).
文摘Research on discrete memristor-based neural networks has received much attention.However,current research mainly focuses on memristor–based discrete homogeneous neuron networks,while memristor-coupled discrete heterogeneous neuron networks are rarely reported.In this study,a new four-stable discrete locally active memristor is proposed and its nonvolatile and locally active properties are verified by its power-off plot and DC V–I diagram.Based on two-dimensional(2D)discrete Izhikevich neuron and 2D discrete Chialvo neuron,a heterogeneous discrete neuron network is constructed by using the proposed discrete memristor as a coupling synapse connecting the two heterogeneous neurons.Considering the coupling strength as the control parameter,chaotic firing,periodic firing,and hyperchaotic firing patterns are revealed.In particular,multiple coexisting firing patterns are observed,which are induced by different initial values of the memristor.Phase synchronization between the two heterogeneous neurons is discussed and it is found that they can achieve perfect synchronous at large coupling strength.Furthermore,the effect of Gaussian white noise on synchronization behaviors is also explored.We demonstrate that the presence of noise not only leads to the transition of firing patterns,but also achieves the phase synchronization between two heterogeneous neurons under low coupling strength.
基金Project supported by the Key Projects of Hunan Provincial Department of Education(Grant No.23A0133)the Natural Science Foundation of Hunan Province(Grant No.2022JJ30572)the National Natural Science Foundations of China(Grant No.62171401).
文摘Synaptic crosstalk is a prevalent phenomenon among neuronal synapses,playing a crucial role in the transmission of neural signals.Therefore,considering synaptic crosstalk behavior and investigating the dynamical behavior of discrete neural networks are highly necessary.In this paper,we propose a heterogeneous discrete neural network(HDNN)consisting of a three-dimensional KTz discrete neuron and a Chialvo discrete neuron.These two neurons are coupled mutually by two discrete memristors and the synaptic crosstalk is considered.The impact of crosstalk strength on the firing behavior of the HDNN is explored through bifurcation diagrams and Lyapunov exponents.It is observed that the HDNN exhibits different coexisting attractors under varying crosstalk strengths.Furthermore,the influence of different crosstalk strengths on the synchronized firing of the HDNN is investigated,revealing a gradual attainment of phase synchronization between the two discrete neurons as the crosstalk strength decreases.
基金Project supported by the Key Area Research and Development Program of Guangdong Province,China(Grant No.2022B0701180001)the National Natural Science Foundation of China(Grant No.61801127)+1 种基金the Science Technology Planning Project of Guangdong Province,China(Grant Nos.2019B010140002 and 2020B111110002)the Guangdong–Hong Kong–Macao Joint Innovation Field Project(Grant No.2021A0505080006).
文摘In the era of big data,the number of images transmitted over the public channel increases exponentially.As a result,it is crucial to devise the efficient and highly secure encryption method to safeguard the sensitive image.In this paper,an improved sine map(ISM)possessing a larger chaotic region,more complex chaotic behavior and greater unpredictability is proposed and extensively tested.Drawing upon the strengths of ISM,we introduce a lightweight symmetric image encryption cryptosystem in wavelet domain(WDLIC).The WDLIC employs selective encryption to strike a satisfactory balance between security and speed.Initially,only the low-frequency-low-frequency component is chosen to encrypt utilizing classic permutation and diffusion.Then leveraging the statistical properties in wavelet domain,Gaussianization operation which opens the minds of encrypting image information in wavelet domain is first proposed and employed to all sub-bands.Simulations and theoretical analysis demonstrate the high speed and the remarkable effectiveness of WDLIC.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12262005,11962003,and 11602062)the Postgraduate Education Reform and Quality Improvement Project of Henan Province(Grant No.YJS2024AL138)the Graduate Education Reform Project of Henan Province(Grant No.2023SJGLX096Y).
文摘Internal polyhedral structures of a granular system can be investigated using the Voronoi tessellations.This technique has gained increasing recognition in research of kinetic properties of granular flows.For systems with mono-sized spherical particles,Voronoi tessellations can be utilized,while radial Voronoi tessellations are necessary for analyzing systems with multi-sized spherical particles.However,research about polyhedral structures of non-spherical particle systems is limited.We utilize the discrete element method to simulate a system of ellipsoidal particles,defined by the equation(x/a)^(2)+(y/1)^(2)+(z/1/a)^(2)=1,where a ranges from 1.1 to 2.0.The system is then dissected by using tangent planes at the contact points,and the geometric quantities of the resulting polyhedra in different shaped systems,such as surface area,volume,number of vertices,number of edges,and number of faces,are calculated.Meanwhile,the longitudinal and transverse wave velocities within the system are calculated with the time-of-flight method.The results demonstrate a strong correlation between the sound velocity of the system and the geometry of the dissected polyhedra.The sound velocity of the system increases with the increase in a,peaking at a=1.3,and then decreases as a continues to increase.The average volume,surface area,number of vertices,number of edges,and number of faces of the polyhedra decrease with the increase in sound velocity.That is,these quantities initially decrease with the increase in a,reaching minima at a=1.3,and then increase with further increase of a.The relationship between sound velocity and the geometric quantities of the dissected polyhedra can serve as a reference for acoustic material design.
基金Project supported by the National Natural Science Foundation of China(Grant No.62375140)the Natural Science Foundation of Suqian,Jiangsu Province,China(Grant No.S202108)+1 种基金the Open Research Fund of the National Laboratory of Solid State Microstructures(Grant No.M36055)the Postgraduate Research&Practice Innovation Program of Jiangsu Province,China(Grant No.KYCX21-0745)。
文摘We propose an optical image watermarking scheme based on orbital angular momentum(OAM)holography.Multiple topological charges(TCs,l)of OAM,as multiple cryptographic sub-keys,are embedded into the host image along with the watermark information.Moreover,the Arnold transformation is employed to further enhance the security and the scrambling time(m)is also served as another cryptographic key.The watermark image is embedded into the host image by using the discrete wavelet transformation(DWT)and singular value decomposition(SVD)methods.Importantly,the interference image is utilized to further enhance security.The imperceptibility of our proposed method is analyzed by using the peak signal-to-noise ratio(PSNR)and the histogram of the watermarked host image.To demonstrate robustness,a series of attack tests,including Gaussian noise,Poisson noise,salt-and-pepper noise,JPEG compression,Gaussian lowpass filtering,cropping,and rotation,are conducted.The experimental results show that our proposed method has advanced security,imperceptibility,and robustness,making it a promising option for optical image watermarking applications.
基金partially supported by the National Natural Science Foundation of China(11871244)the Fundamental Research Funds for the Central Universities,JLU。
文摘We study the Nadaraya-Watson estimators for the drift function of two-sided reflected stochastic differential equations.The estimates,based on either the continuously observed process or the discretely observed process,are considered.Under certain conditions,we prove the strong consistency and the asymptotic normality of the two estimators.Our method is also suitable for one-sided reflected stochastic differential equations.Simulation results demonstrate that the performance of our estimator is superior to that of the estimator proposed by Cholaquidis et al.(Stat Sin,2021,31:29-51).Several real data sets of the currency exchange rate are used to illustrate our proposed methodology.
基金This study has been funded by the National Science Fund for Distinguished Young Scholars(No.52204063)Science Foundation of China University of Petroleum,Beijing(No.2462023BJRC025).Moreover,we would like to express our heartfelt appreciation to the Computational Geosciences group in the Department of Mathematics and Cybernetics at SINTEF Digital for developing and providing the free open-source MATLAB Reservoir Simulation Toolbox(MRST)used in this research.
文摘This study introduces a novel method integrating CO_(2)flooding with radial borehole fracturing for enhanced oil recovery and CO_(2)underground storage,a solution to the limited vertical stimulation reservoir volume in horizontal well fracturing.A numerical model is established to investigate the production rate,reservoir pressure field,and CO_(2)saturation distribution corresponding to changing time of CO_(2)flooding with radial borehole fracturing.A sensitivity analysis on the influence of CO_(2)injection location,layer spacing,pressure difference,borehole number,and hydraulic fractures on oil production and CO_(2)storage is conducted.The CO_(2)flooding process is divided into four stages.Reductions in layer spacing will significantly improve oil production rate and gas storage capacity.However,serious gas channeling can occur when the spacing is lower than 20 m.Increasing the pressure difference between the producer and injector,the borehole number,the hydraulic fracture height,and the fracture width can also increase the oil production rate and gas storage rate.Sensitivity analysis shows that layer spacing and fracture height greatly influence gas storage and oil production.Research outcomes are expected to provide a theoretical basis for the efficient development of shale oil reservoirs in the vertical direction.
基金Project supported by the China Postdoctoral Science Foundation (Grant No.2023M732745)the National Natural Science Foundations of China (Grant Nos.61974110 and 62104177)+1 种基金the Fundamental Research Funds for the Central Universities,China (Grant Nos.QTZX23022 and JBF211103)the Cooperation Program of XDU– Chongqing IC Innovation Research Institute (Grant No.CQ IRI-2022CXY-Z07)。
文摘Acoustic holograms can recover wavefront stored acoustic field information and produce high-fidelity complex acoustic fields. Benefiting from the huge spatial information that traditional acoustic elements cannot match, acoustic holograms pursue the realization of high-resolution complex acoustic fields and gradually tend to high-frequency ultrasound applications. However, conventional continuous phase holograms are limited by three-dimensional(3D) printing size, and the presence of unavoidable small printing errors makes it difficult to achieve acoustic field reconstruction at high frequency accuracy. Here, we present an optimized discrete multi-step phase hologram. It can ensure the reconstruction quality of image with high robustness, and properly lower the requirement for the 3D printing accuracy. Meanwhile, the concept of reconstruction similarity is proposed to refine a measure of acoustic field quality. In addition, the realized complex acoustic field at 20 MHz promotes the application of acoustic holograms at high frequencies and provides a new way to generate high-fidelity acoustic fields.
基金supported by the National Natural Science Foundation of China(61833005)the Humanities and Social Science Fund of Ministry of Education of China(23YJAZH031)+1 种基金the Natural Science Foundation of Hebei Province of China(A2023209002,A2019209005)the Tangshan Science and Technology Bureau Program of Hebei Province of China(19130222g)。
文摘Discrete feedback control was designed to stabilize an unstable hybrid neutral stochastic differential delay system(HNSDDS) under a highly nonlinear constraint in the H_∞ and exponential forms.Nevertheless,the existing work just adapted to autonomous cases,and the obtained results were mainly on exponential stabilization.In comparison with autonomous cases,non-autonomous systems are of great interest and represent an important challenge.Accordingly,discrete feedback control has here been adjusted with a time factor to stabilize an unstable non-autonomous HNSDDS,in which new Lyapunov-Krasovskii functionals and some novel technologies are adopted.It should be noted,in particular,that the stabilization can be achieved not only in the routine H_∞ and exponential forms,but also the polynomial form and even a general form.
文摘In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61971348 and 61201194)。
文摘In order to avoid the complexity of Gaussian modulation and the problem that the traditional point-to-point communication DM-CVQKD protocol cannot meet the demand for multi-user key sharing at the same time, we propose a multi-ring discrete modulation continuous variable quantum key sharing scheme(MR-DM-CVQSS). In this paper, we primarily compare single-ring and multi-ring M-symbol amplitude and phase-shift keying modulations. We analyze their asymptotic key rates against collective attacks and consider the security key rates under finite-size effects. Leveraging the characteristics of discrete modulation, we improve the quantum secret sharing scheme. Non-dealer participants only require simple phase shifters to complete quantum secret sharing. We also provide the general design of the MR-DM-CVQSS protocol.We conduct a comprehensive analysis of the improved protocol's performance, confirming that the enhancement through multi-ring M-PSK allows for longer-distance quantum key distribution. Additionally, it reduces the deployment complexity of the system, thereby increasing the practical value.
基金supported by Undergraduate Training Program for Innovation and Entrepreneurship of China University of Mining and Technology (CUMT)(Grant No. 202110290059Z)Fundamental Research Funds for the Central Universities of CUMT (Grant No. 2020ZDPYMS33)。
文摘Extensive numerical simulations and scaling analysis are performed to investigate competitive growth between the linear and nonlinear stochastic dynamic growth systems, which belong to the Edwards–Wilkinson(EW) and Kardar–Parisi–Zhang(KPZ) universality classes, respectively. The linear growth systems include the EW equation and the model of random deposition with surface relaxation(RDSR), the nonlinear growth systems involve the KPZ equation and typical discrete models including ballistic deposition(BD), etching, and restricted solid on solid(RSOS). The scaling exponents are obtained in both the(1 + 1)-and(2 + 1)-dimensional competitive growth with the nonlinear growth probability p and the linear proportion 1-p. Our results show that, when p changes from 0 to 1, there exist non-trivial crossover effects from EW to KPZ universality classes based on different competitive growth rules. Furthermore, the growth rate and the porosity are also estimated within various linear and nonlinear growths of cooperation and competition.
文摘The cross-level and twist irregularities are the most dangerous irregularity types that could cause wheel unloading with the risk of derailments and additional maintenance expenses.However,the mechanism of the irregularities initiation and development is unclear.The motivation of the present study was the previous experimental studies on the application of wide sleepers in the ballasted track.The long-term track geometry measurements with wide sleepers show an enormous reduction of the vertical longitudinal irregularities compared to the conventional track.However,wide sleepers had higher twist and cross-section level irregularities.The present paper aims to explain the phenomenon by discrete element method(DEM)modeling the development process of sleeper inhomogeneous support at cross-level depending on the sleeper form.The DEM simulations show that the maximal settlement intensity is up to 3.5 times lower for a wide sleeper in comparison with the conventional one.Nevertheless,the cross-level differential settlements are almost the same for both sleepers.The particle loading distribution after all loading cycles is concentrated on the smaller area,up to the half sleeper length,with fully unloaded zones under sleeper ends.Ballast flow limitation under the central part of the sleeper could improve the resilience of wide sleepers to the development of cross-level irregularities.The mechanism of initiation of the cross-level irregularity is proposed,which assumes the loss of sleeper support under sleeper ends.The further growth of inhomogeneous settlements along the sleeper is assumed as a result of the interaction of two processes:ballast flow due to dynamic impact during void closing and on the other side high pressure due to the concentration of the pressure under the middle part of the sleeper.The DEM simulation results support the assumption of the mechanism and agree with the experimental studies.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.62105004 and 52174141)the College Student Innovation and Entrepreneurship Fund Project(Grant No.202210361053)+1 种基金Anhui Mining Machinery and Electrical Equipment Coordination Innovation Center,Anhui University of Science&Technology(Grant No.KSJD202304)the Anhui Province Digital Agricultural Engineering Technology Research Center Open Project(Grant No.AHSZNYGC-ZXKF021)。
文摘A novel color image encryption scheme is developed to enhance the security of encryption without increasing the complexity. Firstly, the plain color image is decomposed into three grayscale plain images, which are converted into the frequency domain coefficient matrices(FDCM) with discrete cosine transform(DCT) operation. After that, a twodimensional(2D) coupled chaotic system is developed and used to generate one group of embedded matrices and another group of encryption matrices, respectively. The embedded matrices are integrated with the FDCM to fulfill the frequency domain encryption, and then the inverse DCT processing is implemented to recover the spatial domain signal. Eventually,under the function of the encryption matrices and the proposed diagonal scrambling algorithm, the final color ciphertext is obtained. The experimental results show that the proposed method can not only ensure efficient encryption but also satisfy various sizes of image encryption. Besides, it has better performance than other similar techniques in statistical feature analysis, such as key space, key sensitivity, anti-differential attack, information entropy, noise attack, etc.
基金the project supported by the National Natural Science Foundation of China(Grant No.52372425)the Fundamental Research Funds for the Central Universities(Science and technology leading talent team project)(Grant No.2022JBXT010).
文摘Laying the under-sleeper pad(USP)is one of the effective measures commonly used to delay ballast degradation and reduce maintenance workload.To explore the impact of application of the USP on the dynamic and static mechanical behavior of the ballast track in the heavy-haul railway system,numerical simulation models of the ballast bed with USP and without USP are presented in this paper by using the discrete element method(DEM)-multi-flexible body dynamic(MFBD)coupling analysis method.The ballast bed support stiffness test and dynamic displacement tests were carried out on the actual operation of a heavy-haul railway line to verify the validity of the models.The results show that using the USP results in a 43.01%reduction in the ballast bed support stiffness and achieves a more uniform distribution of track loads on the sleepers.It effectively reduces the load borne by the sleeper directly under the wheel load,with a 7.89%reduction in the pressure on the sleeper.Furthermore,the laying of the USP changes the lateral resistance sharing ratio of the ballast bed,significantly reducing the stress level of the ballast bed under train loads,with an average stress reduction of 42.19 kPa.It also reduces the plastic displacement of ballast particles and lowers the peak value of rotational angular velocity by about 50%to 70%,which is conducive to slowing down ballast bed settlement deformation and reducing maintenance costs.In summary,laying the USP has a potential value in enhancing the stability and extending the lifespan of the ballast bed in heavy-haul railway systems.
