Application of Newton Krylov methods to theMoMaS benchmark problem1,2 2Laila AMIR & Michel KERN1 ITASCA Consultants, France2 INRIA Rocquencourt, FranceInternational workshopon Modelling Reactive Transport In Porous MediaStrasbourg, FranceJanuary 23, 2008Laila AMIR & Michel KERN Reactive transport 1 / 24Outline1 Motivations2 Chemistry and transportChemical lawsConvection diffusion3 Coupled modelsMathematical formulationAlgorithm for global approach4 Numerical results5 Conclusion and perspectivesLaila AMIR & Michel KERN Reactive transport 2 / 24Outline1 Motivations2 Chemistry and transportChemical lawsConvection diffusion3 Coupled modelsMathematical formulationAlgorithm for global approach4 Numerical results5 Conclusion and perspectivesLaila AMIR & Michel KERN Reactive transport 3 / 24MotivationsO.S methods⇒* Errors caused by separation of operators* slow convergence* small time stepDSA methods⇒* store Jacobian* solving for the individualconcentrations of the componentsGlobal approach with NKM⇒* keep transport and chemistry modules separate* Jacobian free methodLaila AMIR & Michel KERN Reactive transport 4 / 24Outline1 Motivations2 Chemistry and transportChemical lawsConvection diffusion3 Coupled modelsMathematical formulationAlgorithm for global approach4 Numerical results5 Conclusion and perspectivesLaila AMIR & Michel KERN Reactive transport 5 / 24T TDissolved total : C = c+S x, Fixed total : F = A y.C =Φ(T), F ...