In the research of fractal cities, the fractal dimension is very important. It is used to describe the fractal character of the city. The authors have designed two approaches to calculate the fractal dimension by the ...In the research of fractal cities, the fractal dimension is very important. It is used to describe the fractal character of the city. The authors have designed two approaches to calculate the fractal dimension by the box-counting method through an example of Beijing, which are called the vector method and the grid method, respectively. The former calculates the fractal dimension through an intersecting analysis in ArcView; and the latter is carried out by programming in Matlab. They are compared from three aspects： the calculating process, the limits in use, and the results. As a result, the conclusion is made that there are merits and faults on both methods, and they should be chosen to use properly in practical situation.展开更多
Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This ...Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This paper proposes a successive shift box-counting method,in which the studied object is divided into small sub-objects that are composed of a series of gridsaccording to its characteristic scaling. The terrain fractal dimensions in the grids are calculatedwith the successive shift box-counting method and the scattered points with values of fractaldimensions are obtained. The present research shows that the planar variation of fractal dimensionsis well consistent with fault traces and geological boundaries.展开更多
Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media withi...Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.展开更多
In this paper, the box-counting dimension is used to derive an explicit formula for the dimension of a fractal constructed using several contractions or by combining fractals. This dimension agrees with the Hausdorff ...In this paper, the box-counting dimension is used to derive an explicit formula for the dimension of a fractal constructed using several contractions or by combining fractals. This dimension agrees with the Hausdorff dimension in the particular case when the scales factors considered are all the same. A more general sufficient condition for the box-counting dimension and the Hausdorff dimension to be the same is given. It is also shown that the dimension of the fractal obtained by combining two fractals is the weighted average of the dimensions of the two fractals.展开更多
This study successfully deals with the inhomogeneous dimension problem of load separation assumption, which is the theoretical basis of the normalization method. According to the dimensionless load separation principl...This study successfully deals with the inhomogeneous dimension problem of load separation assumption, which is the theoretical basis of the normalization method. According to the dimensionless load separation principle, the normalization method has been improved by intro- ducing a forcible blunting correction. With the improved normalization method, the J-resistance curves of five different metallic materials of CT and SEB specimens are estimated. The forcible blunting correction of initial crack size plays an important role in the J-resistance curve estima- tion, which is closely related to the strain hardening level of material. The higher level of strain hardening leads to a greater difference in JQ determined by different slopes of the blunting line. If the blunting line coefficient recommended by ASTM E1820-11 is used in the improved nor- realization method, it will lead to greater fracture resistance than that processed by the blunting line coefficient recommended by ISO 12135-2002. Therefore, the influence of the blunting line on the determination of JQ must be taken into full account in the fracture toughness assessment of metallic materials.展开更多
This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
In the reliability analysis of slope, the performance functions derived from the most available stability analysis procedures of slopes are usually implicit and cannot be solved by first-order second-moment approach. ...In the reliability analysis of slope, the performance functions derived from the most available stability analysis procedures of slopes are usually implicit and cannot be solved by first-order second-moment approach. A new reliability analysis approach was presented based on three-dimensional Morgenstem-Price method to investigate three-dimensional effect of landslide in stability analyses. To obtain the reliability index, Support Vector Machine （SVM） was applied to approximate the performance function. The time-consuming of this approach is only 0.028% of that using Monte-Carlo method at the same computation accuracy. Also, the influence of time effect of shearing strength parameters of slope soils on the long-term reliability of three-dimensional slopes was investigated by this new approach. It is found that the reliability index of the slope would decrease by 52.54% and the failure probability would increase from 0.000 705% to 1.966%. In the end, the impact of variation coefficients of c andfon reliability index of slopes was taken into discussion and the changing trend was observed.展开更多
We describe how the Unit-Feature Spatial Classification Method(UFSCM) can be used operationally to classify cloud types in satellite imagery efficiently and conveniently.By using a combination of Interactive Data Lang...We describe how the Unit-Feature Spatial Classification Method(UFSCM) can be used operationally to classify cloud types in satellite imagery efficiently and conveniently.By using a combination of Interactive Data Language(IDL) and Visual C++(VC) code in combination to extend the technique in three dimensions(3-D),this paper provides an efficient method to implement interactive computer visualization of the 3-D discrimination matrix modification,so as to deal with the bi-spectral limitations of traditional two dimensional(2-D) UFSCM.The case study of cloud-type classification based on FY-2C satellite data (0600 UTC 18 and 0000 UTC 10 September 2007) is conducted by comparison with ground station data, and indicates that 3-D UFSCM makes more use of the pattern recognition information in multi-spectral imagery,resulting in more reasonable results and an improvement over the 2-D method.展开更多
By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is propose...By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.展开更多
In this paper,we shall establish the superconvergence properties of the Runge-Kutta dis-continuous Galerkin method for solving two-dimensional linear constant hyperbolic equa-tion,where the upwind-biased numerical flu...In this paper,we shall establish the superconvergence properties of the Runge-Kutta dis-continuous Galerkin method for solving two-dimensional linear constant hyperbolic equa-tion,where the upwind-biased numerical flux is used.By suitably defining the correction function and deeply understanding the mechanisms when the spatial derivatives and the correction manipulations are carried out along the same or different directions,we obtain the superconvergence results on the node averages,the numerical fluxes,the cell averages,the solution and the spatial derivatives.The superconvergence properties in space are pre-served as the semi-discrete method,and time discretization solely produces an optimal order error in time.Some numerical experiments also are given.展开更多
In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on...In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on the stream surface and 2-D Navier-Stokes equations on a two dimensional manifold. After introducing stream function on the stream surface,a nonlinear initial-boundary value problem satisfies by stream function is obtained,existence and uniqueness of its solution are proven.Based this theory we proposal a new method called"dimension split method"to solve 3D NSE.展开更多
In order to overcome many problems and shortcomings of our three-dimensional design education in mechanical drawing, make the students master the modem design method and technology easily, modem modeling method and sh...In order to overcome many problems and shortcomings of our three-dimensional design education in mechanical drawing, make the students master the modem design method and technology easily, modem modeling method and shaping method are introduced, and the shaping designing methods of assembled-body are analyzed including basic requirements for shaping design, common principles of shaping and shaping design method. It will be helpful to improve the spatial imagination and modelling ability of student.展开更多
Fracture surface contour study is one of the important requirements for characterization and evaluation of the microstructure of rocks.Based on the improved cube covering method and the 3D contour digital reconstructi...Fracture surface contour study is one of the important requirements for characterization and evaluation of the microstructure of rocks.Based on the improved cube covering method and the 3D contour digital reconstruction model,this study proposes a quantitative microstructure characterization method combining the roughness evaluation index and the 3D fractal dimension to study the change rule of the fracture surface morphology after blasting.This method was applied and validated in the study of the fracture microstructure of the rock after blasting.The results show that the fracture morphology characteristics of the 3D contour digital reconstruction model have good correlation with the changes of the blasting action.The undulation rate of the three-dimensional surface profile of the rock is more prone to dramatic rise and dramatic fall morphology.In terms of tilting trend,the tilting direction also shows gradual disorder,with the tilting angle increasing correspondingly.All the roughness evaluation indexes of the rock fissure surface after blasting show a linear and gradually increasing trend as the distance to the bursting center increases;the difference between the two-dimensional roughness evaluation indexes and the three-dimensional ones of the same micro-area rock samples also becomes increasingly larger,among which the three-dimensional fissure roughness coefficient JRC and the surface roughness ratio Rs display better correlation.Compared with the linear fitting formula of the power function relationship,the three-dimensional fractal dimension of the postblast fissure surface is fitted with the values of JRC and Rs,which renders higher correlation coefficients,and the degree of linear fitting of JRC to the three-dimensional fractal dimension is higher.The fractal characteristics of the blast-affected region form a unity with the three-dimensional roughness evaluation of the fissure surface.展开更多
Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such...Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such as detonation and penetration, the dynamic parallel method （DPM） is designed to adjust the computational domain dynamically to get better load balance. Dynamic parallel method can be separated into two parts： one is division of initial computational domain and location of the data, the other is expansion of the computational domain and adjustment of the data location. DPM program can greatly shorten computational time and be preferable in simulating actual problems. The speedup of the DPM program is linear in parallel test. DPM can be popularized to parallel program of other multi-component high dimension Eulerian methods naturally.展开更多
The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are est...The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.展开更多
This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing t...This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing the order of approximation in fixed number of subdomains,rather than by resorting to a further partitioning.The stability and the convergence of this method are proved.展开更多
The article gives a semi-discrete method for solving high-dimension wave equationBy the method, high-dimension wave equation is converted by, means of diseretizationinto I-D wave equation system which is well-posed. T...The article gives a semi-discrete method for solving high-dimension wave equationBy the method, high-dimension wave equation is converted by, means of diseretizationinto I-D wave equation system which is well-posed. The convergence of the semidijcrete method is given. The numerical calculating resulis show that the speed of convergence is high.展开更多
Finite element method(FEM) simulations were employed to investigate the quenching residual stress distributions of 7085 aluminum alloy plates.The effect of dimensional variation on the quenching residual stress distri...Finite element method(FEM) simulations were employed to investigate the quenching residual stress distributions of 7085 aluminum alloy plates.The effect of dimensional variation on the quenching residual stress distributions was studied and discussed by using models with different dimensions(length,width,and thickness).The accuracy and efficiency of the models were verified by other numerical examples.The order of the dimension effects on the quenching residual stress distributions is:thickness> width=length.The maximum tensile stress and compressive stress increase from 33 to 190 and 39 to 270 MPa,respectively,as the thickness increases from 30 to 150 mm.The ultimate maximum tensile stress(about190 MPa) is equivalent to half of the quenching yield strength at 20℃,while the ultimate maximum compressive stress(about 300 MPa) is equivalent to 80 % of the quenching yield strength at 20℃.There are stress fluctuations at the edge of the large plate both in rolling and in transverse directions.The ratio of the fluctuation region along the rolling direction and transverse direction increases as the thickness increases,while it decreases as the length or width increases.The actual length of the fluctuation region is almost a constant value for the plates with a thickness of 115 mm(about 500 mm in length and 300 mm in width).展开更多
The invariant subspace method is used to construct the explicit solution of a nonlinear evolution equation. The second-order nonlinear differential operators that possess invariant subspaces of submaximal dimension ar...The invariant subspace method is used to construct the explicit solution of a nonlinear evolution equation. The second-order nonlinear differential operators that possess invariant subspaces of submaximal dimension are described. There are second-order nonlinear differential operators, including cubic operators and quadratic operators, which preserve an invariant subspace of submaximal dimension. A full. description, of the second-order cubic operators with constant coefficients admitting a four-dimensional invariant subspace is given. It is shown that the maximal dimension of invaxiant subspaces preserved by a second-order cubic operator is four. Several examples are given for the construction of the exact solutions to nonlinear evolution equations with cubic nonlinearities. These solutions blow up in a finite展开更多
Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the h...Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the help of intrinsically safe borehole video instrument.And according to the geology of this working face,a discrete element model was calculated by UDEC.Combining in situ experimental data with numerical results,the relationship between the fractal dimension of boreholes'wall and the distribution of advanced abutment pressure was studied under the condition of mining advance.The results show that the variation tendency of fractal dimension and the abutment pressure has the same characteristic value.The distance between working face and the peak value of the abutment pressure has a slight increasing trend with the advancing of mining-face.When the working face is set as the original point,the trend of fractal dimension from the far place to the origin can be divided into three phases:constant,steady increasing and constant.And the turning points of these phases are the max-influencing distance(50 m)and peak value(15 m)of abutment pressure.展开更多
文摘In the research of fractal cities, the fractal dimension is very important. It is used to describe the fractal character of the city. The authors have designed two approaches to calculate the fractal dimension by the box-counting method through an example of Beijing, which are called the vector method and the grid method, respectively. The former calculates the fractal dimension through an intersecting analysis in ArcView; and the latter is carried out by programming in Matlab. They are compared from three aspects： the calculating process, the limits in use, and the results. As a result, the conclusion is made that there are merits and faults on both methods, and they should be chosen to use properly in practical situation.
文摘Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This paper proposes a successive shift box-counting method,in which the studied object is divided into small sub-objects that are composed of a series of gridsaccording to its characteristic scaling. The terrain fractal dimensions in the grids are calculatedwith the successive shift box-counting method and the scattered points with values of fractaldimensions are obtained. The present research shows that the planar variation of fractal dimensionsis well consistent with fault traces and geological boundaries.
基金supported by the National Natural Science Foundation of China (Nos.52374078 and 52074043)the Fundamental Research Funds for the Central Universities (No.2023CDJKYJH021)。
文摘Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.
文摘In this paper, the box-counting dimension is used to derive an explicit formula for the dimension of a fractal constructed using several contractions or by combining fractals. This dimension agrees with the Hausdorff dimension in the particular case when the scales factors considered are all the same. A more general sufficient condition for the box-counting dimension and the Hausdorff dimension to be the same is given. It is also shown that the dimension of the fractal obtained by combining two fractals is the weighted average of the dimensions of the two fractals.
基金supported by the National Natural Science Foundation of China(Nos.11472228 and 11202174)the Sichuan Provincial Youth Science and Technology Innovation Team(No.2013TD0004)
文摘This study successfully deals with the inhomogeneous dimension problem of load separation assumption, which is the theoretical basis of the normalization method. According to the dimensionless load separation principle, the normalization method has been improved by intro- ducing a forcible blunting correction. With the improved normalization method, the J-resistance curves of five different metallic materials of CT and SEB specimens are estimated. The forcible blunting correction of initial crack size plays an important role in the J-resistance curve estima- tion, which is closely related to the strain hardening level of material. The higher level of strain hardening leads to a greater difference in JQ determined by different slopes of the blunting line. If the blunting line coefficient recommended by ASTM E1820-11 is used in the improved nor- realization method, it will lead to greater fracture resistance than that processed by the blunting line coefficient recommended by ISO 12135-2002. Therefore, the influence of the blunting line on the determination of JQ must be taken into full account in the fracture toughness assessment of metallic materials.
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
基金Project(50878082) supported by the National Natural Science Foundation of ChinaProject(200631880237) supported by the Science and Technology Program of West Transportation of the Ministry of Transportation of ChinaKey Project(09JJ3104) supported by the Natural Science Foundation of Hunan Province, China
文摘In the reliability analysis of slope, the performance functions derived from the most available stability analysis procedures of slopes are usually implicit and cannot be solved by first-order second-moment approach. A new reliability analysis approach was presented based on three-dimensional Morgenstem-Price method to investigate three-dimensional effect of landslide in stability analyses. To obtain the reliability index, Support Vector Machine （SVM） was applied to approximate the performance function. The time-consuming of this approach is only 0.028% of that using Monte-Carlo method at the same computation accuracy. Also, the influence of time effect of shearing strength parameters of slope soils on the long-term reliability of three-dimensional slopes was investigated by this new approach. It is found that the reliability index of the slope would decrease by 52.54% and the failure probability would increase from 0.000 705% to 1.966%. In the end, the impact of variation coefficients of c andfon reliability index of slopes was taken into discussion and the changing trend was observed.
基金supported by the National Natural Science Foundation of China(Grant No.40875012)the National Basic Research Program of China(Grant No.2009CB421502)the Meteorology Open Fund of Huaihe River Basin(HRM200704).
文摘We describe how the Unit-Feature Spatial Classification Method(UFSCM) can be used operationally to classify cloud types in satellite imagery efficiently and conveniently.By using a combination of Interactive Data Language(IDL) and Visual C++(VC) code in combination to extend the technique in three dimensions(3-D),this paper provides an efficient method to implement interactive computer visualization of the 3-D discrimination matrix modification,so as to deal with the bi-spectral limitations of traditional two dimensional(2-D) UFSCM.The case study of cloud-type classification based on FY-2C satellite data (0600 UTC 18 and 0000 UTC 10 September 2007) is conducted by comparison with ground station data, and indicates that 3-D UFSCM makes more use of the pattern recognition information in multi-spectral imagery,resulting in more reasonable results and an improvement over the 2-D method.
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.
基金Yuan Xu is supported by the NSFC Grant 11671199Qiang Zhang is supported by the NSFC Grant 11671199.
文摘In this paper,we shall establish the superconvergence properties of the Runge-Kutta dis-continuous Galerkin method for solving two-dimensional linear constant hyperbolic equa-tion,where the upwind-biased numerical flux is used.By suitably defining the correction function and deeply understanding the mechanisms when the spatial derivatives and the correction manipulations are carried out along the same or different directions,we obtain the superconvergence results on the node averages,the numerical fluxes,the cell averages,the solution and the spatial derivatives.The superconvergence properties in space are pre-served as the semi-discrete method,and time discretization solely produces an optimal order error in time.Some numerical experiments also are given.
文摘In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on the stream surface and 2-D Navier-Stokes equations on a two dimensional manifold. After introducing stream function on the stream surface,a nonlinear initial-boundary value problem satisfies by stream function is obtained,existence and uniqueness of its solution are proven.Based this theory we proposal a new method called"dimension split method"to solve 3D NSE.
基金Supported by 2015 Undergraduate Education and Teaching Reform and Research Project of University of Science and Technology Beijing(JG2015M15)Fundamental Research Funds for the Central Universities(FRF-BR-15-047A)
文摘In order to overcome many problems and shortcomings of our three-dimensional design education in mechanical drawing, make the students master the modem design method and technology easily, modem modeling method and shaping method are introduced, and the shaping designing methods of assembled-body are analyzed including basic requirements for shaping design, common principles of shaping and shaping design method. It will be helpful to improve the spatial imagination and modelling ability of student.
基金National Key Research and Development Program of China,Grant/Award Number:2021YFC2902103National Natural Science Foundation of China,Grant/Award Number:51934001Fundamental Research Funds for the Central Universities,Grant/Award Number:2023JCCXLJ02。
文摘Fracture surface contour study is one of the important requirements for characterization and evaluation of the microstructure of rocks.Based on the improved cube covering method and the 3D contour digital reconstruction model,this study proposes a quantitative microstructure characterization method combining the roughness evaluation index and the 3D fractal dimension to study the change rule of the fracture surface morphology after blasting.This method was applied and validated in the study of the fracture microstructure of the rock after blasting.The results show that the fracture morphology characteristics of the 3D contour digital reconstruction model have good correlation with the changes of the blasting action.The undulation rate of the three-dimensional surface profile of the rock is more prone to dramatic rise and dramatic fall morphology.In terms of tilting trend,the tilting direction also shows gradual disorder,with the tilting angle increasing correspondingly.All the roughness evaluation indexes of the rock fissure surface after blasting show a linear and gradually increasing trend as the distance to the bursting center increases;the difference between the two-dimensional roughness evaluation indexes and the three-dimensional ones of the same micro-area rock samples also becomes increasingly larger,among which the three-dimensional fissure roughness coefficient JRC and the surface roughness ratio Rs display better correlation.Compared with the linear fitting formula of the power function relationship,the three-dimensional fractal dimension of the postblast fissure surface is fitted with the values of JRC and Rs,which renders higher correlation coefficients,and the degree of linear fitting of JRC to the three-dimensional fractal dimension is higher.The fractal characteristics of the blast-affected region form a unity with the three-dimensional roughness evaluation of the fissure surface.
基金Sponsored by State Key Laboratory of Computational Physics Fundation(9140C690101070C69)
文摘Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such as detonation and penetration, the dynamic parallel method （DPM） is designed to adjust the computational domain dynamically to get better load balance. Dynamic parallel method can be separated into two parts： one is division of initial computational domain and location of the data, the other is expansion of the computational domain and adjustment of the data location. DPM program can greatly shorten computational time and be preferable in simulating actual problems. The speedup of the DPM program is linear in parallel test. DPM can be popularized to parallel program of other multi-component high dimension Eulerian methods naturally.
基金the Science Foundation of the Science and Technology Commission of Shanghai Municipality(No.075105118)the Shanghai Leading Academic Discipline Project(No.T0401)the Fund for E-institute of Shanghai Universities(No.E03004)
文摘The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.
文摘This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing the order of approximation in fixed number of subdomains,rather than by resorting to a further partitioning.The stability and the convergence of this method are proved.
文摘The article gives a semi-discrete method for solving high-dimension wave equationBy the method, high-dimension wave equation is converted by, means of diseretizationinto I-D wave equation system which is well-posed. The convergence of the semidijcrete method is given. The numerical calculating resulis show that the speed of convergence is high.
基金financially supported by the National Basic Research Program of China(No.2012CB619504)the National Natural Science Foundation of China(No.51274046)
文摘Finite element method(FEM) simulations were employed to investigate the quenching residual stress distributions of 7085 aluminum alloy plates.The effect of dimensional variation on the quenching residual stress distributions was studied and discussed by using models with different dimensions(length,width,and thickness).The accuracy and efficiency of the models were verified by other numerical examples.The order of the dimension effects on the quenching residual stress distributions is:thickness> width=length.The maximum tensile stress and compressive stress increase from 33 to 190 and 39 to 270 MPa,respectively,as the thickness increases from 30 to 150 mm.The ultimate maximum tensile stress(about190 MPa) is equivalent to half of the quenching yield strength at 20℃,while the ultimate maximum compressive stress(about 300 MPa) is equivalent to 80 % of the quenching yield strength at 20℃.There are stress fluctuations at the edge of the large plate both in rolling and in transverse directions.The ratio of the fluctuation region along the rolling direction and transverse direction increases as the thickness increases,while it decreases as the length or width increases.The actual length of the fluctuation region is almost a constant value for the plates with a thickness of 115 mm(about 500 mm in length and 300 mm in width).
基金Project supported by the National Natural Science Foundation of China(Grant No.10926082)the Natural Science Foundation of Anhui Province of China(Grant No.KJ2010A128)the Fund for Youth of Anhui Normal University,China(Grant No.2009xqn55)
文摘The invariant subspace method is used to construct the explicit solution of a nonlinear evolution equation. The second-order nonlinear differential operators that possess invariant subspaces of submaximal dimension are described. There are second-order nonlinear differential operators, including cubic operators and quadratic operators, which preserve an invariant subspace of submaximal dimension. A full. description, of the second-order cubic operators with constant coefficients admitting a four-dimensional invariant subspace is given. It is shown that the maximal dimension of invaxiant subspaces preserved by a second-order cubic operator is four. Several examples are given for the construction of the exact solutions to nonlinear evolution equations with cubic nonlinearities. These solutions blow up in a finite
基金financial support from the State Key Basic Research Program of China(Nos.2011CB201201and 2010CB226802)the National Natural Science Foundation of China(No.51204112)
文摘Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the help of intrinsically safe borehole video instrument.And according to the geology of this working face,a discrete element model was calculated by UDEC.Combining in situ experimental data with numerical results,the relationship between the fractal dimension of boreholes'wall and the distribution of advanced abutment pressure was studied under the condition of mining advance.The results show that the variation tendency of fractal dimension and the abutment pressure has the same characteristic value.The distance between working face and the peak value of the abutment pressure has a slight increasing trend with the advancing of mining-face.When the working face is set as the original point,the trend of fractal dimension from the far place to the origin can be divided into three phases:constant,steady increasing and constant.And the turning points of these phases are the max-influencing distance(50 m)and peak value(15 m)of abutment pressure.
