In the framework of elastoplastic theory,by introducing dissipative plastic energy(instead of cumulative plastic strain)and dissipative plastic energy rate(instead of cumulative plastic strain rate)into the ratchettin...In the framework of elastoplastic theory,by introducing dissipative plastic energy(instead of cumulative plastic strain)and dissipative plastic energy rate(instead of cumulative plastic strain rate)into the ratchetting parameter evolution equation and isotropic evolution rules respectively,a cyclic elastoplastic constitutive model based on dissipative plastic energy is established.This model,termed the WDP model,describes the physical meaning and evolution rule of the unclosed stress–strain hysteresis loop using an energy method.A comparison of numerical implementation results with experimental data demonstrates the capability of the WDP model to predict the cyclic deformation of EA4T steel,effectively capturing the cyclic softening characteristics and ratchetting behaviors of axle steel EA4T.展开更多
The electrospray thruster supplied by ionic liquid is a promising micro-propulsion thruster with small size and precise thrust, which can emit both cations and anions to achieve self-neutralization. In order to furthe...The electrospray thruster supplied by ionic liquid is a promising micro-propulsion thruster with small size and precise thrust, which can emit both cations and anions to achieve self-neutralization. In order to further investigate the effect of ion solvation energy on the evaporation of cations and anions from ionic liquid under the action of a uniform electric field, this paper establishes a transient Electrohydrodynamic (EHD) model for free ionic liquid droplets undergoing ion evaporation. The dynamic processes of droplet deformation and ion evaporation are simulated. And the study further focuses on the influence of different ion solvation energies for cations on the droplet morphology and the ion evaporation characteristics at the positively charged end and negatively charged end of the droplet. The results indicate that, when the ion solvation energy for cations is higher than that of anions, it will cause the ion evaporation at the positively charged end of the droplet to lag behind the ion evaporation at the negatively charged end. And the higher the ion solvation energy for the cations, the longer the evaporation lag time at the positively charged end of the droplet, which will lead to a higher peak of surface charge density that can be reached, resulting in a larger evaporation current and sharper droplet stretching deformation. Additionally, the peak surface charge density of the positively charged end of the droplet is linearly related to the ion solvation energy for cations, while the peak surface charge density of the negatively charged end remains almost unchanged and is not significantly affected by the ion solvation energy for cations.展开更多
Coastal wetlands are crucial for the‘blue carbon sink’,significantly contributing to regulating climate change.This study util-ized 160 soil samples,35 remote sensing features,and 5 geo-climatic data to accurately e...Coastal wetlands are crucial for the‘blue carbon sink’,significantly contributing to regulating climate change.This study util-ized 160 soil samples,35 remote sensing features,and 5 geo-climatic data to accurately estimate the soil organic carbon stocks(SOCS)in the coastal wetlands of Tianjin and Hebei,China.To reduce data redundancy,simplify model complexity,and improve model inter-pretability,Pearson correlation analysis(PsCA),Boruta,and recursive feature elimination(RFE)were employed to optimize features.Combined with the optimized features,the soil organic carbon density(SOCD)prediction model was constructed by using multivariate adaptive regression splines(MARS),extreme gradient boosting(XGBoost),and random forest(RF)algorithms and applied to predict the spatial distribution of SOCD and estimate the SOCS of different wetland types in 2020.The results show that:1)different feature combinations have a significant influence on the model performance.Better prediction performance was attained by building a model using RFE-based feature combinations.RF has the best prediction accuracy(R^(2)=0.587,RMSE=0.798 kg/m^(2),MAE=0.660 kg/m^(2)).2)Optical features are more important than radar and geo-climatic features in the MARS,XGBoost,and RF algorithms.3)The size of SOCS is related to SOCD and the area of each wetland type,aquaculture pond has the highest SOCS,followed by marsh,salt pan,mud-flat,and sand shore.展开更多
Atomic-scale strain mapping has become increasingly vital for investigating deformation mechanisms and the governing principles of solid materials.This is due to the significant impact of atomic-scale strain on the ph...Atomic-scale strain mapping has become increasingly vital for investigating deformation mechanisms and the governing principles of solid materials.This is due to the significant impact of atomic-scale strain on the physical,chemical,and mechanical properties of nanomaterials that comprise functional devices such as nanoelectronics,communication devices,electromechanical systems,and sensors.The advent of advanced electron microscopes has enabled the acquisition of high-magnification images with atomic resolution,providing an exceptional platform for measuring the atomic-scale strain of solid materials.However,accurate and unified strain mapping methods and standards for evaluating atomic-scale strain distribution remain scarce.Consequently,a unified strain mapping framework is proposed for atomic-scale strain measurement.Utilizing finite deformation analysis and the least-squares mathematical method,two types of atomic-scale strain field mapping methods have been developed,including the phase analysis-based methods(PAD and PAS)and the peak matching-based strain mapping method(PMS)for high-resolution scanning transmission electron microscope images.The prototypical 2D materials,graphene and molybdenum disulfide,serve as the subjects for the strain field mapping research,conducted through both simulation and experimentation.Upon comparing the theoretical strain mapping results of single-layer graphene and molybdenum disulfide with and without defects,it is demonstrated that the proposed strain mapping methods,particularly the PMS method,can accurately describe the large deformation surrounding a significant strain gradient.展开更多
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.展开更多
With rapid advancements in Infra-Red (IR) detection techniques, the range from where the IR-guided missiles are able to lock the target aircraft has increased. To avoid the detection and tracking by modern IR-guided m...With rapid advancements in Infra-Red (IR) detection techniques, the range from where the IR-guided missiles are able to lock the target aircraft has increased. To avoid the detection and tracking by modern IR-guided missiles, the aircraft and helicopters also demand progress in its stealth techniques. Hence, study of Infra-Red Signature Suppression (IRSS) systems in aircraft and helicopters has become vital even in design stage. Optical blocking (masking) is one of the effective IRSS techniques used to block the Line- Of-Sight (LOS) of the hot engine parts of the exhaust geometry. This paper reviews the various patents on IR signature suppression systems based on the optical blocking method or a combination of IRSS techniques. The performance penalties generated due to installation of various IRSS methods in aircraft and helicopters are also discussed.展开更多
Endocrine disrupting chemicals(EDCs) in the secondary effluent discharged from wastewater treatment plants are of great concern when water reuse is intended. Ozonation and ultrafiltration(UF) are powerful technologies...Endocrine disrupting chemicals(EDCs) in the secondary effluent discharged from wastewater treatment plants are of great concern when water reuse is intended. Ozonation and ultrafiltration(UF) are powerful technologies reported to eliminate EDCs. Due to the importance of effluent organic matters(EfOMs) in secondary effluent, the effects of three kinds of EfOM on the treatment of five EDCs using ozonation and UF were investigated. The three kinds of EfOM studied were humic acid sodium salt(NaAH), bovine serum albumin(BSA)and sodium alginate(NaAg); and the five EDCs were estrone, 17β-estradiol, estriol, 17α-ethynyl estradiol and bisphenol A. The results showed that EfOM accelerated the decay rate of ozone and inhibited the degradation efficiency of EDCs by ozonation in the order NaAH > BSA > NaAg.The ultraviolet absorbance at 280 nm(UVA_(280)) has potential for use as a surrogate indicator to assess EDC removal via ozonation without conducting difficult EDC analyses. When the decline in UVA_(280) exceeded 18%, the five EDCs had been completely removed. The UF behavior of NaAH,BSA and NaAg was found to follow the cake filtration law. The fouling potential of EfOM followed the order NaAg > NaAH > BSA; while EfOM on the membrane surface enhanced EDC removal in the order NaAH > BSA > NaAg. The mean retention rate of the membrane was increased by 24%, 10% and 8%, respectively. The properties of EDCs and EfOM cakes both influenced the EDC removal rates due to adsorption, size exclusion and charge attraction.展开更多
The interplay between chemistry and interfacial-tension-driven hydrodynamic instabilities has been studied experimentally. The system on hand consists of two immiscible liquids separated along an initially plane inter...The interplay between chemistry and interfacial-tension-driven hydrodynamic instabilities has been studied experimentally. The system on hand consists of two immiscible liquids separated along an initially plane interface at which an interfacial reaction takes place to produce in situ a surfactant. It is identified that the dynamics of the system depends on the orientation of the Hele-Shaw cell with respect to the vector of gravity. If the nele-Shaw cell is placed vertically, Marangoni cells with vigorous convection develop in both phases along a nearly planar interface. However, if the Hele-Shaw cell is tilted off the gravity, the instabilities in the system are characterized by the large scale interracial deformation with a spatio-temporal periodicity together with the chemo-Marangoni convection. The focus is on the exploration of the transition from the cellular mode to the large scale interfacial deformation.展开更多
Single thermal cycle simulation tests were carried out for X80 high strength steel pipes from three steel mills by a Gleeble 3500HS thermal simulation test machine,and coincidence degree of the thermal simulation curv...Single thermal cycle simulation tests were carried out for X80 high strength steel pipes from three steel mills by a Gleeble 3500HS thermal simulation test machine,and coincidence degree of the thermal simulation curve with the set curve under heat inputs of 6–30 kJ/cm was observed;The relationship between different heat inputs and microstructure,impact toughness and hardness of steel pipe CGHAZ(coarse grain heat affected zone)was studied by metallographic examination,impact test and hardness test.The results show that with the increase of heat input,original austenite grain size increases gradually,the lath bainite ratio decreases and the granular bainite ratio increases.The impact toughness of C steel pipe is lower than those of A and B steel pipe,and the impact toughness of CGHAZ from the three steel pipes show different trends:for A steel pipe CGHAZ,impact toughness increases first and then decreases,with the highest value of 270–320 J under 20–25 kJ/cm;for B steel pipe CGHAZ,impact toughness decreases slightly;for C steel pipe CGHAZ,impact toughness increases,with the highest value of 260–300 J under 25 kJ/cm.As the heat input increases,the hardness of three X80 steel pipes CGHAZ shows a decreasing trhighend,and C steel pipe has the largest decreasing range.展开更多
The recently developed discrete Boltzmann method(DBM), which is based on a set of uniform linear evolution equations and has high parallel efficiency, is employed to investigate the dynamic nonequilibrium process of K...The recently developed discrete Boltzmann method(DBM), which is based on a set of uniform linear evolution equations and has high parallel efficiency, is employed to investigate the dynamic nonequilibrium process of Kelvin-Helmholtz instability(KHI). It is found that, the relaxation time always strengthens the global nonequilibrium(GNE), entropy of mixing, and free enthalpy of mixing. Specifically, as a combined effect of physical gradients and nonequilibrium area, the GNE intensity first increases but decreases during the whole life-cycle of KHI. The growth rate of entropy of mixing shows firstly reducing, then increasing, and finally decreasing trends during the KHI process. The trend of the free enthalpy of mixing is opposite to that of the entropy of mixing. Detailed explanations are:(i) Initially,binary diffusion smooths quickly the sharp gradient in the mole fraction, which results in a steeply decreasing mixing rate.(ii) Afterwards, the mixing process is significantly promoted by the increasing length of material interface in the evolution of the KHI.(iii) As physical gradients are smoothed due to the binary diffusion and dissipation, the mixing rate reduces and approaches zero in the final stage. Moreover, with the increasing Atwood number, the global strength of viscous stresses on the heavy(light) medium reduces(increases), because the heavy(light) medium has a relatively small(large) velocity change. Furthermore, for a smaller Atwood number, the peaks of nonequilibrium manifestations emerge earlier, the entropy of mixing and free enthalpy of mixing change faster, because the KHI initiates a higher growth rate.展开更多
Amorphous Ni-P coating was plated on AZ31 magnesium alloy via the electroless plating technique, and the plated alloy was subsequently annealed. X-ray dif- fraction (XRD), scanning electron microscopy (SEM), energ...Amorphous Ni-P coating was plated on AZ31 magnesium alloy via the electroless plating technique, and the plated alloy was subsequently annealed. X-ray dif- fraction (XRD), scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS), and differential scanning calorimetry (DSC) were used to characterize the coating. The results show that the hardness of the coating is much higher than that of bare magnesium alloy, which further increases after crystallization. The electrochemical polarization and salt spray tests show that the coating exhibits a much higher corrosion resistance than that of the bare magnesium alloy. Moreover, the crystallized coating still exhibits a much stronger corrosion resistance than that of the bare magnesium alloy, although its corrosion resistance is lower than that of the as-plated one.展开更多
We present an improved lattice Boltzmann(LB) model for thermal liquid-vapor system.In the new model,the Windowed Fast Fourier Transform(WFFT) and its inverse are used to calculate both the convection term and the exte...We present an improved lattice Boltzmann(LB) model for thermal liquid-vapor system.In the new model,the Windowed Fast Fourier Transform(WFFT) and its inverse are used to calculate both the convection term and the external force term of the LB equation.By adopting the WFFT scheme,Gibbs oscillations can be damped effectively in unsmooth regions while high resolution feature of the spectral method can be retained in smooth regions.As a result,spatial discretization errors are dramatically decreased,conservation of the total energy is much better preserved,and the spurious velocities near the liquid-vapor interface are significantly reduced.The high resolution,together with the low complexity of the WFFT approach,endows the proposed method with considerable potential for studying a wide class of problems in the field of multiphase flows.展开更多
2024 aluminum alloys were consolidated by using spark plasma sintering(SPS) method,and then heat treated by solid solution treatment(SST) and aging treatment(AT) procedures.The average grain size of bulk samples sinte...2024 aluminum alloys were consolidated by using spark plasma sintering(SPS) method,and then heat treated by solid solution treatment(SST) and aging treatment(AT) procedures.The average grain size of bulk samples sintered with 5,20 and 50 μm powders was 3.72,4.73 and 8.11 μm,respectively.The difference between the average grain size and original powder size was attributed to the recrystallization during short sintering process.The number of the inclusion phases in these samples decreased after SST and increased after subsequent AT.Besides,it was observed that intergranular corrosion(IGC) cracks initiated from stable pits due to the electrochemical inhomogeneity between the intermetallic particles(IMPs) and the aluminum matrix(176.02,110.83 and 164.80 mV for as-SPS,as-SST and as-AT samples,respectively).Besides,the cracks would propagate along the grain boundaries(GBs) and bypass the IMPs at GBs during propagation.It was revealed that the sample after SST presented the best IGC resistance,and this was ascribed to the reduce of IMPs,both in size and number.展开更多
We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsu...We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.展开更多
Fine zirconium diboride (ZrB2) powders with high purity were successfully prepared by combustion synthesis through magnesiothermic reduction process in Mg-B2O3-ZrO2 system. The reaction mechanism was investigated by...Fine zirconium diboride (ZrB2) powders with high purity were successfully prepared by combustion synthesis through magnesiothermic reduction process in Mg-B2O3-ZrO2 system. The reaction mechanism was investigated by differential thermal analysis and quenching experiment. The results show that the whole magnesio-thermic reduction process includes three stages: first, molten B2O3 and Mg formed above the temperature of 650 ℃, and glassy B2O3 and solid ZrO2 particles were coated on the surface of the molten Mg; thus, the hollow balls can be formed when the molten Mg was exuded under capillary function. Second, ZrO2 particles reacted with molten Mg to form Zr and MgO with dissolution-precip-itation mechanism, which released a large amount of heat to induce the diffusion reaction between B203 and Mg to form B and MgO. Last, Zr reacted with B to form ZrB2 grains. The preparation of ZrB2 by self-propagating syn-thesis in Mg-B2O3-ZrO2 system is a solid-liquid-liquid reaction.展开更多
Potassium ion capacitors(PICs)are regarded as promising large-scale aqueous energy storage systems.However,due to the poor K^(+)transport kinetics and the structural instability of the cathode materials,the key issues...Potassium ion capacitors(PICs)are regarded as promising large-scale aqueous energy storage systems.However,due to the poor K^(+)transport kinetics and the structural instability of the cathode materials,the key issues of limited energy density and poor cyclic stability are obstacles to the in-depth growth of PICs.Herein,a novel O-doped perovskite fluoride is demonstrated via an in-situ electrochemical oxidation strategy as the cathode for PICs,introducing additional defects that improve the capacitance and facilitate the reaction kinetics of the electrode.During the electrochemical oxidation process,it is discovered that the perovskite fluoride crystal tends to transform into disordered O-doped KMnF 3(K_(x)MnF_(y)O_(z)),realizing a structural reconstruction at the electrode material/electrolyte interface.The First-principles calculations based on density functional theory(DFT)are performed to confirm that the improved electrical conduc-tivity and low ionic adsorption energy may be ascribed to the substitution of oxygen for fluorine.The obtained K_(1.14)MnF_(1.17)O_(1.26) cathode achieves a high specific capacitance of 694 F g^(-1) at 1 A g^(-1),as well as high capacitance retention of 91.3%after 10,000 charge/discharge cycles in mild K_(2)SO_(4) electrolyte.This study provides an effective strategy to improve the capacitive performance of perovskite fluoride cathode materials in electrochemical energy storage.展开更多
In order to establish the quantitative relationship between equivalent strain and the performance index of the deformed material within the range of certain passes for equal channel angular processing (ECAP), a new ...In order to establish the quantitative relationship between equivalent strain and the performance index of the deformed material within the range of certain passes for equal channel angular processing (ECAP), a new approach to characterize the equivalent strain was proposed. The results show that there exists better accordance between mechanical property (such as hardness or strength) and equivalent strain after rolling and ECAP in a certain range of deformation amount, and Gauss equation can be satisfied among the equivalent strain and the mechanical properties for ECAP. Through regression analysis on the data of hardness and strength after the deformation, a more generalized expression of equivalent strain for ECAP is proposed as:ε=k0exp[-(k1M-k2)^2], where M is the strength or hardness of the material, k1 is the modified coefficient (k1∈ (0, 1)), ko and k2 are two parameters dependent on the critical strain and mechanical property that reaches saturation state for the material, respectively. In this expression the equivalent strain for ECAP is characterized novelly through the mechanical parameter relating to material property rather than the classical geometry equation.展开更多
A new dual-polarized staggered and stacked patches antenna with wide impedance band-width and high isolation is proposed.The antenna consists of two groups of radiation patches,in 7 layers,and uses the orthogonal adja...A new dual-polarized staggered and stacked patches antenna with wide impedance band-width and high isolation is proposed.The antenna consists of two groups of radiation patches,in 7 layers,and uses the orthogonal adjacent coupling structure on staggered layer to excite a pair of linear polarization modes.Thanks to the staggered feeder mode,it has increased the isolation performance be-tween ports and compressed the transverse size of the antenna.As a result of the combination of staggered stack-up between the patches and the stepped gradient shape of the main radiating patches,it has effectively expanded the impedance bandwidth of the antenna.The proposed antenna is simulated,fabricated and measured.The staggered feeding structure effectively reduces the cross-sectional area of the antenna,and greatly improves the isolation between feeding ports.The measurement results show that the impedance bandwidths for vertical and horizontal polarization modes are 40.2%(638-960 MHz)and 40.0%(645-968 MHz)respectively when the return loss is lower than-10 dB,and the isolation between feeding ports is better than-30 dB.Meanwhile,the antenna has a stable and symmetrical radiation pattern across the working band,therefore making it suitable to be used as antenna and antenna array element of mobile wireless communication base stations.展开更多
Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow...Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from 0 to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.展开更多
In order to analyze the origin of carbon monoxide(CO) in coal seams, stress–strain experiments under temperature of 50, 150 and 250 °C were conducted using lignite from Kailuan mining area. Fourier transform inf...In order to analyze the origin of carbon monoxide(CO) in coal seams, stress–strain experiments under temperature of 50, 150 and 250 °C were conducted using lignite from Kailuan mining area. Fourier transform infrared spectroscopy and elemental analysis were carried out before and after deformation of the samples. The results indicated that CO generated at 150 and 250 °C; the gas component was mostly oxygen(O_2), with small amount of carbon dioxide(CO_2), methane(CH_4) and hydrogen(H_2). At 50 °C, O_2 and a little CO_2 were observed and no CO was found. The carbon content of the coal samples increased slightly after deformation, and the oxygen content, H/C ratio, and O/C ratio decreased. The molecular structure of coal displayed different evolution characteristics at various temperatures. At 50 and 150 °C, the falling off of side chains, broken of ether bond and directional realignment of the aliphatic chains resulting in the formation of long chains were the main performance of coal molecular structure evolution. While at 250 °C, the side chains fell off and short chains formed. Furthermore, at both 150 and 250 °C, condensed degree of aromatic ring increased. Under the action of temperature and pressure, CO forms in two ways.The first is that ether bond breaks, oxygen and carbon atoms combine together and forms CO, or O_2 forming in the broken of ether–oxygen bond leads to the oxidation of free radicals and resulting in the formation of CO. And the second is that CO derives from falling off of C=O group.展开更多
基金supported by the Science&Technology Development Fund of Tianjin Education Commission for Higher Education(No.2023KJ250).
文摘In the framework of elastoplastic theory,by introducing dissipative plastic energy(instead of cumulative plastic strain)and dissipative plastic energy rate(instead of cumulative plastic strain rate)into the ratchetting parameter evolution equation and isotropic evolution rules respectively,a cyclic elastoplastic constitutive model based on dissipative plastic energy is established.This model,termed the WDP model,describes the physical meaning and evolution rule of the unclosed stress–strain hysteresis loop using an energy method.A comparison of numerical implementation results with experimental data demonstrates the capability of the WDP model to predict the cyclic deformation of EA4T steel,effectively capturing the cyclic softening characteristics and ratchetting behaviors of axle steel EA4T.
基金supported by the National Key R&D Program of China(No.2020YFC2201100)the National Natural Science Foundation of China(Nos.12175032,12102082,12275044,12402327,12405290 and 12211530449)+4 种基金the Joint Program of the Science and Technology Program of Liaoning,China(No.2023JH2/101700285)the Fundamental Research Funds for the Central Universities of China(Nos.DUT22RC(3)078,DUT23RC(3)040 and DUT24ZD106)the S&T Program of Hebei,China(No.246Z2301G)the S&T Innovation Program of Hebei,China(Nos.SJMYF2022X18 and SJMYF2022X06)the Beijing Engineering Research Center of Efficient and Green Aerospace Propulsion Technology and Advanced Space Propulsion Laboratory of BICE,China(No.LabASP-2023-07).
文摘The electrospray thruster supplied by ionic liquid is a promising micro-propulsion thruster with small size and precise thrust, which can emit both cations and anions to achieve self-neutralization. In order to further investigate the effect of ion solvation energy on the evaporation of cations and anions from ionic liquid under the action of a uniform electric field, this paper establishes a transient Electrohydrodynamic (EHD) model for free ionic liquid droplets undergoing ion evaporation. The dynamic processes of droplet deformation and ion evaporation are simulated. And the study further focuses on the influence of different ion solvation energies for cations on the droplet morphology and the ion evaporation characteristics at the positively charged end and negatively charged end of the droplet. The results indicate that, when the ion solvation energy for cations is higher than that of anions, it will cause the ion evaporation at the positively charged end of the droplet to lag behind the ion evaporation at the negatively charged end. And the higher the ion solvation energy for the cations, the longer the evaporation lag time at the positively charged end of the droplet, which will lead to a higher peak of surface charge density that can be reached, resulting in a larger evaporation current and sharper droplet stretching deformation. Additionally, the peak surface charge density of the positively charged end of the droplet is linearly related to the ion solvation energy for cations, while the peak surface charge density of the negatively charged end remains almost unchanged and is not significantly affected by the ion solvation energy for cations.
基金Under the auspices of National Natural Science Foundation of China(No.42101393,41901375,52274166)Hebei Natural Science Foundation(No.D2022209005,D2023209008)Central Guided Local Science and Technology Development Fund Project of Hebei Province(No.236Z3305G,246Z4201G)Key Research and Development Program of Science and Technology Plan of Tangshan,China(No.22150221J)。
文摘Coastal wetlands are crucial for the‘blue carbon sink’,significantly contributing to regulating climate change.This study util-ized 160 soil samples,35 remote sensing features,and 5 geo-climatic data to accurately estimate the soil organic carbon stocks(SOCS)in the coastal wetlands of Tianjin and Hebei,China.To reduce data redundancy,simplify model complexity,and improve model inter-pretability,Pearson correlation analysis(PsCA),Boruta,and recursive feature elimination(RFE)were employed to optimize features.Combined with the optimized features,the soil organic carbon density(SOCD)prediction model was constructed by using multivariate adaptive regression splines(MARS),extreme gradient boosting(XGBoost),and random forest(RF)algorithms and applied to predict the spatial distribution of SOCD and estimate the SOCS of different wetland types in 2020.The results show that:1)different feature combinations have a significant influence on the model performance.Better prediction performance was attained by building a model using RFE-based feature combinations.RF has the best prediction accuracy(R^(2)=0.587,RMSE=0.798 kg/m^(2),MAE=0.660 kg/m^(2)).2)Optical features are more important than radar and geo-climatic features in the MARS,XGBoost,and RF algorithms.3)The size of SOCS is related to SOCD and the area of each wetland type,aquaculture pond has the highest SOCS,followed by marsh,salt pan,mud-flat,and sand shore.
基金support from the National Natural Science Foundation of China through Grants 12172190,11872035,11632010,and 12302236。
文摘Atomic-scale strain mapping has become increasingly vital for investigating deformation mechanisms and the governing principles of solid materials.This is due to the significant impact of atomic-scale strain on the physical,chemical,and mechanical properties of nanomaterials that comprise functional devices such as nanoelectronics,communication devices,electromechanical systems,and sensors.The advent of advanced electron microscopes has enabled the acquisition of high-magnification images with atomic resolution,providing an exceptional platform for measuring the atomic-scale strain of solid materials.However,accurate and unified strain mapping methods and standards for evaluating atomic-scale strain distribution remain scarce.Consequently,a unified strain mapping framework is proposed for atomic-scale strain measurement.Utilizing finite deformation analysis and the least-squares mathematical method,two types of atomic-scale strain field mapping methods have been developed,including the phase analysis-based methods(PAD and PAS)and the peak matching-based strain mapping method(PMS)for high-resolution scanning transmission electron microscope images.The prototypical 2D materials,graphene and molybdenum disulfide,serve as the subjects for the strain field mapping research,conducted through both simulation and experimentation.Upon comparing the theoretical strain mapping results of single-layer graphene and molybdenum disulfide with and without defects,it is demonstrated that the proposed strain mapping methods,particularly the PMS method,can accurately describe the large deformation surrounding a significant strain gradient.
基金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.
基金the Indian Institute of Technology Bombay’s Post-Doctoral Research Program, vide appointment no. AO/Admn1/33/2018 dated 10.Aug’2018 for providing funding
文摘With rapid advancements in Infra-Red (IR) detection techniques, the range from where the IR-guided missiles are able to lock the target aircraft has increased. To avoid the detection and tracking by modern IR-guided missiles, the aircraft and helicopters also demand progress in its stealth techniques. Hence, study of Infra-Red Signature Suppression (IRSS) systems in aircraft and helicopters has become vital even in design stage. Optical blocking (masking) is one of the effective IRSS techniques used to block the Line- Of-Sight (LOS) of the hot engine parts of the exhaust geometry. This paper reviews the various patents on IR signature suppression systems based on the optical blocking method or a combination of IRSS techniques. The performance penalties generated due to installation of various IRSS methods in aircraft and helicopters are also discussed.
基金supported by the National Natural Science Foundation of China (No.51808008)the Natural Science Foundation of Hebei Province (No.E2017409004)the Youth Foundation of Hebei Province Education Department (No.QN2016224)
文摘Endocrine disrupting chemicals(EDCs) in the secondary effluent discharged from wastewater treatment plants are of great concern when water reuse is intended. Ozonation and ultrafiltration(UF) are powerful technologies reported to eliminate EDCs. Due to the importance of effluent organic matters(EfOMs) in secondary effluent, the effects of three kinds of EfOM on the treatment of five EDCs using ozonation and UF were investigated. The three kinds of EfOM studied were humic acid sodium salt(NaAH), bovine serum albumin(BSA)and sodium alginate(NaAg); and the five EDCs were estrone, 17β-estradiol, estriol, 17α-ethynyl estradiol and bisphenol A. The results showed that EfOM accelerated the decay rate of ozone and inhibited the degradation efficiency of EDCs by ozonation in the order NaAH > BSA > NaAg.The ultraviolet absorbance at 280 nm(UVA_(280)) has potential for use as a surrogate indicator to assess EDC removal via ozonation without conducting difficult EDC analyses. When the decline in UVA_(280) exceeded 18%, the five EDCs had been completely removed. The UF behavior of NaAH,BSA and NaAg was found to follow the cake filtration law. The fouling potential of EfOM followed the order NaAg > NaAH > BSA; while EfOM on the membrane surface enhanced EDC removal in the order NaAH > BSA > NaAg. The mean retention rate of the membrane was increased by 24%, 10% and 8%, respectively. The properties of EDCs and EfOM cakes both influenced the EDC removal rates due to adsorption, size exclusion and charge attraction.
基金Deutsche Forschung Gemainschaft(Ec/201/1-5)Deutsches Zentrum fuer Luft und Raumfahrt(50WM0058).
文摘The interplay between chemistry and interfacial-tension-driven hydrodynamic instabilities has been studied experimentally. The system on hand consists of two immiscible liquids separated along an initially plane interface at which an interfacial reaction takes place to produce in situ a surfactant. It is identified that the dynamics of the system depends on the orientation of the Hele-Shaw cell with respect to the vector of gravity. If the nele-Shaw cell is placed vertically, Marangoni cells with vigorous convection develop in both phases along a nearly planar interface. However, if the Hele-Shaw cell is tilted off the gravity, the instabilities in the system are characterized by the large scale interracial deformation with a spatio-temporal periodicity together with the chemo-Marangoni convection. The focus is on the exploration of the transition from the cellular mode to the large scale interfacial deformation.
基金supported by Science and Technology Research Project of Universities of Hebei Province(No.QN201.221)。
文摘Single thermal cycle simulation tests were carried out for X80 high strength steel pipes from three steel mills by a Gleeble 3500HS thermal simulation test machine,and coincidence degree of the thermal simulation curve with the set curve under heat inputs of 6–30 kJ/cm was observed;The relationship between different heat inputs and microstructure,impact toughness and hardness of steel pipe CGHAZ(coarse grain heat affected zone)was studied by metallographic examination,impact test and hardness test.The results show that with the increase of heat input,original austenite grain size increases gradually,the lath bainite ratio decreases and the granular bainite ratio increases.The impact toughness of C steel pipe is lower than those of A and B steel pipe,and the impact toughness of CGHAZ from the three steel pipes show different trends:for A steel pipe CGHAZ,impact toughness increases first and then decreases,with the highest value of 270–320 J under 20–25 kJ/cm;for B steel pipe CGHAZ,impact toughness decreases slightly;for C steel pipe CGHAZ,impact toughness increases,with the highest value of 260–300 J under 25 kJ/cm.As the heat input increases,the hardness of three X80 steel pipes CGHAZ shows a decreasing trhighend,and C steel pipe has the largest decreasing range.
基金Supported by the Natural Science Foundation of China under Grant Nos.91441120,51806116,11875001,and 11602162the China Postdoctoral Science Foundation under Grant No.2017M620757+2 种基金the Center for Combustion Energy at Tsinghua Universitythe Natural Science Foundation of Hebei Province under Grant Nos.A2017409014,ZD2017001,and A201500111the UK Engineering and Physical Sciences Research Council under Project UK Consortium on Mesoscale Engineering Sciences(UKCOMES)under Grant Nos.EP/L00030X/1 and EP/R029598/1
文摘The recently developed discrete Boltzmann method(DBM), which is based on a set of uniform linear evolution equations and has high parallel efficiency, is employed to investigate the dynamic nonequilibrium process of Kelvin-Helmholtz instability(KHI). It is found that, the relaxation time always strengthens the global nonequilibrium(GNE), entropy of mixing, and free enthalpy of mixing. Specifically, as a combined effect of physical gradients and nonequilibrium area, the GNE intensity first increases but decreases during the whole life-cycle of KHI. The growth rate of entropy of mixing shows firstly reducing, then increasing, and finally decreasing trends during the KHI process. The trend of the free enthalpy of mixing is opposite to that of the entropy of mixing. Detailed explanations are:(i) Initially,binary diffusion smooths quickly the sharp gradient in the mole fraction, which results in a steeply decreasing mixing rate.(ii) Afterwards, the mixing process is significantly promoted by the increasing length of material interface in the evolution of the KHI.(iii) As physical gradients are smoothed due to the binary diffusion and dissipation, the mixing rate reduces and approaches zero in the final stage. Moreover, with the increasing Atwood number, the global strength of viscous stresses on the heavy(light) medium reduces(increases), because the heavy(light) medium has a relatively small(large) velocity change. Furthermore, for a smaller Atwood number, the peaks of nonequilibrium manifestations emerge earlier, the entropy of mixing and free enthalpy of mixing change faster, because the KHI initiates a higher growth rate.
基金supported by the National Basic Research Program of China (No. 2010CB731600)the National Nature Science Foundation of China (No. 51121061/51171163)the Appling Basic Research Program of Hebei (No. 10965124D)
文摘Amorphous Ni-P coating was plated on AZ31 magnesium alloy via the electroless plating technique, and the plated alloy was subsequently annealed. X-ray dif- fraction (XRD), scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS), and differential scanning calorimetry (DSC) were used to characterize the coating. The results show that the hardness of the coating is much higher than that of bare magnesium alloy, which further increases after crystallization. The electrochemical polarization and salt spray tests show that the coating exhibits a much higher corrosion resistance than that of the bare magnesium alloy. Moreover, the crystallized coating still exhibits a much stronger corrosion resistance than that of the bare magnesium alloy, although its corrosion resistance is lower than that of the as-plated one.
基金Supported by the Science Foundations of Laboratory of Computational Physics and China Academy of Engineering Physics under Grant Nos.2009A0102005 and 2011A0201002National Natural Science Foundation of China under Grant Nos.11075021,91130020,11074300, and 10874242+2 种基金Fundamental research funds for the central university under Grant No.2010YS03Technology Support Program of LangFang under Grant Nos.2010011030,201101118/21/23/24Teaching and Research Foundation of NCIAE under Grant Nos.JY-2011-027-Y and JY-2011-028-Y
文摘We present an improved lattice Boltzmann(LB) model for thermal liquid-vapor system.In the new model,the Windowed Fast Fourier Transform(WFFT) and its inverse are used to calculate both the convection term and the external force term of the LB equation.By adopting the WFFT scheme,Gibbs oscillations can be damped effectively in unsmooth regions while high resolution feature of the spectral method can be retained in smooth regions.As a result,spatial discretization errors are dramatically decreased,conservation of the total energy is much better preserved,and the spurious velocities near the liquid-vapor interface are significantly reduced.The high resolution,together with the low complexity of the WFFT approach,endows the proposed method with considerable potential for studying a wide class of problems in the field of multiphase flows.
基金financially supported by the National Nature Science Foundations of China (Nos.51271012 and 51671013)。
文摘2024 aluminum alloys were consolidated by using spark plasma sintering(SPS) method,and then heat treated by solid solution treatment(SST) and aging treatment(AT) procedures.The average grain size of bulk samples sintered with 5,20 and 50 μm powders was 3.72,4.73 and 8.11 μm,respectively.The difference between the average grain size and original powder size was attributed to the recrystallization during short sintering process.The number of the inclusion phases in these samples decreased after SST and increased after subsequent AT.Besides,it was observed that intergranular corrosion(IGC) cracks initiated from stable pits due to the electrochemical inhomogeneity between the intermetallic particles(IMPs) and the aluminum matrix(176.02,110.83 and 164.80 mV for as-SPS,as-SST and as-AT samples,respectively).Besides,the cracks would propagate along the grain boundaries(GBs) and bypass the IMPs at GBs during propagation.It was revealed that the sample after SST presented the best IGC resistance,and this was ascribed to the reduce of IMPs,both in size and number.
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos.2009A0102005 and 2009B0101012the National Basic Research Program (973 Program) under Grant No.2007CB815105the National Natural Science Foundation under Grant Nos.10775018,10702010,and 10775088
文摘We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.
基金financially supported by the National Natural Science Foundation of China(No.91016014)
文摘Fine zirconium diboride (ZrB2) powders with high purity were successfully prepared by combustion synthesis through magnesiothermic reduction process in Mg-B2O3-ZrO2 system. The reaction mechanism was investigated by differential thermal analysis and quenching experiment. The results show that the whole magnesio-thermic reduction process includes three stages: first, molten B2O3 and Mg formed above the temperature of 650 ℃, and glassy B2O3 and solid ZrO2 particles were coated on the surface of the molten Mg; thus, the hollow balls can be formed when the molten Mg was exuded under capillary function. Second, ZrO2 particles reacted with molten Mg to form Zr and MgO with dissolution-precip-itation mechanism, which released a large amount of heat to induce the diffusion reaction between B203 and Mg to form B and MgO. Last, Zr reacted with B to form ZrB2 grains. The preparation of ZrB2 by self-propagating syn-thesis in Mg-B2O3-ZrO2 system is a solid-liquid-liquid reaction.
基金the financial support from Liaoning Sci-ence and Technology Development Foundation Guided by Cen-tral Government(No.2021JH6/10500139)the Fundamental Research Funds for the Central Universities(No.N2205003)+2 种基金the financial support from the National Natural Science Foundation of China(No.52003007)Nat-ural Science Foundation of Hebei Province(No.E2019409063)Langfang top-notch talent(No.LFBJ202004).
文摘Potassium ion capacitors(PICs)are regarded as promising large-scale aqueous energy storage systems.However,due to the poor K^(+)transport kinetics and the structural instability of the cathode materials,the key issues of limited energy density and poor cyclic stability are obstacles to the in-depth growth of PICs.Herein,a novel O-doped perovskite fluoride is demonstrated via an in-situ electrochemical oxidation strategy as the cathode for PICs,introducing additional defects that improve the capacitance and facilitate the reaction kinetics of the electrode.During the electrochemical oxidation process,it is discovered that the perovskite fluoride crystal tends to transform into disordered O-doped KMnF 3(K_(x)MnF_(y)O_(z)),realizing a structural reconstruction at the electrode material/electrolyte interface.The First-principles calculations based on density functional theory(DFT)are performed to confirm that the improved electrical conduc-tivity and low ionic adsorption energy may be ascribed to the substitution of oxygen for fluorine.The obtained K_(1.14)MnF_(1.17)O_(1.26) cathode achieves a high specific capacitance of 694 F g^(-1) at 1 A g^(-1),as well as high capacitance retention of 91.3%after 10,000 charge/discharge cycles in mild K_(2)SO_(4) electrolyte.This study provides an effective strategy to improve the capacitive performance of perovskite fluoride cathode materials in electrochemical energy storage.
基金Projects(50471102,50671089) supported by the National Natural Science Foundation of China
文摘In order to establish the quantitative relationship between equivalent strain and the performance index of the deformed material within the range of certain passes for equal channel angular processing (ECAP), a new approach to characterize the equivalent strain was proposed. The results show that there exists better accordance between mechanical property (such as hardness or strength) and equivalent strain after rolling and ECAP in a certain range of deformation amount, and Gauss equation can be satisfied among the equivalent strain and the mechanical properties for ECAP. Through regression analysis on the data of hardness and strength after the deformation, a more generalized expression of equivalent strain for ECAP is proposed as:ε=k0exp[-(k1M-k2)^2], where M is the strength or hardness of the material, k1 is the modified coefficient (k1∈ (0, 1)), ko and k2 are two parameters dependent on the critical strain and mechanical property that reaches saturation state for the material, respectively. In this expression the equivalent strain for ECAP is characterized novelly through the mechanical parameter relating to material property rather than the classical geometry equation.
基金Natural Science Founda-tion of China(Grant 61171051 and 61771063)Plan-ning projects of Hebei provincial department(Grant 15210410)Research Fund Project of North China Institute of Aerospace Engineering(Grant BKY-2021-13)。
文摘A new dual-polarized staggered and stacked patches antenna with wide impedance band-width and high isolation is proposed.The antenna consists of two groups of radiation patches,in 7 layers,and uses the orthogonal adjacent coupling structure on staggered layer to excite a pair of linear polarization modes.Thanks to the staggered feeder mode,it has increased the isolation performance be-tween ports and compressed the transverse size of the antenna.As a result of the combination of staggered stack-up between the patches and the stepped gradient shape of the main radiating patches,it has effectively expanded the impedance bandwidth of the antenna.The proposed antenna is simulated,fabricated and measured.The staggered feeding structure effectively reduces the cross-sectional area of the antenna,and greatly improves the isolation between feeding ports.The measurement results show that the impedance bandwidths for vertical and horizontal polarization modes are 40.2%(638-960 MHz)and 40.0%(645-968 MHz)respectively when the return loss is lower than-10 dB,and the isolation between feeding ports is better than-30 dB.Meanwhile,the antenna has a stable and symmetrical radiation pattern across the working band,therefore making it suitable to be used as antenna and antenna array element of mobile wireless communication base stations.
基金the Science Foundation of Laboratory of Computational Physics,the State Key Basic Research Program (973 Program) under Grant No.2007CB815105National Natural Science Foundation of China under Grant Nos.10775018,10474137,10702010,and 10604010
文摘Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from 0 to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.
基金financial support from the National Science foundation of China(No.41430317)the Discipline Innovative Engineering Plan sponsored by the Ministry of Education of China+1 种基金the State Administration of Foreign Experts Affairs of China(No.13023)the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘In order to analyze the origin of carbon monoxide(CO) in coal seams, stress–strain experiments under temperature of 50, 150 and 250 °C were conducted using lignite from Kailuan mining area. Fourier transform infrared spectroscopy and elemental analysis were carried out before and after deformation of the samples. The results indicated that CO generated at 150 and 250 °C; the gas component was mostly oxygen(O_2), with small amount of carbon dioxide(CO_2), methane(CH_4) and hydrogen(H_2). At 50 °C, O_2 and a little CO_2 were observed and no CO was found. The carbon content of the coal samples increased slightly after deformation, and the oxygen content, H/C ratio, and O/C ratio decreased. The molecular structure of coal displayed different evolution characteristics at various temperatures. At 50 and 150 °C, the falling off of side chains, broken of ether bond and directional realignment of the aliphatic chains resulting in the formation of long chains were the main performance of coal molecular structure evolution. While at 250 °C, the side chains fell off and short chains formed. Furthermore, at both 150 and 250 °C, condensed degree of aromatic ring increased. Under the action of temperature and pressure, CO forms in two ways.The first is that ether bond breaks, oxygen and carbon atoms combine together and forms CO, or O_2 forming in the broken of ether–oxygen bond leads to the oxidation of free radicals and resulting in the formation of CO. And the second is that CO derives from falling off of C=O group.