The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experi...The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.展开更多
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.展开更多
In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation...In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.展开更多
To investigate the spatiotemporal variations in the mixed layer depth(MLD)in the Arctic basins,a new criterion to determine the MLD,called the improved maximum angle method(IMAM),was developed.A total of 45123 potenti...To investigate the spatiotemporal variations in the mixed layer depth(MLD)in the Arctic basins,a new criterion to determine the MLD,called the improved maximum angle method(IMAM),was developed.A total of 45123 potential density profiles collected using Ice-Tethered Profilers(ITPs)in the Arctic basins during 2005-2021 were used to demonstrate the method’s effectiveness.By comparing the results obtained by the fixed threshold method(FTM),percentage threshold method(PTM),and maximum gradient method(MGM)for profiles in the Canada Basin,Makarov Basin,and Eurasian Basin,we determined that the quality index(1.0 for perfect identification of the MLD)of the IMAM regarding the assessment of the MLD determination method reached 0.94,which is much greater than those of other criteria.Moreover,two types of the density profiles were identified based on the mixed layer development stage.The MLDs of the typical profiles determined using the IMAM were found to have better consistency with the original definition.By utilizing the new mixed layer criterion,the seasonal variations and regional differences in the MLD in the Arctic basins were analyzed.Spatially,the summer and winter MLDs in the Canada Basin were the shallowest(13.55 m in summer,26.76 m in winter)than those in the Makarov(29.51 m in summer,49.08 m in winter)and Eurasian(20.36 m in summer,46.81 m in winter)basins due to the stable stratification in the upper ocean and the subsequent small effects of dynamic and thermodynamic processes(wind-driven stirring and brine rejection)in the Canada Basin.Seasonally,in the three Arctic basins,the average MLD was shallowest(22.77 m)in summer;it deepened through autumn and reached a winter maximum(41.12 m).展开更多
In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square (MLS) approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin (IIEFG) method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin (EFG) method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (B...The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.展开更多
The improved element-free Galerkin (IEFG) method of elasticity is used to solve the topology optimization problems. In this method, the improved moving least-squares approximation is used to form the shape function....The improved element-free Galerkin (IEFG) method of elasticity is used to solve the topology optimization problems. In this method, the improved moving least-squares approximation is used to form the shape function. In a topology opti- mization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin (EFG) method. The central processing unit (CPU) time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.展开更多
By employing the improved moving least-square (IMLS) approximation, the improved element-free Galerkin (IEFG) method is presented for the unsteady Schrodinger equation. In the IEFG method, the two-dimensional (2D...By employing the improved moving least-square (IMLS) approximation, the improved element-free Galerkin (IEFG) method is presented for the unsteady Schrodinger equation. In the IEFG method, the two-dimensional (2D) trial function is approximated by the IMLS approximation, the variation method is used to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. Because the number of coefficients in the IMLS approximation is less than in the moving least-square (MLS) approximation, fewer nodes are needed in the entire domain when the IMLS approximation is used than when the MLS approximation is adopted. Then the IEFG method has high computational efficiency and accuracy. Several numerical examples are given to verify the accuracy and efficiency of the IEFG method in this paper.展开更多
In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means t...In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.展开更多
In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved c...In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin(ICVEFG) method is presented for two-dimensional(2D) elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.展开更多
An experimental method of measuring the real contact area of transparent blocks based on the principle of total internal reflection is presented, intending to support the investigation of friction characteristics, hea...An experimental method of measuring the real contact area of transparent blocks based on the principle of total internal reflection is presented, intending to support the investigation of friction characteristics, heat conduction, and energy dissipation at the contact interface. A laser sheet illuminates the contact interface, and the transmitted laser sheet is projected onto a screen. Then the contact information is acquired from the screen by a camera. An improved Otsu method is proposed to process the data of experimental images. It can compute the threshold of the overall image and filter out all the pixels one by one. Through analyzing the experimental results, we describe the relationship between the real contact area and the positive pressure during a continuous loading process, at different loading rates, with the polymethyl methacrylate(PMMA)material. A hysteresis phenomenon in the relationship between the real contact area and the positive pressure is found and explained.展开更多
The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, w...The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.展开更多
In this paper, the improved complex variable moving least-squares (ICVMLS) approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squar...In this paper, the improved complex variable moving least-squares (ICVMLS) approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares (CVMLS) approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin (EFG) method and the ICVMLS approximation, the improved complex variable element-free Galerkin (ICVEFG) method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.展开更多
The present study proposes an improved method for generating better typical meteorological years( TMYs) used in building energy simulation in China. Modifications are made to the commonly used Sandia method,by optimiz...The present study proposes an improved method for generating better typical meteorological years( TMYs) used in building energy simulation in China. Modifications are made to the commonly used Sandia method,by optimizing the weightings of indices and thus giving more emphasis to dry bulb temperature and relative humidity and less to wind velocity. After analyzing the solar heat gain on the vertical envelop,an index of diffuse radiation rather than direct normal radiation is added for solar radiation. Using the improved method proposed,TMYs for 5 representative cities of 5 major climate zones in China are generated from the meteorological data recorded during the period 1981—2010. The results show that,compared with previous studies,the monthly diffuse solar radiation of typical meteorological months generated by the improved method are the"closest"to the 30-year average,and the comparison between annual diffuse radiation for the TMY and the 30-year annual average is improved,while little adverse effect is produced on global horizontal radiation comparisons,indicating that the improved method is more suitable to generate the TMY data than previous studies.展开更多
Based on the improved interpolating moving least-squares (ⅡMLS) method and the Galerkin weak form, an improved interpolating element-free Galerkin (ⅡEFG) method is presented for two-dimensional elasticity proble...Based on the improved interpolating moving least-squares (ⅡMLS) method and the Galerkin weak form, an improved interpolating element-free Galerkin (ⅡEFG) method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares (IMLS) method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin (IEFG) method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.展开更多
In this paper, an improved splitting method, based on the completely square-conservative explicit difference schemes, is established. Not only can the time-direction precision of this method be higher than that of the...In this paper, an improved splitting method, based on the completely square-conservative explicit difference schemes, is established. Not only can the time-direction precision of this method be higher than that of the traditional splitting methods but also can the physical feature of mutual dependence of the fast and the slow stages that are calculated separately and splittingly be kept as well. Moreover, the method owns an universality, it can be generalized to other square-conservative difference schemes, such as the implicit and complete ones and the explicit and instantaneous ones. Good time benefits can be acquired when it is applied in the numerical simulations of the monthly mean currents of the South China Sea.展开更多
In this paper, we analyze the generalized Camassa and Holm (CH) equation by the improved element-free Galerkin (IEFG) method. By employing the improved moving least-square (IMLS) approximation, we derive the for...In this paper, we analyze the generalized Camassa and Holm (CH) equation by the improved element-free Galerkin (IEFG) method. By employing the improved moving least-square (IMLS) approximation, we derive the formulas for the generalized CH equation with the IEFG method. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Because there are fewer coefficients in the IMLS approximation than in the MLS approximation, and in the IEFG method, fewer nodes are selected in the entire domain than in the conventional EFG method, the IEFG method should result in a higher computing speed. The effectiveness of the IEFG method for the generalized CH equation is investigated by numerical examples in this paper.展开更多
Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering (CAE) and computer aided manufacturing ...Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering (CAE) and computer aided manufacturing (CAM). This paper presents a high-efficiency improved symmetric successive over-relaxation (ISSOR) preconditioned conjugate gradient (PCG) method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.展开更多
The accuracy of unstructured finite volume methods is greatly influenced by the gradient reconstruction, for which the stencil selection plays a critical role. Compared with the commonly used face-neighbor and vertex-...The accuracy of unstructured finite volume methods is greatly influenced by the gradient reconstruction, for which the stencil selection plays a critical role. Compared with the commonly used face-neighbor and vertex-neighbor stencils, the global-direction stencil is independent of the mesh topology, and characteristics of the flow field can be well reflected by this novel stencil. However, for a high-aspect-ratio triangular grid, the grid skewness is evident, which is one of the most important grid-quality measures known to affect the accuracy and stability of finite volume solvers. On this basis and inspired by an approach of using face-area-weighted centroid to reduce the grid skewness, we explore a method by combining the global-direction stencil and face-area-weighted centroid on high-aspect-ratio triangular grids, so as to improve the computational accuracy. Four representative numerical cases are simulated on high-aspect-ratio triangular grids to examine the validity of the improved global-direction stencil. Results illustrate that errors of this improved methods are the lowest among all methods we tested, and in high-mach-number flow, with the increase of cell aspect ratio, the improved global-direction stencil always has a better stability than commonly used face-neighbor and vertex-neighbor stencils. Therefore, the computational accuracy as well as stability is greatly improved, and superiorities of this novel method are verified.展开更多
A variational method is developed to retrieve winds in the first step and then thermodynamic fields in the second step from Doppler radar observations. In the first step, wind fields are retrieved at two time levels: ...A variational method is developed to retrieve winds in the first step and then thermodynamic fields in the second step from Doppler radar observations. In the first step, wind fields are retrieved at two time levels: the beginning and ending times of the data assimilation period, simultaneously from two successive volume scans by using the weak form constraints provided by the mass continuity and vorticity equations. As the retrieved wind fields are expressed by Legendre polynomial expansions at the beginning and ending times, the time tendency term in the vorticity equation can be conveniently formulated, and the retrieved winds can be compared with the radar observed radial winds in the cost function at the precise time and position of each radar beam. In the second step, the perturbation pressure and temperature fields at the middle time are then derived from the retrieved wind fields and the velocity time tendency by using the weak form constraints provided by the three momentum equations. The merits of the new method are demonstrated by numerical experiments with simulated radar observations and compared with the traditional least squares methods which consider neither the precise observation times and positions nor the velocity time tendency. The new method is also applied to real radar data for a heavy rainfall event during the 2001 Meiyu season in China.展开更多
基金supported by National Natural Science Foundation of China(Nos.92066108 and 51277061)。
文摘The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.
基金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.
基金supported by the Key Projects of Natural Science Foundation of China(No.41931284)the Scientific Research Start-Up Fund for High-Level Introduced Talents of Anhui University of Science and Technology(No.2022yjrc21).
文摘In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.
基金Supported by the National Key R&D Program of China(Nos.2018 YFA 0605903,2019 YFC 1509101)the National Natural Science Foundation of China(No.41976218)the Fundamental Research Funds for the Central Universities(No.202165005)。
文摘To investigate the spatiotemporal variations in the mixed layer depth(MLD)in the Arctic basins,a new criterion to determine the MLD,called the improved maximum angle method(IMAM),was developed.A total of 45123 potential density profiles collected using Ice-Tethered Profilers(ITPs)in the Arctic basins during 2005-2021 were used to demonstrate the method’s effectiveness.By comparing the results obtained by the fixed threshold method(FTM),percentage threshold method(PTM),and maximum gradient method(MGM)for profiles in the Canada Basin,Makarov Basin,and Eurasian Basin,we determined that the quality index(1.0 for perfect identification of the MLD)of the IMAM regarding the assessment of the MLD determination method reached 0.94,which is much greater than those of other criteria.Moreover,two types of the density profiles were identified based on the mixed layer development stage.The MLDs of the typical profiles determined using the IMAM were found to have better consistency with the original definition.By utilizing the new mixed layer criterion,the seasonal variations and regional differences in the MLD in the Arctic basins were analyzed.Spatially,the summer and winter MLDs in the Canada Basin were the shallowest(13.55 m in summer,26.76 m in winter)than those in the Makarov(29.51 m in summer,49.08 m in winter)and Eurasian(20.36 m in summer,46.81 m in winter)basins due to the stable stratification in the upper ocean and the subsequent small effects of dynamic and thermodynamic processes(wind-driven stirring and brine rejection)in the Canada Basin.Seasonally,in the three Arctic basins,the average MLD was shallowest(22.77 m)in summer;it deepened through autumn and reached a winter maximum(41.12 m).
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square (MLS) approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin (IIEFG) method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin (EFG) method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金Project supported by the National Natural Science Foundation of China (Grant No 10871124)Innovation Program of Shanghai Municipal Education Commission (Grant No 09ZZ99)Shanghai Leading Academic Discipline Project (Grant No J50103)
文摘The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.
基金supported by the National Natural Science Foundation of China(Grant Nos.11571223 and U1433104)
文摘The improved element-free Galerkin (IEFG) method of elasticity is used to solve the topology optimization problems. In this method, the improved moving least-squares approximation is used to form the shape function. In a topology opti- mization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin (EFG) method. The central processing unit (CPU) time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Natural Science Foundation of Zhejiang Province,China(Grant No.LY15A020007)+1 种基金the Natural Science Foundation of Ningbo City(Grant No.2014A610028)the K.C.Wong Magna Fund in Ningbo University,China
文摘By employing the improved moving least-square (IMLS) approximation, the improved element-free Galerkin (IEFG) method is presented for the unsteady Schrodinger equation. In the IEFG method, the two-dimensional (2D) trial function is approximated by the IMLS approximation, the variation method is used to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. Because the number of coefficients in the IMLS approximation is less than in the moving least-square (MLS) approximation, fewer nodes are needed in the entire domain when the IMLS approximation is used than when the MLS approximation is adopted. Then the IEFG method has high computational efficiency and accuracy. Several numerical examples are given to verify the accuracy and efficiency of the IEFG method in this paper.
基金supported by the National Key Technology R&D Program (Grant 2011BAJ02B01-02)the National Natural Science Foundation of China (Grant 11602065)
文摘In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11171208 and U1433104)
文摘In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin(ICVEFG) method is presented for two-dimensional(2D) elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.
基金Project supported by the National Natural Science Foundation of China(Grant No.11272171)the Beijing Natural Science Foundation,China(Grant No.3132030)the Education Ministry Doctoral Fund of China(Grant No.20120002110070)
文摘An experimental method of measuring the real contact area of transparent blocks based on the principle of total internal reflection is presented, intending to support the investigation of friction characteristics, heat conduction, and energy dissipation at the contact interface. A laser sheet illuminates the contact interface, and the transmitted laser sheet is projected onto a screen. Then the contact information is acquired from the screen by a camera. An improved Otsu method is proposed to process the data of experimental images. It can compute the threshold of the overall image and filter out all the pixels one by one. Through analyzing the experimental results, we describe the relationship between the real contact area and the positive pressure during a continuous loading process, at different loading rates, with the polymethyl methacrylate(PMMA)material. A hysteresis phenomenon in the relationship between the real contact area and the positive pressure is found and explained.
基金supported by the National Natural Science Foundation of China (11171208)Shanghai Leading Academic Discipline Project (S30106)
文摘The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.
基金supported by the National Natural Science Foundation of China (Grant No.11026223)the Shanghai Leading Academic Discipline Project,China (Grant No.S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No.SHUCX112359)
文摘In this paper, the improved complex variable moving least-squares (ICVMLS) approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares (CVMLS) approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin (EFG) method and the ICVMLS approximation, the improved complex variable element-free Galerkin (ICVEFG) method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51278349)Foundation of Fujian Provincial Education Department(Grant No.JA13185)
文摘The present study proposes an improved method for generating better typical meteorological years( TMYs) used in building energy simulation in China. Modifications are made to the commonly used Sandia method,by optimizing the weightings of indices and thus giving more emphasis to dry bulb temperature and relative humidity and less to wind velocity. After analyzing the solar heat gain on the vertical envelop,an index of diffuse radiation rather than direct normal radiation is added for solar radiation. Using the improved method proposed,TMYs for 5 representative cities of 5 major climate zones in China are generated from the meteorological data recorded during the period 1981—2010. The results show that,compared with previous studies,the monthly diffuse solar radiation of typical meteorological months generated by the improved method are the"closest"to the 30-year average,and the comparison between annual diffuse radiation for the TMY and the 30-year annual average is improved,while little adverse effect is produced on global horizontal radiation comparisons,indicating that the improved method is more suitable to generate the TMY data than previous studies.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)
文摘Based on the improved interpolating moving least-squares (ⅡMLS) method and the Galerkin weak form, an improved interpolating element-free Galerkin (ⅡEFG) method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares (IMLS) method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin (IEFG) method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.
基金Partly supported by the State Major Key Project for Basic Researches
文摘In this paper, an improved splitting method, based on the completely square-conservative explicit difference schemes, is established. Not only can the time-direction precision of this method be higher than that of the traditional splitting methods but also can the physical feature of mutual dependence of the fast and the slow stages that are calculated separately and splittingly be kept as well. Moreover, the method owns an universality, it can be generalized to other square-conservative difference schemes, such as the implicit and complete ones and the explicit and instantaneous ones. Good time benefits can be acquired when it is applied in the numerical simulations of the monthly mean currents of the South China Sea.
基金supported by the Natural Science Foundation of Ningbo City,Zhejiang Province,China (Grant Nos. 2012A610038 and 2012A610023)the Natural Science Foundation of Zhejiang Province,China (Grant No. Y6110007)
文摘In this paper, we analyze the generalized Camassa and Holm (CH) equation by the improved element-free Galerkin (IEFG) method. By employing the improved moving least-square (IMLS) approximation, we derive the formulas for the generalized CH equation with the IEFG method. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Because there are fewer coefficients in the IMLS approximation than in the MLS approximation, and in the IEFG method, fewer nodes are selected in the entire domain than in the conventional EFG method, the IEFG method should result in a higher computing speed. The effectiveness of the IEFG method for the generalized CH equation is investigated by numerical examples in this paper.
基金Project supported by the National Natural Science Foundation of China(Nos.5130926141030747+3 种基金41102181and 51121005)the National Basic Research Program of China(973 Program)(No.2011CB013503)the Young Teachers’ Initial Funding Scheme of Sun Yat-sen University(No.39000-1188140)
文摘Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering (CAE) and computer aided manufacturing (CAM). This paper presents a high-efficiency improved symmetric successive over-relaxation (ISSOR) preconditioned conjugate gradient (PCG) method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.
基金Project supported by the National Key Project, China (Grant No. GJXM92579).
文摘The accuracy of unstructured finite volume methods is greatly influenced by the gradient reconstruction, for which the stencil selection plays a critical role. Compared with the commonly used face-neighbor and vertex-neighbor stencils, the global-direction stencil is independent of the mesh topology, and characteristics of the flow field can be well reflected by this novel stencil. However, for a high-aspect-ratio triangular grid, the grid skewness is evident, which is one of the most important grid-quality measures known to affect the accuracy and stability of finite volume solvers. On this basis and inspired by an approach of using face-area-weighted centroid to reduce the grid skewness, we explore a method by combining the global-direction stencil and face-area-weighted centroid on high-aspect-ratio triangular grids, so as to improve the computational accuracy. Four representative numerical cases are simulated on high-aspect-ratio triangular grids to examine the validity of the improved global-direction stencil. Results illustrate that errors of this improved methods are the lowest among all methods we tested, and in high-mach-number flow, with the increase of cell aspect ratio, the improved global-direction stencil always has a better stability than commonly used face-neighbor and vertex-neighbor stencils. Therefore, the computational accuracy as well as stability is greatly improved, and superiorities of this novel method are verified.
文摘A variational method is developed to retrieve winds in the first step and then thermodynamic fields in the second step from Doppler radar observations. In the first step, wind fields are retrieved at two time levels: the beginning and ending times of the data assimilation period, simultaneously from two successive volume scans by using the weak form constraints provided by the mass continuity and vorticity equations. As the retrieved wind fields are expressed by Legendre polynomial expansions at the beginning and ending times, the time tendency term in the vorticity equation can be conveniently formulated, and the retrieved winds can be compared with the radar observed radial winds in the cost function at the precise time and position of each radar beam. In the second step, the perturbation pressure and temperature fields at the middle time are then derived from the retrieved wind fields and the velocity time tendency by using the weak form constraints provided by the three momentum equations. The merits of the new method are demonstrated by numerical experiments with simulated radar observations and compared with the traditional least squares methods which consider neither the precise observation times and positions nor the velocity time tendency. The new method is also applied to real radar data for a heavy rainfall event during the 2001 Meiyu season in China.