discontinuous flux
Recently Published Documents


TOTAL DOCUMENTS

99
(FIVE YEARS 14)

H-INDEX

19
(FIVE YEARS 2)

Mathematics ◽  
2021 ◽  
Vol 9 (10) ◽  
pp. 1079
Author(s):  
Tingting Xiang ◽  
Guodong Wang ◽  
Suping Zhang

A new modified Engquist–Osher-type flux-splitting scheme is proposed to approximate the scalar conservation laws with discontinuous flux function in space. The fact that the discontinuity of the fluxes in space results in the jump of the unknown function may be the reason why it is difficult to design a high-order scheme to solve this hyperbolic conservation law. In order to implement the WENO flux reconstruction, we apply the new modified Engquist–Osher-type flux to compensate for the discontinuity of fluxes in space. Together the third-order TVD Runge–Kutta time discretization, we can obtain the high-order accurate scheme, which keeps equilibrium state across the discontinuity in space, to solve the scalar conservation laws with discontinuous flux function. Some examples are given to demonstrate the good performance of the new high-order accurate scheme.


Author(s):  
Adrian M Ruf

Abstract We prove that adapted entropy solutions of scalar conservation laws with discontinuous flux are stable with respect to changes in the flux under the assumption that the flux is strictly monotone in $u$ and the spatial dependency is piecewise constant with finitely many discontinuities. We use this stability result to prove a convergence rate for the front tracking method—a numerical method that is widely used in the field of conservation laws with discontinuous flux. To the best of our knowledge, both of these results are the first of their kind in the literature on conservation laws with discontinuous flux. We also present numerical experiments verifying the convergence rate results and comparing numerical solutions computed with the front tracking method to finite volume approximations.


Author(s):  
Jayesh Badwaik ◽  
Christian Klingenberg ◽  
Nils Henrik Risebro ◽  
Adrian M Ruf

We consider conservation laws with discontinuous flux where the initial datum, the flux function, and the discontinuous spatial dependency coefficient are subject to randomness. We establish a notion of random adapted entropy solutions to these equations and prove well-posedness provided that the spatial dependency coefficient is piecewise constant with finitely many discontinuities. In particular, the setting under consideration allows the flux to change across finitely many points in space whose positions are uncertain. We propose a single- and multilevel Monte Carlo method based on a finite volume approximation for each sample. Our analysis includes convergence rate estimates of the resulting Monte Carlo and multilevel Monte Carlo finite volume methods as well as error versus work rates showing that the multilevel variant outperforms the single-level method in terms of efficiency. We present numerical experiments motivated by two-phase reservoir simulations for reservoirs with varying geological properties.


2020 ◽  
Vol 52 (3) ◽  
pp. 322-338
Author(s):  
Nasrin Okhovati ◽  
Mohammad Izadi

In this paper we propose an explicit predictor-corrector finite difference scheme to numerically solve one-dimensional conservation laws with discontinuous flux function appearing in various physical model problems, such as traffic flow and two-phase flow in porous media. The proposed method is based on the second-order MacCormack finite difference scheme and the solution is obtained by correcting first-order schemes. It is shown that the order of convergence is quadratic in the grid spacing for uniform grids when applied to problems with discontinuity. To illustrate some properties of the proposed scheme, numerical results applied to linear as well as non-linear problems are presented.


Author(s):  
M. Graf ◽  
M. Kunzinger ◽  
D. Mitrovic ◽  
D. Vujadinovic

AbstractWe prove existence and uniqueness of a solution to the Cauchy problem corresponding to the dynamics capillarity equation $$\begin{aligned} {\left\{ \begin{array}{ll} \partial _t u_{\varepsilon ,\delta } +\mathrm {div} {\mathfrak f}_{\varepsilon ,\delta }(\mathbf{x}, u_{\varepsilon ,\delta })=\varepsilon \Delta u_{\varepsilon ,\delta }+\delta (\varepsilon ) \partial _t \Delta u_{\varepsilon ,\delta }, \ \ \mathbf{x} \in M, \ \ t\ge 0\\ u|_{t=0}=u_0(\mathbf{x}). \end{array}\right. } \end{aligned}$$ ∂ t u ε , δ + div f ε , δ ( x , u ε , δ ) = ε Δ u ε , δ + δ ( ε ) ∂ t Δ u ε , δ , x ∈ M , t ≥ 0 u | t = 0 = u 0 ( x ) . Here, $${{\mathfrak {f}}}_{\varepsilon ,\delta }$$ f ε , δ and $$u_0$$ u 0 are smooth functions while $$\varepsilon $$ ε and $$\delta =\delta (\varepsilon )$$ δ = δ ( ε ) are fixed constants. Assuming $${{\mathfrak {f}}}_{\varepsilon ,\delta } \rightarrow {{\mathfrak {f}}}\in L^p( {\mathbb {R}}^d\times {\mathbb {R}};{\mathbb {R}}^d)$$ f ε , δ → f ∈ L p ( R d × R ; R d ) for some $$1<p<\infty $$ 1 < p < ∞ , strongly as $$\varepsilon \rightarrow 0$$ ε → 0 , we prove that, under an appropriate relationship between $$\varepsilon $$ ε and $$\delta (\varepsilon )$$ δ ( ε ) depending on the regularity of the flux $${{\mathfrak {f}}}$$ f , the sequence of solutions $$(u_{\varepsilon ,\delta })$$ ( u ε , δ ) strongly converges in $$L^1_{loc}({\mathbb {R}}^+\times {\mathbb {R}}^d)$$ L loc 1 ( R + × R d ) toward a solution to the conservation law $$\begin{aligned} \partial _t u +\mathrm {div} {{\mathfrak {f}}}(\mathbf{x}, u)=0. \end{aligned}$$ ∂ t u + div f ( x , u ) = 0 . The main tools employed in the proof are the Leray–Schauder fixed point theorem for the first part and reduction to the kinetic formulation combined with recent results in the velocity averaging theory for the second. These results have the potential to generate a stable semigroup of solutions to the underlying scalar conservation laws different from the Kruzhkov entropy solutions concept.


2020 ◽  
Vol 21 (1) ◽  
pp. 21
Author(s):  
Isamara L. N. Araujo ◽  
Panters Rodríguez-Bermúdez ◽  
Yoisell Rodríguez-Núñez

In this work we study two-phase flow with gravity either in 1-rock homogeneous media or 2-rocks composed media, this phenomenon can be modeled by a non-linear scalar conservation law with continuous flux function or discontinuous flux function, respectively. Our study is essentially from a numerical point of view, we apply the new Lagrangian-Eulerian finite difference method developed by Abreu and Pérez  and the Lax-Friedrichs classic method to obtain numerical entropic solutions. Comparisons between numerical and analytical solutions show the efficiency of the methods even for discontinuous flux function.


Sign in / Sign up

Export Citation Format

Share Document