Igor Cialenco. Department of Applied Mathematics Illinois Institute of Technology, USA joint with N.
|
|
- Amy Townsend
- 6 years ago
- Views:
Transcription
1 Parameter Estimation for Stochastic Navier-Stokes Equations Igor Cialenco Department of Applied Mathematics Illinois Institute of Technology, USA joint with N. Glatt-Holtz (IU) Asymptotical Statistics of Stochastic Processes VIII Universite du Maine, Le Mans, March, 211 Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 1
2 Chicago, Summer 21 c Ig.C.
3 The Problem SPDE du(t,x) = θau(t,x)dt+fu(t,x)dt+σm(u(t,x))dw(t,x) u in some suitable Hilbert spaces and stochastic basis (Ω,F,F,P) A - linear operator F, M - some (could be nonlinear) operators W(t, x) - cylindrical Brownian motion u is known for all x G, t [,T] - continuous observations θ, σ - parameters (scalars), unknown Goal: Find estimators ˆθ(ω), ˆσ(ω), ω Ω, for parameters θ, σ by observing a single outcome u = u(ω) over a finite time horizon [,T]. Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 3
4 Stochastic ODE The Drift Girsanov Th, maximize Log-Likelihood Ratio w.r.t. θ du = θu(t)dt+σu(t)dw(t), θt = 1 t t du(s) u(s) = 1 u(t) log t u() σ2 2 Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 4
5 dy = y(t)dt+σdw(t), t T; σ > y t = σ 2 t σ = Stochastic ODE Volatility y T T ; In practice: σ2 N k=1 ( ) 2 y( kt N ) y((k 1)T N ) T the drift θ - approximated, the volatility σ - exactly Regular model 1) dp θ dq exists; 2) has a special form (LAN) Same procedure for all Find MLE by maximizing likelihood ratio WHY? Singular model otherwise Individual approach In particular, if P θ1 P θ2 for θ 1 θ 2, then one may find θ exactly Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 5
6 The Heat Equation (simulated by Euler) du = ν u xx dt+ ε σ γ dw, T=1, ν=.1, γ=, ε= What do we have for SPDEs? Mostly Singular Models We will try to understand the singularity and explore it to find the exact θ.
7 SPDE drift additive noise: Huebner-Khasminskii-Rozovskii 92, 95 Bayesian: Bishwal ( 2) Several parameters: Huebner ( 97) Discrete-time observations: Piterbarg-Rozovskii ( 97) q = 2(m 1 2m) d 1 θ(t)-random: Lototsky ( 4) Small noise: Huebner ( 97), Ibragimov-Khasminskii ( 98, 99) almost commutative case or almost diagonalizable model- Rozovskii-Lototsky ( 97), Lototsky ( 1) additive fractional noise: IgC, Lototsky, Pospisil ( 9) multiplicative noise: IgC and Lototsky ( 8), IgC ( 1) nonlinear SPDE: IgC and Glatt-Holtz ( 11) Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 7
8 SPDE General Framework du(t) = θau(t)dt+f(u)dt+σdw(t), U() = U ( A) a linear, selfadjoint, positive-defined (think Laplace) in H with eigenvalues {Φ k } k 1 CONS in H F maybe nonlinear σdw(t) = k 1 σ kφ k dw k (t), W k,k N ind. Brownian Motions σ known θ - parameter/scalar of interest U(t) H is observed/measured/known for all t [,T] but one ω Ω; Observe one path during t [,T] Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 8
9 SPDE General Framework Formal Procedure to Derive an Estimator Project the full system down to N dimensions P N (H) = H N R N du N = (θau N +Ψ N )dt+p N σdw, U() = U Formally treat Ψ N = P N F(U) as an external and known quantity (independent of θ) Assume that P N σ is invertible on H N and take G = (P N σ) 1 Consider P N,T ν ( ) = P(U N ), the measure on C([,T];R N ) generated by U N. For a reference values θ, apply (formally) Girsanov Theorem and get the likelihood ratio (Radon-Nikodym derivative) dp N,T θ /dp N,T θ Maximize the Log-Likelihood Ratio ˆθ(ω) := argmax θ dp N,T θ /dp N,T θ (ω) Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 9
10 SPDE General Framework dp N,T θ dp N,T θ [ =exp 1 2 (θ θ ) AU N,G 2 du N (t) (θ 2 θ 2 ) AU N,G 2 AU N dt ] (θ θ ) AU N,G 2 ψ N dt, θ N = AUN,G 2 du N + AUN,G 2 P N F(U) dt GAUN 2 dt Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 1
11 SPDE General Framework Motivated by MLE type estimator θ N = ˇθ N = AUN,G 2 du N + AUN,G 2 P N F(U) dt GAUN 2 dt AUN,G 2 du N + AUN,G 2 F N (U N ) dt GAUN 2 dt θ N = AUN,G 2 du N GAUN 2 dt Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 11
12 Successfully applied to: Stochastic Navier Stokes Equations, 2D, additive noise Stochastic Reaction-Diffusion Equation, additive noise Stochastic Fractional Burgers Equation, additive noise Work in progress to derive general conditions on A,F,G that guarantee consistency and asymptotic normality of ˇθ N and θ N as number of modes N ˇθ N = θ + + AUN,G 2 N j=1 σ j(u)φ j dw j (t) GAUN 2 dt AUN,G 2 (F N (U) F N (U N )) dt GAUN 2 dt θ N = θ + + AUN,G 2 N j=1 σ j(u)φ j dw j (t) GAUN 2 dt AUN,G 2 F N (U) dt GAUN 2 dt
13 Stochastic Navier Stokes, 2D du +((U )U ν U)dt = k λ γ k Φ kdw k U =, U() = U Models a viscous 2D incompressible fluid. U velocity field, ν kinematic viscosity Periodic or Dirichlet boundary conditions Φ k, λ k λ 1 k eigenfunctions, eigenvalues of the Stokes operator A (P N σ) 1 = A γ on H N = Span{Φ 1,...,Φ N } When γ > 1,! Strong, Pathwise solutions i.e. U C([, ),H 1 ) L 2 loc ([, )H2 ). Higher regularity for larger γ. Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 13
14 Stochastic Navier Stokes, 2D Put G = A α γ, for some α R (determined later). Then, ˇν N = = A1+2α U N,dU N +P N B(U N )dt A1+α U N 2 dt N k=1 λ 1+2α [ k u k(t)du k (t)+ u k(t)b k (t)dt ] N k=1 λ 2(1+α) k u k(t) 2 dt ˆν N = A1+2α U N,dU N = A1+α U N 2 dt = N k=1 λ1+2α k N k=1 λ1+2α k N k=1 λ2(1+α) k [u 2 k (T) u2 k () Tλ 2γ T u k(t) 2 dt 2 N k=1 λ2(1+α) k k ] u k(t)du k (t) u k(t) 2 dt Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 14
15 Stochastic Navier Stokes, 2D Main Result Theorem (Ig.C. & N. Glatt-Holtz 211) Let U = U(ω,t,x) is a solution of the 2D SNSE with Periodic or Dirichlet BCs. If γ > 1 (also γ 1 < 1/4 in the Dirichlet case) then: (i) Consistency: If α > γ 1, then N ˇν lim N = ν, N ˆν lim N ν in Probability. (ii) Asymptotic normality (rate N): If α > γ 1 2, then N( ν N ν) ( d 2ν(α γ +1) 2 ) η N, λ 1 T(α γ +1/2) Note: In particular we can take α = γ. This corresponds to formal MLE. Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 15
16 Sketch of the proof A Splitting Argument ˆν N ν = A1+2α γ(ūn +R N ), k Φ kdw k +A γ P N B(U)dt A1+α (ŪN +R N ) 2 dt = E A1+α Ū N 2 dt A1+α (ŪN +R N ) 2 dt A1+2α γ (ŪN +R N ), k Φ kdw k +A γ P N B(U)dt E A1+αŪN 2 dt Decompose U = Ū +R = linear + nonlinear dū +νaūdt = σdw, Ū() =. t R+νAR = B(U), R() = U. Find explicit and exact rates for moments of the linear part R is more regular in comparison to Ū Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 16
17 Sketch of the proof A Splitting Argument N δ 1 A1+α (ŪN +R N ) 2 dt E 1 T A1+αŪN 2 dt A1+2α γ U N, N k=1 Φ kdw k A1+αŪN 2 dt N δ 2 A1+2α γ R N, N k=1 Φ kdw k A1+αŪN 2 dt N δ 3 A1+2α U N,P N (B(U)) A1+αŪN 2 dt N A1+2α γ U N, N k=1 Φ kdw k E 2γ(α γ +1) 2 T A1+α U N N(, 2 dt λ 1 T(α γ +1/2) ) Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 17
18 1 The Heat Equation (simulated by Euler) du = ν u xx dt+ ε σ γ dw, T=1, ν=.2, γ=.1, ε= du = νu xx dt+dw(t,x), u(,x) = u, u(t,) = u(t,1) =
19 The Burgers Equation (simulated by Euler) du = ν u xx dt β u u x dt + ε σ γ dw, T=1, nu=.2, γ=.1, ε= x t du = νu xx dt uu x dt+dw(t,x), u(,x) = u, u(t,) = u(t,1) =
20 The Burgers Equation (simulated by Euler) du = ν u xx dt β u u x dt + ε σ γ dw, T=1, nu=.2, γ=.1, ε= du = νu xx dt uu x dt+dw(t,x), u(,x) = u, u(t,) = u(t,1) =
21 .23 Heat Equation du = ν u xx dt + σ γ dw, ν=.2, T=1, α=.1, γ= nucheck, Method 1, u(s)du(s) nucheck, Method 2,.5(u 2 (T) u 2 () T λ 2γ nucheckexp, Method 1, u(s)du(s) nucheckexp, Method 2,.5(u 2 (T) u 2 () T λ 2γ True Parameter du = νu xx dt+dw(t,x), u(,x) = u, u(t,) = u(t,1) =
22 Burgers Equation du = ν u xx dt β u u x dt + σ γ dw, ν=.2, T=1, α=.1, γ=.1 nucheck, Method 1, u(s)du(s) nucheck, Method 2,.5(u 2 (T) u 2 () T λ 2γ nucheckexp, Method 1, u(s)du(s) nucheckexp, Method 2,.5(u 2 (T) u 2 () T λ 2γ True Parameter du = νu xx dt uu x dt+dw(t,x), u(,x) = u, u(t,) = u(t,1) =
23 x nuhat for α=γ=.1 nuhat for α = du = νu xx dt uu x dt+dw(t,x), u(,x) = u, u(t,) = u(t,1) =
24 Future Work
25 Future Work... a lot
26 Thank You! The end of the talk But not of the story... Ig.Cialenco, Applied Math, IIT Le Mans, March 22 Slide # 25
Hypothesis testing for Stochastic PDEs. Igor Cialenco
Hypothesis testing for Stochastic PDEs Igor Cialenco Department of Applied Mathematics Illinois Institute of Technology igor@math.iit.edu Joint work with Liaosha Xu Research partially funded by NSF grants
More informationParameter Estimation for the Stochastically Perturbed Navier-Stokes Equations
Parameter Estimation for the Stochastically Perturbed Navier-Stoes Equations Igor Cialenco and Nathan Glatt-Holtz January 6, 211 Abstract We consider a parameter estimation problem to determine the viscosity
More informationParameter Estimation for the Stochastically Perturbed Navier-Stokes Equations
Parameter Estimation for the Stochastically Perturbed Navier-Stoes Equations Igor Cialenco and Nathan Glatt-Holtz June 1, 21 Abstract We consider a parameter estimation problem to determine the viscosity
More informationAsymptotic Properties of an Approximate Maximum Likelihood Estimator for Stochastic PDEs
Asymptotic Properties of an Approximate Maximum Likelihood Estimator for Stochastic PDEs M. Huebner S. Lototsky B.L. Rozovskii In: Yu. M. Kabanov, B. L. Rozovskii, and A. N. Shiryaev editors, Statistics
More informationParameter Estimation for Stochastic Evolution Equations with Non-commuting Operators
Parameter Estimation for Stochastic Evolution Equations with Non-commuting Operators Sergey V. Lototsky and Boris L. Rosovskii In: V. Korolyuk, N. Portenko, and H. Syta (editors), Skorokhod s Ideas in
More informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
More informationPARAMETER ESTIMATION FOR SPDES WITH MULTIPLICATIVE FRACTIONAL NOISE
PARAMETER ESTIMATION FOR SPDES WITH MULTIPLICATIVE FRACTIONAL NOISE IGOR CIALENCO Abstract. We study parameter estimation problem for diagonalizable stochastic partial differential equations driven by
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) LMS
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationBorel Summability in PDE initial value problems
Borel Summability in PDE initial value problems Saleh Tanveer (Ohio State University) Collaborator Ovidiu Costin & Guo Luo Research supported in part by Institute for Math Sciences (IC), EPSRC & NSF. Main
More informationAsymptotic analysis of the sieve estimator for a class of parabolic SPDEs
Asymptotic analysis of the sieve estimator for a class of parabolic SPDEs M. Huebner Michigan State University S. Lototsky Massachusetts Institute of Technology Published in Scandinavian Journal of Statistics,
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationRepresenting Gaussian Processes with Martingales
Representing Gaussian Processes with Martingales with Application to MLE of Langevin Equation Tommi Sottinen University of Vaasa Based on ongoing joint work with Lauri Viitasaari, University of Saarland
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationA note on error estimation for hypothesis testing problems for some linear SPDEs
Stoch PDE: Anal Comp 214 2:48 431 DOI 1.17/s472-14-36-4 A note on error estimation for hypothesis testing problems for some linear SPDEs Igor Cialenco Liaosha Xu Received: 17 December 213 / Published online:
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More informationA STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
A STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS PETER CONSTANTIN AND GAUTAM IYER Abstract. In this paper we derive a probabilistic representation of the
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationRecent Advances in the Numerical Approximation of Stochastic Partial Differential Equations
NSF/CBMS Regional Conference in the Mathematical Sciences Department of Applied Mathematics Illinois Institute of Technology August 9-13 2010 Chicago USA Recent Advances in the Numerical Approximation
More informationBayesian Estimations for Diagonalizable Bilinear SPDEs
ayesian Estimations for Diagonalizable ilinear SPDEs Ziteng Cheng zcheng7@haw.iit.edu Igor Cialenco cialenco@iit.edu http://math.iit.edu/~igor uoting Gong rgong2@iit.edu http://mypages.iit.edu/~rgong2
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationCandidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2011 2012 FLUID DYNAMICS MTH-3D41 Time allowed: 3 hours Attempt FIVE questions. Candidates must show on each answer book the type
More informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationQuestion 1. The correct answers are: (a) (2) (b) (1) (c) (2) (d) (3) (e) (2) (f) (1) (g) (2) (h) (1)
Question 1 The correct answers are: a 2 b 1 c 2 d 3 e 2 f 1 g 2 h 1 Question 2 a Any probability measure Q equivalent to P on F 2 can be described by Q[{x 1, x 2 }] := q x1 q x1,x 2, 1 where q x1, q x1,x
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationExit times of diffusions with incompressible drifts
Exit times of diffusions with incompressible drifts Andrej Zlatoš University of Chicago Joint work with: Gautam Iyer (Carnegie Mellon University) Alexei Novikov (Pennylvania State University) Lenya Ryzhik
More informationA STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
A STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS PETER CONSTANTIN AND GAUTAM IYER Abstract. In this paper we derive a representation of the deterministic
More informationWiener Chaos Solution of Stochastic Evolution Equations
Wiener Chaos Solution of Stochastic Evolution Equations Sergey V. Lototsky Department of Mathematics University of Southern California August 2003 Joint work with Boris Rozovskii The Wiener Chaos (Ω, F,
More informationCross-diffusion models in Ecology
Cross-diffusion models in Ecology Gonzalo Galiano Dpt. of Mathematics -University of Oviedo (University of Oviedo) Review on cross-diffusion 1 / 42 Outline 1 Introduction: The SKT and BT models The BT
More informationOn some relations between Optimal Transport and Stochastic Geometric Mechanics
Title On some relations between Optimal Transport and Stochastic Geometric Mechanics Banff, December 218 Ana Bela Cruzeiro Dep. Mathematics IST and Grupo de Física-Matemática Univ. Lisboa 1 / 2 Title Based
More informationSome Topics in Stochastic Partial Differential Equations
Some Topics in Stochastic Partial Differential Equations November 26, 2015 L Héritage de Kiyosi Itô en perspective Franco-Japonaise, Ambassade de France au Japon Plan of talk 1 Itô s SPDE 2 TDGL equation
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
More informationA Sufficient Condition for the Kolmogorov 4/5 Law for Stationary Martingale Solutions to the 3D Navier-Stokes Equations
A Sufficient Condition for the Kolmogorov 4/5 Law for Stationary Martingale Solutions to the 3D Navier-Stokes Equations Franziska Weber Joint work with: Jacob Bedrossian, Michele Coti Zelati, Samuel Punshon-Smith
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationRough Burgers-like equations with multiplicative noise
Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationTranslation Invariant Experiments with Independent Increments
Translation Invariant Statistical Experiments with Independent Increments (joint work with Nino Kordzakhia and Alex Novikov Steklov Mathematical Institute St.Petersburg, June 10, 2013 Outline 1 Introduction
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More informationNavier-Stokes equations in thin domains with Navier friction boundary conditions
Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
More informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationDensities for the Navier Stokes equations with noise
Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other
More informationwith Current Status Data
Estimation and Testing with Current Status Data Jon A. Wellner University of Washington Estimation and Testing p. 1/4 joint work with Moulinath Banerjee, University of Michigan Talk at Université Paul
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationNonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University
Nonlinear stability of time-periodic viscous shocks Margaret Beck Brown University Motivation Time-periodic patterns in reaction-diffusion systems: t x Experiment: chemical reaction chlorite-iodite-malonic-acid
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationOn Some Variational Optimization Problems in Classical Fluids and Superfluids
On Some Variational Optimization Problems in Classical Fluids and Superfluids Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas
More informationHomogeneous Stochastic Differential Equations
WDS'1 Proceedings of Contributed Papers, Part I, 195 2, 21. ISBN 978-8-7378-139-2 MATFYZPRESS Homogeneous Stochastic Differential Equations J. Bártek Charles University, Faculty of Mathematics and Physics,
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationMalliavin Calculus in Finance
Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationOn the well-posedness of the Prandtl boundary layer equation
On the well-posedness of the Prandtl boundary layer equation Vlad Vicol Department of Mathematics, The University of Chicago Incompressible Fluids, Turbulence and Mixing In honor of Peter Constantin s
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFractional Laplacian
Fractional Laplacian Grzegorz Karch 5ème Ecole de printemps EDP Non-linéaire Mathématiques et Interactions: modèles non locaux et applications Ecole Supérieure de Technologie d Essaouira, du 27 au 30 Avril
More informationMultiscale Hydrodynamic Phenomena
M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is
More informationExit times of diffusions with incompressible drift
Exit times of diffusions with incompressible drift Gautam Iyer, Carnegie Mellon University gautam@math.cmu.edu Collaborators: Alexei Novikov, Penn. State Lenya Ryzhik, Stanford University Andrej Zlatoš,
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More informationEQUITY MARKET STABILITY
EQUITY MARKET STABILITY Adrian Banner INTECH Investment Technologies LLC, Princeton (Joint work with E. Robert Fernholz, Ioannis Karatzas, Vassilios Papathanakos and Phillip Whitman.) Talk at WCMF6 conference,
More informationContinuous Time Finance
Continuous Time Finance Lisbon 2013 Tomas Björk Stockholm School of Economics Tomas Björk, 2013 Contents Stochastic Calculus (Ch 4-5). Black-Scholes (Ch 6-7. Completeness and hedging (Ch 8-9. The martingale
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationThomas Knispel Leibniz Universität Hannover
Optimal long term investment under model ambiguity Optimal long term investment under model ambiguity homas Knispel Leibniz Universität Hannover knispel@stochastik.uni-hannover.de AnStAp0 Vienna, July
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationSeparation for the stationary Prandtl equation
Separation for the stationary Prandtl equation Anne-Laure Dalibard (UPMC) with Nader Masmoudi (Courant Institute, NYU) February 13th-17th, 217 Dynamics of Small Scales in Fluids ICERM, Brown University
More informationTHE POSTPROCESSING GALERKIN AND NONLINEAR GALERKIN METHODS A TRUNCATION ANALYSIS POINT OF VIEW
SIAM J. NUMER. ANAL. Vol. 41, No. 2, pp. 695 714 c 2003 Society for Industrial and Applied Mathematics THE POSTPROCESSING GALERKIN AND NONLINEAR GALERKIN METHODS A TRUNCATION ANALYSIS POINT OF VIEW LEN
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationA Class of Fractional Stochastic Differential Equations
Vietnam Journal of Mathematics 36:38) 71 79 Vietnam Journal of MATHEMATICS VAST 8 A Class of Fractional Stochastic Differential Equations Nguyen Tien Dung Department of Mathematics, Vietnam National University,
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationImplicit-explicit exponential integrators
Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators
More informationHypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations.
Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint works with Michael Salins 2 and Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationMath 575-Lecture 26. KdV equation. Derivation of KdV
Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
More information