A Dirichlet Process Characterization of RBM in a Wedge. Josh Reed. Joint work with Peter Lakner and Bert Zwart. Queueing and Networks Workshop

Similar documents
Obliquely Reflected Brownian motions (ORBMs) in Non-Smooth Domains

Reflected diffusions defined via the extended Skorokhod map

The Skorokhod problem in a time-dependent interval

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 22 12/09/2013. Skorokhod Mapping Theorem. Reflected Brownian Motion

The Generalized Drift Skorokhod Problem in One Dimension: Its solution and application to the GI/GI/1 + GI and. M/M/N/N transient distributions

STEADY-STATE ANALYSIS OF REFLECTED BROWNIAN MOTIONS: CHARACTERIZATION, NUMERICAL METHODS AND QUEUEING APPLICATIONS

Verona Course April Lecture 1. Review of probability

Nonnegativity of solutions to the basic adjoint relationship for some diffusion processes

Deterministic and Stochastic Differential Inclusions with Multiple Surfaces of Discontinuity

Richard F. Bass Krzysztof Burdzy University of Washington

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues

On the p-laplacian and p-fluids

Imaginary Geometry and the Gaussian Free Field

THE SKOROKHOD OBLIQUE REFLECTION PROBLEM IN A CONVEX POLYHEDRON

Optimal investment strategies for an index-linked insurance payment process with stochastic intensity

OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS

Long Time Stability and Control Problems for Stochastic Networks in Heavy Traffic

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

MFE6516 Stochastic Calculus for Finance

On continuous time contract theory

MODELING WEBCHAT SERVICE CENTER WITH MANY LPS SERVERS

The G/GI/N Queue in the Halfin-Whitt Regime I: Infinite Server Queue System Equations

Applications of Ito s Formula

The Cameron-Martin-Girsanov (CMG) Theorem

A Diffusion Approximation for Stationary Distribution of Many-Server Queueing System In Halfin-Whitt Regime

The Skorokhod reflection problem for functions with discontinuities (contractive case)

Lecture 17 Brownian motion as a Markov process

Nonlinear Dynamical Systems Eighth Class

Load Balancing in Distributed Service System: A Survey

Estimation of Dynamic Regression Models

ANALYSIS QUALIFYING EXAM FALL 2016: SOLUTIONS. = lim. F n

Point Process Control

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p

n E(X t T n = lim X s Tn = X s

Lecture 22 Girsanov s Theorem

On Reflecting Brownian Motion with Drift

Translation Invariant Experiments with Independent Increments

BV functions in a Gelfand triple and the stochastic reflection problem on a convex set

Nonlinear representation, backward SDEs, and application to the Principal-Agent problem

OPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS. w 1. v 2. v 3. Ruth J. Williams University of California, San Diego

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.

Robust Mean-Variance Hedging via G-Expectation

HOPF S DECOMPOSITION AND RECURRENT SEMIGROUPS. Josef Teichmann

Stochastic Differential Equations.

and A L T O S O L O LOWELL, MICHIGAN, THURSDAY, OCTCBER Mrs. Thomas' Young Men Good Bye 66 Long Illness Have Sport in

Reflected Brownian Motion

Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2

Particle models for Wasserstein type diffusion

Explicit Solutions for Variational Problems in the Quadrant

P (A G) dp G P (A G)

Optimal stopping for non-linear expectations Part I

Directional Derivatives of Oblique Reflection Maps

Exercises. T 2T. e ita φ(t)dt.

Stability of optimization problems with stochastic dominance constraints

Asymptotic behaviour of solutions of third order nonlinear differential equations

Information and Credit Risk

Hölder continuity of the solution to the 3-dimensional stochastic wave equation

A Concise Course on Stochastic Partial Differential Equations

Solving Linear Rational Expectations Models

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals

1. Stochastic Processes and filtrations

Stochastic integration. P.J.C. Spreij

US01CMTH02 UNIT-3. exists, then it is called the partial derivative of f with respect to y at (a, b) and is denoted by f. f(a, b + b) f(a, b) lim

Assignment #09 - Solution Manual

Brownian Motion on Manifold

CMS winter meeting 2008, Ottawa. The heat kernel on connected sums

Krzysztof Burdzy University of Washington. = X(Y (t)), t 0}

The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge

An Operator Theoretical Approach to Nonlocal Differential Equations

On the stochastic nonlinear Schrödinger equation

Diffuison processes on CR-manifolds

Exercises: Brunn, Minkowski and convex pie

Differentiable Functions

Convexity in R n. The following lemma will be needed in a while. Lemma 1 Let x E, u R n. If τ I(x, u), τ 0, define. f(x + τu) f(x). τ.

Proofs of the martingale FCLT

DIFFERENT KINDS OF ESTIMATORS OF THE MEAN DENSITY OF RANDOM CLOSED SETS: THEORETICAL RESULTS AND NUMERICAL EXPERIMENTS.

Stochastic calculus without probability: Pathwise integration and functional calculus for functionals of paths with arbitrary Hölder regularity

On Stopping Times and Impulse Control with Constraint

The Calculus of Vec- tors

1 Brownian Local Time

Theoretical Tutorial Session 2

Discretization of Stochastic Differential Systems With Singular Coefficients Part II

Stochastic domination in space-time for the supercritical contact process

Random measures, intersections, and applications

Exercises from other sources REAL NUMBERS 2,...,

A class of globally solvable systems of BSDE and applications

Positive recurrence of reflecting Brownian motion in three dimensions

SDE Coefficients. March 4, 2008

A set C R n is convex if and only if δ C is convex if and only if. f : R n R is strictly convex if and only if f is convex and the inequality (1.

Chapter 17: Double and Triple Integrals

LECTURE 2: LOCAL TIME FOR BROWNIAN MOTION

Fundamental Inequalities, Convergence and the Optional Stopping Theorem for Continuous-Time Martingales

ELEMENTS OF STOCHASTIC CALCULUS VIA REGULARISATION. A la mémoire de Paul-André Meyer

Asymptotic inference for a nonstationary double ar(1) model

Lecture 10: Semi-Markov Type Processes

On countably skewed Brownian motion

Control and Observation for Stochastic Partial Differential Equations

Risk Measures in non-dominated Models

Proving the Regularity of the Minimal Probability of Ruin via a Game of Stopping and Control

Transcription:

A Dirichlet Process Characterization of RBM in a Wedge Josh Reed Joint work with Peter Lakner and Bert Zwart Queueing and Networks Workshop Institute for Mathematics and its Applications May 15th, 2018

RBM in a Wedge Let S be a wedge in R 2 whose state space in polar coordinates is given by for some 0 < ξ < 2π. S = {(r,θ) : r 0,0 θ ξ} ξ

RBM in a Wedge Consider a standard 2-d Brownian motion X started from a point z S and reflected back into the interior of the wedge S whenever it hits the boundary of S. The method of reflection is as follows. Let S 1 = {(r,θ) : r > 0,θ = 0} denote the lower boundary of the wedge, and let denote its upper boundary. S 2 = {(r,θ) : r > 0,θ = ξ}

RBM in a Wedge The angle of reflection off of the edge S 1 is denoted by π/2 < θ 1 < π/2, and the angle of reflection off of the edge S 2 = {(r,θ) : r > 0,θ = ξ} is denoted by π/2 < θ 2 < π/2. Bothoftheseanglesaremeasuredwithrespecttotheinward facing normal off their respective edge, with angles directed toward the origin assumed to be positive. We denote by n 1 the inward facing unit normal off of the edge S 1, and by v 1 the unit vector of reflection off of S 1. Similar quantities are defined for S 2.

RBM in a Wedge S 2 n 2 v 2 θ 2 ξ θ 1 n 1 v 1 S 1

RBM in a Wedge X Z In many cases, it turns out to be a difficult problem to obtain samplepathwise such a reflected process Z when starting from a Brownian motion X. However, a process with the probablistic characteristics of the reflected process Z was rigorously defined by Varadhan and Williams (1985).

Definition Let C S denote the space of continuous functions with domain R + = [0, ) and range S. For each t 0 and ω C S, denote by Z(t) : C S S the coordinate map Z(t)(ω) = Z(t,ω) = ω(t), and also define the coordinate mapping process Z = {Z(t),t 0}. Then, for each t 0, set M t = σ{z(s),0 s t}, and set M = σ{z(s),s 0}. Next, for each n 1 and F R 2, denote by C n (F) the set of n- times continuously differentiable functions in some domain containing F, and let C n b be the set of functions in C n (F) that have bounded partial derivatives up to and including order n on F. Finally, define the differential operators D j = v j for j = 1,2, and denote by the Laplacian operator.

Definition Definition 1.[Varadhan and Williams (1985)] A family of probability measures {P z,z S} on (C S,M) is said to solve the submartingale problem if for each z S, the following 3 conditions hold, 1. P z (Z(0) = z) = 1, 2. For each f Cb 2 (S), the process { f(z(t)) 1 t } f(z(s))ds, t 0 2 0 is a submartingale on (C S,M,M t,p z ) whenever f is constant in a neighborhood of the origin and satisfies D i f 0 on S i for i = 1,2, 3. E z [ 0 ] 1{Z(t) = 0}dt = 0.

Definition Reflected Brownian motion arises as the heavy traffic limit of many queueing networks. For instance, the multi-dimensional queue length process of open queueing networks with homogeneous customer populations may be approximated by RBM in a wedge. See Reiman (1984) and Harrison and Williams (1987). The generalized processor sharing queue may also be approximated by RBM in a wedge. See Ramanan and Reiman (2003).

The Quantity α Let α = θ 1+θ 2. ξ α < 1 α = 1 α > 1 Much of the behavior of Z is determined by the value of α.

The Quantity α Varadhan and Williams (1985) proved that if α < 2, then there exists a unique solution to the submartingale probelm. They also showed that if α 2, then there is a solution satisfying only Conditions 1 and 2 of the submartingale problem. In this case, the process Z is absorbed at the origin upon reaching it. The results of Williams (1985) state that Z is a semi-martingale if and only if α < 1 or α 2. The results of Kang and Ramanan (2010) imply that Z is a Dirichlet process if α = 1.

Dirichlet Process Definition Definition 2.Let z S. A continuousprocess Y defined on (C S,F,F t,p z ) is said to be of zero energy if for each T > 0, t i π n Y(t i ) Y(t i 1 ) 2 P 0 as n, (1) for any sequence {π n,n 1} of partitions of [0,T] with π n 0 as n. The notion of a zero-energy process can be used to define a Dirichlet process.

Dirichlet Process Definition Definition 3. Let z S. The stochastic process Z is said to be a Dirichlet process on (C S,F,F t,p z ) if we may write Z = X +Y, (2) where X is an F t -adapted local martingale and Y is a continuous, F t - adapted zero energy process with Y(0) = 0. The class of Dirichlet processes is larger than the class of semi-martingales. Moreover, there exists a change-of-variable formula for the class of Dirichlet processes. The decomposition Z = X+Y appearing in Definition 3 may be shown to be unique.

Dirichlet Processes Result Theorem 1.Suppose that 1 < α < 2. Then, Z has the decomposition Z = X +Y, where (X,Y) is a pair of processes on (C S,F,F t ) such that for each z S, X is a standard Brownian motion started from z and Y is a process of zero energy on (C S,F,F t,p z ). In particular, for each z S, the process Z is a Dirichlet process on (C S,F,F t,p z ). The pair of processes (X,Y) appearing in Theorem 1 do not depend on z S. That is, there is a single pair of processes (X,Y) on (C S,F,F t ) such that the statement of the above theorem hold for each z S.

p-variation Definition Definition 4. Let f : R + R d and p > 0. Then, f is said to be of finite strong p-variation if for each T 0, V p (f,t) = sup f(t i ) f(t i 1 ) p : π π(t) <, t i π where the supremum is taken over all partitions π of [0,T]. Strong p-variation is a generalization of the concept of bounded variation. If V p (f,t) <, then V q (f,t) < for q p > 0. Moreover, functions of finite strong p-variation are, up to a time-change, locally Hölder continuous with exponent 1/p.

p-variation Result Theorem 2. Suppose that 1 < α < 2. Then, for each p > α and z S, P z (V p (Y,[0,T]) < + ) = 1, T 0. Furthermore, for each 0 < p α, P 0 (V p (Y,[0,T]) < + ) = 0, T 0. The sample paths of Y become rougher as α increases from 1 to 2.

The Extended Skorokhod Problem For our last result, we provide a rigorous statement of the fact that Z is the reflected version of the Brownian motion X on all of R +. Let D(R +,R d ) denote the space of R d -valued functions, with domain R +, that are right-continuous with left limits. Also, let D S (R +,R 2 ) be the set of f D(R +,R 2 ) such that f(0) S. Next, let d( ) be a set-valued mapping defined on S such that d(z) is a closed convex cone in R 2 for every z S. In particular, we define {αv 1,α 0}, for z S 1, {αv 2,α 0}, for z S 2, d(z) = (3) V, for z = 0, {0}, for z Int(S).

The Extended Skorokhod Problem Definition 5.[Ramanan(2006)] Thepair of processes (φ,η) D S (R +,R 2 ) D(R +,R 2 ) solve the ESP (S,d( )) for ψ D(R +,R 2 ) if φ(0) = ψ(0), and if for all t R +, the following properties hold, 1. φ(t) = ψ(t)+η(t), 2. φ(t) S, 3. For every s [0,t], 4. η(t) η(t ) Conv[d(φ(t))]. η(t) η(s) Conv[ u (s,t] d(φ(u))], Note that the pushing function η is not required to be of bounded variation.

The Extended Skorokhod Problem Theorem 3.Suppose that 1 < α < 2. Then, for each z S, the ESP (S,d( )) for the Brownian motion X on (C S,F,F t,p z ) has a solution P z -a.s. if and only if Conv(V {αv 1,α 0} {αv 2,α 0}) = R 2. (4) In this case, (Z,Y) solves the ESP (S,d( )) for X. By setting V = R 2, one may find a V such that (4) holds. However, using the fact that 1 < α < 2, it may be verified that the smaller set V = {αv 0,α 0} for any v 0 in the interior of S satisfies (4) as well. Note also that we do not claim in Theorem 3 that (Z,Y) is the unique solution to the ESP (S,d( )) for X.

Thank You THANK YOU!