To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic...In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.展开更多
The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(R...The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.展开更多
Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of...Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.展开更多
We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projec...We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.展开更多
The ^(12)C+^(12)C fusion is one of the most important reactions in modern nuclear astrophysics.The trend and magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes.Ho...The ^(12)C+^(12)C fusion is one of the most important reactions in modern nuclear astrophysics.The trend and magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes.However,direct measurement of this reaction is extremely difficult,which makes it necessary to develop indirect methods.In this study,the ^(23)Na+p reaction system was used to study the compound nucleus ^(24)Mg.We employed a thick-target inverse kinematics method combined with theγ-charged-particle coincidence technique to measure the proton andα exit channels of ^(24)Mg.Technical details of the ^(23)Na+p thick-target inverse kinematics experiment and analysis are presented herein.展开更多
The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the co...The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.展开更多
Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analy...Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.展开更多
We consider the persistence of affine periodic solutions for perturbed affine periodic systems.Such(Q,T)-affine periodic solutions have the form x(t+T)=Qx(t)for all t∈R,where T>0 is fixed and Q is a nonsingular ma...We consider the persistence of affine periodic solutions for perturbed affine periodic systems.Such(Q,T)-affine periodic solutions have the form x(t+T)=Qx(t)for all t∈R,where T>0 is fixed and Q is a nonsingular matrix.These are a kind of spatiotemporal symmetric solutions,e.g.spiral waves.We give the averaging method for the existence of affine periodic solutions in two situations:one in which the initial values of the affine periodic solutions of the unperturbed system form a manifold,and another that does not rely on the structure of the initial values of the unperturbed system's affine periodic solutions.The transversal condition is determined using the Brouwer degree.We also present a higher order averaging method for general degenerate systems by means of the Brouwer degree and a Lyapunov-Schmidt reduction.展开更多
A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential ...A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.展开更多
This paper presents a comprehensive framework for analyzing phase transitions in collective models such as theVicsek model under various noise types. The Vicsek model, focusing on understanding the collective behavior...This paper presents a comprehensive framework for analyzing phase transitions in collective models such as theVicsek model under various noise types. The Vicsek model, focusing on understanding the collective behaviors of socialanimals, is known due to its discontinuous phase transitions under vector noise. However, its behavior under scalar noiseremains less conclusive. Renowned for its efficacy in the analysis of complex systems under both equilibrium and nonequilibriumstates, the eigen microstate method is employed here for a quantitative examination of the phase transitions inthe Vicsek model under both vector and scalar noises. The study finds that the Vicsek model exhibits discontinuous phasetransitions regardless of noise type. Furthermore, the dichotomy method is utilized to identify the critical points for thesephase transitions. A significant finding is the observed increase in the critical point for discontinuous phase transitions withescalation of population density.展开更多
The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state...The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state. The calculation results indicate that the motion of the droplet on the cylinder can be divided into three stages: spreading, sliding, and aggregating.The contact length and contact time of a droplet on a cylindrical surface can be affected by factors such as the wettability gradient of the cylindrical wall, the Bond number, and droplet size. Furthermore, phase diagrams showing the relationship between Bond number, cylinder wall wettability gradient, and contact time as well as maximum contact length for three different droplet sizes are given. A theoretical foundation for additional research into the heat and mass transfer process between the droplet and the cylinder can be established by comprehending the variable rules of maximum contact length and contact time.展开更多
For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of ...For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.展开更多
The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-...The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.展开更多
The selection and coordinated application of government innovation policies are crucial for guiding the direction of enterprise innovation and unleashing their innovation potential.However,due to the lengthy,voluminou...The selection and coordinated application of government innovation policies are crucial for guiding the direction of enterprise innovation and unleashing their innovation potential.However,due to the lengthy,voluminous,complex,and unstructured nature of regional innovation policy texts,traditional policy classification methods often overlook the reality that these texts cover multiple policy topics,leading to lack of objectivity.In contrast,topic mining technology can handle large-scale textual data,overcoming challenges such as the abundance of policy content and difficulty in classification.Although topic models can partition numerous policy texts into topics,they cannot analyze the interplay among policy topics and the impact of policy topic coordination on enterprise innovation in detail.Therefore,we propose a big data analysis scheme for policy coordination paths based on the latent Dirichlet allocation(LDA)model and the fuzzyset qualitative comparative analysis(fsQCA)method by combining topic models with qualitative comparative analysis.The LDA model was employed to derive the topic distribution of each document and the word distribution of each topic and enable automatic classi-fication through algorithms,providing reliable and objective textual classification results.Subsequently,the fsQCA method was used to analyze the coordination paths and dynamic characteristics.Finally,experimental analysis was conducted using innovation policy text data from 31 provincial-level administrative regions in China from 2012 to 2021 as research samples.The results suggest that the proposed method effectively partitions innovation policy topics and analyzes the policy configuration,driving enterprise innovation in different regions.展开更多
Copper-based azide(Cu(N_(3))2 or CuN_(3),CA)chips synthesized by in-situ azide reaction and utilized in miniaturized explosive systems has become a hot research topic in recent years.However,the advantages of in-situ ...Copper-based azide(Cu(N_(3))2 or CuN_(3),CA)chips synthesized by in-situ azide reaction and utilized in miniaturized explosive systems has become a hot research topic in recent years.However,the advantages of in-situ synthesis method,including small size and low dosage,bring about difficulties in quantitative analysis and differences in ignition capabilities of CA chips.The aim of present work is to develop a simplified quantitative analysis method for accurate and safe analysis of components in CA chips to evaluate and investigate the corresponding ignition ability.In this work,Cu(N_(3))2 and CuN_(3)components in CA chips were separated through dissolution and distillation by utilizing the difference in solubility and corresponding content was obtained by measuring N_(3)-concentration through spectrophotometry.The spectrophotometry method was optimized by studying influencing factors and the recovery rate of different separation methods was studied,ensuring the accuracy and reproducibility of test results.The optimized method is linear in range from 1.0-25.0 mg/L,with a correlation coefficient R^(2)=0.9998,which meets the requirements of CA chips with a milligram-level content test.Compared with the existing ICP method,component analysis results of CA chips obtained by spectrophotometry are closer to real component content in samples and have satisfactory accuracy.Moreover,as its application in miniaturized explosive systems,the ignition ability of CA chips with different component contents for direct ink writing CL-20 and the corresponding mechanism was studied.This study provided a basis and idea for the design and performance evaluation of CA chips in miniaturized explosive systems.展开更多
This article introduces the concept of load aggregation,which involves a comprehensive analysis of loads to acquire their external characteristics for the purpose of modeling and analyzing power systems.The online ide...This article introduces the concept of load aggregation,which involves a comprehensive analysis of loads to acquire their external characteristics for the purpose of modeling and analyzing power systems.The online identification method is a computer-involved approach for data collection,processing,and system identification,commonly used for adaptive control and prediction.This paper proposes a method for dynamically aggregating large-scale adjustable loads to support high proportions of new energy integration,aiming to study the aggregation characteristics of regional large-scale adjustable loads using online identification techniques and feature extraction methods.The experiment selected 300 central air conditioners as the research subject and analyzed their regulation characteristics,economic efficiency,and comfort.The experimental results show that as the adjustment time of the air conditioner increases from 5 minutes to 35 minutes,the stable adjustment quantity during the adjustment period decreases from 28.46 to 3.57,indicating that air conditioning loads can be controlled over a long period and have better adjustment effects in the short term.Overall,the experimental results of this paper demonstrate that analyzing the aggregation characteristics of regional large-scale adjustable loads using online identification techniques and feature extraction algorithms is effective.展开更多
As an important material for manufacturing resonant components of musical instruments,Paulownia has an important influence on the sound quality of Ruan.In this paper,a model for evaluating the sound quality of Ruan ba...As an important material for manufacturing resonant components of musical instruments,Paulownia has an important influence on the sound quality of Ruan.In this paper,a model for evaluating the sound quality of Ruan based on the vibration characteristics of wood is developed using machine learning methods.Generally,the selection of materials for Ruan manufacturing relies primarily on manually weighing,observing,striking,and listening by the instrument technician.Deficiencies in scientific theory have hindered the quality of the finished Ruan.In this study,nine Ruans were manufactured,and a prediction model of Ruan sound quality was proposed based on the raw material information of Ruans.Out of a total of 180 data sets,145 and 45 sets were chosen for training and validation,respec-tively.In this paper,typical correlation analysis was used to determine the correlation between two single indicators in two adjacent pairwise combinations of the measured objects in each stage of the production process in Ruan.The vibra-tion characteristics of the wood were tested,and a model for predicting the evaluation of Ruan’s acoustic qualities was developed by measuring the vibration characteristics of the resonating plate material.The acoustic quality of the Ruan sound board wood was evaluated and predicted using machine learning model generalized regression neural net-work.The results show that the prediction of Ruan sound quality can be achieved using Matlab simulation based on the vibration characteristics of the soundboard wood.When the model-predicted values were compared with the tradi-tional predicted results,it was found that the generalized regression neural network had good performance,achieving an accuracy of 93.8%which was highly consistent with the experimental results.It was concluded that the model can accurately predict the acoustic quality of the Ruan based on the vibration performance of the soundboards.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金The Guangdong Basic and Applied Basic Research Foundation(2022A1515010730)National Natural Science Foundation of China(32001647)+2 种基金National Natural Science Foundation of China(31972022)Financial and moral assistance supported by the Guangdong Basic and Applied Basic Research Foundation(2019A1515011996)111 Project(B17018)。
文摘In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.
基金supported by the National Key Research and Development Program of China(No.2023YFC2907600)the National Natural Science Foundation of China(Nos.42077267,42277174 and 52074164)+2 种基金the Natural Science Foundation of Shandong Province,China(No.ZR2020JQ23)the Opening Project of State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology(No.KFJJ21-02Z)the Fundamental Research Funds for the Central Universities,China(No.2022JCCXSB03).
文摘The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.
基金Supported by the National Natural Science Foundation of China(12001395)the special fund for Science and Technology Innovation Teams of Shanxi Province(202204051002018)+1 种基金Research Project Supported by Shanxi Scholarship Council of China(2022-169)Graduate Education Innovation Project of Taiyuan Normal University(SYYJSYC-2314)。
文摘Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.
基金financial support from the National Natural Science Foundation of China (Grant No. 12227901)the financial support from the National Natural Science Foundation of China (Grant Nos. 11974263 and 12174291)。
文摘We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.
基金supported by the National Key Research and Development Project (No. 2022YFA1602301)the National Natural Science Foundation of China (Nos. U2267205, 12275361, 12125509, 12222514, 11961141003, and 12005304)+2 种基金the CAST Young Talent Support Planthe CNNC Science Fund for Talented Young Scholarsthe Continuous-Support Basic Scientific Research Project
文摘The ^(12)C+^(12)C fusion is one of the most important reactions in modern nuclear astrophysics.The trend and magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes.However,direct measurement of this reaction is extremely difficult,which makes it necessary to develop indirect methods.In this study,the ^(23)Na+p reaction system was used to study the compound nucleus ^(24)Mg.We employed a thick-target inverse kinematics method combined with theγ-charged-particle coincidence technique to measure the proton andα exit channels of ^(24)Mg.Technical details of the ^(23)Na+p thick-target inverse kinematics experiment and analysis are presented herein.
文摘The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.
基金supported by the Na-tional Natural Science Foundation of China(No.52272369).
文摘Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.
基金supported by the National Natural Science Foundation of China(1237119112071175)+4 种基金supported by the NSFC(1207117511901080)supported by the NSFC(12071175)the Fundamental Research Funds For the Central Universities(2412023YQ003)the Natural Science Foundation of Jilin Province(20200201253JC)。
文摘We consider the persistence of affine periodic solutions for perturbed affine periodic systems.Such(Q,T)-affine periodic solutions have the form x(t+T)=Qx(t)for all t∈R,where T>0 is fixed and Q is a nonsingular matrix.These are a kind of spatiotemporal symmetric solutions,e.g.spiral waves.We give the averaging method for the existence of affine periodic solutions in two situations:one in which the initial values of the affine periodic solutions of the unperturbed system form a manifold,and another that does not rely on the structure of the initial values of the unperturbed system's affine periodic solutions.The transversal condition is determined using the Brouwer degree.We also present a higher order averaging method for general degenerate systems by means of the Brouwer degree and a Lyapunov-Schmidt reduction.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.52171251,U2106225,and 52231011)Dalian Science and Technology Innovation Fund (Grant No.2022JJ12GX036)。
文摘A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.
基金the National Natural Science Foundation of China(Grant No.62273033).
文摘This paper presents a comprehensive framework for analyzing phase transitions in collective models such as theVicsek model under various noise types. The Vicsek model, focusing on understanding the collective behaviors of socialanimals, is known due to its discontinuous phase transitions under vector noise. However, its behavior under scalar noiseremains less conclusive. Renowned for its efficacy in the analysis of complex systems under both equilibrium and nonequilibriumstates, the eigen microstate method is employed here for a quantitative examination of the phase transitions inthe Vicsek model under both vector and scalar noises. The study finds that the Vicsek model exhibits discontinuous phasetransitions regardless of noise type. Furthermore, the dichotomy method is utilized to identify the critical points for thesephase transitions. A significant finding is the observed increase in the critical point for discontinuous phase transitions withescalation of population density.
文摘The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state. The calculation results indicate that the motion of the droplet on the cylinder can be divided into three stages: spreading, sliding, and aggregating.The contact length and contact time of a droplet on a cylindrical surface can be affected by factors such as the wettability gradient of the cylindrical wall, the Bond number, and droplet size. Furthermore, phase diagrams showing the relationship between Bond number, cylinder wall wettability gradient, and contact time as well as maximum contact length for three different droplet sizes are given. A theoretical foundation for additional research into the heat and mass transfer process between the droplet and the cylinder can be established by comprehending the variable rules of maximum contact length and contact time.
基金supported by National Natural Science Foundation of China(11771257)the Shandong Provincial Natural Science Foundation of China(ZR2023YQ002,ZR2023MA007,ZR2021MA004)。
文摘For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.
基金supported by the National Key Research and Development Program of China(No.2020YFB1901900)the National Natural Science Foundation of China(Nos.U20B2011,12175138)the Shanghai Rising-Star Program。
文摘The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.
文摘The selection and coordinated application of government innovation policies are crucial for guiding the direction of enterprise innovation and unleashing their innovation potential.However,due to the lengthy,voluminous,complex,and unstructured nature of regional innovation policy texts,traditional policy classification methods often overlook the reality that these texts cover multiple policy topics,leading to lack of objectivity.In contrast,topic mining technology can handle large-scale textual data,overcoming challenges such as the abundance of policy content and difficulty in classification.Although topic models can partition numerous policy texts into topics,they cannot analyze the interplay among policy topics and the impact of policy topic coordination on enterprise innovation in detail.Therefore,we propose a big data analysis scheme for policy coordination paths based on the latent Dirichlet allocation(LDA)model and the fuzzyset qualitative comparative analysis(fsQCA)method by combining topic models with qualitative comparative analysis.The LDA model was employed to derive the topic distribution of each document and the word distribution of each topic and enable automatic classi-fication through algorithms,providing reliable and objective textual classification results.Subsequently,the fsQCA method was used to analyze the coordination paths and dynamic characteristics.Finally,experimental analysis was conducted using innovation policy text data from 31 provincial-level administrative regions in China from 2012 to 2021 as research samples.The results suggest that the proposed method effectively partitions innovation policy topics and analyzes the policy configuration,driving enterprise innovation in different regions.
基金the financial support provided by the National Natural Science Foundation of China(Grant No.11872013).
文摘Copper-based azide(Cu(N_(3))2 or CuN_(3),CA)chips synthesized by in-situ azide reaction and utilized in miniaturized explosive systems has become a hot research topic in recent years.However,the advantages of in-situ synthesis method,including small size and low dosage,bring about difficulties in quantitative analysis and differences in ignition capabilities of CA chips.The aim of present work is to develop a simplified quantitative analysis method for accurate and safe analysis of components in CA chips to evaluate and investigate the corresponding ignition ability.In this work,Cu(N_(3))2 and CuN_(3)components in CA chips were separated through dissolution and distillation by utilizing the difference in solubility and corresponding content was obtained by measuring N_(3)-concentration through spectrophotometry.The spectrophotometry method was optimized by studying influencing factors and the recovery rate of different separation methods was studied,ensuring the accuracy and reproducibility of test results.The optimized method is linear in range from 1.0-25.0 mg/L,with a correlation coefficient R^(2)=0.9998,which meets the requirements of CA chips with a milligram-level content test.Compared with the existing ICP method,component analysis results of CA chips obtained by spectrophotometry are closer to real component content in samples and have satisfactory accuracy.Moreover,as its application in miniaturized explosive systems,the ignition ability of CA chips with different component contents for direct ink writing CL-20 and the corresponding mechanism was studied.This study provided a basis and idea for the design and performance evaluation of CA chips in miniaturized explosive systems.
基金supported by the State Grid Science&Technology Project(5100-202114296A-0-0-00).
文摘This article introduces the concept of load aggregation,which involves a comprehensive analysis of loads to acquire their external characteristics for the purpose of modeling and analyzing power systems.The online identification method is a computer-involved approach for data collection,processing,and system identification,commonly used for adaptive control and prediction.This paper proposes a method for dynamically aggregating large-scale adjustable loads to support high proportions of new energy integration,aiming to study the aggregation characteristics of regional large-scale adjustable loads using online identification techniques and feature extraction methods.The experiment selected 300 central air conditioners as the research subject and analyzed their regulation characteristics,economic efficiency,and comfort.The experimental results show that as the adjustment time of the air conditioner increases from 5 minutes to 35 minutes,the stable adjustment quantity during the adjustment period decreases from 28.46 to 3.57,indicating that air conditioning loads can be controlled over a long period and have better adjustment effects in the short term.Overall,the experimental results of this paper demonstrate that analyzing the aggregation characteristics of regional large-scale adjustable loads using online identification techniques and feature extraction algorithms is effective.
基金supported by China Postdoctoral Science Foundation(2019M651240)National Natural Science Foundation of China(31670559).
文摘As an important material for manufacturing resonant components of musical instruments,Paulownia has an important influence on the sound quality of Ruan.In this paper,a model for evaluating the sound quality of Ruan based on the vibration characteristics of wood is developed using machine learning methods.Generally,the selection of materials for Ruan manufacturing relies primarily on manually weighing,observing,striking,and listening by the instrument technician.Deficiencies in scientific theory have hindered the quality of the finished Ruan.In this study,nine Ruans were manufactured,and a prediction model of Ruan sound quality was proposed based on the raw material information of Ruans.Out of a total of 180 data sets,145 and 45 sets were chosen for training and validation,respec-tively.In this paper,typical correlation analysis was used to determine the correlation between two single indicators in two adjacent pairwise combinations of the measured objects in each stage of the production process in Ruan.The vibra-tion characteristics of the wood were tested,and a model for predicting the evaluation of Ruan’s acoustic qualities was developed by measuring the vibration characteristics of the resonating plate material.The acoustic quality of the Ruan sound board wood was evaluated and predicted using machine learning model generalized regression neural net-work.The results show that the prediction of Ruan sound quality can be achieved using Matlab simulation based on the vibration characteristics of the soundboard wood.When the model-predicted values were compared with the tradi-tional predicted results,it was found that the generalized regression neural network had good performance,achieving an accuracy of 93.8%which was highly consistent with the experimental results.It was concluded that the model can accurately predict the acoustic quality of the Ruan based on the vibration performance of the soundboards.