Bench blasting is commonly used in open-pit mining.Some design parameters such as positions of hole packing and caving holes have great influences on the blasting effects.In this work,with a hybrid discrete-finite ele...Bench blasting is commonly used in open-pit mining.Some design parameters such as positions of hole packing and caving holes have great influences on the blasting effects.In this work,with a hybrid discrete-finite element method,numerical simulations of bench blasting are conducted,capturing the whole continuous-discontinuous processes.Considering two engineering cases,the influences of hole packing and caving holes are evaluated.The numerical results not only lead to some improved designs by relocating the packing positions and caving holes but also indicate the reliability of the adopted numerical tools.展开更多
The rigid body limit equilibrium method （LEM） and the nonlinear finite element method （NFEM） are often used in the analysis of anti-sliding stability of gravity dam. But LEM cannot reflect the process of progressi...The rigid body limit equilibrium method （LEM） and the nonlinear finite element method （NFEM） are often used in the analysis of anti-sliding stability of gravity dam. But LEM cannot reflect the process of progressive instability and mechanical mecha- nism on failure for rock mass while NFEM is difficult to use to solve the displacement discontinuity of weak structural plane. Combining the research with Xiangjiaba Hydropower Station project, the analysis of anti-sliding stability for segment 12# of the dam has been carried out using interface stress element method （ISEM）. The results can reflect the most dangerous location, the scope and distribution of failure zone in weak structural plane, and present the process of progressive failure in dam foun- dation as well as the safety coefficient of possible sliding body. These achievements provide an important technical reference for dam foundation treatment measures. The computational results show that ISEM can naturally describe discontinuous de- formation of rock mass such as dislocation, openness and sliding. Besides, this method is characterized by good adaptability, convenient calculation and high compatibility, thus it is regarded as an effective way to make an analysis of anti-sliding stabil- ity of gravity dam展开更多
We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations ...We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy o...In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.展开更多
Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle in...Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle induced areas.It is a scenario of strong coupled thermomechanical process involving discontinuity behaviours of rocks.In this context,a numerical model was proposed to investigate the thermal cracking of rocks,in a framework of the continuous-discontinuous element method(CDEM)for efficiently capturing the initiation and propagation of multiple cracks.A simplex integration strategy was adopted to account for the influences of temperature-dependent material properties.Several benchmark tests were considered and the obtained results were compared with analytical solutions and numerical results from the literature.The results show that the fracture degree of the cases when considering temperature-dependent material parameters had 10%differences approximately compared with the cases with constant parameters.展开更多
We deal with the numerical solution of the Navier-Stokes equations describing a motion of viscous compressible fluids.In order to obtain a sufficiently stable higher order scheme with respect to the time and space coo...We deal with the numerical solution of the Navier-Stokes equations describing a motion of viscous compressible fluids.In order to obtain a sufficiently stable higher order scheme with respect to the time and space coordinates,we develop a combination of the discontinuous Galerkin finite element(DGFE)method for the space discretization and the backward difference formulae(BDF)for the time discretization.Since the resulting discrete problem leads to a system of nonlinear algebraic equations at each time step,we employ suitable linearizations of inviscid as well as viscous fluxes which give a linear algebraic problem at each time step.Finally,the resulting BDF-DGFE scheme is applied to steady as well as unsteady flows and achieved results are compared with reference data.展开更多
The fatigue damage model based on theory of damage mechanics is capable of predicting the fatigue life under multiaxial loading. Meanwhile, the application of critical plane method in the prediction of multiaxial fati...The fatigue damage model based on theory of damage mechanics is capable of predicting the fatigue life under multiaxial loading. Meanwhile, the application of critical plane method in the prediction of multiaxial fatigue life has made certain progress. According to the law of thermodynamics, a new damage evolution equation is developed in the present study to predict the fatigue life of geometrically discontinuous structure under tension-torsion loading based on damage mechanics and the critical plane method. The essence of this approach is tha t the st rain parame ter of the uniaxial nonlinear fatigue damage model is replaced with the equivalent strain, which consists of the releva nt parame ters of the critical plane. However, it is difficult to calculate the stress-strain status and the critical plane position of geometrically dis? continuous structure by theoretical methods because of the existence of stress concentration and the multiaxial nonproportional characteristics. Therefore, a new numerical simulation method is proposed to determine the critical plane of geometrically discontinuous structure under multiaxial loading by means of the finite element method and MATLAB software. The fatigue life of notched specimens subjected to combined bending and torsion is predicted using the proposed met hod, and the result is compared with t hose from the experimen ts and the Manson-Cfiffin law. The comparisons show that the proposed method is superior to the Manson-Coffin law and is capable of reproducing the experimental results reasonably when the geometry of the structure is complex. It completely meets the needs of engineering practice.展开更多
The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in...The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in the S N method to automatically optimize the angular distribution and minimize angular discretization errors with lower expenses.The proposed method enables linear dis-continuous finite element quadrature sets over an icosahe-dron to vary their quadrature orders in a one-twentieth sphere so that fine resolutions can be applied to the angular domains that are important.An error estimation that operates in conjunction with the spherical harmonics method is developed to determine the locations where more refinement is required.The adaptive quadrature sets are applied to three duct problems,including the Kobayashi benchmarks and the IRI-TUB research reactor,which emphasize the ability of this method to resolve neutron streaming through ducts with voids.The results indicate that the performance of the adaptive method is more effi-cient than that of uniform quadrature sets for duct transport problems.Our adaptive method offers an appropriate placement of angular unknowns to accurately integrate angular fluxes while reducing the computational costs in terms of unknowns and run times.展开更多
A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations....A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.展开更多
In this paper we investigate the existence, uniqueness and regularity of the solution of semilinear parabolic equations with coefficients that are discontinuous across the interface, some prior estimates are obtained....In this paper we investigate the existence, uniqueness and regularity of the solution of semilinear parabolic equations with coefficients that are discontinuous across the interface, some prior estimates are obtained. A net shape of the finite elements around the singular points was designed in [7] to solve the linear elliptic problems, by means of that net, we prove that the approximate solution has the same convergence rate as that without singularity.展开更多
This paper deals with the shape and influenced factors of surface non-continuous deformation due to mining. With finite element method, analysis are made to derive the relations between discontinuous deformation and m...This paper deals with the shape and influenced factors of surface non-continuous deformation due to mining. With finite element method, analysis are made to derive the relations between discontinuous deformation and mining affection, weak plane’s position & thickness, and mechanical property of weak-plane medium. The mutual affection of multiple weak-planes is also discussed. The results of the paper lay a foundation for constructing the calculation method of surface discontinuous deformation.展开更多
In automotive industries,panel acoustic contribution analysis(PACA)is used to investigate the contributions of the body panels to the acoustic pressure at a certain point of interest.Currently,PACA is implementedmostl...In automotive industries,panel acoustic contribution analysis(PACA)is used to investigate the contributions of the body panels to the acoustic pressure at a certain point of interest.Currently,PACA is implementedmostly by either experiment-based methods or traditional numerical methods.However,these schemes are effort-consuming and inefficient in solving engineering problems,thereby restraining the further development of PACA in automotive acoustics.In this work,we propose a PACA scheme using discontinuous isogeometric boundary element method(IGABEM)to build an easily implementable and efficient method to identify the relative acoustic contributions of each automotive body panel.Discontinuous IGABEMis more accurate and converges faster than continuous BEM and IGABEM in the interior sound pressure evaluation of automotive compartments.In this work,a contribution ratio is defined to estimate the relative acoustic contribution of the structure panels;it can be calculated by reusing the coefficient matrix that has already been generated in the sound pressure evaluation process.The utilization of the parallel technique enables the proposed method to be more efficient than conventional methods;it is validated in two numerical examples,including a car passenger compartment subjected to realistic boundary conditions.A sound pressure response experiment based on a steel box is conducted to verify the accuracy of the interior sound pressure calculation using discontinuous IGABEM.This work is expected to promote the practical process of IGABEM for application in automotive acoustic problems.展开更多
基金the financial support by the National Natural Science Foundation of China(NSFC)(52178324).
文摘Bench blasting is commonly used in open-pit mining.Some design parameters such as positions of hole packing and caving holes have great influences on the blasting effects.In this work,with a hybrid discrete-finite element method,numerical simulations of bench blasting are conducted,capturing the whole continuous-discontinuous processes.Considering two engineering cases,the influences of hole packing and caving holes are evaluated.The numerical results not only lead to some improved designs by relocating the packing positions and caving holes but also indicate the reliability of the adopted numerical tools.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51179064, 11132003 and 10972072)the National Science and Technology Supporting Plan (Grant No. 2008BAB29B03)
文摘The rigid body limit equilibrium method （LEM） and the nonlinear finite element method （NFEM） are often used in the analysis of anti-sliding stability of gravity dam. But LEM cannot reflect the process of progressive instability and mechanical mecha- nism on failure for rock mass while NFEM is difficult to use to solve the displacement discontinuity of weak structural plane. Combining the research with Xiangjiaba Hydropower Station project, the analysis of anti-sliding stability for segment 12# of the dam has been carried out using interface stress element method （ISEM）. The results can reflect the most dangerous location, the scope and distribution of failure zone in weak structural plane, and present the process of progressive failure in dam foun- dation as well as the safety coefficient of possible sliding body. These achievements provide an important technical reference for dam foundation treatment measures. The computational results show that ISEM can naturally describe discontinuous de- formation of rock mass such as dislocation, openness and sliding. Besides, this method is characterized by good adaptability, convenient calculation and high compatibility, thus it is regarded as an effective way to make an analysis of anti-sliding stabil- ity of gravity dam
文摘We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
文摘In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.
基金the financial support from the Natural Science Foundation of Hebei Province(Grant No.E2020050012)the National Natural Science Foundation of China(NSFC)(Grant No.52178324)the National Key Research and Development Project of China,the Ministry of Science and Technology of China(Grant No.2018YFC1505504).
文摘Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle induced areas.It is a scenario of strong coupled thermomechanical process involving discontinuity behaviours of rocks.In this context,a numerical model was proposed to investigate the thermal cracking of rocks,in a framework of the continuous-discontinuous element method(CDEM)for efficiently capturing the initiation and propagation of multiple cracks.A simplex integration strategy was adopted to account for the influences of temperature-dependent material properties.Several benchmark tests were considered and the obtained results were compared with analytical solutions and numerical results from the literature.The results show that the fracture degree of the cases when considering temperature-dependent material parameters had 10%differences approximately compared with the cases with constant parameters.
文摘We deal with the numerical solution of the Navier-Stokes equations describing a motion of viscous compressible fluids.In order to obtain a sufficiently stable higher order scheme with respect to the time and space coordinates,we develop a combination of the discontinuous Galerkin finite element(DGFE)method for the space discretization and the backward difference formulae(BDF)for the time discretization.Since the resulting discrete problem leads to a system of nonlinear algebraic equations at each time step,we employ suitable linearizations of inviscid as well as viscous fluxes which give a linear algebraic problem at each time step.Finally,the resulting BDF-DGFE scheme is applied to steady as well as unsteady flows and achieved results are compared with reference data.
基金the National Natural Science Foundation of China (Grant No. 51605212)the Natural Science Foundation of Gansu Province (Grant No. 17JR5RA122)the Project of Hongliu First-class Disciplines Development Program of Lanzhou University of Technology.
文摘The fatigue damage model based on theory of damage mechanics is capable of predicting the fatigue life under multiaxial loading. Meanwhile, the application of critical plane method in the prediction of multiaxial fatigue life has made certain progress. According to the law of thermodynamics, a new damage evolution equation is developed in the present study to predict the fatigue life of geometrically discontinuous structure under tension-torsion loading based on damage mechanics and the critical plane method. The essence of this approach is tha t the st rain parame ter of the uniaxial nonlinear fatigue damage model is replaced with the equivalent strain, which consists of the releva nt parame ters of the critical plane. However, it is difficult to calculate the stress-strain status and the critical plane position of geometrically dis? continuous structure by theoretical methods because of the existence of stress concentration and the multiaxial nonproportional characteristics. Therefore, a new numerical simulation method is proposed to determine the critical plane of geometrically discontinuous structure under multiaxial loading by means of the finite element method and MATLAB software. The fatigue life of notched specimens subjected to combined bending and torsion is predicted using the proposed met hod, and the result is compared with t hose from the experimen ts and the Manson-Cfiffin law. The comparisons show that the proposed method is superior to the Manson-Coffin law and is capable of reproducing the experimental results reasonably when the geometry of the structure is complex. It completely meets the needs of engineering practice.
基金supported by the National Natural Science Foundation of China(No.11975097)the Fundamental Research Funds for the Central Universities(No.2019MS038).
文摘The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in the S N method to automatically optimize the angular distribution and minimize angular discretization errors with lower expenses.The proposed method enables linear dis-continuous finite element quadrature sets over an icosahe-dron to vary their quadrature orders in a one-twentieth sphere so that fine resolutions can be applied to the angular domains that are important.An error estimation that operates in conjunction with the spherical harmonics method is developed to determine the locations where more refinement is required.The adaptive quadrature sets are applied to three duct problems,including the Kobayashi benchmarks and the IRI-TUB research reactor,which emphasize the ability of this method to resolve neutron streaming through ducts with voids.The results indicate that the performance of the adaptive method is more effi-cient than that of uniform quadrature sets for duct transport problems.Our adaptive method offers an appropriate placement of angular unknowns to accurately integrate angular fluxes while reducing the computational costs in terms of unknowns and run times.
基金Supported by the National Natural Science Foundation of China(50976072,51106099,10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)the Science Foundation for the Excellent Youth Scholar of Higher Education of Shanghai(slg09003)~~
文摘A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.
文摘In this paper we investigate the existence, uniqueness and regularity of the solution of semilinear parabolic equations with coefficients that are discontinuous across the interface, some prior estimates are obtained. A net shape of the finite elements around the singular points was designed in [7] to solve the linear elliptic problems, by means of that net, we prove that the approximate solution has the same convergence rate as that without singularity.
文摘This paper deals with the shape and influenced factors of surface non-continuous deformation due to mining. With finite element method, analysis are made to derive the relations between discontinuous deformation and mining affection, weak plane’s position & thickness, and mechanical property of weak-plane medium. The mutual affection of multiple weak-planes is also discussed. The results of the paper lay a foundation for constructing the calculation method of surface discontinuous deformation.
基金funded by the National Natural Science Foundation of China (Grant No.52175111)。
文摘In automotive industries,panel acoustic contribution analysis(PACA)is used to investigate the contributions of the body panels to the acoustic pressure at a certain point of interest.Currently,PACA is implementedmostly by either experiment-based methods or traditional numerical methods.However,these schemes are effort-consuming and inefficient in solving engineering problems,thereby restraining the further development of PACA in automotive acoustics.In this work,we propose a PACA scheme using discontinuous isogeometric boundary element method(IGABEM)to build an easily implementable and efficient method to identify the relative acoustic contributions of each automotive body panel.Discontinuous IGABEMis more accurate and converges faster than continuous BEM and IGABEM in the interior sound pressure evaluation of automotive compartments.In this work,a contribution ratio is defined to estimate the relative acoustic contribution of the structure panels;it can be calculated by reusing the coefficient matrix that has already been generated in the sound pressure evaluation process.The utilization of the parallel technique enables the proposed method to be more efficient than conventional methods;it is validated in two numerical examples,including a car passenger compartment subjected to realistic boundary conditions.A sound pressure response experiment based on a steel box is conducted to verify the accuracy of the interior sound pressure calculation using discontinuous IGABEM.This work is expected to promote the practical process of IGABEM for application in automotive acoustic problems.
