NUMERICAL GRADIENT FLOW DISCRETIZATION OF VISCOUS THIN FILMS ON CURVED GEOMETRIES

2013 ◽  
Vol 23 (05) ◽  
pp. 917-947 ◽  
Author(s):  
MARTIN RUMPF ◽  
ORESTIS VANTZOS
Keyword(s):  
Thin Films ◽  
Transport Equation ◽  
Degrees Of Freedom ◽  
Gradient Flow ◽  
Time Discretization ◽  
Time Step ◽  
Time Stepping ◽  
Exterior Calculus ◽  
Natural Time ◽  
Space Discretization

The evolution of a viscous thin film on a curved geometry is numerically approximated based on the natural time discretization of the underlying gradient flow. This discretization leads to a variational problem to be solved at each time step, which reflects the balance between the decay of the free (gravitational and surface) energy and the viscous dissipation. Both dissipation and energy are derived from a lubrication approximation for a small ratio between the characteristic film height and the characteristic length scale of the surface. The dissipation is formulated in terms of a corresponding flux field, whereas the energy primarily depends on the fluid volume per unit surface, which is a conserved quantity. These two degrees of freedom are coupled by the underlying transport equation. Hence, one is naturally led to a PDE-constrained optimization problem, where the variational time stepping problem has to be solved under the constraint described by the transport equation. For the space discretization a discrete exterior calculus approach is investigated. Various applications demonstrate the qualitative and quantitative behavior of one- and two-dimensional thin films on curved geometries.

scholarly journals Contrast-Independent, Partially-Explicit Time Discretizations for Nonlinear Multiscale Problems

Mathematics ◽  
2021 ◽  
Vol 9 (23) ◽  
pp. 3000
Author(s):  
Eric T. Chung ◽  
Yalchin Efendiev ◽  
Wing Tat Leung ◽  
Wenyuan Li
Keyword(s):  
Linear Equations ◽  
Nonlinear Problems ◽  
Time Discretization ◽  
Numerical Results ◽  
Explicit Methods ◽  
Implicit Methods ◽  
Time Step ◽  
Multiscale Problems ◽  
Time Stepping ◽  
Space Decomposition

This work continues a line of work on developing partially explicit methods for multiscale problems. In our previous works, we considered linear multiscale problems where the spatial heterogeneities are at the subgrid level and are not resolved. In these works, we have introduced contrast-independent, partially explicit time discretizations for linear equations. The contrast-independent, partially explicit time discretization divides the spatial space into two components: contrast dependent (fast) and contrast independent (slow) spaces defined via multiscale space decomposition. Following this decomposition, temporal splitting was proposed, which treats fast components implicitly and slow components explicitly. The space decomposition and temporal splitting are chosen such that they guarantees stability, and we formulated a condition for the time stepping. This condition was formulated as a condition on slow spaces. In this paper, we extend this approach to nonlinear problems. We propose a splitting approach and derive a condition that guarantees stability. This condition requires some type of contrast-independent spaces for slow components of the solution. We present numerical results and show that the proposed methods provide results similar to implicit methods with a time step that is independent of the contrast.

scholarly journals Error Estimates of a Continuous Galerkin Time Stepping Method for Subdiffusion Problem

2021 ◽  
Vol 88 (3) ◽  
Author(s):  
Yuyuan Yan ◽  
Bernard A. Egwu ◽  
Zongqi Liang ◽  
Yubin Yan
Keyword(s):  
Piecewise Linear ◽  
Time Discretization ◽  
Convergence Order ◽  
Step Size ◽  
Time Step ◽  
Transform Method ◽  
Nonsmooth Data ◽  
Time Step Size ◽  
Time Stepping ◽  
Continuous Galerkin

AbstractA continuous Galerkin time stepping method is introduced and analyzed for subdiffusion problem in an abstract setting. The approximate solution will be sought as a continuous piecewise linear function in time t and the test space is based on the discontinuous piecewise constant functions. We prove that the proposed time stepping method has the convergence order $$O(\tau ^{1+ \alpha }), \, \alpha \in (0, 1)$$ O ( τ 1 + α ) , α ∈ ( 0 , 1 ) for general sectorial elliptic operators for nonsmooth data by using the Laplace transform method, where $$\tau $$ τ is the time step size. This convergence order is higher than the convergence orders of the popular convolution quadrature methods (e.g., Lubich’s convolution methods) and L-type methods (e.g., L1 method), which have only $$O(\tau )$$ O ( τ ) convergence for the nonsmooth data. Numerical examples are given to verify the robustness of the time discretization schemes with respect to data regularity.

LINEARLY IMPLICIT ENERGY-PRESERVING FOURIER PSEUDOSPECTRAL SCHEMES FOR THE COMPLEX MODIFIED KORTEWEG–DE VRIES EQUATION

2020 ◽  
Vol 62 (3) ◽  
pp. 256-273
Author(s):  
J. L. YAN ◽  
L. H. ZHENG ◽  
L. ZHU ◽  
F. Q. LU
Keyword(s):  
Vries Equation ◽  
Pseudospectral Method ◽  
Time Discretization ◽  
Time Step ◽  
Fourier Pseudospectral Method ◽  
Fully Discrete ◽  
De Vries ◽  
Space Discretization ◽  
Invariant Energy Quadratization ◽  
Fourier Pseudospectral

AbstractWe propose two linearly implicit energy-preserving schemes for the complex modified Korteweg–de Vries equation, based on the invariant energy quadratization method. First, a new variable is introduced and a new Hamiltonian system is constructed for this equation. Then the Fourier pseudospectral method is used for the space discretization and the Crank–Nicolson leap-frog schemes for the time discretization. The proposed schemes are linearly implicit, which is only needed to solve a linear system at each time step. The fully discrete schemes can be shown to conserve both mass and energy in the discrete setting. Some numerical examples are also presented to validate the effectiveness of the proposed schemes.

scholarly journals Artificial Compressibility Methods for the Incompressible Navier–Stokes Equations Using Lowest-Order Face-Based Schemes on Polytopal Meshes

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Riccardo Milani ◽  
Jérôme Bonelle ◽  
Alexandre Ern
Keyword(s):  
Stokes Equations ◽  
Time Discretization ◽  
Second Order ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Artificial Compressibility ◽  
Time Stepping ◽  
Monolithic Approach ◽  
Polytopal Meshes

Abstract We investigate artificial compressibility (AC) techniques for the time discretization of the incompressible Navier–Stokes equations. The space discretization is based on a lowest-order face-based scheme supporting polytopal meshes, namely discrete velocities are attached to the mesh faces and cells, whereas discrete pressures are attached to the mesh cells. This face-based scheme can be embedded into the framework of hybrid mixed mimetic schemes and gradient schemes, and has close links to the lowest-order version of hybrid high-order methods devised for the steady incompressible Navier–Stokes equations. The AC time-stepping uncouples at each time step the velocity update from the pressure update. The performances of this approach are compared against those of the more traditional monolithic approach which maintains the velocity-pressure coupling at each time step. We consider both first-order and second-order time schemes and either an implicit or an explicit treatment of the nonlinear convection term. We investigate numerically the CFL stability restriction resulting from an explicit treatment, both on Cartesian and polytopal meshes. Finally, numerical tests on large 3D polytopal meshes highlight the efficiency of the AC approach and the benefits of using second-order schemes whenever accurate discrete solutions are to be attained.

Linearly implicit energy-preserving Fourier pseudospectral schemes for the complex modified Korteweg-de Vries equation

2021 ◽  
Vol 62 ◽  
pp. 256-273
Author(s):  
J. L. Yan ◽  
L. H. Zheng ◽  
L. Zhu ◽  
F. Q. Lu
Keyword(s):  
Vries Equation ◽  
Pseudospectral Method ◽  
Time Discretization ◽  
Time Step ◽  
Fourier Pseudospectral Method ◽  
Fully Discrete ◽  
De Vries ◽  
Space Discretization ◽  
Invariant Energy Quadratization ◽  
Fourier Pseudospectral

We propose two linearly implicit energy-preserving schemes for the complex modified Korteweg–de Vries equation, based on the invariant energy quadratization method. First, a new variable is introduced and a new Hamiltonian system is constructed for this equation. Then the Fourier pseudospectral method is used for the space discretization and the Crank–Nicolson leap-frog schemes for the time discretization. The proposed schemes are linearly implicit, which is only needed to solve a linear system at each time step. The fully discrete schemes can be shown to conserve both mass and energy in the discrete setting. Some numerical examples are also presented to validate the effectiveness of the proposed schemes.   doi:10.1017/S1446181120000218

Efficient Time-Stepping/Spectral Methods for the Navier-Stokes-Nernst-Planck-Poisson Equations

2017 ◽  
Vol 21 (5) ◽  
pp. 1408-1428 ◽  
Author(s):  
Xiaoling Liu ◽  
Chuanju Xu
Keyword(s):  
Stability Condition ◽  
Time Discretization ◽  
Equation System ◽  
Second Order ◽  
Navier Stokes ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
First Order ◽  
Time Stepping

AbstractThis paper is concerned with numerical methods for the Navier-Stokes-Nernst-Planck-Poisson equation system. The main goal is to construct and analyze some stable time stepping schemes for the time discretization and use a spectral method for the spatial discretization. The main contribution of the paper includes: 1) an useful stability inequality for the weak solution is derived; 2) a first order time stepping scheme is constructed, and the non-negativity of the concentration components of the discrete solution is proved. This is an important property since the exact solution shares the same property. Moreover, the stability of the scheme is established, together with a stability condition on the time step size; 3) a modified first order scheme is proposed in order to decouple the calculation of the velocity and pressure in the fluid field. This new scheme equally preserves the non-negativity of the discrete concentration solution, and is stable under a similar stability condition; 4) a stabilization technique is introduced to make the above mentioned schemes stable without restriction condition on the time step size; 5) finally we construct a second order finite difference scheme in time and spectral discretization in space. The numerical tests carried out in the paper show that all the proposed schemes possess some desirable properties, such as conditionally/unconditionally stability, first/second order convergence, non-negativity of the discrete concentrations, and so on.

scholarly journals Stability Analysis of Decoupled Time-stepping Schemes for the Specialized Conduction System/myocardium Coupled Problem in Cardiology

2017 ◽  
Vol 12 (5) ◽  
pp. 208-239 ◽  
Author(s):  
S. Mani Aouadi ◽  
W. Mbarki ◽  
N. Zemzemi
Keyword(s):  
Time Discretization ◽  
Conduction System ◽  
Splitting Schemes ◽  
Time Step ◽  
Time Stepping ◽  
Time Splitting ◽  
Coupling Terms ◽  
The Stability ◽  
Fully Coupled ◽  
Purkinje Network

The Purkinje network is the rapid conduction system in the heart. It ensures the physiological spread of the electrical wave in the ventricles. In this work, we consider a problem that models the coupling between the Purkinje network and the myocardium. We first prove the stability of the space semi-discretized problem. Then we present four different strategies for solving the Purkinje/ myocardium coupling. The strategies are based on different time discretization of the coupling terms. The first scheme is fully coupled, where the coupling terms are considered implicit. The second and the third schemes are based on Gauss-Seidel time-splitting schemes where one coupling term is considered explicit and the other is implicit. The last is a Jacobi-like time-splitting scheme where both coupling terms are considered explicit. Our main result is the proof of the stability of the three considered schemes under the same restriction on the time step. Moreover, we show that the energy of the problem is slightly affected by the time-splitting schemes. We illustrate the theoretical result by different numerical simulations in 2D. We also conduct 3D simulations using physiologically detailed ionic models.

Fast Conjugate Heat Transfer Simulation of Long Transient Flexible Operations Using Adaptive Time Stepping

2018 ◽  
Vol 140 (9) ◽  
Author(s):  
R. Maffulli ◽  
L. He ◽  
P. Stein ◽  
G. Marinescu
Keyword(s):  
Heat Transfer ◽  
Conjugate Heat Transfer ◽  
Time Derivative ◽  
Computational Cost ◽  
Strong Impact ◽  
Time Step ◽  
Time Stepping ◽  
Adaptive Time ◽  
Fully Coupled ◽  
Adaptive Time Stepping

The emerging renewable energy market calls for more advanced prediction tools for turbine transient operations in fast startup/shutdown cycles. Reliable numerical analysis of such transient cycles is complicated by the disparity in time scales of the thermal responses in fluid and solid domains. Obtaining fully coupled time-accurate unsteady conjugate heat transfer (CHT) results under these conditions would require to march in both domains using the time-step dictated by the fluid domain: typically, several orders of magnitude smaller than the one required by the solid. This requirement has strong impact on the computational cost of the simulation as well as being potentially detrimental to the accuracy of the solution due to accumulation of round-off errors in the solid. A novel loosely coupled CHT methodology has been recently proposed, and successfully applied to both natural and forced convection cases that remove these requirements through a source-term based modeling (STM) approach of the physical time derivative terms in the relevant equations. The method has been shown to be numerically stable for very large time steps with adequate accuracy. The present effort is aimed at further exploiting the potential of the methodology through a new adaptive time stepping approach. The proposed method allows for automatic time-step adjustment based on estimating the magnitude of the truncation error of the time discretization. The developed automatic time stepping strategy is applied to natural convection cases under long (2000 s) transients: relevant to the prediction of turbine thermal loads during fast startups/shutdowns. The results of the method are compared with fully coupled unsteady simulations showing comparable accuracy with a significant reduction of the computational costs.

scholarly journals Goal-Oriented Error Estimation for the Fractional Step Theta Scheme

2014 ◽  
Vol 14 (2) ◽  
pp. 203-230 ◽  
Author(s):  
Dominik Meidner ◽  
Thomas Richter
Keyword(s):  
Time Discretization ◽  
Posteriori Error ◽  
Numerical Dissipation ◽  
Error Estimator ◽  
Fractional Step ◽  
Weighted Residual Method ◽  
Time Stepping ◽  
Backward Euler ◽  
Galerkin Formulation ◽  
Theta Method

Abstract. In this work, we derive a goal-oriented a posteriori error estimator for the error due to time-discretization of nonlinear parabolic partial differential equations by the fractional step theta method. This time-stepping scheme is assembled by three steps of the general theta method, that also unifies simple schemes like forward and backward Euler as well as the Crank–Nicolson method. Further, by combining three substeps of the theta time-stepping scheme, the fractional step theta time-stepping scheme is derived. It possesses highly desired stability and numerical dissipation properties and is second order accurate. The derived error estimator is based on a Petrov–Galerkin formulation that is up to a numerical quadrature error equivalent to the theta time-stepping scheme. The error estimator is assembled as one weighted residual term given by the dual weighted residual method and one additional residual estimating the Galerkin error between time-stepping scheme and Petrov–Galerkin formulation.

Influence of Clearances and Thermal Effects on the Dynamic Behavior of Gear-Hydrodynamic Journal Bearing Systems

2013 ◽  
Vol 135 (6) ◽  
Author(s):  
R. Fargère ◽  
P. Velex
Keyword(s):  
Degrees Of Freedom ◽  
Integration Scheme ◽  
Time Step ◽  
Specific Element ◽  
Normal Contact ◽  
Contact Algorithm ◽  
Gear Teeth ◽  
Proposed Model ◽  
Elastic Couplings ◽  
Definition Of

A global model of mechanical transmissions is introduced which deals with most of the possible interactions between gears, shafts, and hydrodynamic journal bearings. A specific element for wide-faced gears with nonlinear time-varying mesh stiffness and tooth shape deviations is combined with shaft finite elements, whereas the bearing contributions are introduced based on the direct solution of Reynolds' equation. Because of the large bearing clearances, particular attention has been paid to the definition of the degrees-of-freedom and their datum. Solutions are derived by combining a time step integration scheme, a Newton–Raphson method, and a normal contact algorithm in such a way that the contact conditions in the bearings and on the gear teeth are simultaneously dealt with. A series of comparisons with the experimental results obtained on a test rig are given which prove that the proposed model is sound. Finally, a number of results are presented which show that parameters often discarded in global models such as the location of the oil inlet area, the oil temperature in the bearings, the clearance/elastic couplings interactions, etc. can be influential on static and dynamic tooth loading.

Sign in / Sign up

Export Citation Format

Share Document