Inverse problems in tumor growth modelling
|
|
- Philip Turner
- 5 years ago
- Views:
Transcription
1 MC2 Damiano Lombardi Angelo Iollo Thierry Colin Olivier Saut Marseille,
2 Outline 1) Introduction: phenomenological model for tumor growth models 2) Models based on mixture theory: a simplified model 3) Inverse problems: formulation and numerical procedure 4) Realistic Applications: Prognosis tool for lung tumors 5) Conclusions: analysis of results and perspectives
3 Introduction Tumor growth is an intrinsic multi-scale phenomenon: Microscopic scale: description of genetic mutations and molecular pathways; sub cellular level. Mesoscopic scale: behavior of a cell in its social context; biochemical signals and transport phenomena; cellular level. Macroscopic scale: interactions between cellular phenotypes in a tissue, mechanical propagation; tissue and organ level. It is impossible to model tumor growth using first principles
4 Introduction The role of mathematical modeling: Understand quantitatively the complex relationships involved Use models in realistic applications: prognosis tools, treatment efficacy evaluations, therapy impact Different type of models: Individual based: discrete models, agent based models and cellular automata. PDE models: continuous type models, based on reaction diffusion equations, transport equations.
5 Continuous type models Population dynamics: mass conservation P i t + (J) =Γ(P j,c k ) N(P j,c k ) Variation of the cell density for the species i: J is the flux Production and Reduction, depending on other species and nutrients The choice of J determines the model nature: J depending on the gradient of P: reaction-diffusion model J depending on P and a vector field v: transport model P i t + (v ip i )=Γ(P j,c k ) N(P j,c k ),
6 Continuous type models The meaning of v: v(x; t) = X(x; t) t Evolution of the Lagrangian Coordinate v is unknown: conservation of momentum (P i v i ) t + (P i v i v i )= σ i + j U i,j (P i,p j ; v i, v j ), Variation of momentum in a given point, at time t Divergence of Stress Tensor: Interaction between species: The stress tensor is linked to v through constitutive equations
7 Continuous type models Nutrients evolution: a reaction diffusion model C k t (D(P i ) C k )= ( ) S j (P j ) M j (P j ) j C k Time variation of the nutrients concentration; Diffusion: in general it is non-uniform and non-isotropic; Production and Consumption. The nutrient equation is coupled with the population dynamics Coefficients depend on tissue properties
8 Simplified 2D Darcy model 2 cellular species: P t Q t + (vp ) = (2γ 1)P, + (vq) = (1 γ)p. P is the proliferating cell density Q is the dead cell density Healthy tissue is denoted by S, and P+Q+S=1: ( ) v = γp. v = k(p, Q) Π. { } Equation for the divergence Darcy-Law for the velocity The concentration of oxygen: C t (D(Y ) C) =α S Ω 1 Y αpc λc 1 Y dω
9 Simplified 2D Darcy model Equation for the porosity: k = k 1 +(k 2 k 1 )(P + Q) Equation for the diffusivity: D = D max K(P + Q) Hypoxia threshold: γ = 1 + tanh(r(c C hyp)) 2, Non-uniform, isotropic, linear Non-uniform, isotropic, linear Regularization of unit step Hypoxia function realizes the coupling between oxygen diffusion and population dynamics.
10 The need of inverse problems: Inverse Problems; The model is a system of coupled non-linear parametric PDE. Microscopic scales are not modeled: their effects is lumped into parameters that are not known a priori In order to apply models, parameters have to be determined via inverse problems, the data coming from medical imagery
11 Inverse problems: Inverse Problems; Observables have to be defined: Y = P + Q With Scan different tumor phenotypes may not be distinguished: the sum is the only observable. The time evolution of the observable is derived: Y t + (vy )=γp The other equations become: v = γp k(y ) v = k(y ) v C t (D(Y ) C) = αpc λc,
12 Inverse problems: Inverse Problems; Unknowns: k 2 /k 1,D max, K, α, λ,c hyp v, P, C Parameters Fields that are not observables From a medical standpoint P and C are of the highest importance Normally only two-three scans are available: the problem is greatly under-determined Infinite solutions are possible and fit the images in an optimal sense.
13 Inverse problems: Inverse Problems; In order to find one solution a regularization approach is defined. It consists in choosing a priori a particular functional space, and looking for the solution that belongs to this space. The base of the space is constructed by means of semiempirical eigenfunctions, defined using the results of numerical simulations. It is equivalent to add information imposing a functional constrain to the unknown fields.
14 Inverse problems: Inverse Problems; P = a (P ) i φ (P ) i,i=1,..., N P, C = a (C) i φ (C) i,i=1,..., N C, v = a (v) i φ (v) i,i=1,..., N v, γp = a (γp ) i φ (γp ) i,i=1,..., N γp, Eigenfunctions are extracted from a database of simulations by means of Proper Orthogonal Decomposition (POD). POD given the dimension N maximizes the energy representation of a given signal (P) φ i = j bi jp j λ 1/2 i. λ and b are the eigenvalues and the eigenvectors of the Autocorrelation matrix A (P ) ij =< P k (t h ),P m (t n ) > h, n =1,..., N; k, m =1,..., M
15 Inverse problems: Inverse Problems; The POD expansions are substituted into the equations written for the observable Y Ẏ + a v i (Y φv i )=aγp j φ (γp ) j a v i (φv i )=aγp j φ γp j k(y )a v i φ v i = k(y ) a v i φ v i a C i φc i a C i (K(Y ) φ C i )= αa p i ac i φ C i φ P i λa C i φ C i 2a γp i φ γp i = 1 + tanh(r(a C i φc i C hyp )) Unknowns: k 2 /k 1,D max, K, α, λ,c hyp a P i,a C i,a v i,a γp i Parameters Expansion Coefficients
16 Inverse problems: procedure Direct equations written at certain time instants Are data sufficient to perform a complete identification? yes Inverse Problems; no We build a data base varying all the parameters We need to add information We extract an eigenfunction basis We identify the system using a non-linear Newton solver
17 Realistic cases: Tumour in lung: 2D scan (Institut Bergonié) first scan: second scan: third scan: Starting from the first and the second scan we try to reconstruct the third one
18 Realistic cases: 1) Identification: all the parameters and the unknown quantities are identified in such a way that we are able to simulate the entire tumour growth; 2) Simulation: we take the first image as initial condition and we simulate the tumour (125 simulations in the database); we observe the differences between the third image and the tumour simulated at the same time. The simulation is rather realistic: 1) Tumour grows up in the same way; 2) It moves towards the boundary; 3) The dynamics seem to be well captured.
19 Realistic cases: Simulations of the second and third scans second and third scans
20 Realistic cases: Volume curve: - Simulation (continuous line); - scan (circle) Errors: Simulations reach the volume of the third scan a little bit after; Volume T(days) Why? 1) We have only 2D partial information; 2) The model is approximated, not explicitly designed for a lung; 3) We do not have considered angiogenesis.
21 Realistic cases: Another case, from the other side of the biological spectrum: slow, quasi-stable lung nodules The amount of information is higher, we try to solve the same problem
22 Realistic cases: Simuation databases: 125 and 800 simulations. The offline stage is parallelized. Given the first two images, we tried to reconstruct all the others. Initial condition is the first image, for all the simulations. The simulation is realistic since a slow tumor growth is represented. The shape is not so good. Irregularity in the boundary are introduced mainly from initial condition
23 Realistic cases: Simulations of the second and third scans second and third scans
24 Realistic cases: 0.09 Volume curve: what about non-logistic behavior? 0.08 Predictions Vol Data months There is an intrinsic non-linearity in the model
25 Conclusions: An efficient numerical procedure have been set up: The numerical cost is concentrated in the construction of the database and the eigenmodes, but this stage is parallelized Prognose is realiable on the time scale of medical exams Errors and open questions: Insufficient representation of tissue anisotropies: the shape of the tumor is not yet well captured. High dimensional database imply more difficult optimization with local methods: study the sensitivity to initial guess. The reconstruction of the initial conditions is a critical stage for several reasons. Regularization process and time scale of the observed phenomena.
26 Perspectives: It is difficult to improve the accuracy, but we have to improve the reliability. Study the sensitivity to initial conditions, the variance of the results. Identification in several cases: a population approach. Models: Model for lungs: tissue anisotropies and particular geometry Angiogenesis and therapy model Inverse problems: General database construction, 3D cases Optimal trasport theory to approximate evolution of observables using semiempirical eigenfunctions
27 Thank you for your attention!
Numerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationTwo-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationIntroduction to some topics in Mathematical Oncology
Introduction to some topics in Mathematical Oncology Franco Flandoli, University of Pisa y, Finance and Physics, Berlin 2014 The field received considerable attentions in the past 10 years One of the plenary
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 information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationModelling Drug Delivery and Tumour Growth
Modelling Drug Delivery and Tumour Growth Matthew Hubbard and Christian Groh, Pamela Jones, Brian Sleeman, Steven Smye, Christopher Twelves (Leeds) Paul Loadman, Nagaraja Periasamy, Roger Phillips (Bradford)
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationInvestigating the role of tortuosity in the Kozeny-Carman equation
Investigating the role of tortuosity in the Kozeny-Carman equation Rebecca Allen, Shuyu Sun King Abdullah University of Science and Technology rebecca.allen@kaust.edu.sa, shuyu.sun@kaust.edu.sa Sept 30,
More informationLecture 6. Four postulates of quantum mechanics. The eigenvalue equation. Momentum and energy operators. Dirac delta function. Expectation values
Lecture 6 Four postulates of quantum mechanics The eigenvalue equation Momentum and energy operators Dirac delta function Expectation values Objectives Learn about eigenvalue equations and operators. Learn
More informationLinear vector spaces and subspaces.
Math 2051 W2008 Margo Kondratieva Week 1 Linear vector spaces and subspaces. Section 1.1 The notion of a linear vector space. For the purpose of these notes we regard (m 1)-matrices as m-dimensional vectors,
More informationEfficient Observation of Random Phenomena
Lecture 9 Efficient Observation of Random Phenomena Tokyo Polytechnic University The 1st Century Center of Excellence Program Yukio Tamura POD Proper Orthogonal Decomposition Stochastic Representation
More informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationSimulating Solid Tumor Growth Using Multigrid Algorithms
Simulating Solid Tumor Growth Using Multigrid Algorithms Asia Wyatt Applied Mathematics, Statistics, and Scientific Computation Program Advisor: Doron Levy Department of Mathematics/CSCAMM Abstract In
More informationImprovement of Reduced Order Modeling based on Proper Orthogonal Decomposition
ICCFD5, Seoul, Korea, July 7-11, 28 p. 1 Improvement of Reduced Order Modeling based on Proper Orthogonal Decomposition Michel Bergmann, Charles-Henri Bruneau & Angelo Iollo Michel.Bergmann@inria.fr http://www.math.u-bordeaux.fr/
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationMacroscopic quasi-stationary distribution and microscopic particle systems
Macroscopic quasi-stationary distribution and microscopic particle systems Matthieu Jonckheere, UBA-CONICET, BCAM visiting fellow Coauthors: A. Asselah, P. Ferrari, P. Groisman, J. Martinez, S. Saglietti.
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationJuan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a
Juan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a Universidad de Buenos Aires, Fac. Ing., IGPUBA, ARGENTINA b Department of Mathematics, Purdue University, USA c Universidad
More informationOn a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model
On a Data Assimilation Method coupling, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model 2016 SIAM Conference on Uncertainty Quantification Basile Marchand 1, Ludovic
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationIntroduction to stochastic multiscale modelling of tumour growth
Introduction to stochastic multiscale modelling of tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 1 CIMPA, Santiago de Cuba, June 2016
More informationMaximum variance formulation
12.1. Principal Component Analysis 561 Figure 12.2 Principal component analysis seeks a space of lower dimensionality, known as the principal subspace and denoted by the magenta line, such that the orthogonal
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationNUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING
NUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING Prof. Joško Ožbolt 1 Josipa Bošnjak 1, Goran Periškić 1, Akanshu Sharma 2 1 Institute of Construction Materials,
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationSolving a non-linear partial differential equation for the simulation of tumour oxygenation
Solving a non-linear partial differential equation for the simulation of tumour oxygenation Julian Köllermeier, Lisa Kusch, Thorsten Lajewski MathCCES, RWTH Aachen Lunch Talk, 26.05.2011 J. Köllermeier,
More informationA PDE approach to studying evolutionary and spatial dynamics in cancer cell populations
A PDE approach to studying evolutionary and spatial dynamics in cancer cell populations Tommaso Lorenzi LJLL, 30th June 2017 T. Lorenzi (University of St Andrews) June 30th 2 / 43 Synopsis Object of study
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
More informationAN ITERATIVE PROPER ORTHOGONAL DECOMPOSITION (I-POD) TECHNIQUE WITH APPLICATION TO THE FILTERING OF PARTIAL DIFFERENTIAL EQUATIONS
AAS 2-622 AN ITERATIVE PROPER ORTHOGONAL DECOMPOSITION (I-POD) TECHNIQUE WITH APPLICATION TO THE FILTERING OF PARTIAL DIFFERENTIAL EQUATIONS Dan Yu, Suman Chakravorty INTRODUCTION In this paper, we consider
More informationApplications of Diagonalization
Applications of Diagonalization Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw Apr 2, 200 The notes cover applications of matrix diagonalization Boas 3.2. Quadratic curves Consider the quadratic curve, 5x 2 4xy
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationPrincipal Component Analysis
Principal Component Analysis Laurenz Wiskott Institute for Theoretical Biology Humboldt-University Berlin Invalidenstraße 43 D-10115 Berlin, Germany 11 March 2004 1 Intuition Problem Statement Experimental
More information1 Principal Components Analysis
Lecture 3 and 4 Sept. 18 and Sept.20-2006 Data Visualization STAT 442 / 890, CM 462 Lecture: Ali Ghodsi 1 Principal Components Analysis Principal components analysis (PCA) is a very popular technique for
More informationREDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION
REDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION Lindaura Maria Steffens Dara Liandra Lanznaster lindaura.steffens@udesc.br daraliandra@gmail.com Santa Catarina State University CESFI, Av. Central,
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
More informationUncorrelated Multilinear Principal Component Analysis through Successive Variance Maximization
Uncorrelated Multilinear Principal Component Analysis through Successive Variance Maximization Haiping Lu 1 K. N. Plataniotis 1 A. N. Venetsanopoulos 1,2 1 Department of Electrical & Computer Engineering,
More informationFundamental equations of relativistic fluid dynamics
CHAPTER VI Fundamental equations of relativistic fluid dynamics When the energy density becomes large as may happen for instance in compact astrophysical objects, in the early Universe, or in high-energy
More informationSTA G. Conformal Field Theory in Momentum space. Kostas Skenderis Southampton Theory Astrophysics and Gravity research centre.
Southampton Theory Astrophysics and Gravity research centre STA G Research Centre Oxford University 3 March 2015 Outline 1 Introduction 2 3 4 5 Introduction Conformal invariance imposes strong constraints
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationProper Orthogonal Decomposition
Proper Orthogonal Decomposition Kameswararao Anupindi School of Mechanical Engineering Purdue University October 15, 2010 Kameswararao Anupindi (Purdue University) ME611, Principles of Turbulence October
More informationComputational Modelling of Mechanics and Transport in Growing Tissue
Computational Modelling of Mechanics and Transport in Growing Tissue H. Narayanan, K. Garikipati, E. M. Arruda & K. Grosh University of Michigan Eighth U.S. National Congress on Computational Mechanics
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationThermodynamics for fluid flow in porous structures
Communications to SIMAI Congress, ISSN 1827-9015, Vol. 1 (2006) DOI: 10.1685/CSC06105 Thermodynamics for fluid flow in porous structures M.E. Malaspina University of Messina, Department of Mathematics
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationMaximal Entropy for Reconstruction of Back Projection Images
Maximal Entropy for Reconstruction of Back Projection Images Tryphon Georgiou Department of Electrical and Computer Engineering University of Minnesota Minneapolis, MN 5545 Peter J Olver Department of
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationSolving a non-linear partial differential equation for the simulation of tumour oxygenation
Solving a non-linear partial differential equation for the simulation of tumour oxygenation Julian Köllermeier, Lisa Kusch, Thorsten Lajewski MathCCES, RWTH Aachen Talk at Karolinska Institute, 19.09.2011
More informationOptimal control in diffuse interface models of tumor growth
Optimal control in diffuse interface models of tumor growth E. Rocca Università degli Studi di Pavia 37ème Colloque de la Société Francophone de Biologie Théorique, Poitiers joint work with Harald Garcke
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationA model order reduction technique for speeding up computational homogenisation
A model order reduction technique for speeding up computational homogenisation Olivier Goury, Pierre Kerfriden, Wing Kam Liu, Stéphane Bordas Cardiff University Outline Introduction Heterogeneous materials
More information2.4 Eigenvalue problems
2.4 Eigenvalue problems Associated with the boundary problem (2.1) (Poisson eq.), we call λ an eigenvalue if Lu = λu (2.36) for a nonzero function u C 2 0 ((0, 1)). Recall Lu = u. Then u is called an eigenfunction.
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationContinuum Modeling of Transportation Networks with Differential Equations
with Differential Equations King Abdullah University of Science and Technology Thuwal, KSA Examples of transportation networks The Silk Road Examples of transportation networks Painting by Latifa Echakhch
More informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationUnsupervised learning: beyond simple clustering and PCA
Unsupervised learning: beyond simple clustering and PCA Liza Rebrova Self organizing maps (SOM) Goal: approximate data points in R p by a low-dimensional manifold Unlike PCA, the manifold does not have
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationNon-equilibrium mixtures of gases: modeling and computation
Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present
More informationCovariance and Correlation Matrix
Covariance and Correlation Matrix Given sample {x n } N 1, where x Rd, x n = x 1n x 2n. x dn sample mean x = 1 N N n=1 x n, and entries of sample mean are x i = 1 N N n=1 x in sample covariance matrix
More informationTensor Visualization. CSC 7443: Scientific Information Visualization
Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its
More informationSolving Nonhomogeneous PDEs (Eigenfunction Expansions)
Chapter 12 Solving Nonhomogeneous PDEs (Eigenfunction Expansions) 12.1 Goal We know how to solve diffusion problems for which both the PDE and the s are homogeneous using the separation of variables method.
More informationPrincipal component analysis
Principal component analysis Angela Montanari 1 Introduction Principal component analysis (PCA) is one of the most popular multivariate statistical methods. It was first introduced by Pearson (1901) and
More informationLecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1
L1.P1 Lecture #1 Review Postulates of quantum mechanics (1-3) Postulate 1 The state of a system at any instant of time may be represented by a wave function which is continuous and differentiable. Specifically,
More informationP The Entropy Trajectory: A Perspective to Classify Complex Systems. Tomoaki SUZUDO Japan Atomic Energy Research Institute, JAERI
P 0 0 8 The Entropy Trajectory: A Perspective to Classify Complex Systems Tomoaki SUZUDO Japan Atomic Energy Research Institute, JAERI What are complex systems? Having macroscopic properties which are
More information(U c. t)/b (U t)/b
DYNAMICAL MODELING OF THE LARGE-SCALE MOTION OF A PLANAR TURBULENT JET USING POD MODES. S. Gordeyev 1 and F. O. Thomas 1 University of Notre Dame, Notre Dame, USA University of Notre Dame, Notre Dame,
More informationStochastic optimization - how to improve computational efficiency?
Stochastic optimization - how to improve computational efficiency? Christian Bucher Center of Mechanics and Structural Dynamics Vienna University of Technology & DYNARDO GmbH, Vienna Presentation at Czech
More informationFast and Accurate HARDI and its Application to Neurological Diagnosis
Fast and Accurate HARDI and its Application to Neurological Diagnosis Dr. Oleg Michailovich Department of Electrical and Computer Engineering University of Waterloo June 21, 2011 Outline 1 Diffusion imaging
More informationDimension Reduction Techniques. Presented by Jie (Jerry) Yu
Dimension Reduction Techniques Presented by Jie (Jerry) Yu Outline Problem Modeling Review of PCA and MDS Isomap Local Linear Embedding (LLE) Charting Background Advances in data collection and storage
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationSingular Value Decomposition and Principal Component Analysis (PCA) I
Singular Value Decomposition and Principal Component Analysis (PCA) I Prof Ned Wingreen MOL 40/50 Microarray review Data per array: 0000 genes, I (green) i,i (red) i 000 000+ data points! The expression
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationLaplacian transport towards irregular interfaces: the mathematics
Laplacian transport towards irregular interfaces: the mathematics Denis S GREBENKOV Laboratoire de Physique de la Matière Condensée CNRS Ecole Polytechnique, France Email: denisgrebenkov@polytechniqueedu
More informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
More informationOn Computing Reactive Flows Using Slow Manifolds
On Computing Reactive Flows Using Slow Manifolds by S. Paolucci 1, J. M. Powers 2, and S. Singh 3 Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, Indiana 46556-5637,
More informationDiscriminative Direction for Kernel Classifiers
Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering
More informationFluctuating Hydrodynamics for Transport in Neutral and Charged Fluid Mixtures
Fluctuating Hydrodynamics for Transport in Neutral and Charged Fluid Mixtures Alejandro L. Garcia San Jose State University Lawrence Berkeley Nat. Lab. Hydrodynamic Fluctuations in Soft Matter Simulations
More informationDetailed Program of the Workshop
Detailed Program of the Workshop Inverse Problems: Theoretical and Numerical Aspects Laboratoire de Mathématiques de Reims, December 17-19 2018 December 17 14h30 Opening of the Workhop 15h00-16h00 Mourad
More informationDEIM-BASED PGD FOR PARAMETRIC NONLINEAR MODEL ORDER REDUCTION
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) DEIM-BASED PGD FOR PARAMETRIC NONLINEAR MODEL ORDER REDUCTION JOSE
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationLecture: Computational Systems Biology Universität des Saarlandes, SS Introduction. Dr. Jürgen Pahle
Lecture: Computational Systems Biology Universität des Saarlandes, SS 2012 01 Introduction Dr. Jürgen Pahle 24.4.2012 Who am I? Dr. Jürgen Pahle 2009-2012 Manchester Interdisciplinary Biocentre, The University
More informationV15 Flux Balance Analysis Extreme Pathways
V15 Flux Balance Analysis Extreme Pathways Stoichiometric matrix S: m n matrix with stochiometries of the n reactions as columns and participations of m metabolites as rows. The stochiometric matrix is
More informationStability of Functionally Graded Plate under In-Plane Time-Dependent Compression
Mechanics and Mechanical Engineering Vol. 7, No. 2 (2004) 5 12 c Technical University of Lodz Stability of Functionally Graded Plate under In-Plane Time-Dependent Compression Andrzej TYLIKOWSKI Warsaw
More informationMini-project in scientific computing
Mini-project in scientific computing Eran Treister Computer Science Department, Ben-Gurion University of the Negev, Israel. March 7, 2018 1 / 30 Scientific computing Involves the solution of large computational
More informationAn introduction to PDE-constrained optimization
An introduction to PDE-constrained optimization Wolfgang Bangerth Department of Mathematics Texas A&M University 1 Overview Why partial differential equations? Why optimization? Examples of PDE optimization
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More informationPart 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA
Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA Review of Fundamentals displacement-strain relation stress-strain relation balance of momentum (deformation) (constitutive equation) (Newton's Law)
More informationTwo-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix
Two-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix René Vidal Stefano Soatto Shankar Sastry Department of EECS, UC Berkeley Department of Computer Sciences, UCLA 30 Cory Hall,
More informationStructures in the turbulent wake of an Ahmed body: LES vs Experiments using POD and some ideas about control
LES vs Experiments using POD and some ideas about control Stéphanie PELLERIN and Bérengère PODVIN Laboratoire d Informatique pour la Mécanique et les Sciences de l Ingénieur Context Université Paris Sud,
More informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationChapter 0. Preliminaries. 0.1 Things you should already know
Chapter 0 Preliminaries These notes cover the course MATH45061 (Continuum Mechanics) and are intended to supplement the lectures. The course does not follow any particular text, so you do not need to buy
More informationLocal Anisotropy In Globally Isotropic Granular Packings. Kamran Karimi Craig E Maloney
Local Anisotropy In Globally Isotropic Granular Packings Kamran Karimi Craig E Maloney Granular Materials 2 A Granular Material Is A Conglomeration Of Discrete Solid, Natural Macroscopic Particles Characterized
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More information