A Newton Krylov method for coupling transport with chemistry in porous media
|
|
- Theresa Willis
- 5 years ago
- Views:
Transcription
1 A Newton Krylov method for coupling transport with chemistry in porous media Laila Amir 1,2 Michel Kern 2 1 Itasca Consultants, Lyon 2 INRIA Rocquencourt Joint work with J. Erhel (INRIA Rennes /IRISA) Funded by GDR MoMaS Groupe de travail «Méthodes numériques» Laboratoire Jacques-Louis Lions 4 décembre 2006 L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 1 / 41
2 Outline 1 Modeling Chemical phenomena Chemical phenomena in aqueous chemistry Modeling equilibrium systems Sorption models 2 Multispecies equilibrium reactive transport Chemical problem Transport equations The coupled system Coupling algorithms 3 Numerical results Pyrite test case Chromatography Simplified MoMaS benchmark L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 2 / 41
3 Outline 1 Modeling Chemical phenomena Chemical phenomena in aqueous chemistry Modeling equilibrium systems Sorption models 2 Multispecies equilibrium reactive transport Chemical problem Transport equations The coupled system Coupling algorithms 3 Numerical results Pyrite test case Chromatography Simplified MoMaS benchmark L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 3 / 41
4 Reactive transport in a porous medium L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 4 / 41
5 Classification of chemical reactions According to nature of reaction Homogeneous In the same phase (aqueous, gaseous,...) Examples: Acid base, oxydo reduction According to speed of reaction L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 5 / 41
6 Classification of chemical reactions According to nature of reaction Homogeneous In the same phase (aqueous, gaseous,...) Examples: Acid base, oxydo reduction Heterogeneous Involve different phases Examples: Sorption, precipitation / dissolution,... According to speed of reaction L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 5 / 41
7 Classification of chemical reactions According to nature of reaction Homogeneous In the same phase (aqueous, gaseous,...) Examples: Acid base, oxydo reduction Heterogeneous Involve different phases Examples: Sorption, precipitation / dissolution,... According to speed of reaction Slow reactions Irreversible, modeled using kinetic law L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 5 / 41
8 Classification of chemical reactions According to nature of reaction Homogeneous In the same phase (aqueous, gaseous,...) Examples: Acid base, oxydo reduction Heterogeneous Involve different phases Examples: Sorption, precipitation / dissolution,... According to speed of reaction Slow reactions Irreversible, modeled using kinetic law Fast reactions Reversible, modeled using equilibrium L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 5 / 41
9 Classification of chemical reactions According to nature of reaction Homogeneous In the same phase (aqueous, gaseous,...) Examples: Acid base, oxydo reduction Heterogeneous Involve different phases Examples: Sorption, precipitation / dissolution,... According to speed of reaction Slow reactions Irreversible, modeled using kinetic law Fast reactions Reversible, modeled using equilibrium Depends on relative speed of reactions and transport. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 5 / 41
10 Classification of chemical reactions According to nature of reaction Homogeneous In the same phase (aqueous, gaseous,...) Examples: Acid base, oxydo reduction Heterogeneous Involve different phases Examples: Sorption, precipitation / dissolution,... According to speed of reaction Slow reactions Irreversible, modeled using kinetic law Fast reactions Reversible, modeled using equilibrium Depends on relative speed of reactions and transport. In this talk: Equilibrium reactions, with sorption. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 5 / 41
11 Aqueous reactions: chemical equilibrium Chemical reaction: α 1 C 1 + α 2 C 2 S, in general α i C i 0 i Thermodynamic equilibrium: minimize (change in) Gibb s free energy G = G 0 + RT i α i ln(c i ), R = 8.31 J/K/mol Leads to mass action law ( c α i i = K, K = exp G ) 0 RT i L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 6 / 41
12 Modeling general equilibrium General chemical reactions : N s species, N r reactions N s j=1 ν ij Y j 0, i = 1,..., N r ν ij stoichiometric coefficients. Matrix equation νy = 0 L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 7 / 41
13 Modeling general equilibrium General chemical reactions : N s species, N r reactions N s j=1 ν ij Y j 0, i = 1,..., N r ν ij stoichiometric coefficients. Matrix equation νy = 0 Assumption ν has full rank : Rank ν = N r. Basis for null-space of ν has dimensions N c = N s N r. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 7 / 41
14 Modeling general equilibrium General chemical reactions : N s species, N r reactions N s j=1 ν ij Y j 0, i = 1,..., N r ν ij stoichiometric coefficients. Matrix equation νy = 0 Assumption ν has full rank : Rank ν = N r. Basis for null-space of ν has dimensions N c = N s N r. Partition ν = ( B N ), B R Nr Nr invertible, N R Nc Nr, R = B 1 N ( ) X General solution of νy = 0: Y =, X = RC. C R C Nc, X R Nr. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 7 / 41
15 Example: carbonic acid dissociation Species H 2 O, H 3 O +, OH, H 2 CO 3, HCO 3, CO2 3 Reactions Possible components H 3 O + + OH 2H 2 O, H 2 O + H 2 CO 3 H 3 O + + HCO 3, HCO 3 + H 2O CO H 3 O +. H 2 O, H 3 O +, H 2 CO 3, H 2 O, H 3 O +, HCO 3, H 2 O, H 3 O +, CO 2 3, H 2 O, OH, H 2 CO 3, H 2 O, OH, HCO 3, H 2 O, OH, CO 2 3. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 8 / 41
16 Components: the Morel tableau Species in c: components, in x: secondary species. Rewrite chemical system as (TDB give components, then species) N c j=1 r ij C j X i Components Equ. C 1 C j C Nc constants X 1 R 11 R 1Nc K 1.. X i R i1 R ij R inc K i.. X Ns R 1Ns R NsNc K Ns Total conc. T 1 T i T Nc L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 9 / 41
17 From chemistry to mathematics Each reaction, mass action law ({X i } = activity of X i ) N c {X i } = K i {C j } r ij, i = 1,..., N r j=1 Each component, mass conservation ([X i ] = concentration of X i ) N r T j = [C j ] + r ij [X i ], i=1 j = 1,..., N c L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 10 / 41
18 From chemistry to mathematics Each reaction, mass action law ({X i } = activity of X i ) N c {X i } = K i {C j } r ij, i = 1,..., N r j=1 Each component, mass conservation ([X i ] = concentration of X i ) N r T j = [C j ] + r ij [X i ], j = 1,..., N c i=1 Assume ideal solution [X i ] = {X i } := x (for dissolved species). System of nonlinear algebraic equations log x = R log c + log K T = c + R T x L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 10 / 41
19 Sorption processes Definition Sorption is the accumulation of a fluid on a solid at the fluid solid interface. Main mechanism for exchanges between dissolved species and solid surfaces. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 11 / 41
20 Sorption processes Definition Sorption is the accumulation of a fluid on a solid at the fluid solid interface. Main mechanism for exchanges between dissolved species and solid surfaces. Several possible mechanisms Surface complexation Formation of bond between surface and aqueous species, due to electrostatic interactions. Depends on surface potential. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 11 / 41
21 Sorption processes Definition Sorption is the accumulation of a fluid on a solid at the fluid solid interface. Main mechanism for exchanges between dissolved species and solid surfaces. Several possible mechanisms Surface complexation Formation of bond between surface and aqueous species, due to electrostatic interactions. Depends on surface potential. Ion exchange Ions are exchanged between sorption sites on the surface. Depends on Cationic Exchange Capacity. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 11 / 41
22 Sorption processes Definition Sorption is the accumulation of a fluid on a solid at the fluid solid interface. Main mechanism for exchanges between dissolved species and solid surfaces. Several possible mechanisms Surface complexation Formation of bond between surface and aqueous species, due to electrostatic interactions. Depends on surface potential. Ion exchange Ions are exchanged between sorption sites on the surface. Depends on Cationic Exchange Capacity. Can be modeled as mass action law L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 11 / 41
23 Sorptions isotherms Definition An adsorption isotherm relates F (mol/g) quantity of adsorbed component to its concentration C (mol/l) in the fluid L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 12 / 41
24 Sorptions isotherms Definition An adsorption isotherm relates F (mol/g) quantity of adsorbed component to its concentration C (mol/l) in the fluid Transport of one solute in saturated porous medium, interaction with solid phase : with C t + ρ F t Non-equilibrium df = k(ψ(c) F ) dt Equilibrium F = Ψ(C) + div ( uc D C ) = 0 L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 12 / 41
25 Common equilibrium isotherms Linear Ψ(C) = K d C Langmuir Ψ(C) = κ 1C 1 + κ 2 C Concave near C = 0 Freundlich Ψ(C) = κc p 1 Concave near C = 0, Ψ (0+) = for 0 < p < 1, 0.5 convex near C = 0 for 1 < p,. 0 Analysis for general function Ψ : Knabner, van Duijn Langmuir freundlich Linear L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 13 / 41
26 Examples for linear, Langmuir and Freundlich isotherm Concentration profiles at tmax/3, 2*tmax/3, and tmax Concentration profiles at tmax/3, 2*tmax/3, and tmax Concentration profiles at tmax/3, 2*tmax/3, and tmax Linear isotherm, Pe = Langmuir isotherm 0.9 Freundlich isotherm K L = 6.1/C 0 n F = Pe = Pe = C/C C/C C/C x (cm) x (cm) x (cm) J. A. Cunningham (Texas A&M) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 14 / 41
27 Outline 1 Modeling Chemical phenomena Chemical phenomena in aqueous chemistry Modeling equilibrium systems Sorption models 2 Multispecies equilibrium reactive transport Chemical problem Transport equations The coupled system Coupling algorithms 3 Numerical results Pyrite test case Chromatography Simplified MoMaS benchmark L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 15 / 41
28 Reactions for aqueous/solid system c j aqueous components, s j sorbant components, x i aqueous secondary species, fixed y i secondary species. N c x i S ij c j, j=1 i = 1,..., N x N c N s y j A ij c j + B ij s j, i = 1,..., N y, j=1 j=1 L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 16 / 41
29 Mass action law N c x i = K x i c S ij j, i = 1,..., N x j=1 N c y i = K y i j=1 Use logarithm: linear algebra Mass conservation c A ij j N s j=1 s B ij j, i = 1,..., N y, c + S t x + A T y = T s + B T y = W, T, W : total concentration in components L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 17 / 41
30 The chemical problem System of non-linear equations c + S T x + A T y = T, s + B T y = W, log x = S log c + log K x, log y = A log c + B log s + log K y. Dissolved total: C = c + S T x, Fixed total: F = A T y. Role of chemical model Given totals T (and W, known), split into mobile and immobile total concentrations. C = Φ(T ), F = Ψ(T ) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 18 / 41
31 The chemical problem (2) Take concentration logarithms as main unknowns Nonlinear system H(z) = ( ) T W z = (log c, log s), K = (log K x, log K y ) H(z) = exp(z + S T exp(k + Sz) Jacobian matrix H (z) = diag exp(z) + S T diag(exp(k + Sz)) S L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 19 / 41
32 Numerical solution of nonlinear problem Use globalized Newton s method (line search, trust region). Ion exchange: 6 species, 4 components (vary initial guess) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 20 / 41
33 Transport in a porous medium Diffusion convection equation ω c t D 2 c 2 x + u c x = f for 0 < x < L c = c d at x = 0 D c x = x 0 at x = L c(x, 0) = c 0 (x), 0 < x < L. ω : porosity Assumption D : dispersion coefficient Let L(c) = D 2 c 2 x + u c x. Dispersion tensor independent of species u : Darcy velocity L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 21 / 41
34 Numerical method for transport Space-time finite difference method c n j t x 2 c n+θ = θc n+1 + (1 θ)c n. ω cn+1 j D cn+θ j+1 2cn+θ j + c n+θ j 1 + u cn+θ j c n+θ j 1 x = f n+θ j Implicit scheme Unconditionally stable Upwind scheme (first order in space) θ = 1/2 : Crank Nicolson scheme (2nd order in time) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 22 / 41
35 The coupled system Transport for each species and component x i t y i t + L(x i ) = r x i, = r y i, c j t s j t + L(c j ) = r c j, = r s j, L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 23 / 41
36 The coupled system Transport for each species and component x i t y i t + L(x i ) = r x i, = r y i, c j t s j t + L(c j ) = r c j, = r s j, Eliminate (unknown) reaction rates by using conservation laws: CD equations for totals (T = C + F ) T ic + L(C ic ) = 0, ic = 1,..., N c t T ic ix = Cic ix + F ic ix ic = 1,.., N c and ix = 1,.., N x F ix = Ψ(T ix ) ix = 1,..., N x. Number of transport equations reduced from N x + N y to N c + N s L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 23 / 41
37 Different formulations (1) CC formulation, explicit chemistry (J. Erhel) dc + df + LC = 0 dt dt ( ) C + F H(z) = 0 W F F (z) = 0. + Explicit Jacobian + Chemistry function, no chemical solve Intrusive approach (chemistry not a black box) Precipitation not easy to include L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 24 / 41
38 Different formulations (2) TC formulation, implicit chemistry dt + LC = 0 dt T C F = 0 F Ψ(T ) = 0 + Non-intrusive approach (chemistry as black box) + Precipitation can (probably) be included No explicitjacobian (finite differences) One chemical solve for each function evaluation L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 25 / 41
39 Standard iterative algorithm Block Gauss Seidel method Transport C n+1,k+1 + F n+1,k T n + L(C n+1,k+1 ) = 0, t T n+1,k+1 = C n+1,k+1 + F n+1,k, Chemistry F n+1,k+1 = Ψ(T n+1,k+1 ) Yeh Tripathi (1989), Saaltink et al. (2001), Carrayrou (2001), Dimier, Montarnal et al. (2004). L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 26 / 41
40 DAE approach for CC Coupled system is index 1 DAE (J. Erhel, C. de Dieuleveult) M dy dt + f (y) = 0 L C I 0 I ( C ) y = z, M = 0 0 0, f (y) = C + F H(z) W F F F (z) Use standard DAE software L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 27 / 41
41 Global method for TC C n+1 C n + F n+1 F n + L(C n+1 ) = 0 t t T n+1 = C n+1 + F n+1 F n+1 = Ψ(T n+1 ) C n+1 F n+1 + (I + tl)c n+1 C n F n f T = F n+1 T n+1 C n+1 F n+1 F n+1 Ψ(T n+1 ) = 0. Solve by Newton s method L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 28 / 41
42 Structure of Jacobian matrix (I + tl) 0 I Jacobian : f (C, T, F ) = I I I 0 Ψ (T ) I Ψ (T ) jacobian of chemistry Storage of jacobian matrix is expensive, size of matrix is 3N x N c 3N x N c L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 29 / 41
43 Newton Krylov method Solve the linear system by an iterative method GMRES, TFQMR and BiCGStab require only jacobian matrix by vector products. Can be approximated by finite differences or computed analytically. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 30 / 41
44 Newton Krylov method Solve the linear system by an iterative method GMRES, TFQMR and BiCGStab require only jacobian matrix by vector products. Can be approximated by finite differences or computed analytically. Inexact Newton Approximation of the Newton s direction: f (x k )d + f (x k ) η f (x k ) (0<η<1) Choice of the forcing term η? Keep quadratic convergence (locally) Avoid oversolving the linear system η = γ f (x k ) 2 / f (x k 1 ) 2 (Kelley, Eisenstat and Walker) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 30 / 41
45 Newton Krylov method (2) Computing the jacobian-vector product f f (y + hw) f (y) (y)w h Choice of h? h = 10 7 x w (Kelley). Outstanding issue: preconditioning References Hammond, Valocchi, Lichtner (CMWR, 2002) Knoll, Keyes (JCP, 2004) Mousseau, Knoll (JCP, 2004) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 31 / 41
46 Outline 1 Modeling Chemical phenomena Chemical phenomena in aqueous chemistry Modeling equilibrium systems Sorption models 2 Multispecies equilibrium reactive transport Chemical problem Transport equations The coupled system Coupling algorithms 3 Numerical results Pyrite test case Chromatography Simplified MoMaS benchmark L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 32 / 41
47 Pyrite test case 4 components, 39 aqueuous and 13 fixed species Somewhat artificial without precipitation Use DAE software (Variable time step and order. In Matlab, function ode15s modified to use UMFPACK). C. de Dieuleveult s thesis (INRIA/Andra) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 33 / 41
48 Example: ion exchange Column experiment (Phreeqc ex. 11, Alliances ex. 3) Column contains a solution with 1mmol Na, 0.2mmol K and 1.2mmol NO 3. Inject solution with 1.2mmolCaCl 2. CEC = Chemical system: Ca X CaX 2, log K 1 = 0.8, Na + + X NaX, log K 2 = 0, K + + X KX, log K 3 = 0.7 L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 34 / 41
49 Chromatography: result Concentration at end of column Concentrations (mmol/l) Ca Na Cl K Without diffusion, can be solved semi-analytically (Appelo et al.). Extension to diffusive case? Temps (h) L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 35 / 41
50 Algorithm performance 2000 Comparaison of CPU time 1800 CPU time (second) Global approach Splitting appraoch Iterative Splitting appraoch As mesh is refined Number of Newton iterations remains stable Number of GMRES iterations grows number of nodes For each time step: Block Gauss-Seidel iterations, Newton Krylov 4-7 Newton steps, GMRES steps L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 36 / 41
51 Simplified MoMaS benchmark One dimension, diffusive regime, shorter time period Easy chemistry, Morel tableau X 1 X 2 X 3 X 4 S K C C C C C CS CS Initial: X 2, X 4, inject X 1, X 2, X 3, flush with X 2, X 4. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 37 / 41
52 Simplified benchmark : results Split x = 1/30, t = 2. Split x = 1/50, t = 0.5. Global x = 1/30, t = 2. L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 38 / 41
53 Benchmark: performance Nonlinear residual Nonlinear convergence Residual norm Nonlinear convergence Outer iteration For each time step: 8-9 Newton iterations Function evaluations Number of inner iterations increases for each further Newton iteration function evaluations L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 39 / 41
54 MoMaS benchmark proposal Written by J. Carrayrou (IMFS). International Scientific Committee 1D and 2D geometries 3 levels of chemical difficulty Easy Equilibrium, ion exchange Medium More species, kinetics Hard Precipitation dissolution L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 40 / 41
55 Conclusions Perspectives Formulation of reactive transport within mathematical framework Implementation of Newton Krylov algorithm Preliminary performance tests L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 41 / 41
56 Conclusions Perspectives Formulation of reactive transport within mathematical framework Implementation of Newton Krylov algorithm Preliminary performance tests Areas for further work Numerical analysis of algorithms? Comparison of different formulations, and algorithms Handling of precipitation reactions Handling of per species diffusion coefficient Extension to multi-phase flows L. Amir, M. Kern Coupling transport and chemistry GT Méth. numériques LJLL 41 / 41
Reactive transport in porous media: formulations, non-linear solvers and preconditioners
Reactive transport in porous media: formulations, non-linear solvers and preconditioners Michel Kern with L. Amir, B. Gueslin, A. Taakili Institut National de Recherche en Informatique et Automatique High
More informationCoupling transport with geochemistry Applications to CO2 sequestration
Coupling transport with geochemistry Applications to CO2 sequestration Michel Kern with L. Amir, B. Gueslin, A. Taakili INRIA Paris Rocquencourt 3ème Conférence Internationale de la Société Marocaine de
More informationPreconditioning the Newton-Krylov method for reactive transport
Preconditioning the Newton-Krylov method for reactive transport Michel Kern, L. Amir, A. Taakili INRIA Paris-Rocquencourt Maison de la Simulation MoMaS workshop Reactive Transport Modeling in the Geological
More informationPRECONDITIONING A COUPLED MODEL FOR REACTIVE TRANSPORT IN POROUS MEDIA
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 16, Number 1, Pages 18 48 c 2019 Institute for Scientific Computing and Information PRECONDITIONING A COUPLED MODEL FOR REACTIVE TRANSPORT
More informationThree-dimensional Modelling of Reactive Solutes Transport in Porous Media
151 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 41, 214 Guest Editors: Simonetta Palmas, Michele Mascia, Annalisa Vacca Copyright 214, AIDIC Servizi S.r.l., ISBN 978-88-9568-32-7; ISSN 2283-9216
More informationReactive transport modeling of carbonate diagenesis on unstructured grids
Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität Leoben Dolomitization and geochemical modeling Hollywood may never
More informationTIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA
XVIII International Conference on Water Resources CMWR 2 J. Carrera (Ed) c CIMNE, Barcelona, 2 TIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA Florian Haeberlein and Anthony Michel
More informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
More informationReactive Transport in Porous Media
Reactive Transport in Porous Media Jean-Baptiste Apoung-Kamga, Pascal Have, Jean Houot, Michel Kern, Adrien Semin To cite this version: Jean-Baptiste Apoung-Kamga, Pascal Have, Jean Houot, Michel Kern,
More informationComputational modelling of reactive transport in hydrogeological systems
Water Resources Management III 239 Computational modelling of reactive transport in hydrogeological systems N. J. Kiani, M. K. Patel & C.-H. Lai School of Computing and Mathematical Sciences, University
More informationA posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media
A posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media D. A. Di Pietro, M. Vohraĺık, and S. Yousef Université Montpellier 2 Marseille,
More informationCoupled Thermo-Hydro-Mechanical and Chemical Analysis of Expansive Clay Subjected to Heating and Hydration
Transport in Porous Media (2007) 66:341 372 Springer 2007 DOI 10.1007/s11242-006-0014-z Coupled Thermo-Hydro-Mechanical and Chemical Analysis of Expansive Clay Subjected to Heating and Hydration L. DO
More informationReactive-transport modelling of electrokinetic extraction of heavy metals from marine sediments
EREM63 Reactive-transport modelling of electrokinetic extraction of heavy metals from marine sediments Matteo Masi a, *, Alessio Ceccarini b, Renato Iannelli a a University of Pisa, Department of Energy,
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationChemical Hydrogeology
Physical hydrogeology: study of movement and occurrence of groundwater Chemical hydrogeology: study of chemical constituents in groundwater Chemical Hydrogeology Relevant courses General geochemistry [Donahoe]
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationColloid and colloid-facilitated transport modeling using HPx
Colloid and colloid-facilitated transport modeling using HPx D. Jacques 1, D. Zhou 2, J. Makselon 3, I. Engelhardt 2, and S. Thiele-Bruh 4 1 SCK-CEN, 2 TU Berlin, 3 Forschungszentrum Jülich GmbH, 4 University
More informationContents Preface Introduction Model Concepts
Preface xi 1 Introduction 1 1.1 Environmental Problems and Geochemical Modeling 1 1.1.1 High-Level Radioactive Waste Disposal 1 1.1.2 Mining Related Environmental Issues 4 1.1.3 Landfills 8 1.1.4 Deep
More informationComparison of numerical schemes for multiphase reactive transport
Énergies renouvelables Production éco-responsable Transports innovants Procédés éco-efficients Ressources durables Comparison of numerical schemes for multiphase reactive transport T. Faney, A. Michel,
More informationEvaporation-driven soil salinization
Evaporation-driven soil salinization Vishal Jambhekar 1 Karen Schmid 1, Rainer Helmig 1,Sorin Pop 2 and Nima Shokri 3 1 Department of Hydromechanics and Hydrosystemmodeling, University of Stuttgart 2 Department
More informationELEMENT-BASED FORMULATIONS FOR COUPLED FLOW, TRANSPORT, AND CHEMICAL REACTIONS
ELEMENT-BASED FORMULATIONS FOR COUPLED FLOW, TRANSPORT, AND CHEMICAL REACTIONS A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE
More informationAn Accelerated Block-Parallel Newton Method via Overlapped Partitioning
An Accelerated Block-Parallel Newton Method via Overlapped Partitioning Yurong Chen Lab. of Parallel Computing, Institute of Software, CAS (http://www.rdcps.ac.cn/~ychen/english.htm) Summary. This paper
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationNumerical Methods for Large-Scale Nonlinear Equations
Slide 1 Numerical Methods for Large-Scale Nonlinear Equations Homer Walker MA 512 April 28, 2005 Inexact Newton and Newton Krylov Methods a. Newton-iterative and inexact Newton methods. Slide 2 i. Formulation
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationC. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2
On iterative solvers combined with projected Newton methods for reacting flow problems C. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2 1 Delft University of Technology J.M. Burgerscentrum 2 Delft University
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationDepartment of Computer Science, University of Illinois at Urbana-Champaign
Department of Computer Science, University of Illinois at Urbana-Champaign Probing for Schur Complements and Preconditioning Generalized Saddle-Point Problems Eric de Sturler, sturler@cs.uiuc.edu, http://www-faculty.cs.uiuc.edu/~sturler
More informationElectrical double layer
Electrical double layer Márta Berka és István Bányai, University of Debrecen Dept of Colloid and Environmental Chemistry http://dragon.unideb.hu/~kolloid/ 7. lecture Adsorption of strong electrolytes from
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationThermohaline Flow and Reactive Solute Transport in Fractured Rock
Thermohaline Flow and Reactive Solute Transport in Fractured Rock Thomas GRAF Institute of Fluid Mechanics in Civil Engineering Alberta Geological Survey, Edmonton 2008-2009 Basin Hydrogeologist Sasdf
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationImplicit-explicit exponential integrators
Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators
More informationGEOCHEMISTRY, GROUNDWATER AND POLLUTION,
GEOCHEMISTRY, GROUNDWATER AND POLLUTION, 2 ND EDITION C.A.J. APPELO Hydrochemical Consultant, Amsterdam, the Netherlands D. POSTMA Environment & Resources DTU, Technical University of Denmark, Kgs. Lyngby,
More informationThe Phase Field Method
The Phase Field Method To simulate microstructural evolutions in materials Nele Moelans Group meeting 8 June 2004 Outline Introduction Phase Field equations Phase Field simulations Grain growth Diffusion
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationEstimation of Flow Geometry, Swept Volume, and Surface Area from Tracer Tests
Estimation of Flow Geometry, Swept Volume, and Surface Area from Tracer Tests Paul W. Reimus, Los Alamos National Laboratory G. Michael Shook, Chevron Energy Technology Company 28 th Oil Shale Symposium
More informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationC. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2
s: application to chemical vapor deposition C. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2 1 Delft University of Technology J.M. Burgerscentrum 2 Delft University of Technology Department of Multi Scale Physics
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 informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationELECTROCHEMICAL SYSTEMS
ELECTROCHEMICAL SYSTEMS Third Edition JOHN NEWMAN and KAREN E. THOMAS-ALYEA University of California, Berkeley ELECTROCHEMICAL SOCIETY SERIES WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC PUBLICATION PREFACE
More informationParallel Iterative Methods for Sparse Linear Systems. H. Martin Bücker Lehrstuhl für Hochleistungsrechnen
Parallel Iterative Methods for Sparse Linear Systems Lehrstuhl für Hochleistungsrechnen www.sc.rwth-aachen.de RWTH Aachen Large and Sparse Small and Dense Outline Problem with Direct Methods Iterative
More informationAltered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems
Altered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems Sanjay K Khattri Abstract We present an Altered Jacobian Newton Iterative Method for solving nonlinear elliptic problems Effectiveness
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationReactive Transport in Porous Media
Reactive Transport in Porous Media Daniel M. Tartakovsky Department of Mechanical & Aerospace Engineering University of California, San Diego Alberto Guadagnini, Politecnico di Milano, Italy Peter C. Lichtner,
More informationRosenbrock time integration combined with Krylov subspace enrichment for unsteady flow simulations
Master of Science Thesis Rosenbrock time integration combined with Krylov subspace enrichment for unsteady flow simulations Unsteady aerodynamics David Blom January 11, 2013 Ad Rosenbrock time integration
More informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More informationSpace-Time Domain Decomposition Methods for Transport Problems in Porous Media
1 / 49 Space-Time Domain Decomposition Methods for Transport Problems in Porous Media Thi-Thao-Phuong Hoang, Elyes Ahmed, Jérôme Jaffré, Caroline Japhet, Michel Kern, Jean Roberts INRIA Paris-Rocquencourt
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 informationTRANSPORT 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 information1D Verification Examples
1 Introduction 1D Verification Examples Software verification involves comparing the numerical solution with an analytical solution. The objective of this example is to compare the results from CTRAN/W
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More information7th AIAA/USAF/NASA/ISSMO Symposium on Multidisciplinary Analysis and Optimization September 2 4, 1998/ St. Louis, MO
AIAA 98 4712 Acceleration of Multidisciplinary Analysis Solvers by Inexact Subsystem Simulations Matthias Heinkenschloss Mary Beth Hribar Michael Kokkolaras Rice University Houston, Texas 77005-1892 7th
More informationNumerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
More informationNumerical solution of Least Squares Problems 1/32
Numerical solution of Least Squares Problems 1/32 Linear Least Squares Problems Suppose that we have a matrix A of the size m n and the vector b of the size m 1. The linear least square problem is to find
More information2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 00:1 6 [Version: 2000/07/27 v1.0] Nonlinear Additive Schwarz Preconditioners and Application in Computational Fluid
More informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
More informationGYPSUM PRECIPITATION/DISSOLUTION DURING SEAWATER INTRUSION
18 SWIM. Cartagena 2004, Spain. (Ed. Araguás, Custodio and Manzano). IGME GYPSUM PRECIPITATION/DISSOLUTION DURING SEAWATER INTRUSION N. BOLUDA-BOTELLA, V. GOMIS-YAGÜES, F. RUIZ-BEVIÁ and Mª. D. SAQUETE-FERRÁNDIZ
More informationChapter 7: Anion and molecular retention
I. Anions and molecules of importance in soils Anions of major importance to agricultural soils and soil chemistry are: H 2 PO - 4, HPO 2-4, SO 2-4, HCO - 3, NO - 3, Cl -, F - and OH -. Also, micronutrients
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationDownloaded 08/04/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 38, No. 4, pp. B593 B618 c 2016 Society for Industrial and Applied Mathematics ACTIVE-SET REDUCED-SPACE METHODS WITH NONLINEAR ELIMINATION FOR TWO-PHASE FLOW PROBLEMS IN POROUS
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationInexact Newton Methods and Nonlinear Constrained Optimization
Inexact Newton Methods and Nonlinear Constrained Optimization Frank E. Curtis EPSRC Symposium Capstone Conference Warwick Mathematics Institute July 2, 2009 Outline PDE-Constrained Optimization Newton
More informationHow to simulate in Simulink: DAE, DAE-EKF, MPC & MHE
How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE Tamal Das PhD Candidate IKP, NTNU 24th May 2017 Outline 1 Differential algebraic equations (DAE) 2 Extended Kalman filter (EKF) for DAE 3 Moving horizon
More informationMéthodes de raffinement espace-temps
Méthodes de raffinement espace-temps Caroline Japhet Université Paris 13 & University of Maryland Laurence Halpern Université Paris 13 Jérémie Szeftel Ecole Normale Supérieure, Paris Pascal Omnes CEA,
More informationNumerical simulation of the Gross-Pitaevskii equation by pseudo-spectral and finite element methods comparison of GPS code and FreeFem++
. p.1/13 10 ec. 2014, LJLL, Paris FreeFem++ workshop : BECASIM session Numerical simulation of the Gross-Pitaevskii equation by pseudo-spectral and finite element methods comparison of GPS code and FreeFem++
More informationA Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin
A Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin Jiang M.S. Candidate Joined Fall 2013 1 Main Points Advanced
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 informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationLecture 11: CMSC 878R/AMSC698R. Iterative Methods An introduction. Outline. Inverse, LU decomposition, Cholesky, SVD, etc.
Lecture 11: CMSC 878R/AMSC698R Iterative Methods An introduction Outline Direct Solution of Linear Systems Inverse, LU decomposition, Cholesky, SVD, etc. Iterative methods for linear systems Why? Matrix
More informationPreconditioned inverse iteration and shift-invert Arnoldi method
Preconditioned inverse iteration and shift-invert Arnoldi method Melina Freitag Department of Mathematical Sciences University of Bath CSC Seminar Max-Planck-Institute for Dynamics of Complex Technical
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
More informationMultiscale hp-adaptive quiescent reactors
Multiscale hp-adaptive quiescent reactors C.E.Michoski Institute for Computational Engineering and Sciences (ICES) University of Texas (UT) at Austin Group Meeting Nov. 7, 2011 Outline 1 2 Equation (or
More informationHigh Performance Nonlinear Solvers
What is a nonlinear system? High Performance Nonlinear Solvers Michael McCourt Division Argonne National Laboratory IIT Meshfree Seminar September 19, 2011 Every nonlinear system of equations can be described
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationAn Inexact Newton Method for Nonlinear Constrained Optimization
An Inexact Newton Method for Nonlinear Constrained Optimization Frank E. Curtis Numerical Analysis Seminar, January 23, 2009 Outline Motivation and background Algorithm development and theoretical results
More informationModeling the next battery generation: Lithium-sulfur and lithium-air cells
Modeling the next battery generation: Lithium-sulfur and lithium-air cells D. N. Fronczek, T. Danner, B. Horstmann, Wolfgang G. Bessler German Aerospace Center (DLR) University Stuttgart (ITW) Helmholtz
More informationIntroduction to Heat and Mass Transfer. Week 7
Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms
More informationCOUPLING METHODOLOGY WITHIN THE SOFTWARE PLATFORM ALLIANCES
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2005 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) CIMNE, Barcelona, 2005 COUPLING METHODOLOGY WITHIN
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Catalysis of environmental reactions Dr. Zifei Liu Catalysis and catalysts Catalysis is the increase in the rate of a chemical reaction due to the participation
More informationDynamic Optimization in Air Quality Modeling
Dynamic Optimization in Air Quality Modeling A. Caboussat Department of Mathematics, University of Houston Houston, Texas 204-3008 caboussat@math.uh.edu http://aero.math.uh.edu Project supported by the
More informationGroundwater chemistry
Read: Ch. 3, sections 1, 2, 3, 5, 7, 9; Ch. 7, sections 2, 3 PART 14 Groundwater chemistry Introduction Matter present in water can be divided into three categories: (1) Suspended solids (finest among
More informationANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow
ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow Michel Kern Inria Paris Rocquencourt Maison de la Simulation C2S@Exa Days, Inria Paris Centre, Novembre 2016 M. Kern
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationAnswers to Exercises Computational Fluid Dynamics
Answers to Exercises Computational Fluid Dynamics Exercise - Artificial diffusion upwind computations.9 k=. exact.8 k=..7 k=..6 k=..5.4.3.2...2.3.4.5.6.7.8.9 x For k =.,. and., and for N = 2, the discrete
More informationSimplifying Chemical Kinetics
Texas A&M University September 1, 2010 Chemically Reactive Flows: In gas-phase chemical reaction systems, the governing porcesses such as flow, chemical reactions and molecular transport, occur at time
More informationIterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems Luca Bergamaschi e-mail: berga@dmsa.unipd.it - http://www.dmsa.unipd.it/ berga Department of Mathematical Methods and Models for Scientific Applications University
More informationModeling and numerical approximation of multi-component anisothermal flows in porous media
Modeling and numerical approximation of multi-component anisothermal flows in porous media M. Amara, D. Capatina, L. Lizaik,, P. Terpolilli Laboratory of Applied Mathematics, CNRS UMR 5142, University
More informationDepartment of Civil Engineering-I.I.T. Delhi CVL722 1st Semester HW Set2. Adsorption
Department of Civil Engineering-I.I.T. Delhi CVL722 1st Semester 2016-17 HW Set2 Adsorption Q1: For the following information, determine Langmuir and Freundlich model constants? Also plot Q of these models
More informationDepartment of Ocean Engineering, Indian Institute of Technology-Madras, Chennai , India. *Corresponding author.
J. Earth Syst. Sci. (2018) 127:53 c Indian Academy of Sciences https://doi.org/10.1007/s12040-018-0950-3 Interaction of dissolution, sorption and biodegradation on transport of BTEX in a saturated groundwater
More informationSplitting Iteration Methods for Positive Definite Linear Systems
Splitting Iteration Methods for Positive Definite Linear Systems Zhong-Zhi Bai a State Key Lab. of Sci./Engrg. Computing Inst. of Comput. Math. & Sci./Engrg. Computing Academy of Mathematics and System
More information