INTRODUCTION The articles in the“Atomic and molecular physics for controlled fusion and astrophysics”special issue cover a wide range of topics in atomic and molecular physics in the context of hot plasmas.Basic ato...INTRODUCTION The articles in the“Atomic and molecular physics for controlled fusion and astrophysics”special issue cover a wide range of topics in atomic and molecular physics in the context of hot plasmas.Basic atomic processes are of fundamental importance in confinement fusion and astrophysical environments,and also for ultrahigh–intensity interaction of lasers with matter.Atomic physics in extreme environments such as high pressures and hot or dense plasmas^(1,2)presents new challenges to the community,and these have to be addressed by both theoretical and experimental studies.Several extreme configurations are investigated in this special issue,which should be understood as an initiative to draw the attention of the community to important ongoing work in the context of extreme states of matter.This special issue presents eight articles from scientists actively working in this field and shows the important advances that have been made in basic atomic processes and related areas of plasma properties and plasma diagnosis over the last few years.展开更多
The concept of local shock strength and a quantitative measure index str of local shock strength are proposed,derived from the oblique shock relation and the monotonic relationship between total pressure loss ratio an...The concept of local shock strength and a quantitative measure index str of local shock strength are proposed,derived from the oblique shock relation and the monotonic relationship between total pressure loss ratio and normal Mach number.Utilizing the high density gradient characteristic of shock waves and the oblique shock relation,a post-processing algorithm for two-dimensional flow field data is developed.The objective of the post-processing algorithm is to obtain specific shock wave location coordinates and calculate the corresponding str from flow filed data under the calibration of the oblique shock relation.Valida-tion of this post-processing algorithm is conducted using a standard model example that can be solved analytically.Combining the concept of local shock strength with the post-processing algorithm,a local shock strength quantitative mapping approach is established for the first time.This approach enables a quantitative measure and visualization of local shock strength at distinct locations,represented by color mapping on the shock structures.The approach can be applied to post-processing numerical sim-ulation data of two-dimensional flows.Applications to the intersection of two left-running oblique shock waves(straight shock waves),the bow shock in front of a cylinder(curved shock wave),and Mach reflection(mixed straight and curved shock waves)demonstrate the accuracy,and effectiveness of the mapping approach in investigating diverse shock wave phenomena.The quan-titative mapping approach of str may be a valuable tool in the design of supersonic/hypersonic vehicles and the exploration of shock wave evolution.展开更多
This study uses nonequilibrium molecular dynamics simulations to explore the dynamic failures and deformation mechanisms of a cylindrical shell composed of nanocrystalline nickel-titanium alloy under implosion loading...This study uses nonequilibrium molecular dynamics simulations to explore the dynamic failures and deformation mechanisms of a cylindrical shell composed of nanocrystalline nickel-titanium alloy under implosion loading.We discover that some individual spall planes are sequentially generated in the material along the propagation of a radial stress wave,indicative of the formation of multiple spallation.For larger grain sizes,void nucleation at the first spallation occurs in a coexisting intergranular/transgranular manner,whereas with decreasing grain size,voids tend to nucleate along the grain boundaries.Correspondingly,the spall strength exhibits a transition from an inverse Hall-Petch to a Hall-Petch relationship.For larger grain sizes,at the secondary spallation,localized shearing zones and grain boundaries provide potential void-nucleated sites.Importantly,the formation of shear deformation bands promotes grain refinement,contributing to a reduction in the dislocation-induced strengthening effect.Consequently,a lower spall strength is produced,in contrast to the first spallation.As the grain size becomes smaller,voids nucleate mostly along grain boundaries,and plastic deformation is dominated by dense grain boundaries.Overall,the high temperature caused by shear localization leads to material weakening,and in turn there is a significant decrease in the spall strength for the secondary spallation,compared with the first.Finally,significant penetration between two spall planes is observed for large grain size,which can be attributed to the nucleation of voids on linking grain boundaries,with temperatures exceeding the melting point of the material.展开更多
Abstract The continuously rotating detonation engine （CRDE） is a new concept of engines for air- craft and spacecraft. Quasi-stable continuously rotating detonation （CRD） can be observed in an annular combustion c...Abstract The continuously rotating detonation engine （CRDE） is a new concept of engines for air- craft and spacecraft. Quasi-stable continuously rotating detonation （CRD） can be observed in an annular combustion chamber, but the sustaining, stabilizing and adjusting mechanisms are not yet clear. To learn more deeply into the CRDE, experimental studies have been carried out to inves- tigate hydrogen-oxygen CRDE. Pressure histories are obtained during each shot, which show that stable CRD waves are generated in the combustor, when feeding pressures are higher than 0.5 MPa for fuel and oxidizer, respectively. Each shot can keep running as long as fresh gas feeding main- tains. Close-up of the pressure history shows the repeatability of pressure peaks and indicates the detonation velocity in hydrogen-oxygen CRD, which proves the success of forming a stable CRD in the annular chamber. Spectrum of the pressure history matches the close-up analysis and confirms the CRD. It also shows multi-wave phenomenon and affirms the fact that in this case a single detonation wave is rotating in the annulus. Moreover, oscillation phenomenon is found in pressure peaks and a self-adjusting mechanism is proposed to explain the phenomenon.展开更多
A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving speci...A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.展开更多
The physics of laser-plasma interaction is studied on the Shenguang III prototype laser facility under conditions relevant to inertial confinement fusion designs.A sub-millimeter-size underdense hot plasma is created ...The physics of laser-plasma interaction is studied on the Shenguang III prototype laser facility under conditions relevant to inertial confinement fusion designs.A sub-millimeter-size underdense hot plasma is created by ionization of a low-density plastic foam by four high-energy(3.2 kJ)laser beams.An interaction beam is fired with a delay permitting evaluation of the excitation of parametric instabilities at different stages of plasma evolution.Multiple diagnostics are used for plasma characterization,scattered radiation,and accelerated electrons.The experimental results are analyzed with radiation hydrodynamic simulations that take account of foam ionization and homogenization.The measured level of stimulated Raman scattering is almost one order of magnitude larger than that measured in experiments with gasbags and hohlraums on the same installation,possibly because of a greater plasma density.Notable amplification is achieved in high-intensity speckles,indicating the importance of implementing laser temporal smoothing techniques with a large bandwidth for controlling laser propagation and absorption.展开更多
The droplet formation dynamics of a Newtonian liquid in a drop-on-demand (DOD) inkjet process is numerically investigated by using a volume-of-fluid (VOF) method. We focus on the nozzle geometry, wettability of the in...The droplet formation dynamics of a Newtonian liquid in a drop-on-demand (DOD) inkjet process is numerically investigated by using a volume-of-fluid (VOF) method. We focus on the nozzle geometry, wettability of the interior surface, and the fluid properties to achieve the stable droplet formation with higher velocity. It is found that a nozzle with contracting angle of 45° generates the most stable and fastest single droplet, which is beneficial for the enhanced printing quality and high-throughput printing rate. For this nozzle with the optimal geometry, we systematically change the wettability of the interior surface, i.e., different contact angles. As the contact angle increases, pinch-off time increases and the droplet speed reduces. Finally, fluids with different properties are investigated to identify the printability range.展开更多
The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical sym- metry properties of the numerical schemes are discus...The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical sym- metry properties of the numerical schemes are discussed, and an area weighted scheme is extended from a Lagrangian method to an arbitrary Lagrangian and Eulerian （ALE） method. Numerical results are presented to compare three discrete configurations, i.e., the control volume scheme, the area weighted scheme, and the plane scheme with the addition of a geometrical source. The fact that the singularity arises from the geometri- cal source term in the plane scheme is illustrated. A suggestion for choosing the discrete formulation is given when the strong shock wave problems are simulated.展开更多
We propose an efficient and robust way to design absorbing boundary conditions in atomistic computations. An optimal discrete boundary condition is obtained by minimizing a functional of a reflection coefficient integ...We propose an efficient and robust way to design absorbing boundary conditions in atomistic computations. An optimal discrete boundary condition is obtained by minimizing a functional of a reflection coefficient integral over a range of wave numbers. The minimization is performed with respect to a set of wave numbers, at which transparent absorption is reached. Compared with the optimization with respect to the boundary condition coefficients suggested by E and Huang [Phys. Rev. Lett. 87 （2001） 133501], we reduce considerably the number of independent variables and the computing cost. We further demonstrate with numerical examples that both the optimization and the wave absorption are more robust in the proposed design.展开更多
In the present paper, we investigate the instability, adiabaticity, and controlling effects of external fields for a dark state in a homonuclear atom-tetramer conversion that is implemented by a generalized stimulated...In the present paper, we investigate the instability, adiabaticity, and controlling effects of external fields for a dark state in a homonuclear atom-tetramer conversion that is implemented by a generalized stimulated Raman adiabatic passage. We analytically obtain the regions for the appearance of dynamical instability and study the adiabatic evolution by a newly defined adiabatic fidelity. Moreover, the effects of the external field parameters and the spontaneous emissions on the conversion efficiency are also investigated.展开更多
A simplified theoretical model for the linear Rayleigh-Taylor instability of finite thickness elastic-plastic solid constantly accelerated by finite thickness viscous fluid is performed.With the irrotational assumptio...A simplified theoretical model for the linear Rayleigh-Taylor instability of finite thickness elastic-plastic solid constantly accelerated by finite thickness viscous fluid is performed.With the irrotational assumption,it is possible to consider viscosity,surface tension,elasticity or plasticity effects simultaneously.The model considers thicknesses at rigid wall boundary conditions with the velocity potentials,and deals with solid elastic-plastic transition and fluid viscosity based on the velocity continuity and force equilibrium at contact interface.The complete analytical expressions of the amplitude motion equation,the growth rate,and the instability boundary are obtained for arbitrary Atwood number,viscosity,thicknesses of solid and fluid.The thicknesses effects of two materials on the growth rate and the instability boundary are discussed.展开更多
The photon Doppler velocimetry(PDV) spectrum is investigated in an attempt to reveal the particle parameters of ejecta from shock-loaded samples in a vacuum. A GPU-accelerated Monte–Carlo algorithm, which considers t...The photon Doppler velocimetry(PDV) spectrum is investigated in an attempt to reveal the particle parameters of ejecta from shock-loaded samples in a vacuum. A GPU-accelerated Monte–Carlo algorithm, which considers the multiplescattering effects of light, is applied to reconstruct the light field of the ejecta and simulate the corresponding PDV spectrum.The influence of the velocity profile, total area mass, and particle size of the ejecta on the simulated spectra is discussed qualitatively. To facilitate a quantitative discussion, a novel theoretical optical model is proposed in which the singlescattering assumption is applied. With this model, the relationships between the particle parameters of ejecta and the peak information of the PDV spectrum are derived, enabling direct extraction of the particle parameters from the PDV spectrum.The values of the ejecta parameters estimated from the experimental spectrum are in good agreement with those measured by a piezoelectric probe.展开更多
The interaction between shock waves and multiple cylinders,referred to as shock–cylinder interaction(SCI),is an important phenomenon in science and engineering.However,its underlying physical mechanisms remain unclea...The interaction between shock waves and multiple cylinders,referred to as shock–cylinder interaction(SCI),is an important phenomenon in science and engineering.However,its underlying physical mechanisms remain unclear.This study entailed the numerical simulation of the aerobreakup of two tandem water columns subjected to a high-speed gas flow by using an adaptive mesh refinement(AMR)-based diffusion-interface model.The objective was to elucidate the changes in water–column deformation patterns over a wide range of Weber numbers.Statistical analysis was performed to examine the deformation of the water columns in vertical directions.Results reveal distinct deformation patterns between the two columns as the Weber number increases.Additionally,an extended exponential stretching law model was devised,and its improved capability to predict the deformation patterns was demonstrated.展开更多
We are intrigued by the issues of shock instability,with a particular emphasis on numerical schemes that address the carbuncle phenomenon by reducing dissipation rather than increasing it.For a specific class of plana...We are intrigued by the issues of shock instability,with a particular emphasis on numerical schemes that address the carbuncle phenomenon by reducing dissipation rather than increasing it.For a specific class of planar flow fields where the transverse direction exhibits vanishing but non-zero velocity components,such as a disturbed onedimensional(1D)steady shock wave,we conduct a formal asymptotic analysis for the Euler system and associated numerical methods.This analysis aims to illustrate the discrepancies among various low-dissipative numerical algorithms.Furthermore,a numerical stability analysis of steady shock is undertaken to identify the key factors underlying shock-stable algorithms.To verify the stability mechanism,a consistent,low-dissipation,and shock-stable HLLC-type Riemann solver is presented.展开更多
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.展开更多
This paper extends the previous work[1]for the three-temperature gray radiative transfer equations to the frequency-dependent case.Since the additional frequency variable is considered,the equations are more complicat...This paper extends the previous work[1]for the three-temperature gray radiative transfer equations to the frequency-dependent case.Since the additional frequency variable is considered,the equations are more complicated than those in the gray case.Moreover,opacity may be typically a decreasing function of the frequency variable in applications.At the same spatial location,the equations can be in the optically thick case for low frequency photons,while in the optically thin case for high frequency ones.Thus,the resulting discrete equations can significantly increase the computational cost for opacity having the multi-scale property in multiple frequency radiation.Due to the presence of the radiation-electron coupling,electronion coupling,and electron-ion diffusion terms,the model under consideration exhibits strong nonlinearity and strong coupling properties.In this paper,the multigroup method is used to discretize the frequency variable and the H_(N)^(T)method to discretize the angular variable first.Then,within the framework of a unified gas kinetic scheme(UGKS),a multigroup H_(N)^(T)-UGKS method is constructed to solve this complex model iteratively.Furthermore,it can be shown that as the Knudsen number tends to zero,with variations in the electron-ion coupling,absorption,and scattering coefficients,the multigroup H_(N)^(T)-UGKS scheme can converge to numerical schemes for the single-temperature,two-temperature,and the frequency-dependent three-temperature,two-temperature diffusion limit equations,respectively.Finally,several numerical examples are provided to validate the effectiveness and stability of the proposed scheme.展开更多
For the observed line at 799.23°A in tungsten EBIT experiment,which was assigned to be^(3)F_(4)^(o)−^(3)F_(3)^(o)([Ar]4s^(2)4p^(5)4d)of W^(38+)ion,there were noticeable deviations for most calculated wavelengths ...For the observed line at 799.23°A in tungsten EBIT experiment,which was assigned to be^(3)F_(4)^(o)−^(3)F_(3)^(o)([Ar]4s^(2)4p^(5)4d)of W^(38+)ion,there were noticeable deviations for most calculated wavelengths from the measured value.To clarify this issue,we carry out an extensive calculation for energy levels and transition properties of W^(38+)ion using the multi-configuration Dirac–Hartree–Fock and relativistic configuration interaction method,in which more deeper inner core electron correlations are included,and different forms of Breit interaction as well as quantum electrodynamics corrections are investigated.It is found that the inner core electron correlations can affect the total energy of levels,while only slightly modify the excited energy of levels in 4s^(2)4p^(5)4d complex.The present calculated wavelengths agree with the corresponding measured values excellently except the line at 799.23Å.Thus we are strongly suspicious this line should be misidentified,and suggest that new experiment with higher resolution and spectra analysis based on more accurate atomic data should be performed for W^(38+)ion.展开更多
FeO_(2)is proposed to be a kind of substance in the Earth’s lower mantle in recent years.In this pa-per,the equation of state,elastic properties and sound velocities are obtained based on the first principle calcula-...FeO_(2)is proposed to be a kind of substance in the Earth’s lower mantle in recent years.In this pa-per,the equation of state,elastic properties and sound velocities are obtained based on the first principle calcula-tions.By solving the Boltzmann transport equations,we investigated the lattice thermal con-ductivity of FeO_(2)under high pressure and high temperature.The calculated compressional and shear sound velocities of FeO_(2)agree with the data of preliminary reference Earth model.The lattice thermal conductivity of FeO_(2)at core-mantleboundary(~135 GPa,~3500 K)is 1.77 W/mK,and the total thermal conductivity is 135.10 W/mK.The influence of lattice thermal conductivity can be ignored above 3000 K.展开更多
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.展开更多
As an emerging multifunctional metal with the lowest melting point except for mercury,gallium combines a wide range of metallic and non-metallic elements to form advanced semiconductors critically important in cutting...As an emerging multifunctional metal with the lowest melting point except for mercury,gallium combines a wide range of metallic and non-metallic elements to form advanced semiconductors critically important in cutting-edge technologies.However,due to its low melting point and poor machinability,it is quite difficult to simultaneously characterize gallium’s elastic properties and damping characteristics using conventional methods,which is es-sential in designing and evaluating gallium-based structures.Therefore,developing effective methods to achieve accurate and efficient measurements of Young’s modulus and corresponding internal friction of gallium is of great significance.This letter studies simultaneous measurements of the variations in Young’s modulus and internal friction of gallium at varying temperatures by employing the modified piezoelectric ultrasonic composite oscil-lator technique.Combining the explicit theoretical formulas with the measured resonance and anti-resonance frequencies,it has been discovered that Young’s modulus undergoes an approximately linear decrease as the temperature rises,declining from 83.84 GPa at -70℃ to 79.37 GPa at 20℃.Moreover,like aluminum in the same Group ⅢA of the Periodic Table of Elements and exhibits a grain-boundary internal friction peak,gallium displays a longitudinal internal friction peak at approximately-12°C,with the peak value reaching 1.77×10^(-3).This basic research on gallium’s elastic properties and damping characteristics under low-temperature condi-tions will inspire further explorations of the mechanical properties of a diverse spectrum of low-melting-point functional materials and facilitate applications of gallium-based structures under complex conditions.展开更多
文摘INTRODUCTION The articles in the“Atomic and molecular physics for controlled fusion and astrophysics”special issue cover a wide range of topics in atomic and molecular physics in the context of hot plasmas.Basic atomic processes are of fundamental importance in confinement fusion and astrophysical environments,and also for ultrahigh–intensity interaction of lasers with matter.Atomic physics in extreme environments such as high pressures and hot or dense plasmas^(1,2)presents new challenges to the community,and these have to be addressed by both theoretical and experimental studies.Several extreme configurations are investigated in this special issue,which should be understood as an initiative to draw the attention of the community to important ongoing work in the context of extreme states of matter.This special issue presents eight articles from scientists actively working in this field and shows the important advances that have been made in basic atomic processes and related areas of plasma properties and plasma diagnosis over the last few years.
基金supported by the National Natural Science Foundation of China(Grant No.12372233)the Fund of NPU-Duke China Seed Program(Grant No.119003067)the“111 Project”of China(Grant No.B17037-106).
文摘The concept of local shock strength and a quantitative measure index str of local shock strength are proposed,derived from the oblique shock relation and the monotonic relationship between total pressure loss ratio and normal Mach number.Utilizing the high density gradient characteristic of shock waves and the oblique shock relation,a post-processing algorithm for two-dimensional flow field data is developed.The objective of the post-processing algorithm is to obtain specific shock wave location coordinates and calculate the corresponding str from flow filed data under the calibration of the oblique shock relation.Valida-tion of this post-processing algorithm is conducted using a standard model example that can be solved analytically.Combining the concept of local shock strength with the post-processing algorithm,a local shock strength quantitative mapping approach is established for the first time.This approach enables a quantitative measure and visualization of local shock strength at distinct locations,represented by color mapping on the shock structures.The approach can be applied to post-processing numerical sim-ulation data of two-dimensional flows.Applications to the intersection of two left-running oblique shock waves(straight shock waves),the bow shock in front of a cylinder(curved shock wave),and Mach reflection(mixed straight and curved shock waves)demonstrate the accuracy,and effectiveness of the mapping approach in investigating diverse shock wave phenomena.The quan-titative mapping approach of str may be a valuable tool in the design of supersonic/hypersonic vehicles and the exploration of shock wave evolution.
基金support of the National Natural Science Foundation of China under Grant Nos.12372367 and 12202081the Special Foundation from the Institute of Fluid Physics of CAEP under Grant No.2022-YCHT-0641.
文摘This study uses nonequilibrium molecular dynamics simulations to explore the dynamic failures and deformation mechanisms of a cylindrical shell composed of nanocrystalline nickel-titanium alloy under implosion loading.We discover that some individual spall planes are sequentially generated in the material along the propagation of a radial stress wave,indicative of the formation of multiple spallation.For larger grain sizes,void nucleation at the first spallation occurs in a coexisting intergranular/transgranular manner,whereas with decreasing grain size,voids tend to nucleate along the grain boundaries.Correspondingly,the spall strength exhibits a transition from an inverse Hall-Petch to a Hall-Petch relationship.For larger grain sizes,at the secondary spallation,localized shearing zones and grain boundaries provide potential void-nucleated sites.Importantly,the formation of shear deformation bands promotes grain refinement,contributing to a reduction in the dislocation-induced strengthening effect.Consequently,a lower spall strength is produced,in contrast to the first spallation.As the grain size becomes smaller,voids nucleate mostly along grain boundaries,and plastic deformation is dominated by dense grain boundaries.Overall,the high temperature caused by shear localization leads to material weakening,and in turn there is a significant decrease in the spall strength for the secondary spallation,compared with the first.Finally,significant penetration between two spall planes is observed for large grain size,which can be attributed to the nucleation of voids on linking grain boundaries,with temperatures exceeding the melting point of the material.
基金supported by the National Natural Science Foundation of China(No.91441110)
文摘Abstract The continuously rotating detonation engine （CRDE） is a new concept of engines for air- craft and spacecraft. Quasi-stable continuously rotating detonation （CRD） can be observed in an annular combustion chamber, but the sustaining, stabilizing and adjusting mechanisms are not yet clear. To learn more deeply into the CRDE, experimental studies have been carried out to inves- tigate hydrogen-oxygen CRDE. Pressure histories are obtained during each shot, which show that stable CRD waves are generated in the combustor, when feeding pressures are higher than 0.5 MPa for fuel and oxidizer, respectively. Each shot can keep running as long as fresh gas feeding main- tains. Close-up of the pressure history shows the repeatability of pressure peaks and indicates the detonation velocity in hydrogen-oxygen CRD, which proves the success of forming a stable CRD in the annular chamber. Spectrum of the pressure history matches the close-up analysis and confirms the CRD. It also shows multi-wave phenomenon and affirms the fact that in this case a single detonation wave is rotating in the annulus. Moreover, oscillation phenomenon is found in pressure peaks and a self-adjusting mechanism is proposed to explain the phenomenon.
基金Supported by National Natural Science Foundation of China under Grant Nos.11772064,and 11502117CAEP Foundation under Grant No.CX2019033+1 种基金the opening project of State Key Laboratory of Explosion Science and Technology(Beijing Institute of Technology)Science Challenge Project under Grant No.JCKY2016212A501
文摘A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.
基金This project was partially supported by the Advanced Research Using High Intensity Laser Produced Photons and Particles(ADONIS)project(Grant No.CZ.02.1.01/0.0/0.0/16_019/0000789)the CAAS project(Grant No.CZ.02.1.01/0.0/0.0/16_019/0000778)+3 种基金both from the European Regional Development FundThe results of the LQ1606 project were partially obtained with the financial support from the Ministry of Education,Youth and Sports as part of targeted support from the National Programme of Sustainability IIThe authors acknowledge support from the National Natural Science Foundation of China(Grant Nos.11775033,11875241,11975215,11905204,12035002)the Laser Fusion Research Center Funds for Young Talents(Grant No.RCFPD3-2019-6).
文摘The physics of laser-plasma interaction is studied on the Shenguang III prototype laser facility under conditions relevant to inertial confinement fusion designs.A sub-millimeter-size underdense hot plasma is created by ionization of a low-density plastic foam by four high-energy(3.2 kJ)laser beams.An interaction beam is fired with a delay permitting evaluation of the excitation of parametric instabilities at different stages of plasma evolution.Multiple diagnostics are used for plasma characterization,scattered radiation,and accelerated electrons.The experimental results are analyzed with radiation hydrodynamic simulations that take account of foam ionization and homogenization.The measured level of stimulated Raman scattering is almost one order of magnitude larger than that measured in experiments with gasbags and hohlraums on the same installation,possibly because of a greater plasma density.Notable amplification is achieved in high-intensity speckles,indicating the importance of implementing laser temporal smoothing techniques with a large bandwidth for controlling laser propagation and absorption.
基金Project supported by the National Natural Science Foundation of China(No.11802004)
文摘The droplet formation dynamics of a Newtonian liquid in a drop-on-demand (DOD) inkjet process is numerically investigated by using a volume-of-fluid (VOF) method. We focus on the nozzle geometry, wettability of the interior surface, and the fluid properties to achieve the stable droplet formation with higher velocity. It is found that a nozzle with contracting angle of 45° generates the most stable and fastest single droplet, which is beneficial for the enhanced printing quality and high-throughput printing rate. For this nozzle with the optimal geometry, we systematically change the wettability of the interior surface, i.e., different contact angles. As the contact angle increases, pinch-off time increases and the droplet speed reduces. Finally, fluids with different properties are investigated to identify the printability range.
基金Project supported by the National Natural Science Foundation of China(Nos.11471048 and U1630249)the Foundation of Chinese Academy of Engineering Physics(No.2014A0202010)the Foundation of Laboratory of Computational Physics(No.9140C690202140C69293)
文摘The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical sym- metry properties of the numerical schemes are discussed, and an area weighted scheme is extended from a Lagrangian method to an arbitrary Lagrangian and Eulerian （ALE） method. Numerical results are presented to compare three discrete configurations, i.e., the control volume scheme, the area weighted scheme, and the plane scheme with the addition of a geometrical source. The fact that the singularity arises from the geometri- cal source term in the plane scheme is illustrated. A suggestion for choosing the discrete formulation is given when the strong shock wave problems are simulated.
基金Supported in part by the National Natural Science Foundation of China under Grant No 10872004, the National Basic Research Program of China under Grant No 2007CB814800, and the Ministry of Education of China under Grant Nos NCET-06-0011 and 200800010013.
文摘We propose an efficient and robust way to design absorbing boundary conditions in atomistic computations. An optimal discrete boundary condition is obtained by minimizing a functional of a reflection coefficient integral over a range of wave numbers. The minimization is performed with respect to a set of wave numbers, at which transparent absorption is reached. Compared with the optimization with respect to the boundary condition coefficients suggested by E and Huang [Phys. Rev. Lett. 87 （2001） 133501], we reduce considerably the number of independent variables and the computing cost. We further demonstrate with numerical examples that both the optimization and the wave absorption are more robust in the proposed design.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11005055,11075020,and 11204117)the National Fundamental Research Programme of China(Grant No.2011CB921503)+1 种基金the Ph.D.Programs Foundation of Liaoning Provincial Science and Technology Bureau(GrantNo.201103778)the Higher School Excellent Researcher Award Program from the Educational Department of Liaoning Province of China(GrantNo.LJQ2011005)
文摘In the present paper, we investigate the instability, adiabaticity, and controlling effects of external fields for a dark state in a homonuclear atom-tetramer conversion that is implemented by a generalized stimulated Raman adiabatic passage. We analytically obtain the regions for the appearance of dynamical instability and study the adiabatic evolution by a newly defined adiabatic fidelity. Moreover, the effects of the external field parameters and the spontaneous emissions on the conversion efficiency are also investigated.
基金Project supported by of the Science Challenge Project of China(Grant No.TZ2018001)。
文摘A simplified theoretical model for the linear Rayleigh-Taylor instability of finite thickness elastic-plastic solid constantly accelerated by finite thickness viscous fluid is performed.With the irrotational assumption,it is possible to consider viscosity,surface tension,elasticity or plasticity effects simultaneously.The model considers thicknesses at rigid wall boundary conditions with the velocity potentials,and deals with solid elastic-plastic transition and fluid viscosity based on the velocity continuity and force equilibrium at contact interface.The complete analytical expressions of the amplitude motion equation,the growth rate,and the instability boundary are obtained for arbitrary Atwood number,viscosity,thicknesses of solid and fluid.The thicknesses effects of two materials on the growth rate and the instability boundary are discussed.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11902043 and 11772065)the Science Challenge Project (Grant No. TZ2016001)。
文摘The photon Doppler velocimetry(PDV) spectrum is investigated in an attempt to reveal the particle parameters of ejecta from shock-loaded samples in a vacuum. A GPU-accelerated Monte–Carlo algorithm, which considers the multiplescattering effects of light, is applied to reconstruct the light field of the ejecta and simulate the corresponding PDV spectrum.The influence of the velocity profile, total area mass, and particle size of the ejecta on the simulated spectra is discussed qualitatively. To facilitate a quantitative discussion, a novel theoretical optical model is proposed in which the singlescattering assumption is applied. With this model, the relationships between the particle parameters of ejecta and the peak information of the PDV spectrum are derived, enabling direct extraction of the particle parameters from the PDV spectrum.The values of the ejecta parameters estimated from the experimental spectrum are in good agreement with those measured by a piezoelectric probe.
基金supported by the National Natural Science Foundation of China (Grant Nos.12202070 and 11772065)the Foundation of National Key Laboratory of Computational Physics.
文摘The interaction between shock waves and multiple cylinders,referred to as shock–cylinder interaction(SCI),is an important phenomenon in science and engineering.However,its underlying physical mechanisms remain unclear.This study entailed the numerical simulation of the aerobreakup of two tandem water columns subjected to a high-speed gas flow by using an adaptive mesh refinement(AMR)-based diffusion-interface model.The objective was to elucidate the changes in water–column deformation patterns over a wide range of Weber numbers.Statistical analysis was performed to examine the deformation of the water columns in vertical directions.Results reveal distinct deformation patterns between the two columns as the Weber number increases.Additionally,an extended exponential stretching law model was devised,and its improved capability to predict the deformation patterns was demonstrated.
基金Project supported by the National Natural Science Foundation of China(Nos.12471367 and12361076)the Research Program of Science and Technology at Universities of Inner Mongolia Autonomous Region(Nos.NJZY19186,NJZY22036,and NJZY23003)。
文摘We are intrigued by the issues of shock instability,with a particular emphasis on numerical schemes that address the carbuncle phenomenon by reducing dissipation rather than increasing it.For a specific class of planar flow fields where the transverse direction exhibits vanishing but non-zero velocity components,such as a disturbed onedimensional(1D)steady shock wave,we conduct a formal asymptotic analysis for the Euler system and associated numerical methods.This analysis aims to illustrate the discrepancies among various low-dissipative numerical algorithms.Furthermore,a numerical stability analysis of steady shock is undertaken to identify the key factors underlying shock-stable algorithms.To verify the stability mechanism,a consistent,low-dissipation,and shock-stable HLLC-type Riemann solver is presented.
基金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.
基金supported by the Beijing Natural Science Foundation(Z230003)for Sunby the National Key R&D Program(2020YFA0712200)+1 种基金the National Key Project(GJXM92579)the Sino-German Science Center(GZ 1465)for Jiang。
文摘This paper extends the previous work[1]for the three-temperature gray radiative transfer equations to the frequency-dependent case.Since the additional frequency variable is considered,the equations are more complicated than those in the gray case.Moreover,opacity may be typically a decreasing function of the frequency variable in applications.At the same spatial location,the equations can be in the optically thick case for low frequency photons,while in the optically thin case for high frequency ones.Thus,the resulting discrete equations can significantly increase the computational cost for opacity having the multi-scale property in multiple frequency radiation.Due to the presence of the radiation-electron coupling,electronion coupling,and electron-ion diffusion terms,the model under consideration exhibits strong nonlinearity and strong coupling properties.In this paper,the multigroup method is used to discretize the frequency variable and the H_(N)^(T)method to discretize the angular variable first.Then,within the framework of a unified gas kinetic scheme(UGKS),a multigroup H_(N)^(T)-UGKS method is constructed to solve this complex model iteratively.Furthermore,it can be shown that as the Knudsen number tends to zero,with variations in the electron-ion coupling,absorption,and scattering coefficients,the multigroup H_(N)^(T)-UGKS scheme can converge to numerical schemes for the single-temperature,two-temperature,and the frequency-dependent three-temperature,two-temperature diffusion limit equations,respectively.Finally,several numerical examples are provided to validate the effectiveness and stability of the proposed scheme.
基金supported by the Science Challenge Project of China Academy of Engineering Physics(CAEP)(Grant No.TZ2018005)the National Natural Science Foundation of China(Grant Nos.12474277,12374259,12104095,12074081,and 12074082).
文摘For the observed line at 799.23°A in tungsten EBIT experiment,which was assigned to be^(3)F_(4)^(o)−^(3)F_(3)^(o)([Ar]4s^(2)4p^(5)4d)of W^(38+)ion,there were noticeable deviations for most calculated wavelengths from the measured value.To clarify this issue,we carry out an extensive calculation for energy levels and transition properties of W^(38+)ion using the multi-configuration Dirac–Hartree–Fock and relativistic configuration interaction method,in which more deeper inner core electron correlations are included,and different forms of Breit interaction as well as quantum electrodynamics corrections are investigated.It is found that the inner core electron correlations can affect the total energy of levels,while only slightly modify the excited energy of levels in 4s^(2)4p^(5)4d complex.The present calculated wavelengths agree with the corresponding measured values excellently except the line at 799.23Å.Thus we are strongly suspicious this line should be misidentified,and suggest that new experiment with higher resolution and spectra analysis based on more accurate atomic data should be performed for W^(38+)ion.
基金supported by the National Natural Science Foundation of China(No.12072044)the Natu-ral Science Foundation of Chongqing City(No.cstc2020jcyj-msxmX0616).
文摘FeO_(2)is proposed to be a kind of substance in the Earth’s lower mantle in recent years.In this pa-per,the equation of state,elastic properties and sound velocities are obtained based on the first principle calcula-tions.By solving the Boltzmann transport equations,we investigated the lattice thermal con-ductivity of FeO_(2)under high pressure and high temperature.The calculated compressional and shear sound velocities of FeO_(2)agree with the data of preliminary reference Earth model.The lattice thermal conductivity of FeO_(2)at core-mantleboundary(~135 GPa,~3500 K)is 1.77 W/mK,and the total thermal conductivity is 135.10 W/mK.The influence of lattice thermal conductivity can be ignored above 3000 K.
基金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.
基金financially supported by the National Key R&D Program of China(Grant No.2023YFF0716800)the National Natural Science Foundation of China(Grant No.12074160)the Natural Science Foundation of Liaoning Province of China(Grant No.2024-MS-181).
文摘As an emerging multifunctional metal with the lowest melting point except for mercury,gallium combines a wide range of metallic and non-metallic elements to form advanced semiconductors critically important in cutting-edge technologies.However,due to its low melting point and poor machinability,it is quite difficult to simultaneously characterize gallium’s elastic properties and damping characteristics using conventional methods,which is es-sential in designing and evaluating gallium-based structures.Therefore,developing effective methods to achieve accurate and efficient measurements of Young’s modulus and corresponding internal friction of gallium is of great significance.This letter studies simultaneous measurements of the variations in Young’s modulus and internal friction of gallium at varying temperatures by employing the modified piezoelectric ultrasonic composite oscil-lator technique.Combining the explicit theoretical formulas with the measured resonance and anti-resonance frequencies,it has been discovered that Young’s modulus undergoes an approximately linear decrease as the temperature rises,declining from 83.84 GPa at -70℃ to 79.37 GPa at 20℃.Moreover,like aluminum in the same Group ⅢA of the Periodic Table of Elements and exhibits a grain-boundary internal friction peak,gallium displays a longitudinal internal friction peak at approximately-12°C,with the peak value reaching 1.77×10^(-3).This basic research on gallium’s elastic properties and damping characteristics under low-temperature condi-tions will inspire further explorations of the mechanical properties of a diverse spectrum of low-melting-point functional materials and facilitate applications of gallium-based structures under complex conditions.
