A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Inelastic Scattering of Dark Matter with Heavy Cosmic Rays

We investigate the impact of inelastic collisions between dark matter(DM)and heavy cosmic ray(CR)nuclei on CR propagation.We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections,allowing us to assess how these collisions affect the spectra of CR boron and carbon.We derive new CR spectra from DM-CR collisions by incorporating their cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species.In a specific example with a coupling strength of b_(χ)=0.1 and a DM mass of m_(χ)=0.1 GeV,considering a simplified scenario where DM interacts exclusively with oxygen,a notable modification in the boron-to-carbon spectrum due to the DM-CR interaction is observed.Particularly,the peak within the spectrum,spanning from 0.1 to 10 GeV,experiences an enhancement of approximately 1.5 times.However,in a more realistic scenario where DM particles interact with all CRs,this peak can be amplified to twice its original value.Utilizing the latest data from AMS-02 and DAMPE on the boron-to-carbon ratio,we estimate a 95%upper limit for the effective inelastic cross-section of DM-proton as a function of DM mass.Our findings reveal that at m_(χ)?2 MeV,the effective inelastic cross-section between DM and protons must be less than O(10^(-32))cm^(2).

Time-resolved ultra-small-angle X-ray scattering beamline(BL10U1)at SSRF

The construction of a new beamline,BL10U1,was completed at the Shanghai synchrotron radiation facility in 2020.This multipurpose beamline was designed to provide X-ray scattering techniques such as ultra-small-angle X-ray scattering(USAXS),small-angle X-ray scattering(SAXS),wide-angle X-ray scattering,and microfocus SAXS(μSAXS)for a broad user community.To realize fast time-resolved USAXS experiments,the beamline adopted an in-vacuum undulator with a total length of 1.6 m as the photon source.An in-house cryogenic-cooled double multilayer monochromator was installed to deliver a photon flux of approximately 10^(13) photons/s at a photon energy of 10 keV.The three-year successful operation of this beamline demonstrated that the monochromator operated smoothly,as expected.BL10U1 has three end stations in succession:USAXS end station,μSAXS end station,and end station for industrial applications.The minimum scattering vector q~0.0042 nm^(-1) at 10 keV can be achieved at the USAXS end station equipped with a 28 m-long and 1.8 m-diameter vacuum flight tube.At theμSAXS end station,a beam spot of less than 10×8μm was achieved for micro-SAXS experiments.In contrast,in situ experimental instruments up to 5 m high and 8 m wide can be mounted at the industrial application end station,which offers industrial scientists the opportunity to use their large industrial equipment.BL10U1 opens up a new capability to investigate phenomena such as non-equilibrium and dynamic processes of materials with a wide length scale from angstroms to micrometers with millisecond time resolution.In this paper,we also report beamline design considerations and commissioning results.

Study of deep transportation and plugging performance of deformable gel particles in porous media

Deformable gel particles(DGPs) possess the capability of deep profile control and flooding. However, the deep migration behavior and plugging mechanism along their path remain unclear. Breakage, an inevitable phenomenon during particle migration, significantly impacts the deep plugging effect. Due to the complexity of the process, few studies have been conducted on this subject. In this paper, we conducted DGP flow experiments using a physical model of a multi-point sandpack under various injection rates and particle sizes. Particle size and concentration tests were performed at each measurement point to investigate the transportation behavior of particles in the deep part of the reservoir. The residual resistance coefficient and concentration changes along the porous media were combined to analyze the plugging performance of DGPs. Furthermore, the particle breakage along their path was revealed by analyzing the changes in particle size along the way. A mathematical model of breakage and concentration changes along the path was established. The results showed that the passage after breakage is a significant migration behavior of particles in porous media. The particles were reduced to less than half of their initial size at the front of the porous media. Breakage is an essential reason for the continuous decreases in particle concentration, size, and residual resistance coefficient. However, the particles can remain in porous media after breakage and play a significant role in deep plugging. Higher injection rates or larger particle sizes resulted in faster breakage along the injection direction, higher degrees of breakage, and faster decreases in residual resistance coefficient along the path. These conditions also led to a weaker deep plugging ability. Smaller particles were more evenly retained along the path, but more particles flowed out of the porous media, resulting in a poor deep plugging effect. The particle size is a function of particle size before injection, transport distance, and different injection parameters(injection rate or the diameter ratio of DGP to throat). Likewise, the particle concentration is a function of initial concentration, transport distance, and different injection parameters. These models can be utilized to optimize particle injection parameters, thereby achieving the goal of fine-tuning oil displacement.

Suppression of stimulated Brillouin and Raman scatterings using an alternating frequency laser and transverse magnetic fields

A novel scheme to suppress both stimulated Brillouin scattering(SBS) and stimulated Raman scattering(SRS) by combining an alternating frequency(AF) laser and a transverse magnetic field is proposed. The AF laser allows the laser frequency to change discretely and alternately over time. The suppression of SBS is significant as long as the AF difference is greater than the linear growth rate of SBS or the alternating time of the laser frequency is shorter than the linear growth time of SBS. However, the AF laser proves ineffective in suppressing SRS, which usually has a much higher linear growth rate than SBS. To remedy that, a transverse magnetic field is included to suppress the SRS instability. The electrons trapped in the electron plasma waves(EPWs) of SRS can be accelerated by the surfatron mechanism in a transverse magnetic field and eventually detrapped. While continuously extracting energy from EPWs, the EPWs are dissipated and the kinetic inflation of SRS is suppressed. The one-dimensional particle-in-cell simulation results show that both SBS and SRS can be effectively suppressed by combining the AF laser with a transverse magnetic field with tens of Tesla. The total reflectivity can be dramatically reduced by more than one order of magnitude. These results provide a potential reference for controlling SBS and SRS under the related parameters of inertial confinement fusion.

Computer vision-aided DEM study on the compaction characteristics of graded subgrade filler considering realistic coarse particle shapes

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.

Mass transfer enhancement and hydrodynamic performance with wire mesh coupling solid particles in bubble column reactor

It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on bubble motion and gas-liquid mass transfer process in a bubble column.Particle image velocimetry(PIV)technology was employed to analyze the flow field and bubble motion behavior,and dynamic oxygen absorption technology was used to measure the gas-liquid volumetric mass transfer coefficient(kLa).The effect of embedding wire mesh,adding solid particles,and wire mesh coupling solid particles on the flow characteristic and kLa were analyzed and compared.The results show that the gas-liquid interface area increases by 33%-72%when using the wire mesh coupling solid particles strategy compared to the gas-liquid two-phase flow,which is superior to the other two strengthening methods.Compared with the system without reinforcement,kLa in the bubble column increased by 0.5-1.8 times with wire mesh coupling solid particles method,which is higher than the sum of kLa increases with inserting wire mesh and adding particles,and the coupling reinforcement mechanism for affecting gas-liquid mass transfer process was discussed to provide a new idea for enhancing gas-liquid mass transfer.

Study of the spatial growth of stimulated Brillouin scattering in a gas-filled hohlraum via detecting the driven ion acoustic wave

In an experiment performed on the Shenguang-III prototype laser facility, collective Thomson scattering (TS) is used to study the spatialgrowth of stimulated Brillouin scattering (SBS) in a gas-filled hohlraum by detecting the SBS-driven ion acoustic wave. High-quality timeresolved SBS and TS spectra are obtained simultaneously in the experiment, and these are analyzed by a steady-state code based on theray-tracing model. The analysis indicates that ion–ion collisions may play an important role in suppressing SBS growth in the Au plasma;as aresult, the SBS excited in the filled gas region is dominant. In the early phase of the laser pulse, SBS originates primarily from the high-densityplasma at the edges of the interaction beam channel, which is piled up by the heating of the interaction beam. Throughout the duration of thelaser pulse, the presence of the TS probe beam might mitigate SBS by perturbing the density distribution around the region overlapping withthe interaction beam.

Fully polarized Compton scattering in plane waves and its polarization transfer

Fully polarized Compton scattering from a beam of spin-polarized electrons is investigated in plane-wave backgrounds in a broad intensity region from the perturbative to the nonperturbative regimes.In the perturbative regime,polarized linear Compton scattering is considered for investigating polarization transfer from a single laser photon to a scattered photon,and in the high-intensity region,the polarized locally monochromatic approximation and locally constant field approximation are established and are employed to study polarization transfer from an incoming electron to a scattered photon.The numerical results suggest an appreciable improvement of about 10%in the scattering probability in the intermediate-intensity region if the electron’s longitudinal spin is parallel to the laser rotation.The longitudinal spin of the incoming electron can be transferred to the scattered photon with an efficiency that increases with laser intensity and collisional energy.For collision between an optical laser with frequency1 eV and a 10 GeV electron,this polarization transfer efficiency can increase from about 20%in the perturbative regime to about 50%in the nonperturbative regime for scattered photons with relatively high energy.

Backward scattering of laser plasma interactions from hundreds-of-joules broadband laser on thick target

The use of broadband laser technology is a novel approach for inhibiting processes related to laser plasma interactions(LPIs).In this study,several preliminary experiments into broadband-laser-driven LPIs are carried out using a newly established hundreds-of-joules broadband second-harmonic-generation laser facility.Through direct comparison with LPI results for a traditional narrowband laser,the actual LPI-suppression effect of the broadband laser is shown.The broadband laser had a clear suppressive effect on both back-stimulated Raman scattering and back-stimulated Brillouin scattering at laser intensities below 1×10^(15) W cm^(−2).An abnormal hot-electron phenomenon is also investigated,using targets of different thicknesses.

Electron kinetic effects in back-stimulated Raman scattering bursts driven by broadband laser pulses

We examine electron kinetic effects in broadband-laser-driven back-stimulated Raman scattering(BSRS)bursts using particle-in-cell simulations.These bursts occur during the nonlinear stage,causing reflectivity spikes and generating large numbers of hot electrons.Long-duration simulations are performed to observe burst events,and a simplified model is developed to eliminate the interference of the broadband laser’s random intensity fluctuations.Using the simplified model,we isolate and characterize the spectrum of electron plasma waves.The spectrum changes from a sideband structure to a turbulence-like structure during the burst.A significant asymmetry in the spectrum is observed.This asymmetry is amplified and transferred to electron phase space by high-intensity broadband laser pulses,leading to violent vortex-merging and generation of hot electrons.The proportion of hot electrons increases from 6.76%to 14.7%during a single violent burst event.We demonstrate that kinetic effects profoundly influence the BSRS evolution driven by broadband lasers.

Correction:Brillouin scattering spectrum for liquid detection and applications in oceanography

After publication of this article1,it was brought to our at-tention that the mathematical expressions‘‰’were mis-takenly replaced by‘%’for salinities.Details are listed below.1.In the last sentence in abstract,“approximately 0.1℃and 0.5%”should be“approximately 0.1℃and 0.5‰”.

State Estimation of Drive-by-Wire Chassis Vehicle Based on Dual Unscented Particle Filter Algorithm

Accurate vehicle dynamic information plays an important role in vehicle driving safety.However,due to the characteristics of high mobility and multiple controllable degrees of freedom of drive-by-wire chassis vehicles,the current mature application of traditional vehicle state estimation algorithms can not meet the requirements of drive-by-wire chassis vehicle state estimation.This paper proposes a state estimation method for drive-by-wire chassis vehicle based on the dual unscented particle filter algorithm,which make full use of the known advantages of the four-wheel drive torque and steer angle parameters of the drive-by-wire chassis vehicle.In the dual unscented particle filter algorithm,two unscented particle filter transfer information to each other,observe the vehicle state information and the tire force parameter information of the four wheels respectively,which reduce the influence of parameter uncertainty and model parameter changes on the estimation accuracy during driving.The performance with the dual unscented particle filter algorithm,which is analyzed in terms of the time-average square error,is superior of the unscented Kalman filter algorithm.The effectiveness of the algorithm is further verified by driving simulator test.In this paper,a vehicle state estimator based on dual unscented particle filter algorithm was proposed for the first time to improve the estimation accuracy of vehicle parameters and states.

Optimizing laser coupling,matter heating,and particle acceleration from solids using multiplexed ultraintense lasers

Realizing the full potential of ultrahigh-intensity lasers for particle and radiation generation will require multi-beam arrangements due to technology limitations.Here,we investigate how to optimize their coupling with solid targets.Experimentally,we show that overlapping two intense lasers in a mirror-like configuration onto a solid with a large preplasma can greatly improve the generation of hot electrons at the target front and ion acceleration at the target backside.The underlying mechanisms are analyzed through multidimensional particle-in-cell simulations,revealing that the self-induced magnetic fields driven by the two laser beams at the target front are susceptible to reconnection,which is one possible mechanism to boost electron energization.In addition,the resistive magnetic field generated during the transport of the hot electrons in the target bulk tends to improve their collimation.Our simulations also indicate that such effects can be further enhanced by overlapping more than two laser beams.

Significance of including lid thickness and particle shape factor in numerical modeling for prediction of particle trap efficiency of invert trap

Sediment accumulation on the bed of open sewers and drains reduces hydraulic efficiency and can cause localized flooding.Slotted invert traps installed underneath the bed of open sewers and drains can eliminate sediment build-up by catching sediment load.Previous three-dimensional(3D)computational studies have examined the particle trapping performance of invert traps of different shapes and depths under varied sediment and flow conditions,considering particles as spheres.For two-dimensional and 3D numerical modeling,researchers assumed the lid geometry to be a thin line and a plane,respectively.In this 3D numerical study,the particle trapping efficiency of a slotted irregular hexagonal invert trap fitted at the flume bottom was examined by incorporating the particle shape factor of non-spherical sewage solid particles and the thicknesses of upstream and downstream lids over the trap in the discrete phase model of the ANSYS Fluent 2020 R1 software.The volume of fluid(VOF)and the realizable k-turbulence models were used to predict the velocity field.The two-dimensional particle image velocimetry(PIV)was used to measure the velocity field inside the invert trap.The results showed that the thicknesses of upstream and downstream lids affected the velocity field and turbulent kinetic energy at all flow depths.The joint impact of the particle shape factor and lid thickness on the trap efficiency was significant.When both the lid thickness and particle shape factor were considered in the numerical modeling,trap efficiencies were underestimated,with relative errors of-8.66%to-0.65%in comparison to the experimental values of Mohsin and Kaushal(2017).They were also lower than the values predicted by Mohsin and Kaushal(2017),which showed an overall overestimation with errors of-2.3%to 17.4%.

Numerical study of alpha particle loss with toroidal field ripple based on CFETR steady-state scenario

Effects of plasma equilibrium parameters on the alpha particle loss with the toroidal field ripple based on the CFETR steady-state scenario have been numerically investigated by the orbit-following code GYCAVA. It is found that alpha particle losses decrease and loss regions become narrower with the plasma current increasing or with the magnetic field decreasing. It is because the ripple stochastic transport and the ripple well loss of alpha particle are reduced with the safety factor decreasing. Decrease of the plasma density and temperature can reduce alpha particle losses due to enhancement of the slowing-down effect. The direction of the toroidal magnetic field can significantly affect heat loads induced by lost alpha particle. The vertical asymmetry of heat loads induced by the clockwise and counter-clockwise toroidal magnetic fields are due to the fact that the ripple distribution is asymmetric about the mid-plane, which can be explained by the typical orbits of alpha particle. The maximal heat load of alpha particle for the clockwise toroidal magnetic field is much smaller than that for the counter-clockwise one.

Single-photon scattering and quantum entanglement of two giant atoms with azimuthal angle differences in a waveguide system

We theoretically investigate coherent scattering of single photons and quantum entanglement of two giant atoms with azimuthal angle differences in a waveguide system.Using the real-space Hamiltonian,analytical expressions are derived for the transport spectra scattered by these two giant atoms with four azimuthal angles.Fano-like resonance can be exhibited in the scattering spectra by adjusting the azimuthal angle difference.High concurrence of the entangled state for two atoms can be implemented in a wide angle-difference range,and the entanglement of the atomic states can be switched on/off by modulating the additional azimuthal angle differences from the giant atoms.This suggests a novel handle to effectively control the single-photon scattering and quantum entanglement.

Numerical simulation for the initial state of avalanche in polydisperse particle systems

Numerical simulation is employed to investigate the initial state of avalanche in polydisperse particle systems.Nucleation and propagation processes are illustrated for pentadisperse and triadisperse particle systems,respectively.In these processes,particles involved in the avalanche grow slowly in the early stage and explosively in the later stage,which is clearly different from the continuous and steady growth trend in the monodisperse system.By examining the avalanche propagation,the number growth of particles involved in the avalanche and the slope of the number growth,the initial state can be divided into three stages:T1(nucleation stage),T2(propagation stage),T3(overall avalanche stage).We focus on the characteristics of the avalanche in the T2 stage,and find that propagation distances increase almost linearly in both axial and radial directions in polydisperse systems.We also consider the distribution characteristics of the average coordination number and average velocity for the moving particles.The results support that the polydisperse particle systems are more stable in the T2 stage.

Particle aggregation and breakage kinetics in cemented paste backfill

The macroscopic flow behavior and rheological properties of cemented paste backfill(CPB)are highly impacted by the inherent structure of the paste matrix.In this study,the effects of shear-induced forces and proportioning parameters on the microstructure of fresh CPB were studied.The size evolution and distribution of floc/agglomerate/particles of paste were monitored by focused beam reflection measuring(FBRM)technique,and the influencing factors of aggregation and breakage kinetics of CPB were discussed.The results indicate that influenced by both internal and external factors,the paste kinetics evolution covers the dynamic phase and the stable phase.Increasing the mass content or the cement-tailings ratio can accelerate aggregation kinetics,which is advantageous for the rise of average floc size.Besides,the admixture and high shear can improve breaking kinetics,which is beneficial to reduce the average floc size.The chord length resembles a normal distribution somewhat,with a peak value of approximate 20μm.The particle disaggregation con-stant(k_(2))is positively correlated with the agitation rate,and k_(2) is five orders of magnitude greater than the particle aggregation constant(k1).The kinetics model depicts the evolution law of particles over time quantitatively and provides a theoretical foundation for the micromechanics of complicated rheological behavior of paste.

High-Volume Mineral Admixtures Cement: The Effects of Particle Size Distribution

The effects of high-volume slag-fly ash cement with different particle sizes on hydration degree,microstructure and mechanical properties were systematically studied,by means of laser particle size(DLS),X-ray diffraction (XRD),comprehensive thermal analysis (TG-DTA),scanning electron microscopy(SEM) and mechanical properties tests.The results show that suitable particle size distribution of cementitious material has significantly promoting effects on hydration reaction rate and mechanical properties.Compared with slag without further grinding,the slag after ball milling for 4 h has an obvious improvement in reactivity,which also provides a faster hydration rate and higher compressive strength for the cementitious material.When the slag milled for 1 and 4 h is mixed at a mass ratio of 2:1 (i e,slag with D_(50) of 7.4μm and average size of 9.9μm,and slag with D_(50) value of 2.6μm and average size of 5.3μm),and a certain amount of fly ash is added in,the most obvious improvement of compressive strength of cement is achieved.