2D Coupled 3D:A New Numerical Model for Dual - Structured - Aquifer System

This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.

3D Finite Element Analysis of TBM Water Diversion Tunnel Segment Coupled with Seepage Field

In most studies of tunnel boring machine(TBM)tunnelling, the groundwater pressure was not considered, or was simplified and exerted on the boundary of lining structure. Meanwhile, the leakage, which mainly occurs in the segment joints, was often ignored in the relevant studies of TBM tunnelling. Additionally, the geological models in these studies were simplified to different extents, and mostly were simplified as homogenous bodies. Considering the deficiencies above, a 3D refined model of the surrounding rock of a tunnel is firstly established using NURBS-TIN-BRe P hybrid data structure in this paper. Then the seepage field of the surrounding rock considering the leakage in the segment joints is simulated. Finally, the stability of TBM water diversion tunnel is studied coupled with the seepage simulation, to analyze the stress-strain conditions, the axial force and the bending moment of tunnel segment considering the leakage in the segment joints. The results illustrate that the maximum radial displacement, the minimum principal stress, the maximum principal stress and the axial force of segment lining considering the seepage effect are all larger than those disregarding the seepage effect.

Coupled thermo-hydro-mechanical simulation of CO2 enhanced gas recovery with an extended equation of state module for TOUGH2MP-FLAC3D

As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.

Hydrodynamic Response of A Fully Coupled TLP Hull-TTR System with Detailed Modeling of A Hydraulic Pneumatic Tensioner and Riser Joints

Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull.Owing to the complicated configuration of the tensioners,the hull and TTRs form a strong coupled system.Traditionally,some simplified tensioner models are applied to analyze the TLP structures.There is a large discrepancy between their analysis results and the actual mechanism behaviors of a tensioner.It is very necessary to develop a more detailed tensioner model to consider the coupling effects between TLP and TTRs.In the present study,a fully coupled TLP hull-TTR system for hydrodynamic numerical simulation is established.A specific hydraulic pneumatic tensioner is modeled by considering 4 cylinders.The production TTR model is stacked up by specific riser joints.The simulation is also extended to analyze an array of TTRs.Different regular and irregular waves are considered.The behaviors of different cylinders are presented.The results show that it is important to consider the specific configurations of the tensioner and TTRs,which may lead to obviously different response behaviors,compared with those from a simplified model.

3D numerical simulation in acoustic-elastic coupled media with staggered-grid finite-difference method

Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.

Shear behavior of intact granite under thermo-mechanical coupling and three-dimensional morphology of shear-formed fractures

The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.

Micro Mechanical Model of 3D Woven Composites

A combined beam model representing the periodicity of the microstructure and micro deformation of 3D woven composites is developed for predicting mechanical properties. The model considers the effects of off axial tension/compression and bending/shearing couplings as well as the mutual reactions of fiber yarns. The method determining microstructure by using woven parameters is described for a typical 3D woven composite material. An analytical cell, constructed by a minimum periodic section of yarn and interlayer matrix, is adopted. Micro stresses in the cell under in-plane tensile loading are obtained by using the proposed beam model and macro modulus is then obtained by the averaging method. Material tests and a 2D micro FEM analysis are made to evaluate this model. Analyses reveal that micro stress caused by tensile/bending coupling effect is not negligible in the stress analysis.

Bicharacteristic Curves in Modeling of the Earth Crust and Upper Mantle

The aim of this paper is the introduction of a new approach to 3D modelling of elastic piecewise homogeneous media, in particular Earth crust and upper Mantle. The method is based on the principle of tomography with Earthquake as a source of the signal and receiver stations on the surface. The wave propagation in solid media is described by a system of three strongly coupled hyperbolic equations with piece - wise constant coefitients. The characteristic set and hi-characteristic curves of this system are computed in a homogeneous half-space with free boundary and the formulae of reflection and diffraction of the hi-characteristics on the internal boundaries of the media. Applications of the characteristic set and bi-eharacteristic curves for the inverse problem in geophysics and Earth modelling are given.

Coupled 3D discrete-continuum numerical modeling of pile penetration in sand

A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numerical simulation program was used which considered sand near the pile as interacted particles using a discrete element method;the sand away from the pile was simulated as a continuous medium exhibiting linear elastic behaviors.The domain analyzed was divided into two zones.Contact forces at the interface between the two zones were obtained from a discrete zone and applied to the continuum boundaries as nodal forces,while the interface velocities were obtained from the continuum zone and applied to the discrete boundaries.We show that the coupled discrete-continuum simulation can give a microscopic description of the pile penetration process without losing the discrete nature of the zone concerned,and may significantly improve computational efficiency.

Study on coupled 3D seepage and stress fields of the complex channel project

Seepage is a vital reason that may bring many adverse consequences such as subsidence,inclination and fracture to the channel,and is harmful to the safety and stability of the channel.Thus,a 3D visualization model is established for the engineering geological information of the research channel section by using NURBS-TIN-BRep hybrid data structure.Coupled with the VOF(volume of fluid)method,the N-S(Navier-Stokes)equations are applied to seepage simulation of the research channel section.Then the stability of the channel is studied coupled with the seepage simulation results,to comprehensively analyze the stress and displacement conditions of the channel under the impact of different factors such as seepage and underground goafs.The results of this study illustrate that the channel seepage has great influence on its stability,especially on the displacement field:it will lead to a significant sedimentation to the foundation.Therefore,during the practical construction,it is suggested that the certain part of the channel should be reinforced and effective seepage control measures should be taken.

A three-dimensional coupled numerical model of nonlinear waves in a harbor

A 3-D time-domain numerical coupled model for nonlinear waves acting on a ship in a harbor has been developed in the present study.The whole domain is divided into the inner domain and the outer domain.The inner domain is the area around the ship,where the flow is expressed by the Laplace equation and numerically solved by the finite element method.The other area is the outer domain,where the flow is described by the higher-order Boussinesq equations and numerically solved by the finite difference method.The matching conditions on the interfaces between the inner domain and the outer domain,the procedure of coupled solution,the length of common domain and the mesh generation in the inner domain are discussed in detail.The other coupled model with the flow in the inner domain governed by the simplified linear Euler equations and relevant physical experiment are adopted to validate the present coupled model,and it is shown that the numerical results of the present model agree with the experimental data,so the present model can be used for the study on the effect of nonlinear waves acting on a fixed ship in a large area and provide a reference for the time-domain simulation of nonlinear wave forces on an arbitrary object in a large harbor and the 3-D district computation in the future.

PEM电解槽流场的多物理场耦合建模与分析

基于传热、传质和电化学动力学等多物理场耦合理论,建立PEMWE单电池三维模型,研究流场结构对PEMWE性能的影响。对比分析了4种流场结构电解槽内流体流速、物质分布和极化曲线的情况。研究结果表明:在相同的操作条件下,叶栅形流场PEMWE的槽内物质流速和压力分布更加均匀,电化学反应更加充分,而传统的平行流场结构电化学性能较差;在80℃、0.1 MPa的工作条件下,叶栅形流场PEMWE的极化电压更低,当电流密度为2 A/cm^(2)时,单电池极化电压仅为1.686 V,表明采用叶栅形流场结构相较其余几种方式能更有效地提升PEMWE单电池的性能。

SIMULATION OF HYDRAULIC TRANSIENTS IN HYDROPOWER SYSTEMS USING THE 1-D-3-D COUPLING APPROACH

Although the hydraulic transients in pipe systems are usually simulated by using a one-dimensional （l-D） approach, local three-dimensional （3-D） simulations are necessary because of obvious 3-D flow features in some local regions of the hydropower systems. This paper combines the 1-D method with a 3-D fluid flow model to simulate the Multi-Dimensional （MD） hydraulic transients in hydropower systems and proposes two methods for modeling the compressible water with the correct wave speed, and two strategies for efficiently coupling the 1-D and 3-D computational domains. The methods are validated by simulating the water hammer waves and the oscillations of the water level in a surge tank, and comparing the results ~with the 1-D solution data. An MD study is conducted for the transient flows in a realistic water conveying system that consists of a draft tube, a tailrace surge tank and a tailrace tunnel. It is shown that the 1-D-3-D coupling approach is an efficient and promising way to simulate the hydraulic transients in the hydropower systems in which the interactions between 1-D hydraulic fluctuations of the pipeline systems and the local 3-D flow patterns should be considered.

地表河网-地下水流系统耦合模拟Ⅰ:模型

本文通过采用联立求解地表河流一维明渠非恒定渐变流与地下水拟三维非稳定流运动方程的技术路线,提出了地表河网与地下水流耦合模拟算法。该算法中地表河网流动分别由质量和动量守恒定律而建立的连续性方程和运动方程所组成;地下水流采用多层"拟三维"非稳定流模型,通过彼此间的水力联系,建立了地表河网-地下水流耦合模拟评价模型,该模型可用于地表水与地下水转化关系较为密切地区的水资源综合评价预测中。

车辆-轨道耦合系统中基于变参数的三维图形仿真研究

建立一个随参数改变而改变的三维可视模型 ,是机车车辆 -轨道耦合系统动力学可视化的一个急待解决的重要问题 .文中将可视化、三维建模技术引入机车车辆 -轨道耦合动力学仿真中 ,系统地研究了机车车辆轨道耦合系统的基于变参数的三维图形仿真的方法 .综合建立了车辆模型的动态装配三维图形系统 ,直、曲线轨道的的三维离散结构等 ,为机车车辆 -轨道耦合动力学仿真建立了一个有效的可视化综合仿真平台 .

机车车辆轨道耦合几何结构基于面向对象的建模分析研究

将可视化、三维建模等技术引入机车车辆轨道耦合动力学仿真中 ,从系统工程的观点 ,在机车车辆—轨道耦合动力学计算核心模块的基础上 ,提出了基于面向对象的机车车辆—轨道耦合动力学的计算机仿真方法及相应的模型解释和标识机制 ,成功地形成了一个以三维可视输入输出 ,并与核心模块有机结合的机车车辆轨道动力学分析综合仿真平台。

河流-含水层系统数值模拟方法探讨

地表水与地下水转化量的模拟预测一直是区域水资源评价和管理中的难题。提出了一种基于河流越流系数的处理方法,即将转化量计算模型同三维地下水数值模型完全耦合,能够预测在有人工干预条件下,地表水与地下水转化量的变化趋势。提出的方法适合多含水层系统的情形。实例研究表明,提出的河流-含水层系统模拟方法原理简单,易于实现,可以推广应用。

基于组合式场路耦合法的多绕组变压器建模与阻抗参数设计

为解决低励磁阻抗多绕组变压器建模、多绕组变压器复合短路阻抗及多并联支路间电流分配计算的问题,基于组合式场路耦合分析,给出一种高效、通用的多绕组变压器建模和阻抗参数设计方法。根据初步设计的变压器结构,采用磁场能量法,二次开发出基于ANSYS参数化设计语言(ANSYS parametric design language,APDL)的三维静态磁场专用分析程序,计算出多绕组电感矩阵,并以此为依据构建多绕组变压器的线性仿真模型,进行相关仿真研究。通过仿真结果反馈的信息适当修改变压器结构以满足要求。实验结果表明,所提方法能够准确计算出多并联支路间电流分配,各种短路工况下复合短路阻抗计算结果可满足工程实践的需要。

珠江三角洲一维河网与三维河口耦合水质模型模拟与验证

在一维与三维耦合水动力、悬沙模型的基础上,构建了珠江三角洲一维与三维耦合水质模型。模型中包括浮游植物-氮-磷-碳-溶解氧循环,并考虑颗粒物在底泥与水体界面的再悬浮过程。选择了三个代表性时段(1999年1月、1999年7月与2000年)进行水质模拟,并采用相应水期的实测资料进行验证。结果表明,耦合水质模型计算效果良好,较好地再现了整个研究区域内各水质要素的时空分布规律,成功描述了珠江三角洲河网与河口区的动力与生物化学变化过程。