Water Artificial Circulation For Eutrophication Control
|
|
- Ursula Barton
- 6 years ago
- Views:
Transcription
1 Water Artificial Circulation For Eutrophication Control Francisco Javier Fernández Centro Universitario de la Defensa, Escuela Naval Militar, Marn, Spain Joint work with Lino J. Álvarez-Vázquez and Aurea Martínez Depto. Matemática Aplicada II, Universidade de Vigo, E.I. Telecomunicación Vigo, Spain
2 Table of contents. The environmental problem Mathematical formulation Mathematical analysis Numerical resolution Numerical tests
3 The environmental problem
4 The environmental problem Eutrophication (environmental problem) Eutrophication is one of the most important problems of large masses of water and it is caused by high levels of pollutants that reach the waters. This pollutants come mainly from human activities and can cause an excessive phytoplankton growth that lead to undesirable effects like algal blooms. 1 Artificial circulation (main idea) Its is a management technique for oxygenating eutrophic water bodies subject to quality problems, such as loss of oxygen, sediment accumulation and algal blooms. It disrupts stratification and minimizes the development of stagnant zones that may be subject to water quality problems (low levels of oxygen). In our case we are interested in increase the dissolved oxygen content in the bottom layers.
5 The environmental problem 1 Artificial circulation (how it works) A flow pump takes water from the well aerated upper layers by means of a collector and injecting it into the poorly oxygenated bottom layers, setting up a circulation pattern that prevents stratification. Oxygen-poor water from the bottom is circulated to the surface, where oxygenation from the atmosphere and photosynthesis can occur. Objective We want to optimize the volume of pumped water in order to minimize the energy cost and ensure that the concentration level of dissolved oxygen in the bottom layer is over a threshold.
6 Mathematical formulation
7 Mathematical formulation 1 N C 1 C 2 C T 2 T 1 S T 3 The control variables C 3 C 4 T 4 N N The physical domain R 3 (for instance a lake). C, control domain. N CT number of pairs collector-injector. = Γ S Γ N Γ C Γ T, where: Γ S, surface of the lake, Γ N, bottom and walls, Γ C, collectors, Γ C = N CT k=1 C k, Γ T, injectors, Γ T = N CT k=1 T k. g(t) = (g 1 (t),..., g N CT (t)), where g k (t) is the volume of water displaced by pump k. The state equations The hydrodynamic model: modified Navier-Stokes equations (Smagorinsky). ( ) v/ t + vv (2ν + 2ν tur [ɛ(v) : ɛ(v)] 1/2 )ɛ(v) + p = F in (0, T ), v = 0 in (0, T ), v = 0 on (Γ S Γ N ) (0, T ), v = +g k (t)/µ(c k ) n on C k (0, T ), for k = 1,..., N CT, v = g k (t)/µ(t k ) n on T k (0, T ), for k = 1,..., N CT, v(0) = v 0 in,
8 Mathematical formulation The eutrophication model: Michaelis-Menten kinetics. u 1 generic nutrient (usually, nitrogen and/or phosphorus), u 2 phytoplankton u 3 zooplankton, u 4 organic detritus, and u 5 dissolved oxygen. Interactions between five above species can be modelled by the following system: u i / t + v u i (µ i u i ) = A i (u) + f i in (0, T ), u i n = 0 on (Γ S Γ N Γ C ) (0, T ), u i = C k u i dγ/µ(c k ) on T k (0, T ), k = 1,..., N CT, u i (0) = u i 0 in, for i = 1,..., 5, where the reaction term A = (A i ) : (0, T ) [R +] 5 R 5 is: ) u C nc (L(x, 1 t) K N + u 1 u2 K r u 2 + C nc K rd Θ θ(x,t) 20 u 4 ( ) u 1 L(x, t) K N + u 1 u2 K r u 2 K mf u 2 u 2 K z K F + u 2 u3 A(x, t, u) = u 2 C fz K z K F + u 2 u3 K mz u 3 K mf u 2 + K mz u 3 K rd Θ θ(x,t) 20 u 4 ) u C oc (L(x, 1 t) K N + u 1 u2 K r u 2 C oc K rd Θ θ(x,t) 20 u 4
9 Mathematical formulation The optimal control problem where: (P) mín{j(g) : g U ad and 1 u(t) dx K C, t [0, T ]}, µ( C ) C U ad = {g U : 0 < c 1 g k (t) c 2, t [0, T ], k = 1,..., N CT }, with U a suitable functional space and the constants c 1, c 2 related to technological constraints on the pumps. We will have to redefine it in order to perform the mathematical analysis of the problem. J(g) is the cost functional given by: J(g) = 1 N CT T g k (t) 2 dt, 2 k=1 0 representing the energy saving related to pumps operation. K C R 5 is the rectangle: K C = [λ m 1, λm 1 ] [λm 5, λm 5 ], with λ m k and λm k the minimal and maximal mean concentrations allowed for the species in C. (v, u) are the solutions of the state system associated to control g. We have a optimal control problem with nonlinear state constraints.
10 Mathematical analysis
11 Mathematical analysis The space for the control variables The state equations U ad = { g H s 1/2 ( (0, T )) L 2 (0, T ; H 1 ( )) : c 1 g(x, t) 0, a.e. (x, t) T k (0, T ), k = 1,..., N CT, +0 g(x, t) c 2, a.e. (x, t) C k (0, T ), k = 1,..., N CT, g(x, t) = 0 a.e. (x, t) (Γ N Γ S ) (0, T ), k = 1,..., N CT, T k C k g(t) dγ = 0, a.e. t (0, T ), k = 1,..., N CT, g H s 1/2 ( (0,T )) L 2 (0,T ; H1 ( )) c 3 } The hydrodynamic model. If s 7/2, there exists an unique solution: v W 1,,2 (0, T ; W, [L 2 ()] 3 ) C([0, T ]; W), v(0) = v 0 a.e. in. v = g n, a.e. (x, t) (0, T ). v verifies the variational formulation: v t z dx + vv z dx + β(ɛ(v))ɛ(v) : ɛ(z) dx = F z dx, a.e. t (0, T ), z V. where: W = {v [W 1,3 ()] 3 : v = 0 and v (ΓN Γ S ) = 0}, V = {v [W 1,3 ()] 3 : v = 0 and v = 0} W 1,p,q (I ; X, X ) = { u L p (I ; X ) : u/ t L q (I ; X ) }.
12 Mathematical analysis The eutrophication model. If v W 1,,2 (0, T ; W, [L 2 ()] 3 ) C([0, T ]; W) is the solution of the hydrodynamic model, then there exists a unique solution for the eutrophication model (we also need compatibility conditions for oefficients and data): u W 1,2,2 (0, T ; [H 1 ()] 5, [H 1 ()] 5 ), u i (0) = u i 0 a.e. x, u i T k = 1 u i dγ a.e. (x, t) T k (0, T ), µ(c k ) C k u i verifies the variational formulation: u i t η dx + v u i η dx + µ i u i η dx = f i η dx + A i (u)η dx, η X, a.e. t (0, T ). The control problem Existence. Assuming that the coefficients and data of the problem satisfy the necessary conditions for the existence of a solution of the state systems, there exists, at least, a solution of the optimal control problem (P). Optimality conditions. We can prove the existence of first order necessary conditions that must be verified by the optimal control.
13 Numerical resolution
14 Numerical resolution Space-Time discretization A partition {t 0, t 1,..., t N } of the time interval [0, T ] for a time step t = T /N, where t n = n t for n = 0, 1,..., N, will be considered. We will use the Characteristics Method for the discretization of the material derivatives. We will consider a regular triangular mesh τ h of a polygonal domain and we will use the Finite Element Method for the space discretization, with: W h /V h (P 1b FEM space) for the water velocity (u), M h (P 1 FEM space) for the water pressure (p), Zh/X h (P 1 FEM space) for the species concentration of the eutrophication model (u). Software development The space and time discretizations will be implemented in the scientific software Freefem++. The numerical resolution of the constrained optimization problem (obtained from the space-time discretization of the continuous control problem) has been developed with the help of the interior point algorithm IPOPT, interfaced with Freefem++
15 Numerical resolution The discretized control problem After the space-time discretization of the continuous problem (P), we obtain the following constrained nonlinear minimization problem: where: (Q) mín{j(g) : g U ad }, Discretized control: We will consider the following set of admissible discrete controls: U ad = {g R N N CT : c 1 g k,n c 2, k = 1,..., N CT, n = 1,..., N}, with c 1 and c 2 are bounds related to pumps technical characteristics. Discretized cost functional: We consider the following discretized cost function: J(g) = 1 N N CT (g k,n ) 2. 2 n=1 k=1 Discretized state constraints: We define the real-valued mapping: G = (G 1,..., G N ) : U ad R 5N g G(g), where, for each n = 1,..., N, G n 1 (g) = u n+1 dx, µ( C ) C with u n+1 Z h the solution of the discretized eutrophication model. Thus: U ad = {g U ad : G n (g) K C, n = 1,..., N}.
16 Numerical resolution Discretized state equations: Given an admissible control g = (g 1,1, g 2,1,..., g N CT,1,..., g 1,N, g 2,N,..., g N CT,N ) U ad, }{{}}{{} g 1 g N The hydrodynamic model: For n = 0, 1,..., N 1, (v n+1, p n+1 ) W h M h, with v n+1 T k = g k,n+1 /µ(t k ) n, v n+1 C k = g k,n+1 /µ(c k ) n, k = 1,..., N CT, is the solution of: α v n+1 z dx + 2ν ɛ(v n+1 ) : ɛ(z) dx +2ν tur [ɛ(v n+1 ) : ɛ(v n+1 )] 1/2 ɛ(v n+1 ) : ɛ(z) dx p n+1 z dx v n+1 q dx λ p n+1 q dx = α (v n X n ) z dx + F n+1 z dx, z V h, q M h. The eutrophication model: For n = 0,..., N, u n+1 Z h, with u n+1 T k = u n+1 dγ/µ(c k ), k = 1,..., N CT, C k is the solution of: α u n+1 η dx + Λ µ u n+1 : η dx = α (u n X n ) η dx + f n+1 η dx + A n+1 (u n+1 ) η dx, η X h.
17 Numerical resolution Remark. It should be noted that we employ a different number of time steps in both models: N for the hydrodynamic model and N + 1 for the species model. This shift is motivated by the following dependency scheme v 0 u 0 g 1 v 1 u 1 g 2 v 2 u 2 G 1 (g) = G 1 (g 1 ) g 3 v 3 u 3 G 2 (g) = G 2 (g 1, g 2 ).... g N v N u N G N 1 (g) = G N 1 (g 1, g 2,..., g N 1 ) u N+1 G N (g) = G N (g 1, g 2,..., g N 1, g N ) The Jacobian associated to state constraints is a block lower triangular matrix: δ g1 G 1 (g) δ g1 G 2 (g) δ g2 G 2 (g) δ gg(g) = δ g1 G N 1 (g) δ g2 G N 1 (g)... δ gn 1 G N 1 (g) 0 δ g1 G N (g) δ g2 G N (g)... δ gn 1 G N (g) δ gn G N (g).
18 Numerical resolution Computing the Jacobian associated to state constraints In order to use the IPOPT algorithm in the numerical resolution of the problem (Q), we need a function that evaluates the Jacobian associated to the state constraints. We can use the linearized state equations or the adjoint state equations. The choice of one method or another depends on the relationship between the size of the control and the constraints. If we want to use the linearized equations, we must solve N CT N times this equations. If we employ the adjoint state equations, we must solve 5 N times this equations. We have that if N CT < 5 it is more convenient to use the linearized equations, whereas if N CT > 5, it is more advantageous to employ the adjoint state. In the example presented in the next section, we have considered N CT = 4, therefore, we will use the linearized equations for computing the Jacobian of the discretized state constraints: δg n 1 (g)(δg) = δu n+1 dx, n = 1,..., N, µ( C ) C where δv 0 = 0, δu 0 = 0 and, given a direction δg, the linearized state is obtained by:
19 Numerical resolution The linearized hydrodynamic model: For n = 0,..., N 1, (δv n+1 h, δp n+1 h ) W h M h, with δv n+1 k,n+1 δg T k = µ(t k ) n, k,n+1 δg δvn+1 Ck = µ(c k ) n, k = 1,..., N CT, is the solution of: α δv n+1 z dx + 2ν ɛ(δv n+1 ) : ɛ(z) dx +2ν tur [ɛ(v n+1 ) : ɛ(v n+1 )] 1/2 ɛ(δv n+1 ) : ɛ(v n+1 ) ɛ(v n+1 ) : ɛ(z) dx +2ν tur [ɛ(v n+1 ) : ɛ(v n+1 )] 1/2 ɛ(δv n+1 ) : ɛ(z) dx δp n+1 z dx δv n+1 q dx λ δp n+1 q dx = α (δv n X n ) z dx ( v n X n )δv n z dx, v V h, q M h. The linearized eutrophication model: For n = 0,..., N, the element δu n+1 Z h, with δu n+1 T k = 1 δu n+1 dγ, k = 1,..., N CT, µ(c k ) C k is the solution of: α δu n+1 η dx + Λ µ u n+1 : η dx = α (δu n X n ) η dx ( u n X n )δv n η dx + δ ua n+1 (u n+1 )(δu n+1 ) z dx, η X h.
20 Numerical tests
21 Numerical tests The following data have been considered for testing the methodology proposed (2D example): Computation domain Time discretization (4 scenarios) Test 1: 12h, t = 3600s. Test 2: 12h, t = 1800s. C 1 C 2 S C 3 C 4 Test 3: 24h, t = 3600s. Test 4: 48h, t = 3600s. N N Space discretization = [0, 14] [0, 16] (meters). N CT = 4 pairs collector-injector. C T 2 T 1 T 3 T 4 N C = [0, 14] [0, 2]. We will consider a regular mesh of 3075 vertices.
22 Numerical tests Test 1 (12h, t = 3600s) Test 2 (12h, t = 1800s)
23 Numerical tests Test 1 (12h, t = 3600s) Test 2 (12h, t = 1800s) Dissolved Oxygen, time step 12 IsoValue Dissolved Oxygen, time step 25 IsoValue Vec Value 0 Vec Value e e e e e
24 Numerical tests Test 3 (24h, t = 3600s) Test 4 (48h, t = 3600s)
25 Numerical tests Test 3 (24h, t = 3600s) Test 4 (48h, t = 3600s) Dissolved Oxygen, time step 25 IsoValue Dissolved Oxygen, time step 49 IsoValue Vec Value 0 Vec Value e e e e e e e e e e e e e e e e e
26 Thank you for your attention. Any Questions?
CONTROLLING EUTROPHICATION IN A MOVING DOMAIN
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationUrban heat islands: An optimal control approach. Departamento de Matemática Aplicada II Universidad de Vigo Spain
Urban heat islands: An optimal control approach L.J. Alvarez-Vázquez, F.J.Fernández, N. García-Chan, A. Martínez, M.E. Vázquez-Méndez. Departamento de Matemática Aplicada II Universidad de Vigo Spain The
More informationJournal of Mathematical Analysis and Applications. Mathematical analysis of a three-dimensional eutrophication model
J Math Anal Appl 349 9) 35 55 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications wwwelseviercom/locate/jmaa Mathematical analysis of a three-dimensional eutrophication
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationStability Analysis of Plankton Ecosystem Model. Affected by Oxygen Deficit
Applied Mathematical Sciences Vol 9 2015 no 81 4043-4052 HIKARI Ltd wwwm-hikaricom http://dxdoiorg/1012988/ams201553255 Stability Analysis of Plankton Ecosystem Model Affected by Oxygen Deficit Yuriska
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 informationA Finite Element Approximation of a Coupled Model for Two Turbulent Fluids
A Finite Element Approximation of a Coupled Model for Two Turbulent Fluids Driss Yakoubi: GIREF, Université Laval, Québec. Canada joint works with Tomas Chacon Rebollo: Laboratoire Jacques-Louis Lions,
More informationMeasurements of In Situ Pick-up Rate of Nutrients on Riverbed
1th International Conference on Integrated Diffuse Pollution Management (IWA DIPCON 8). Research Center for Environmental and Hazardous Substance Management (EHSM), Khon Kaen University, Thailand; 5-9
More informationMathematical Modeling and Analysis of Eutrophication of Water Bodies Caused by Nutrients
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 511 524 Mathematical Modeling and Analysis of Eutrophication of Water Bodies Caused by Nutrients A. K. Misra Department of Mathematics,
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
More informationOn pollution and R&D-based growth in a trade model between countries concerned and unconcerned about climate change
On pollution and R&D-based growth in a trade model between countries concerned and unconcerned about climate change Francisco Cabo 1 Guiomar Martín-Herrán 1 M.Pilar Martínez-García 2 1 IMUVa, Dpt. Economía
More informationOptimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints
Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints
More informationModelling tools for assessing lagoon water quality and biological resources
Modelling tools for assessing lagoon water quality and biological resources Arturas Razinkovas-Baziukas Petras Zemlys 2011-03 Spatial heterogeneity and morphological complexity requires not traditional
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) LMS
More informationLecture Pure Twist
Lecture 4-2003 Pure Twist pure twist around center of rotation D => neither axial (σ) nor bending forces (Mx, My) act on section; as previously, D is fixed, but (for now) arbitrary point. as before: a)
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationDesign of optimal RF pulses for NMR as a discrete-valued control problem
Design of optimal RF pulses for NMR as a discrete-valued control problem Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Carla Tameling (Göttingen) and Benedikt Wirth
More informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationFinite element analysis of the Reynolds lubrication equation with cavitation
Rakenteiden Mekaniikka (Journal of Structural Mechanics) Vol. 45, No 4, 22, pp. 88 2 Finite element analysis of the Reynolds lubrication equation with cavitation Arto Sorsimo, Mika Juntunen, Rolf Stenberg
More informationPartitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers
Copyright c 2007 ICCES ICCES, vol.2, no.1, pp.21-27, 2007 Partitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers L. Rodríguez-Tembleque 1, J.A.
More informationConstrained optimization. Unconstrained optimization. One-dimensional. Multi-dimensional. Newton with equality constraints. Active-set method.
Optimization Unconstrained optimization One-dimensional Multi-dimensional Newton s method Basic Newton Gauss- Newton Quasi- Newton Descent methods Gradient descent Conjugate gradient Constrained optimization
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
More informationAlgae and Dissolved Oxygen Dynamics of Landa Lake and the Upper Spring Run
Algae and Dissolved Oxygen Dynamics of Landa Lake and the Upper Spring Run Why study algae and dissolved oxygen dynamics of Landa Lake and the Upper Spring Run? During low-flow conditions, extensive algal
More informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationTheoretical and numerical results for a chemo-repulsion model with quadratic production
Theoretical and numerical results for a chemo-repulsion model with quadratic production F. Guillén-Gonzalez, M. A. Rodríguez-Bellido & and D. A. Rueda-Gómez Dpto. Ecuaciones Diferenciales y Análisis Numérico
More informationGlobal Carleman inequalities and theoretical and numerical control results for systems governed by PDEs
Global Carleman inequalities and theoretical and numerical control results for systems governed by PDEs Enrique FERNÁNDEZ-CARA Dpto. E.D.A.N. - Univ. of Sevilla joint work with A. MÜNCH Lab. Mathématiques,
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationSome inverse problems with applications in Elastography: theoretical and numerical aspects
Some inverse problems with applications in Elastography: theoretical and numerical aspects Anna DOUBOVA Dpto. E.D.A.N. - Univ. of Sevilla joint work with E. Fernández-Cara (Universidad de Sevilla) J. Rocha
More informationPeristaltic Pump. Introduction. Model Definition
Peristaltic Pump Introduction In a peristaltic pump, rotating rollers are squeezing a flexible tube. As the pushed down rollers move along the tube, the fluid in the tube follows the motion. The main advantage
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationControl Theory Approach to Aero Shape Optimization
Control Theory Approach to Aero Shape Optimization. Overview. Finite Volume Method 3. Formulation of Optimal Design Problem 4. Continuous Adjoint Approach 5. Discrete Adjoint Approach Overview of Adjoint
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More information3D DLM/FD METHODS FOR SIMULATING THE MOTION OF SPHERES IN BOUNDED SHEAR FLOWS OF OLDROYD-B FLUIDS
3D DLM/FD METHODS FOR SIMULATING THE MOTION OF SPHERES IN BOUNDED SHEAR FLOWS OF OLDROYD-B FLUIDS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationNon-Differentiable Embedding of Lagrangian structures
Non-Differentiable Embedding of Lagrangian structures Isabelle Greff Joint work with J. Cresson Université de Pau et des Pays de l Adour CNAM, Paris, April, 22nd 2010 Position of the problem 1. Example
More informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More informationarxiv: v2 [math.ap] 4 Mar 2016
Existence of a weak solution to a fluid-elastic structure interaction problem with the Navier slip boundary condition Boris Muha Sunčica Čanić arxiv:155.446v [math.ap] 4 Mar 16 Abstract We study a nonlinear,
More information14. *14.8 CASTIGLIANO S THEOREM
*14.8 CASTIGLIANO S THEOREM Consider a body of arbitrary shape subjected to a series of n forces P 1, P 2, P n. Since external work done by forces is equal to internal strain energy stored in body, by
More informationFINITE ELEMENT METHODS IN LOCAL ACTIVE CONTROL OF SOUND
FINITE ELEMENT METHODS IN LOCAL ACTIVE CONTROL OF SOUND ALFREDO BERMÚDEZ, PABLO GAMALLO, AND RODOLFO RODRÍGUEZ Abstract. The active control of sound is analyzed in the framework of the mathematical theory
More informationPAPER 331 HYDRODYNAMIC STABILITY
MATHEMATICAL TRIPOS Part III Thursday, 6 May, 016 1:30 pm to 4:30 pm PAPER 331 HYDRODYNAMIC STABILITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationThermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space
1/29 Thermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space Jungho Park Department of Mathematics New York Institute of Technology SIAM
More informationNumerical Methods for PDE-Constrained Optimization
Numerical Methods for PDE-Constrained Optimization Richard H. Byrd 1 Frank E. Curtis 2 Jorge Nocedal 2 1 University of Colorado at Boulder 2 Northwestern University Courant Institute of Mathematical Sciences,
More informationAGITATION AND AERATION
AGITATION AND AERATION Although in many aerobic cultures, gas sparging provides the method for both mixing and aeration - it is important that these two aspects of fermenter design be considered separately.
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationLinear Transformations: Standard Matrix
Linear Transformations: Standard Matrix Linear Algebra Josh Engwer TTU November 5 Josh Engwer (TTU) Linear Transformations: Standard Matrix November 5 / 9 PART I PART I: STANDARD MATRIX (THE EASY CASE)
More informationVorticity and Dynamics
Vorticity and Dynamics In Navier-Stokes equation Nonlinear term ω u the Lamb vector is related to the nonlinear term u 2 (u ) u = + ω u 2 Sort of Coriolis force in a rotation frame Viscous term ν u = ν
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationOPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL
OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL EDUARDO CASAS, ROLAND HERZOG, AND GERD WACHSMUTH Abstract. Semilinear elliptic optimal control
More informationLake Pyykösjärvi water quality model Timo Huttula and Ninni Liukko Freshwater Centre, Finnish Environment Institute, Finland
Lake Pyykösjärvi water quality model Timo Huttula and Ninni Liukko Freshwater Centre, Finnish Environment Institute, Finland timo.huttula@ymparisto.fi Abstract: Lake Pyykösjärvi waters are eutrophied with
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 informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationOn Some Variational Optimization Problems in Classical Fluids and Superfluids
On Some Variational Optimization Problems in Classical Fluids and Superfluids Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas
More informationWork with a partner. Read Section page 60 in Section 2.4, and discuss answers to questions C F. Discuss your responses with the class. Any Questions?
Work with a partner Read Section page 60 in Section 2.4, and discuss answers to questions C F. Discuss your responses with the class. LokenTimer2.swf Any Questions? Title: Jun 7 8:32 AM (1 of 25) Title:
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationTurbulence Modeling of Air Flow in the Heat Accumulator Layer
PIERS ONLINE, VOL. 2, NO. 6, 2006 662 Turbulence Modeling of Air Flow in the Heat Accumulator Layer I. Behunek and P. Fiala Department of Theoretical and Experimental Electrical Engineering Faculty of
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 informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationChapter 2. Formulation of Finite Element Method by Variational Principle
Chapter 2 Formulation of Finite Element Method by Variational Principle The Concept of Variation of FUNCTIONALS Variation Principle: Is to keep the DIFFERENCE between a REAL situation and an APPROXIMATE
More informationNumerical control and inverse problems and Carleman estimates
Numerical control and inverse problems and Carleman estimates Enrique FERNÁNDEZ-CARA Dpto. E.D.A.N. - Univ. of Sevilla from some joint works with A. MÜNCH Lab. Mathématiques, Univ. Blaise Pascal, C-F 2,
More informationInverse problems in lithospheric flexure and viscoelasticity
Inverse problems in lithospheric flexure and viscoelasticity Axel Osses in coll with M. de Buhan (CNRS, France), E. Contreras (Geophysics Department, Chile), B. Palacios (DIM, Chile) DIM - Departamento
More informationApplications of a Local Energy Method to Systems of PDE s Involving Free Boundaries
Applications of a Local Energy Method to Systems of PDE s Involving Free Boundaries Gonzalo Galiano Abstract. We present a method of analysis for free boundary problems which is based on local energy estimates.
More informationAn Eulerian-Lagrangian approach for fluid-structure coupled systems using X-FEM Antoine Legay / CNAM-Paris
An Eulerian-Lagrangian approach for fluid-structure coupled systems using X-FEM Antoine Legay / CNAM-Paris Collaborations: A. Tralli, P. Gaudenzi, Università di Roma La Sapienza, Italia 8th USNCCM, 25-27th
More informationc 2002 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 4 No. pp. 507 5 c 00 Society for Industrial and Applied Mathematics WEAK SECOND ORDER CONDITIONS FOR STOCHASTIC RUNGE KUTTA METHODS A. TOCINO AND J. VIGO-AGUIAR Abstract. A general
More informationFigure 1: Surface waves
4 Surface Waves on Liquids 1 4 Surface Waves on Liquids 4.1 Introduction We consider waves on the surface of liquids, e.g. waves on the sea or a lake or a river. These can be generated by the wind, by
More informationStochastic Modeling of Chemical Reactions
Stochastic Modeling of Chemical Reactions Eric Haseltine Department of Chemical Engineering University of Wisconsin-Madison TWMCC 27 September Texas-Wisconsin Modeling and Control Consortium 22TWMCC Outline
More informationOPTIMAL TRANSPORTATION PLANS AND CONVERGENCE IN DISTRIBUTION
OPTIMAL TRANSPORTATION PLANS AND CONVERGENCE IN DISTRIBUTION J.A. Cuesta-Albertos 1, C. Matrán 2 and A. Tuero-Díaz 1 1 Departamento de Matemáticas, Estadística y Computación. Universidad de Cantabria.
More informationChapter 3 Steady-State, ne- mens onal C on uction
Chapter 3 Steady-State, One-Dimensional i Conduction 3.1 The Plane Wall 3.1.1 Temperature Distribution For one-dimensional, steady-state conduction in a plane wall with no heat generation, the differential
More informationChanges in phytoplankton and zooplankton production in the Nordic and the Arctic Seas under a warmer climatic regime.
Changes in phytoplankton and zooplankton production in the Nordic and the Arctic Seas under a warmer climatic regime. Dag Slagstad, Ingrid Ellingsen and Morten Alver, Trondheim, Norway 2 nd ESSAS Open
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationAssignment 2 Continuous Matter 4335/8191(B)
Assignment 2 Continuous Matter 4335/8191(B Tyler Shendruk October 29, 2012 1 Problem 1 Given the stream function ψ = x x 3 y/2, determine if mass is conserved the flow is irrotational and the velocity
More informationError Analyses of Stabilized Pseudo-Derivative MLS Methods
Department of Mathematical Sciences Revision: April 11, 2015 Error Analyses of Stabilized Pseudo-Derivative MLS Methods Don French, M. Osorio and J. Clack Moving least square (MLS) methods are a popular
More informationLevel Set Tumor Growth Model
Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth
More informationMultiple Choice Review Problems
Multiple Choice Review Problems 1. (NC) Which graph best represents the position of a particle, st ( ), as a function of time, if the particle's velocity is negative and the particle's acceleration is
More informationCHAPTER 5 WARM UPS. Mrs. Hilliard
CHAPTER 5 WARM UPS Mrs. Hilliard CHAPTER 5 VOCABULARY 1. Photosynthesis 2. Cellular respiration 3. Producer 4. Consumer 5. Decomposer 6. Food chain 7. Food web 8. Trophic level 9. Carbon cycle 10. Nitrogen-fixing
More informationWATER QUALITY MANAGEMENT AND CONTROL OF POLLUTION IN LAKES
WATER QUALITY MANAGEMENT AND CONTROL OF POLLUTION IN LAKES FACULTY OF MATHEMATICS XI Modelling Week UCM Giulia Del Veccia Ignacio Ferrero Melchor Víctor Manuel González Vicente Karina Haydeé Molina Martínez
More informationLecture Kollbruner Section 5.2 Characteristics of Thin Walled Sections and.. Kollbruner Section 5.3 Bending without Twist
Lecture 3-003 Kollbruner Section 5. Characteristics of Thin Walled Sections and.. Kollbruner Section 5.3 Bending without Twist thin walled => (cross section shape arbitrary and thickness can vary) axial
More informationBayesian inverse problems with Laplacian noise
Bayesian inverse problems with Laplacian noise Remo Kretschmann Faculty of Mathematics, University of Duisburg-Essen Applied Inverse Problems 2017, M27 Hangzhou, 1 June 2017 1 / 33 Outline 1 Inverse heat
More informationIntroduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationA C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani
A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes Marco Verani MOX, Department of Mathematics, Politecnico di Milano Joint work with: P. F. Antonietti (MOX Politecnico di
More information1. The accumulated net change function or area-so-far function
Name: Section: Names of collaborators: Main Points: 1. The accumulated net change function ( area-so-far function) 2. Connection to antiderivative functions: the Fundamental Theorem of Calculus 3. Evaluating
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationCH 24. Electric Potential
CH 24 Electric Potential [SHIVOK SP212] January 8, 2016 I. Electric Potential Energy A. Experimentally, physicists and engineers discovered that the electric force is conservative and thus has an associated
More informationTurbulent spectra generated by singularities
Turbulent spectra generated by singularities p. Turbulent spectra generated by singularities E.A. Kuznetsov Landau Institute for Theoretical Physics, Moscow Russia In collaboration with: V. Naulin A.H.
More informationStochastic calculus of variations on the diffeomorphisms group
Title Stochastic calculus of variations on the diffeomorphisms group Warwick, April 2012 Ana Bela Cruzeiro Dep. Mathematics IST (TUL) Grupo de Física-Matemática Univ. Lisboa ( ) 1 / 19 The diffeomorphisms
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationFaraday Tensor & Maxwell Spinor (Part I)
February 2015 Volume 6 Issue 2 pp. 88-97 Faraday Tensor & Maxwell Spinor (Part I) A. Hernández-Galeana #, R. López-Vázquez #, J. López-Bonilla * & G. R. Pérez-Teruel & 88 Article # Escuela Superior de
More informationGradient estimates and global existence of smooth solutions to a system of reaction-diffusion equations with cross-diffusion
Gradient estimates and global existence of smooth solutions to a system of reaction-diffusion equations with cross-diffusion Tuoc V. Phan University of Tennessee - Knoxville, TN Workshop in nonlinear PDES
More informationReynolds number scaling of inertial particle statistics in turbulent channel flows
Reynolds number scaling of inertial particle statistics in turbulent channel flows Matteo Bernardini Dipartimento di Ingegneria Meccanica e Aerospaziale Università di Roma La Sapienza Paolo Orlandi s 70th
More informationBoundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids
Boundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids Wendy Kress and Jonas Nilsson Department of Scientific Computing Uppsala University Box S-75 4 Uppsala Sweden
More informationNX: Limnological and Engineering Analysis of a Polluted Urban Lake. Prelude to Environmental Management of Onondaga Lake, New York.
Steven W. Effler Editor Limnological and Engineering Analysis of a Polluted Urban Lake Prelude to Environmental Management of Onondaga Lake, New York With 557 Illustrations NX: Springer Series Preface
More informationGeneral Curvilinear Ocean Model (GCOM): Enabling Thermodynamics
General Curvilinear Ocean Model (GCOM): Enabling Thermodynamics M. Abouali, C. Torres, R. Walls, G. Larrazabal, M. Stramska, D. Decchis, and J.E. Castillo AP0901 09 General Curvilinear Ocean Model (GCOM):
More informationSolving Partial Differential Equations Numerically. Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah
Solving Partial Dierential Equations Numerically Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah Overview What are partial dierential equations? How do we solve them? (Example)
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
More informationA PRELIMINARY STUDY OF MAGNETIC PROPERTIES IN SEDIMENTS FROM THE ANLLÓNS RIVERBED, AND THEIR POSSIBLE ASSOCIATION TO ANTHROPOGENIC POLLUTION
A PRELIMINARY STUDY OF MAGNETIC PROPERTIES IN SEDIMENTS FROM THE ANLLÓNS RIVERBED, AND THEIR POSSIBLE ASSOCIATION TO ANTHROPOGENIC POLLUTION Vincenzo Costanzo-Álvarez 1, Diego López-Rodríguez 2, Rosa Devesa-Rey
More information