The crack tip strain gauge method in the compliance technique was used to determine the opening load of notched crack of axle steel,and the nonlinear finite element ADINA program, to which the cyclic stress-strain cur...The crack tip strain gauge method in the compliance technique was used to determine the opening load of notched crack of axle steel,and the nonlinear finite element ADINA program, to which the cyclic stress-strain curve of axle steel was applied,was used to analyze the stress-strain field ahead of the crack tip and the opening load of notched crack.The results of both the compliance technique and the numerical method were in good agreement.In this pa- per,the concept of the sensitive point is proposed and the key to the determination of the crack opening load in the experiment is to place a strain gauge at sensitive point.It is certified by both experimental and numerical methods that the sensitive point has the best linear relation- ship character and the value of strain is much greater.展开更多
A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. F...A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. First, the principle of TRL technique is described. Then, the availability of TRL solver for nonlinear application is particularly discussed demonstrating that this solver can be easily used while still remaining great efficiency. The implementation on how to apply this technique in FEM based on magnetic vector potential （MVP） is also introduced. Finally, a numerical example of 3-D magnetostatic modeling using the TRL solver and FEMLAB is given. It shows that a huge computer resource can be saved by employing the new solver.展开更多
Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harm...Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.展开更多
Two kinds of fractures can be observed in the SPH （smoothed particle hydrodynamics） simulations, which are the physical fracture and the numerical fracture. The physical one exists in reality, while the numerical on...Two kinds of fractures can be observed in the SPH （smoothed particle hydrodynamics） simulations, which are the physical fracture and the numerical fracture. The physical one exists in reality, while the numerical one is fictitious. This paper presents the effects of both fractures and proposes a simple adding particle technique to avoid the numerical fracture. The real physical fracture is then figured out by using an applicable fracture criterion. Firstly, the effect of the numerical fracture on the computational accuracy is investigated by introducing the artificial fracture in a model of wave propagation. Secondly, a simple adding particle technique is proposed and validated by a three dimensional bending test. Finally, the experiments of penetration on the skin of aircrafts are simulated by both the initial SPH method and the improved method with the adding particle technique. The results show that the improved SPH method can describe the physical fracture very well with better accuracy.展开更多
In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier ser...In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.展开更多
Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressur...Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.展开更多
The success of LISA Pathfinder in demonstrating the LISA drag-free requirement paved the way for using space interferometers to detect low-frequency and middle-frequency gravitational waves(GWs). The TAIJI GW mission ...The success of LISA Pathfinder in demonstrating the LISA drag-free requirement paved the way for using space interferometers to detect low-frequency and middle-frequency gravitational waves(GWs). The TAIJI GW mission and the new LISA GW mission propose using an arm length of 3 Gm(1 Gm = 10~6 km) and an arm length of 2.5 Gm respectively. For a space laser-interferometric GW antenna,due to astrodynamical orbit variation, time delay interferometry(TDI) is needed to achieve nearly equivalent equal-arms for suppressing the laser frequency noise below the level of optical path noise, acceleration noise, etc in order to attain the requisite sensitivity. In this paper, we simulate TDI numerically for the TAIJI mission and the new LISA mission. To do this, we work out a set of 2200-day(6-year) optimized science orbits for each mission starting on 2028 March 22 using the CGC 2.7.1 ephemeris framework. Then we use the numerical method to calculate the residual optical path differences of the first-generation TDI configurations and the selected second-generation TDI configurations. The resulting optical path differences of the second-generation TDI configurations calculated for TAIJI, new LISA and eLISA are well below their respective requirements for laser frequency noise cancelation. However, for the first-generation TDI configurations, the original requirements need to be relaxed by 3 to 30 fold to be satisfied. For TAIJI and the new LISA, about one order of magnitude relaxation would be good and recommended; this could be borne on the laser stability requirement in view of recent progress in laser stability, or the GW detection sensitivities of the second-generation TDIs have to be used in the diagnosis of the observed data instead of the commonly used X, Y and Z TDIs.展开更多
The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the m...The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the motion of the free surface, the COF method with geometric reconstruction is employed. The velocity of the flow is measured by means of LDV technique. Four types of flow patterns, the position of the separation and reattachment point, the distribution of shear stress on the bed at downstream of the weir are presented and discussed. The numerical results agree well with the experiment data.展开更多
The 1927 Gulang M_S8.0 earthquake is a severe earthquake that followed the Haiyuan M_S8.5 earthquake of 1920 in the Qilian Mt._Hexi Corridor earthquake zone. There are divergences of opinion in the previous studies ab...The 1927 Gulang M_S8.0 earthquake is a severe earthquake that followed the Haiyuan M_S8.5 earthquake of 1920 in the Qilian Mt._Hexi Corridor earthquake zone. There are divergences of opinion in the previous studies about the rupture properties of the earthquake. Based on trenching and field investigation, and analysis of historical data, we hold that the earthquake resulted from the joint process of the Tianqiaogou_Huangyangchuan fault, Dongqingding segment of the Huangcheng_Shuangta fault and the Wuwei_Tianzhu buried fault, which constitute the Gulang nappe. By finite_element numerical simulation on the deformation mechanism of Gulang nappe, it is found that the stress and strain mainly concentrate in the western segment of the Tianqiaogou_Huangyangchuan fault, the Dongqingding segment of the Huangcheng_Shuangta fault, and the Gulangxia segment of the Wuwei_Tianzhu buried fault and the Gulang_Shuangta fault. The stress concentration coincides with the distribution of the earthquake surface rupture. It also proves that the earthquake is an outcome of the Gulang nappe activity as a whole.展开更多
In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Comb...In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.展开更多
The 21 cm radiation of neutral hydrogen provides crucial information for studying the early universe and its evolution.To advance this research,countries have made significant investments in constructing large lowfreq...The 21 cm radiation of neutral hydrogen provides crucial information for studying the early universe and its evolution.To advance this research,countries have made significant investments in constructing large lowfrequency radio telescope arrays,such as the Low Frequency Array and the Square Kilometre Array Phase 1 Low Frequency.These instruments are pivotal for radio astronomy research.However,challenges such as ionospheric plasma interference,ambient radio noise,and instrument-related effects have become increasingly prominent,posing major obstacles in cosmology research.To address these issues,this paper proposes an efficient signal processing method that combines wavelet transform and mathematical morphology.The method involves the following steps:Background Subtraction:Background interference in radio observation signals is eliminated.Wavelet Transform:The signal,after removing background noise,undergoes a two-dimensional discrete wavelet transform.Threshold processing is then applied to the wavelet coefficients to effectively remove interference components.Wavelet Inversion:The processed signal is reconstructed using wavelet inversion.Mathematical Morphology:The reconstructed signal is further optimized using mathematical morphology to refine the results.Experimental verification was conducted using solar observation data from the Xinjiang Observatory and the Yunnan Observatory.The results demonstrate that this method successfully removes interference signals while preserving useful signals,thus improving the accuracy of radio astronomy observations and reducing the impact of radio frequency interference.展开更多
A new numerical scheme for solving the tidal flow in an opening channel using the advective-diffusion shallow-water equations as the governing equations is proposed based on the combination of the MacCormack and the f...A new numerical scheme for solving the tidal flow in an opening channel using the advective-diffusion shallow-water equations as the governing equations is proposed based on the combination of the MacCormack and the finite analysis methods. In the present scheme, the finite analysis method is used to discretize the momentum equation and the MacCormack technique is used to discretize the continuity equation in a single grid system. The matrix of the discretized momentum equation is characterized by predominantly main diagonal elements, which ensures favorable convergence and stability for the numerical simulation by the combined method. To verify the present method, hydraulics simulation is carried out for a section down mainstream of the Huangpu River. The computational results agree with the measured data. By use of orthogonal curvilinear coordinate system, the methods can be easily extended to the numerical simulation of the tidal flow in a tortuous channel.展开更多
The wave-based method （WBM） has been applied for the prediction of mid-frequency vibrations of fiat plates. The scaling factors, Gauss point selection rule and truncation rule are introduced to insure the wave model...The wave-based method （WBM） has been applied for the prediction of mid-frequency vibrations of fiat plates. The scaling factors, Gauss point selection rule and truncation rule are introduced to insure the wave model to converge. Numerical results show that the prediction tech- nique based on WBM is with higher accuracy and smaller computational effort than the one on FEM, which implies that this new technique on WBM can be applied to higher-frequency range.展开更多
Geosynthetic-reinforced and pile-supported （GRPS） embankment has been increasingly constructed in a large number of regions and for a wide range of projects in the past decades. However, many disadvantages are expos...Geosynthetic-reinforced and pile-supported （GRPS） embankment has been increasingly constructed in a large number of regions and for a wide range of projects in the past decades. However, many disadvantages are exposed through a lot of applications on conventional technique of GRPS embankment （called CT embankment）, i.e., intolerable settlement and lateral displacement, low geosynthetic efficiency, etc. In view of these disadvantages, the fixed geosynthetic technique of GRPS embankment （called FGT embankment） is developed in this work. In this system, the geosynthetic is fixed on the pile head by the steel bar fulcrum and concrete fixed top. The principles and construction techniques involved in the FGT embankment are described firstly. Then, the numerical analysis method and two-stage analysis method are used to study the performance of FGT embankment, respectively. It is shown that the FGT embankment can provide a better improvement technique to construct a high embankment over soft ground.展开更多
A new Runge-Kutta （PK） fourth order with four stages embedded method with error control is presentea m this paper for raster simulation in cellular neural network （CNN） environment. Through versatile algorithm, si...A new Runge-Kutta （PK） fourth order with four stages embedded method with error control is presentea m this paper for raster simulation in cellular neural network （CNN） environment. Through versatile algorithm, single layer/raster CNN array is implemented by incorporating the proposed technique. Simulation results have been obtained, and comparison has also been carried out to show the efficiency of the proposed numerical integration algorithm. The analytic expressions for local truncation error and global truncation error are derived. It is seen that the RK-embedded root mean square outperforms the RK-embedded Heronian mean and RK-embedded harmonic mean.展开更多
Spectrum denoising is an important procedure for large-scale spectroscopical surveys. This work proposes a novel stellar spectrum denoising method based on deep Bayesian modeling. The construction of our model include...Spectrum denoising is an important procedure for large-scale spectroscopical surveys. This work proposes a novel stellar spectrum denoising method based on deep Bayesian modeling. The construction of our model includes a prior distribution for each stellar subclass, a spectrum generator and a flow-based noise model. Our method takes into account the noise correlation structure, and it is not susceptible to strong sky emission lines and cosmic rays. Moreover, it is able to naturally handle spectra with missing flux values without ad-hoc imputation. The proposed method is evaluated on real stellar spectra from the Sloan Digital Sky Survey(SDSS) with a comprehensive list of common stellar subclasses and compared to the standard denoising auto-encoder. Our denoising method demonstrates a superior performance to the standard denoising auto-encoder, in respect of denoising quality and missing flux imputation. It may be potentially helpful in improving the accuracy of the classification and physical parameter measurement of stars when applying our method during data preprocessing.展开更多
The hard X-ray modulation telescope （HXMT） mission is mainly devoted to performing an all-sky survey at 1- 250 keV with both high sensitivity and high spatial resolution. The observed data reduction as well as the i...The hard X-ray modulation telescope （HXMT） mission is mainly devoted to performing an all-sky survey at 1- 250 keV with both high sensitivity and high spatial resolution. The observed data reduction as well as the image reconstruction for HXMT can be achieved by using the direct demodulation method （DDM）. However the original DDM is too computationally expensive for multi-dimensional data with high resolution to be employed for HXMT data. We propose an accelerated direct demodulation method especially adapted for data from HXMT. Simulations are also presented to demonstrate this method.展开更多
A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling ...A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling variational formulation and the convergence of the approximate solution corresponding to the algorithm are proved. If the fine mesh h is chosen as coarse mesh H-square, the nonlinear Galerkin method, nonlinearity is only treated on the coarse grid and linearity is treated on the fine grid. Hence, the new algorithm can save a large amount of computational time.展开更多
In this paper TVD explicit schemes developed recently in CFD are used to simulate five models of bores due to dam-break. Numerical experiments are presented to demonstrate that they achieve high resolution and generat...In this paper TVD explicit schemes developed recently in CFD are used to simulate five models of bores due to dam-break. Numerical experiments are presented to demonstrate that they achieve high resolution and generate no spurious oscillations and to show the versatility of the schemes.展开更多
The perturbation method and finite strip method are combined to solve the large deflection bending problems of rectangular plates. Perturbation method is used to reduce the nonlinear differential equations into a seri...The perturbation method and finite strip method are combined to solve the large deflection bending problems of rectangular plates. Perturbation method is used to reduce the nonlinear differential equations into a series of linear differential equations. The finite strip method is then employed to tackle these linear equations. Some calculation examples are compared with those got by other methods.展开更多
文摘The crack tip strain gauge method in the compliance technique was used to determine the opening load of notched crack of axle steel,and the nonlinear finite element ADINA program, to which the cyclic stress-strain curve of axle steel was applied,was used to analyze the stress-strain field ahead of the crack tip and the opening load of notched crack.The results of both the compliance technique and the numerical method were in good agreement.In this pa- per,the concept of the sensitive point is proposed and the key to the determination of the crack opening load in the experiment is to place a strain gauge at sensitive point.It is certified by both experimental and numerical methods that the sensitive point has the best linear relation- ship character and the value of strain is much greater.
基金Open Funds of State Key Laboratory of MillimeterWaves,China (No. K200401), Outstanding Teaching and ResearchAwards for Young Teachers of Nanjing Normal University (No.1320BL51)
文摘A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. First, the principle of TRL technique is described. Then, the availability of TRL solver for nonlinear application is particularly discussed demonstrating that this solver can be easily used while still remaining great efficiency. The implementation on how to apply this technique in FEM based on magnetic vector potential （MVP） is also introduced. Finally, a numerical example of 3-D magnetostatic modeling using the TRL solver and FEMLAB is given. It shows that a huge computer resource can be saved by employing the new solver.
基金Project supported by the National Natural Science Foundation of China (No.10172038)
文摘Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.
基金supported by the National Natural Science Foundation of China (No.10577016)the 111 Project (No. B07050)the program for 2008 New Century Excellent Talents in University (No.NCET080454)
文摘Two kinds of fractures can be observed in the SPH （smoothed particle hydrodynamics） simulations, which are the physical fracture and the numerical fracture. The physical one exists in reality, while the numerical one is fictitious. This paper presents the effects of both fractures and proposes a simple adding particle technique to avoid the numerical fracture. The real physical fracture is then figured out by using an applicable fracture criterion. Firstly, the effect of the numerical fracture on the computational accuracy is investigated by introducing the artificial fracture in a model of wave propagation. Secondly, a simple adding particle technique is proposed and validated by a three dimensional bending test. Finally, the experiments of penetration on the skin of aircrafts are simulated by both the initial SPH method and the improved method with the adding particle technique. The results show that the improved SPH method can describe the physical fracture very well with better accuracy.
文摘In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.
基金Project supported by the National Natural Science Foundation of China(No.5140232)the National Science and Technology Major Project(No.2011ZX05038003)the China Postdoctoral Science Foundation(No.2014M561074)
文摘Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.
基金funding in support of his work leading to these results from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme FP7/2007-2013/ (PEOPLE-2013-ITN) under REA grant agreement n [606176]
文摘The success of LISA Pathfinder in demonstrating the LISA drag-free requirement paved the way for using space interferometers to detect low-frequency and middle-frequency gravitational waves(GWs). The TAIJI GW mission and the new LISA GW mission propose using an arm length of 3 Gm(1 Gm = 10~6 km) and an arm length of 2.5 Gm respectively. For a space laser-interferometric GW antenna,due to astrodynamical orbit variation, time delay interferometry(TDI) is needed to achieve nearly equivalent equal-arms for suppressing the laser frequency noise below the level of optical path noise, acceleration noise, etc in order to attain the requisite sensitivity. In this paper, we simulate TDI numerically for the TAIJI mission and the new LISA mission. To do this, we work out a set of 2200-day(6-year) optimized science orbits for each mission starting on 2028 March 22 using the CGC 2.7.1 ephemeris framework. Then we use the numerical method to calculate the residual optical path differences of the first-generation TDI configurations and the selected second-generation TDI configurations. The resulting optical path differences of the second-generation TDI configurations calculated for TAIJI, new LISA and eLISA are well below their respective requirements for laser frequency noise cancelation. However, for the first-generation TDI configurations, the original requirements need to be relaxed by 3 to 30 fold to be satisfied. For TAIJI and the new LISA, about one order of magnitude relaxation would be good and recommended; this could be borne on the laser stability requirement in view of recent progress in laser stability, or the GW detection sensitivities of the second-generation TDIs have to be used in the diagnosis of the observed data instead of the commonly used X, Y and Z TDIs.
文摘The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the motion of the free surface, the COF method with geometric reconstruction is employed. The velocity of the flow is measured by means of LDV technique. Four types of flow patterns, the position of the separation and reattachment point, the distribution of shear stress on the bed at downstream of the weir are presented and discussed. The numerical results agree well with the experiment data.
文摘The 1927 Gulang M_S8.0 earthquake is a severe earthquake that followed the Haiyuan M_S8.5 earthquake of 1920 in the Qilian Mt._Hexi Corridor earthquake zone. There are divergences of opinion in the previous studies about the rupture properties of the earthquake. Based on trenching and field investigation, and analysis of historical data, we hold that the earthquake resulted from the joint process of the Tianqiaogou_Huangyangchuan fault, Dongqingding segment of the Huangcheng_Shuangta fault and the Wuwei_Tianzhu buried fault, which constitute the Gulang nappe. By finite_element numerical simulation on the deformation mechanism of Gulang nappe, it is found that the stress and strain mainly concentrate in the western segment of the Tianqiaogou_Huangyangchuan fault, the Dongqingding segment of the Huangcheng_Shuangta fault, and the Gulangxia segment of the Wuwei_Tianzhu buried fault and the Gulang_Shuangta fault. The stress concentration coincides with the distribution of the earthquake surface rupture. It also proves that the earthquake is an outcome of the Gulang nappe activity as a whole.
文摘In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.
基金funded by the National Key Research and Development Program’s intergovernmental International Science and Technology Innovation Cooperation project,titled Remote Sensing and Radio Astronomy Observation of Space Weather in Low and Middle Latitudes(project number:2022YFE0140000)Supported by International Partnership Program of Chinese Academy of Sciences,grant No.114A11KYSB20200001。
文摘The 21 cm radiation of neutral hydrogen provides crucial information for studying the early universe and its evolution.To advance this research,countries have made significant investments in constructing large lowfrequency radio telescope arrays,such as the Low Frequency Array and the Square Kilometre Array Phase 1 Low Frequency.These instruments are pivotal for radio astronomy research.However,challenges such as ionospheric plasma interference,ambient radio noise,and instrument-related effects have become increasingly prominent,posing major obstacles in cosmology research.To address these issues,this paper proposes an efficient signal processing method that combines wavelet transform and mathematical morphology.The method involves the following steps:Background Subtraction:Background interference in radio observation signals is eliminated.Wavelet Transform:The signal,after removing background noise,undergoes a two-dimensional discrete wavelet transform.Threshold processing is then applied to the wavelet coefficients to effectively remove interference components.Wavelet Inversion:The processed signal is reconstructed using wavelet inversion.Mathematical Morphology:The reconstructed signal is further optimized using mathematical morphology to refine the results.Experimental verification was conducted using solar observation data from the Xinjiang Observatory and the Yunnan Observatory.The results demonstrate that this method successfully removes interference signals while preserving useful signals,thus improving the accuracy of radio astronomy observations and reducing the impact of radio frequency interference.
基金supported by the National Natural Science Foundation of China (Grant No.10572084)the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No.20050280008)
文摘A new numerical scheme for solving the tidal flow in an opening channel using the advective-diffusion shallow-water equations as the governing equations is proposed based on the combination of the MacCormack and the finite analysis methods. In the present scheme, the finite analysis method is used to discretize the momentum equation and the MacCormack technique is used to discretize the continuity equation in a single grid system. The matrix of the discretized momentum equation is characterized by predominantly main diagonal elements, which ensures favorable convergence and stability for the numerical simulation by the combined method. To verify the present method, hydraulics simulation is carried out for a section down mainstream of the Huangpu River. The computational results agree with the measured data. By use of orthogonal curvilinear coordinate system, the methods can be easily extended to the numerical simulation of the tidal flow in a tortuous channel.
基金Project supported by the National Natural Science Foundation of China (No.10472035).
文摘The wave-based method （WBM） has been applied for the prediction of mid-frequency vibrations of fiat plates. The scaling factors, Gauss point selection rule and truncation rule are introduced to insure the wave model to converge. Numerical results show that the prediction tech- nique based on WBM is with higher accuracy and smaller computational effort than the one on FEM, which implies that this new technique on WBM can be applied to higher-frequency range.
基金Foundation item： Project（51278216） supported by the National Natural Science Foundation of China Project（11-2-05） supported by the Scientific and Technological Project for Shanxi Communication Construction, China Project（HF-08-01-2011-240） supported by the Graduates＇ Innovation Fund of Huazhong University of Science and Technology, China
文摘Geosynthetic-reinforced and pile-supported （GRPS） embankment has been increasingly constructed in a large number of regions and for a wide range of projects in the past decades. However, many disadvantages are exposed through a lot of applications on conventional technique of GRPS embankment （called CT embankment）, i.e., intolerable settlement and lateral displacement, low geosynthetic efficiency, etc. In view of these disadvantages, the fixed geosynthetic technique of GRPS embankment （called FGT embankment） is developed in this work. In this system, the geosynthetic is fixed on the pile head by the steel bar fulcrum and concrete fixed top. The principles and construction techniques involved in the FGT embankment are described firstly. Then, the numerical analysis method and two-stage analysis method are used to study the performance of FGT embankment, respectively. It is shown that the FGT embankment can provide a better improvement technique to construct a high embankment over soft ground.
基金supported as a part of Technical Quality Improvement Programme (TEQIP)
文摘A new Runge-Kutta （PK） fourth order with four stages embedded method with error control is presentea m this paper for raster simulation in cellular neural network （CNN） environment. Through versatile algorithm, single layer/raster CNN array is implemented by incorporating the proposed technique. Simulation results have been obtained, and comparison has also been carried out to show the efficiency of the proposed numerical integration algorithm. The analytic expressions for local truncation error and global truncation error are derived. It is seen that the RK-embedded root mean square outperforms the RK-embedded Heronian mean and RK-embedded harmonic mean.
基金funded by the National Natural Science Foundation of China(Grant Nos.11873066 and U1731109)。
文摘Spectrum denoising is an important procedure for large-scale spectroscopical surveys. This work proposes a novel stellar spectrum denoising method based on deep Bayesian modeling. The construction of our model includes a prior distribution for each stellar subclass, a spectrum generator and a flow-based noise model. Our method takes into account the noise correlation structure, and it is not susceptible to strong sky emission lines and cosmic rays. Moreover, it is able to naturally handle spectra with missing flux values without ad-hoc imputation. The proposed method is evaluated on real stellar spectra from the Sloan Digital Sky Survey(SDSS) with a comprehensive list of common stellar subclasses and compared to the standard denoising auto-encoder. Our denoising method demonstrates a superior performance to the standard denoising auto-encoder, in respect of denoising quality and missing flux imputation. It may be potentially helpful in improving the accuracy of the classification and physical parameter measurement of stars when applying our method during data preprocessing.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11173038 and 11103022)the Tsinghua University Initiative Scientific Research Program (Grant No. 20111081102)
文摘The hard X-ray modulation telescope （HXMT） mission is mainly devoted to performing an all-sky survey at 1- 250 keV with both high sensitivity and high spatial resolution. The observed data reduction as well as the image reconstruction for HXMT can be achieved by using the direct demodulation method （DDM）. However the original DDM is too computationally expensive for multi-dimensional data with high resolution to be employed for HXMT data. We propose an accelerated direct demodulation method especially adapted for data from HXMT. Simulations are also presented to demonstrate this method.
文摘A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling variational formulation and the convergence of the approximate solution corresponding to the algorithm are proved. If the fine mesh h is chosen as coarse mesh H-square, the nonlinear Galerkin method, nonlinearity is only treated on the coarse grid and linearity is treated on the fine grid. Hence, the new algorithm can save a large amount of computational time.
文摘In this paper TVD explicit schemes developed recently in CFD are used to simulate five models of bores due to dam-break. Numerical experiments are presented to demonstrate that they achieve high resolution and generate no spurious oscillations and to show the versatility of the schemes.
文摘The perturbation method and finite strip method are combined to solve the large deflection bending problems of rectangular plates. Perturbation method is used to reduce the nonlinear differential equations into a series of linear differential equations. The finite strip method is then employed to tackle these linear equations. Some calculation examples are compared with those got by other methods.
