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.展开更多
Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the meta...Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.展开更多
The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical ...The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical properties of the resulting fibre to deviate from the designed properties. Therefore, a method of evaluating the optical properties of the actual fibre is necessary for the purpose of application. Up to now, the methods employed to measure the properties of the actual photonic crystal fibre often require long fibre samples or complex expensive equipments. To our knowledge, there are few studies of modeling an actual photonic crystal fibre and evaluating its properties rapidly. In this paper, a novel method, based on the combination model of digital image processing and the finite element method, is proposed to rapidly model the optical properties of the actual photonic crystal fibre. Two kinds of photonic crystal fibres made by Crystal Fiber A/S are modeled. It is confirmed from numerical results that the proposed method is simple, rapid and accurate for evaluating the optical properties of the actual photonic crystal fibre without requiring complex equipment.展开更多
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 two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictio...The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.展开更多
On the basis of a detailed discussion of the development of total ionizing dose （TID） effect model, a new commercial-model-independent TID modeling approach for partially depleted silicon-on-insulator metal-oxide- s...On the basis of a detailed discussion of the development of total ionizing dose （TID） effect model, a new commercial-model-independent TID modeling approach for partially depleted silicon-on-insulator metal-oxide- semiconductor field effect transistors is developed. An exponential approximation is proposed to simplify the trap charge calculation. Irradiation experiments with 60Co gamma rays for IO and core devices are performed to validate the simulation results. An excellent agreement of measurement with the simulation results is observed.展开更多
On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate an...On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.展开更多
In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite ...In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.展开更多
We proposed a practical way for mapping the results of coarse-grained molecular simulations to the observables in hydrogen change experiments.By combining an atomic-interaction based coarse-grained model with an all-a...We proposed a practical way for mapping the results of coarse-grained molecular simulations to the observables in hydrogen change experiments.By combining an atomic-interaction based coarse-grained model with an all-atom structure reconstruction algorithm,we reproduced the experimental hydrogen exchange data with reasonable accuracy using molecular dynamics simulations.We also showed that the coarse-grained model can be further improved by imposing experimental restraints from hydrogen exchange data via an iterative optimization strategy.These results suggest that it is feasible to develop an integrative molecular simulation scheme by incorporating the hydrogen exchange data into the coarse-grained molecular dynamics simulations and therefore help to overcome the accuracy bottleneck of coarse-grained models.展开更多
This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of com...This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.展开更多
In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the ...In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.展开更多
Currently,the cranes used at sea do not have enough flexibility,efficiency,and safety.Thus,this study proposed a floating multirobot coordinated towing system to meet the demands for offshore towing.Because of the fle...Currently,the cranes used at sea do not have enough flexibility,efficiency,and safety.Thus,this study proposed a floating multirobot coordinated towing system to meet the demands for offshore towing.Because of the flexibility of rope-driven robots,the one-way pulling characteristics of the rope,and the floating characteristics of the base,towing robots are easily overturned.First,the spatial configuration of the towing system was established according to the towing task,and the kinematic model of the towing system was established using the coordinate transformation.Then,the dynamic model of the towing system was established according to the rigid-body dynamics and hydrodynamic theory.Finally,the stability of the towing system was analyzed using the stability cone method.The simulation experiments provide a reference for the practical application of the floating multirobot coordinated towing system,which can improve the stability of towing systems by changing the configuration of the towing robot.展开更多
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.展开更多
Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an...Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.展开更多
Accurate acquisition of the rock stress is crucial for various rock engineering applications.The hollow inclusion (HI) technique is widely used for measuring in-situ rock stress.This technique calculates the stress te...Accurate acquisition of the rock stress is crucial for various rock engineering applications.The hollow inclusion (HI) technique is widely used for measuring in-situ rock stress.This technique calculates the stress tensor by measuring strain using an HI strain cell.However,existing analytical solutions for stress calculation based on an HI strain cell in a double-layer medium are not applicable when an HI strain cell is used in a three-layer medium,leading to erroneous stress calculations.To address this issue,this paper presents a method for calculating stress tensors in a three-layer medium using numerical simulations,specifically by obtaining a constitutive matrix that relates strain measurements to stress tensors in a three-layer medium.Furthermore,using Latin hypercube sampling (LHS) and orthogonal experimental design strategies,764 groups of numerical models encompassing various stress measurement scenarios have been established and calculated using FLAC^(3D)software.Finally,a surrogate model based on artificial neural network (ANN) was developed to predict constitutive matrices,achieving a goodness of fit (R^(2)) of 0.999 and a mean squared error (MSE) of 1.254.A software program has been developed from this surrogate model for ease of use in practical engineering applications.The method’s accuracy was verified through numerical simulations,analytical solution and laboratory experiment,demonstrating its effectiveness in calculating stress in a three-layer medium.The surrogate model was applied to calculate mining-induced stress in the roadway roof rock of a coal mine,a typical case for stress measurement in a three-layer medium.Errors in stress calculations arising from the use of existing analytical solutions were corrected.The study also highlights the significant errors associated with using double-layer analytical solutions in a three-layer medium,which could lead to inappropriate engineering design.展开更多
We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-or...We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.展开更多
We propose a fractional-order improved Fitz Hugh–Nagumo(FHN)neuron model in terms of a generalized Caputo fractional derivative.Following the existence of a unique solution for the proposed model,we derive the numeri...We propose a fractional-order improved Fitz Hugh–Nagumo(FHN)neuron model in terms of a generalized Caputo fractional derivative.Following the existence of a unique solution for the proposed model,we derive the numerical solution using a recently proposed L1 predictor–corrector method.The given method is based on the L1-type discretization algorithm and the spline interpolation scheme.We perform the error and stability analyses for the given method.We perform graphical simulations demonstrating that the proposed FHN neuron model generates rich electrical activities of periodic spiking patterns,chaotic patterns,and quasi-periodic patterns.The motivation behind proposing a fractional-order improved FHN neuron model is that such a system can provide a more nuanced description of the process with better understanding and simulation of the neuronal responses by incorporating memory effects and non-local dynamics,which are inherent to many biological systems.展开更多
An optical emission spectroscopy(OES)method with a non-invasive measurement capability,without inducing disturbance to the discharge,represents an effective method for material monitoring.However,when the OES method i...An optical emission spectroscopy(OES)method with a non-invasive measurement capability,without inducing disturbance to the discharge,represents an effective method for material monitoring.However,when the OES method is employed to monitor the trace erosion product within the ceramic channel of a Hall thruster,it becomes challenging to distinguish between signal and noise.In this study,we propose a model filtering method based on the signal characteristics of the Hall thruster plume spectrometer.This method integrates the slit imaging and spectral resolution features of the spectrometer.Employing this method,we extract the spectral signals of the erosion product and working gas from the Hall thruster under different operating conditions.The results indicate that our new method performs comparably to the traditional method without model filtering when extracting atom signals from strong xenon working gas.However,for trace amounts of the erosion product,our approach significantly enhances the signal-to-noise ratio(SNR),enabling the identification of extremely weak spectral signals even under low mass flow rate and low-voltage conditions.We obtain boron atom concentration of 3.91×10^(-3) kg/m^(3) at a mass flow rate of 4×10^(-7) kg/s and voltage of 200 V while monitoring a wider range of thruster operating conditions.The new method proposed in this study is suitable for monitoring other low-concentration elements,making it valuable for materials processing,environmental monitoring and space propulsion applications.展开更多
On the basis of introducing the fundamental theory and the basic analysis steps of the sub model method, the strength of the new engine complex assembly structure was analyzed according to the properties of the engin...On the basis of introducing the fundamental theory and the basic analysis steps of the sub model method, the strength of the new engine complex assembly structure was analyzed according to the properties of the engine structures, some of the key parts of the engine were analyzed with refined mesh by sub model method and the error of the FEM solution was estimated by the extrapolation method. The example showed that the sub model can not only analyze the comlex structures without the restriction of the software and hardware of the computers, but get the more precise analysis result also. This method is more suitable for the strength analysis of the complex assembly structure.展开更多
文摘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.
基金supported jointly by the National Basic Research Program of China("973"Program)(No2014CB046200)the National Science Foundation of Jiangsu Province(No.BK2014059)+1 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe National Natural Science Foundation of China(No.11172135)
文摘Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.
基金Project supported by the National Basic Research Program of China(Grant No.2010CB328206)the National Natural Science Foundation of China(Grant No.60977033)the Science and Technology Innovation Foundation for Excellent Doctors of Beijing Jiaotong University,China(Grant Nos.141055522 and 141060522)
文摘The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical properties of the resulting fibre to deviate from the designed properties. Therefore, a method of evaluating the optical properties of the actual fibre is necessary for the purpose of application. Up to now, the methods employed to measure the properties of the actual photonic crystal fibre often require long fibre samples or complex expensive equipments. To our knowledge, there are few studies of modeling an actual photonic crystal fibre and evaluating its properties rapidly. In this paper, a novel method, based on the combination model of digital image processing and the finite element method, is proposed to rapidly model the optical properties of the actual photonic crystal fibre. Two kinds of photonic crystal fibres made by Crystal Fiber A/S are modeled. It is confirmed from numerical results that the proposed method is simple, rapid and accurate for evaluating the optical properties of the actual photonic crystal fibre without requiring complex equipment.
基金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.
基金Project supported by the Hefei National Research Center for Physical Sciences at the Microscale (Grant No.KF2021002)the Natural Science Foundation of Shanxi Province,China (Grant Nos.202303021221029 and 202103021224051)+2 种基金the National Natural Science Foundation of China (Grant Nos.11975024,12047503,and 12275263)the Anhui Provincial Supporting Program for Excellent Young Talents in Colleges and Universities (Grant No.gxyq ZD2019023)the National Key Research and Development Program of China (Grant No.2018YFA0306501)。
文摘The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.
基金Supported by the National Natural Science Foundation of China under Grant Nos 61404151 and 61574153
文摘On the basis of a detailed discussion of the development of total ionizing dose （TID） effect model, a new commercial-model-independent TID modeling approach for partially depleted silicon-on-insulator metal-oxide- semiconductor field effect transistors is developed. An exponential approximation is proposed to simplify the trap charge calculation. Irradiation experiments with 60Co gamma rays for IO and core devices are performed to validate the simulation results. An excellent agreement of measurement with the simulation results is observed.
基金supported by the Key Projects of the National Science Foundation of China(Nos.52178190,52078250,11832013)
文摘On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.
基金the National Natural Science Foundation of China (5963140) Doctor Point Fund of National Education Committee Parent Company Fund of Aviation Industry
文摘In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.
基金the National Natural Science Foundation of China(Grant Nos.11974173 and 11934008)the HPC Center of Nanjing University。
文摘We proposed a practical way for mapping the results of coarse-grained molecular simulations to the observables in hydrogen change experiments.By combining an atomic-interaction based coarse-grained model with an all-atom structure reconstruction algorithm,we reproduced the experimental hydrogen exchange data with reasonable accuracy using molecular dynamics simulations.We also showed that the coarse-grained model can be further improved by imposing experimental restraints from hydrogen exchange data via an iterative optimization strategy.These results suggest that it is feasible to develop an integrative molecular simulation scheme by incorporating the hydrogen exchange data into the coarse-grained molecular dynamics simulations and therefore help to overcome the accuracy bottleneck of coarse-grained models.
基金This work was supported by China Postdoctoral Science Foundation(No.2022M723391)the Science and Technology Innovation Project of Higher Education in Shanxi Province(No.2019L0754)+1 种基金the Central Guiding Local Science and Technology Development Fund Project(No.YDZJSX2021B021)Shanxi Province Basic Research Plan General Project(No.202203021221294).
文摘This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.
基金Supported by the Beijing Municipal Science&Technology Commission(Z211100004421012),the Key Reaserch and Development Pro⁃gram of China(2022YFF0605902)。
文摘In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.
基金Supported by the National Natural Science Foundation of China under Grant No.51965032the Natural Science Foundation of Gansu Province of China under Grant No.22JR5RA319+2 种基金the Excellent Doctoral Student Foundation of Gansu Province of China under Grant No.23JRRA842the Sichuan Province Engineering Technology Research Center of General Aircraft Maintenance under Grant No.GAMRC2023YB05the Key Research and Development Project of Lanzhou Jiaotong University under Grant No.LZJTUZDYF2302.
文摘Currently,the cranes used at sea do not have enough flexibility,efficiency,and safety.Thus,this study proposed a floating multirobot coordinated towing system to meet the demands for offshore towing.Because of the flexibility of rope-driven robots,the one-way pulling characteristics of the rope,and the floating characteristics of the base,towing robots are easily overturned.First,the spatial configuration of the towing system was established according to the towing task,and the kinematic model of the towing system was established using the coordinate transformation.Then,the dynamic model of the towing system was established according to the rigid-body dynamics and hydrodynamic theory.Finally,the stability of the towing system was analyzed using the stability cone method.The simulation experiments provide a reference for the practical application of the floating multirobot coordinated towing system,which can improve the stability of towing systems by changing the configuration of the towing robot.
基金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.
基金Supported by the National Natural Science Foundation of China (Grant Nos.52088102 and 51879287)National Key Research and Development Program of China (Grant No.2022YFB2602301)。
文摘Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.
基金funding support from the National Natural Science Foundation of China (Nos. 42477208 and 52079134)the Natural Science Foundation of Hubei Province, China (No. 2024AFA072)+2 种基金the Youth Innovation Promotion Association CAS (No. 2022332)the National Key R&D Program of China (No. 2024YFF0508203)the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering Safety (Nos. SKLGME-JBGS2402 and SKLGME022022)。
文摘Accurate acquisition of the rock stress is crucial for various rock engineering applications.The hollow inclusion (HI) technique is widely used for measuring in-situ rock stress.This technique calculates the stress tensor by measuring strain using an HI strain cell.However,existing analytical solutions for stress calculation based on an HI strain cell in a double-layer medium are not applicable when an HI strain cell is used in a three-layer medium,leading to erroneous stress calculations.To address this issue,this paper presents a method for calculating stress tensors in a three-layer medium using numerical simulations,specifically by obtaining a constitutive matrix that relates strain measurements to stress tensors in a three-layer medium.Furthermore,using Latin hypercube sampling (LHS) and orthogonal experimental design strategies,764 groups of numerical models encompassing various stress measurement scenarios have been established and calculated using FLAC^(3D)software.Finally,a surrogate model based on artificial neural network (ANN) was developed to predict constitutive matrices,achieving a goodness of fit (R^(2)) of 0.999 and a mean squared error (MSE) of 1.254.A software program has been developed from this surrogate model for ease of use in practical engineering applications.The method’s accuracy was verified through numerical simulations,analytical solution and laboratory experiment,demonstrating its effectiveness in calculating stress in a three-layer medium.The surrogate model was applied to calculate mining-induced stress in the roadway roof rock of a coal mine,a typical case for stress measurement in a three-layer medium.Errors in stress calculations arising from the use of existing analytical solutions were corrected.The study also highlights the significant errors associated with using double-layer analytical solutions in a three-layer medium,which could lead to inappropriate engineering design.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11571366 and 11501570the Open Foundation of State Key Laboratory of High Performance Computing of China+1 种基金the Research Fund of National University of Defense Technology under Grant No JC15-02-02the Fund from HPCL
文摘We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.
文摘We propose a fractional-order improved Fitz Hugh–Nagumo(FHN)neuron model in terms of a generalized Caputo fractional derivative.Following the existence of a unique solution for the proposed model,we derive the numerical solution using a recently proposed L1 predictor–corrector method.The given method is based on the L1-type discretization algorithm and the spline interpolation scheme.We perform the error and stability analyses for the given method.We perform graphical simulations demonstrating that the proposed FHN neuron model generates rich electrical activities of periodic spiking patterns,chaotic patterns,and quasi-periodic patterns.The motivation behind proposing a fractional-order improved FHN neuron model is that such a system can provide a more nuanced description of the process with better understanding and simulation of the neuronal responses by incorporating memory effects and non-local dynamics,which are inherent to many biological systems.
基金financially supported by National Natural Science Foundation of China(No.U22B2094)。
文摘An optical emission spectroscopy(OES)method with a non-invasive measurement capability,without inducing disturbance to the discharge,represents an effective method for material monitoring.However,when the OES method is employed to monitor the trace erosion product within the ceramic channel of a Hall thruster,it becomes challenging to distinguish between signal and noise.In this study,we propose a model filtering method based on the signal characteristics of the Hall thruster plume spectrometer.This method integrates the slit imaging and spectral resolution features of the spectrometer.Employing this method,we extract the spectral signals of the erosion product and working gas from the Hall thruster under different operating conditions.The results indicate that our new method performs comparably to the traditional method without model filtering when extracting atom signals from strong xenon working gas.However,for trace amounts of the erosion product,our approach significantly enhances the signal-to-noise ratio(SNR),enabling the identification of extremely weak spectral signals even under low mass flow rate and low-voltage conditions.We obtain boron atom concentration of 3.91×10^(-3) kg/m^(3) at a mass flow rate of 4×10^(-7) kg/s and voltage of 200 V while monitoring a wider range of thruster operating conditions.The new method proposed in this study is suitable for monitoring other low-concentration elements,making it valuable for materials processing,environmental monitoring and space propulsion applications.
文摘On the basis of introducing the fundamental theory and the basic analysis steps of the sub model method, the strength of the new engine complex assembly structure was analyzed according to the properties of the engine structures, some of the key parts of the engine were analyzed with refined mesh by sub model method and the error of the FEM solution was estimated by the extrapolation method. The example showed that the sub model can not only analyze the comlex structures without the restriction of the software and hardware of the computers, but get the more precise analysis result also. This method is more suitable for the strength analysis of the complex assembly structure.
