UNITIVE ANALYSIS SCHEMES ON PROBLEMS OF MULTIPLE MOVING BOUNDARIES WITH 3-D LIQUID-SOLID MULTIPLE NONLINEAR COUPLING FOR UPLIFT OF ANCHOREDLIQUID STORAGE TANKS

In the paper, 3-D analysis method with unitive schemes is set up, which is used to resolve the uplift with multiple moving boundaries and multiple nonlinear coupling for anchored liquid storage tanks. hi it, an algorithm of quasi-harmonious finite elements for arbitrary quadrilateral of thin plates and shells is built up to analyze the multiple coupling problems of general thin plates and shells structures with three dimensions, the complementary equations for analyzing uplifting moving boundary problems are deduced. The axial symmetry and presumption of beam type mode are not used. In it, an algorithm is put forward for analyzing the Navier-Stokes problems of unsteady, three-dimensional, and viscous liquid with sloshing of moving boundary surfaces in large amplitude under ALE frame by scheme of time-split-steps to which linear potential theory is not applied. The algorithms can be used to analyze the solid-liquid multiple nonlinear coupling problems with 3-D moving boundary with friction in multiple places.

ASecond-Order Finite-Difference Method for Compressible Fluids in Domains with Moving Boundaries

In this paper,we describe how to construct a finite-difference shockcapturing method for the numerical solution of the Euler equation of gas dynamics on arbitrary two-dimensional domainΩ,possibly with moving boundary.The boundaries of the domain are assumed to be changing due to the movement of solid objects/obstacles/walls.Although the motion of the boundary could be coupled with the fluid,all of the numerical tests are performed assuming that such a motion is prescribed and independent of the fluid flow.The method is based on discretizing the equation on a regular Cartesian grid in a rectangular domainΩ_(R)⊃Ω.Ωe identify inner and ghost points.The inner points are the grid points located insideΩ,while the ghost points are the grid points that are outsideΩbut have at least one neighbor insideΩ.The evolution equations for inner points data are obtained from the discretization of the governing equation,while the data at the ghost points are obtained by a suitable extrapolation of the primitive variables(density,velocities and pressure).Particular care is devoted to a proper description of the boundary conditions for both fixed and time dependent domains.Several numerical experiments are conducted to illustrate the validity of themethod.Ωe demonstrate that the second order of accuracy is numerically assessed on genuinely two-dimensional problems.

Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.

Stabilization for a Fourth Order Nonlinear Schrödinger Equation in Domains with Moving Boundaries

In the present paper we study the well-posedness using the Galerkin method and the stabilization considering multiplier techniques for a fourth-order nonlinear Schrödinger equation in domains with moving boundaries.We consider two situations for the stabilization:the conservative case and the dissipative case.

An MPI parallel DEM-IMB-LBM framework for simulating fluid-solid interaction problems

The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.

Approach of moving boundary and its application in 3D tidal current simulation of Changjiang (Yangtze) River Estuary based on ECOMSED model

Considering three-dimensional model ECOMSED can not simulate wetting-drying of shoal with its fixed boundary, an approach to represent moving boundary in the model is introduced here. This approach smoothly joints the internal and external mode by making use of wetting and drying technique and is verified by a numerical test which presents a good agreement with the previous test results obtained by other researchers. A three dimensional numerical model is established to simulate the hydrodynamics in spring tide in the Changjiang （Yangtze） River estuary by this modified ECOMSED model which is also validated through the observed field data, the simulation presents a good periodic tidal change. It also successfully simulates the tidal current of computational areas and reproduces the tidal flat intermittent appearance.

Impacts of proppant distribution on development of tight oil reservoirs with threshold pressure gradient

Field evidence indicates that proppant distribution and threshold pressure gradient have great impacts on well productivity.Aiming at the development of unconventional oil reservoirs in Triassic Chang-7 Unit,Ordos Basin of China,we presented an integrated workflow to investigate how(1)proppant placement in induced fracture and(2)non-linear flow in reservoir matrix would affect well productivity and fluid flow in the reservoir.Compared with our research before(Yue et al.,2020),here we extended this study into the development of multi-stage fractured horizontal wells(MFHWs)with large-scale complicated fracture geometry.The integrated workflow is based on the finite element method and consists of simulation models for proppant-laden fluid flow,fracture flow,and non-linear seepage flow,respectively.Simulation results indicate that the distribution of proppant inside the induced cracks significantly affects the productivity of the MFHW.When we assign an idealized proppant distribution instead of the real distribution,there will be an overestimation of 44.98%in daily oil rate and 30.63%in cumulative oil production after continuous development of 1000 days.Besides,threshold pressure gradient(TPG)also significantly affects the well performance in tight oil reservoirs.If we simply apply linear Darcy’s law to the reservoir matrix,the overall cumulative oil production can be overrated by 77%after 1000 days of development.In general,this research provides new insights into the development of tight oil reservoirs with TPG and meanwhile reveals the significance of proppant distribution and non-linear fluid flow in the production scenario design.

Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy's law

A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equa- tions with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be sta- bly numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy＇s flow problem, the exact analytical solution for the case of Darcy＇s flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclu- sion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equa- tions; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pres- sure gradient increasing for the one-dimensional problem.

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.

A New Wetting and Drying Method for Moving Boundary in Shallow Water Flow Models

To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying （WD） method is proposed. In the new method, a ＂predicted water depth＂ is evaluated explicitly based on the simplified shallow water equations and used to determine the status （wet or dry） together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in beth typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.

Characteristic finite difference method and application for moving boundary value problem of coupled system

The 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.A kind of characteristic finite difference schemes is put forward,from which optimal order estimates in l~2 norm are derived for the error in the approximate solutions.The research is important both theoretically and practically for the model analysis in the field,the model numerical method and software development.

A high order boundary scheme to simulate complex moving rigid body under impingement of shock wave

In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary under the impingement of an inviscid shock wave. Based on the high order inverse Lax-Wendroff(ILW) procedure developed in the previous work(TAN, S. and SHU, C. W. A high order moving boundary treatment for compressible inviscid flows. Journal of Computational Physics, 230(15),6023–6036(2011)), in which the authors only considered the translation of the rigid body,we consider both translation and rotation of the body in this paper. In particular, we reformulate the material derivative on the moving boundary with no-penetration condition, and the newly obtained formula plays a key role in the proposed algorithm. Several numerical examples, including cylinder, elliptic cylinder, and NACA0012 airfoil, are given to indicate the effectiveness and robustness of the present method.

Linear Complementary Method for Moving Boundary Problems with Phase Transformation

In this paper, the linear complementary method for moving boundary problems with phase transformation is presented, in which a pair of unknown vectors of heat source with phase transforming and the temperature field can be solved exactly, and a large amount of iterative calculations can be avoided.

A Novel and Optimized Method for Electro-focusing and Moving Neutralization Reaction Boundary Formed by HCl and NaOH

A novel method is developed for electro-focusing and moving neutralization reaction boundary (MNRB) created with HCl and NaOH. The optimized conditions are screened out. By using this method, the experiments are performed on MNRB formed with HCl and NaOH in agarose gel. The experiments are quantitatively in coincidence with the predictions with the theory of moving chemical reaction boundary (MCRB).

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.

Modified characteristic finite difference fractional step method for moving boundary value problem of nonlinear percolation system

A fractional step scheme with modified characteristic finite differences run- ning in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium with moving boundary values. With the help of theoretical techniques including the change of regions, piecewise threefold quadratic interpolation, calculus of variations, multiplicative commutation rule of differ- ence operators, multiplicative commutation rule of difference operators, decomposition of high order difference operators, induction hypothesis, and prior estimates, an optimal order in 12 norm is displayed to complete the convergence analysis of the numerical algo- rithm. Some numerical results arising in the actual simulation of migration-accumulation of oil resources by this method are listed in the last section.

MOVING BOUNDARY PROBLEM FOR DIFFUSION RELEASE OF DRUG FROM A CYLINDER POLYMERIC MATRIX

An approximate analytical solution of moving boundary problem for diffusion release of drug from a cylinder polymeric matrix was obtained by use of refined integral method. The release kinetics has been analyzed for non-erodible matrices with perfect sink condition. The formulas of the moving boundary and the fractional drug release were given. The moving boundary and the fractional drug release have been calculated at various drug loading levels, mid the calculated results were in good agreement with those of experiments. The comparison of the moving boundary in spherical, cylinder, planar matrices has been completed. An approximate formula for estimating the available release time was presented. These results are useful for the clinic experiments. This investigation provides a new theoretical tool for studying the diffusion release of drug from a cylinder polymeric matrix and designing the controlled released drug.

Upwind finite difference method for miscible oil and water displacement problem with moving boundary values

The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational evaluation in 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. For the twodimensional bounded region, the upwind finite difference schemes are proposed. Some techniques, such as the calculus of variations, the change of variables, and the theory of a priori estimates, are used. The optimal orderl2-norm estimates are derived for the errors in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, the model numerical method, and the software development.

Adaptive thermo-fluid moving boundary computations for interfacial dynamics

In this study, we present adaptive moving boundary computation technique with parallel implementation on a distributed memory multi-processor system for large scale thermo-fluid and interfacial flow computations. The solver utilizes Eulerian-Lagrangian method to track moving （Lagrangian） interfaces explicitly on the stationary （Eulerian） Cartesian grid where the flow fields are computed. We address the domain decomposition strategies of Eulerian- Lagrangian method by illustrating its intricate complexity of the computation involved on two different spaces inter- actively and consequently, and then propose a trade-off ap- proach aiming for parallel scalability. Spatial domain decomposition is adopted for both Eulerian and Lagrangian do- main due to easy load balancing and data locality for mini- mum communication between processors. In addition, parallel cell-based unstructured adaptive mesh refinement （AMR） technique is implemented for the flexible local refinement and even-distributed computational workload among processors. Selected cases are presented to highlight the computa- tional capabilities, including Faraday type interfacial waves with capillary and gravitational forcing, flows around varied geometric configurations and induced by boundary conditions and/or body forces, and thermo-fluid dynamics with phase change. With the aid of the present techniques, large scale challenging moving boundary problems can be effectively addressed.

Hybrid Cartesian Grid Method for Moving Boundary Problems

A hybrid Cartesian structured grid method is proposed for solving moving boundary unsteady problems.The near body region is discretized by using the body-fitted structured grids,while the remaining computational domain is tessellated with the generated Cartesian grids.As the body moves,the structured grids move with the body and the outer boundaries of inside grids are used to generate new holes in the outside adaptive Cartesian grid to facilitate data communication.By using the alternating digital tree(ADT)algorithm,the computational time of hole-cutting and identification of donor cells can be reduced significantly.A compressible solver for unsteady flow problems is developed.A cell-centered,second-order accurate finite volume method is employed in spatial discretization and an implicit dual-time stepping low-upper symmetric Gauss-Seidel(LU-SGS)approach is employed in temporal discretization.Geometrybased adaptation is used during unsteady simulation time steps when boundary moves and the flow solution is interpolated from the old Cartesian grids to the new one with inverse distance weighting interpolation formula.Both laminar and turbulent unsteady cases are tested to demonstrate the accuracy and efficiency of the proposed method.Then,a 2-D store separation problem is simulated.The result shows that the hybrid Cartesian grid method can handle the unsteady flow problems involving large-scale moving boundaries.