Entropy-dissipating finite-difference schemes for nonlinear fourth-order parabolic equations

1. Introduction

The design of numerical schemes that preserve the structure of the associated partial differential equations is an important task in numerical mathematics. In this paper, we develop new finite-difference approximations conserving the mass, preserving the positivity, and dissipating the entropy of nonlinear fourth-order parabolic equations of the type

(1)

\partial_{t} u = - J_{x}, J = u^{β} u_{x x x} + a u^{β - 1} u_{x x} u_{x} + b u^{β - 2} u_{x}^{3} in T,

where $T = R / Z$ is the one-dimensional torus, $a$ , $b \in R$ , and $β \geq 0$ , together with the initial condition $u (0) = u^{0}$ in $T$ . We also discuss briefly the discretization of multi-dimensional equations; see Section 4.

A special case of (1) (with $a = b = 0$ ) is the thin-film equation

(2)

\partial_{t} u = - (u^{β} u_{x x x})_{x},

which models the flow of a thin liquid along a solid surface with film height $u (x, t)$ or the thin neck of a Hele-Shaw flow in the lubrication approximation [1]. The multi-dimensional version is given by $\partial_{t} u = - div (u^{β} \nabla Δ u)$ and its discretization will be discussed in Section 4. Another example (with $a = - 2$ , $b = 1$ , $β = 0$ ) is the Derrida–Lebowitz–Speer–Spohn (DLSS) equation

(3)

\partial_{t} u = - 2 (u (\frac{(\sqrt{u})_{x x}}{\sqrt{u}})_{x})_{x} = - (u_{x x x} - 2 \frac{u_{x x} u_{x}}{u} + \frac{u_{x}^{3}}{u^{2}})_{x},

which arises as a scaling limit of interface fluctuations in spin systems [7] and describes the evolution of the electron density $u (x, t)$ in a quantum semiconductor under simplifying assumptions.

Equations (2) and (3) conserve the mass $\int_{T} u (x, t) d x$ and preserve the positivity of the solution, even for the corresponding multi-dimensional versions [6, 14]. They also dissipate the entropy¹¹1Strictly speaking, equation (1) produces entropy and dissipates energy, but the notion entropy dissipation seems to be common in numerical schemes. in the sense that there exists a Lyapunov functional $S (t) = \int_{T} s (u (x, t)) d x$ with entropy density $s (u)$ such that the entropy production $- d S / d t$

provides gradient estimates. For instance, (

2) and (3) dissipate the entropy density

s (u) = s_{α} (u)

, where

	$s_{α} (u) = \frac{u^{α}}{α (α - 1)}$	$for α > 0, α \neq 0, 1,$
(4)		$s_{0} (u) = - log u$	$for α = 0,$
	$s_{1} (u) = u (log u - 1)$	$for α = 1,$

if $\frac{3}{2} \leq α + β \leq 3$ (thin-film equation) and $0 \leq α \leq \frac{3}{2}$ (DLSS equation); see [15, Section 4]. These bounds are optimal. We call $s_{1} (u)$ the Shannon entropy (density), $s_{0} (u)$ the logarithmic entropy, and $s_{α} (u)$ for $α \neq 0, 1$ a Rényi entropy. It holds that $s_{α} (u) \to s_{0} (u)$ pointwise for $α \to 0$ and $s_{α} (u) \to u log u$ pointwise for $α \to 1$ . We prefer to define $s_{1} (u)$ as in (4) to avoid the additional term in $d (u log u) / d u = log u + 1$ which would complicate the subsequent computations.

The proof of the dissipation of the entropy is based on suitable integrations by parts and the chain rule [15]. On the discrete level, we face the problem that the chain rule is generally not available. We are aware of two general strategies.

The first strategy is to exploit the gradient-flow structure of the parabolic equation (if it exists). It involves only one integration by parts, and the discrete chain rule can be formulated by means of suitable mean functions. This idea was elaborated as the Discrete Variational Derivative Method for finite-difference approximations [11]. The gradient-flow formulation (with respect to the $L^{2}$ -Wasserstein metric) yields a natural semi-discretization in time of the evolution using the minimizing movement scheme in finite-dimensional spaces from finite-volume or finite-difference approximations. These techniques were used in [8, 13, 19, 20] for the multi-dimensional DLSS equation. It allows for the proof of entropy dissipation of the Shannon entropy and the Fisher information $\int_{T} (\sqrt{u})_{x}^{2} d x$ , but not for general Rényi entropies, since no Wasserstein gradient-flow formulation seems to be available for these functionals. The thin-film equation with $0 < β < 1$ is shown in [18] to be a gradient flow with respect to a weighted Wasserstein metric. In the work [22], a finite-difference scheme that dissipates the discrete $H^{1}$ norm of the solution to the one-dimensional thin-film equation was analyzed.

The minimizing movement scheme is based on the implicit Euler method. We mention that higher-order time discretizations were investigated too, in the framework of semi-discrete problems; see [4] using the two-step BDF method and [17] using one-leg multistep generalizations. A generic framework for Galerkin methods in space and discontinuous Galerkin methods in time was presented in [9].

The second strategy uses time-continuous Markov chains on finite state spaces. Birth-death processes that define the Markov chain can be interpreted as a finite-volume discretization of a one-dimensional Fokker–Planck equation, and the dissipation of the discrete Shannon entropy can be proved. The nonlinear integrations by parts are reduced to a discrete Bochner-type inequality

[5, 10, 16], which is obtained by identifying the Radon–Nikodym derivative of a measure involving the jump rates of the Markov chain [3, Section 2]. It seems that this idea is restricted to linear diffusion equations.

In this paper, we suggest a third strategy. The idea is to write the flux $J$ as a combination of derivatives of the function $s^{'} (u)$ . This allows for integrations by parts that can be extended to the discrete level and it avoids the application of the chain rule. More precisely, we determine two functions $A$ and $B$ depending on $v := s^{'} (u)$ and its derivatives such that $J = A_{x} - v B_{x}$ . The function $s^{'} (u)$ is known in thermodynamics as the chemical potential, and the formulation of the flux in terms of the chemical potential seems to be natural from a thermodynamic viewpoint. We apply this idea to fourth-order parabolic equations for the first time. It turns out that for $s = s_{α}$ with $α \neq 1$ , we can write

A = \frac{u^{α + β}}{(α - 1)^{2} v} (λ_{1} ξ_{2} + λ_{2} ξ_{1}^{2}), B = \frac{u^{α + β}}{(α - 1)^{2} v^{2}} (λ_{3} ξ_{2} + λ_{4} ξ_{1}^{2}),

where $ξ_{1} = v_{x} / v$ , $ξ_{2} = v_{x x} / v$ , and the coefficients $λ_{i}$ depend on $a$ , $b$ , $α$ , and $β$ ; see (10) and (11). Integrating by parts twice and using equation (1) gives for $S_{α} = \int_{T} s_{α} (u) d x$ :

\frac{d S_{α}}{d t} = \int_{T} J v_{x} d x = - \int_{T} (v_{x x} A - (v v_{x})_{x} B) d x .

The task is to show that the integrand, written as a polynomial in $(ξ_{1}, ξ_{2})$ , is nonnegative for all values of $(ξ_{1}, ξ_{2}) \in R^{2}$ . It follows from the product rule $(v v_{x})_{x} = v v_{x x} + v_{x}^{2}$ that, for $α \neq 1$ ,

	$\frac{d S_{α}}{d t}$	$= - \int_{T} (ξ_{2} v A - (ξ_{2} + ξ_{1}^{2}) v^{2} B) d x$
(5)			$= - \int_{T} \frac{u^{α + β}}{(α - 1)^{2}} ((λ_{1} - λ_{3}) ξ_{2}^{2} + (λ_{2} - λ_{3} - λ_{4}) ξ_{2} ξ_{1}^{2} - λ_{4} ξ_{1}^{4}) d x .$

Under certain conditions on the parameters, we expect that the integrand is bounded from below by $u^{α + β} (ξ_{2}^{2} + ξ_{1}^{4})$ , up to a factor, which yields some gradient estimates.

On the discrete level, we imitate this idea: The flux $J = A_{x} - v B_{x}$ and the variables $ξ_{1}$ , $ξ_{2}$ of the polynomials $A$ and $B$ is discretized by central finite differences. For this, let $T_{h} = T / (h Z)$ be a discrete torus with grid size $h > 0$ and define the scheme

(6)

\partial_{t} u_{i} = - \frac{1}{h} (J_{i + 1 / 2} - J_{i - 1 / 2}), J_{i + 1 / 2} = \frac{1}{h} (A_{i + 1} - A_{i}) - \frac{1}{2 h} (v_{i + 1} + v_{i}) (B_{i + 1} - B_{i}),

with the initial condition $u_{i} (0) = u^{0} (i)$ for $i \in T_{h}$ , where $u_{i} = u (i)$ , $v_{i} = s^{'} (u_{i})$ , and $A_{i}$ and $B_{i}$ are the polynomials $A$ and $B$ , evaluated at $i \in T_{h}$ , respectively; see (19). We show below that with the discrete entropy $S_{α}^{h} = h \sum_{i \in T_{h}} s_{α} (u)$ , the discrete analog of (5) becomes

\frac{d S_{α}^{h}}{d t} = - h \sum i \in T_{h} (ξ_{2, i} v_{i} A_{i} - (v_{i} ξ_{2, i} + ξ_{1, i}^{2}) v_{i}^{2} B_{i}),

where $ξ_{1, i}$ , $ξ_{2, i}$ approximate $v_{x} (i) / v (i)$ , $v_{x x} (i) / v (i)$ , respectively. This yields exactly the polynomial of the continuous case. Thus, we obtain the same conditions on the parameters as for the continuous equation.

We still need a discrete analog of the product rule $(v v_{x})_{x} = v v_{x x} + v_{x}^{2}$ to conclude. This is done by carefully choosing $ξ_{1, i}$ and $ξ_{2, i}$ . Definition (20) ensures that $v_{i}^{2} (ξ_{2, i} + ξ_{1, i}^{2}) = (v_{i + 1}^{2} - 2 v_{i} + v_{i - 1}^{2}) / (2 h^{2})$ which approximates $\frac{1}{2} (v^{2})_{x x} = v v_{x x} + v_{x}^{2}$ . This choice is used in the central scheme (6) for $J_{i + 1 / 2}$ . Noncentral schemes require different definitions of $ξ_{1, i}$ and $ξ_{2, i}$ ; see Remark 7.

A drawback of our technique is that scheme (6) depends on the entropy to be dissipated. The scheme does not dissipate all admissible entropies. In applications, however, one usually wants to dissipate only that entropy which is physically relevant.

Our main results can be sketched as follows:

Lyapunov functional: Let $α \geq 0$ , $α \neq 1$ and assume that

(7) $K (α, β) := - 2 α^{2} + (3 a - 4 β + 9) α - 2 β^{2} + (3 a + 9) β - 9 (a + b + 1) \geq 0.$

Then the continuous entropy $S_{α}$ and the discrete entropy $S_{α}^{h}$ are dissipated in the sense that $d S_{α} / d t \leq 0$ and $d S_{α}^{h} / d t \leq 0$ , i.e., $S_{α}$ and $S_{α}^{h}$ are Lyapunov functionals for $u (t)$ and $u_{i} (t)$ , respectively; see Theorems 2 and 5. Condition (7) is optimal for the thin-film and DLSS equations.
Entropy dissipation: Let $α \geq 0$ , $α \neq 1$ and assume that $K (α, β) > 0$ . Then there exists a constant $c (α, β) > 0$ such that for all $t > 0$ ,

(8) $\frac{d S_{α}}{d t} + c (α, β) \int_{T} u^{α + β} (ξ_{2}^{2} + ξ_{1}^{4}) d x$ $\leq 0,$

(9) $\frac{d S_{α}^{h}}{d t} + c (α, β) h \sum i \in T_{h} {¯ u}_{i}^{α + β} (ξ_{2, i}^{2} + ξ_{1, i}^{4})$ $\leq 0,$

where $ξ_{1} = v_{x} / v$ , $ξ_{2} = v_{x x} / v$ , $v = s_{α} (u)$ , ${¯ u}_{i}$ is an arbitrary average of $(u_{i})$ , and $ξ_{1, i}$ , $ξ_{2, i}$ are defined in (20). This result is proved in Theorems 2 and 5.
Case $α = 1$ : For the case $α = 1$ , we need the formulation $J = A_{x} - w B_{x}$ instead of $J = A_{x} - v B_{x}$ , where $w = s_{0}^{'} (u) = - u^{- 1}$ and $v = log u$ , since $J$ generally does not depend on the logarithm. For details, see Proposition 8.
Case $α = 0$ : We show that scheme (6) with $α = 0$ possesses global positive solutions. This result is a consequence of the discrete entropy inequality and mass conservation, which imply that $h \sum_{i \in T_{h}} (u_{i} (t) - log u_{i} (t))$ is bounded for all $t > 0$ . Consequently, $u_{i} (t) - log u_{i} (t)$ is bounded for all $i \in T_{h}$ and $t > 0$ , and since the function $s \mapsto s - log s$ tends to infinity if either $s \to 0$ or $s \to \infty$ , this proves that $u_{i} (t)$ is bounded from below and above. We refer to Proposition 6 for details.
Multi-dimensional case: In principle, the multi-dimensional case can be treated using functions $A$ and $B$ with many variables. Practically, however, the computations are becoming too involved and it may be unclear how to discretize mixed derivatives. One idea to overcome this issue is to use scalar variables only, like $ξ_{1} = | \nabla v | / v$ and $ξ_{2} = Δ v / v$ . This allows us to treat the multi-dimensional thin-film equation; see Proposition 10.

The paper is organized as follows. We prove the continuous entropy inequality (8) in Section 2 and the discrete entropy inequality (9) in Section 3. A scheme for the multi-dimensional thin-film equation is proposed and analyzed in Section 4. Numerical simulations for the thin-film and DLSS equations in one space dimension and for a thin-film equation in two space dimensions are presented in Section 5.

2. General continuous equation

To prepare the discretization, we need to analyze the entropy dissipation properties of the continuous equation (1). We show first that $J$ can be written as $A_{x} - v B_{x}$ with $v = s_{α}^{'} (u)$ and functions $A$ and $B$ which depend on $v$ , $v_{x}$ , and $v_{x x}$ .

Lemma 1.

Let $J$ be given as in (1) and $s_{α}$ as in (4) and let $α \geq 0$ satisfy $α \neq 1$ . Then $J = A_{x} - v B_{x}$ , where

	$A$	$= u^{α + β} (λ_{1} u^{2 - 2 α} (s_{α}^{'} (u))_{x x} + λ_{2} (α - 1) u^{3 - 3 α} (s_{α}^{'} (u))_{x}^{2})$
		$= \frac{u^{α + β}}{(α - 1)^{2} v} (λ_{1} \frac{}{v_{x x}} v + λ_{2} (\frac{v_{x}}{v})^{2}),$
	$B$	$= u^{α + β} (λ_{3} (α - 1) u^{3 - 3 α} (s_{α}^{'} (u))_{x x} + λ_{4} (α - 1)^{2} u^{4 - 4 α} (s_{α}^{'} (u))_{x}^{2})$
		$= \frac{u^{α + β}}{(α - 1)^{2} v^{2}} (λ_{3} \frac{v_{x x}}{v} + λ_{4} (\frac{v_{x}}{v})^{2})$

and

	$λ_{1}$	$= \frac{- 2 α^{2} + (β + 5) α + a β - β^{2} - a - 2 b - 3}{(α - 1) (β - 2 α + 3)} + \frac{2 (α - 1)}{β - 2 α + 3} λ_{4},$
(10)		$λ_{2}$	$= \frac{2 α^{2} - (a + 7) α + 2 a + b + 6}{(α - 1) (β - 2 α + 3)} + \frac{β - 3 α + 4}{β - 2 α + 3} λ_{4},$
	$λ_{3}$	$= \frac{(a + 1) β - β^{2} - a - 2 b}{(1 - α) (β - 2 α + 3)} + \frac{2 (α - 1)}{β - 2 α + 3} λ_{4},$

with $λ_{4} \in R$ being a free parameter.

The lemma shows in particular that the formulation $J = A_{x} - v B_{x}$ is not unique. This fact is used to optimize later the range of admissible parameters $α$ , $β$ , $a$ , and $b$ .

Proof.

Let $α > 0$ with $α \neq 1$ . A direct computation yields

	$A$	$= λ_{1} u^{β} u_{x x} + ((λ_{1} + λ_{2}) α - 2 λ_{1} - λ_{2}) u^{β - 1} u_{x}^{2},$
	$B$	$= λ_{3} (α - 1) u^{β - α + 1} u_{x x} + ((λ_{3} + λ_{4}) α - 2 λ_{3} - λ_{4}) u^{β - α} u_{x}^{2} .$

Inserting these expressions into $A_{x} - v B_{x}$ and simplifying leads to

	$A_{x}$	$- \frac{u^{α - 1}}{α - 1} B_{x} = (λ_{1} - λ_{3}) u^{β} u_{x x x}$
		$+ ((2 λ_{1} + 2 λ_{2} - λ_{3} - 2 λ_{4}) α + (λ_{1} - λ_{3}) β - 4 λ_{1} - 2 λ_{2} + 3 λ_{3} + 2 λ_{4}) u^{β - 1} u_{x x} u_{x}$
		$+ ((λ_{3} + λ_{4}) α^{2} + (λ_{1} + λ_{2} - λ_{3} - λ_{4}) α β - (λ_{1} + λ_{2} + 2 λ_{3} + λ_{4}) α$
		$+ (- 2 λ_{1} - λ_{2} + 2 λ_{3} + λ_{4}) β + 2 λ_{1} + λ_{2}) u^{β - 2} u_{x}^{3} .$

We identify the coefficients with those in the expression (1) for $J$ :

	$1$	$= λ_{1} - λ_{3},$
	$a$	$= (2 λ_{1} + 2 λ_{2} - λ_{3} - 2 λ_{4}) α + (λ_{1} - λ_{3}) β - 4 λ_{1} - 2 λ_{2} + 3 λ_{3} + 2 λ_{4},$
	$b$	$= (λ_{3} + λ_{4}) α^{2} + (λ_{1} + λ_{2} - λ_{3} - λ_{4}) α β - (λ_{1} + λ_{2} + 2 λ_{3} + λ_{4}) α$
		$+ (- 2 λ_{1} - λ_{2} + 2 λ_{3} + λ_{4}) β + 2 λ_{1} + λ_{2} .$

The general solution of this linear system for $(λ_{1}, λ_{2}, λ_{3}, λ_{4})$ , with free parameter $λ_{4} \in R$ , gives (10).

Next, let $α = 0$ . The ansatz for $A$ and $B$ becomes

	$A$	$= u^{β} (λ_{1} u^{2} (- u^{- 1})_{x x} - λ_{2} u^{3} (- u^{- 1})_{x}^{2}),$
	$B$	$= u^{β} (- λ_{3} u^{3} (- u^{- 1})_{x x} + λ_{4} u^{4} (- u^{- 1})_{x}^{2}) .$

Then

	$A_{x} - v B_{x}$	$= (λ_{1} - λ_{3}) u^{β} u_{x x x} + ((λ_{1} - λ_{3}) β - 4 λ_{1} - 2 λ_{2} + 3 λ_{3} + 2 λ_{4}) u^{β - 1} u_{x x} u_{x}$
		$+ ((- 2 λ_{1} - λ_{2} + 2 λ_{3} + λ_{4}) β + 2 λ_{1} + λ_{2}) u^{β - 2} u_{x}^{3} .$

Identifying the coefficients with those in (1) again gives a linear system for the parameters $(λ_{1}, λ_{2}, λ_{3}, λ_{4})$ . The general solution reads as

	$λ_{1}$	$= \frac{β^{2} - a β + a + 2 b + 3}{β + 3} - \frac{2}{β + 3} λ_{4},$
(11)		$λ_{2}$	$= - \frac{2 a + b + 6}{β + 3} + \frac{β + 4}{β + 3} λ_{4},$
	$λ_{3}$	$= \frac{β^{2} - a β - β + a + 2 b}{β + 3} - \frac{2}{β + 3} λ_{4},$

These expressions are the same as (10) with $α = 0$ . ∎

For the following theorem, we recall definition (4) of $s_{α}$ and set $S_{α} (u) = \int_{T} s_{α} (u) d x$ .

Theorem 2.

Let $u$ be a smooth positive solution to (1) and let $α \geq 0$ . If $K (α, β) \geq 0$ (see definition (7)), then $S_{α}$ is a Lyapunov functional, i.e. $d S_{α} / d t \leq 0$ for all $t > 0$ . If $K (α, β) > 0$ then there exists $c (α, β) > 0$ such that for all $t > 0$ ,

\frac{d S_{α}}{d t} + c (α, β) \int_{T} u^{α + β} (u^{2 (1 - α)} (u^{α - 1})_{x x}^{2} + u^{4 (1 - α)} (u^{α - 1})_{x}^{4}) d x \leq 0

and there exists another constant $C (α, β) > 0$ such that for all $t > 0$ ,

\frac{d S_{α}}{d t} + C (α, β) \int_{T} u^{α + β} (u^{- 2} u_{x x}^{2} + u^{- 4} u_{x}^{4}) d x \leq 0.

Proof.

Let first $α \geq 0$ with $α \neq 1$ . We calculate the time derivative of the entropy, using integration by parts twice:

	$\frac{d S_{α}}{d t}$	$= \int_{T} s_{α}^{'} (u) \partial_{t} u d x = - \int_{T} v J_{x} d x = \int_{T} J v_{x} d x$
		$= \int_{T} (A_{x} - v B_{x}) v_{x} d x = - \int_{T} (A v_{x x} - (v v_{x x} + v_{x}^{2}) B) d x$

(12)

where

(13)

P (ξ_{1}, ξ_{2}) = (λ_{1} - λ_{3}) ξ_{2}^{2} + (λ_{2} - λ_{3} - λ_{4}) ξ_{2} ξ_{1} - λ_{4} ξ_{1}^{4} .

The right-hand side of (12) is nonpositive if $P (ξ_{1}, ξ_{2}) \geq 0$ for all $(ξ_{1}, ξ_{2}) \in R^{2}$ . Taking into account that $λ_{1} - λ_{3} = 1$ , this is the case if and only if

- 4 λ_{4} - (λ_{2} - λ_{3} - λ_{4})^{2} \geq 0.

In view of definition (10) of $λ_{2}$ and $λ_{3}$ , we may interpret $λ_{4}$ as a free parameter und $λ_{2} = λ_{2} (λ_{4})$ and $λ_{3} = λ_{3} (λ_{4})$ as affine functions in $λ_{4}$ . Therefore, we require that

(14)

f (λ_{4}) := - 4 λ_{4} - (λ_{2} (λ_{4}) - λ_{3} (λ_{4}) - λ_{4})^{2} \geq 0.

We choose the optimal value of $λ_{4}$ by computing the critical value of $f$ :

	$0 = f^{'} (λ_{4})$	$= \frac{2 (- 2 α^{2} + (- 3 a + 8 β + 3) α - 3 a β + β^{2} + 9 (a + b) - 15 β)}{(β - 2 α + 3)^{2}}$
		$- \frac{18 (α - 1)^{2}}{(β - 2 α + 3)^{2}} λ_{4} .$

This yields

(15)

λ_{4} = \frac{- 2 α^{2} + (- 3 a + 8 β + 3) α - 3 a β + β^{2} + 9 (a + b) - 15 β}{9 (α - 1)^{2}} .

Inserting $λ_{4}$ in (14) leads to

0 \leq f (λ_{4}) = \frac{4}{9 (α - 1)^{2}} (- 2 α^{2} + (3 a - 4 β + 9) α - 2 β^{2} + (3 a + 9) β - 9 (a + b + 1)),

which is equivalent to $K (α, β) \geq 0$ , see (7).

If $K (α, β) > 0$ , there exists $c_{0} (α, β) > 0$ such that for all $(ξ_{1}, ξ_{2}) \in R^{2}$ , $P (ξ_{1}, ξ_{2}) \geq c_{0} (α, β) (ξ_{1}^{2} + ξ_{1}^{4})$ . Inserting this information in (12), we infer that

	$\frac{d S_{α}}{d t}$	$\leq - c_{0} (α, β) \int_{T} \frac{u^{α + β}}{(α - 1)^{2}} [(\frac{v_{x x}}{v})^{2} + (\frac{v_{x}}{v})^{4}] d x$
		$= - c_{0} (α, β) \int_{T} u^{α + β - 2} [(\frac{u_{x x}}{u})^{2} + 2 (α - 2) \frac{u_{x x}}{u} (\frac{u_{x}}{u})^{2} + (2 α^{2} - 6 α + 5) (\frac{u_{x}}{u})^{4}] d x .$

The discriminant equals

1 \cdot (2 α^{2} - 6 α + 5) - (α - 2)^{2} = (α - 1)^{2},

and it is positive for all $α \neq 1$ . Therefore, there exists $k (α) > 0$ such that

\frac{d S_{α}}{d t} \leq - c_{0} (α, β) k (α) \int_{T} u^{α + β} ((\frac{u_{x x}}{u})^{2} + (\frac{u_{x}}{} u)^{4}) d x,

and this gives the conclusion for $c (α, β) = c_{0} (α, β) / (α - 1)^{2}$ and $C (α, β) = c_{0} (α, β) k (α)$ when $α \neq 1$ .

It remains to analyze the case $α = 1$ . Here, we cannot formulate the flux as $J = A_{x} - v B_{x}$ with $v = s_{1}^{'} (u) = log u$ , since $J$ does not contain logarithmic terms. Our idea is to write $J = A_{x} - w B_{x}$ with $w = - u^{- 1}$ and functions $A$ and $B$ that depend on $w$ , $w_{x}$ , and $w_{x x}$ . The time derivative of the entropy $S_{1}$ can be written in terms of $w$ and its derivatives only, since the logarithmic term $v = log u$ only appears with its derivatives.

The formulation $J = A_{x} - w B_{x}$ corresponds to the expression used for $α = 0$ . In fact, we have

	$A$	$= u^{β} (λ_{1} u^{2} (- u^{- 1})_{x x} - λ_{2} u^{3} (- u^{- 1})_{x}^{2}) = \frac{u^{β}}{w} (λ_{1} \frac{w_{x x}}{w} + λ_{2} (\frac{w_{x}}{w})^{2}),$
	$B$	$= u^{β} (- λ_{3} u^{3} (- u^{- 1})_{x x} + λ_{4} u^{4} (- u^{- 1})_{x}^{2}) = \frac{u^{β}}{w^{2}} (λ_{3} \frac{w_{x x}}{w} + λ_{4} (\frac{w_{x}}{w})^{2}),$

where $λ_{1}, λ_{2}, λ_{3}$ are given by (11) and $λ_{4}$ is a free parameter. As before, the time derivative becomes

\frac{d S_{1}}{d t} = \int_{T} \partial_{t} u log u d x = - \int_{T} J_{x} v d x = \int_{T} J v_{x} d x = - \int_{T} (A v_{x x} - B (v_{x} w)_{x}) d x .

Set $ξ_{1} = w_{x} / w$ and $ξ_{2} = w_{x x} / w$ . Since $v_{x x} = - w_{x x} / w + (w_{x} / w)^{2} = - ξ_{2} + ξ_{1}^{2}$ and $(v_{x} w)_{x} = - w_{x x} = - w ξ_{2}$ , we obtain

	$\frac{d S_{1}}{d t}$	$= - \int_{T} \frac{u^{β}}{w} ((λ_{1} ξ_{2} + λ_{2} ξ_{1}^{2}) (- ξ_{2} + ξ_{1}^{2}) + (λ_{3} ξ_{2} + λ_{4} ξ_{1}^{2}) ξ_{2}) d x$
		$= - \int_{T} u^{β + 1} (- (λ_{1} ξ_{2} + λ_{2} ξ_{1}^{2}) (- ξ_{2} + ξ_{1}^{2}) - (λ_{3} ξ_{2} + λ_{4} ξ_{1}^{2}) ξ_{2}) d x$
		$= - \int_{T} u^{β + 1} P_{1} (ξ_{1}, ξ_{2}) d x,$

where

(16)

P_{1} (ξ_{1}, ξ_{2}) = (λ_{1} - λ_{3}) ξ_{2}^{2} + (- λ_{1} + λ_{2} - λ_{4}) ξ_{2} ξ_{1}^{2} - λ_{2} ξ_{1}^{4} .

Observe that $λ_{1} - λ_{3} = 1$ . Thus, $d S_{1} / d t \geq 0$ if the polynomial $P_{1}$ is nonnegative for all $(ξ_{1}, ξ_{2}) \in R^{2}$ , which is the case if and only if

f_{1} (λ_{4}) := - 4 λ_{2} (λ_{4}) - (- λ_{1} (λ_{4}) + λ_{2} (λ_{4}) - λ_{4})^{2} \geq 0.

We choose the optimal value $λ_{4}$ from $f_{1}^{'} (λ_{4}) = 0$ , i.e.

λ_{4} = \frac{1}{9} (β^{2} - (3 a + 14) β + 9 (a + b) + 3) .

Then $f (λ_{4}) \geq 0$ is equivalent to $K (1, β) \geq 0$ .

Finally, if $K (1, β) > 0$ , we have $P_{1} (ξ_{1}, ξ_{2}) \geq c_{0} (1, β) (ξ_{2}^{2} + ξ_{1}^{4})$ and consequently,

	$\frac{d S_{1}}{d t}$
		$= - c_{0} (1, β) \int_{T} u^{β + 1} (u^{- 2} u_{x x}^{2} - 4 u^{- 3} u_{x x} u_{x}^{2} + 5 u^{- 4} u_{x}^{4}) d x .$

The discriminant is positive and there exists $k (1) > 0$ such that $ξ_{2}^{2} - 4 ξ_{2} ξ_{1}^{2} + 5 ξ_{1}^{4} \geq k (1) (ξ_{2}^{2} + ξ_{1}^{4})$ . We infer that

\frac{d S_{1}}{d t} \leq - c_{0} (1, β) k (1) \int_{T} u^{β + 1} (u^{- 2} u_{x x}^{2} + u^{- 4} u_{x}^{4}) d x,

which concludes the proof with $c (1, β) = c_{0} (1, β)$ and $C (1, β) = c_{0} (1, β) k (1)$ . ∎

Remark 3 (Examples).

The DLSS equation corresponds to (1) if $a = - 2$ , $b = 0$ , and $β = 0$ . Then condition (7) becomes $0 < α \leq \frac{3}{2}$ , which is the optimal interval. Choosing $a = b = 0$ , we obtain the thin-film equation, and condition (7) is equivalent to $\frac{3}{2} \leq α + β \leq 3$ , which again is the optimal parameter range. ∎

Remark 4 (Systematic integration by parts).

The result of Theorem 2 can be also derived by the method of systematic integration by parts of [15]. Our proof is taylored in such a way that it can be directly “translated” to the discrete level. Indeed, the method of [15] needs several chain rules that are not available on the discrete level and our technique needs to be used.

Still, systematic integration by parts and our strategy are strongly related. Systematic integration by parts means that we are adding so-called integration-by-parts formulas whose derivatives vanishe. In this way, we can derive (12),

\frac{d S_{α}}{d t} = \int_{T} (A_{x} - v B_{x}) v_{x} d x + c_{1} \int_{T} (A v_{x} - B v v_{x})_{x} d x = - \int_{T} (A v_{x x} - B (v v_{x})_{x}) d x,

by choosing $c_{1} =$

Entropy-dissipating finite-difference schemes for nonlinear fourth-order parabolic equations

Authors

Energy-consistent finite difference schemes for compressible hydrodynamics and magnetohydrodynamics using nonlinear filtering

Extending nonstandard finite difference schemes rules to systems of nonlinear ODEs with constant coefficients

A generalized optimal fourth-order finite difference scheme for a 2D Helmholtz equation with the perfectly matched layer boundary condition

Sixth order weighted essentially non-oscillatory schemes with Z-type nonlinear weighting procedure for nonlinear degenerate parabolic equations

Adaptive Third Order Adams-Bashforth Time Stepping Scheme for 2D Extended Boussinesq Equations

OPESCI-FD: Automatic Code Generation Package for Finite Difference Models

Truly multi-dimensional all-speed schemes for the Euler equations on Cartesian grids

1. Introduction

2. General continuous equation

Lemma 1.

Proof.

Theorem 2.

Proof.

Remark 3 (Examples).

Remark 4 (Systematic integration by parts).

(8)		$\frac{d S_{α}}{d t} + c (α, β) \int_{T} u^{α + β} (ξ_{2}^{2} + ξ_{1}^{4}) d x$	$\leq 0,$
(9)		$\frac{d S_{α}^{h}}{d t} + c (α, β) h \sum i \in T_{h} {¯ u}_{i}^{α + β} (ξ_{2, i}^{2} + ξ_{1, i}^{4})$	$\leq 0,$

Entropy-dissipating finite-difference schemes for nonlinear fourth-order parabolic equations

Authors

Related Research

Energy-consistent finite difference schemes for compressible hydrodynamics and magnetohydrodynamics using nonlinear filtering

Extending nonstandard finite difference schemes rules to systems of nonlinear ODEs with constant coefficients

A generalized optimal fourth-order finite difference scheme for a 2D Helmholtz equation with the perfectly matched layer boundary condition

Sixth order weighted essentially non-oscillatory schemes with Z-type nonlinear weighting procedure for nonlinear degenerate parabolic equations

Adaptive Third Order Adams-Bashforth Time Stepping Scheme for 2D Extended Boussinesq Equations

OPESCI-FD: Automatic Code Generation Package for Finite Difference Models

Truly multi-dimensional all-speed schemes for the Euler equations on Cartesian grids

This week in AI

1. Introduction

2. General continuous equation

Lemma 1.

Proof.

Theorem 2.

Proof.

Remark 3 (Examples).

Remark 4 (Systematic integration by parts).