Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under sat...Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit(ADE)scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in nonuniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourthorder finite difference(FD) approximation to the spatial derivatives of the axisymmetric fluid-diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps,giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua(FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%-50% that of FLAC’s basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%-1.8%.展开更多
An explicitly coupled two-dimensional (2D) multiphysics finite element method (FEM) framework comprised of thermal, phase field, mechanical and electromagnetic (TPME) equations was developed to simulate the conversion...An explicitly coupled two-dimensional (2D) multiphysics finite element method (FEM) framework comprised of thermal, phase field, mechanical and electromagnetic (TPME) equations was developed to simulate the conversion of solid kerogen in oil shale to liquid oil through </span><i><span style="font-family:Verdana;font-size:12px;">in-situ</span></i><span style="font-family:Verdana;font-size:12px;"> pyrolysis by radio frequency heating. Radio frequency heating as a method of <i></span><i><span style="font-family:Verdana;font-size:12px;">in-situ</span></i><span style="font-family:Verdana;font-size:12px;"></i> pyrolysis represents a tenable enhanced oil recovery method, whereby an applied electrical potential difference across a target oil shale formation is converted to thermal energy, heating the oil shale and causing it to liquify to become liquid oil. A number of <i></span><i><span style="font-family:Verdana;font-size:12px;">in-situ</span></i><span style="font-family:Verdana;font-size:12px;"></i> pyrolysis methods are reviewed but the focus of this work is on the verification of the TPME numerical framework to model radio frequency heating as a potential dielectric heating process for enhanced oil recovery.</span></span><span style="font-size:10pt;font-family:""> </span><span style="font-family:Verdana;">Very few studies exist which describe production from oil shale;furthermore, there are none that specifically address the verification of numerical models describing radio frequency heating. As a result, the Method of Manufactured Solutions (MMS) was used as an analytical verification method of the developed numerical code. Results show that the multiphysics finite element framework was adequately modeled enabling the simulation of kerogen conversion to oil as a part of the analysis of a TPME numerical model.展开更多
In this work we considered bi-domain partial differential equations(PDEs)with two coupling interface conditions.The one domain is corresponding to the ocean and the second is to the atmosphere.The two coupling conditi...In this work we considered bi-domain partial differential equations(PDEs)with two coupling interface conditions.The one domain is corresponding to the ocean and the second is to the atmosphere.The two coupling conditions are used to linked the interaction between these two regions.As we know that almost every engineering problem modeled via PDEs.The analytical solutions of these kind of problems are not easy,so we use numerical approximations.In this study we discuss the two essential properties,namely mass conservation and stability analysis of two types of coupling interface conditions for the oceanatmosphere model.The coupling conditions arise in general circulation models used in climate simulations.The two coupling conditions are the Dirichlet-Neumann and bulk interface conditions.For the stability analysis,we use the Godunov-Ryabenkii theory of normal-mode analysis.The main empha-sis of this work is to study the numerical properties of coupling conditions and an important point is to maintain conservativity of the overall scheme.Furthermore,for the numerical approximation we use two methods,an explicit and implicit couplings.The implicit coupling have further two algorithms,monolithic algorithm and partitioned iterative algorithm.The partitioned iterative approach is complex as compared to the monolithic approach.In addition,the comparison of the numerical results are exhibited through graphical illustration and simulation results are drafted in tabular form to validate our theoretical investigation.The novel characteristics of the findings from this paper can be of great importance in science and ocean engineering.展开更多
基金the support from the University Transportation Center for Underground Transportation Infrastructure at the Colorado School of Mines for partially funding this research under Grant No. 69A3551747118 of the Fixing America's Surface Transportation Act (FAST Act) of U.S. DoT FY2016
文摘Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit(ADE)scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in nonuniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourthorder finite difference(FD) approximation to the spatial derivatives of the axisymmetric fluid-diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps,giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua(FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%-50% that of FLAC’s basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%-1.8%.
文摘An explicitly coupled two-dimensional (2D) multiphysics finite element method (FEM) framework comprised of thermal, phase field, mechanical and electromagnetic (TPME) equations was developed to simulate the conversion of solid kerogen in oil shale to liquid oil through </span><i><span style="font-family:Verdana;font-size:12px;">in-situ</span></i><span style="font-family:Verdana;font-size:12px;"> pyrolysis by radio frequency heating. Radio frequency heating as a method of <i></span><i><span style="font-family:Verdana;font-size:12px;">in-situ</span></i><span style="font-family:Verdana;font-size:12px;"></i> pyrolysis represents a tenable enhanced oil recovery method, whereby an applied electrical potential difference across a target oil shale formation is converted to thermal energy, heating the oil shale and causing it to liquify to become liquid oil. A number of <i></span><i><span style="font-family:Verdana;font-size:12px;">in-situ</span></i><span style="font-family:Verdana;font-size:12px;"></i> pyrolysis methods are reviewed but the focus of this work is on the verification of the TPME numerical framework to model radio frequency heating as a potential dielectric heating process for enhanced oil recovery.</span></span><span style="font-size:10pt;font-family:""> </span><span style="font-family:Verdana;">Very few studies exist which describe production from oil shale;furthermore, there are none that specifically address the verification of numerical models describing radio frequency heating. As a result, the Method of Manufactured Solutions (MMS) was used as an analytical verification method of the developed numerical code. Results show that the multiphysics finite element framework was adequately modeled enabling the simulation of kerogen conversion to oil as a part of the analysis of a TPME numerical model.
基金the Deans of Scientific Research at King Khalid University,Abha,Saudi Arabia for fund-ing this work through research group program under grant number GRP-216/1443.
文摘In this work we considered bi-domain partial differential equations(PDEs)with two coupling interface conditions.The one domain is corresponding to the ocean and the second is to the atmosphere.The two coupling conditions are used to linked the interaction between these two regions.As we know that almost every engineering problem modeled via PDEs.The analytical solutions of these kind of problems are not easy,so we use numerical approximations.In this study we discuss the two essential properties,namely mass conservation and stability analysis of two types of coupling interface conditions for the oceanatmosphere model.The coupling conditions arise in general circulation models used in climate simulations.The two coupling conditions are the Dirichlet-Neumann and bulk interface conditions.For the stability analysis,we use the Godunov-Ryabenkii theory of normal-mode analysis.The main empha-sis of this work is to study the numerical properties of coupling conditions and an important point is to maintain conservativity of the overall scheme.Furthermore,for the numerical approximation we use two methods,an explicit and implicit couplings.The implicit coupling have further two algorithms,monolithic algorithm and partitioned iterative algorithm.The partitioned iterative approach is complex as compared to the monolithic approach.In addition,the comparison of the numerical results are exhibited through graphical illustration and simulation results are drafted in tabular form to validate our theoretical investigation.The novel characteristics of the findings from this paper can be of great importance in science and ocean engineering.
