A finite volume method for scalar conservation laws with stochastic time–space dependent flux functions

We propose a new finite volume method for scalar conservation laws with stochastic time–space dependent flux functions. The stochastic effects appear in the flux function and can be interpreted as a random manner to localize the discontinuity in the time–space dependent flux function. The location of the interface between the fluxes can be obtained by solving a system of stochastic differential equations for the velocity fluctuation and displacement variable. In this paper we develop a modified Rusanov method for the reconstruction of numerical fluxes in the finite volume discretization. To solve the system of stochastic differential equations for the interface we apply a second-order Runge–Kutta scheme. Numerical results are presented for stochastic problems in traffic flow and two-phase flow applications. It is found that the proposed finite volume method offers a robust and accurate approach for solving scalar conservation laws with stochastic time–space dependent flux functions.     

Two positivity preserving flux limited,second‐order numerical methods for a haptotaxis model

Two numerical methods for a one‐dimensional haptotaxis model, which exploit the use of van Leer flux limiter, are developed and analyzed. Sufficient conditions time step size and flux limiting are given for such formulation to ensure the non‐negativity of the discrete solution and second‐order accuracy in space. Another advantage is that we avoid solving large nonlinear systems of algebraic equations. The discrete preservation of total conservation of cell density, concentration, and logarithmic density is also verified for the numerical solution. Numerical results concerning accuracy, convergence rate, positivity, and conservation properties are presented and discussed. Similar approach could be applied efficiently in the corresponding two‐ and three‐dimensional problems. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

More on P-Stable Convex Sets in Banach Spaces

We study the asymptotic behavior and limit distributions for sums S _n =b_n ^-1 _i=1 ⁿ _i,where _i, i 1, are i.i.d. random convex compact (cc) sets in a given separable Banach space B and summation is defined in a sense of Minkowski. The following results are obtained: (i) Series (LePage type) and Poisson integral representations of random stable cc sets in B are established; (ii) The invariance principle for processes S _n(t) =b_n ^-1 _i=1 ^[nt] _i, t[0, 1], and the existence of p-stable cc Levy motion are proved; (iii) In the case, where _i are segments, the limit of S _n is proved to be countable zonotope. Furthermore, if B = R ^d , the singularity of distributions of two countable zonotopes Y_{p 1}, ₁,Y_{p 2}, ₂, corresponding to values of exponents p ₁, p ₂ and spectral measures ₁, ₂, is proved if either p ₁ p ₂ or _{1 2}; (iv) Some new simple estimates of parameters of stable laws in R ^d , based on these results are suggested.     

Marcinkiewicz–Zygmund Strong Laws for Infinite Variance Time Series

We study rate of convergence in the strong law of large numbers for finite and infinite variance time series in both contexts of weak and strong dependence.     

Pathwise superhedging for time-dependent barrier options on càdlàg paths—Finite or infinite tradeable European,One-Touch,lookback or forward starting options

We establish pathwise duality using simple predictable trading strategies for the robust hedging problem associated with a barrier option whose payoff depends on the terminal level and the infimum of a càdlàg strictly positive stock price process, given tradeable European options at all strikes at a single maturity. The result allows for a significant dimension reduction in the computation of the superhedging cost, via an alternate lower-dimensional formulation of the primal problem as a convex optimization problem, which is qualitatively similar to the duality which was formally sketched using linear programming arguments in Duembgen and Rogers [10] for the case where we only consider continuous sample paths. The proof exploits a simplification of a classical result by Rogers (1993) which characterizes the attainable joint laws for the supremum and the drawdown of a uniformly integrable martingale (not necessarily continuous), combined with classical convex duality results from Rockefellar (1974) using paired spaces with compatible locally convex topologies and the Hahn–Banach theorem. We later adapt this result to include additional tradeable One-Touch options using the Kertz and Rösler (1990) condition. We also compute the superhedging cost when in the more realistic situation where there is only finite tradeable European options; for this case we obtain the full duality in the sense of quantile hedging as in Soner (2015), where the superhedge works with probability $1?\epsilon$ where $\epsilon$ can be arbitrarily small), and we obtain an upper bound for the true pathwise superhedging cost. In Section 5, we extend our analysis to include time-dependent barrier options using martingale coupling arguments, where we now have tradeable European options at both maturities at all strikes and tradeable forward starting options at all strikes. This set up is designed to approximate the more realistic situation where we have a finite number of tradeable Europeans at both maturities plus a finite number of tradeable forward starting options.¹     

A priori estimates for rough PDEs with application to rough conservation laws

We introduce a general weak formulation for PDEs driven by rough paths, as well as a new strategy to prove well-posedness. Our procedure is based on a combination of fundamental a priori estimates with (rough) Gronwall-type arguments. In particular this approach does not rely on any sort of transformation formula (flow transformation, Feynman–Kac representation formula etc.) and is therefore rather flexible. As an application, we study conservation laws driven by rough paths establishing well–posedness for the corresponding kinetic formulation.     

RELAXATION SOLVERS FOR IDEAL MHD EQUATIONS-A REVIEW

We have developed approximate Riemann solvers for ideal MHD equations based on a relaxation approach in [4], [5]. These lead to entropy consistent solutions with good properties like guaranteed positive density. We describe the extension to higher order and multiple space dimensions. Finally we show our implementation of all this into the astrophysics code FLASH.     

NONLINEAR EVOLUTION SYSTEMS AND GREEN’S FUNCTION

In this paper, we will introduce how to apply Green’s function method to get the pointwise estimates for the solutions of Cauchy problem of nonlinear evolution equations with dissipative structure. First of all, we introduce the pointwise estimates of the time-asymptotic shape of the solutions of the isentropic Navier-Stokes equations and show to exhibit the generalized Huygen’s principle. Then, for other nonlinear dissipative evolution equations, we will only introduce the result and give some brief explanations. Our approach is based on the detailed analysis of the Green’s function of the linearized system and micro-local analysis, such as frequency decomposition and so on.     

一类二维单个守恒律方程的Riemann问题

本文研究一类二维单个守恒律方程的Riemann问题.用广义特征分析方法研究了这类方程,给出了基本波的分类,解决了初值为两片常数的二维Riemann问题,给出了Riemann解.