An Adaptive SDG Method for the Stokes System

Staggered grid techniques are attractive ideas for flow problems due to their more enhanced conservation properties. Recently, a staggered discontinuous Galerkin method is developed for the Stokes system. This method has several distinctive advantages, namely high order optimal convergence as well as local and global conservation properties. In addition, a local postprocessing technique is developed, and the postprocessed velocity is superconvergent and pointwisely divergence-free. Thus, the staggered discontinuous Galerkin method provides a convincing alternative to existing schemes. For problems with corner singularities and flows in porous media, adaptive mesh refinement is crucial in order to reduce the computational cost. In this paper, we will derive a computable error indicator for the staggered discontinuous Galerkin method and prove that this indicator is both efficient and reliable. Moreover, we will present some numerical results with corner singularities and flows in porous media to show that the proposed error indicator gives a good performance.     

On efficient simulation of non-Newtonian flow in saturated porous media with a multigrid adaptive refinement solver

Flow of non-Newtonian fluid in saturated porous media can be described by the continuity equation and the generalized Darcy law. Here we discuss the efficient solution of the resulting second order nonlinear elliptic equation. The equation is discretized by the finite volume method on a cell-centered grid. Local adaptive refinement of the grid is introduced in order to reduce the number of unknowns. We develop a special implementation, that allows us to perform unstructured local refinement in conjunction with the finite volume discretization. Two residual based error indicators are exploited in the adaptive refinement criterion. Second order accurate discretization of the fluxes on the interfaces between refined and non-refined subdomains, as well as on the boundaries with Dirichlet boundary condition, are presented here as an essential part of an accurate and efficient algorithm. A nonlinear full approximation storage multigrid algorithm is developed especially for the above described composite (coarse plus locally refined) grid approach. In particular, second order approximation of the fluxes around interfaces is a result of a quadratic approximation of slave nodes in the multigrid-adaptive refinement (MG-AR) algorithm. Results from numerical solution of various academic and practice-induced problems are presented and the performance of the solver is discussed.     

A fully 3D adaptive finite element simulation and prototype for gas flow in a porous media

This paper presents a three-dimensional (3D) adaptive finite element solution for gas flow in a porous media. The solution obtained is truly 3D and employs a dynamic h-adaptive refinement technique for efficiency. The adaptive procedure uses a new node-based storage mechanism [Wang GH, Tyler JM, Weltman JS, Callahan JD. Advances in Engineering Software including Computing Systems in Engineering 1999;30:31–41]. This node-based structure substantially reduces the memory necessary to store the finite element mesh. A prototype simulator, written in C++, has been implemented for Eugene Island block 305, a multi-well condensate reservoir off the coast of Louisiana to demonstrate the use of this adaptive procedure. Results, presented in this paper, show dramatic agreement with the actual production data. This prototype simulator uses Windows workstations to support a fully dynamic 3D mesh plus mesh generation.     

基于自适应深度迁移学习的多孔介质重构

目前用于多孔介质重构的多点统计法（MPS）等传统方法需要多次扫描训练图像,然后进行后续复杂的概率计算得到模拟结果,导致重构效率较低,模拟过程复杂,因此提出一种基于自适应深度迁移学习的重构方法。首先利用深度神经网络从多孔介质的训练图像中提取复杂特征,然后在深度迁移学习中添加自适应层以减少训练数据和预测数据之间的数据分布差异,最后使用自适应迁移学习复制这些特征来获得与真实训练数据结构相似的重构结果。通过与典型的多孔介质重构方法MPS的比较实验,结果显示在多点连通曲线、变差函数曲线和孔隙度方面,该方法重构质量更好,平均重构耗时从840 s减少到166 s,平均CPU占用率从98%下降到20%,平均内存占用下降了69%。所提方法在保证重构结果质量更好的前提下,显著提高了多孔介质重构的效率。     

Adaptive finite volume methods for displacement problems in porous media

In this paper we consider adaptive numerical simulation of miscible and immiscible displacement problems in porous media, which are modeled by single and two phase flow equations. Using the IMPES formulation of the two phase flow equation both problems can be treated in the same numerical framework. We discretise the equations by an operator splitting technique where the flow equation is approximated by Raviart-Thomas mixed finite elements and the saturation or concentration equation by vertex centered finite volume methods. Using a posteriori error estimates for both approximation schemes we deduce an adaptive solution algorithm for the system of equations and show the applicability in several examples.     

Numerical methods for diffusion-reaction-transport processes in unsaturated porous media

This paper introduces a new grid refinement and coarsening technique for the approximation of partial differential equations including a first order time derivative. This hierarchical movement algorithm is based on the nested iteration method. The combination of this algorithm, a quasi Newton method and the Schur-complement multi-grid method leads to an efficient method for the solution of partial differential equations describing complex real life problems. As a test case, diffusion-reaction-transport processes in heterogeneous unsaturated porous media are considered. Some simulation results are presented. Received: 10 March 1997 / Accepted: 16 June 1997     

多介质拉氏自适应流体动力学软件LAD2D研制及其应用

拉氏方法是内爆动力学过程数值模拟的主要方法.针对高温、高压、多介质和大变形等内爆问题,采用非结构任意多边形网格底层管理、计算过程中网格邻域可变技术,以及拉氏自适应网格加密方法和层次化、模块化程序设计思想,自主研发非结构拉氏自适应网格流体动力学软件LAD2D.从物理模型、计算方法、程序设计、程序验证与确认、大变形问题数值模拟等方面系统地介绍LAD2D.LAD2D对多介质爆轰弹塑性流体大变形问题有很强的适应能力.     

Degenerate Two-Phase Incompressible Flow IV: Local Refinement and Domain Decomposition

This is the fourth paper of a series in which we analyze mathematical properties and develop numerical methods for a degenerate elliptic-parabolic partial differential system which describes the flow of two incompressible, immiscible fluids in porous media. In this paper we describe a finite element approximation for this system on locally refined grids. This adaptive approximation is based on a mixed finite element method for the elliptic pressure equation and a Galerkin finite element method for the degenerate parabolic saturation equation. Both discrete stability and sharp a priori error estimates are established for this approximation. Iterative techniques of domain decomposition type for solving it are discussed, and numerical results are presented.     

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart–Thomas(–Nédélec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters.     

A posteriori error estimation and adaptive mesh refinement for combined thermal-stress finite element analysis

Local and global error estimators and an associated h-based adaptive mesh refinement schemes are proposed for coupled thermal-stress problems. The error estimators are based on the “flux smoothing” technique of Zienkiewicz and Zhu with important modifications to improve convergence performance and computational efficiency. Adaptive mesh refinement is based on the concept of adaptive accuracy criteria, previously presented by the authors for stress-based problems and extended here for coupled thermal-stress problems. Three methods of mesh refinement are presented and numerical results indicate that the proposed method is the most efficient in terms of number of adaptive mesh refinements required for convergence in both the thermal and stress solutions. Also, the proposed method required a smaller number of active degrees of freedom to obtain an accurate solution.     

An adaptive edge-based smoothed point interpolation method for mechanics problems

This paper develops a smoothing domain-based energy (SDE) error indicator and an efficient adaptive procedure using edge-based smoothed point interpolation methods (ES-PIM), in which the strain field is constructed via the generalized smoothing operation over smoothing domains associated with edges of three-node triangular background cells. Because the ES-PIM can produce a close-to-exact stiffness and achieve ‘super-convergence’ and ‘ultra-accurate’ solutions, it is an ideal candidate for adaptive analysis. A SDE error indicator is first devised to make use of the features of the ES-PIM. A local refinement technique based on the Delaunay algorithm is then implemented to achieve high efficiency. The refinement of nodal neighbourhood is accomplished simply by adjusting a scaling factor assigned to control local nodal density. Intensive numerical studies, including the problems with stress concentration and solution singularity, demonstrate that the proposed adaptive procedure is effective and efficient in producing solutions with desired accuracy.     

Observer‐based fault‐tolerant control of hypersonic scramjet vehicles in the presence of actuator faults and saturation

An adaptive sliding mode observer (SMO)–based fault‐tolerant control method taking into consideration of actuator saturation is proposed for a hypersonic scramjet vehicle (HSV) under a class of time‐varying actuator faults. The SMO is designed to robustly estimate the HSV states and reconstruct the fault signals. The adaptive technique is integrated into the SMO to approximate the unknown bounds of system uncertainties, actuator faults, and estimation errors. The robust SMO synthesis condition, which can be formulated as a set of linear matrix inequalities, is improved by relaxing structure constraints to the Lyapunov matrix. An anti‐windup feedback control law, which utilizes the estimated HSV states and the fault signals, is designed to counteract the negative effects of actuator saturation induced by actuator faults. Simulation results demonstrate that the proposed approach can guarantee stability and maintain performance of the closed‐loop system in the presence of HSV actuator faults and saturation.     

A new treatment of capillarity to improve the stability of IMPES two-phase flow formulation

In this paper, we present an efficient numerical method for two-phase immiscible flow in porous media with different capillarity pressures. In highly heterogeneous permeable media, the saturation is discontinuous due to different capillary pressure functions. One popular scheme is to split the system into a pressure and a saturation equation, and to apply IMplicit Pressure Explicit Saturation (IMPES) approach for time stepping. One disadvantage of IMPES is instability resulting from the explicit treatment for capillary pressure. To improve stability, the capillary pressure is usually incorporated in the saturation equation which gradients of saturation appear. This approach, however, does not apply to the case of different capillary pressure functions for multiple rock-types, because of the discontinuity of saturation across rock interfaces. In this paper, we present a new treatment of capillary pressure, which appears implicitly in the pressure equation. Using an approximation of capillary function, we substitute the implicit saturation equation into the pressure equation. The coupled pressure equation will be solved implicitly and followed by the explicit saturation equation. Five numerical examples are provided to demonstrate the advantages of our approach. Comparison shows that our proposed method is more efficient and stable than the classical IMPES approach.     

Controlled cost of adaptive mesh refinement in practical 3D finite element analysis

In this paper, attention is restricted to mesh adaptivity. Traditionally, the most common mesh adaptive strategies for linear problems are used to reach a prescribed accuracy. This goal is best met with an h-adaptive scheme in combination with an error estimator. In an industrial context, the aim of the mechanical simulations in engineering design is not only to obtain greatest quality but more often a compromise between the desired quality and the computation cost (CPU time, storage, software, competence, human cost, computer used). In this paper we propose the use of alternative mesh refinement with an h-adaptive procedure for 3D elastic problems. The alternative mesh refinement criteria allow to obtain the maximum of accuracy for a prescribed cost. These adaptive strategies are based on a technique of error in constitutive relation (the process could be used with other error estimators) and an efficient adaptive technique which automatically takes into account the steep gradient areas. This work proposes a 3D method of adaptivity with the latest version of the INRIA automatic mesh generator GAMHIC3D.     

Robust Tracking Control of Uncertain MIMO Nonlinear Systems with Application to UAVs

In this paper, we consider the robust adaptive tracking control of uncertain multi-input and multi-output (MIMO) nonlinear systems with input saturation and unknown external disturbance. The nonlinear disturbance observer (NDO) is employed to tackle the system uncertainty as well as the external disturbance. To handle the input saturation, an auxiliary system is constructed as a saturation compensator. By using the backstepping technique and the dynamic surface method, a robust adaptive tracking control scheme is developed. The closed-loop system is proved to be uniformly ultimately bounded thorough Lyapunov stability analysis. Simulation results with application to an unmanned aerial vehicle (UAV) demonstrate the effectiveness of the proposed robust control scheme.      

Direct adaptive fuzzy backstepping control of uncertain nonlinear systems in the presence of input saturation

In this paper, a novel direct adaptive fuzzy control approach is presented for uncertain nonlinear systems in the presence of input saturation. Fuzzy logic systems are directly used to tackle unknown nonlinear functions, and the adaptive fuzzy tracking controller is constructed by using the backstepping recursive design techniques. To overcome the problem of input saturation, a new auxiliary design system and Nussbaum gain functions are incorporated into the control scheme, respectively. It is proved that the proposed control approach can guarantee that all the signals of the resulting closed-loop system are semi-globally uniformly ultimately bounded (SGUUB), and the tracking error converges to a small neighborhood of the origin. A simulation example is included to illustrate the effectiveness of the proposed approach. Two key advantages of the scheme are that (i) the direct adaptive fuzzy control method is proposed for uncertain nonlinear system with input saturation by using Nussbaum function technique and (ii) The number of the online adaptive learning parameters is reduced.     

Adaptive Fuzzy Backstepping Tracking Control for Flexible Robotic Manipulator

In this paper, an adaptive fuzzy state feedback control method is proposed for the single-link robotic manipulator system. The considered system contains unknown nonlinear function and actuator saturation. Fuzzy logic systems (FLSs) and a smooth function are used to approximate the unknown nonlinearities and the actuator saturation, respectively. By combining the command-filter technique with the backstepping design algorithm, a novel adaptive fuzzy tracking backstepping control method is developed. It is proved that the adaptive fuzzy control scheme can guarantee that all the variables in the closed-loop system are bounded, and the system output can track the given reference signal as close as possible. Simulation results are provided to illustrate the effectiveness of the proposed approach.      

注蒸汽热采稠油油藏数值模拟自适应网法计算软件的设计开发

注蒸汽热采稠油是一个带有相变的多孔介质中多相流动问题。利用自适应网格法动态追踪油藏中的温度和饱和度锋面,以大幅度提高计算速度。在对渗透率空间分布不均匀、油藏区域分布不同种类输运性质岩石、裂缝-孔隙双重介质、非规则边界及内部含有断层油藏等复杂地层注蒸汽热采稠油问题开展自适应网格法的研究基础上,自主开发了注蒸汽热采稠油油藏数值模拟自适应网格法计算软件。该软件具有较高的计算精度和稳定性,可满足科研生产需求。     

Automatic nested refinement: a technique for the generation of high quality multi-block structured grids for multi-scale problems using GridPro

A technique for the generation of conformal adaptive refinement of hex meshes is presented. Automatic nested refinement is a technique for generating recursively nested topology automatically. It can be applied inside GridPro’s topological paradigm to generate block structured grids, which can resolve tiny features in the problem while providing for a smooth and consistent way to transition to a larger scale. The selection of self-similar, ‘fractal-like’ topological templates makes the selection of number of levels easy, while making the technique feasible for infinite levels of adaptive refinement. The technique has been programmed and integrated into the GUI of GridPro, making it very accessible and easy to use. This method illustrates a way of generating structured grids in an unstructured way, made possible because of the topological paradigm of GridPro. Because the technique uses the topological paradigm, it inherits all the advantages the paradigm offers, including the ability to handle very complex geometries, parametric variation of surfaces, and the ease of use, speed and quality of GridPro. The technique has been illustrated using a variety of applications. This method has proved to be a fast, efficient, automatic and reliable means to perform physical simulations that have a disparity in scale.     

Nonlinear Robust Adaptive Deterministic Control for Flexible Hypersonic Vehicles in the Presence of Input Constraint

A nonlinear deterministic robust control scheme is developed for a flexible hypersonic vehicle with input saturation. Firstly, the model analysis is conducted for the hypersonic vehicle model via the input‐output linearized technique. Secondly, the sliding mode manifold is designed based on homogeneity theory. Then an adaptive high order sliding mode control scheme is proposed to achieve tracking for the step change in altitude and velocity for hypersonic vehicles where the uncertainty boundary is unknown. Furthermore, the control input constraint is investigated and another new adaptive law is proposed to estimate the uncertainties and to guarantee the stability of the system with input saturation. Finally, the simulation results are provided to demonstrate the effectiveness of the proposed method.