We construct optimal k-step, 5- to 10-stage, explicit, strong-stability-preserving Hermite-Birkhoff (SSP HB) methods of order 12 with nonnegative coefficients by combining linear k-step methods of order 9 with 5- to 1...We construct optimal k-step, 5- to 10-stage, explicit, strong-stability-preserving Hermite-Birkhoff (SSP HB) methods of order 12 with nonnegative coefficients by combining linear k-step methods of order 9 with 5- to 10-stage Runge-Kutta (RK) methods of order 4. Since these methods maintain the monotonicity property, they are well suited for solving hyperbolic PDEs by the method of lines after a spatial discretization. It is seen that the 8-step 7-stage HB methods have largest effective SSP coefficient among the HB methods of order 12 on hand. On Burgers’ equations, some of the new HB methods have larger maximum effective CFL numbers than Huang’s 7-step hybrid method of order 7, thus allowing larger step size.展开更多
In this paper,we extend the work of Brenner and Sung[Math.Comp.59,321–338(1992)]and present a regularity estimate for the elastic equations in concave domains.Based on the regularity estimate we prove that the consta...In this paper,we extend the work of Brenner and Sung[Math.Comp.59,321–338(1992)]and present a regularity estimate for the elastic equations in concave domains.Based on the regularity estimate we prove that the constants in the error estimates of the nonconforming Crouzeix-Raviart element approximations for the elastic equations/eigenvalue problem are independent of Laméconstant,which means the nonconforming Crouzeix-Raviart element approximations are locking-free.We also establish two kinds of two-grid discretization schemes for the elastic eigenvalue problem,and analyze that when the mesh sizes of coarse grid and fine grid satisfy some relationship,the resulting solutions can achieve the optimal accuracy.Numerical examples are provided to show the efficiency of two-grid schemes for the elastic eigenvalue problem.展开更多
A 32-channel charge-sensitive amplifier(CSA)is designed for fast timing in the delay-line readout of a parallel plate avalanche counter(PPAC)array.It is realized on a PCB with operational amplifiers and other discrete...A 32-channel charge-sensitive amplifier(CSA)is designed for fast timing in the delay-line readout of a parallel plate avalanche counter(PPAC)array.It is realized on a PCB with operational amplifiers and other discrete components.Each channel consists of an integrator,a pole-zero cancellation net,and a linear amplification stage,which can be adapted to accommodate either positive or negative input signals.The RMS equivalent input noise charges are 3.3 fC,the conversion gains are approximately±2 mV∕fC,and the intrinsic time resolution reaches 32 ps.In the prototype PPAC application,the CSA performs as well as the commercial FTA820A amplifier,providing a position resolution as good as 0.17 mm,and exhibiting reliable stability during several hours of continuous data acquisition.展开更多
In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical m...In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical model,and its validity was verified using a simple impact test.A crushable discrete element method(DEM)framework is built based on the previously established theoretical model.The tensile strength,which considers the fractal theory,size effect,and Weibull variation,was assigned to each generated particle.The assigned strength is then used for crush detection by comparing it with its maximum tensile stress.Mass conservation is ensured by inserting a series of sub-particles whose total mass was equal to the quality loss.Based on the crushable DEM framework,a numerical simulation of the crushing behavior of a pebble bed with hollow cylindrical geometry under a uniaxial compression test was performed.The results of this investigation showed that the particle withstands the external load by contact and sliding at the beginning of the compression process,and the results confirmed that crushing can be considered an important method of resisting the increasing external load.A relatively regular particle arrangement aids in resisting the load and reduces the occurrence of particle crushing.However,a limit exists to the promotion of resistance.When the strain increases beyond this limit,the distribution of the crushing position tends to be isotropic over the entire pebble bed.The theoretical model and crushable DEM framework provide a new method for exploring the pebble bed in a fusion reactor,considering particle crushing.展开更多
In contrast to cyclic polymers with ring-like backbones,side-chain cyclization is another intriguing structural feature that has not been extensively studied.In this study,a library of orthogonally protected monomers ...In contrast to cyclic polymers with ring-like backbones,side-chain cyclization is another intriguing structural feature that has not been extensively studied.In this study,a library of orthogonally protected monomers featuring monocyclic,dicyclic,or tricyclic pendant motifs was designed and prepared based on malic acid derivatives.Polyesters with precise chemical structures and uniform chain lengths were prepared modularly through iterative growth.Meticulous control over the chemical details allows for a close investigation of the topological effects on the polymer properties.Compared to their linear side chain counterparts,the presence of cyclic pendant groups has a significant impact on chain conformation,leading to a reduction in hydrodynamic volume and an enhancement in the glass transition temperature.These results underscore the potential of tailoring polymer properties through rational engineering of side chain topology.展开更多
The pseudo-two-dimensional(P2D)model plays an important role in exploring physicochemical mechanisms,predicting the state of health,and improving the fast charge capability for Li-ion batteries(LIBs).However,the fast ...The pseudo-two-dimensional(P2D)model plays an important role in exploring physicochemical mechanisms,predicting the state of health,and improving the fast charge capability for Li-ion batteries(LIBs).However,the fast charge leads to the lithium concentration gradient in the solid and electrolyte phases and the non-uniform electrochemical reaction at the solid/electrolyte interface.In order to decouple charge transfer reactions in LIBs under dynamic conditions,understanding the spatio-temporal resolution of the P2D model is urgently required.Till now,the study of this aspect is still insufficient.This work studies the spatio-temporal resolution for dynamic/static electrochemical impedance spectroscopy(DEIS/SEIS)on multiple scales.In detail,DEIS and SEIS with spatio-temporal resolutions are used to decouple charge transfer reactions in LIBs based on the numerical solution of the P2D model in the frequency domain.The calculated results indicate that decoupling solid diffusion requires a high spatial resolution along the r-direction in particles,decoupling electrolyte diffusion and interfacial transfer reaction requires a high spatial resolution along the x-direction,and decoupling charge transfer reactions in LIBs at an extremely low state of charge(SOC)requires an extremely high temporal resolution along the t-direction.Finally,the optimal range of spatio-temporal resolutions for DEIS/SEIS is derived,and the method to decouple charge transfer reactions with spatio-temporal resolutions is developed.展开更多
Seabed mining operations have been found to induce significant movement and deformation in overlying rock strata,posing serious threats to mining safety.The presence of geological faults further complicates these defo...Seabed mining operations have been found to induce significant movement and deformation in overlying rock strata,posing serious threats to mining safety.The presence of geological faults further complicates these deformation patterns.This study utilized geophysical surveys and the continuum-based discrete element method(CDEM)to investigate how fault activity influences rock deformation and failure.The results demonstrate that:1)Acting in mechanically weak zones,faults exerted a pronounced barrier effect on deformation propagation and stress redistribution within the surrounding rock,leading to markedly divergent displacement patterns on either side of the fault plane.Comparative analyses between single-fault and double-fault models revealed an 18%−22%expansion of the damage zone under the latter,together with significantly intensified deformation and failure;2)The double-fault model exhibited a larger maximum cumulative vertical displacement and a spatial shift in the location of peak deformation,thereby posing a heightened threat to mine safety;3)Acting in an orebody substitute,backfill effectively constrained surrounding rock deformation,enhanced its load-bearing capacity,and delayed the overburden subsidence.Nevertheless,backfill only reduced the amplitude of deformation;it could not entirely prevent settlement.These findings provide essential theoretical insights and foundational knowledge for safer submarine mining practices.展开更多
Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture...Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.展开更多
The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearitie...The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearities.To discretize the continuous problem in space,we utilize a spectral Galerkin method.Subsequently,we introduce a nonlinear-tamed exponential integrator scheme,resulting in a fully discrete scheme.Within the framework of semigroup theory,this study provides precise estimations of the Sobolev regularity,L^(∞) regularity in space,and Hölder continuity in time for the mild solution,as well as for its semi-discrete and full-discrete approximations.Building upon these results,we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions.A numerical example is presented to validate the theoretical findings.展开更多
It is well known that the inhomogeneous Calderón-Zygmund convolution operators are bounded on the local Hardy spaces.In this paper,we prove that these operators are bounded on the local product Hardy spaces and t...It is well known that the inhomogeneous Calderón-Zygmund convolution operators are bounded on the local Hardy spaces.In this paper,we prove that these operators are bounded on the local product Hardy spaces and the Lipschitz spaces.The key ideas used here are the discrete local Calderón identity and a density argument for the inhomogeneous product Lipschitz spaces in the weak sense.展开更多
Fractures are typically characterized by roughness that significantlyaffects the mechanical and hydraulic characteristics of reservoirs.However,hydraulic fracturing mechanisms under the influenceof fracture morphology...Fractures are typically characterized by roughness that significantlyaffects the mechanical and hydraulic characteristics of reservoirs.However,hydraulic fracturing mechanisms under the influenceof fracture morphology remain largely unexplored.Leveraging the advantages of the finite-discrete element method(FDEM)for explicitly simulating fracture propagation and the strengths of the unifiedpipe model(UPM)for efficientlymodeling dual-permeability seepage,we propose a new hydromechanical(HM)coupling approach for modeling hydraulic fracturing.Validated against benchmark examples,the proposed FDEM-UPM model is further augmented by incorporating a Fourier-based methodology for reconstructing non-planar fractures,enabling quantitative analysis of hydraulic fracturing behavior within rough discrete fracture networks(DFNs).The FDEM-UPM model demonstrates computational advantages in accurately capturing transient hydraulic seepage phenomena,while the asynchronous time-stepping schemes between hydraulic and mechanical analyses substantially enhanced computational efficiencywithout compromising computational accuracy.Our results show that fracture morphology can affect both macroscopic fracture networks and microscopic interaction types between hydraulic fractures(HFs)and natural fractures(NFs).In an isotropic stress field,the initiation azimuth,propagation direction and microcracking mechanism are significantly influencedby fracture roughness.In an anisotropic stress field,HFs invariably propagate parallel to the direction of the maximum principal stress,reducing the overall complexity of the stimulated fracture networks.Additionally,stress concentration and perturbation attributed to fracture morphology tend to be compromised as the leak-off increases,while the breakdown and propagation pressures remain insensitive to fracture morphology.These findingsprovide new insights into the hydraulic fracturing mechanisms of fractured reservoirs containing complex rough DFNs.展开更多
Zero-dimensional(0D)organic-inorganic metal halide perovskite is one of the hot research topics in the field of optoelectronic materials.Their structure generally consists of discrete metal halide octahedra entirely i...Zero-dimensional(0D)organic-inorganic metal halide perovskite is one of the hot research topics in the field of optoelectronic materials.Their structure generally consists of discrete metal halide octahedra entirely isolated by surrounding organic cations,forming independent luminescent centers[1,2].Such a configuration results in high exciton binding energy and exceptional luminescence efficiency,due to strong quantum confinement[3,4].展开更多
Some two-scale finite element discretizations are introduced for a class of linear partial differential equations. Both boundary value and eigenvalue problems are studied. Based on the two-scale error resolution techn...Some two-scale finite element discretizations are introduced for a class of linear partial differential equations. Both boundary value and eigenvalue problems are studied. Based on the two-scale error resolution techniques, several two-scale finite element algorithms are proposed and analyzed. It is shown that this type of two-scale algorithms not only significantly reduces the number of degrees of freedom but also produces very accurate approximations.展开更多
In this paper, we use Hermite weighted essentially non-oscillatory （HWENO） schemes with a Lax-Wendroff time discretization procedure, termed HWENO-LW schemes, to solve Hamilton-Jacobi equations. The idea of the reco...In this paper, we use Hermite weighted essentially non-oscillatory （HWENO） schemes with a Lax-Wendroff time discretization procedure, termed HWENO-LW schemes, to solve Hamilton-Jacobi equations. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first derivative values are evolved in time and are used in the reconstruction. One major advantage of HWENO schemes is its compactness in the reconstruction. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with HWENO with Runge-Kutta time discretizations schemes （HWENO-RK） of Qiu and Shu [19] for Hamilton-Jacobi equations, the major advantages of HWENO-LW schemes are their saving of computational cost and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method.展开更多
In this paper,we discuss the local discontinuous Galerkin methods coupled with two specific explicitimplicit-null time discretizations for solving one-dimensional nonlinear diffusion problems Ut=(a(U)Ux)x.The basic id...In this paper,we discuss the local discontinuous Galerkin methods coupled with two specific explicitimplicit-null time discretizations for solving one-dimensional nonlinear diffusion problems Ut=(a(U)Ux)x.The basic idea is to add and subtract two equal terms a0 Uxx the right-hand side of the partial differential equation,then to treat the term a0 Uxx implicitly and the other terms(a(U)Ux)x-a0 Uxx explicitly.We give stability analysis for the method on a simplified model by the aid of energy analysis,which gives a guidance for the choice of a0,i.e.,a0≥max{a(u)}/2 to ensure the unconditional stability of the first order and second order schemes.The optimal error estimate is also derived for the simplified model,and numerical experiments are given to demonstrate the stability,accuracy and performance of the schemes for nonlinear diffusion equations.展开更多
This paper studies the two-stage fourth-order accurate time discretization[J.Q.Li and Z.F.Du,SIAM J.Sci.Comput.,38(2016)]and its application to the special relativistic hydrodynamical equations.Our analysis reveals th...This paper studies the two-stage fourth-order accurate time discretization[J.Q.Li and Z.F.Du,SIAM J.Sci.Comput.,38(2016)]and its application to the special relativistic hydrodynamical equations.Our analysis reveals that the new two-stage fourth-order accurate time discretizations can be proposed.With the aid of the direct Eulerian GRP(generalized Riemann problem)methods and the analytical resolution of the local“quasi 1D”GRP,the two-stage fourth-order accurate time discretizations are successfully implemented for the 1D and 2D special relativistic hydrodynamical equations.Several numerical experiments demonstrate the performance and accuracy as well as robustness of our schemes.展开更多
Upwind algorithms are becoming progressively popular for river flood routing due to their capability of resolving trans-critical flow regimes. For consistency, these algorithms suggest natural upwind discretization of...Upwind algorithms are becoming progressively popular for river flood routing due to their capability of resolving trans-critical flow regimes. For consistency, these algorithms suggest natural upwind discretization of the source term, which may be essential for natural channels with irregular geometry. Yet applications of these upwind algorithms to natural river flows are rare, and in such applications the traditional and simpler pointwise, rather than upwind discretization of the source term is used. Within the framework of a first-order upwind algorithm, this paper presents a comparison of upwind and pointwise discretizations of the source term. Numerical simulations were carried out for a selected irregular channel comprising a pool-riffle sequence Jn the River Lune, England with observed data. It is Shown that the impact of pointwise discretization, compared to the upwind, is appreciable mainly in flow zones with the Froude number closer to or larger than unity. The discrepancy due to pointwise and upwind discretizations of the source term is negligible in flow depth and hence in water surface elevation, but well manifested in mean velocity and derived flow quantities. Also the occurrence of flow reversal and equalisation over the pool-riffle sequence in response to increasing discharges is demonstrated.展开更多
Integrable discretizations of are proposed. N-soliton solutions for analogues of the complex and real Dym the complex and real Dym equations both semi-discrete and fully discrete equations are also presented.
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hyd...Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.展开更多
文摘We construct optimal k-step, 5- to 10-stage, explicit, strong-stability-preserving Hermite-Birkhoff (SSP HB) methods of order 12 with nonnegative coefficients by combining linear k-step methods of order 9 with 5- to 10-stage Runge-Kutta (RK) methods of order 4. Since these methods maintain the monotonicity property, they are well suited for solving hyperbolic PDEs by the method of lines after a spatial discretization. It is seen that the 8-step 7-stage HB methods have largest effective SSP coefficient among the HB methods of order 12 on hand. On Burgers’ equations, some of the new HB methods have larger maximum effective CFL numbers than Huang’s 7-step hybrid method of order 7, thus allowing larger step size.
基金supported by the National Natural Science Foundation of China (Grant No.11761022)。
文摘In this paper,we extend the work of Brenner and Sung[Math.Comp.59,321–338(1992)]and present a regularity estimate for the elastic equations in concave domains.Based on the regularity estimate we prove that the constants in the error estimates of the nonconforming Crouzeix-Raviart element approximations for the elastic equations/eigenvalue problem are independent of Laméconstant,which means the nonconforming Crouzeix-Raviart element approximations are locking-free.We also establish two kinds of two-grid discretization schemes for the elastic eigenvalue problem,and analyze that when the mesh sizes of coarse grid and fine grid satisfy some relationship,the resulting solutions can achieve the optimal accuracy.Numerical examples are provided to show the efficiency of two-grid schemes for the elastic eigenvalue problem.
基金supported by the National Natural Science Foundation of China(Nos.U2167202,12225504,12005276)the Natural Science Foundation of Shandong Province(No.ZR2024QA172)the Fundamental Research Funds of Shandong University.
文摘A 32-channel charge-sensitive amplifier(CSA)is designed for fast timing in the delay-line readout of a parallel plate avalanche counter(PPAC)array.It is realized on a PCB with operational amplifiers and other discrete components.Each channel consists of an integrator,a pole-zero cancellation net,and a linear amplification stage,which can be adapted to accommodate either positive or negative input signals.The RMS equivalent input noise charges are 3.3 fC,the conversion gains are approximately±2 mV∕fC,and the intrinsic time resolution reaches 32 ps.In the prototype PPAC application,the CSA performs as well as the commercial FTA820A amplifier,providing a position resolution as good as 0.17 mm,and exhibiting reliable stability during several hours of continuous data acquisition.
基金supported by Anhui Provincial Natural Science Foundation(2408085QA030)Natural Science Research Project of Anhui Educational Committee,China(2022AH050825)+3 种基金Medical Special Cultivation Project of Anhui University of Science and Technology(YZ2023H2C008)the Excellent Research and Innovation Team of Anhui Province,China(2022AH010052)the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology,China(2021yjrc51)Collaborative Innovation Program of Hefei Science Center,CAS,China(2019HSC-CIP006).
文摘In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical model,and its validity was verified using a simple impact test.A crushable discrete element method(DEM)framework is built based on the previously established theoretical model.The tensile strength,which considers the fractal theory,size effect,and Weibull variation,was assigned to each generated particle.The assigned strength is then used for crush detection by comparing it with its maximum tensile stress.Mass conservation is ensured by inserting a series of sub-particles whose total mass was equal to the quality loss.Based on the crushable DEM framework,a numerical simulation of the crushing behavior of a pebble bed with hollow cylindrical geometry under a uniaxial compression test was performed.The results of this investigation showed that the particle withstands the external load by contact and sliding at the beginning of the compression process,and the results confirmed that crushing can be considered an important method of resisting the increasing external load.A relatively regular particle arrangement aids in resisting the load and reduces the occurrence of particle crushing.However,a limit exists to the promotion of resistance.When the strain increases beyond this limit,the distribution of the crushing position tends to be isotropic over the entire pebble bed.The theoretical model and crushable DEM framework provide a new method for exploring the pebble bed in a fusion reactor,considering particle crushing.
基金financially supported by the National Natural Science Foundation of China(No.22273026)Scientific Research Innovation Capability Support Project for Young Faculty(No.ZYGXQNJSKYCXNLZCXM-I15)+3 种基金Basic and Applied Basic Research Foundation of Guangdong Province(2024A1515012401)GJYC program of Guangzhou(No.2024D03J0002)the China Postdoctoral Science Foundation(No.2024M750938)Postdoctoral Fellowship Program of CPSF(No.GZC20240492)for their financial support。
文摘In contrast to cyclic polymers with ring-like backbones,side-chain cyclization is another intriguing structural feature that has not been extensively studied.In this study,a library of orthogonally protected monomers featuring monocyclic,dicyclic,or tricyclic pendant motifs was designed and prepared based on malic acid derivatives.Polyesters with precise chemical structures and uniform chain lengths were prepared modularly through iterative growth.Meticulous control over the chemical details allows for a close investigation of the topological effects on the polymer properties.Compared to their linear side chain counterparts,the presence of cyclic pendant groups has a significant impact on chain conformation,leading to a reduction in hydrodynamic volume and an enhancement in the glass transition temperature.These results underscore the potential of tailoring polymer properties through rational engineering of side chain topology.
基金supported by the National Natural Science Foundation of China(Nos.22479092 and 22078190)。
文摘The pseudo-two-dimensional(P2D)model plays an important role in exploring physicochemical mechanisms,predicting the state of health,and improving the fast charge capability for Li-ion batteries(LIBs).However,the fast charge leads to the lithium concentration gradient in the solid and electrolyte phases and the non-uniform electrochemical reaction at the solid/electrolyte interface.In order to decouple charge transfer reactions in LIBs under dynamic conditions,understanding the spatio-temporal resolution of the P2D model is urgently required.Till now,the study of this aspect is still insufficient.This work studies the spatio-temporal resolution for dynamic/static electrochemical impedance spectroscopy(DEIS/SEIS)on multiple scales.In detail,DEIS and SEIS with spatio-temporal resolutions are used to decouple charge transfer reactions in LIBs based on the numerical solution of the P2D model in the frequency domain.The calculated results indicate that decoupling solid diffusion requires a high spatial resolution along the r-direction in particles,decoupling electrolyte diffusion and interfacial transfer reaction requires a high spatial resolution along the x-direction,and decoupling charge transfer reactions in LIBs at an extremely low state of charge(SOC)requires an extremely high temporal resolution along the t-direction.Finally,the optimal range of spatio-temporal resolutions for DEIS/SEIS is derived,and the method to decouple charge transfer reactions with spatio-temporal resolutions is developed.
基金Project(42072305)supported by the National Natural Science Foundation of China。
文摘Seabed mining operations have been found to induce significant movement and deformation in overlying rock strata,posing serious threats to mining safety.The presence of geological faults further complicates these deformation patterns.This study utilized geophysical surveys and the continuum-based discrete element method(CDEM)to investigate how fault activity influences rock deformation and failure.The results demonstrate that:1)Acting in mechanically weak zones,faults exerted a pronounced barrier effect on deformation propagation and stress redistribution within the surrounding rock,leading to markedly divergent displacement patterns on either side of the fault plane.Comparative analyses between single-fault and double-fault models revealed an 18%−22%expansion of the damage zone under the latter,together with significantly intensified deformation and failure;2)The double-fault model exhibited a larger maximum cumulative vertical displacement and a spatial shift in the location of peak deformation,thereby posing a heightened threat to mine safety;3)Acting in an orebody substitute,backfill effectively constrained surrounding rock deformation,enhanced its load-bearing capacity,and delayed the overburden subsidence.Nevertheless,backfill only reduced the amplitude of deformation;it could not entirely prevent settlement.These findings provide essential theoretical insights and foundational knowledge for safer submarine mining practices.
基金supported by the Japan Society for the Promotion of Science(JSPS)KAKENHI(Grant Nos.JP23KK0182,JP23K26356,and JP24K00971).
文摘Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.
基金partially supported by the National Natural Science Foundation of China(Grant No.12071073)financial support by the Jiangsu Provincial Scientific Research Center of Applied Mathematics(Grant No.BK20233002).
文摘The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearities.To discretize the continuous problem in space,we utilize a spectral Galerkin method.Subsequently,we introduce a nonlinear-tamed exponential integrator scheme,resulting in a fully discrete scheme.Within the framework of semigroup theory,this study provides precise estimations of the Sobolev regularity,L^(∞) regularity in space,and Hölder continuity in time for the mild solution,as well as for its semi-discrete and full-discrete approximations.Building upon these results,we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions.A numerical example is presented to validate the theoretical findings.
基金supported by the NSFC(12301115)the Natural Science Foundation of Huzhou(2023YZ11,2024YZ37)the second author was supported by the NSFC(12071437).
文摘It is well known that the inhomogeneous Calderón-Zygmund convolution operators are bounded on the local Hardy spaces.In this paper,we prove that these operators are bounded on the local product Hardy spaces and the Lipschitz spaces.The key ideas used here are the discrete local Calderón identity and a density argument for the inhomogeneous product Lipschitz spaces in the weak sense.
基金supported by the National Natural Science Foundation of China(Grant Nos.52574103 and 42277150).
文摘Fractures are typically characterized by roughness that significantlyaffects the mechanical and hydraulic characteristics of reservoirs.However,hydraulic fracturing mechanisms under the influenceof fracture morphology remain largely unexplored.Leveraging the advantages of the finite-discrete element method(FDEM)for explicitly simulating fracture propagation and the strengths of the unifiedpipe model(UPM)for efficientlymodeling dual-permeability seepage,we propose a new hydromechanical(HM)coupling approach for modeling hydraulic fracturing.Validated against benchmark examples,the proposed FDEM-UPM model is further augmented by incorporating a Fourier-based methodology for reconstructing non-planar fractures,enabling quantitative analysis of hydraulic fracturing behavior within rough discrete fracture networks(DFNs).The FDEM-UPM model demonstrates computational advantages in accurately capturing transient hydraulic seepage phenomena,while the asynchronous time-stepping schemes between hydraulic and mechanical analyses substantially enhanced computational efficiencywithout compromising computational accuracy.Our results show that fracture morphology can affect both macroscopic fracture networks and microscopic interaction types between hydraulic fractures(HFs)and natural fractures(NFs).In an isotropic stress field,the initiation azimuth,propagation direction and microcracking mechanism are significantly influencedby fracture roughness.In an anisotropic stress field,HFs invariably propagate parallel to the direction of the maximum principal stress,reducing the overall complexity of the stimulated fracture networks.Additionally,stress concentration and perturbation attributed to fracture morphology tend to be compromised as the leak-off increases,while the breakdown and propagation pressures remain insensitive to fracture morphology.These findingsprovide new insights into the hydraulic fracturing mechanisms of fractured reservoirs containing complex rough DFNs.
文摘Zero-dimensional(0D)organic-inorganic metal halide perovskite is one of the hot research topics in the field of optoelectronic materials.Their structure generally consists of discrete metal halide octahedra entirely isolated by surrounding organic cations,forming independent luminescent centers[1,2].Such a configuration results in high exciton binding energy and exceptional luminescence efficiency,due to strong quantum confinement[3,4].
基金This work was supported by the National Natural Science Foundation of China （Grant No. 10425105） and subsidized by the Special Funds for Major State Basic Research Projects （Grant No. 2005CB321704）.
文摘Some two-scale finite element discretizations are introduced for a class of linear partial differential equations. Both boundary value and eigenvalue problems are studied. Based on the two-scale error resolution techniques, several two-scale finite element algorithms are proposed and analyzed. It is shown that this type of two-scale algorithms not only significantly reduces the number of degrees of freedom but also produces very accurate approximations.
基金Research partially supported by NNSFC grant 10371118,SRF for ROCS,SEM and Nanjing University Talent Development Foundation.
文摘In this paper, we use Hermite weighted essentially non-oscillatory （HWENO） schemes with a Lax-Wendroff time discretization procedure, termed HWENO-LW schemes, to solve Hamilton-Jacobi equations. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first derivative values are evolved in time and are used in the reconstruction. One major advantage of HWENO schemes is its compactness in the reconstruction. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with HWENO with Runge-Kutta time discretizations schemes （HWENO-RK） of Qiu and Shu [19] for Hamilton-Jacobi equations, the major advantages of HWENO-LW schemes are their saving of computational cost and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method.
基金supported by National Natural Science Foundation of China(Grant Nos.11601241,11671199,11571290 and 11672082)Natural Science Foundation of Jiangsu Province(Grant No.BK20160877)+1 种基金ARO(Grant No.W911NF-15-1-0226)National Science Foundation of USA(Grant No.DMS-1719410)
文摘In this paper,we discuss the local discontinuous Galerkin methods coupled with two specific explicitimplicit-null time discretizations for solving one-dimensional nonlinear diffusion problems Ut=(a(U)Ux)x.The basic idea is to add and subtract two equal terms a0 Uxx the right-hand side of the partial differential equation,then to treat the term a0 Uxx implicitly and the other terms(a(U)Ux)x-a0 Uxx explicitly.We give stability analysis for the method on a simplified model by the aid of energy analysis,which gives a guidance for the choice of a0,i.e.,a0≥max{a(u)}/2 to ensure the unconditional stability of the first order and second order schemes.The optimal error estimate is also derived for the simplified model,and numerical experiments are given to demonstrate the stability,accuracy and performance of the schemes for nonlinear diffusion equations.
基金The authors were partially supported by the Special Project on High-performance Computing under the National Key R&D Program(No.2016YF B0200603)Sci-ence Challenge Project(No.JCK Y2016212A502)the National Natural Science Foundation of China(Nos.91630310&11421101).
文摘This paper studies the two-stage fourth-order accurate time discretization[J.Q.Li and Z.F.Du,SIAM J.Sci.Comput.,38(2016)]and its application to the special relativistic hydrodynamical equations.Our analysis reveals that the new two-stage fourth-order accurate time discretizations can be proposed.With the aid of the direct Eulerian GRP(generalized Riemann problem)methods and the analytical resolution of the local“quasi 1D”GRP,the two-stage fourth-order accurate time discretizations are successfully implemented for the 1D and 2D special relativistic hydrodynamical equations.Several numerical experiments demonstrate the performance and accuracy as well as robustness of our schemes.
基金Project supported by the Open Research Fund of the State Key Laboratory of Water Resources and Hydropower Engineering Science, and the National Natural Science Foundation of China (Grant No: 50459001).
文摘Upwind algorithms are becoming progressively popular for river flood routing due to their capability of resolving trans-critical flow regimes. For consistency, these algorithms suggest natural upwind discretization of the source term, which may be essential for natural channels with irregular geometry. Yet applications of these upwind algorithms to natural river flows are rare, and in such applications the traditional and simpler pointwise, rather than upwind discretization of the source term is used. Within the framework of a first-order upwind algorithm, this paper presents a comparison of upwind and pointwise discretizations of the source term. Numerical simulations were carried out for a selected irregular channel comprising a pool-riffle sequence Jn the River Lune, England with observed data. It is Shown that the impact of pointwise discretization, compared to the upwind, is appreciable mainly in flow zones with the Froude number closer to or larger than unity. The discrepancy due to pointwise and upwind discretizations of the source term is negligible in flow depth and hence in water surface elevation, but well manifested in mean velocity and derived flow quantities. Also the occurrence of flow reversal and equalisation over the pool-riffle sequence in response to increasing discharges is demonstrated.
文摘Integrable discretizations of are proposed. N-soliton solutions for analogues of the complex and real Dym the complex and real Dym equations both semi-discrete and fully discrete equations are also presented.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金Australian Research Council Linkage Program(LP200301404)for sponsoring this researchthe financial support provided by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology,SKLGP2021K002)National Natural Science Foundation of China(52374101,32111530138).
文摘Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.
