Quantifying the post-earthquake functional recovery of railway stations presents significant challenges.This paper first establishes a post-earthquake function calculation method for railway stations,encompassing the ...Quantifying the post-earthquake functional recovery of railway stations presents significant challenges.This paper first establishes a post-earthquake function calculation method for railway stations,encompassing the establishment of relationships between the station’s function and the damage state,function loss,and failure probability of components and professional equipment in each layer.Also,the“4 stages-6 sequences”post-earthquake repair method is present,taking into account the functional and structural characteristics of railway stations.Additionally,a novel piecewise function for the post-earthquake functional dynamic recovery of railway stations is developed.A case study is conducted on a typical railway station to demonstrate the analysis procedure.Results indicate that under fortification,rare,and extremely rare earthquake scenarios,the interlayer drift ratio(IDR)of the railway station were 1/276,1/143,and 1/52,respectively,and corresponding peak floor acceleration(PFA)were 6.31 m/s^(2),7.82 m/s^(2),and 8.57 m/s^(2),respectively.The post-earthquake function of the railway station was 93.21%,82.33%,and 64.16%of its initial function.The repair times were 6.66 days,18.65 days,and 37.42 days.The displacement-sensitive,non-structural components were identified as the most vulnerable to damage.And the first repair stage(R_(1))which was mainly used to repair structural components and non-structural transport components,accounted for the highest proportion of total repair time.展开更多
Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of thes...Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.展开更多
This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the acc...This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.展开更多
A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the re...A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the reference functions method to describe the solution to the depth-separated wave equation approximately using parabolic cylinder functions.The conditions for the validity of this approximation are also discussed.Furthermore,a formula that incorporates waveguide effects for the modal group velocity is derived,revealing that boundary effects at very low frequencies can have a significant impact on the propagation characteristics of even low-order normal modes.The present method not only offers improved accuracy compared to the classical WKB approximation and the uniform asymptotic approximation based on Airy functions,but also provides a wider range of depth applicability.Additionally,this method exhibits strong agreement with numerical methods and offers valuable physical insights.Finally,the method is applied to the study of very-low-frequency sound propagation in the South China Sea,leading to sound transmission loss predictions that closely align with experimental observations.展开更多
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.展开更多
Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instr...Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instrument function,is measured in X-ray Thom-son scattering(XRTS)experiments,which allow the study of electronic structure properties at the microscopic level.Among the various ab initio methods,linear-response time-dependent density-functional theory(LR-TDDFT)is a key framework for simulating the DSF.The standard approach in LR-TDDFT for computing the DSF relies on the orbital representation.A significant drawback of this method is the unfavorable scaling of the number of required empty bands as the wavenumber increases,making LR-TDDFT impractical for modeling XRTS measurements over large energy scales,such as in backward scattering geometry.In this work,we consider and test an alternative approach to LR-TDDFT that employs the Liouville–Lanczos(LL)method for simulating the DSF of WDM.This approach does not require empty states and allows the DSF at large momentum transfer values and over a broad frequency range to be accessed.We compare the results obtained from the LL method with those from the solution of Dyson’s equation using the standard LR-TDDFT within the projector augmented-wave formalism for isochorically heated aluminum and warm dense hydrogen.Additionally,we utilize exact path integral Monte Carlo results for the imaginary-time density-density correlation function(ITCF)of warm dense hydrogen to rigorously benchmark the LL approach.We discuss the application of the LL method for calculating DSFs and ITCFs at different wavenumbers,the effects of pseudopotentials,and the role of Lorentzian smearing.The successful validation of the LL method under WDM conditions makes it a valuable addition to the ab initio simulation landscape,supporting experimental efforts and advancing WDM theory.展开更多
Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct dou...Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct doubly-periodic solutions of the Zakharov-Kuznetsov equation, which has been derived by Gottwald as a two-dimensional model for nonlinear Rossby waves. When the modulus k →1, these solutions reduce to the solitary wave solutions of the equation.展开更多
We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projec...We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.展开更多
This paper elucidates the effectiveness of combining the Poincare-Lighthill-Kuo method (PLK method, for short) and symbolic computation. Firstly, the idea and history of the PLK method are briefly introduced. Then, th...This paper elucidates the effectiveness of combining the Poincare-Lighthill-Kuo method (PLK method, for short) and symbolic computation. Firstly, the idea and history of the PLK method are briefly introduced. Then, the difficulty of intermediate expression swell, often encountered in symbolic computation, is outlined. For overcoming the difficulty, a semi-inverse algorithm was proposed by the author, with which the lengthy ports of intermediate expressions are first frozen in the form of symbols till the Fnal stage of seeking perturbation solutions. Tn discuss the applications of the above algorithm, the related work of the author and his research group on nonlinear oscillations and waves is concisely reviewed. The computer-extended perturbation solution of the Duffing equation shows that the asymptotic solution obtained with the PLK method possesses the convergence radius of 1 and thus the range of validity of the solution is considerably enlarged. The studies on internal solitary waves in stratified fluid and on the head-on collision between two solitary waves in a hyperelastic rod indicate that by means of the presented methods, very complicated manipulation, unconceivable in hand calculation, can be conducted and thus result in higher-order evolution equations and asymptotic solutions. The examples illustrate that the algorithm helps to realize the symbolic computation on micro-commputers. Finally, it is concluded that,vith the aid of symbolic computation, the vitality of the PLK method is greatly. Strengthened and at least for the solutions to conservative systems of oscillations and waves, it is a powerful tool.展开更多
This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node...This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.展开更多
The objective in this presentation is to introduce some of the unique properties and applications of nullors in active circuit analysis and designs. The emphasis is to discuss the role nullors can play in symbolic rep...The objective in this presentation is to introduce some of the unique properties and applications of nullors in active circuit analysis and designs. The emphasis is to discuss the role nullors can play in symbolic representation of transfer functions. To show this we adopt the topological platform for the circuit analysis and use a recently developed Admittance Method (AM) to achieve the Sum of Tree Products (STP), replacing the determinant and cofactors of the Nodal Admittance Matrix (NAM) of the circuit. To construct a transfer function, we start with a given active circuit and convert all its controlled sources and I/O-ports to nullors. Now, with a solid nullor circuit (passive elements and nullors) we first eliminate the passive elements through AM operations. This produces the STPs. Second, the all-nullor circuit is then used to find the signs or the STPs. Finally, the transfer function (in symbolic, if chosen) is obtained from the ratio between the STPs.展开更多
With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying micr...With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.展开更多
Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore th...Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.展开更多
In this paper, by using symbolic and algebra computation, Chen and Wang＇s multiple R/ccati equations rational expansion method was further extended. Many double soliton-like and other novel combined forms of exact so...In this paper, by using symbolic and algebra computation, Chen and Wang＇s multiple R/ccati equations rational expansion method was further extended. Many double soliton-like and other novel combined forms of exact solutions of the （2＋1）-dimensional Breaking soliton equation are derived by using the extended multiple Riccatl equations expansion method.展开更多
To quantify the seismic resilience of buildings,a method for evaluating functional loss from the component level to the overall building is proposed,and the dual-parameter seismic resilience assessment method based on...To quantify the seismic resilience of buildings,a method for evaluating functional loss from the component level to the overall building is proposed,and the dual-parameter seismic resilience assessment method based on postearthquake loss and recovery time is improved.A threelevel function tree model is established,which can consider the dynamic changes in weight coefficients of different category of components relative to their functional losses.Bayesian networks are utilized to quantify the impact of weather conditions,construction technology levels,and worker skill levels on component repair time.A method for determining the real-time functional recovery curve of buildings based on the component repair process is proposed.Taking a three-story teaching building as an example,the seismic resilience indices under basic earthquakes and rare earthquakes are calculated.The results show that the seismic resilience grade of the teaching building is comprehensively judged as GradeⅢ,and its resilience grade is more significantly affected by postearthquake loss.The proposed method can be used to predict the seismic resilience of buildings prior to earthquakes,identify weak components within buildings,and provide guidance for taking measures to enhance the seismic resilience of buildings.展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
The effects of plasma screening on the ^(1)P^(o) resonance states of H-and He below the n=3 and n=4 thresholds of the respective subsystemsare investigated using the stabilization method and correlated exponential wav...The effects of plasma screening on the ^(1)P^(o) resonance states of H-and He below the n=3 and n=4 thresholds of the respective subsystemsare investigated using the stabilization method and correlated exponential wave functions.Two plasma mediums,namely,the Debye plasma and quantum plasma environments are considered.The screened Coulomb potential(SCP)obtained from Debye-Hückel model is used to represent Debye plasma environments and the exponential cosine screened Coulomb potential(ECSCP)obtained from a modified Debye-Hückel model is used to represent quantum plasma environments.The resonance parameters(resonance positions and widths)are presented in terms of the screening parameters.展开更多
The concepts of complementary cofactor pairs, normal double-graphs and feasible torn vertex seta are introduced. By using them a decomposition theorem for first-order cofactor C(Y) is derived. Combining it with the mo...The concepts of complementary cofactor pairs, normal double-graphs and feasible torn vertex seta are introduced. By using them a decomposition theorem for first-order cofactor C(Y) is derived. Combining it with the modified double-graph method, a new decomposition analysis-modified double-graph decomposition analysis is presented for finding symbolic network functions. Its advantages are that the resultant symbolic expressions are compact and contain no cancellation terms, and its sign evaluation is very simple.展开更多
In this paper, symbolic code matrix ,constant matrix and count matrix are defined .The first twomatrices are used to describe the elemental expression of augmented matrix and the nede admittance equa-tion is thus obta...In this paper, symbolic code matrix ,constant matrix and count matrix are defined .The first twomatrices are used to describe the elemental expression of augmented matrix and the nede admittance equa-tion is thus obtained. The third matrix is used to obtain the incoming degree matrix, and according to thematrix all the 1- factors of the Coates graph are given. By using the data code, the determinant is expandedand the same items in the expansion are merged. Thus the symbolic network function in which no term can-cellation occurs is generated.展开更多
In a recent article [Commun. Theor. Phys. （Beijing, China） 47 （2007） 270], Cao et al. gave some nontrivial solutions of a Riccati equation by using symbolic and algebra computation. They took these solutions, whic...In a recent article [Commun. Theor. Phys. （Beijing, China） 47 （2007） 270], Cao et al. gave some nontrivial solutions of a Riccati equation by using symbolic and algebra computation. They took these solutions, which are in the form of q-deformed hyperbolic and triangular functions as new solutions. In this comment, we will show that these solutions are just the special cases of some known solutions of the Riccati equation and thus they are not new solutions.展开更多
基金National Natural Science Foundation of China under Grant No.52278534the Sichuan Provincial Natural Science Foundation of China under Grant No.2022NSFSC0423。
文摘Quantifying the post-earthquake functional recovery of railway stations presents significant challenges.This paper first establishes a post-earthquake function calculation method for railway stations,encompassing the establishment of relationships between the station’s function and the damage state,function loss,and failure probability of components and professional equipment in each layer.Also,the“4 stages-6 sequences”post-earthquake repair method is present,taking into account the functional and structural characteristics of railway stations.Additionally,a novel piecewise function for the post-earthquake functional dynamic recovery of railway stations is developed.A case study is conducted on a typical railway station to demonstrate the analysis procedure.Results indicate that under fortification,rare,and extremely rare earthquake scenarios,the interlayer drift ratio(IDR)of the railway station were 1/276,1/143,and 1/52,respectively,and corresponding peak floor acceleration(PFA)were 6.31 m/s^(2),7.82 m/s^(2),and 8.57 m/s^(2),respectively.The post-earthquake function of the railway station was 93.21%,82.33%,and 64.16%of its initial function.The repair times were 6.66 days,18.65 days,and 37.42 days.The displacement-sensitive,non-structural components were identified as the most vulnerable to damage.And the first repair stage(R_(1))which was mainly used to repair structural components and non-structural transport components,accounted for the highest proportion of total repair time.
基金supported by the Open Funds for Hubei Key Laboratory of Marine Geological Resources,China University of Geosciences(No.MGR202308)the Natural Science Foundation of Shandong Province(No.ZR2020MD085)+3 种基金the National Natural Science Foundation of China(No.41821004)the Taishan Scholar Program(No.tstp2022114)the Shandong Provincial Natural Science Foundation(No.DKXZZ202206)the National Key Research and Development Program of China(No.2016YFC1402404).
文摘Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.
基金supported by the Advance Research Project of Civil Aerospace Technology(Grant No.D020304)National Nat-ural Science Foundation of China(Grant Nos.52205257 and U22B2083).
文摘This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12174048 and 12204128)。
文摘A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the reference functions method to describe the solution to the depth-separated wave equation approximately using parabolic cylinder functions.The conditions for the validity of this approximation are also discussed.Furthermore,a formula that incorporates waveguide effects for the modal group velocity is derived,revealing that boundary effects at very low frequencies can have a significant impact on the propagation characteristics of even low-order normal modes.The present method not only offers improved accuracy compared to the classical WKB approximation and the uniform asymptotic approximation based on Airy functions,but also provides a wider range of depth applicability.Additionally,this method exhibits strong agreement with numerical methods and offers valuable physical insights.Finally,the method is applied to the study of very-low-frequency sound propagation in the South China Sea,leading to sound transmission loss predictions that closely align with experimental observations.
基金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.
基金supported by the Center for Advanced Systems Understanding(CASUS),financed by Germany’s Federal Ministry of Education and Research(BMBF)and the Saxon State Government out of the State Budget approved by the Saxon State Parliamentfunding from the European Research Council(ERC)under the European Union’s Horizon 2022 research and innovation programme(Grant Agreement No.101076233,“PREXTREME”)funding from the European Union’s Just Transition Fund(JTF)within the project Röntgenlaser-Optimierung der Laserfusion(ROLF),Contract No.5086999001,co-financed by the Saxon State Government out of the State Budget approved by the Saxon State Parliament.
文摘Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instrument function,is measured in X-ray Thom-son scattering(XRTS)experiments,which allow the study of electronic structure properties at the microscopic level.Among the various ab initio methods,linear-response time-dependent density-functional theory(LR-TDDFT)is a key framework for simulating the DSF.The standard approach in LR-TDDFT for computing the DSF relies on the orbital representation.A significant drawback of this method is the unfavorable scaling of the number of required empty bands as the wavenumber increases,making LR-TDDFT impractical for modeling XRTS measurements over large energy scales,such as in backward scattering geometry.In this work,we consider and test an alternative approach to LR-TDDFT that employs the Liouville–Lanczos(LL)method for simulating the DSF of WDM.This approach does not require empty states and allows the DSF at large momentum transfer values and over a broad frequency range to be accessed.We compare the results obtained from the LL method with those from the solution of Dyson’s equation using the standard LR-TDDFT within the projector augmented-wave formalism for isochorically heated aluminum and warm dense hydrogen.Additionally,we utilize exact path integral Monte Carlo results for the imaginary-time density-density correlation function(ITCF)of warm dense hydrogen to rigorously benchmark the LL approach.We discuss the application of the LL method for calculating DSFs and ITCFs at different wavenumbers,the effects of pseudopotentials,and the role of Lorentzian smearing.The successful validation of the LL method under WDM conditions makes it a valuable addition to the ab initio simulation landscape,supporting experimental efforts and advancing WDM theory.
文摘Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct doubly-periodic solutions of the Zakharov-Kuznetsov equation, which has been derived by Gottwald as a two-dimensional model for nonlinear Rossby waves. When the modulus k →1, these solutions reduce to the solitary wave solutions of the equation.
基金financial support from the National Natural Science Foundation of China (Grant No. 12227901)the financial support from the National Natural Science Foundation of China (Grant Nos. 11974263 and 12174291)。
文摘We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.
文摘This paper elucidates the effectiveness of combining the Poincare-Lighthill-Kuo method (PLK method, for short) and symbolic computation. Firstly, the idea and history of the PLK method are briefly introduced. Then, the difficulty of intermediate expression swell, often encountered in symbolic computation, is outlined. For overcoming the difficulty, a semi-inverse algorithm was proposed by the author, with which the lengthy ports of intermediate expressions are first frozen in the form of symbols till the Fnal stage of seeking perturbation solutions. Tn discuss the applications of the above algorithm, the related work of the author and his research group on nonlinear oscillations and waves is concisely reviewed. The computer-extended perturbation solution of the Duffing equation shows that the asymptotic solution obtained with the PLK method possesses the convergence radius of 1 and thus the range of validity of the solution is considerably enlarged. The studies on internal solitary waves in stratified fluid and on the head-on collision between two solitary waves in a hyperelastic rod indicate that by means of the presented methods, very complicated manipulation, unconceivable in hand calculation, can be conducted and thus result in higher-order evolution equations and asymptotic solutions. The examples illustrate that the algorithm helps to realize the symbolic computation on micro-commputers. Finally, it is concluded that,vith the aid of symbolic computation, the vitality of the PLK method is greatly. Strengthened and at least for the solutions to conservative systems of oscillations and waves, it is a powerful tool.
基金Anhui Provincial Natural Science Foundation(2308085QD124)Anhui Province University Natural Science Research Project(GrantNo.2023AH050918)The University Outstanding Youth Talent Support Program of Anhui Province.
文摘This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.
文摘The objective in this presentation is to introduce some of the unique properties and applications of nullors in active circuit analysis and designs. The emphasis is to discuss the role nullors can play in symbolic representation of transfer functions. To show this we adopt the topological platform for the circuit analysis and use a recently developed Admittance Method (AM) to achieve the Sum of Tree Products (STP), replacing the determinant and cofactors of the Nodal Admittance Matrix (NAM) of the circuit. To construct a transfer function, we start with a given active circuit and convert all its controlled sources and I/O-ports to nullors. Now, with a solid nullor circuit (passive elements and nullors) we first eliminate the passive elements through AM operations. This produces the STPs. Second, the all-nullor circuit is then used to find the signs or the STPs. Finally, the transfer function (in symbolic, if chosen) is obtained from the ratio between the STPs.
基金the National Key Research and Development Program of China(Grant Number 2021YFB1714600)the National Natural Science Foundation of China(Grant Number 52075195)the Fundamental Research Funds for the Central Universities,China through Program No.2172019kfyXJJS078.
文摘With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.
基金the National Natural Science Foundation of China(Nos.12302007,12372006,and 12202109)the Specific Research Project of Guangxi for Research Bases and Talents(No.AD23026051)。
文摘Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.
基金The project partially supported by National Natural Science Foundation of China under Grant No. 10471143 and the State 973 Project under Grant No. 2004CB318001 The authors are very grateful to Prof. Hong-Bo Li, Yong Chen, Zhen-Ya Yan, and Zhuo-Sheng Lii for their kind help and valuable suggestions. They also thank Prof. En-Gui Fan and Prof. Chun-Ping Liu for their constructive suggestions about the solutions of Riccati equation.
文摘In this paper, by using symbolic and algebra computation, Chen and Wang＇s multiple R/ccati equations rational expansion method was further extended. Many double soliton-like and other novel combined forms of exact solutions of the （2＋1）-dimensional Breaking soliton equation are derived by using the extended multiple Riccatl equations expansion method.
基金The National Key Research and Development Program of China(No.2023YFC3805003)。
文摘To quantify the seismic resilience of buildings,a method for evaluating functional loss from the component level to the overall building is proposed,and the dual-parameter seismic resilience assessment method based on postearthquake loss and recovery time is improved.A threelevel function tree model is established,which can consider the dynamic changes in weight coefficients of different category of components relative to their functional losses.Bayesian networks are utilized to quantify the impact of weather conditions,construction technology levels,and worker skill levels on component repair time.A method for determining the real-time functional recovery curve of buildings based on the component repair process is proposed.Taking a three-story teaching building as an example,the seismic resilience indices under basic earthquakes and rare earthquakes are calculated.The results show that the seismic resilience grade of the teaching building is comprehensively judged as GradeⅢ,and its resilience grade is more significantly affected by postearthquake loss.The proposed method can be used to predict the seismic resilience of buildings prior to earthquakes,identify weak components within buildings,and provide guidance for taking measures to enhance the seismic resilience of buildings.
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
基金Supported by the Natural Science Foundation of Heilongjiang Province(LH2024A025)。
文摘The effects of plasma screening on the ^(1)P^(o) resonance states of H-and He below the n=3 and n=4 thresholds of the respective subsystemsare investigated using the stabilization method and correlated exponential wave functions.Two plasma mediums,namely,the Debye plasma and quantum plasma environments are considered.The screened Coulomb potential(SCP)obtained from Debye-Hückel model is used to represent Debye plasma environments and the exponential cosine screened Coulomb potential(ECSCP)obtained from a modified Debye-Hückel model is used to represent quantum plasma environments.The resonance parameters(resonance positions and widths)are presented in terms of the screening parameters.
文摘The concepts of complementary cofactor pairs, normal double-graphs and feasible torn vertex seta are introduced. By using them a decomposition theorem for first-order cofactor C(Y) is derived. Combining it with the modified double-graph method, a new decomposition analysis-modified double-graph decomposition analysis is presented for finding symbolic network functions. Its advantages are that the resultant symbolic expressions are compact and contain no cancellation terms, and its sign evaluation is very simple.
基金The Project Supported by National Natural Science Foundation of China
文摘In this paper, symbolic code matrix ,constant matrix and count matrix are defined .The first twomatrices are used to describe the elemental expression of augmented matrix and the nede admittance equa-tion is thus obtained. The third matrix is used to obtain the incoming degree matrix, and according to thematrix all the 1- factors of the Coates graph are given. By using the data code, the determinant is expandedand the same items in the expansion are merged. Thus the symbolic network function in which no term can-cellation occurs is generated.
基金supported by National Natural Science Foundation of China under Grant No.10671172the Natural Science Foundation of Jiangsu Province under Grant No.BK2006064
文摘In a recent article [Commun. Theor. Phys. （Beijing, China） 47 （2007） 270], Cao et al. gave some nontrivial solutions of a Riccati equation by using symbolic and algebra computation. They took these solutions, which are in the form of q-deformed hyperbolic and triangular functions as new solutions. In this comment, we will show that these solutions are just the special cases of some known solutions of the Riccati equation and thus they are not new solutions.
