The influence of inter-band rock on rib spalling in longwall panel with large mining height

In order to improve rib stability,failure criteria and instability mode of a thick coal seam with inter-band rock layer are analysed in this study.A three-dimensional mechanical model is established for the rib by considering the rock layer.A safety factor is defined foy the rib,and it is observed that the safety factor exhibits a positive correlation with the thickness and strength of the inter-band rock.A calculation method for determining critical parameters of the rock layer is presented to ensure the rib stability.It is revealed that incomplete propagation of the fracture at the hard rock constitutes a fundamental prerequisite for ensuring the rib stability.The influence of the position of the inter-band rock in the coal seam on failure mechanism of the rib was thoroughly investigated by developing a series of physical models for the rib at the face area.The best position for the inter-band rock in the coal seam is at a height of 1.5 m away from the roof line,which tends to provide a good stability state for the rib.For different inter-band rock positions,two ways of controlling rib by increasing supports stiffness and flexible grouting reinforcement are proposed.

Relative Impacts of Sea Ice Loss and Atmospheric Internal Variability on the Winter Arctic to East Asian Surface Air Temperature Based on Large-Ensemble Simulations with NorESM2

To quantify the relative contributions of Arctic sea ice and unforced atmospheric internal variability to the “warm Arctic, cold East Asia”(WACE) teleconnection, this study analyses three sets of large-ensemble simulations carried out by the Norwegian Earth System Model with a coupled atmosphere–land surface model, forced by seasonal sea ice conditions from preindustrial, present-day, and future periods. Each ensemble member within the same set uses the same forcing but with small perturbations to the atmospheric initial state. Hence, the difference between the present-day(or future) ensemble mean and the preindustrial ensemble mean provides the ice-loss-induced response, while the difference of the individual members within the present-day(or future) set is the effect of atmospheric internal variability. Results indicate that both present-day and future sea ice loss can force a negative phase of the Arctic Oscillation with a WACE pattern in winter. The magnitude of ice-induced Arctic warming is over four(ten) times larger than the ice-induced East Asian cooling in the present-day(future) experiment;the latter having a magnitude that is about 30% of the observed cooling. Sea ice loss contributes about 60%(80%) to the Arctic winter warming in the present-day(future) experiment. Atmospheric internal variability can also induce a WACE pattern with comparable magnitudes between the Arctic and East Asia. Ice-lossinduced East Asian cooling can easily be masked by atmospheric internal variability effects because random atmospheric internal variability may induce a larger magnitude warming. The observed WACE pattern occurs as a result of both Arctic sea ice loss and atmospheric internal variability, with the former dominating Arctic warming and the latter dominating East Asian cooling.

Breeding design in wheat by combining the QTL information in a GWAS panel with a general genetic map and computer simulation

A large amount of genome-wide association study(GWAS)panels together with quantitative-trait locus(QTL)information associated with breeding-targeted traits have been described in wheat(Triticum aestivum L.).However,the application of mapping results from a GWAS panel to conventional wheat breeding remains a challenge.In this study,we first report a general genetic map which was constructed from 44 published linkage maps.It permits the estimation of genetic distances between any two genetic loci with physical map positions,thereby unifying the linkage relationships between QTL,genes,and genomic markers from multiple genetic populations.Second,we describe QTL mapping in a wheat GWAS panel of 688 accessions,identifying 77 QTL associated with 12 yield and grain-quality traits.Because these QTL have known physical map positions,they could be mapped onto the general map.Finally,we present a design approach to wheat breeding by using known QTL information and computer simulation.Potential crosses between parents in the GWAS panel may be evaluated by the relative frequency of the target genotype,trait correlations in simulated progeny populations,and genetic gain of selected progenies.It is possible to simultaneously improve yield and grain quality by suitable parental selection,progeny population size,and progeny selection scheme.Applying the design approach will allow identifying the most promising crosses and selection schemes in advance of the field experiment,increasing predictability and efficiency in wheat breeding.

MicroMagnetic.jl:A Julia package for micromagnetic and atomistic simulations with GPU support

MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.

Generative optimization of bistable plates with deep learning

Bistate plates have found extensive applications in the domains of smart structures and energy harvesting devices.Most bistable curved plates are characterized by a constant thickness profile.Regrettably,due to the inherent complexity of this problem,relatively little attention has been devoted to this area.In this study,we demonstrate how deep learning can facilitate the discovery of novel plate profiles that cater to multiple objectives,including maximizing stiffness,forward snapping force,and backward snapping force.Our proposed approach is distinguished by its efficiency in terms of low computational energy consumption and high effectiveness.It holds promise for future applications in the design and optimization of multistable structures with diverse objectives,addressing the requirements of various fields.

Numerical Simulation and Entropy Production Analysis of Centrifugal Pump with Various Viscosity

The fluid’s viscosity significantly affects the performance of a centrifugal pump.The entropy production method and leakage are employed to analyze the performance changes under various viscosities by numerical simulation and validated by experiments.The results showed that increasing viscosity reduces both the pump head and efficiency.In addition,the optimal operating point shifts to the left.Leakage is influenced by vortex distribution in the front chamber and boundary layer thickness in wear-ring clearance,leading to an initial increase and subsequent decrease in leakage with increasing viscosity.The total entropy production Spro,Total inside the pump rises with increasing viscosity.The different mechanisms dominate under varying conditions:Turbulent dissipation dominates at low viscosity.Under high-viscosity conditions,energy loss is primarily caused by direct dissipation Spro,D and wall entropy production Spro,W.This study provides a deeper and more objective understanding of the energy characteristics of centrifugal pumps handling fluids of various viscosity,potentially aiding in optimizing pump design and improving energy conversion efficiency.

Experimental and numerical study of buckling-restrained braces configured with out-of-plane eccentricity under cyclic loading

This study focuses on variations in the hysteretic behavior of buckling-restrained braces(BRBs)configured with or without out-of-plane eccentricity under cyclic loading.Quasi-static experiments and numerical simulations were carried out on concentrically and eccentrically loaded BRB specimens to investigate the mechanical properties,energy dissipation performance,stress distribution,and high-order deformation pattern.The experimental and numerical results showed that compared to the concentrically loaded BRBs,the stiffness,yield force,cumulated plastic ductility(CPD)coefficient,equivalent viscous damping coefficient and energy dissipation decreased,and the yield displacement and compression strength adjustment factor increased for the eccentrically loaded BRBs.With the existence of the out-of-plane eccentricity,the initial yield position changes from the yield segment to the junction between the yield segment and transition segment under a tensile load,while the initial high-order buckling pattern changes from a first-order C-shape to a secondorder S-shape under a compressive load.

Numerical Investigation on the Harbour Resonance Inside A Harbour with Lateral Cavities

Wave-induced harbour resonance is numerically investigated inside a harbour with lateral cavities.The theoretical solutions for the amplification parameter are compared with the simulated results under varying dimensionless wave numbers in order to verify the simulation model in a rectangular harbour at a constant depth.The results indicate that the numerical model can correctly calculate the natural frequency and the natural wave height.A range of calculations are performed for harbour resonance with one pair of lateral cavities,two pairs of lateral cavities and three pairs of lateral cavities,respectively.The simulated results indicate that the amplitude of the amplification parameter decreases both at the primary natural oscillation and the secondary natural oscillation,as the number of lateral cavities increases.The dimensionless wave number reduces as the number of lateral cavities increases both at the primary natural oscillation and the secondary natural oscillation as well.

Overall Assessment of Heat Transfer for a Rarefied Flow in a Microchannel with Obstacles Using Lattice Boltzmann Method

The objective of this investigation is to assess the effect of obstacles on numerical heat transfer and fluid flow momentum in a rectangular microchannel(MC).Two distinct configurations were studied:one without obstacles and the other with alternating obstacles placed on the upper and lower walls.The research utilized the thermal lattice Boltzmann method(LBM),which solves the energy and momentum equations of fluids with the BGK approximation,implemented in a Python coding environment.Temperature jump and slip velocity conditions were utilized in the simulation for the MC and extended to all obstacle boundaries.The study aims to analyze the rarefaction effect,with Knudsen numbers(Kn)of 0.012,0.02,and 0.05.The outcomes indicate that rarefaction has a significant impact on the velocity and temperature distribution.The presence of nine obstacles led to slower fluid movement inside the microchannel MC,resulting in faster cooling at the outlet.In MCs with obstacles,the rarefaction effect plays a crucial role in decreasing the Nusselt number(Nu)and skin friction coefficient(Cf).Furthermore,the study demonstrated that the obstacles played a crucial role in boosting fluid flow and heat transfer in the MC.The findings suggest that the examined configurations could have potential applications as cooling technologies in micro-electro-mechanical systems and microdevice applications.

Creation and annihilation of artificial magnetic skyrmions with the electric field

Recent theory and experiments show that artificial magnetic skyrmions can be stabilized at room temperature without the need for the external magnetic field,casting strong potentials for the device applications.In this work,we study the electric field manipulation of artificial magnetic skyrmions imprinted by Co disks on CoPt multilayers utilizing the micromagnetic simulations.We find that the reversible annihilation and creation of skyrmions can be realized with the electric field via the strain mediated magnetoelastic coupling.In addition,we also demonstrate controllable manipulation of individual skyrmion,which opens a new platform for constructing magnetic field-free and low-energy dissipation skyrmion based media.

Toward Next-Generation Heterogeneous Catalysts:Empowering Surface Reactivity Prediction with Machine Learning

Heterogeneous catalysis remains at the core of various bulk chemical manufacturing and energy conversion processes,and its revolution necessitates the hunt for new materials with ideal catalytic activities and economic feasibility.Computational high-throughput screening presents a viable solution to this challenge,as machine learning(ML)has demonstrated its great potential in accelerating such processes by providing satisfactory estimations of surface reactivity with relatively low-cost information.This review focuses on recent progress in applying ML in adsorption energy prediction,which predominantly quantifies the catalytic potential of a solid catalyst.ML models that leverage inputs from different categories and exhibit various levels of complexity are classified and discussed.At the end of the review,an outlook on the current challenges and future opportunities of ML-assisted catalyst screening is supplied.We believe that this review summarizes major achievements in accelerating catalyst discovery through ML and can inspire researchers to further devise novel strategies to accelerate materials design and,ultimately,reshape the chemical industry and energy landscape.

Ground settlement and tunnel response due to twin-curved shield tunnelling in soft ground with small clear distance

Twin curved tunnels are often encountered in shield tunnelling,where significant complexities in densely exploited underground space are observed.In this study,the ground settlement and tunnel deformation due to twin-curved shield tunnelling in soft ground were investigated using numerical simulation and field monitoring.Different curvature radii of twin curved tunnels and subsequent effects of tunnel construction were considered to reveal the tunnelling effect on ground surface settlement and tunnel deformation.The results show that the settlement trough yields one offset towards inside of curved shield tunnelling.The location of settlement trough and maximum settlement were affected by curvature radius but except for the shape and width of settlement trough.Adjacent parallel twin-curved shield tunnelling could increase the offset of existing settlement trough and maximum settlement.Then,an empirical prediction of surface settlement trough due to twin-curved shield tunnelling with same tunnel diameters in soft clay was proposed,which was applicable to curvature radius less than 800 m.Finally,a minimum radius of 600 m of curvature tunnel was proposed in terms of allowable convergence deformation of tunnel.The result could provide guidance on safety evaluation for twin curved shield tunnelling construction.

A NOVEL STOCHASTIC HEPATITIS B VIRUS EPIDEMIC MODEL WITH SECOND-ORDER MULTIPLICATIVE α-STABLE NOISE AND REAL DATA

This work presents an advanced and detailed analysis of the mechanisms of hepatitis B virus(HBV)propagation in an environment characterized by variability and stochas-ticity.Based on some biological features of the virus and the assumptions,the corresponding deterministic model is formulated,which takes into consideration the effect of vaccination.This deterministic model is extended to a stochastic framework by considering a new form of disturbance which makes it possible to simulate strong and significant fluctuations.The long-term behaviors of the virus are predicted by using stochastic differential equations with second-order multiplicative α-stable jumps.By developing the assumptions and employing the novel theoretical tools,the threshold parameter responsible for ergodicity(persistence)and extinction is provided.The theoretical results of the current study are validated by numerical simulations and parameters estimation is also performed.Moreover,we obtain the following new interesting findings:(a)in each class,the average time depends on the value ofα;(b)the second-order noise has an inverse effect on the spread of the virus;(c)the shapes of population densities at stationary level quickly changes at certain values of α.The last three conclusions can provide a solid research base for further investigation in the field of biological and ecological modeling.

Recent advances in protein conformation sampling by combining machine learning with molecular simulation

The rapid advancement and broad application of machine learning(ML)have driven a groundbreaking revolution in computational biology.One of the most cutting-edge and important applications of ML is its integration with molecular simulations to improve the sampling efficiency of the vast conformational space of large biomolecules.This review focuses on recent studies that utilize ML-based techniques in the exploration of protein conformational landscape.We first highlight the recent development of ML-aided enhanced sampling methods,including heuristic algorithms and neural networks that are designed to refine the selection of reaction coordinates for the construction of bias potential,or facilitate the exploration of the unsampled region of the energy landscape.Further,we review the development of autoencoder based methods that combine molecular simulations and deep learning to expand the search for protein conformations.Lastly,we discuss the cutting-edge methodologies for the one-shot generation of protein conformations with precise Boltzmann weights.Collectively,this review demonstrates the promising potential of machine learning in revolutionizing our insight into the complex conformational ensembles of proteins.

Experimental Study on the Variation of Optical Imaging Characteristics with Zenith Angle due to Internal Solitary Waves in Sunglint

Internal solitary waves(ISWs)change the roughness of the sea surface,thus producing dark and bright bands in optical images.However,reasons for changes in imaging characteristics with the solar zenith angle remain unclear.In this paper,the optical imaging pattern of ISWs in sunglint under different zenith angles of the light source is investigated by collecting optical images of ISWs through physical simulation.The experiment involves setting 10 zenith angles of the light source,which are divided into area a the optical images of ISWs in the three areas show dark-bright mode,single bright band,and bright-dark mode,which are consistent with those observed by optical remote sensing.In addition,this study analyzed the percentage of the dark and bright areas of the bands and the change in the relative gray difference and found changes in both areas under different zenith angles of the light source.The MODIS and ASAR images display a similar brightness-darkness distance of the same ISWs.Therefore,the relationship between the brightness-darkness distance and the characteristic half-width of ISWs is determined in accordance with the eKdV theory and the imaging mechanism of ISWs of the SAR image.Overall,the relationship between them in the experiment is almost consistent with the theoretical result.

Damage Characteristics Analysis and Fractal Study of Shale With Prefabricated Fractures under Thermal-mechanical Coupling

To study the damage and failure of shale with different fracture inclination angles under uniaxial compression loading,in this work,RFPA2D-Thermal,a two-dimensional real failure process analysis software,was used for numerical simulation.Numerical simulation results show that quartz in shale mainly affects the tensile and compressive strength of shale by increasing rock brittleness.The coupling of temperature and pressure will cause lateral and volume destruction of shale,which enables the shale body to be more easily broken.Fracture inclination is the key factor affecting shale damage patterns.The failure mode of shale with low-and high-angle fractures is mainly shear failure,and the compressive strength does not vary with crack inclination.The damage mode of obliquely intersecting fractured shale is slip damage along the fracture face,the compressive strength decreases and then increases with the fracture inclination,and a minimum value exists.The acoustic emission simulation results of the damage process effectively reflect the accumulated internal damage and macroscopic crack appearance until fracture instability when the prefabricated fractured shale is subjected to uniaxial compressive loading.The crack inclinations of 0°and 120℃ corresponds to the most complex"N"shape damage mode.The crack inclinations of 30°and 60°,and the damage mode is an inverted"λ"shape.

Simulating the evolution of focused waves by a two-layer Boussinesq-type model

Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simulate deep-water focused waves of a two-layer Boussinesq-type(BT)model,which has been shown to have excellent linear and nonlinear performance.To further improve the numerical accuracy and stability,the internal wavegenerated method is introduced into the two-layer Boussinesq-type model.Firstly,the sensitivity of the numerical results to the grid resolution is analyzed to verify the convergence of the model;secondly,the focused wave propagating in two opposite directions is simulated to prove the symmetry of the numerical results and the feasibility of the internal wave-generated method;thirdly,the limiting focused wave condition is simulated to compare and analyze the wave surface and the horizontal velocity of the profile at the focusing position,which is in good agreement with the measured values.Meanwhile the simulation of focused waves in very deep waters agrees well with the measured values,which further demonstrates the capability of the two-layer BT model in simulating focused waves in deep waters.

Local Torrential Rainfall Event within a Mei-Yu Season Mesoscale Convective System:Importance of Back-Building Processes

An extreme rainfall event occurred over Hangzhou,China,during the afternoon hours on 24 June 2013.This event occurred under suitable synoptic conditions and the maximum 4-h cumulative rainfall amount was over 150 mm.This rainfall event had two major rainbands.One was caused by a quasi-stationary convective line,and the other by a backbuilding convective line related to the interaction of the outflow boundary from the first rainband and an existing low-level mesoscale convergence line associated with a mei-yu frontal system.The rainfall event lasted 4 h,while the back-building process occurred in 2 h when the extreme rainfall center formed.So far,few studies have examined the back-building processes in the mei-yu season that are caused by the interaction of a mesoscale convergence line and a convective cold pool.The two rainbands are successfully reproduced by the Weather Research and Forecasting(WRF)model with fourlevel,two-way interactive nesting.In the model,new cells repeatedly occur at the west side of older cells,and the backbuilding process occurs in an environment with large CAPE,a low LFC,and plenty of water vapor.Outflows from older cells enhance the low-level convergence that forces new cells.High precipitation efficiency of the back-building training cells leads to accumulated precipitation of over 150 mm.Sensitivity experiments without evaporation of rainwater show that the convective cold pool plays an important role in the organization of the back-building process in the current extreme precipitation case.

Numerical simulations for radon migration and exhalation behavior during measuring radon exhalation rate with closed-loop method

Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.

Growth mechanism and characteristics of electron drift instability in Hall thruster with different propellant types

The existence of a significant electron drift instability(EDI) in the Hall thruster is considered as one of the possible causes of the abnormal increase in axial electron mobility near the outlet of the channel. In recent years, extensive simulation research on the characteristics of EDI has been conducted, but the excitation mechanism and growth mechanism of EDI in linear stage and nonlinear stage remain unclear. In this work, a one-dimensional PIC model in the azimuthal direction of the thruster near-exit region is established to gain further insights into the mechanism of the EDI in detail, and the effects of different types of propellants on EDI characteristics are discussed. The changes in axial electron transport caused by EDI under different types of propellants and electromagnetic field strengths are also examined. The results indicate that EDI undergoes a short linear growth phase before transitioning to the nonlinear phase and finally reaching saturation through the ion Landau damping. The EDI drives a significant ion heating in the azimuthal direction through electron–ion friction before entering the quasi-steady state, which increases the axial mobility of the electrons. Using lighter atomic weight propellant can effectively suppress the oscillation amplitude of EDI, but it will increase the linear growth rate, frequency, and phase velocity of EDI. Compared with the classical mobility, the axial electron mobility under the EDI increases by three orders of magnitude, which is consistent with experimental phenomena. The change of propellant type is insufficient to significantly change the axial electron mobility. It is also found that the collisions between electrons and neutral gasescan significantly affect the axial electron mobility under the influence of EDI, and lead the strength of the electric field to increase and the strength of the magnetic field to decrease, thereby both effectively suppressing the axial transport of electrons.