An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approxi...An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.展开更多
In recent years, super high-rise buildings (>500 m) are developing very quickly and become an important frontier of civil engineering. The collapse resistance of super high-rise buildings subjected to extremely str...In recent years, super high-rise buildings (>500 m) are developing very quickly and become an important frontier of civil engineering. The collapse resistance of super high-rise buildings subjected to extremely strong earthquake is a critical problem that must be intensively studied. This paper builds up a nonlinear finite element model of the tallest building in China, Shang- hai Tower (632 m), and proposes the modeling method and failure criteria for different structural elements. The dynamic char- acters of this building are then analyzed, and the possible failure modes and collapse processes due to earthquakes are pre- dicted, as well as the corresponding collapse mechanism. This work will be helpful in collapse prevention and the seismic design of super high-rise buildings.展开更多
Thermoelectric effect is the most efficient way to convert electric energy directly from the temperature gradient. Thermoelectric effect-based power generation, cooling and heating devices are solid-stated, environmen...Thermoelectric effect is the most efficient way to convert electric energy directly from the temperature gradient. Thermoelectric effect-based power generation, cooling and heating devices are solid-stated, environmentally friendly, reliable, long-lived, easily maintainable, and easy to achieve miniaturization and integration. So they have unparalleled advantages in the aerospace, vehicle industry, waste heat recovery, electronic cooling, etc. This paper reviews the progress in thermodynamic analyses and optimizations for single- and multiple-element, single- and multiple-stage, and combined thermoelectric generators, thermoelectric refrigerators and thermoelectric heat pumps, especially in the aspects of non-equilibrium thermodynamics and finite time thermodynamics. It also discusses the developing trends of thermoelectric devices, such as the heat sources of thermoelectric generators, multi-stage thermoelectric devices, combined thermoelectric devices, and heat transfer enhancement of thermoelectric devices.展开更多
Regular assumption of finite element meshes is a basic condition of most analysis of finite element approximations both for conventional conforming elements and nonconforming elements. The aim of this paper is to pres...Regular assumption of finite element meshes is a basic condition of most analysis of finite element approximations both for conventional conforming elements and nonconforming elements. The aim of this paper is to present a novel approach of dealing with the approximation of a four-degree nonconforming finite element for the second order elliptic problems on the anisotropic meshes. The optimal error estimates of energy norm and L^2-norm without the regular assumption or quasi-uniform assumption are obtained based on some new special features of this element discovered herein. Numerical results are given to demonstrate validity of our theoretical analysis.展开更多
The liquid-solid phase change process of a simple one-dimensional slab is studied in this paper.By taking entransy dissipation minimization as optimization objective,the optimal external reservoir temperature profiles...The liquid-solid phase change process of a simple one-dimensional slab is studied in this paper.By taking entransy dissipation minimization as optimization objective,the optimal external reservoir temperature profiles are derived by using optimal control theory under the condition of a fixed freezing or melting time.The entransy dissipation corresponding to the optimal heat exchange strategies of minimum entransy dissipation is 8/9 of that corresponding to constant reservoir temperature operations,which is independent of all system parameters.The obtained results for entransy dissipation minimization are also compared with those obtained for the optimal heat exchange strategies of minimum entropy generation and constant reservoir temperature operations by numerical examples.The obtained results can provide some theoretical guidelines for the choice of optimal cooling or heating strategy in practical liquid-solid phase change processes.展开更多
As one kind of key anti-fatigue manufacture approaches with simplicity and effectiveness, the hole cold expansion technology satisfies the increasing needs for light weight and durability of aircraft structures.It can...As one kind of key anti-fatigue manufacture approaches with simplicity and effectiveness, the hole cold expansion technology satisfies the increasing needs for light weight and durability of aircraft structures.It can improve the fatigue life by several times at no additional weight conditions.The hole cold expansion technology has been widely used in manufacturing and repairing of both fighters and commercial aircraft, and has become a research hotspot in the strengthening technology.In recent years, hole cold expansion process methods, residual stress around expanded holes, the behavior of fatigue crack initiation and propagation, and fatigue lives after cold expansion are researched extensively through lots of experiments and finite element simulations.A review on the hole cold expansion technology research status in the last twenty years is presented in this paper.Via the analysis of the current characteristics and defects of the hole cold expansion technology, combined with the actual needs in design and manufacture of new-generation aircraft, development trends and novel research directions are presented for realizing precise and high-efficiency anti-fatigue manufacture.展开更多
Underground utility tunnels are widely used in urban areas throughout the world for lifeline networks due to their easy maintenance and environmental protection capabilities. However, knowledge about their seismic per...Underground utility tunnels are widely used in urban areas throughout the world for lifeline networks due to their easy maintenance and environmental protection capabilities. However, knowledge about their seismic performance is still quite limited and seismic design procedures are not included in current design codes. This paper describes a series of shaking table tests the authors performed on a scaled utility tunnel model to explore its performance under earthquake excitation. Details of the experimental setup are first presented focusing on aspects such as the design of the soil container, scaled structural model, sensor array arrangement and test procedure. The main observations from the test program, including structural response, soil response, soil-structure interaction and earth pressure, are summarized and discussed. Further, a finite element model （FEM） of the test utility tunnel is established where the nonlinear soil properties are modeled by the Drucker- Prager constitutive model; the master-slave surface mechanism is employed to simulate the soil-structure dynamic interaction; and the confining effect of the laminar shear box to soil is considered by proper boundary modeling. The results from the numerical model are compared with experiment measurements in terms of displacement, acceleration and amplification factor of the structural model and the soil. The comparison shows that the numerical results match the experimental measurements quite well. The validated numerical model can be adopted for further analysis.展开更多
The mass entransy and its dissipation extremum principle have opened up a new direction for the mass transfer optimization. Firstly, the emergence and development process of both the mass entransy and its dissipalion ...The mass entransy and its dissipation extremum principle have opened up a new direction for the mass transfer optimization. Firstly, the emergence and development process of both the mass entransy and its dissipalion extremum principle are reviewed. Secondly, the combination of the mass entransy dissipation extremum principle and the finite-time thermodynamics for opti- mizing the mass transfer processes of one-way isothermal mass transfer, two-way isothermal equimolar mass transfer, and iso- thermal throttling and isothermal crystallization are summarized. Thirdly, the combination of the mass entransy dissipation ex- tremum principle and the constructal theory for optimizing the mass transfer processes of disc-to-point and volume-to-point problems are summarized. The scientific features of the mass entransy dissipation extremam principle are emphasized.展开更多
This paper is devoted to the five parameters nonconforming finite element schemes with moving grids for velocity-pressure mixed formulations of the nonstationary Stokes problem in 2-D. We show that this element has an...This paper is devoted to the five parameters nonconforming finite element schemes with moving grids for velocity-pressure mixed formulations of the nonstationary Stokes problem in 2-D. We show that this element has anisotropic behavior and derive anisotropic error estimations in some certain norms of the velocity and the pressure based on some novel techniques. Especially through careful analysis we get an interesting result on consistency error estimation, which has never been seen for mixed finite element methods in the previously literatures.展开更多
提出一种基于确定的有穷状态自动机(deterministic finite automaton,简称DFA)的正则表达式压缩算法.首先,定义了膨胀率DR(distending rate)来描述正则表达式的膨胀特性.然后基于DR提出一种分片的算法RECCADR(regular expressions cut a...提出一种基于确定的有穷状态自动机(deterministic finite automaton,简称DFA)的正则表达式压缩算法.首先,定义了膨胀率DR(distending rate)来描述正则表达式的膨胀特性.然后基于DR提出一种分片的算法RECCADR(regular expressions cut and combine algorithm based on DR),有效地选择出导致DFA状态膨胀的片段并隔离,降低了单个正则表达式存储需求.同时,基于正则表达式的组合关系提出一种选择性分群算法REGADR(regular expressions group algorithm based on DR),在可以接受的存储需求总量下,通过选择性分群大幅度减少了状态机的个数,有效地降低了匹配算法的复杂性.展开更多
This paper deals with the high accuracy analysis of bilinear finite element on the class of anisotropic rectangular meshes. The inverse inequalities on anisotropic meshes are established. The superclose and the superc...This paper deals with the high accuracy analysis of bilinear finite element on the class of anisotropic rectangular meshes. The inverse inequalities on anisotropic meshes are established. The superclose and the superconvergence are obtained for the second order elliptic problem. A numerical test is given, which coincides with our theoretical analysis.展开更多
The research of rock properties based on its inherent microscopic to mesoscopic porous structure has drawn great attention for its potential in predicting the macroscopic behavior of rocks.An accurate reconstruction o...The research of rock properties based on its inherent microscopic to mesoscopic porous structure has drawn great attention for its potential in predicting the macroscopic behavior of rocks.An accurate reconstruction of the threedimensional porous structure is a premise for the related studies of hydraulic and mechanical properties of rocks,such as the transport properties and mechanical responses under pressures.In this paper,we present a computer procedure for reconstructing the 3D porous structure of low-permeability sandstone.Two large-size 3D models are reconstructed based on the information of a reference model which is established from computed tomography(CT)images.A self-developed finite element method is applied to analyze the nonlinear mechanical behavior of the sandstone based on its reconstructed model and to compare the results with those based on the reference model.The good consistency of the obtained mechanical responses indicates the potential of using reconstruction models to predict the influences of porous structure on the mechanical properties of low-permeability sandstone.展开更多
Viscous fluid model and potential flow model with and without artificial damping force(f=-μV,μ the damping coefficient and V the local averaging flow velocity) are employed in this work to investigate the phenomenon...Viscous fluid model and potential flow model with and without artificial damping force(f=-μV,μ the damping coefficient and V the local averaging flow velocity) are employed in this work to investigate the phenomenon of fluid resonance in narrow gaps between multi-bodies in close proximity under water waves.The numerical results are compared with experimental data available in the literature.The comparison demonstrates that both the viscous fluid model and the potential flow model are able to predict the resonant frequency reasonably well.However the conventional potential flow model(without artificial damping term) significantly over-predicts the wave height in narrow gaps around the resonant frequency.In order to calibrate the appropriate damping coefficient used for the potential model and make it work as well as the viscous fluid model in predicting the resonant wave height in narrow gaps but with little computational efforts,the dependence of damping coefficient μ on the body geometric dimensions is examined considering the parameters of gap width Bg,body draft D,body breadth ratio Br and body number n(n = 2,3),where Br = BB/BA for the case of two bodies(Body A and Body B) with different breadths of BA and BB,respectively.It was confirmed that the damping coefficient used for the potential flow model is not sensitive to the geometric dimensions and spatial arrangement.It was found that μ∈ [0.4,0.5] may guarantee the variation of Hg/H0 with kh to be generally in good agreement with the experimental data and the results of viscous fluid model,where Hg is the excited wave height in narrow gaps under various dimensionless incident wave frequencies kh,H0 is the incident wave height,k = 2π/L is the wave number and h is the water depth.展开更多
The pore characteristics,mineral compositions,physical and mechanical properties of the subarkose sandstones were acquired by means of CT scan,X-ray diffraction and physical tests.A few physical models possessing the ...The pore characteristics,mineral compositions,physical and mechanical properties of the subarkose sandstones were acquired by means of CT scan,X-ray diffraction and physical tests.A few physical models possessing the same pore characteristics and matrix properties but different porosities compared to the natural sandstones were developed.The 3D finite element models of the rock media with varied porosities were established based on the CT image processing of the physical models and the MIMICS software platform.The failure processes of the porous rock media loaded by the split Hopkinson pressure bar(SHPB) were simulated by satisfying the elastic wave propagation theory.The dynamic responses,stress transition,deformation and failure mechanisms of the porous rock media subjected to the wave stresses were analyzed.It is shown that an explicit and quantitative analysis of the stress,strain and deformation and failure mechanisms of porous rocks under the wave stresses can be achieved by using the developed 3D finite element models.With applied wave stresses of certain amplitude and velocity,no evident pore deformation was observed for the rock media with a porosity less than 15%.The deformation is dominantly the combination of microplasticity(shear strain),cracking(tensile strain) of matrix and coalescence of the cracked regions around pores.Shear stresses lead to microplasticity,while tensile stresses result in cracking of the matrix.Cracking and coalescence of the matrix elements in the neighborhood of pores resulted from the high transverse tensile stress or tensile strain which exceeded the threshold values.The simulation results of stress wave propagation,deformation and failure mechanisms and energy dissipation in porous rock media were in good agreement with the physical tests.The present study provides a reference for analyzing the intrinsic mechanisms of the complex dynamic response,stress transit mode,deformation and failure mechanisms and the disaster mechanisms of rock media.展开更多
Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient applicat...Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.展开更多
The extended finite element method (XFEM) is a new numerical method for modeling discontinuity. Research about numerical modeling for concrete hydraulic fracturing by XFEM is explored. By building the virtual work pri...The extended finite element method (XFEM) is a new numerical method for modeling discontinuity. Research about numerical modeling for concrete hydraulic fracturing by XFEM is explored. By building the virtual work principle of the fracture problem considering water pressure on the crack surface, the governing equations of XFEM for hydraulic fracture modeling are derived. Implementation of the XFEM for hydraulic fracturing is presented. Finally, the method is verified by two examples and the advan- tages of the XFEM for hydraulic fracturing analysis are displayed.展开更多
An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same t...An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.展开更多
A nonconforming H^1-Calerkin mixed finite element method is analyzed for Sobolev equations on anisotropic meshes. The error estimates are obtained without using Ritz-Volterra projection.
The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space ...The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).展开更多
基金supported by the National Natural Science Foundation of China No.10671184
文摘An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.
基金supported by the National Natural Science Foundation of China (Grant No. 90815025)the Tsinghua University Research Funds (Grant No. 2010THZ02-1)the "Program for New Century Excellent Talents in University"
文摘In recent years, super high-rise buildings (>500 m) are developing very quickly and become an important frontier of civil engineering. The collapse resistance of super high-rise buildings subjected to extremely strong earthquake is a critical problem that must be intensively studied. This paper builds up a nonlinear finite element model of the tallest building in China, Shang- hai Tower (632 m), and proposes the modeling method and failure criteria for different structural elements. The dynamic char- acters of this building are then analyzed, and the possible failure modes and collapse processes due to earthquakes are pre- dicted, as well as the corresponding collapse mechanism. This work will be helpful in collapse prevention and the seismic design of super high-rise buildings.
基金supported by the National Natural Science Foundation of China(Grant Nos.11305266&51576207)the National Basic Research Program of China("973"Project)(Grant No.2012CB720405)
文摘Thermoelectric effect is the most efficient way to convert electric energy directly from the temperature gradient. Thermoelectric effect-based power generation, cooling and heating devices are solid-stated, environmentally friendly, reliable, long-lived, easily maintainable, and easy to achieve miniaturization and integration. So they have unparalleled advantages in the aerospace, vehicle industry, waste heat recovery, electronic cooling, etc. This paper reviews the progress in thermodynamic analyses and optimizations for single- and multiple-element, single- and multiple-stage, and combined thermoelectric generators, thermoelectric refrigerators and thermoelectric heat pumps, especially in the aspects of non-equilibrium thermodynamics and finite time thermodynamics. It also discusses the developing trends of thermoelectric devices, such as the heat sources of thermoelectric generators, multi-stage thermoelectric devices, combined thermoelectric devices, and heat transfer enhancement of thermoelectric devices.
文摘Regular assumption of finite element meshes is a basic condition of most analysis of finite element approximations both for conventional conforming elements and nonconforming elements. The aim of this paper is to present a novel approach of dealing with the approximation of a four-degree nonconforming finite element for the second order elliptic problems on the anisotropic meshes. The optimal error estimates of energy norm and L^2-norm without the regular assumption or quasi-uniform assumption are obtained based on some new special features of this element discovered herein. Numerical results are given to demonstrate validity of our theoretical analysis.
基金supported by the Program for New Century Excellent Talents in Universities of China (Grant No 20041006)the Foundation for the Authors of National Excellent Doctoral Dissertation of China (Grant No 200136)
文摘The liquid-solid phase change process of a simple one-dimensional slab is studied in this paper.By taking entransy dissipation minimization as optimization objective,the optimal external reservoir temperature profiles are derived by using optimal control theory under the condition of a fixed freezing or melting time.The entransy dissipation corresponding to the optimal heat exchange strategies of minimum entransy dissipation is 8/9 of that corresponding to constant reservoir temperature operations,which is independent of all system parameters.The obtained results for entransy dissipation minimization are also compared with those obtained for the optimal heat exchange strategies of minimum entropy generation and constant reservoir temperature operations by numerical examples.The obtained results can provide some theoretical guidelines for the choice of optimal cooling or heating strategy in practical liquid-solid phase change processes.
基金the Fund of National Engineering and Research Center for Commercial Aircraft Manufacturing of China (No.SAMC12-JS-15-021)the Funding of Jiangsu Innovation Program for Graduate Education of China (No.CXLX12_0137)the Fundamental Research Funds for the Central Universities of China
文摘As one kind of key anti-fatigue manufacture approaches with simplicity and effectiveness, the hole cold expansion technology satisfies the increasing needs for light weight and durability of aircraft structures.It can improve the fatigue life by several times at no additional weight conditions.The hole cold expansion technology has been widely used in manufacturing and repairing of both fighters and commercial aircraft, and has become a research hotspot in the strengthening technology.In recent years, hole cold expansion process methods, residual stress around expanded holes, the behavior of fatigue crack initiation and propagation, and fatigue lives after cold expansion are researched extensively through lots of experiments and finite element simulations.A review on the hole cold expansion technology research status in the last twenty years is presented in this paper.Via the analysis of the current characteristics and defects of the hole cold expansion technology, combined with the actual needs in design and manufacture of new-generation aircraft, development trends and novel research directions are presented for realizing precise and high-efficiency anti-fatigue manufacture.
基金Key Project in the National Science & Technology Pillar Program Under Grant No. 2006BAJ03B03Research Fund for Young Teacher Supported by State Key Laboratory for Disaster Reduction in Civil Engineering Under Grant No. SLDRCE08-C-03
文摘Underground utility tunnels are widely used in urban areas throughout the world for lifeline networks due to their easy maintenance and environmental protection capabilities. However, knowledge about their seismic performance is still quite limited and seismic design procedures are not included in current design codes. This paper describes a series of shaking table tests the authors performed on a scaled utility tunnel model to explore its performance under earthquake excitation. Details of the experimental setup are first presented focusing on aspects such as the design of the soil container, scaled structural model, sensor array arrangement and test procedure. The main observations from the test program, including structural response, soil response, soil-structure interaction and earth pressure, are summarized and discussed. Further, a finite element model （FEM） of the test utility tunnel is established where the nonlinear soil properties are modeled by the Drucker- Prager constitutive model; the master-slave surface mechanism is employed to simulate the soil-structure dynamic interaction; and the confining effect of the laminar shear box to soil is considered by proper boundary modeling. The results from the numerical model are compared with experiment measurements in terms of displacement, acceleration and amplification factor of the structural model and the soil. The comparison shows that the numerical results match the experimental measurements quite well. The validated numerical model can be adopted for further analysis.
基金supported by the National Natural Science Foundation China(Grant Nos.51176203 and 10905093)
文摘The mass entransy and its dissipation extremum principle have opened up a new direction for the mass transfer optimization. Firstly, the emergence and development process of both the mass entransy and its dissipalion extremum principle are reviewed. Secondly, the combination of the mass entransy dissipation extremum principle and the finite-time thermodynamics for opti- mizing the mass transfer processes of one-way isothermal mass transfer, two-way isothermal equimolar mass transfer, and iso- thermal throttling and isothermal crystallization are summarized. Thirdly, the combination of the mass entransy dissipation ex- tremum principle and the constructal theory for optimizing the mass transfer processes of disc-to-point and volume-to-point problems are summarized. The scientific features of the mass entransy dissipation extremam principle are emphasized.
基金This research is supported by the National Science Foundation of China（No.10371113）.The authors would like to thank the anonymous referees for their helpful suggestions.
文摘This paper is devoted to the five parameters nonconforming finite element schemes with moving grids for velocity-pressure mixed formulations of the nonstationary Stokes problem in 2-D. We show that this element has anisotropic behavior and derive anisotropic error estimations in some certain norms of the velocity and the pressure based on some novel techniques. Especially through careful analysis we get an interesting result on consistency error estimation, which has never been seen for mixed finite element methods in the previously literatures.
文摘提出一种基于确定的有穷状态自动机(deterministic finite automaton,简称DFA)的正则表达式压缩算法.首先,定义了膨胀率DR(distending rate)来描述正则表达式的膨胀特性.然后基于DR提出一种分片的算法RECCADR(regular expressions cut and combine algorithm based on DR),有效地选择出导致DFA状态膨胀的片段并隔离,降低了单个正则表达式存储需求.同时,基于正则表达式的组合关系提出一种选择性分群算法REGADR(regular expressions group algorithm based on DR),在可以接受的存储需求总量下,通过选择性分群大幅度减少了状态机的个数,有效地降低了匹配算法的复杂性.
基金This research is supported by the National Natural Science Foundation of China(No.10371113) Foundation of Oversea Scholar of China(No.2001(119)) the Project of Creative Engineering of Henan Province of China 2002(219) NSF of Henan Province of China.
文摘This paper deals with the high accuracy analysis of bilinear finite element on the class of anisotropic rectangular meshes. The inverse inequalities on anisotropic meshes are established. The superclose and the superconvergence are obtained for the second order elliptic problem. A numerical test is given, which coincides with our theoretical analysis.
基金the financial supports from the National Natural Science Foundation for Distinguished Young Scholars of China(51125017)the National Natural Science Foundation of China(51374213)the National Basic Research Pro-gram of China(Grant 2010CB226804,2011CB201201).
文摘The research of rock properties based on its inherent microscopic to mesoscopic porous structure has drawn great attention for its potential in predicting the macroscopic behavior of rocks.An accurate reconstruction of the threedimensional porous structure is a premise for the related studies of hydraulic and mechanical properties of rocks,such as the transport properties and mechanical responses under pressures.In this paper,we present a computer procedure for reconstructing the 3D porous structure of low-permeability sandstone.Two large-size 3D models are reconstructed based on the information of a reference model which is established from computed tomography(CT)images.A self-developed finite element method is applied to analyze the nonlinear mechanical behavior of the sandstone based on its reconstructed model and to compare the results with those based on the reference model.The good consistency of the obtained mechanical responses indicates the potential of using reconstruction models to predict the influences of porous structure on the mechanical properties of low-permeability sandstone.
基金supports from the Natural National Science Foundation of China (Grant Nos.50909016,50921001 and 10802014)support of ARC Discovery Project Program (Grant No. DP0557060)supported by the Open Fund from the State Key Laboratory of Structural Analysis for Industrial Equipment (Grant No. GZ0909)
文摘Viscous fluid model and potential flow model with and without artificial damping force(f=-μV,μ the damping coefficient and V the local averaging flow velocity) are employed in this work to investigate the phenomenon of fluid resonance in narrow gaps between multi-bodies in close proximity under water waves.The numerical results are compared with experimental data available in the literature.The comparison demonstrates that both the viscous fluid model and the potential flow model are able to predict the resonant frequency reasonably well.However the conventional potential flow model(without artificial damping term) significantly over-predicts the wave height in narrow gaps around the resonant frequency.In order to calibrate the appropriate damping coefficient used for the potential model and make it work as well as the viscous fluid model in predicting the resonant wave height in narrow gaps but with little computational efforts,the dependence of damping coefficient μ on the body geometric dimensions is examined considering the parameters of gap width Bg,body draft D,body breadth ratio Br and body number n(n = 2,3),where Br = BB/BA for the case of two bodies(Body A and Body B) with different breadths of BA and BB,respectively.It was confirmed that the damping coefficient used for the potential flow model is not sensitive to the geometric dimensions and spatial arrangement.It was found that μ∈ [0.4,0.5] may guarantee the variation of Hg/H0 with kh to be generally in good agreement with the experimental data and the results of viscous fluid model,where Hg is the excited wave height in narrow gaps under various dimensionless incident wave frequencies kh,H0 is the incident wave height,k = 2π/L is the wave number and h is the water depth.
基金supported by the National Basic Research Program of China ("973" Project) (Grant Nos 2010CB226804, 2002CB412705)the National Natural Science Foundation of China (Grant No 50974125)+1 种基金the Natural Sciences & Engineering Research Council of Canada (PGS-D2-2006)Beijing Key Laboratory Projects
文摘The pore characteristics,mineral compositions,physical and mechanical properties of the subarkose sandstones were acquired by means of CT scan,X-ray diffraction and physical tests.A few physical models possessing the same pore characteristics and matrix properties but different porosities compared to the natural sandstones were developed.The 3D finite element models of the rock media with varied porosities were established based on the CT image processing of the physical models and the MIMICS software platform.The failure processes of the porous rock media loaded by the split Hopkinson pressure bar(SHPB) were simulated by satisfying the elastic wave propagation theory.The dynamic responses,stress transition,deformation and failure mechanisms of the porous rock media subjected to the wave stresses were analyzed.It is shown that an explicit and quantitative analysis of the stress,strain and deformation and failure mechanisms of porous rocks under the wave stresses can be achieved by using the developed 3D finite element models.With applied wave stresses of certain amplitude and velocity,no evident pore deformation was observed for the rock media with a porosity less than 15%.The deformation is dominantly the combination of microplasticity(shear strain),cracking(tensile strain) of matrix and coalescence of the cracked regions around pores.Shear stresses lead to microplasticity,while tensile stresses result in cracking of the matrix.Cracking and coalescence of the matrix elements in the neighborhood of pores resulted from the high transverse tensile stress or tensile strain which exceeded the threshold values.The simulation results of stress wave propagation,deformation and failure mechanisms and energy dissipation in porous rock media were in good agreement with the physical tests.The present study provides a reference for analyzing the intrinsic mechanisms of the complex dynamic response,stress transit mode,deformation and failure mechanisms and the disaster mechanisms of rock media.
文摘Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.
基金Supported by the National Natural Science Foundation of China (Grant Nos. 50539030, 50609004)the National Basic Research Program of China ("973" Program) (Grant No. 2007CB714104)
文摘The extended finite element method (XFEM) is a new numerical method for modeling discontinuity. Research about numerical modeling for concrete hydraulic fracturing by XFEM is explored. By building the virtual work principle of the fracture problem considering water pressure on the crack surface, the governing equations of XFEM for hydraulic fracture modeling are derived. Implementation of the XFEM for hydraulic fracturing is presented. Finally, the method is verified by two examples and the advan- tages of the XFEM for hydraulic fracturing analysis are displayed.
基金Supported by the National Natural Science Foundation of China (10601022)Natural Science Foundation of Inner Mongolia Autonomous Region (200607010106)Youth Science Foundation of Inner Mongolia University(ND0702)
文摘An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.
基金Supported by the National Natural Science Foundation of China(No.10671184).
文摘A nonconforming H^1-Calerkin mixed finite element method is analyzed for Sobolev equations on anisotropic meshes. The error estimates are obtained without using Ritz-Volterra projection.
文摘The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).
