THE CHARACTERISTIC FINITE ELEMENT ALTERNATING-DIRECTION METHOD AND ANALYSIS FOR THREE-DIMENSIONAL NUMERICAL RESERVOIR SIMULATION

Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model can be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the three-dimensional characteristic of large-scale science-engineering computation, we put forward a kind of characteristic finite element alternating-direction schemes and obtain optimal order estimates in L~2 norm for the error in the approximate assumption.     

A Spatial Geometric Nonlinearity Spline Beam Element With Nodal Parameters Containing Strains北大核心CSCD

Many slender rods in engineering can be modeled as Euler-Bernoulli beams. For the analysis of their dynamic behaviors, it is necessary to establish the dynamic models for the flexible multi-body systems. Geometric nonlinear elements with absolute nodal coordinates help solve a large number of dynamic problems of flexible beams, but they still face such problems as shear locking, nodal stress discontinuity and low computation efficiency. Based on the theory of large deformation beams’ virtual power equations, the functional formulas between displacements and rotation angles at the nodes were established, which can satisfy the deformation coupling relationships. The generalized strains to describe geometric nonlinear effects in this case were derived. Some parameters of boundary nodes were replaced by axial strains and sectional curvatures to obtain a more accurate and concise constraint method for applying external forces. To improve the numerical efficiency and stability of the system’s motion equations, a model-smoothing method was used to filter high frequencies out of the model. The numerical examples verify the rationality and effectiveness of the proposed element. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

The Theory of Finite Models without Equal Sign

In this paper, it is the first time ever to suggest that we study the model theory of all finite structures and to put the equal sign in the same situtation as the other relations. Using formulas of infinite lengths we obtain new theorems for the preservation of model extensions, submodels, model homomorphisms and inverse homomorphisms. These kinds of theorems were discussed in Chang and Keisler＇s Model Theory, systematically for general models, but Gurevich obtained some different theorems in this direction for finite models. In our paper the old theorems manage to survive in the finite model theory. There are some differences between into homomorphisms and onto homomorphisms in preservation theorems too. We also study reduced models and minimum models. The characterization sentence of a model is given, which derives a general result for any theory T to be equivalent to a set of existential-universal sentences. Some results about completeness and model completeness are also given.     

THE CHARACTERISTIC FINITE DIFFERENCE METHOD FOR THREE-DIMENSIONAL MOVING BOUNDARY VALUE PROBLEM

The software for oil-gas transport and accumulation is to describe the history of oil-gas transport and accumulation in basin evolution. It is of great value in rational evaluation of prospecting and exploiting oil-gas resources. The mathematical model can be discribed as a coupled system of nonlinear partial differential equations with moving boundary value problem. For a generic case of the three-demensional bounded region, bi this thesis, the effects of gravitation、buoyancy and capillary pressure are considered, we put forward a kind of characteristic finite difference schemes and make use thick and thin grids to form a complete set, and of calculus of vaviations, the change of variable, the theory of prior estimates and techniques, Optimal order estimates in l~2 norm are derived for the error in approximate assumption, Thus we have completely solved the well-known theoretical problem proposed by J. Douglas, Jr.     

CONSTRAINED RATIONAL CUBIC SPLINE AND ITS APPLICATION

1. IntroductionDesign of high quality, manufacturable surfaces, such as the outer shape of a ship, car oraeroplane, is an important yet challenging task in today's manufacturing industries. Althoughsignificam progress has been made in the last decade in developing and commercializing pro--duction quality CAD tools, demand for more effective tools is still high due to the ever increajsein model complexity and the needs to address and incorporate manufacturing requirements inthe early stage of …     

Finite Integral Transform Solutions for Free Vibrations of Rectangular Thin Plates With Mixed Boundary Constraints北大核心CSCD

Analytical solutions, with unique research value, can serve as benchmarks for empirical formulas and numerical methods, a tool for rapid parameter analysis and optimization, and a theoretical basis for experimental designs. Conventional analytical methods, e.g., the Lévy solution method, are only applicable to mechanical problems of plates and shells with opposite simply-supported edges, which, however, may fail to obtain analytical solutions for the issues with complex boundary constraints. In recent years, the finite integral transform method for plate and shell problems was developed to deal with non-Lévy-type plates and shells, but it is still infeasible to solve the mixed boundary constrains-induced complex boundary value problems of higher-order partial differential equations. Herein, for the first time, the finite integral transform method was combined with the sub-domain decomposition technique to solve the free vibrations of rectangular thin plates with mixed boundary constraints. The rectangular plate was first divided into 2 sub-domains according to the mixed boundary constraints, and the 2 sub-domains were solved analytically with the finite integral transform method. Finally, the continuity conditions were introduced to obtain the analytical solution of the original problem. Based on the side spot-welded cantilever plates commonly used in engineering, the free vibration problem of a rectangular thin plate with 1 edge subjected to clamped-simply supported constraints and the other 3 edges free, was analyzed. The obtained natural frequencies and mode shapes are in good agreement with those from the finite element method as well as the solutions in literature, thus verifying the accuracy of the proposed method. The solution procedure of the finite integral transform method can be implemented based on the governing equations without any assumption of the solution form. Therefore, this strict analytical method is widely applicable to complex boundary value problems of higher-order partial differential equations for such mechanical problems of plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

UPWIND FINITE VOLUME SCHEMES FOR SEMICONDUCTOR DEVICE

The mathematical model of semiconductor devices is described by the initial boundary value problem of a system of three nonlinear partial differential equations. One equation in elliptic form is for the electrostatic potential; two equations of convection-dominated diffusion type are for the electron and hole concentrations. Finite volume element procedure are put forward for the electrostatic potential, while upwind     

Numerical Investigation of Particle Focusing Patterns in Laminar Pipe Flow With High Reynolds Numbers北大核心CSCD

The inertial focusing characteristics of particles in laminar flow pipes with high Re numbers were studied based on the “relative motion model”. In order to solve the problem of long pipes with high Re number flow, periodic boundary conditions were imposed on the inlet and outlet of the pipe. The research results show that the use of periodic boundary conditions can effectively reduce the computational, and the mechanical properties of particles in high Re flow can be calculated by using L=4D pipe. The difference from the low Re number is that as the Re number continues to increase,the lift force of the particles in the radial direction is no longer distributed as a parabola. The lift curve has a concave area between r+ =0.5 ~ 0.7, and there is a tendency for a new inertial focus point to appear in this section. By means of particles of a+ =1/17 for Re > 1 000, this new focus point position is solvable. In addition, in the analysis of the flow field, a secondary flow occurs around the particle, and its intensity gradually increases with the Re number and the closeness of the particle to the wall. The generation of the secondary flow affects the spatial distribution of the particle lift. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

A Variational Finite Element Model for Large-Eddy Simulations of Turbulent Flows

The authors introduce a new Large Eddy Simulation model in a channel,based on the projection on finite element spaces as filtering operation in its variational form,for a given triangulation{Th}h>0.The eddy viscosity is expressed in terms of the friction velocity in the boundary layer due to the wall,and is of a standard sub grid-model form outside the boundary layer.The mixing length scale is locally equal to the grid size.The computational domain is the channel without the linear sub-layer of the boundary layer.The no-slip boundary condition(or BC for short)is replaced by a Navier(BC)at the computational wall.Considering the steady state case,the authors show that the variational finite element model they have introduced,has a solution(vh,ph)h>0that converges to a solution of the steady state Navier-Stokes equation with Navier BC.     

THE UPWIND FINITE DIFFERENCE METHOD FOR MOVING BOUNDARY VALUE PROBLEM OF COUPLED SYSTEM

Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.     

Evaluation of the approximated diffusion flamelet concept using fuels with different chemical complexity

The ability of flamelet models to reproduce turbulent combustion in devices such as diesel engines or gas turbines has enhanced the usage of these approaches in Computational Fluid Dynamics (CFD) simulations. The models based on turbulent look-up tables generated from counterflow laminar diffusion flames (DF model) permit drastic reduction of the computational cost of the CFD calculation. Nevertheless, for complex molecular fuels, such as n-heptane, the oxidation process involves hundreds of species and the calculation of the transport equations together with the ODE system that models the chemical kinetics for the DF solution becomes unaffordable for industrial devices where hundreds of flamelets are required. In this context, new hypotheses have to be introduced in order to reduce the computational cost maintaining the coherence of the combustion process. Recently, a new model known as Approximated Diffusion Flamelet (ADF) has been proposed with the aim of solving the turbulent combustion for complex fuels in a reduced time. However, the validity of this model is still an open question and has to be verified in order to justify subsequent CFD calculations. This work assesses the ADF model and its ability to reproduce accurately the combustion process and its main parameters for three fuels with different chemical complexity and boundary conditions by its comparison with the DF model. Results show that although some discrepancies arise, the ADF model has the ability to correctly describe the ignition delay and the combustion structure in the auto-ignition zone that is the most relevant one for industrial processes.     

Monotone iterative method of upper and lower solutions applied to a multilayer combustion model in porous media

This paper presents a new nonlinear reaction–diffusion–convection system coupled with a system of ordinary differential equations that models a combustion front in a multilayer porous medium. The model includes heat transfer between the layers and heat loss to the external environment. A few assumptions are made to simplify the model, such as incompressibility; then, the unknowns are determined to be the temperature and fuel concentration in each layer. When the fuel concentration in each layer is a known function, we prove the existence and uniqueness of a classical solution for the initial and boundary value problem for the corresponding system. The proof uses a new approach for combustion problems in porous media. We construct monotone iterations of upper and lower solutions and prove that these iterations converge to a unique solution for the problem, first locally and then, in time, globally.     

Einfluss eines Verbrennungsvorganges auf den Wärme- und Stoffaustausch in einer laminaren Grenzschicht

Summary The influence of combustion on heat and mass transfer is investigated on the following model. A mixture of an inert with a combustible gas (air) flows in steady, laminar flow over a flat plate. A mass flux of gaseous fuel away from the plate surface is produced by some means. Combustion is assumed to occur with very fast reaction rate so that the process is purely controlled by diffusion and the equilibrium is assumed as very close to complete combustion. It is studied under which conditions the combustion occurs at the surface or when the flame is displaced into the boundary layer. The influence of combustion on the heat transfer from a hot gas to the plate surface is calculated, for the condition that combustion occurs at the surface.      

The Rotation of a Gravitating Sphere in a Monatomic Gas,II

The purpose of this work is to study the fluid motion caused by the high speed rotation of a gravitating sphere in a monatomic gas. It has been possible to find a stable steady solution only for very small Prandtl number, which can be interpreted to mean an optically thick gas. The flow is characterized by a flat radial jet in the equatorial plane and a viscous boundary layer on the spherical surface which, in some cases, lies beneath a thermal boundary layer. That the outer region must be hydrostatic puts very stringent constraints on the associated velocity field which necessitate still another boundary layer on the sphere. This last layer is shown to be unstable to small disturbances in certain temperature ranges. Finally, a similar solution that exists for order one Prandtl number must be disregarded because this last boundary layer is always unstable.     

The second‐order upwind finite difference fractional steps method for moving boundary value problem of oil‐water percolation

The research of the three‐dimensional (3D) compressible miscible (oil and water) displacement problem with moving boundary values is of great value to the history of oil‐gas transport and accumulation in basin evolution, as well as to the rational evaluation in prospecting and exploiting oil‐gas resources, and numerical simulation of seawater intrusion. The mathematical model can be described as a 3D‐coupled system of nonlinear partial differential equations with moving boundary values. For a generic case of 3D‐bounded region, a kind of second‐order upwind finite difference fractional steps schemes applicable to parallel arithmetic is put forward. Some techniques, such as the change of variables, calculus of variations, and the theory of a priori estimates, are adopted. Optimal order estimates in l² norm are derived for the errors in approximate solutions. The research is important both theoretically and practically for model analysis in the field, for model numerical method and for software development. Thus, the well‐known problem has been solved.Copyright © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1103–1129, 2014     

Nonclassical Transonic Boundary Layers: Toward Overcoming Dead-End Situations in High-Speed Aerodynamics

Analytical models of unsteady free viscous-inviscid interaction of gas flows at transonic speeds, i.e., a transonic boundary layer with self-induced pressure (nonclassical boundary layer) are considered. It is shown that an adequate flow model can be constructed by applying methods of singular perturbations. The results of a comparative analysis of classical and regularized stability models for a boundary layer with self-induced pressure in the case of interaction at transonic speeds are overviewed.     

The applicability of continuum models in the transitional regime of hypersonic flow over blunt bodies

Hypersonic rarefied gas flow over blunt bodies in the transitional flow regime (from continuum to free-molecule) is investigated. Asymptotically correct boundary conditions on the body surface are derived for the full and thin viscous shock layer models. The effect of taking into account the slip velocity and the temperature jump in the boundary condition along the surface on the extension of the limits of applicability of continuum models to high free-stream Knudsen numbers is investigated. Analytic relations are obtained, by an asymptotic method, for the heat transfer coefficient, the skin friction coefficient and the pressure as functions of the free-stream parameters and the geometry of the body in the flow field at low Reynolds number; the values of these coefficients approach their values in free-molecule flow (for unit accommodation coefficient) as the Reynolds number approaches zero. Numerical solutions of the thin viscous shock layer and full viscous shock layer equations, both with the no-slip boundary conditions and with boundary conditions taking into account the effects slip on the surface are obtained by the implicit finite-difference marching method of high accuracy of approximation. The asymptotic and numerical solutions are compared with the results of calculations by the Direct Simulation Monte Carlo method for flow over bodies of different shape and for the free-stream conditions corresponding to altitudes of 75–150 km of the trajectory of the Space Shuttle, and also with the known solutions for the free-molecule flow regine. The areas of applicability of the thin and full viscous shock layer models for calculating the pressure, skin friction and heat transfer on blunt bodies, in the hypersonic gas flow are estimated for various free-stream Knudsen numbers.     

Lattice Boltzmann simulation for high-speed compressible viscous flows with a boundary layer

In this study, the lattice Boltzmann method is employed for simulating high-speed compressible viscous flows with a boundary layer. The coupled double-distribution-function lattice Boltzmann method proposed by Li et al. (2007) is employed because of its good numerical stability and non-free-parameter feature. The non-uniform mesh construction near the wall boundary in fine grids is combined with an appropriate wall boundary treatment for the finite difference method in order to obtain accurate spatial resolution in the boundary layer problem. Three typical problems in high-speed viscous flows are solved in the lattice Boltzmann simulation, i.e., the compressible boundary layer problem, shock wave problem, and shock boundary layer interaction problem. In addition, in-depth comparisons are made with the non-oscillatory and non-free-parameter dissipation (NND) scheme and second order upwind scheme in the present lattice Boltzmann model. Our simulation results indicate the great potential of the lattice Boltzmann method for simulating high-speed compressible viscous flows with a boundary layer. Further research is needed (e.g., better numerical models and appropriate finite difference schemes) because the lattice Boltzmann method is still immature for high-speed compressible viscous flow applications.     

Existence of nonlinear boundary layer to the Boltzmann equation with reverse reflection boundary condition

Many physical models have boundaries. When the Boltzmann equation is used to study a physical problem with boundary, there usually exists a layer of width of the order of the Knudsen number along the boundary. Hence, the research on the boundary layer problem is important both in mathematics and physics. Based on the previous work, in this paper, we consider the existence of boundary layer solution to the Boltzmann equation for hard sphere model with positive Mach number. The boundary condition is imposed on incoming particles of reverse reflection type, and the solution is assumed to approach to a global Maxwellian in the far field. Similar to the problem with Dirichlet boundary condition studied in [S. Ukai, T. Yang, S.H. Yu, Nonlinear boundary layers of the Boltzmann equation: I. Existence, Comm. Math. Phys. 3 (2003) 373-393], the existence of a solution is shown to depend on the Mach number of the far field Maxwellian. Moreover, there is an implicit solvability condition on the boundary data. According to the solvability condition, the co-dimension of the boundary data related to the number of the positive characteristic speeds is obtained.     

弹性约束边界条件下矩形蜂窝夹芯板的自由振动分析

蜂窝夹芯板在飞行器、高速列车等领域有广泛的用途，对其开展振动分析具有明确的科学价值及工程意义.为区别于诸简支等传统约束边界，提出了弹性约束边界下蜂窝夹芯板结构的自由振动特性分析方法.具体来说，首先通过将蜂窝夹芯层等效为各向异性板，将夹芯板问题转变为三层板结构.进一步地，将板结构的位移场函数由改进的二维Fourier级数表示，并基于能量原理的Rayleigh Ritz法得到结构的固有频率和固有振型，理论预测结果与数值模拟分析吻合较好.提出的理论模型可用于系统讨论约束边界对蜂窝夹芯结构自由振动特性的影响，为此类结构的约束方案设计提供理论依据.