with deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics
|
|
- Richard Bryan
- 5 years ago
- Views:
Transcription
1 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied Mathematics, University of Crete, Heraklion, Greece Institute of Applied and Computational Mathematics, FO.R.T.H., Heraklion, Crete, Greece 1
2 1. The Problem Non-Homogeneous Cahn-Hilliard Equation tu = ( 2 u + f (u)) + F (x, t; ), x Ω R n, n 3, t > 0. Fourth order non-linear p.d.e., u = u(x, t) (mass concentration), > 0 (measure for the inner interfaces width). f (u) = uf(u) = u( 1 4 (u2 1) 2 ) = u(u 2 1). F smooth function on Ω (0, + ). f = 1 4 (u2 1) 2 double equal-well potential of minimum value 0 at u = ±1. 2
3 2. Model for concentration evolution The mass container: Ω bounded domain of R n, n 3. Mass concentration: u(x, t). Total mass in Ω: m(t) = Ω u(x, t)dx. Evolution process Reason: Non-Equilibrium state - Non-minimized Free Energy. (1) Concentration Evolution to equilibrium state lim t t eq u(x, t) = u eq(x). (2) Equilibrium (Free energy minimization) - Evolution stops u(x, t) = ueq(x), for any t teq. 3
4 Critical dynamics: Ginzburg-Landau model Hohenberg, Halperin, (Rev. Mod. Phys., 1977) Define: E := u 2 + f(u), for f(u) = 1 4 (u2 1) 2 double equal-well potential. System with no conservation laws: Model A (Allen-Cahn related) ut = Γ0Eu + Γ0h(x, t) + θ(x, t), h external field, θ Gaussian white noise source. System with conservation: Model B (Cahn-Hilliard related) in Model A replace Γ0 with λ0 and get ut = λ0 (Eu h(x, t)) + θ(x, t). 4
5 The generalized Cahn-Hilliard equation Gurtin, (Physica D, 1996) ut = ( 2 u + f (u) γ(x, t)) + m(x, t), γ external microforce, m external mass supply. Chemical potential: µ := 2 u + f (u) γ(x, t). Time dependent total free energy of concentration (1) Local energy of u: E(x, t; u) = Ek(x, t; u) + Ed(x, t; u), (2) Ek kinetic energy: u(x, t) (speed) concentration change in space. (3) Ed dynamical energy: Ed = F(x, t; u)du Due to a space-time force field F = F(x, t; u). Total Free Energy Functional: FE[u](t) := E(x, t; u)dx. Ω 5
6 Concentration Evolution Equation tu = K L, K const.> 0, L := ue(x, t; u) = uek(x, t; u) + F(x, t; u). Mass evolution: tm(t) = K Ω nl(s, t; u)ds. Mass conservation condition nl(x, t; u) = 0 on Ω n( uek(x, t; u) + F(x, t; u)) = 0 on Ω tm(t) = 0. Evolution Stages (a) Homogenization (Constant concentration) (b) Spinodal (c) Coarsening (Sharp Interface) (d) Equilibrium 6
7 Figure 1: (1) Homogenization Spinodal (2) Spinodal: Deterministic Model VS Stochastic Model 7
8 Figure 2: Coarsening 8
9 Figure 3: Spinodal Coarsening 9
10 3. Homogeneous Cahn-Hilliard equation Fourth order non-linear P.D.E. tu = K L, K := 1, L = uek + ued, with Ek(x, t; u) = u 2 Ed(x, t; u) = f (u)du = f(u) (double equal well potential) (1) Non-Linear force field considered F(x, t; u) = F(u) := f (u) = u(u 2 1) models the tendency of a binary alloy (two species) forced to homogenization to return in a two separated phases equilibrium. (2) The kinetic energy of the mass 2 models the material preference to be as uniform as possible. 10
11 Replace L and get Homogeneous Cahn-Hilliard equation tu = ( 2 u + f (u)), t > 0. Simple phase separation model of a binary alloy at a fixed temperature, > 0 measure for the inner interfaces width (phase transitions). B.C.: nl(x, t; u) = 0 on Ω mass conservation and free energy decreasing condition n( 2 u + f (u)) = 0 on Ω. Existence-Uniqueness of solution Homogeneous Cahn Hilliard with initial condition u0(x) H 2 (Ω), with nu0 = 0 on Ω and boundary conditions: n( 2 u + f (u)) = 0 on Ω, and nu = 0 on Ω, admits a unique solution, Elliott - Zheng, (Arch. Rat. Mech. Anal., 1986). 11
12 4. Non-homogeneous term F in Cahn-Hilliard: F := Fex + Fn Fex external field (source term, production rate in a reaction diffusion process A + B C). Fn Gaussian white noise source (thermal noise). Thermal noise= 0 source brings system to equilibrium. Source = deterministic + noise. Mean-field assumption: source noise = 0. The Equation - General Case Mass conservation condition: tu = ( 2 u + f (u)) + F (x, t; ), t > 0. Ω n( 2 u + f (u))ds = F (x, t; )dx. Ω 12
13 Case 1: Extra dynamics in free energy F := F2 = ψ(x, t; ) may model a free energy extra term (Gurtin-in chemical potential, Langer-ψ = ψ(x) (Annals of Physics, 1971)) tu = ( 2 u + [f (u) + ψ(x, t; )]) = ( 2 u + f (u)) + F2(x, t; ), t > 0. Ed(x, t; u) = f(u) + ψ(x, t; )u ued(x, t; u) = f (u) + ψ(x, t; ). Extra linear potential ψ(x, t; )u, extra space-time force field ψ (reaction - boundary interaction rate). Mass conservation b.c: n( 2 u + [f (u) + ψ(x, t; )]) = 0 on Ω. Case 2: Extra dynamics in free energy, free energy independent terms tu = ( 2 u + [f (u) + ψ(x, t; )]) + F1(x, t; ) = ( 2 u + f (u)) + F1(x, t; ) + F2(x, t; ), t > 0, F := F1 + F2. 13
14 5. Asymptotics 1. Multidimensional case, O() limit problem Alikakos-Bates-Chen (Arch. Rat. Mech. Anal., 1994), Pego (Proc. R. Soc. Lond., 1989) Homogeneous Cahn-Hilliard on R n, n 2 Hele-Shaw Inner Free Boundary Problem. Consider ut(x, t; ) = ( u + f (u) ), x Ω, t > 0, u(x, 0) = u0(x), nu = n( u + f (u) ) = 0 on Ω. 14
15 Set v = lim ). 0 + ( u f (u) Then v satisfies the sharp interface Hele-Shaw Problem: v = 0 in Ω\Γ(t), t > 0, nv = 0 on Ω, v = λh on Γ(t), V = 1 2 ( nv + nv ), Γ(0) = Γ0. Γ(t) closed n 1 dimensional hypersurface (sharp interface) H mean curvature of Γ, and V velocity of Γ Ω = Ω + (t) Ω (t) Γ(t), λ > 0. 15
16 Main result for the Non-Homogeneous case Consider G uniformly bounded as 0 +. tu(x, t; ) = ( u + f (u) ) + G(x, t; ), 1. Split G: G = G1 G2, G1(x, t; ), G2(x, t; ), uniformly bounded as Get: 3. Define: v() := u f (u) 4. Linearize: v() = v0 + v1 + O( 2 ). ut = ( u + f (u) + G2, get the system G2) + G1, ut = v + G1, v = u f (u) + G2. 16
17 Theorem. The order O() approximation of v (i.e. v0), satisfies the non-homogeneous problem v0 = G1 + O( 2 ) in Ω\Γ(t), t > 0, v0 = λh + G2 + O() on Γ(t), V = 1 2 ( Qz( ) Qz(+ )) + O( 2 ), Q(z) := v 1, z := d, d = d(x, t) := inf y Γ(t) x y. (1) Order of convergence to a homogeneous problem: = min{o(), O(G1), O(G2)}. (2) G2 ( Free energy): Local contribution while, G1: Global contribution. Assume that O(G1) = O( λ 1 ), O(G 2) = O( λ 2 ) then 1. If min{λ1, λ2} = 1 order of convergence to homogeneous Hele-Shaw remains the same. 2. If 0 < min{λ1, λ2} < 1 convergence is delayed. 3. If O(G1), O(G2) =const. Non-Homogeneous Cahn-Hilliard Non-Homogeneous Hele-Shaw. 17
18 Conclusion Assumption: Convergence to equilibrium admits that Homogeneous Hele-Shaw problem appears first. Non-Homogeneous terms may slow down equilibrium. Remarks: 1. Rogers, Elder, Desai, (Phys. Rev. B, 1988): Numerical experiments demonstrate that increasing of the noise strength leads to more diffuse interfaces and slower transitions. 2. Physical models: Noise in G. Split G On the limit: noise in G1, noise in G2. Free energy noise terms are local, as they appear at phase transitions Γ(t). Other noises (obtained by G1) influence globally evolution process. 18
19 Mass Conservation and minimal surface conditions Consider: Non-Homogeneous Hele-Shaw Problem v = G1 in Ω\Γ(t), t > 0, nv = 0 on Ω, v = λh + G2 on Γ(t), V = 1 2 ( nv + nv ), Γ(0) = Γ0. Get: d Per(Γ(t)) := dt Γ HV = 1 2λ d dt Vol(Ω (t)) := V = 1 Γ 2 Ω\Γ G1v + 1 λ G1 Ω\Γ Γ G2V 1 2λ Ω\Γ v 2 (1) d dt Per(Γ(t)) 0 minimal surface (2) d dt Vol(Ω (t)) = 0 mass conservation 19
20 7. One dimensional case, O( ) limit problem Ω R, 1 ξ(x, t ξ(x, t ) uniformly bounded as 0+. tu = ( u + f (u) ) 1 + ξ(x, t ), t > 0, ) deterministic free energy independent that follows white noise scale in t, Noise scale: Consider white noise Wt defined as the t derivative of an Itô integral ( L 2 (P )) Wt := d t dt 0 1dB t s, use 0 1dB s = Bt and get by Brownian scale in R Bt = B t = t 1dBs, 0 then set τ = t/ and get Wt = d dt B t = d dt t 0 1dBs = d dτ τ 0 1dBs = 1 Wτ. 20
21 Space and time scalings: x, z := x, and t, τ := t Linearize: u = u0(t, τ, x, z) + u1(t, τ, x, z) + u2(t, τ, x, z) + Outer expansion: u =1 + u 1/2 (τ, x) + v1(t, x)+ u1(τ, x) + 3/2 v 3/2 (t, x) + 3/2 u 3/2 (τ, x) + 2 v2(t, x) + 2 u2(τ, x) + Compute: tu = 1 τ u 1/2 + tv1 + τ u1 + 3/2 tv 3/2 + τ u 3/2 + 2 tv2 + τ u2 + 21
22 and u + f (u) = 1 f (1) + 1 f (1)u 1/2 + f (1)(v 1 + u 1 ) + 1 f (1)u 2 1/2 2 + ( u 1/2 (v 1 + u 1 )f (1)) + f (1)(u 3/2 + v 3/2 ) ) + ( f (1)(v 2 + u 2 ) (v 1 + u 1 ) 2 f (1) + f (1)u 1/2 (v 3/2 + u 3/2 ) ) + ( u 1/2 + f (1)u 1/2 (v 2 + u 2 ) + f (1)(v 1 + u 1 )(v 3/2 + u 3/2 ) ) + 2( v 1 u f (1)(v 3/2 + u 3/2 ) 2 + f (1)(v 1 + u 1 )(v 2 + u 2 ) ) + 5/2( v 3/2 u 3/2 + f (1)(v 3/2 + u 3/2 )(v 2 + u 2 ) ) + 3( v 2 u f (1)(v 2 + u 2 ) 2) + Collect the O( 1 ) terms to obtain τ u 1/2 = f (1) u 1/2 + ξ(x, τ) on Ω\Γ. 22
23 Ω R, expand d (distance from the interface): d(x, t, τ) = x ( x0(t) + y 1/2 (τ) + (x1(t) + y1(τ)) ) + Inner expansion: d(x, t, τ) u = m( ) + d(x, t, τ) m 1/2 ( ) + d(x, t, τ) u 1/2 (τ, x) + m1( ) +v1(t, x) + u1(τ, x) + 3/2 d(x, t, τ) m 3/2 ( ) + 3/2 v 3/2 (t, x) + 3/2 u 3/2 (τ, x) + 2 d(x, t, τ) m2( ) + 2 v2(t, x) + 2 u2(τ, x) +, m unique solution of Euler-Lagrange equation, for z := d m (z) + f ( m(z)) = 0, for z R lim m(z) = ±1, m(0) = 0. z ± 23
24 Expand f (u), use inner expansion, and seek: u + f (u) = 1 ( m + f ( m)) + 1 ( m 1/2 + f ( m)(m 1/2 + u 1/2 ) + ( m 1 + f ( m)(m 1 + v 1 + u 1 ) f ( m)(m 1/2 + u 1/2 ) 2 ) + ( (m 1/2 + u 1/2 )(m 1 + v 1 + u 1 )f ( m))+ ( m 3/2 + f ( m)m 3/2 ) +f ) ( m)(u 3/2 + v 3/2 ) }{{} L mm 3/2 ( + m 2 + f ( m)(m 2 + v 2 + u 2 ) + 1 (m 1 + v 1 + u 1 ) 2 f ( m) 2 +f ) ( m)(m 1/2 + u 1/2 )(m 3/2 + v 3/2 + u 3/2 ) + ( u 1/2 + f ( m)(m 1/2 + u 1/2 )(m 2 + v 2 + u 2 ) +f ) ( m)(m 1 + v 1 + u 1 )(m 3/2 + v 3/2 + u 3/2 ) + 2( v 1 u f ( m)(m 3/2 + v 3/2 + u 3/2 ) 2 2 +f ) ( m)(m 1 + v 1 + u 1 )(m 2 + v 2 + u 2 ) + 5/2( v 3/2 u 3/2 + f ( m)(m 3/2 + v 3/2 + u 3/2 )(m 2 + v 2 + u 2 ) ) + 3( v 2 u f ( m)(m 2 + v 2 + u 2 ) 2) + 24
25 Observe terms of the form: L mm 3/2 = m 3/2 + f ( m)m 3/2, L mm 1/2 = m 1/2 + f ( m)m 1/2 that are linearized operators around m. Apply Fredholm alternative solvability condition to get interface conditions: On the interface all the terms disappear and we end up to principal term Diffusion Equation: τ u 1/2 = f (1) u 1/2 + ξ(x, τ) on Ω\Γ. For ξ = 0, this diffusion is the linearization of the ill-posed problem τ u 1/2 (x, τ) = (f (u 1/2 )(x, τ)). Conclusions 1. Approximate problem Diffusion Equation. 2. No inner boundary conditions lower order approximation ( ) slower time scale. 3. ξ delays the phenomenon, and we capture evolution before coarsening. 25
26 REFERENCES REFERENCES References [1] N. D. Alikakos, P. W. Bates and X. Chen, Convergence of the Cahn-Hilliard Equation to the Hele-Shaw Model, Arch. Rat. Mech. Anal. 128 (1994), pp [2] T. Antal, M. Droz, J. Magnin and Z. Rácz, Formation of Liesengang Patterns: A Spinodal Decomposition Scenarion, Phys. Rev. Lett. 83 (1999), pp [3] J. W. Cahn, J. E. Hilliard, Free Energy of a Nonuniform System. I. Interfacial Free Energy, J. Chem. Phys. 28 (1958), pp [4] C. M. Elliott, S. Zheng, On the Cahn-Hilliard equation, Arch. Rat. Mech. Anal. 96 (1986), pp [5] P. C. Fife, Models for phase separation and their mathematics, El. Journ. Diff. E. 48 (2000), pp [6] M. E. Gurtin, Generalized Ginzburg-Landau and Cahn-Hilliard equations based on a microforce balance, Physica D 92 (1996), pp [7] P. C. Hohenberg and B.I. Halperin, Theory of dynamic critical phenomena, J. Rev. Mod. Phys. 49 (1977), pp
27 REFERENCES REFERENCES [8] J. S. Langer, Theory of Spinodal Decomposition in Alloys, Ann. of Phys. 65 (1971), pp [9] B. csendal, Stochastic Differential Equations, Springer, 2003, New York. [10] R. L. Pego, Front migration in the non-linear Cahn-Hilliard equation, Proc. R. Soc. Lond. A 422 (1989), pp [11] T. M. Rogers, K. R. Elder, and R. C. Desai, Numerical study of the late stages of spinodal decomposition, Phys. Rev. B 37 (1988), pp [12] C. Cardon-Weber, Cahn-Hilliard stochastic equation: existence of the solution and of its density, Bernoulli 7 (2001), pp
Dynamics and stochastic limit for phase transition problems. problems with noise. Dimitra Antonopoulou IACM-FORTH
Dynamics and stochastic limit for phase transition problems with noise IACM-FORTH International Conference on Applied Mathematics, September 16-20, 2013, ACMAC, Heraklion, Greece Table of contents 1 1.1
More informationTHE CAUCHY PROBLEM AND STEADY STATE SOLUTIONS FOR A NONLOCAL CAHN-HILLIARD EQUATION
Electronic Journal of Differential Equations, Vol. 24(24), No. 113, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) THE CAUCHY
More informationStability Analysis of Stationary Solutions for the Cahn Hilliard Equation
Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Peter Howard, Texas A&M University University of Louisville, Oct. 19, 2007 References d = 1: Commun. Math. Phys. 269 (2007) 765
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationTHE SHARP INTERFACE LIMIT FOR THE STOCHASTIC CAHN-HILLIARD EQUATION
THE SHARP INTERFACE LIMIT FOR THE STOCHASTIC CAHN-HILLIARD EQUATION D.C. ANTONOPOULOU $, D. BLÖMKER, AND G.D. KARALI Abstract. We study the -dependent two and three dimensional stochastic Cahn-Hilliard
More informationPREPRINT 2009:20. Finite Element Approximation of the Linearized Cahn-Hilliard-Cook Equation STIG LARSSON ALI MESFORUSH
PREPRINT 9: Finite Element Approximation of the Linearized Cahn-Hilliard-Coo Equation STIG LARSSON ALI MESFORUSH Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationA sharp diffuse interface tracking method for approximating evolving interfaces
A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationStep Bunching in Epitaxial Growth with Elasticity Effects
Step Bunching in Epitaxial Growth with Elasticity Effects Tao Luo Department of Mathematics The Hong Kong University of Science and Technology joint work with Yang Xiang, Aaron Yip 05 Jan 2017 Tao Luo
More informationWell-posedness and asymptotic analysis for a Penrose-Fife type phase-field system
0-0 Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system Buona positura e analisi asintotica per un sistema di phase field di tipo Penrose-Fife Salò, 3-5 Luglio 2003 Riccarda
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationMOTION OF A DROPLET FOR THE STOCHASTIC MASS-CONSERVING ALLEN-CAHN EQUATION
MOTION OF A DROPLET FOR THE STOCHASTIC MASS-CONSERVING ALLEN-CAHN EQUATION D.C. ANTONOPOULOU $, P.W. BATES, D. BLÖMKER, G.D. KARALI Abstract. We study the stochastic mass-conserving Allen-Cahn equation
More informationTHE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,
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 informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
More informationLinear Theory of Evolution to an Unstable State
Chapter 2 Linear Theory of Evolution to an Unstable State c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 2.1 Introduction The simple theory of nucleation that we introduced in
More informationarxiv: v1 [math.na] 15 Aug 2007
A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF THE CAHN-HILLIARD EQUATION AND THE HELE-SHAW FLOW XIAOBING FENG AND HAIJUN WU arxiv:78.116v1 [math.na] 15 Aug 7 Abstract. This paper develops
More informationWeak solutions for the Cahn-Hilliard equation with degenerate mobility
Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationClassical solutions for the quasi-stationary Stefan problem with surface tension
Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationString method for the Cahn-Hilliard dynamics
String method for the Cahn-Hilliard dynamics Tiejun Li School of Mathematical Sciences Peking University tieli@pku.edu.cn Joint work with Wei Zhang and Pingwen Zhang Outline Background Problem Set-up Algorithms
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 informationA fast reaction - slow diffusion limit for propagating redox fronts in mineral rocks
for propagating redox fronts in mineral rocks Centre for Analysis, Scientific computing and Applications (CASA), Department of Mathematics and Computer Science, TU Eindhoven, The Netherlands Joint work
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 informationNumerical Simulations on Two Nonlinear Biharmonic Evolution Equations
Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationHOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS. 1. Introduction
HOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS CARMEN CORTAZAR, MANUEL ELGUETA, ULIO D. ROSSI, AND NOEMI WOLANSKI Abstract. We present a model for
More informationMetastability for interacting particles in double-well potentials and Allen Cahn SPDEs
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ SPA 2017 Contributed Session: Interacting particle systems Metastability for interacting particles in double-well
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationDemixing in viscous fluids: Connection with OT. Felix Otto. joint work with: R. V. Kohn, Y. Brenier, C. Seis, D. Slepcev
Demixing in viscous fluids: Connection with OT Felix Otto Max Planck Institute for Mathematics in the Sciences Leipzig, Germany joint work with: R. V. Kohn, Y. Brenier, C. Seis, D. Slepcev Geometric evolution
More informationSolution of Stochastic Optimal Control Problems and Financial Applications
Journal of Mathematical Extension Vol. 11, No. 4, (2017), 27-44 ISSN: 1735-8299 URL: http://www.ijmex.com Solution of Stochastic Optimal Control Problems and Financial Applications 2 Mat B. Kafash 1 Faculty
More informationTravelling bubbles in a moving boundary problem of Hele-Shaw type
Travelling bubbles in a moving boundary problem of Hele-Shaw type G. Prokert, TU Eindhoven Center for Analysis, Scientific Computing, and Applications (CASA) g.prokert@tue.nl Joint work with M. Günther,
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES XIAOBING FENG AND ANDREAS PROHL Abstract. In this first part
More informationTRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM
TRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM PAUL C. FIFE, Mathematics Department, University of Utah, Salt Lake City, Utah 84112, USA Abstract Existence, uniqueness, and regularity properties
More informationCoupled second order singular perturbations for phase transitions
Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and
More informationMulti-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form
Multi-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form Yipeng Yang * Under the supervision of Dr. Michael Taksar Department of Mathematics University of Missouri-Columbia Oct
More informationTransactions on Modelling and Simulation vol 17, 1997 WIT Press, ISSN X
Thermodynamieally consistent phase-field models of solidification processes C. Charach*, P.C. Fife* "CEEP, J. Blaustein Institute, Ben Gurion University of the Negev, Sede-Boqer Campus, Sede-Boqer, 84990,
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationA global solution curve for a class of free boundary value problems arising in plasma physics
A global solution curve for a class of free boundary value problems arising in plasma physics Philip Korman epartment of Mathematical Sciences University of Cincinnati Cincinnati Ohio 4522-0025 Abstract
More informationAdaptive algorithm for saddle point problem for Phase Field model
Adaptive algorithm for saddle point problem for Phase Field model Jian Zhang Supercomputing Center, CNIC,CAS Collaborators: Qiang Du(PSU), Jingyan Zhang(PSU), Xiaoqiang Wang(FSU), Jiangwei Zhao(SCCAS),
More informationarxiv: v1 [math.ap] 11 Apr 2019
Global well-posedness of the 3D Cahn-Hilliard equations Zhenbang Li Caifeng Liu arxiv:194.6191v1 [math.ap] 11 Apr 219 April 15, 219 Abstract The Cauchy problem of the Cahn-Hilliard equations is studied
More informationComplex Fluids. University of California - Santa Barbara, CA April 26, Abstract
Gravitational Effects on Structure Development in Quenched Complex Fluids V. E. Badalassi, a H. D. Ceniceros b and S. Banerjee c a,c Department of Chemical Engineering, b Department of Mathematics University
More informationACMAC s PrePrint Repository
ACMAC s PrePrint Repository Existence of solution for a generalized stochastic Cahn-Hilliard equation on convex domains Georgia Karali and imitra Antonopoulou Original Citation: Karali, Georgia and Antonopoulou,
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationCoarsening kinetics from a variable-mobility Cahn-Hilliard equation: Application of a semi-implicit Fourier spectral method
PHYSICAL REVIEW E VOLUME 60, NUMBER 4 OCTOBER 1999 Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: Application of a semi-implicit Fourier spectral method Jingzhi Zhu and Long-Qing
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationPartial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula
Partial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula Group 4: Bertan Yilmaz, Richard Oti-Aboagye and Di Liu May, 15 Chapter 1 Proving Dynkin s formula
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationOn some nonlinear parabolic equation involving variable exponents
On some nonlinear parabolic equation involving variable exponents Goro Akagi (Kobe University, Japan) Based on a joint work with Giulio Schimperna (Pavia Univ., Italy) Workshop DIMO-2013 Diffuse Interface
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 informationStochastic Viral Dynamics with Beddington-DeAngelis Functional Response
Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response Junyi Tu, Yuncheng You University of South Florida, USA you@mail.usf.edu IMA Workshop in Memory of George R. Sell June 016 Outline
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationGeneric Singularities of Solutions to some Nonlinear Wave Equations
Generic Singularities of Solutions to some Nonlinear Wave Equations Alberto Bressan Deartment of Mathematics, Penn State University (Oberwolfach, June 2016) Alberto Bressan (Penn State) generic singularities
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationModelling and analysis for contact angle hysteresis on rough surfaces
Modelling and analysis for contact angle hysteresis on rough surfaces Xianmin Xu Institute of Computational Mathematics, Chinese Academy of Sciences Collaborators: Xiaoping Wang(HKUST) Workshop on Modeling
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 informationBSDEs and PDEs with discontinuous coecients Applications to homogenization K. Bahlali, A. Elouain, E. Pardoux. Jena, March 2009
BSDEs and PDEs with discontinuous coecients Applications to homogenization K. Bahlali, A. Elouain, E. Pardoux. Jena, 16-20 March 2009 1 1) L p viscosity solution to 2nd order semilinear parabolic PDEs
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 informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationEXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.
More informationW I A S Uniqueness in nonlinearly coupled PDE systems DICOP 08, Cortona, September 26, 2008
W I A S Weierstrass Institute for Applied Analysis and Stochastics in Forschungsverbund B erlin e.v. Pavel Krejčí Uniqueness in nonlinearly coupled PDE systems joint work with Lucia Panizzi DICOP 8, Cor
More informationStochastic optimal control with rough paths
Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction
More informationFrom nonlocal to local Cahn-Hilliard equation. Stefano Melchionna Helene Ranetbauer Lara Trussardi. Uni Wien (Austria) September 18, 2018
From nonlocal to local Cahn-Hilliard equation Stefano Melchionna Helene Ranetbauer Lara Trussardi Uni Wien (Austria) September 18, 2018 SFB P D ME S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationTHE EXISTENCE OF GLOBAL ATTRACTOR FOR A SIXTH-ORDER PHASE-FIELD EQUATION IN H k SPACE
U.P.B. Sci. Bull., Series A, Vol. 8, Iss. 2, 218 ISSN 1223-727 THE EXISTENCE OF GLOBAL ATTRACTOR FOR A SIXTH-ORDER PHASE-FIELD EQUATION IN H k SPACE Xi Bao 1, Ning Duan 2, Xiaopeng Zhao 3 In this paper,
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
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 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 informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
More informationOn semilinear elliptic equations with measure data
On semilinear elliptic equations with measure data Andrzej Rozkosz (joint work with T. Klimsiak) Nicolaus Copernicus University (Toruń, Poland) Controlled Deterministic and Stochastic Systems Iasi, July
More informationPattern size in Gaussian fields from spinodal decomposition
Pattern size in Gaussian fields from spinodal decomposition Luigi Amedeo Bianchi (joint work with D. Blömker and P. Düren) 8th Random Dynamical Systems Workshop Bielefeld, 6th November 2015 Research partially
More informationA kinetic study of the ordering process in ternary III V semiconductor alloys
J. Phys.: Condens. Matter 9 (1997) 5737 5749. Printed in the UK PII: S0953-8984(97)78738-1 A kinetic study of the ordering process in ternary III V semiconductor alloys Zhi-Feng Huang and Bing-Lin Gu CCAST
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationA C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani
A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes Marco Verani MOX, Department of Mathematics, Politecnico di Milano Joint work with: P. F. Antonietti (MOX Politecnico di
More informationConvergence in Almost Periodic Fisher and Kolmogorov Models
Convergence in Almost Periodic Fisher and Kolmogorov Models Wenxian Shen Department of Mathematics Auburn University Auburn, AL 36849 and Yingfei Yi School of Mathematics Georgia Institute of Technology
More informationSupporting Information. Liesegang rings engineered from charged nanoparticles
Supporting Information Liesegang rings engineered from charged nanoparticles Istvan Lagzi, Bartlomiej Kowalczyk, and Bartosz A. Grzybowski * Department of Chemical and Biological Engineering Northwestern
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationTotal Variation Image Edge Detection
Total Variation Image Edge Detection PETER NDAJAH Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-28, JAPAN ndajah@telecom0.eng.niigata-u.ac.jp
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationWeakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence
Sergei Kuksin Weakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence (Toronto, 10 January, 2014 ) 1 1 Introduction: weak turbulence (WT) (One of) origines: Rudolf Peierls,
More informationLong-time behavior and different shear regimes in quenched binary mixtures
PHYSICAL REVIEW E 75, 011501 2007 Long-time behavior and different shear regimes in quenched binary mixtures G. Gonnella* Dipartimento di Fisica, Università di Bari, and INFN, Sezione di Bari, Via Amendola
More informationA RAPIDLY CONVERGING PHASE FIELD MODEL. Xinfu Chen and Gunduz Caginalp. Christof Eck
DISCETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 15, Number 4, August 2006 pp. 1017 1034 A APIDLY CONVEGING PHASE FIELD MODEL Xinfu Chen and Gunduz Caginalp Department of
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationQualitative behaviour of numerical methods for SDEs and application to homogenization
Qualitative behaviour of numerical methods for SDEs and application to homogenization K. C. Zygalakis Oxford Centre For Collaborative Applied Mathematics, University of Oxford. Center for Nonlinear Analysis,
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationOn Asymptotic Variational Wave Equations
On Asymptotic Variational Wave Equations Alberto Bressan 1, Ping Zhang 2, and Yuxi Zheng 1 1 Department of Mathematics, Penn State University, PA 1682. E-mail: bressan@math.psu.edu; yzheng@math.psu.edu
More informationFor rough surface,wenzel [26] proposed the equation for the effective contact angle θ e in terms of static contact angle θ s
DERIVATION OF WENZEL S AND CASSIE S EQUATIONS FROM A PHASE FIELD MODEL FOR TWO PHASE FLOW ON ROUGH SURFACE XIANMIN XU AND XIAOPING WANG Abstract. In this paper, the equilibrium behavior of an immiscible
More information