Laplacian transport towards irregular interfaces: the mathematics

Size: px
Start display at page:

Download "Laplacian transport towards irregular interfaces: the mathematics"

Transcription

1 Laplacian transport towards irregular interfaces: the mathematics Denis S GREBENKOV Laboratoire de Physique de la Matière Condensée CNRS Ecole Polytechnique, France denisgrebenkov@polytechniqueedu Web: Denis S Grebenkov, CNRS Ecole Polytechnique p1/35

2 Plan Laplacian transport phenomena Theoretical approach Numerical and experimental study Conclusion, questions and perspectives Denis S Grebenkov, CNRS Ecole Polytechnique p2/35

3 Oxygen diffusion towards and through alveolar membranes Denis S Grebenkov, CNRS Ecole Polytechnique p3/35

4 Oxygen diffusion towards and through alveolar membranes Dichotomic tree Volume 4 5 liters Area m 2 Human lung 23 generations : convection 15 diffusion 8 Denis S Grebenkov, CNRS Ecole Polytechnique p3/35

5 Oxygen diffusion towards and through alveolar membranes Laplacian transport C = 0 in the volume C = C 0 at the entrance J = D C Fick s law J n = WC on the membrane Alveoli Alveolated duct Airflow C=C 0 2 C = 0 n C C = 1 Λ Λ = D/W characteristic length of the problem D diffusion coefficient W membrane permeability Denis S Grebenkov, CNRS Ecole Polytechnique p3/35

6 Oxygen diffusion towards and through alveolar membranes Laplacian transport Airflow C=C 0 2 C = 0 C = 0 C = C 0 C n = C Λ in the volume at the entrance on the membrane Alveoli Alveolated duct n C C = 1 Λ Λ = D/W characteristic length of the problem D diffusion coefficient W membrane permeability Denis S Grebenkov, CNRS Ecole Polytechnique p3/35

7 Heterogeneous catalysis Denis S Grebenkov, CNRS Ecole Polytechnique p4/35

8 Heterogeneous catalysis Chemical reaction C = C 0 C n = C Λ A A on the source on the surface Denis S Grebenkov, CNRS Ecole Polytechnique p4/35

9 Heterogeneous catalysis Description C = 0 C = C 0 C n = C Λ Λ = D/K in the volume on the source on the surface characteristic length of the problem D diffusion coefficient K catalyst s reactivity Denis S Grebenkov, CNRS Ecole Polytechnique p4/35

10 Electric transport through an irregular electrode Denis S Grebenkov, CNRS Ecole Polytechnique p5/35

11 Electric transport through an irregular electrode Microroughness (of metallic surface) Irregular geometry on different scales Impedance? 1µm Denis S Grebenkov, CNRS Ecole Polytechnique p5/35

12 Electric transport through an irregular electrode Microroughness (of metallic surface) Irregular geometry on different scales Impedance? 1µm Denis S Grebenkov, CNRS Ecole Polytechnique p5/35

13 Electric transport through an irregular electrode Microroughness (of metallic surface) Irregular geometry on different scales Impedance? 1µm Denis S Grebenkov, CNRS Ecole Polytechnique p5/35

14 Electric transport through an irregular electrode Electric transport: (direct constant potential) V = 0 in the volume V = V 0 on the counter-electrode V n = 1 Λ V on the working electrode 1µm Λ = r/ρ characteristic length of the problem r resistance of the working electrode ρ resistivity of the electrolyte Denis S Grebenkov, CNRS Ecole Polytechnique p5/35

15 Laplacian transport source u = 1 u = 0 Λ u n = u working surface Denis S Grebenkov, CNRS Ecole Polytechnique p6/35

16 Laplacian transport source u = 1 u = 0 Λ u n = u working surface How does a geometrical irregularity influence the transport properties of the system? Denis S Grebenkov, CNRS Ecole Polytechnique p6/35

17 Spectroscopic impedance Mathematical model: mixed boundary value problem for Laplace equation source Diffusive transport C = C 0 C = 0 Λ C n = C working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

18 Spectroscopic impedance Mathematical model: mixed boundary value problem for Laplace equation We are interested in the total flux Φ Λ through the working surface: source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Λ C n = C Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

19 Spectroscopic impedance Mathematical model: mixed boundary value problem for Laplace equation We are interested in the total flux Φ Λ through the working surface: Cell impedance: source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Z cell (Λ) = C 0 Φ Λ Λ C n = C Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

20 Spectroscopic impedance Mathematical model: mixed boundary value problem for Laplace equation We are interested in the total flux Φ Λ through the working surface: Spectroscopic impedance: source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Λ C n = C Λ = D W Z(Λ) = Z cell (Λ) Z cell (0) working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

21 Spectroscopic impedance Dirichlet mode: Λ 0 source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Λ C n = C Z(Λ) = C 0 Φ Λ C 0 Φ 0 Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

22 Spectroscopic impedance Dirichlet mode: Λ 0 Absorbing surface source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Z(Λ) = C 0 Φ Λ C 0 Φ 0 C 0 Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

23 Spectroscopic impedance Dirichlet mode: Λ 0 Absorbing surface Z(Λ) Λ source C = C 0 C = 0 DL dir Φ Λ = D Diffusive transport Ω ds C Λ n Z(Λ) = C 0 Φ Λ C 0 Φ 0 L dir Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

24 Spectroscopic impedance Neumann mode: Λ source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Λ C n = C Z(Λ) = C 0 Φ Λ C 0 Φ 0 Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

25 Spectroscopic impedance Neumann mode: Λ Reflecting surface Z(Λ) source C = C 0 C = 0 Λ DL p Φ Λ = D Diffusive transport Ω ds C Λ n Z(Λ) = C 0 Φ Λ C 0 Φ 0 C n 0 Λ = D W L p working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

26 Spectroscopic impedance Intermediate mode: Constant phase angle behavior Z(Λ) Λ β source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Λ C n = C Z(Λ) = C 0 Φ Λ C 0 Φ 0 Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

27 Spectroscopic impedance Intermediate mode: Constant phase angle behavior Z(Λ) Λ β Fractal boundary β = 1 D f β = τ(2) D f source Φ Λ = D Diffusive transport C = C 0 Ω C = 0 ds C Λ n Λ C n = C Z(Λ) = C 0 Φ Λ C 0 Φ 0 Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

28 Spectroscopic impedance Intermediate mode: Arbitrary surface How does the geometrical irregularity influence C = 0 transport proper- ties? source Φ Λ = D Diffusive transport C = C 0 Ω ds C Λ n Λ C n = C Z(Λ) = C 0 Φ Λ C 0 Φ 0 Λ = D W working surface Denis S Grebenkov, CNRS Ecole Polytechnique p7/35

29 Plan Laplacian transport phenomena Theoretical approach Numerical and experimental study Conclusion, questions and perspectives Denis S Grebenkov, CNRS Ecole Polytechnique p8/35

30 Theoretical description Laplace equation + mixed condition: Description in average Denis S Grebenkov, CNRS Ecole Polytechnique p9/35

31 Theoretical description Laplace equation + mixed condition: Description in average Partially reflected Brownian motion: Random trajectory of each particle Denis S Grebenkov, CNRS Ecole Polytechnique p9/35

32 Theoretical description Laplace equation + mixed condition: Description in average Partially reflected Brownian motion: Random trajectory of each particle Dirichlet-to-Neumann operator: Spectral analysis of the problem Denis S Grebenkov, CNRS Ecole Polytechnique p9/35

33 Definition of the Dirichlet-to-Neumann operator For a given domain Ω, Ω Ω Denis S Grebenkov, CNRS Ecole Polytechnique p10/35

34 Definition of the Dirichlet-to-Neumann operator For a given domain Ω, one takes a function f H 1 ( Ω) on the boundary Ω, Ω Ω f Denis S Grebenkov, CNRS Ecole Polytechnique p10/35

35 Definition of the Dirichlet-to-Neumann operator For a given domain Ω, one takes a function f H 1 ( Ω) on the boundary Ω, one solves the Dirichlet problem in this domain Ω Ω u = f u = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p10/35

36 Definition of the Dirichlet-to-Neumann operator For a given domain Ω, one takes a function f H 1 ( Ω) on the boundary Ω, one solves the Dirichlet problem in this domain, and one applies the normal derivative to the solution Ω Ω u = f u = 0 f g = u n [Mf] Denis S Grebenkov, CNRS Ecole Polytechnique p10/35

37 Properties of the Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator M is dependent solely on the geometry; Denis S Grebenkov, CNRS Ecole Polytechnique p11/35

38 Properties of the Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator M is dependent solely on the geometry; self-adjoint; Denis S Grebenkov, CNRS Ecole Polytechnique p11/35

39 Properties of the Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator M is dependent solely on the geometry; self-adjoint; with discrete spectrum {µ α } bounded below; Denis S Grebenkov, CNRS Ecole Polytechnique p11/35

40 Properties of the Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator M is dependent solely on the geometry; self-adjoint; with discrete spectrum {µ α } bounded below; with smooth eigenfunctions V α (s) (defined on the boundary Ω) which form a complete basis in L 2 ( Ω) Denis S Grebenkov, CNRS Ecole Polytechnique p11/35

41 Representation of the spectroscopic impedance Z(Λ) = Z cell (Λ) Z cell (0) = C 0 Φ Λ C 0 Φ 0 Denis S Grebenkov, CNRS Ecole Polytechnique p12/35

42 Representation of the spectroscopic impedance Z(Λ) = Z cell (Λ) Z cell (0) = C 0 Φ Λ C 0 Φ 0 Z(Λ) = Λ D α=0 F α 1 + µ α Λ Denis S Grebenkov, CNRS Ecole Polytechnique p12/35

43 Representation of the spectroscopic impedance Z(Λ) = Z cell (Λ) Z cell (0) = C 0 Φ Λ C 0 Φ 0 Physics Λ, D Z(Λ) = Λ D Λ = D W, D K, r ρ α=0 F α 1 + µ α Λ Geometry µ α, F α F α = ( V α ω 0 ) 2 Denis S Grebenkov, CNRS Ecole Polytechnique p12/35

44 Plan Laplacian transport phenomena Theoretical approach Numerical and experimental study Conclusion, questions and perspectives Denis S Grebenkov, CNRS Ecole Polytechnique p13/35

45 Koch interface Denis S Grebenkov, CNRS Ecole Polytechnique p14/35

46 Koch interface Denis S Grebenkov, CNRS Ecole Polytechnique p14/35

47 Koch interface Denis S Grebenkov, CNRS Ecole Polytechnique p14/35

48 Koch interface D f = ln 5/ ln 3 Denis S Grebenkov, CNRS Ecole Polytechnique p14/35

49 Localization of the spectrum µ α L tot Denis S Grebenkov, CNRS Ecole Polytechnique p15/35

50 Localization of the spectrum µ α L tot Denis S Grebenkov, CNRS Ecole Polytechnique p15/35

51 Localization of the spectrum µ α L tot Denis S Grebenkov, CNRS Ecole Polytechnique p15/35

52 Localization of the spectrum µ α L tot Denis S Grebenkov, CNRS Ecole Polytechnique p15/35

53 Reduction of eigenmodes µ α L p Z(Λ) = Λ D α=0 F α 1 + µ α Λ Denis S Grebenkov, CNRS Ecole Polytechnique p16/35

54 Reduction of eigenmodes 05 F α L Z(Λ) Λ D ( F0 1 + µ 0 Λ + F µ 1 Λ ) µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p16/35

55 Eigenvectors First generation Plane surface (line) V 0 01 V F 0 L tot = 1 F 0 L tot = 1 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

56 Eigenvectors First generation Plane surface (line) V 1 01 V F 1 L tot = 067 F 1 L tot = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

57 Eigenvectors First generation Plane surface (line) V 2 01 V F 2 L tot = 0 F 2 L tot = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

58 Eigenvectors First generation Plane surface (line) V 3 01 V F 3 L tot = 0 F 3 L tot = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

59 Eigenvectors First generation Plane surface (line) V 4 01 V F 4 L tot = 0017 F 4 L tot = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

60 Eigenvectors First generation Plane surface (line) V 9 01 V F 9 L tot = 0021 F 9 L tot = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

61 Eigenvectors First generation Plane surface (line) V V F 50 L tot = 0 F 50 L tot = 0 Denis S Grebenkov, CNRS Ecole Polytechnique p17/35

62 Reduction of eigenmodes 05 0 g=1 0 S µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p18/35

63 Reduction of eigenmodes 05 g=1 0 g=2 0 S µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p18/35

64 Reduction of eigenmodes g=1 05 g=2 0 g=3 0 S µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p18/35

65 Reduction of eigenmodes g=1 g=2 g= S g= µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p18/35

66 Reduction of eigenmodes g=1 g=2 g= S g= ˆµ k L p 5 k µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p18/35

67 Reduction of eigenmodes g=1 g=2 g= S g= ˆµ k L p 5 k ˆFk L p (5/3) k µ α L p Denis S Grebenkov, CNRS Ecole Polytechnique p18/35

68 Spectroscopic impedance: model and calculations Fourth generation: 100% 50% (Z Z mod )/Z Z mod (Λ) = Λ D 4 k=0 ˆF k 1 + Λˆµ k 0 50% 100% Λ/L Denis S Grebenkov, CNRS Ecole Polytechnique p19/35

69 Analytical model: anomalous behavior Tenth generation: Z mod (Λ) = Λ D 10 k=0 ˆF k 1 + Λˆµ k Denis S Grebenkov, CNRS Ecole Polytechnique p20/35

70 Analytical model: anomalous behavior Tenth generation: 10 5 Z mod (Λ) = Λ D 10 k=0 DZ mod (Λ) Λ/L p ˆF k 1 + Λˆµ k DZ mod (Λ) Λ/L Z mod (Λ) Λ β β Λ/L Denis S Grebenkov, CNRS Ecole Polytechnique p20/35

71 Analytical model: anomalous behavior Tenth generation: β(λ) Z mod (Λ) = Λ D 10 k=0 ˆF k 1 + Λˆµ k g= β 068 1/D f = Λ Denis S Grebenkov, CNRS Ecole Polytechnique p20/35

72 Analytical model: anomalous behavior Tenth generation: β(λ) Z mod (Λ) = Λ D 10 k=0 ˆF k 1 + Λˆµ k g=10 g=5 07 β 068 1/D f = Λ Denis S Grebenkov, CNRS Ecole Polytechnique p20/35

73 Analytical model: anomalous behavior Tenth generation: β(λ) g=10 g=5 Z mod (Λ) = Λ D g=3 10 k=0 ˆF k 1 + Λˆµ k 07 β 068 1/D f = Λ Denis S Grebenkov, CNRS Ecole Polytechnique p20/35

74 Irregular boundaries Random Koch curves Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

75 Irregular boundaries Random Koch curves Koch curves of different dimensions Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

76 Irregular boundaries Random Koch curves Koch curves of different dimensions Koch surfaces Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

77 Irregular boundaries Random Koch curves Koch curves of different dimensions Koch surfaces Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

78 Irregular boundaries Random Koch curves Koch curves of different dimensions Koch surfaces Effect of eigenmodes reduction Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

79 Irregular boundaries Random Koch curves Koch curves of different dimensions Koch surfaces Effect of eigenmodes reduction Hierarchy of eigenvalues Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

80 Irregular boundaries Random Koch curves Koch curves of different dimensions Koch surfaces Effect of eigenmodes reduction Hierarchy of eigenvalues Analytical model of the impedance Denis S Grebenkov, CNRS Ecole Polytechnique p21/35

81 Experimental study How does the geometrical irregularity influence impedance behavior of real electrodes? Denis S Grebenkov, CNRS Ecole Polytechnique p22/35

82 Experimental study How does the geometrical irregularity influence impedance behavior of real electrodes? Denis S Grebenkov, CNRS Ecole Polytechnique p22/35

83 Classical theory of double layer Electric transport: V = V 0 e iωt V = 0 V = 0 V = V 0 V n = 1 Λ V Λ = ζ(ω) ρ in the volume on the counter-electrode on the working electrode where ζ(ω) is the surface impedance Z(Λ) = ρλ α= F α 1 + µ α Λ Denis S Grebenkov, CNRS Ecole Polytechnique p23/35

84 Practical implementation One measures the physical characteristics, ρ and ζ(ω), to obtain Λ = ζ(ω)/ρ; Denis S Grebenkov, CNRS Ecole Polytechnique p24/35

85 Practical implementation One measures the physical characteristics, ρ and ζ(ω), to obtain Λ = ζ(ω)/ρ; One finds the geometrical characteristics, µ α and F α ; Denis S Grebenkov, CNRS Ecole Polytechnique p24/35

86 Practical implementation One measures the physical characteristics, ρ and ζ(ω), to obtain Λ = ζ(ω)/ρ; One finds the geometrical characteristics, µ α and F α ; One substitutes them into the spectral decomposition Z(Λ) = ρλ α=0 F α 1 + µ α Λ Denis S Grebenkov, CNRS Ecole Polytechnique p24/35

87 Second generation: experiment and theory Re[Z(ω)] R (Ohm) Im[Z(ω)] (Ohm) ω (Hz) ω (Hz) Denis S Grebenkov, CNRS Ecole Polytechnique p25/35

88 Second generation: experiment and theory Re[Z(ω)] R (Ohm) Im[Z(ω)] (Ohm) ω (Hz) ω (Hz) Denis S Grebenkov, CNRS Ecole Polytechnique p25/35

89 Second generation: experiment and theory Re[Z(ω)] R (Ohm) Im[Z(ω)] (Ohm) ω (Hz) ω (Hz) In conclusion, without knowledge of the microscopic transport mechanism, the geometrical irregularity can be taken into account by our theoretical approach Denis S Grebenkov, CNRS Ecole Polytechnique p25/35

90 Theoretical results Theoretical description of Laplacian transport towards irregular interfaces based on the Dirichlet-to-Neumann operator and study of the partially reflected Brownian motion Implementation of the (discrete) boundary element method for contruction of the Dirichlet-to-Neumann operator Denis S Grebenkov, CNRS Ecole Polytechnique p26/35

91 Numerical results Effective contribution of a small number of eigenmodes to linear response of the system Hierarchy of these significative eigenmodes (characteristic scales of the boundary) Analytical model of the impedance Denis S Grebenkov, CNRS Ecole Polytechnique p27/35

92 Experimental results Measurement of the impedance of brass and nickel electrodes with boundaries reproducing two generations of the convex Von Koch surface Observation of anomalous behavior of the impedance in the case of planar electrodes Accounting for geometrical irregularity by theoretical approach Denis S Grebenkov, CNRS Ecole Polytechnique p28/35

93 Construction of the Dirichlet-to- Neumann operator Present approach: A kind of the boundary element method based on potential theory Denis S Grebenkov, CNRS Ecole Polytechnique p29/35

94 Construction of the Dirichlet-to- Neumann operator Present approach: A kind of the boundary element method based on potential theory Main problem: one has to calculate all the elements of the matrix representing this operator in a discrete form Denis S Grebenkov, CNRS Ecole Polytechnique p29/35

95 Construction of the Dirichlet-to- Neumann operator Present approach: A kind of the boundary element method based on potential theory Main problem: one has to calculate all the elements of the matrix representing this operator in a discrete form Reality: Only a small fraction of the eigenmodes is in fact needed Denis S Grebenkov, CNRS Ecole Polytechnique p29/35

96 Construction of the Dirichlet-to- Neumann operator We are looking for a numerical method allowing to compute the eigenmodes of the operator M successively, without constructing the whole operator Denis S Grebenkov, CNRS Ecole Polytechnique p30/35

97 Construction of the Dirichlet-to- Neumann operator Denis S Grebenkov, CNRS Ecole Polytechnique p30/35

98 Perspectives Extensions of the theoretical approach: Denis S Grebenkov, CNRS Ecole Polytechnique p31/35

99 Perspectives Extensions of the theoretical approach: other equations (time-dependent diffusion) Denis S Grebenkov, CNRS Ecole Polytechnique p31/35

100 Perspectives Extensions of the theoretical approach: other equations (time-dependent diffusion) other conditions (non linear) Denis S Grebenkov, CNRS Ecole Polytechnique p31/35

101 Perspectives Extensions of the theoretical approach Numerical study of higher orders generations, especially in 3 dimensions Denis S Grebenkov, CNRS Ecole Polytechnique p32/35

102 Perspectives Extensions of the theoretical approach Numerical study of higher orders generations, especially in 3 dimensions Analytical model of the impedance for random interfaces? Denis S Grebenkov, CNRS Ecole Polytechnique p32/35

103 Perspectives Extensions of the theoretical approach Numerical study of higher orders generations, especially in 3 dimensions Analytical model of the impedance for random interfaces? Inverse problem: knowing the spectroscopic impedance, what can one say about the geometry? Denis S Grebenkov, CNRS Ecole Polytechnique p32/35

104 Perspectives Extensions of the theoretical approach Numerical study of higher orders generations, especially in 3 dimensions Analytical model of the impedance for random interfaces? Inverse problem: knowing the spectroscopic impedance, what can one say about the geometry? Experimental study of irregular electrodes Denis S Grebenkov, CNRS Ecole Polytechnique p32/35

105 Perspectives Denis S Grebenkov, CNRS Ecole Polytechnique p33/35

106 Open problems Justification of the approach for really irregular interfaces Denis S Grebenkov, CNRS Ecole Polytechnique p34/35

107 Open problems Justification of the approach for really irregular interfaces Reduction of the number of contributing eigenmodes Denis S Grebenkov, CNRS Ecole Polytechnique p34/35

108 Open problems Justification of the approach for really irregular interfaces Reduction of the number of contributing eigenmodes Relation between the eigenvalues and the characteristic length scales of the boundary Denis S Grebenkov, CNRS Ecole Polytechnique p34/35

109 Open problems Justification of the approach for really irregular interfaces Reduction of the number of contributing eigenmodes Relation between the eigenvalues and the characteristic length scales of the boundary Correction, generalization and mathematical justification for the simplified impedance model β = 1 = β = τ(2) D f D f Denis S Grebenkov, CNRS Ecole Polytechnique p34/35

110 Related papers 1 B Sapoval, Transport Across Irregular Interfaces: Fractal Electrodes, Membranes and Catalysts, in Fractals and Disordered Systems, Eds A Bunde, S Havlin, pp (Springer, 1996) 2 B Sapoval, General Formulation of Laplacian Transfer Across Irregular Surfaces, Phys Rev Lett 73, 3314 (1994) 3 M Filoche and B Sapoval, Can One Hear the Shape of an Electrode? II Theoretical Study of the Laplacian Transfer, Eur Phys J B 9, 755 (1999) 4 D S Grebenkov, M Filoche, B Sapoval, Spectral Properties of the Brownian Self-Transport Operator, Eur Phys J B 36, (2003) 5 D S Grebenkov, M Filoche, B Sapoval, Mathematical Basis for a General Theory of Laplacian Transport towards Irregular Interfaces, Phys Rev E 73, (2006) 6 D S Grebenkov, M Filoche, B Sapoval, A Simplified Analytical Model for Laplacian Transfer Across Deterministic Prefractal Interfaces, Fractals (2007, in press) 7 D S Grebenkov, Partially Reflected Brownian Motion: A Stochastic Approach to Transport Phenomena, in Focus on Probability Theory, Ed L R Velle, pp (Nova Science Publishers, 2006) [Electronic version is available on Denis S Grebenkov, CNRS Ecole Polytechnique p35/35

Outline of the talk How to describe restricted diffusion? How to monitor restricted diffusion? Laplacian eigenfunctions in NMR Other applications Loca

Outline of the talk How to describe restricted diffusion? How to monitor restricted diffusion? Laplacian eigenfunctions in NMR Other applications Loca Laplacian Eigenfunctions in NMR Denis S. Grebenkov Laboratoire de Physique de la Matière Condensée CNRS Ecole Polytechnique, Palaiseau, France IPAM Workshop «Laplacian Eigenvalues and Eigenfunctions» February

More information

INTRODUCTION I. QUALITY FACTOR AND AMPLITUDE DISTRIBUTION

INTRODUCTION I. QUALITY FACTOR AND AMPLITUDE DISTRIBUTION Acoustical properties of irregular and fractal cavities B. Sapoval, O. Haeberlé, and S. Russ Laboratoire de Physique de la Matière Condensée, Ecole Polytechnique, C.N.R.S., 91128 Palaiseau Cédex, France

More information

arxiv: v1 [math-ph] 28 Feb 2013

arxiv: v1 [math-ph] 28 Feb 2013 On the exponential decay of Laplacian eigenfunctions in planar domains with branches Binh T Nguyen 1, Andrey L Delytsin and Denis S Grebenkov 1 November 1, 018 1 Laboratoire de Physique de la Matière Condensée

More information

On the Exponential Decay of Laplacian Eigenfunctions in Planar Domains with Branches

On the Exponential Decay of Laplacian Eigenfunctions in Planar Domains with Branches Contemporary Mathematics Volume 630, 2014 http://dx.doi.org/10.1090/conm/630/12674 On the Exponential Decay of Laplacian Eigenfunctions in Planar Domains with Branches Binh T. Nguyen, Denis S. Grebenkov,

More information

One dimensional steady state diffusion, with and without source. Effective transfer coefficients

One dimensional steady state diffusion, with and without source. Effective transfer coefficients One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport

More information

Detailed Program of the Workshop

Detailed Program of the Workshop Detailed Program of the Workshop Inverse Problems: Theoretical and Numerical Aspects Laboratoire de Mathématiques de Reims, December 17-19 2018 December 17 14h30 Opening of the Workhop 15h00-16h00 Mourad

More information

Dynamics of a tagged monomer: Effects of elastic pinning and harmonic absorption. Shamik Gupta

Dynamics of a tagged monomer: Effects of elastic pinning and harmonic absorption. Shamik Gupta : Effects of elastic pinning and harmonic absorption Laboratoire de Physique Théorique et Modèles Statistiques, Université Paris-Sud, France Joint work with Alberto Rosso Christophe Texier Ref.: Phys.

More information

Modeling as a tool for understanding the MEA. Henrik Ekström Utö Summer School, June 22 nd 2010

Modeling as a tool for understanding the MEA. Henrik Ekström Utö Summer School, June 22 nd 2010 Modeling as a tool for understanding the MEA Henrik Ekström Utö Summer School, June 22 nd 2010 COMSOL Multiphysics and Electrochemistry Modeling The software is based on the finite element method A number

More information

What means dimension?

What means dimension? What means dimension? Christiane ROUSSEAU Universite de Montre al November 2011 How do we measure the size of a geometric object? For subsets of the plane we often use perimeter, length, area, diameter,

More information

Stochastic nonlinear Schrödinger equations and modulation of solitary waves

Stochastic 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 information

Physique de la Matière Condensée, Ecole Polytechnique, Palaiseau, France

Physique de la Matière Condensée, Ecole Polytechnique, Palaiseau, France Denis Grebenkov a and Gleb Oshanin b a Physique de la Matière Condensée, Ecole Polytechnique, Palaiseau, France b Physique Théorique de la Matière Condensée (UMR CNRS) Sorbonne Universités Université Pierre

More information

Modelling crystal growth with the phase-field method Mathis Plapp

Modelling crystal growth with the phase-field method Mathis Plapp Modelling crystal growth with the phase-field method Mathis Plapp Laboratoire de Physique de la Matière Condensée CNRS/Ecole Polytechnique, 91128 Palaiseau, France Solidification microstructures Hexagonal

More information

Boundary value problems

Boundary value problems 1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the

More information

2. The Steady State and the Diffusion Equation

2. The Steady State and the Diffusion Equation 2. The Steady State and the Diffusion Equation The Neutron Field Basic field quantity in reactor physics is the neutron angular flux density distribution: Φ( r r, E, r Ω,t) = v(e)n( r r, E, r Ω,t) -- distribution

More information

Forced Convectional Heat Transfer in Solid Oxide Fuel Cells: An Analytical Treatment.

Forced Convectional Heat Transfer in Solid Oxide Fuel Cells: An Analytical Treatment. Ionics 9 (2003) 83 Forced Convectional Heat Transfer in Solid Oxide Fuel Cells: An Analytical Treatment. F.A. Coutelieris 1, A.K. Demin 2, S.L. Douvartzides 1 and P.E. Tsiakaras 1 1 University of Thessalia,

More information

LAPLACIANS COMPACT METRIC SPACES. Sponsoring. Jean BELLISSARD a. Collaboration:

LAPLACIANS COMPACT METRIC SPACES. Sponsoring. Jean BELLISSARD a. Collaboration: LAPLACIANS on Sponsoring COMPACT METRIC SPACES Jean BELLISSARD a Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics Collaboration: I. PALMER (Georgia Tech, Atlanta) a e-mail:

More information

Fractional Laplacian

Fractional 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 information

State-Space Modeling of Electrochemical Processes. Michel Prestat

State-Space Modeling of Electrochemical Processes. Michel Prestat State-Space Modeling of Electrochemical Processes Who uses up my battery power? Michel Prestat ETH-Zürich Institute for Nonmetallic Materials Head: Prof. L.J. Gauckler Outline Electrochemistry Electrochemical

More information

Pulse propagation in random waveguides with turning points and application to imaging

Pulse propagation in random waveguides with turning points and application to imaging Pulse propagation in random waveguides with turning points and application to imaging Liliana Borcea Mathematics, University of Michigan, Ann Arbor and Josselin Garnier Centre de Mathématiques Appliquées,

More information

ON THE ROBIN EIGENVALUES OF THE LAPLACIAN IN THE EXTERIOR OF A CONVEX POLYGON

ON THE ROBIN EIGENVALUES OF THE LAPLACIAN IN THE EXTERIOR OF A CONVEX POLYGON NANOSYSTES: PHYSICS, CHEISTRY, ATHEATICS, 15, 6 (1, P. 46 56 ON THE ROBIN EIGENVALUES OF THE LAPLACIAN IN THE EXTERIOR OF A CONVEX POLYGON Konstantin Pankrashkin Laboratoire de mathématiques, Université

More information

1 Eigenvalues and eigenvectors of the Laplacian on a parallelepiped

1 Eigenvalues and eigenvectors of the Laplacian on a parallelepiped Eigenvalues and eigenvectors of Laplacian on a parallelepiped I. Ţevy, C. Udrişte, F. Munteanu, D. Zugrăvescu, L. Petrescu Abstract. The purpose of this paper is threefold: (i) to prove new properties

More information

Construction of a New Domain Decomposition Method for the Stokes Equations

Construction of a New Domain Decomposition Method for the Stokes Equations Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,

More information

Sound Propagation in the Nocturnal Boundary Layer. Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert

Sound Propagation in the Nocturnal Boundary Layer. Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert Sound Propagation in the Nocturnal Boundary Layer Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert The Propagation of Sound Outdoors (over flat ground) The atmosphere is a gas under the influence

More information

Reverse-Engineering the Robustness of Mammalian Lungs

Reverse-Engineering the Robustness of Mammalian Lungs 11 Reverse-Engineering the Robustness of Mammalian Lungs Michael Mayo 1, Peter Pfeifer 2 and Chen Hou 3 1 Environmental Laboratory, US Army Engineer Research and Development Center, Vicksburg, MS 2 Department

More information

INFLUENCE OF ROOM GEOMETRY ON MODAL DENSITY, SPATIAL DISTRIBUTION OF MODES AND THEIR DAMPING M. MEISSNER

INFLUENCE OF ROOM GEOMETRY ON MODAL DENSITY, SPATIAL DISTRIBUTION OF MODES AND THEIR DAMPING M. MEISSNER ARCHIVES OF ACOUSTICS 30, 4 (Supplement), 203 208 (2005) INFLUENCE OF ROOM GEOMETRY ON MODAL DENSITY, SPATIAL DISTRIBUTION OF MODES AND THEIR DAMPING M. MEISSNER Institute of Fundamental Technological

More information

ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS

ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS VINCENT FONTAINE 1,2 AND ANIS YOUNES 2 1 Laboratoire de Physique des Bâtiments et des Systèmes, Université de La Réunion, 15 avenue

More information

Cube-based diffusion algorithm

Cube-based diffusion algorithm The Open-Access Journal for the Basic Principles of Diffusion Theory, Experiment and Application Cube-based diffusion algorithm Stephan Zschiegner, 1,2 Stefanie Russ, 3 Armin Bunde, 1 Jörg Kärger 2 1 Universität

More information

Microscopic-Macroscopic connection. Silvana Botti

Microscopic-Macroscopic connection. Silvana Botti relating experiment and theory European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre for Computational

More information

L -uniqueness of Schrödinger operators on a Riemannian manifold

L -uniqueness of Schrödinger operators on a Riemannian manifold L -uniqueness of Schrödinger operators on a Riemannian manifold Ludovic Dan Lemle Abstract. The main purpose of this paper is to study L -uniqueness of Schrödinger operators and generalized Schrödinger

More information

Modeling, simulation and characterization of atomic force microscopy measurements for ionic transport and impedance in PEM fuel cells

Modeling, simulation and characterization of atomic force microscopy measurements for ionic transport and impedance in PEM fuel cells Modeling, simulation and characterization of atomic force microscopy measurements for ionic transport and impedance in PEM fuel cells Peter M. Pinsky David M. Barnett Yongxing Shen Department of Mechanical

More information

Surfaces and Interfaces

Surfaces and Interfaces 1/16 Surfaces and Interfaces Fouad MAROUN Laboratoire de Physique de la Matière Condensée, Ecole Polytechnique CNRS Palaiseau, France 2/16 Outline Definition of surfaces and interfaces and 3 examples in

More information

Review of temperature distribution in cathode of PEMFC

Review of temperature distribution in cathode of PEMFC Project Report 2008 MVK 160 Heat and Mass Transport May 08, 2008, Lund, Sweden Review of temperature distribution in cathode of PEMFC Munir Ahmed Khan Department of Energy Sciences, Lund Institute of Technology,

More information

Shlomo Havlin } Anomalous Transport in Scale-free Networks, López, et al,prl (2005) Bar-Ilan University. Reuven Cohen Tomer Kalisky Shay Carmi

Shlomo Havlin } Anomalous Transport in Scale-free Networks, López, et al,prl (2005) Bar-Ilan University. Reuven Cohen Tomer Kalisky Shay Carmi Anomalous Transport in Complex Networs Reuven Cohen Tomer Kalisy Shay Carmi Edoardo Lopez Gene Stanley Shlomo Havlin } } Bar-Ilan University Boston University Anomalous Transport in Scale-free Networs,

More information

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties

More information

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

BV functions in a Gelfand triple and the stochastic reflection problem on a convex set BV functions in a Gelfand triple and the stochastic reflection problem on a convex set Xiangchan Zhu Joint work with Prof. Michael Röckner and Rongchan Zhu Xiangchan Zhu ( Joint work with Prof. Michael

More information

Statistical mechanics of random billiard systems

Statistical mechanics of random billiard systems Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore

More information

A Shrinking Core Model for Steam Hydration of CaO-Based Sorbents Cycled for CO2 Capture: Supplementary Information

A Shrinking Core Model for Steam Hydration of CaO-Based Sorbents Cycled for CO2 Capture: Supplementary Information A Shrinking Core Model for Steam Hydration of CaO-Based Sorbents Cycled for CO2 Capture: Supplementary Information John Blamey 1, Ming Zhao 2, Vasilije Manovic 3, Edward J. Anthony 3, Denis R. Dugwell

More information

A class of domains with fractal boundaries: Functions spaces and numerical methods

A class of domains with fractal boundaries: Functions spaces and numerical methods A class of domains with fractal boundaries: Functions spaces and numerical methods Yves Achdou joint work with T. Deheuvels and N. Tchou Laboratoire J-L Lions, Université Paris Diderot École Centrale -

More information

Spectral asymptotics for stable trees and the critical random graph

Spectral asymptotics for stable trees and the critical random graph Spectral asymptotics for stable trees and the critical random graph EPSRC SYMPOSIUM WORKSHOP DISORDERED MEDIA UNIVERSITY OF WARWICK, 5-9 SEPTEMBER 2011 David Croydon (University of Warwick) Based on joint

More information

Asymptotic distribution of eigenvalues of Laplace operator

Asymptotic distribution of eigenvalues of Laplace operator Asymptotic distribution of eigenvalues of Laplace operator 23.8.2013 Topics We will talk about: the number of eigenvalues of Laplace operator smaller than some λ as a function of λ asymptotic behaviour

More information

Fundamental Solutions and Green s functions. Simulation Methods in Acoustics

Fundamental Solutions and Green s functions. Simulation Methods in Acoustics Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental

More information

Contact time of a bouncing drop

Contact time of a bouncing drop Contact time of a bouncing drop Denis Richard, Christophe Clanet (*) & David Quéré Laboratoire de Physique de la Matière Condensée, URA 792 du CNRS, Collège de France, 75231 Paris Cedex 05, France (*)

More information

Inverse problems in tumor growth modelling

Inverse problems in tumor growth modelling MC2 Damiano Lombardi Angelo Iollo Thierry Colin Olivier Saut Marseille, 18-09-09 Outline 1) Introduction: phenomenological model for tumor growth models 2) Models based on mixture theory: a simplified

More information

Some Mathematical and Physical Background

Some Mathematical and Physical Background Some Mathematical and Physical Background Linear partial differential operators Let H be a second-order, elliptic, self-adjoint PDO, on scalar functions, in a d-dimensional region Prototypical categories

More information

Lecture 12: Detailed balance and Eigenfunction methods

Lecture 12: Detailed balance and Eigenfunction methods Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner

More information

Electrostatics: Electrostatic Devices

Electrostatics: Electrostatic Devices Electrostatics: Electrostatic Devices EE331 Electromagnetic Field Theory Outline Laplace s Equation Derivation Meaning Solving Laplace s equation Resistors Capacitors Electrostatics -- Devices Slide 1

More information

On the spectrum of the Hodge Laplacian and the John ellipsoid

On the spectrum of the Hodge Laplacian and the John ellipsoid Banff, July 203 On the spectrum of the Hodge Laplacian and the John ellipsoid Alessandro Savo, Sapienza Università di Roma We give upper and lower bounds for the first eigenvalue of the Hodge Laplacian

More information

Quantitative Electrophysiology

Quantitative Electrophysiology ECE 795: Quantitative Electrophysiology Notes for Lecture #1 Wednesday, September 13, 2006 1. INTRODUCTION TO EXCITABLE CELLS Historical perspective: Bioelectricity first discovered by Luigi Galvani in

More information

The fractal nature of a diffusion front and the relation to percolation

The fractal nature of a diffusion front and the relation to percolation The fractal nature of a diffusion front and the relation to percolation Bernard Sapoval, Michel Rosso, JeanFrançois Gouyet To cite this version: Bernard Sapoval, Michel Rosso, JeanFrançois Gouyet. The

More information

Anderson Localization on the Sierpinski Gasket

Anderson Localization on the Sierpinski Gasket Anderson Localization on the Sierpinski Gasket G. Mograby 1 M. Zhang 2 1 Department of Physics Technical University of Berlin, Germany 2 Department of Mathematics Jacobs University, Germany 5th Cornell

More information

Physical Pharmacy. Diffusion

Physical Pharmacy. Diffusion Physical Pharmacy Diffusion Diffusion Diffusion is defined as a process of mass transfer of individual molecules of a substance brought about by random molecular motion and associated with a driving force

More information

Pointwise bounds on the gradient and the spectrum of the Neumann-Poincaré operator: The case of 2 discs

Pointwise bounds on the gradient and the spectrum of the Neumann-Poincaré operator: The case of 2 discs Pointwise bounds on the gradient and the spectrum of the Neumann-Poincaré operator: The case of discs Eric Bonnetier and Faouzi Triki Abstract. We compute the spectrum of the Neumann-Poincaré operator

More information

Data Analysis and Manifold Learning Lecture 3: Graphs, Graph Matrices, and Graph Embeddings

Data Analysis and Manifold Learning Lecture 3: Graphs, Graph Matrices, and Graph Embeddings Data Analysis and Manifold Learning Lecture 3: Graphs, Graph Matrices, and Graph Embeddings Radu Horaud INRIA Grenoble Rhone-Alpes, France Radu.Horaud@inrialpes.fr http://perception.inrialpes.fr/ Outline

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods

More information

A LATTICE GAS MODEL OF ELECTROCHEMICAL CELLS : MEAN-FIELD KINETIC APPROACH

A LATTICE GAS MODEL OF ELECTROCHEMICAL CELLS : MEAN-FIELD KINETIC APPROACH A LATTICE GAS MODEL OF ELECTROCHEMICAL CELLS : MEAN-FIELD KINETIC APPROACH M.-O. BERNARD, M. PLAPP, J.-F. GOUYET Laboratoire de Physique de la Matière Condensée, Ecole Polytechnique, F-91128 Palaiseau,

More information

TRANSPORT IN POROUS MEDIA

TRANSPORT IN POROUS MEDIA 1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case

More information

Signal Loss. A1 A L[Neper] = ln or L[dB] = 20log 1. Proportional loss of signal amplitude with increasing propagation distance: = α d

Signal Loss. A1 A L[Neper] = ln or L[dB] = 20log 1. Proportional loss of signal amplitude with increasing propagation distance: = α d Part 6 ATTENUATION Signal Loss Loss of signal amplitude: A1 A L[Neper] = ln or L[dB] = 0log 1 A A A 1 is the amplitude without loss A is the amplitude with loss Proportional loss of signal amplitude with

More information

Standard Finite Elements and Weighted Regularization

Standard Finite Elements and Weighted Regularization Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the

More information

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 23 Feb 1998

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 23 Feb 1998 Multiplicative processes and power laws arxiv:cond-mat/978231v2 [cond-mat.stat-mech] 23 Feb 1998 Didier Sornette 1,2 1 Laboratoire de Physique de la Matière Condensée, CNRS UMR6632 Université des Sciences,

More information

Elements of linear algebra

Elements of linear algebra Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v

More information

Quantitative Electrophysiology

Quantitative Electrophysiology ECE 795: Quantitative Electrophysiology Notes for Lecture #1 Tuesday, September 18, 2012 1. INTRODUCTION TO EXCITABLE CELLS Historical perspective: Bioelectricity first discovered by Luigi Galvani in 1780s

More information

A Porous Media Approach for Bifurcating Flow and Mass Transfer in a Human Lung

A Porous Media Approach for Bifurcating Flow and Mass Transfer in a Human Lung Fujio Kuwahara Yoshihiko Sano Department of Mechanical Engineering, Shizuoka University, 3-5-1 Johoku, Hamamatsu 43-8561, Japan Jianjun Liu Department of Civil Engineering, Wuhan Polytechnic University,

More information

On continuous time contract theory

On continuous time contract theory Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem

More information

IFT CONTINUOUS LAPLACIAN Mikhail Bessmeltsev

IFT CONTINUOUS LAPLACIAN   Mikhail Bessmeltsev IFT 6112 07 CONTINUOUS LAPLACIAN http://www-labs.iro.umontreal.ca/~bmpix/teaching/6112/2018/ Mikhail Bessmeltsev Famous Motivation An Experiment Unreasonable to Ask? Length of string http://www.takamine.com/templates/default/images/gclassical.png

More information

Some Planar Isospectral Domains. Peter Buser, John Conway, Peter Doyle, and Klaus-Dieter Semmler. 1 Introduction

Some Planar Isospectral Domains. Peter Buser, John Conway, Peter Doyle, and Klaus-Dieter Semmler. 1 Introduction IMRN International Mathematics Research Notices 1994, No. 9 Some Planar Isospectral Domains Peter Buser, John Conway, Peter Doyle, and Klaus-Dieter Semmler 1 Introduction In 1965, Mark Kac [6] asked, Can

More information

MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued).

MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued). MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued). Bessel s differential equation of order m 0: z 2 d2 f dz 2 + z df dz + (z2 m 2 )f = 0 The equation is considered on the

More information

B. LAQUERRIÈRE AND N. PRIVAULT

B. LAQUERRIÈRE AND N. PRIVAULT DEVIAION INEQUALIIES FOR EXPONENIAL JUMP-DIFFUSION PROCESSES B. LAQUERRIÈRE AND N. PRIVAUL Abstract. In this note we obtain deviation inequalities for the law of exponential jump-diffusion processes at

More information

Electrodes MB - JASS 09. Metal in electrolyte

Electrodes MB - JASS 09. Metal in electrolyte Electrodes MB - JASS 09 Metal in electrolyte 1 Helmholtz double layer (1) Helmholtz double layer : simplest approximation surface charge is neutralized by opposite signed counterions placed away from the

More information

Title of communication, titles not fitting in one line will break automatically

Title of communication, titles not fitting in one line will break automatically Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add

More information

A Finite Element Model for Numerical Analysis of Sintering

A Finite Element Model for Numerical Analysis of Sintering A Finite Element Model for Numerical Analysis of Sintering DANIELA CÂRSTEA High-School Group of Railways, Craiova ION CÂRSTEA Department of Computer Engineering and Communication University of Craiova

More information

Chapter 4. RWs on Fractals and Networks.

Chapter 4. RWs on Fractals and Networks. Chapter 4. RWs on Fractals and Networks. 1. RWs on Deterministic Fractals. 2. Linear Excitation on Disordered lattice; Fracton; Spectral dimension 3. RWs on disordered lattice 4. Random Resistor Network

More information

The Truth about diffusion (in liquids)

The Truth about diffusion (in liquids) The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion

More information

Theoretical questions for the final exam ED 2012.

Theoretical questions for the final exam ED 2012. Theoretical questions for the final exam ED 2012. 1. Radiation a) Properties and types of radiation b) Physical parameters of radiation 2. Law of attenuation of radiation a) Experimental interpretation

More information

Absorption of gas by a falling liquid film

Absorption of gas by a falling liquid film Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German

More information

Deterministic chaos, fractals and diffusion: From simple models towards experiments

Deterministic chaos, fractals and diffusion: From simple models towards experiments Deterministic chaos, fractals and diffusion: From simple models towards experiments Rainer Klages Queen Mary University of London, School of Mathematical Sciences Université Pierre et Marie Curie Paris,

More information

Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating

Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating by T.C. Luce in collaboration with J.C. DeBoo, C.C. Petty, J. Pino, J.M. Nelson, and J.C.M. dehaas Presented at 9th EU-US Transport

More information

Electric Field--Definition. Brownian motion and drift velocity

Electric Field--Definition. Brownian motion and drift velocity Electric Field--Definition Definition of electrostatic (electrical) potential, energy diagram and how to remember (visualize) relationships E x Electrons roll downhill (this is a definition ) Holes are

More information

New constructions of domain decomposition methods for systems of PDEs

New constructions of domain decomposition methods for systems of PDEs New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.

More information

Electrochemical methods : Fundamentals and Applications

Electrochemical methods : Fundamentals and Applications Electrochemical methods : Fundamentals and Applications Lecture Note 7 May 19, 2014 Kwang Kim Yonsei University kbkim@yonsei.ac.kr 39 8 7 34 53 Y O N Se I 88.91 16.00 14.01 78.96 126.9 Electrochemical

More information

Introduction to the mathematical modeling of multi-scale phenomena

Introduction to the mathematical modeling of multi-scale phenomena Introduction to the mathematical modeling of multi-scale phenomena Diffusion Brownian motion Brownian motion (named after botanist Robert Brown) refers to the random motion of particles suspended in a

More information

Theory of Charge Transport in Mixed Conductors: Description of Interfacial Contributions Compatible with the Gibbs Thermodynamics

Theory of Charge Transport in Mixed Conductors: Description of Interfacial Contributions Compatible with the Gibbs Thermodynamics Theory of Charge Transport in Mixed Conductors: Description of Interfacial Contributions Compatible with the Gibbs Thermodynamics Mikhail A. Vorotyntsev LSEO-UMR 5188 CNRS, Université de Bourgogne, Dijon,

More information

Dissipative Solitons in Physical Systems

Dissipative Solitons in Physical Systems Dissipative Solitons in Physical Systems Talk given by H.-G. Purwins Institute of Applied Physics University of Münster, Germany Chapter 1 Introduction Cosmology and Pattern Formation 1-1 Complex Behaviour

More information

Applicability of the High Field Model: A Preliminary Numerical Study

Applicability of the High Field Model: A Preliminary Numerical Study Applicability of the High Field Model: A Preliminary Numerical Study Carlo Cercignani,IreneM.Gamba, Joseph W. Jerome, and Chi-Wang Shu Abstract In a companion presentation, we have discussed the theory

More information

The Klein-Gordon Equation Meets the Cauchy Horizon

The Klein-Gordon Equation Meets the Cauchy Horizon Enrico Fermi Institute and Department of Physics University of Chicago University of Mississippi May 10, 2005 Relativistic Wave Equations At the present time, our best theory for describing nature is Quantum

More information

Local vs. Nonlocal Diffusions A Tale of Two Laplacians

Local 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 information

A Multigrid Method for Two Dimensional Maxwell Interface Problems

A Multigrid Method for Two Dimensional Maxwell Interface Problems A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A

More information

HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM

HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM Homogenization of a heat transfer problem 1 G. Allaire HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM Work partially supported by CEA Grégoire ALLAIRE, Ecole Polytechnique

More information

Quantifying Intermittent Transport in Cell Cytoplasm

Quantifying Intermittent Transport in Cell Cytoplasm Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent

More information

Design of Structures for Earthquake Resistance

Design of Structures for Earthquake Resistance NATIONAL TECHNICAL UNIVERSITY OF ATHENS Design of Structures for Earthquake Resistance Basic principles Ioannis N. Psycharis Lecture 3 MDOF systems Equation of motion M u + C u + K u = M r x g(t) where:

More information

Final Exam May 4, 2016

Final 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 information

On semilinear elliptic equations with measure data

On 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 information

Lecture 12: Detailed balance and Eigenfunction methods

Lecture 12: Detailed balance and Eigenfunction methods Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),

More information

Mean-Field Limits for Large Particle Systems Lecture 2: From Schrödinger to Hartree

Mean-Field Limits for Large Particle Systems Lecture 2: From Schrödinger to Hartree for Large Particle Systems Lecture 2: From Schrödinger to Hartree CMLS, École polytechnique & CNRS, Université Paris-Saclay FRUMAM, Marseilles, March 13-15th 2017 A CRASH COURSE ON QUANTUM N-PARTICLE DYNAMICS

More information

Divorcing pressure from viscosity in incompressible Navier-Stokes dynamics

Divorcing pressure from viscosity in incompressible Navier-Stokes dynamics Divorcing pressure from viscosity in incompressible Navier-Stokes dynamics Jian-Guo Liu 1, Jie Liu 2, and Robert L. Pego 3 February 18, 2005 Abstract We show that in bounded domains with no-slip boundary

More information

Fractals and Disordered Systems

Fractals and Disordered Systems Fractals and Disordered Systems Springer Berlin Heidelberg New York Barcelona Budapest Hong Kong London Milan Paris Santa Clara Singapore Tokyo Armin Bunde Shlomo Havlin (Eds.) Fractals and Disordered

More information

Corner Singularities

Corner Singularities Corner Singularities Martin Costabel IRMAR, Université de Rennes 1 Analysis and Numerics of Acoustic and Electromagnetic Problems Linz, 10 15 October 2016 Martin Costabel (Rennes) Corner Singularities

More information

38 Green's Function: Diffusion Equation

38 Green's Function: Diffusion Equation 38 Green's Function: Diffusion Equation The Green's function method to solve the general initialboundary value problem for diffusion equations is given. The Markovian property of the free-space Green's

More information

ES 272 Assignment #2. in,3

ES 272 Assignment #2. in,3 ES 272 Assignment #2 Due: March 14th, 2014; 5pm sharp, in the dropbox outside MD 131 (Donhee Ham office) Instructor: Donhee Ham (copyright c 2014 by D. Ham) (Problem 1) The kt/c Noise (50pt) Imagine an

More information

GEOMETRIC AND FRACTAL PROPERTIES OF SCHRAMM-LOEWNER EVOLUTION (SLE)

GEOMETRIC AND FRACTAL PROPERTIES OF SCHRAMM-LOEWNER EVOLUTION (SLE) GEOMETRIC AND FRACTAL PROPERTIES OF SCHRAMM-LOEWNER EVOLUTION (SLE) Triennial Ahlfors-Bers Colloquium Gregory F. Lawler Department of Mathematics Department of Statistics University of Chicago 5734 S.

More information

Physics 506 Winter 2004

Physics 506 Winter 2004 Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading

More information