Partially explicit splitting scheme with explicit–implicit-null method for nonlinear multiscale flow problems

Fully implicit schemes with second‐order time evolutions have been applied to simulate nonlinear diffusion problems precisely for a long time, but there is seldom theoretical study for either their convergence properties or efficient iterations. Here, a second‐order time evolution fully implicit scheme for two‐dimensional nonlinear divergence diffusion problem is analyzed. The unique existence of its solution is given. Two new methods are provided to prove its convergence, including entire inductive hypothesis reasoning and a two‐step reasoning process. Rigorous analysis shows the scheme is stable; its solution has second‐order convergence in both space and time to the exact solution of the problem. The convergence is applied to analyze a Newton iteration accelerating the computation and show its quadratic convergent speed and second‐order accuracy. The reasoning techniques also adapt to first‐order time accuracy schemes, and can be extended to analyze a wide class of nonlinear schemes for nonlinear problems. Numerical tests highlight the theoretical results and demonstrate the high performance of the algorithms. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 121–140, 2016

The main objective of this paper is to reduce the dimensionality of unknown coefficient vectors of finite-element (FE) solutions in two-grid (CN) FE (TGCNFE) format for the nonlinear unsaturated soil water flow problem by using a proper orthogonal decomposition (POD) and to design a new reduced-dimension iteration TGCNFE (RDITGCNFE). For this objective, a new time semi-discrete CN (TSDCN) scheme for the nonlinear unsaturated soil water flow problem is first designed and the existence, stability, and error estimates of TSDCN solutions are demonstrated. Subsequently, a new TGCNFE format for the nonlinear unsaturated soil water flow problem is designed and the existence, unconditional stability, and error estimates of TGCNFE solutions are demonstrated. Next, a new RDITGCNFE format with the same FE basis functions as the TGCNFE format is built by the POD method and the existence, unconditional stability, and error estimates of RDITGCNFE solutions are discussed. Ultimately, the rightness of theory results and the superiority of the RDITGCNFE format are verified by two sets of numerical tests. It is worth noting that the RDITGCNFE format differs completely from all previous reduced-dimension methods, including the authors’ previous works. Therefore, the study of this paper can not only provide a new theoretical method for the dimensionality reduction of numerical models for nonlinear problems but also provide an algorithm implementation technology for the numerical simulation of practical engineering problems.

Operator-splitting methods for the 2D convective Cahn–Hilliard equation

This work presents the development of an explicit/implicit particle finite element method (PFEM) for the 2D modeling of metal cutting processes. The purpose is to study the efficiency of implicit and explicit time integration schemes in terms of precision, accuracy and computing time. The formulation for implicit and explicit time marching schemes is developed, and a detailed study on the explicit solution steps is presented. The PFEM remeshing procedures for insertion and removal of particles have been improved to model the multiple scales of time and/or space of the solution. The detection and treatment of the rigid tool contact are presented for both, implicit and explicit schemes. The performance of explicit/implicit integration is studied with a set of different two-dimensional orthogonal cutting tests of AISI 4340 steel at cutting speeds ranging from 1 m/s up to 30 m/s. It was shown that if the correct selection of the time integration scheme is made, the computing time can decrease up to 40 times. It allows us to affirm that the computing time of the PFEM simulations can be excessive due to the used time marching scheme independently of the meshing process. As a practical result, a set of recommendations to select the time integration schemes for a given cutting speed are given. This is intended to minimize one of the negative constraints pointed out by the industry when using metal cutting simulators.

Diffusive relaxation systems provide a general framework to approximate nonlinear diffusion problems, also in the degenerate case (Aregba-Driollet et al. in Math. Comput. 73(245):63---94, 2004; Boscarino et al. in Implicit-explicit Runge-Kutta schemes for hyperbolic systems and kinetic equations in the diffusion limit, 2011; Cavalli et al. in SIAM J. Sci. Comput. 34:A137---A160, 2012; SIAM J. Numer. Anal. 45(5):2098---2119, 2007; Naldi and Pareschi in SIAM J. Numer. Anal. 37:1246---1270, 2000; Naldi et al. in Surveys Math. Indust. 10(4):315---343, 2002). Their discretization is usually obtained by explicit schemes in time coupled with a suitable method in space, which inherits the standard stability parabolic constraint. In this paper we combine the effectiveness of the relaxation systems with the computational efficiency and robustness of the implicit approximations, avoiding the need to resolve nonlinear problems and avoiding stability constraints on time step. In particular we consider an implicit scheme for the whole relaxation system except for the nonlinear source term, which is treated though a suitable linearization technique. We give some theoretical stability results in a particular case of linearization and we provide insight on the general case. Several numerical simulations confirm the theoretical results and give evidence of the stability and convergence also in the case of nonlinear degenerate diffusion.

We propose an approach for large-scale non-separable nonlinear multicommodity flow problems by solving a sequence of subproblems which can be addressed by commercial solvers. Using a combination of solution methods such as modified gradient projection, shortest path algorithm and golden section search, the approach can handle general problem instances, including those with (i) non-separable cost, (ii) objective function not available analytically as polynomial but are evaluated using black-boxes, and (iii) additional side constraints not of network flow types. Implemented as a toolbox in commercial solvers, it allows researchers and practitioners, currently conversant with linear instances, to easily manage large-scale convex instances as well. In this article, we compared the proposed algorithm with alternative approaches in the literature, covering both theory and large test cases. New test cases with non-separable convex costs and non-network flow side constraints are also presented and evaluated. The toolbox is available free for academic use upon request.

We address the problem of generating Dynamic Compact Thermal Models for non-linear and non-homogeneous heat diffusion problems. We propose a novel Model Order Reduction technique that is based on i) a specifically designed discretization that properly models the non-linear behaviour of different materials, and ii) the introduction of an improved Hyper Reduction method that dramatically reduces the problem dimension. Our experiments show that our method performs similarly to full Finite Volume simulation in terms of accuracy, while being much more computationally efficient.

Learning-based multi-continuum model for multiscale flow problems

Model reduction for multi-scale transport problems using model-form preserving least-squares projections with variable transformation

A numerical model is presented for the computation of unsteady two‐fluid interfaces in nonlinear porous media flow. The nonlinear Forchheimer equation is included in the Navier‐Stokes equations for porous media flow. The model is based on capturing the interface on a fixed mesh domain. The zero level set of a pseudo‐concentration function, which defines the interface between the two fluids, is governed by a time‐dependent advection equation. The time‐dependent Navier‐Stokes equations and the advection equation are spatially discretized by the finite element (FE) method. The fully coupled implicit time integration scheme and the explicit forward Eulerian scheme are implemented for the advancement in time. The trapezoidal rule is applied to the fully implicit scheme, while the operator‐splitting algorithm is used for the velocity‐pressure segregation in the explicit scheme. The spatial and time discretizations are stabilized using FE stabilization techniques. Numerical examples of unsteady flow of two‐fluid interfaces in an earth dam are investigated.

A finite volume method based on radial basis functions for two-dimensional nonlinear diffusion equations

A qualitative analysis and numerical simulations of a nonlinear second-order anisotropic diffusion problem with non-homogeneous Cauchy–Neumann boundary conditions

Estimation of concentration-dependent diffusion coefficient in pressure-decay experiment of heavy oils and bitumen

Many important initial value problems have the property that energy is nonincreasing in time. Energy stable methods, also referred to as strongly stable methods, guarantee the same property discretely. We investigate requirements for conditional energy stability of explicit Runge--Kutta methods for nonlinear or nonautonomous problems. We provide both necessary and sufficient conditions for energy stability over these classes of problems. Examples of conditionally energy stable schemes are constructed, and an example is given in which unconditional energy stability is obtained with an explicit scheme.

In this paper high order time integration schemes are investigated for the integration of uid-structure interaction (FSI) model problem with deforming domains. The integration algorithm consists of a partitioned scheme in which a combination of implicit and explicit (IMEX) high order Runge-Kutta time integration schemes is used: the uid and structure systems are decoupled at their interface and both systems are integrated using the implicit system under a given boundary condition at the interface. The inuence of this boundary condition (the coupling term) is integrated using the explicit Runge-Kutta scheme. The uid dynamics are written in the arbitrary Lagrange Eulerian formulation to cope with the deforming uid mesh. An expression is proposed for the mesh face velocities of the deforming mesh such that the discrete geometric conservation law is satised. For the model problem numerical investigations are made into the eciency of these high order time integration schemes compared to lower order methods. The IMEX scheme retains the order of the implicit and explicit schemes and obtains the same accuracy as a second order staggered scheme for less computational work.