In this manuscript,we propose an analytical equivalent linear viscoelastic constitutive model for fiber-reinforced composites,bypassing general computational homogenization.The method is based on the reduced-order hom...In this manuscript,we propose an analytical equivalent linear viscoelastic constitutive model for fiber-reinforced composites,bypassing general computational homogenization.The method is based on the reduced-order homogenization(ROH)approach.The ROH method typically involves solving multiple finite element problems under periodic conditions to evaluate elastic strain and eigenstrain influence functions in an‘off-line’stage,which offers substantial cost savings compared to direct computational homogenization methods.Due to the unique structure of the fibrous unit cell,“off-line”stage calculation can be eliminated by influence functions obtained analytically.Introducing the standard solid model to the ROH method enables the creation of a comprehensive analytical homogeneous viscoelastic constitutive model.This method treats fibrous composite materials as homogeneous,anisotropic viscoelastic materials,significantly reducing computational time due to its analytical nature.This approach also enables precise determination of a homogenized anisotropic relaxation modulus and accurate capture of various viscoelastic responses under different loading conditions.Three sets of numerical examples,including unit cell tests,three-point beam bending tests,and torsion tests,are given to demonstrate the predictive performance of the homogenized viscoelastic model.Furthermore,the model is validated against experimental measurements,confirming its accuracy and reliability.展开更多
Hydrodynamic instability growth at the deuterium-tritium(DT)fuel-ablator interface plays a critical role in determining the performance of inertial confinement fusion implosions.During the late stages of implosion,ins...Hydrodynamic instability growth at the deuterium-tritium(DT)fuel-ablator interface plays a critical role in determining the performance of inertial confinement fusion implosions.During the late stages of implosion,insufficient doping of the ablator material can result in highenergy X-ray preheat,which may trigger the development of a classical-like Rayleigh-Taylor instability(RTI)at the fuel-ablator interface.In implosion experiments at the Shenguang 100 kJ-level laser facility,the primary source of perturbation is the roughness of the inner DT ice interface.In this study,we propose an analytical model to describe the feed-out process of the initial roughness of the inner DT ice interface.The perturbation amplitude derived from this model serves as the initial seed for the late-time RTI during the acceleration phase.Our findings confirm the presence of classical-like RTI at the fuel-ablator interface.Numerical simulations conducted using a radiation hydrodynamic code validate the proposed analytical model and demonstrate the existence of a peak mode number in both the feed-out process and the classical-like RTI.It provides an alternative bridge between the current target fabrication limitations and the unexpected implosion performance.展开更多
Inverse bremsstrahlung absorption in laser-heated plasmas is studied using the Fokker–Planck equation in the low-field limit.Compared with the commonly used fitting formulas of Langdon and Matte et al.,our work emplo...Inverse bremsstrahlung absorption in laser-heated plasmas is studied using the Fokker–Planck equation in the low-field limit.Compared with the commonly used fitting formulas of Langdon and Matte et al.,our work employs fewer approximations and provides more accurate predictions for the super-Gaussian orderβand the heating rate.Simulation results show that the super-Gaussian order is generally lower than the fitting results of Matte et al.,which leads to an increase in absorption.However,we find two other factors that reduce absorption:the high-order term of the collision frequency and the effects caused by high laser intensity.Therefore,the final simulated absorption can either be higher or lower,depending on the conditions.These phenomena are theoretically analyzed using the Fokker–Planck equation.Fitting formulas are proposed for the super-Gaussian order and the heating rate,showing a discrepancy within∼10%of the simulation results.We also compare the simulation results with the experimental results from several recent papers.展开更多
We present a fully time-dependent quantum wave packet evolution method for investigating molecular dynamics in intense laser fields.This approach enables the simultaneous treatment of interactions among multiple elect...We present a fully time-dependent quantum wave packet evolution method for investigating molecular dynamics in intense laser fields.This approach enables the simultaneous treatment of interactions among multiple electronic states while simultaneously tracking their time-dependent electronic,vibrational,and rotational dynamics.As an illustrative example,we consider neutral H_(2)molecules and simulate the laser-induced excitation dynamics of electronic and rotational states in strong laser fields,quantitatively distinguishing the respective contributions of electronic dipole transitions(within the classical-field approximation)and non-resonant Raman processes to the overall molecular dynamics.Furthermore,we precisely evaluate the relative contributions of direct tunneling ionization from the ground state and ionization following electronic excitation in the strong-field ionization of H_(2).The developed methodology shows strong potential for performing high-precision theoretical simulations of electronic-vibrational-rotational state excitations,ionization,and dissociation dynamics in molecules and their ions under intense laser fields.展开更多
By adopting stochastic density functional theory(SDFT)and mixed stochastic-deterministic density functional theory(MDFT)methods,we perform first-principles calculations to predict the shock Hugoniot curves of boron(pr...By adopting stochastic density functional theory(SDFT)and mixed stochastic-deterministic density functional theory(MDFT)methods,we perform first-principles calculations to predict the shock Hugoniot curves of boron(pressure P=7.9×10^(3)-1.6×10^(6) GPa and temperature T=25-2800 eV),silicon(P=2.6×10^(3)-7.9×10^(5) GPa and T=21.5-1393 eV),and aluminum(P=5.2×10^(3)-9.0×10^(5) GPa and T=25-1393 eV)over wide ranges of pressure and temperature.In particular,we systematically investigate the impact of different cutoff radii in norm-conserving pseudopotentials on the calculated properties at elevated temperatures,such as pressure,ionization energy,and equation of state.By comparing the SDFT and MDFT results with those of other first-principles methods,such as extended first-principles molecular dynamics and path integral Monte Carlo methods,we find that the SDFT and MDFT methods show satisfactory precision,which advances our understanding of first-principles methods when applied to studies of matter at extremely high pressures and temperatures.展开更多
The multiscale computational method with asymptotic analysis and reduced-order homogenization(ROH)gives a practical numerical solution for engineering problems,especially composite materials.Under the ROH framework,a ...The multiscale computational method with asymptotic analysis and reduced-order homogenization(ROH)gives a practical numerical solution for engineering problems,especially composite materials.Under the ROH framework,a partition-based unitcell structure at the mesoscale is utilized to give a mechanical state at the macro-scale quadrature point with pre-evaluated influence functions.In the past,the“1-phase,1-partition”rule was usually adopted in numerical analysis,where one constituent phase at the mesoscale formed one partition.The numerical cost then is significantly reduced by introducing an assumption that the mechanical responses are the same all the time at the same constituent,while it also introduces numerical inaccuracy.This study proposes a new partitioning method for fibrous unitcells under a reduced-order homogenization methodology.In this method,the fiber phase remains 1 partition,but the matrix phase is divided into 2 partitions,which refers to the“12”partitioning scheme.Analytical elastic influence+functions are derived by introducing the elastic strain energy equivalence(Hill-Mandel condition).This research also obtains the analytical eigenstrain influence functions by alleviating the so-called“inclusion-locking”phenomenon.In addition,a numerical approach to minimize the error of strain energy density is introduced to determine the partitioning of the matrix phase.Several numerical examples are presented to compare the differences among direct numerical simulation(DNS),“11”,and“12”partitioning schemes.The numerical simulations show improved++numerical accuracy by the“12”partitioning scheme.展开更多
This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solv...This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.展开更多
The properties of the non-trivial quantum state in an all-optical environment come mainly from the higher-order quantum electrodynamics effect,which remains one of the few unverified predictions of this theory due to ...The properties of the non-trivial quantum state in an all-optical environment come mainly from the higher-order quantum electrodynamics effect,which remains one of the few unverified predictions of this theory due to its weak signal.Here,we propose a scheme specifically designed to detect this quantum vacuum,where a tightly focused pump laser interacts with an optical frequency comb(OFC)in its resonant cavity.When the OFC pulse passes through the vacuum polarized by the high-intensity pump laser,its carrier frequency and envelope change.This can be intuitively understood as the asymmetric photon acceleration induced by the ponderomotive force of the pump laser.By leveraging the exceptional ultrahigh frequency and temporal resolution of the OFC,this scheme holds the potential to improve the accuracy of quantum vacuum signal.Combining theoretical and simulation results,we discuss possible experimental conditions,and the detectable OFC signal is shown to be orders of magnitude better than the instrumental detection threshold.This shows our scheme can be verified on the forthcoming laser systems.展开更多
Large-angle stimulated Raman scattering(LA-SRS)in a longitudinally inhomogeneous plasma with a transverse density modulation is studied using a three-wave coupled model and numerical simulations.The simulations show t...Large-angle stimulated Raman scattering(LA-SRS)in a longitudinally inhomogeneous plasma with a transverse density modulation is studied using a three-wave coupled model and numerical simulations.The simulations show that the scattering angle of SRS in a longitudinally inhomogeneous plasma can be significantly affected by transverse density modulation.Under transverse density modulation conditions,the laser focuses into underdense regions,owing to the transversely modulated refractive index.The angle of LA-SRS,neither a purely 90° angle side scattering nor purely backscattering,is almost consistent with the specific angle at which the density inhomogeneity vanishes.In modulated plasmas,the nonuniform distribution of laser intensity shifts the regions of scattering and gain compared with those in uniform plasmas,ultimately affecting the laser transmission.SRS is suppressed in weakly modulated regimes,whereas it is enhanced under strong modulation conditions,and a theoretical criterion distinguishing between strong and weak modulation is established.展开更多
We put forward a new design of a compact beam transport system for intense laser-driven proton therapy,where instead of using conventional pulsed solenoids,our design relies on a helical coil irradiated by a nanosecon...We put forward a new design of a compact beam transport system for intense laser-driven proton therapy,where instead of using conventional pulsed solenoids,our design relies on a helical coil irradiated by a nanosecond laser pulse to generate strong magnetic fields for focusing protons.A pair of dipole magnets and apertures are employed to further filter protons with large divergences and low energies.Our numerical studies combine particle-in-cell simulations for laser-plasma interaction to generate high-energy monoenergetic proton beams,finite element analysis for evaluating the magnetic field distribution inside the coil,and MonteCarlo simulations for beam transport and energy deposition.Our results show that with this design,a spread-out Bragg peak in a range of several centimeters to a deep-seated tumor with a dose of approximately 16.5 cGy and fluctuation around 2% can be achieved.The instantaneous dose rate reaches up to 10^(9)Gy/s,holding the potential for future FLASH radiotherapy research.展开更多
Terahertz(THz)radiation is rapidly emerging as a powerful tool with diverse applications,including high-speed imaging,laser-driven particle acceleration,and ultra-high frequency(UHF)communications.However,generating m...Terahertz(THz)radiation is rapidly emerging as a powerful tool with diverse applications,including high-speed imaging,laser-driven particle acceleration,and ultra-high frequency(UHF)communications.However,generating multipulse THz radiation with controllable time intervals remains a significant challenge.This study presents an approach to overcome this hurdle by exploiting the interaction between an electron beam and plasma.Using numerical simulations and theoretical analysis,we investigated the behavior of an electron beam within a plasma and its interaction with the longitudinal sheath field.This interaction resulted in the generation of multiple distinct THz pulses.We demonstrated that the plasma length adjustment allows for precise tuning of the interval between THz pulses.Moreover,the radiation intensity could be controlled by the electron beam energy and the electron bunch duration.The proposed scheme can generate multipulse THz radiation in a flexible and precise manner,paving the way for advancements in applications requiring high temporal resolution.展开更多
Terahertz(THz) radiation, spanning the frequency range 100 GHz to 10 THz, offers diverse applications in spectroscopy, materials characterization, medical diagnostics and environmental monitoring. Despite its potentia...Terahertz(THz) radiation, spanning the frequency range 100 GHz to 10 THz, offers diverse applications in spectroscopy, materials characterization, medical diagnostics and environmental monitoring. Despite its potential, the generation of high-intensity, tunable THz radiation remains a significant challenge. In this work, we explore a novel approach to the efficient generation of THz radiation based on laser–plasma interactions, utilizing the principles of photon deceleration.When a relativistic CO_(2) laser passes through a pre-ionized plasma, the laser induces a nonlinear wakefield, creating a strong refractive index gradient. This gradient, combined with the lower-density region of the wakefield, slows down the laser, facilitating the accumulation of THz radiation. The resulting THz pulse exhibits extreme collimation, high energy efficiency and tunability. Our work shows that this method can achieve up to 10% conversion efficiency with optimal plasma density near the critical density. This technique presents a promising solution for overcoming current limitations in THz source development and offers potential for diverse applications.展开更多
A Discrete Boltzmann Method(DBM)with a Maxwell-type boundary condition is constructed to investigate the influence of rarefaction on laminar Shock Wave/Boundary Layer Interaction(SWBLI).Due to the complexity of compre...A Discrete Boltzmann Method(DBM)with a Maxwell-type boundary condition is constructed to investigate the influence of rarefaction on laminar Shock Wave/Boundary Layer Interaction(SWBLI).Due to the complexity of compressible flow,a Knudsen number vector Kn,whose components include the local Knudsen numbers such as Kn_(ρ)and Kn_(U),is introduced to characterize the local structures,where Kn_(ρ)and Kn_(U)are Knudsen numbers defined in terms of the density and velocity interfaces,respectively.Since first focusing on the steady state of SWBLI,the DBM considers up to the second-order Kn_(ρ)(rarefaction/non-equilibrium)effects.The model is validated using Mach number 2 SWBLI and the necessity of using DBM with sufficient physical accuracy is confirmed by the shock collision problem.Key findings include the following:the leading-edge shock wave increases the local density Knudsen number Kn_(ρ)and eventually leads to the failure of linear constitutive relations in the Navier-Stokes(N-S)model and surely also in the lower-order DBM;the non-equilibrium effect differences in regions behind the leading-edge shock wave are primarily correlated with Kn_(ρ),while in the separation region are primarily correlated with Kn_(U);the non-equilibrium quantities D_(2)and D_(4,2),as well as the viscous entropy production rate S_(NOMF)can be used to identify the separation zone.The findings clarify various effects and main mechanisms in different regions associated with SWBLI,which are concealed in N-S model.展开更多
Controlling terahertz(THz)polarization with high stability and tunability is essential for achieving further progress in ultrafast spectroscopy,structured-light manipulation,and quantum information processing.Here,we ...Controlling terahertz(THz)polarization with high stability and tunability is essential for achieving further progress in ultrafast spectroscopy,structured-light manipulation,and quantum information processing.Here,we propose a magnetized plasma platform for dynamic THz polarization control by exploiting the intrinsic birefringence between extraordinary and ordinary modes.We identify a strong-magnetization,zero-group-velocity-mismatch regime where the two modes share matched group velocities while retaining finite phase birefringence,enabling robust,phase-stable spin angular momentum control.By tuning the plasma length and magnetic field,we realize programmable phase retardation and demonstrate universal single-qubit gates through parameterized unitary operations.Full-wave particle-in-cell simulations validate high-fidelity polarization transformations across the Poincarésphere and demonstrate the potential for generating structured vector beams under spatially varying magnetic fields.The platform offers ultrafast response,resilience to extreme THz intensities,and in situ tunability,positioning magnetized plasmas as a versatile and damage-resilient medium for next-generation THz polarization control and structured-wave applications.展开更多
X-ray free-electron lasers(XFELs)can generate bright X-ray pulses with short durations and narrow bandwidths,leading to extensive applica-tions in many disciplines such as biology,materials science,and ultrafast scien...X-ray free-electron lasers(XFELs)can generate bright X-ray pulses with short durations and narrow bandwidths,leading to extensive applica-tions in many disciplines such as biology,materials science,and ultrafast science.Recently,there has been a growing demand for X-ray pulses with high photon energy,especially from developments in“diffraction-before-destruction”applications and in dynamic mesoscale materials science.Here,we propose utilizing the electron beams at XFELs to drive a meter-scale two-bunch plasma wakefield accelerator and double the energy of the accelerated beam in a compact and inexpensive way.Particle-in-cell simulations are performed to study the beam quality degradation under different beam loading scenarios and nonideal issues,and the results show that more than half of the accelerated beam can meet the requirements of XFELs.After its transport to the undulator,the accelerated beam can improve the photon energy to 22 keV by a factor of around four while maintaining the peak power,thus offering a promising pathway toward high-photon-energy XFELs.展开更多
Correction to:Nuclear Science and Techniques(2025)36:100 https://doi.org/10.1007/s41365-025-01692-6 In this article,Fig.9 appeared incorrectly and have now been corrected in the original publication.For completeness a...Correction to:Nuclear Science and Techniques(2025)36:100 https://doi.org/10.1007/s41365-025-01692-6 In this article,Fig.9 appeared incorrectly and have now been corrected in the original publication.For completeness and transparency,both correct and incorrect versions are displayed below.展开更多
High-power laser pulses interacting with targets can generate intense electromagnetic pulses(EMPs),which can disrupt physical experimental diagnostics and even damage diagnostic equipment,posing a threat to the reliab...High-power laser pulses interacting with targets can generate intense electromagnetic pulses(EMPs),which can disrupt physical experimental diagnostics and even damage diagnostic equipment,posing a threat to the reliable operation of experiments.In this study,EMPs resulting from multi-petawatt laser irradiating nitrogen gas jets were systematically analyzed and investigated.The experimental results revealed that the EMP amplitude is positively correlated with the quantity and energy of the electrons captured and accelerated by the plasma channel.These factors are reflected by parameters such as laser energy and nitrogen gas jet pressure.Additionally,we propose several potential sources of EMPs produced by laser-irradiated gas jets and separately analyzed their spatiotemporal distributions.The findings provide insight into the mechanisms of EMP generation and introduce a new approach to achieve controllable EMPs by regulating the laser energy and gas jet pressure.展开更多
基金support by the National Key R&D Program of China(Grant No.2023YFA1008901)the National Natural Science Foundation of China(Grant Nos.11988102,12172009)is gratefully acknowledged.
文摘In this manuscript,we propose an analytical equivalent linear viscoelastic constitutive model for fiber-reinforced composites,bypassing general computational homogenization.The method is based on the reduced-order homogenization(ROH)approach.The ROH method typically involves solving multiple finite element problems under periodic conditions to evaluate elastic strain and eigenstrain influence functions in an‘off-line’stage,which offers substantial cost savings compared to direct computational homogenization methods.Due to the unique structure of the fibrous unit cell,“off-line”stage calculation can be eliminated by influence functions obtained analytically.Introducing the standard solid model to the ROH method enables the creation of a comprehensive analytical homogeneous viscoelastic constitutive model.This method treats fibrous composite materials as homogeneous,anisotropic viscoelastic materials,significantly reducing computational time due to its analytical nature.This approach also enables precise determination of a homogenized anisotropic relaxation modulus and accurate capture of various viscoelastic responses under different loading conditions.Three sets of numerical examples,including unit cell tests,three-point beam bending tests,and torsion tests,are given to demonstrate the predictive performance of the homogenized viscoelastic model.Furthermore,the model is validated against experimental measurements,confirming its accuracy and reliability.
基金funded by the National Key R&D Program of China(Grant No.2023YFA1608400)the National Natural Science Foundation of China(Grant No.12302281).
文摘Hydrodynamic instability growth at the deuterium-tritium(DT)fuel-ablator interface plays a critical role in determining the performance of inertial confinement fusion implosions.During the late stages of implosion,insufficient doping of the ablator material can result in highenergy X-ray preheat,which may trigger the development of a classical-like Rayleigh-Taylor instability(RTI)at the fuel-ablator interface.In implosion experiments at the Shenguang 100 kJ-level laser facility,the primary source of perturbation is the roughness of the inner DT ice interface.In this study,we propose an analytical model to describe the feed-out process of the initial roughness of the inner DT ice interface.The perturbation amplitude derived from this model serves as the initial seed for the late-time RTI during the acceleration phase.Our findings confirm the presence of classical-like RTI at the fuel-ablator interface.Numerical simulations conducted using a radiation hydrodynamic code validate the proposed analytical model and demonstrate the existence of a peak mode number in both the feed-out process and the classical-like RTI.It provides an alternative bridge between the current target fabrication limitations and the unexpected implosion performance.
基金supported by the Natural Science Foundation of China(Grant Nos.12035002,12405235,and U2430207)the China Postdoctoral Science Foundation(Grant No.2023M740336)the CAEP foundation(Grant No.YZJJZQ2023020).
文摘Inverse bremsstrahlung absorption in laser-heated plasmas is studied using the Fokker–Planck equation in the low-field limit.Compared with the commonly used fitting formulas of Langdon and Matte et al.,our work employs fewer approximations and provides more accurate predictions for the super-Gaussian orderβand the heating rate.Simulation results show that the super-Gaussian order is generally lower than the fitting results of Matte et al.,which leads to an increase in absorption.However,we find two other factors that reduce absorption:the high-order term of the collision frequency and the effects caused by high laser intensity.Therefore,the final simulated absorption can either be higher or lower,depending on the conditions.These phenomena are theoretically analyzed using the Fokker–Planck equation.Fitting formulas are proposed for the super-Gaussian order and the heating rate,showing a discrepancy within∼10%of the simulation results.We also compare the simulation results with the experimental results from several recent papers.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFA1602502)the National Natural Science Foundation of China(Grant No.12450404)。
文摘We present a fully time-dependent quantum wave packet evolution method for investigating molecular dynamics in intense laser fields.This approach enables the simultaneous treatment of interactions among multiple electronic states while simultaneously tracking their time-dependent electronic,vibrational,and rotational dynamics.As an illustrative example,we consider neutral H_(2)molecules and simulate the laser-induced excitation dynamics of electronic and rotational states in strong laser fields,quantitatively distinguishing the respective contributions of electronic dipole transitions(within the classical-field approximation)and non-resonant Raman processes to the overall molecular dynamics.Furthermore,we precisely evaluate the relative contributions of direct tunneling ionization from the ground state and ionization following electronic excitation in the strong-field ionization of H_(2).The developed methodology shows strong potential for performing high-precision theoretical simulations of electronic-vibrational-rotational state excitations,ionization,and dissociation dynamics in molecules and their ions under intense laser fields.
基金supported by the National Key R&D Program of China under Grant No.2025YFB3003603the National Natural Science Foundation of China under Grant Nos.12135002 and 12105209.
文摘By adopting stochastic density functional theory(SDFT)and mixed stochastic-deterministic density functional theory(MDFT)methods,we perform first-principles calculations to predict the shock Hugoniot curves of boron(pressure P=7.9×10^(3)-1.6×10^(6) GPa and temperature T=25-2800 eV),silicon(P=2.6×10^(3)-7.9×10^(5) GPa and T=21.5-1393 eV),and aluminum(P=5.2×10^(3)-9.0×10^(5) GPa and T=25-1393 eV)over wide ranges of pressure and temperature.In particular,we systematically investigate the impact of different cutoff radii in norm-conserving pseudopotentials on the calculated properties at elevated temperatures,such as pressure,ionization energy,and equation of state.By comparing the SDFT and MDFT results with those of other first-principles methods,such as extended first-principles molecular dynamics and path integral Monte Carlo methods,we find that the SDFT and MDFT methods show satisfactory precision,which advances our understanding of first-principles methods when applied to studies of matter at extremely high pressures and temperatures.
基金funded by the National Key R&D Program of China(Grant No.2023YFA1008901)the National Natural Science Foundation of China(Grant Nos.11988102,12172009)“The Fundamental Research Funds for the Central Universities,Peking University”.
文摘The multiscale computational method with asymptotic analysis and reduced-order homogenization(ROH)gives a practical numerical solution for engineering problems,especially composite materials.Under the ROH framework,a partition-based unitcell structure at the mesoscale is utilized to give a mechanical state at the macro-scale quadrature point with pre-evaluated influence functions.In the past,the“1-phase,1-partition”rule was usually adopted in numerical analysis,where one constituent phase at the mesoscale formed one partition.The numerical cost then is significantly reduced by introducing an assumption that the mechanical responses are the same all the time at the same constituent,while it also introduces numerical inaccuracy.This study proposes a new partitioning method for fibrous unitcells under a reduced-order homogenization methodology.In this method,the fiber phase remains 1 partition,but the matrix phase is divided into 2 partitions,which refers to the“12”partitioning scheme.Analytical elastic influence+functions are derived by introducing the elastic strain energy equivalence(Hill-Mandel condition).This research also obtains the analytical eigenstrain influence functions by alleviating the so-called“inclusion-locking”phenomenon.In addition,a numerical approach to minimize the error of strain energy density is introduced to determine the partitioning of the matrix phase.Several numerical examples are presented to compare the differences among direct numerical simulation(DNS),“11”,and“12”partitioning schemes.The numerical simulations show improved++numerical accuracy by the“12”partitioning scheme.
基金supported by the National Natural Science Foundation of China Excellence Research Group Program for“Multiscale Problems in Nonlinear Mechanics”(Grant No.12588201)the National Key R&D Program of China(Grant No.2023YFA1008901)+1 种基金the National Nat-ural Science Foundation of China(Grant No.12172009)supported by“The Fundamental Research Funds for the Central Universities,Peking University”.
文摘This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.
基金supported by the National Key R&D Program of China(Grant Nos.2022YFA1603200,2022YFA1603201,2024YFA1613400)the National Natural Science Foundation of China(Grant Nos.12135001,11825502,12075014,12475243)+2 种基金the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA25050900)the Science and Technology on Plasma Physics Laboratory(Grant No.6142A04210110)the National Natural Science Funds for Distinguished Young Scholars(Grant No.11825502)。
文摘The properties of the non-trivial quantum state in an all-optical environment come mainly from the higher-order quantum electrodynamics effect,which remains one of the few unverified predictions of this theory due to its weak signal.Here,we propose a scheme specifically designed to detect this quantum vacuum,where a tightly focused pump laser interacts with an optical frequency comb(OFC)in its resonant cavity.When the OFC pulse passes through the vacuum polarized by the high-intensity pump laser,its carrier frequency and envelope change.This can be intuitively understood as the asymmetric photon acceleration induced by the ponderomotive force of the pump laser.By leveraging the exceptional ultrahigh frequency and temporal resolution of the OFC,this scheme holds the potential to improve the accuracy of quantum vacuum signal.Combining theoretical and simulation results,we discuss possible experimental conditions,and the detectable OFC signal is shown to be orders of magnitude better than the instrumental detection threshold.This shows our scheme can be verified on the forthcoming laser systems.
基金supported by the National Natural Science Foundation of China under Grant Nos.U2430207,12035002,and 12305258by the CAEP Foundation under Grant No.YZJJZQ2023020.
文摘Large-angle stimulated Raman scattering(LA-SRS)in a longitudinally inhomogeneous plasma with a transverse density modulation is studied using a three-wave coupled model and numerical simulations.The simulations show that the scattering angle of SRS in a longitudinally inhomogeneous plasma can be significantly affected by transverse density modulation.Under transverse density modulation conditions,the laser focuses into underdense regions,owing to the transversely modulated refractive index.The angle of LA-SRS,neither a purely 90° angle side scattering nor purely backscattering,is almost consistent with the specific angle at which the density inhomogeneity vanishes.In modulated plasmas,the nonuniform distribution of laser intensity shifts the regions of scattering and gain compared with those in uniform plasmas,ultimately affecting the laser transmission.SRS is suppressed in weakly modulated regimes,whereas it is enhanced under strong modulation conditions,and a theoretical criterion distinguishing between strong and weak modulation is established.
基金supported by the National Key R&D Program of China(Nos.2022YFA1603200 and 2022YFA1603201)National Natural Science Foundation of China(Nos.12135001,11921006,12475243 and 11825502)+1 种基金Strategic Priority Research Program of CAS(No.XDA25050900)support from the National Natural Science Funds for Distinguished Young Scholar(No.11825502)。
文摘We put forward a new design of a compact beam transport system for intense laser-driven proton therapy,where instead of using conventional pulsed solenoids,our design relies on a helical coil irradiated by a nanosecond laser pulse to generate strong magnetic fields for focusing protons.A pair of dipole magnets and apertures are employed to further filter protons with large divergences and low energies.Our numerical studies combine particle-in-cell simulations for laser-plasma interaction to generate high-energy monoenergetic proton beams,finite element analysis for evaluating the magnetic field distribution inside the coil,and MonteCarlo simulations for beam transport and energy deposition.Our results show that with this design,a spread-out Bragg peak in a range of several centimeters to a deep-seated tumor with a dose of approximately 16.5 cGy and fluctuation around 2% can be achieved.The instantaneous dose rate reaches up to 10^(9)Gy/s,holding the potential for future FLASH radiotherapy research.
基金supported by the National Natural Science Foundation of China(Grant No.12175058)the National Science Fund of Hunan Province for Distinguished Young Scholars(Grant No.2024JJ2009)。
文摘Terahertz(THz)radiation is rapidly emerging as a powerful tool with diverse applications,including high-speed imaging,laser-driven particle acceleration,and ultra-high frequency(UHF)communications.However,generating multipulse THz radiation with controllable time intervals remains a significant challenge.This study presents an approach to overcome this hurdle by exploiting the interaction between an electron beam and plasma.Using numerical simulations and theoretical analysis,we investigated the behavior of an electron beam within a plasma and its interaction with the longitudinal sheath field.This interaction resulted in the generation of multiple distinct THz pulses.We demonstrated that the plasma length adjustment allows for precise tuning of the interval between THz pulses.Moreover,the radiation intensity could be controlled by the electron beam energy and the electron bunch duration.The proposed scheme can generate multipulse THz radiation in a flexible and precise manner,paving the way for advancements in applications requiring high temporal resolution.
基金Project supported by the China Postdoctoral Science Foundation (Grant No. 2024T170021)the Beijing Municipal Science & Technology Commission, Administrative Commission of Zhongguancun Science Park (Grant No. Z231100006023003)+2 种基金the National Natural Science Foundation of China (Grant Nos. 12175058, 12205007, and 11921006)the National Science Fund of Hunan Province for Distinguished Young Scholars (Grant No. 2024JJ2009)The computing was supported by the High-performance Computing Platform of Peking University。
文摘Terahertz(THz) radiation, spanning the frequency range 100 GHz to 10 THz, offers diverse applications in spectroscopy, materials characterization, medical diagnostics and environmental monitoring. Despite its potential, the generation of high-intensity, tunable THz radiation remains a significant challenge. In this work, we explore a novel approach to the efficient generation of THz radiation based on laser–plasma interactions, utilizing the principles of photon deceleration.When a relativistic CO_(2) laser passes through a pre-ionized plasma, the laser induces a nonlinear wakefield, creating a strong refractive index gradient. This gradient, combined with the lower-density region of the wakefield, slows down the laser, facilitating the accumulation of THz radiation. The resulting THz pulse exhibits extreme collimation, high energy efficiency and tunability. Our work shows that this method can achieve up to 10% conversion efficiency with optimal plasma density near the critical density. This technique presents a promising solution for overcoming current limitations in THz source development and offers potential for diverse applications.
基金support from the National Key R&D Program of China(No.2020YFC2201100)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics,China(No.JCKYS2023212003)+1 种基金the National Natural Science Foundation of China(No.12172061)the Opening Project of State Key Laboratory of Explosion Science and Safety Protection(Beijing Institute of Technology)(No.KFJJ25-02M).
文摘A Discrete Boltzmann Method(DBM)with a Maxwell-type boundary condition is constructed to investigate the influence of rarefaction on laminar Shock Wave/Boundary Layer Interaction(SWBLI).Due to the complexity of compressible flow,a Knudsen number vector Kn,whose components include the local Knudsen numbers such as Kn_(ρ)and Kn_(U),is introduced to characterize the local structures,where Kn_(ρ)and Kn_(U)are Knudsen numbers defined in terms of the density and velocity interfaces,respectively.Since first focusing on the steady state of SWBLI,the DBM considers up to the second-order Kn_(ρ)(rarefaction/non-equilibrium)effects.The model is validated using Mach number 2 SWBLI and the necessity of using DBM with sufficient physical accuracy is confirmed by the shock collision problem.Key findings include the following:the leading-edge shock wave increases the local density Knudsen number Kn_(ρ)and eventually leads to the failure of linear constitutive relations in the Navier-Stokes(N-S)model and surely also in the lower-order DBM;the non-equilibrium effect differences in regions behind the leading-edge shock wave are primarily correlated with Kn_(ρ),while in the separation region are primarily correlated with Kn_(U);the non-equilibrium quantities D_(2)and D_(4,2),as well as the viscous entropy production rate S_(NOMF)can be used to identify the separation zone.The findings clarify various effects and main mechanisms in different regions associated with SWBLI,which are concealed in N-S model.
基金supported by the National Natural Science Foundation of China (Grant Nos. 12175058 and 11921006)the National Grand Instrument Project (No. 2019YFF01014402)the Beijing Distinguished Young Scientist Program and National Grand Instrument Project No. SQ2019YFF01014400
文摘Controlling terahertz(THz)polarization with high stability and tunability is essential for achieving further progress in ultrafast spectroscopy,structured-light manipulation,and quantum information processing.Here,we propose a magnetized plasma platform for dynamic THz polarization control by exploiting the intrinsic birefringence between extraordinary and ordinary modes.We identify a strong-magnetization,zero-group-velocity-mismatch regime where the two modes share matched group velocities while retaining finite phase birefringence,enabling robust,phase-stable spin angular momentum control.By tuning the plasma length and magnetic field,we realize programmable phase retardation and demonstrate universal single-qubit gates through parameterized unitary operations.Full-wave particle-in-cell simulations validate high-fidelity polarization transformations across the Poincarésphere and demonstrate the potential for generating structured vector beams under spatially varying magnetic fields.The platform offers ultrafast response,resilience to extreme THz intensities,and in situ tunability,positioning magnetized plasmas as a versatile and damage-resilient medium for next-generation THz polarization control and structured-wave applications.
基金supported by the National Grand Instrument Project No. SQ2019YFF01014400the Natural Science Foundation of China (Grant Nos. 12375147, 12435011, 12075030)+2 种基金the Beijing Outstanding Young Scientist Project, Project for Young Scientists in Basic Research of Chinese Academy of Sciences (YSBR-115)the Beijing Normal University Scientific Research Initiation Fund for Introducing Talents No. 310432104the Fundamental Research Funds for the Central Universities, Peking University
文摘X-ray free-electron lasers(XFELs)can generate bright X-ray pulses with short durations and narrow bandwidths,leading to extensive applica-tions in many disciplines such as biology,materials science,and ultrafast science.Recently,there has been a growing demand for X-ray pulses with high photon energy,especially from developments in“diffraction-before-destruction”applications and in dynamic mesoscale materials science.Here,we propose utilizing the electron beams at XFELs to drive a meter-scale two-bunch plasma wakefield accelerator and double the energy of the accelerated beam in a compact and inexpensive way.Particle-in-cell simulations are performed to study the beam quality degradation under different beam loading scenarios and nonideal issues,and the results show that more than half of the accelerated beam can meet the requirements of XFELs.After its transport to the undulator,the accelerated beam can improve the photon energy to 22 keV by a factor of around four while maintaining the peak power,thus offering a promising pathway toward high-photon-energy XFELs.
文摘Correction to:Nuclear Science and Techniques(2025)36:100 https://doi.org/10.1007/s41365-025-01692-6 In this article,Fig.9 appeared incorrectly and have now been corrected in the original publication.For completeness and transparency,both correct and incorrect versions are displayed below.
基金supported by the National Grand Instrument Project(No.2019YFF01014404)the Natural Science Foundation of China(Nos.12122501,61631001,11921006,U2241281,and 11975037)the Foundation of Science and Technology on Plasma Physics Laboratory(No.6142A04220108)。
文摘High-power laser pulses interacting with targets can generate intense electromagnetic pulses(EMPs),which can disrupt physical experimental diagnostics and even damage diagnostic equipment,posing a threat to the reliable operation of experiments.In this study,EMPs resulting from multi-petawatt laser irradiating nitrogen gas jets were systematically analyzed and investigated.The experimental results revealed that the EMP amplitude is positively correlated with the quantity and energy of the electrons captured and accelerated by the plasma channel.These factors are reflected by parameters such as laser energy and nitrogen gas jet pressure.Additionally,we propose several potential sources of EMPs produced by laser-irradiated gas jets and separately analyzed their spatiotemporal distributions.The findings provide insight into the mechanisms of EMP generation and introduce a new approach to achieve controllable EMPs by regulating the laser energy and gas jet pressure.
