We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground...We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.展开更多
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.展开更多
In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian co...In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.展开更多
In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discreti...In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.展开更多
A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous H...A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.展开更多
For classical Hamiltonian with general form H = 1/2∑ijMijpipj＋1/2∑ijLijqiqj we find a new convenient way to obtain its normal coordinates, namely, let H be quantised and then employ the invariant eigen-operator （...For classical Hamiltonian with general form H = 1/2∑ijMijpipj＋1/2∑ijLijqiqj we find a new convenient way to obtain its normal coordinates, namely, let H be quantised and then employ the invariant eigen-operator （IEO） method （Fan et al. 2004 Phys. Lett. A 321 75） to derive them. The general matrix equation, which relies on M and L, for obtaining the normal coordinates of H is derived.展开更多
Eigenvalue-solution to those Hamiltonians involving non-commutative coordinates is not easily obtained. In this paper we apply the invariant eigen-operator （IEO） method to solving the energy spectrmn of the three-mo...Eigenvalue-solution to those Hamiltonians involving non-commutative coordinates is not easily obtained. In this paper we apply the invariant eigen-operator （IEO） method to solving the energy spectrmn of the three-mode harmonic oscillator in non-commutative space with the coordinate operators satisfying cyclic commutative relations, [X1, X2] = [X2, X3]=[X3, X1] = iθ, and this method seems effective and concise.展开更多
Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved ...Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.展开更多
Based on a total of 16 indicators selected from the tourism and transport industries,an evaluation index system of the coupling and coordination development level of tourism and transport is constructed.The entropy va...Based on a total of 16 indicators selected from the tourism and transport industries,an evaluation index system of the coupling and coordination development level of tourism and transport is constructed.The entropy value method and the coupling coordination degree model are used to conduct an empirical study on the development level and coupling coordination level of the transport and tourism industries in Chengdu City from 2011 to 2020.The results show that,on the whole,the coupling coordination degree of transport and tourism in Chengdu is poor and has been in a state of mild to moderate dysfunction.The development level of tourism lagged behind the development of transport from 2011 to 2012,and the two were in a state of mild dysfunction.However,from 2013 onwards,the development level of tourism was prioritized over the development level of transport.This shift caused the coupling coordination degree of the two industries to decline sharply to 0.23305 in 2013.The development level of the tourism industry increased again,reaching 0.34206 in 2019,which marked an improvement.Consequently,the coordination degree of the transport and tourism industries evolved from moderate dislocation to mild dislocation.Finally,the results of the empirical research are analyzed,and corresponding suggestions are put forward to promote the sustainable growth of the transport and tourism industries in Chengdu City.These suggestions aim to improve the coupled and coordinated development level of the two industries.展开更多
Bond valence method illustrates the relation between valence and length of a particular bond type. This theory has been used to predict structure information, but the effect is very limited. In this paper, two indexes...Bond valence method illustrates the relation between valence and length of a particular bond type. This theory has been used to predict structure information, but the effect is very limited. In this paper, two indexes, i.e., global instability index（GII） and bond strain index（BSI）, are adopted as a judgment of a search-match program for prediction. The results show that with GII and BSI combined as judgment, the predicted atom positions are very close to real ones. The mechanism and validity of this searching program are also discussed. The GII ＆ BSI distribution contour map reveals that the predicted function is a reflection of exponential feature of bond valence formula. This combined searching method may be integrated with other structure-determination method, and may be helpful in refining and testifying light atom positions.展开更多
A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of part...A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of partial differential equations defined in a 2D or 3D geological model. The new approach refers to a “strong version” against the “weak version” of the subspace spectral method based on the variational principle or Galerkin’s weighting scheme. We incorporate local nonlinear transformations and global spline interpolations in a curved coordinate system and make the discrete grid exactly matches geometry of the model so that it is achieved to convert the global domain into subdomains and apply Chebyshev points to locally sampling physical quantities and globally computing the spatial derivatives. This new approach not only remains exponential convergence of the standard spectral method in subdomains, but also yields a sparse assembled matrix when applied for the global domain simulations. We conducted 2D and 3D synthetic experiments and compared accuracies of the numerical differentiations with traditional finite difference approaches. The results show that as the points of differentiation vector are larger than five, the subdomain Chebyshev spectral method significantly improve the accuracies of the finite difference approaches.展开更多
This paper presents an array pattern synthesis algorithm for arbitrary arrays based on coordinate descent method (CDM). With this algorithm, the complex element weights are found to minimize a weighted L2 norm of the ...This paper presents an array pattern synthesis algorithm for arbitrary arrays based on coordinate descent method (CDM). With this algorithm, the complex element weights are found to minimize a weighted L2 norm of the difference between desired and achieved pattern. Compared with traditional optimization techniques, CDM is easy to implement and efficient to reach the optimum solutions. Main advantage is the flexibility. CDM is suitable for linear and planar array with arbitrary array elements on arbitrary positions. With this method, we can configure arbitrary beam pattern, which gives it the ability to solve variety of beam forming problem, e.g. focused beam, shaped beam, nulls at arbitrary direction and with arbitrary beam width. CDM is applicable for phase-only and amplitude-only arrays as well, and furthermore, it is a suitable method to treat the problem of array with element failures.展开更多
In atomic,molecular,and nuclear physics,the method of complex coordinate rotation is a widely used theoretical tool for studying resonant states.Here,we propose a novel implementation of this method based on the gradi...In atomic,molecular,and nuclear physics,the method of complex coordinate rotation is a widely used theoretical tool for studying resonant states.Here,we propose a novel implementation of this method based on the gradient optimization(CCR-GO).The main strength of the CCR-GO method is that it does not require manual adjustment of optimization parameters in the wave function;instead,a mathematically well-defined optimization path can be followed.Our method is proven to be very efficient in searching resonant positions and widths over a variety of few-body atomic systems,and can significantly improve the accuracy of the results.As a special case,the CCR-GO method is equally capable of dealing with bound-state problems with high accuracy,which is traditionally achieved through the usual extreme conditions of energy itself.展开更多
Dirac's method which itself is for constrained Boson fields and particle systems is followed and developed to treat Dirac fields in light-front coordinates.
In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are reg...In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are regulated through sensitivity filtering tomitigate numerical instabilities associatedwith stress concentrations.Ap-norm aggregation function is employed to globalize local stress constraints,and a normalization technique linearly weights strain energy and stress,transforming the multi-objective problem into a single-objective formulation.The sensitivity of the objective function with respect to design variables is rigorously derived.Three numerical examples are presented,comparing the optimized structures in terms of strain energy,mass,and stress across five different mathematical models with varying combinations of optimization objectives.The results validate the effectiveness and feasibility of the proposed method for achieving a balanced design between structural stiffness and strength.This approach offers a new perspective for future research on stiffness-strength coordinated structural optimization.展开更多
The precise acquisition of the quality characteristic parameters of large aircraft directly affects its performance characteristics.For large aircrafts such as missiles and rockets with internal fillings,traditional m...The precise acquisition of the quality characteristic parameters of large aircraft directly affects its performance characteristics.For large aircrafts such as missiles and rockets with internal fillings,traditional measurement methods involving large-angle tilting or rotation may pose safety risks.In light of the characteristics of large aircraft and in combination with existing measurement methods,we design a mass and centroid measurement method based on four-point support and small-angle tilting,and develop a set of mass and centroid testing system.This method obtains the intersection point of the gravity action line in the product coordinate system through coordinate transformation in two postures,thereby obtaining the three-dimensional centroid of the aircraft.We first elaborate on the principle of this method in detail,then introduce the composition of the equipment,and analyze the structural stress of key components.Finally,experimental verification and uncertainty analysis are carried out.Experimental verification shows that the maximum deviation of the mass measurement accuracy is less than 0.02%,the centroid measurement accuracy in the X direction is±0.15 mm,in the Y direction it is±0.21 mm,and in the Z direction it is±0.19 mm.展开更多
The analysis and calculating method of dynamic errors of CMMs during probing are discussed.To relate the dynamic displacement errors with the dynamic rotational errors a method for obtaining the displacement errors at...The analysis and calculating method of dynamic errors of CMMs during probing are discussed.To relate the dynamic displacement errors with the dynamic rotational errors a method for obtaining the displacement errors at the probing position from dynamic rotational errors is presented.It is pointed out that the finite element method might be used for modeling dynamic errors.However,dynamic errors are difficult to be modeled so a combined practical and theoretical approach is needed.In addition,the dynamic errors are measured with inductive position sensors.展开更多
A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution f...A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier–Stokes equations, including contribution of chemical reaction, via the Chapman–Enskog expansion.For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.展开更多
Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water （CDW）, a new hybrid vertical coordinate is designed, which uses σ co...Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water （CDW）, a new hybrid vertical coordinate is designed, which uses σ coordinate for current and σ-z coordinate for salinity. To combine the current and salinity, the Eulerian-Lagrangian method is used for the salinity calculation, and the baroclinic pressure gradient （BPG） is calculated on the salinity sited layers. The new hybrid vertical coordinate is introduced to the widely used model of POM （Princeton Ocean Model） to make a new model of POM-σ-z. The BPG calculations of an ideal case show that POM-σ-z model brings smaller error than POM model does. The simulations of CDW also show that POM-σ-z model is better than POM model on simulating the salinity and its front.展开更多
基金supported by the National Key R&D Program of China(No.2023YFA1606701)the National Natural Science Foundation of China(Nos.12175042,11890710,11890714,12047514,12147101,and 12347106)+1 种基金Guangdong Major Project of Basic and Applied Basic Research(No.2020B0301030008)China National Key R&D Program(No.2022YFA1602402).
文摘We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11261035 and 11171038)the Science Research Foundation of the Institute of Higher Education of Inner Mongolia Autonomous Region, China (Grant No. NJZZ12198)the Natural Science Foundation of Inner Mongolia Autonomous Region, China (Grant No. 2012MS0102)
文摘In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11261035,11171038,and 10771019)the Science Reaearch Foundation of Institute of Higher Education of Inner Mongolia Autonomous Region,China (Grant No. NJZZ12198)the Natural Science Foundation of Inner Mongolia Autonomous Region,China (Grant No. 2012MS0102)
文摘In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.
文摘A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.
基金supported by the National Natural Science Foundation of China (Grant No.10874174)the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No.20070358009)
文摘For classical Hamiltonian with general form H = 1/2∑ijMijpipj＋1/2∑ijLijqiqj we find a new convenient way to obtain its normal coordinates, namely, let H be quantised and then employ the invariant eigen-operator （IEO） method （Fan et al. 2004 Phys. Lett. A 321 75） to derive them. The general matrix equation, which relies on M and L, for obtaining the normal coordinates of H is derived.
基金the President Foundation of the Chinese Academy of Sciencesthe Specialized Research Fund for the Doctoral Program of Higher Education
文摘Eigenvalue-solution to those Hamiltonians involving non-commutative coordinates is not easily obtained. In this paper we apply the invariant eigen-operator （IEO） method to solving the energy spectrmn of the three-mode harmonic oscillator in non-commutative space with the coordinate operators satisfying cyclic commutative relations, [X1, X2] = [X2, X3]=[X3, X1] = iθ, and this method seems effective and concise.
基金support of this work by the National Natural Science Foundation of China(No.51405096)the Fundamental Research Funds for the Central Universities(HEUCF210710).
文摘Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.
基金Chongqing University of Science and Technology Postgraduate Innovation Program Project(Project No.YKJCX2320902)。
文摘Based on a total of 16 indicators selected from the tourism and transport industries,an evaluation index system of the coupling and coordination development level of tourism and transport is constructed.The entropy value method and the coupling coordination degree model are used to conduct an empirical study on the development level and coupling coordination level of the transport and tourism industries in Chengdu City from 2011 to 2020.The results show that,on the whole,the coupling coordination degree of transport and tourism in Chengdu is poor and has been in a state of mild to moderate dysfunction.The development level of tourism lagged behind the development of transport from 2011 to 2012,and the two were in a state of mild dysfunction.However,from 2013 onwards,the development level of tourism was prioritized over the development level of transport.This shift caused the coupling coordination degree of the two industries to decline sharply to 0.23305 in 2013.The development level of the tourism industry increased again,reaching 0.34206 in 2019,which marked an improvement.Consequently,the coordination degree of the transport and tourism industries evolved from moderate dislocation to mild dislocation.Finally,the results of the empirical research are analyzed,and corresponding suggestions are put forward to promote the sustainable growth of the transport and tourism industries in Chengdu City.These suggestions aim to improve the coupled and coordinated development level of the two industries.
基金Project supported by the National Natural Science Foundation of China(Grant No.51272027)
文摘Bond valence method illustrates the relation between valence and length of a particular bond type. This theory has been used to predict structure information, but the effect is very limited. In this paper, two indexes, i.e., global instability index（GII） and bond strain index（BSI）, are adopted as a judgment of a search-match program for prediction. The results show that with GII and BSI combined as judgment, the predicted atom positions are very close to real ones. The mechanism and validity of this searching program are also discussed. The GII ＆ BSI distribution contour map reveals that the predicted function is a reflection of exponential feature of bond valence formula. This combined searching method may be integrated with other structure-determination method, and may be helpful in refining and testifying light atom positions.
文摘A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of partial differential equations defined in a 2D or 3D geological model. The new approach refers to a “strong version” against the “weak version” of the subspace spectral method based on the variational principle or Galerkin’s weighting scheme. We incorporate local nonlinear transformations and global spline interpolations in a curved coordinate system and make the discrete grid exactly matches geometry of the model so that it is achieved to convert the global domain into subdomains and apply Chebyshev points to locally sampling physical quantities and globally computing the spatial derivatives. This new approach not only remains exponential convergence of the standard spectral method in subdomains, but also yields a sparse assembled matrix when applied for the global domain simulations. We conducted 2D and 3D synthetic experiments and compared accuracies of the numerical differentiations with traditional finite difference approaches. The results show that as the points of differentiation vector are larger than five, the subdomain Chebyshev spectral method significantly improve the accuracies of the finite difference approaches.
文摘This paper presents an array pattern synthesis algorithm for arbitrary arrays based on coordinate descent method (CDM). With this algorithm, the complex element weights are found to minimize a weighted L2 norm of the difference between desired and achieved pattern. Compared with traditional optimization techniques, CDM is easy to implement and efficient to reach the optimum solutions. Main advantage is the flexibility. CDM is suitable for linear and planar array with arbitrary array elements on arbitrary positions. With this method, we can configure arbitrary beam pattern, which gives it the ability to solve variety of beam forming problem, e.g. focused beam, shaped beam, nulls at arbitrary direction and with arbitrary beam width. CDM is applicable for phase-only and amplitude-only arrays as well, and furthermore, it is a suitable method to treat the problem of array with element failures.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.91636216,11974382,and 11474316)the Chinese Academy of Sciences Strategic Priority Research Program(Grant No.XDB21020200)+1 种基金by the YIPA Programthe support of NSERC,SHARCnet,ACEnet of Canada。
文摘In atomic,molecular,and nuclear physics,the method of complex coordinate rotation is a widely used theoretical tool for studying resonant states.Here,we propose a novel implementation of this method based on the gradient optimization(CCR-GO).The main strength of the CCR-GO method is that it does not require manual adjustment of optimization parameters in the wave function;instead,a mathematically well-defined optimization path can be followed.Our method is proven to be very efficient in searching resonant positions and widths over a variety of few-body atomic systems,and can significantly improve the accuracy of the results.As a special case,the CCR-GO method is equally capable of dealing with bound-state problems with high accuracy,which is traditionally achieved through the usual extreme conditions of energy itself.
文摘Dirac's method which itself is for constrained Boson fields and particle systems is followed and developed to treat Dirac fields in light-front coordinates.
基金funded by National Nature Science Foundation of China(92266203)National Nature Science Foundation of China(52205278)+1 种基金Key Projects of Shijiazhuang Basic Research Program(241791077A)Central Guide Local Science and Technology Development Fund Project of Hebei Province(246Z1022G).
文摘In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are regulated through sensitivity filtering tomitigate numerical instabilities associatedwith stress concentrations.Ap-norm aggregation function is employed to globalize local stress constraints,and a normalization technique linearly weights strain energy and stress,transforming the multi-objective problem into a single-objective formulation.The sensitivity of the objective function with respect to design variables is rigorously derived.Three numerical examples are presented,comparing the optimized structures in terms of strain energy,mass,and stress across five different mathematical models with varying combinations of optimization objectives.The results validate the effectiveness and feasibility of the proposed method for achieving a balanced design between structural stiffness and strength.This approach offers a new perspective for future research on stiffness-strength coordinated structural optimization.
基金supported by National Natural Science Foundation of China-Youth Program(No.62303420)。
文摘The precise acquisition of the quality characteristic parameters of large aircraft directly affects its performance characteristics.For large aircrafts such as missiles and rockets with internal fillings,traditional measurement methods involving large-angle tilting or rotation may pose safety risks.In light of the characteristics of large aircraft and in combination with existing measurement methods,we design a mass and centroid measurement method based on four-point support and small-angle tilting,and develop a set of mass and centroid testing system.This method obtains the intersection point of the gravity action line in the product coordinate system through coordinate transformation in two postures,thereby obtaining the three-dimensional centroid of the aircraft.We first elaborate on the principle of this method in detail,then introduce the composition of the equipment,and analyze the structural stress of key components.Finally,experimental verification and uncertainty analysis are carried out.Experimental verification shows that the maximum deviation of the mass measurement accuracy is less than 0.02%,the centroid measurement accuracy in the X direction is±0.15 mm,in the Y direction it is±0.21 mm,and in the Z direction it is±0.19 mm.
文摘The analysis and calculating method of dynamic errors of CMMs during probing are discussed.To relate the dynamic displacement errors with the dynamic rotational errors a method for obtaining the displacement errors at the probing position from dynamic rotational errors is presented.It is pointed out that the finite element method might be used for modeling dynamic errors.However,dynamic errors are difficult to be modeled so a combined practical and theoretical approach is needed.In addition,the dynamic errors are measured with inductive position sensors.
基金support of the Science Foundations of National Key Laboratory of Computational Physics,National Natural Science Foundation of China under Grant No.11202003the opening project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology) under Grant No.KFJJ14-1M.YL +1 种基金support of National Natural Science Foundation of China under Grant Nos.11074300 and 41472130National Basic Research Program of China under Grant No.2013CBA01504
文摘A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier–Stokes equations, including contribution of chemical reaction, via the Chapman–Enskog expansion.For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.
文摘Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water （CDW）, a new hybrid vertical coordinate is designed, which uses σ coordinate for current and σ-z coordinate for salinity. To combine the current and salinity, the Eulerian-Lagrangian method is used for the salinity calculation, and the baroclinic pressure gradient （BPG） is calculated on the salinity sited layers. The new hybrid vertical coordinate is introduced to the widely used model of POM （Princeton Ocean Model） to make a new model of POM-σ-z. The BPG calculations of an ideal case show that POM-σ-z model brings smaller error than POM model does. The simulations of CDW also show that POM-σ-z model is better than POM model on simulating the salinity and its front.
