VALIDATION OF A POROELASTIC MODEL
|
|
- Kristian Lindsey
- 5 years ago
- Views:
Transcription
1 VALIDATION OF A POROELASTIC MODEL Tobias Ansali 21/06/2012 Abstract This stuies aim, is to verify from a mathematical point of view the valiity of a moel for the infusion of a ug insie a cancer tissue. Avery important parameter has been optimise in the escription of the phiscal fenomenon the hyaulic conuctivity K. The optimal value has been etermine with Lagrange's approach. The function K was consiere as the optimal value that better represent sperimental ata. The results show how the optimise function has a reasonable tenency from a phisycal point of view an furthermore has a singular tenency to the one obtaine in the moel. The moel In this reasearch we hypothize the soli tumor to have a spherical form.the changes (eformations) that the tumor unergoes are to be consiere innitesimal, for that reason the eformations are governe by Hooke's law (elastic el). We stuy the transfer of the therapeuic agent insie the tumor consiering a porelastic meium carachterize by hyaulic conuctivity k of the tissue (given the relationship between inamic viscosity an permeability of the meium) an by Lame's coecient G an λ. The relationship between the therapeutical agent an bloo vessels is governe by starling's law: L p conuttività vascolare [ S V cm mmhg s ] sup erci e vasc olare p er unità i volume[ cm2 cm 3 ] p e pressione vascolare eettiva[mmhg]. p pressione interstiz iale[mmhg]. Ωcan act both as a well or as a source base on the ierence of pressure insie or outsie the bloo vessel. the tumor is by nature strongly heterogeneous, we only consier it in its raial irection, which is expresse by the hyaulic conuctivity of the tumor K. another stuy has observe that conuctivity k is strongly inuence by the eformation of the tumor. The farmaceutical agent is introuce in the center of the tumor, creating a small raius (a) cavity, which its imensions can be compare to the tip of the neele. The general strength acting on the tumor woul be: T tensor of eective stress T σ pi (2) σ tensor of contact stress p interstizial ui pressure (IFP) consiering that the tension eformation of the tumor is governe by hooke's law, throught the equation of the bon we get: T pi + λ( u)i + 2G[ 1 2 (( u) + ( u)t )] (3) where: Ω L p S V (p e p) (1) where ueformation of soli λ 2ν 1 2ν G 1
2 ove ν is the Poisson's coecient assuming we foun ourselves in stationary conitions, transforming the equation in cilinical coorinates, an an remembering that all variables are hypotetically only expresse in accorance of the raial r coorinate, we obtain: (2G + λ) (u + 2u r ) p (4) To be able to n the istribution of the pressure an eformation we nee another equation. We have to consier the conservation of the mass: q Ω (5) Please note that q has the irection an imention of the velocity, But truthfully is not the actual velocity insie the pores, but it refers to the volumetric carrying for unit area. The secon meaning of55 represent, as alreay sai, the relationship between the vascular net an the pharmaceutical agent uring the evolution of this last one. A ui that evolves insie a pore is escribe by Darcy's law: q K p (6) Taking into consieration the 1 an 5 becomes: ( K p) L p S(p e p) (7) Taking into consieration the 1 r 2 (r2 K p ) L S p V (p e p) (8) Regaring the hyaulic conuctivity of the tumor it has been chosen to consier it as epening from the eformation with a semi-empiric exponential law: u M[α K K 0 e +(1 α) u r ] (9) The equations are close by the bounary conitions. p 0 r R p p inf r a u + 2ν u 1 ν u + 2ν 1 ν r 0 r R u r 0 r a where a an R are respectivly the insie an outsie raiuous after the eformation. Linearization of equations The linearization of the equation is mae through a perturbative analysis of the moel which functions in accorance to a characteristic parameter of the problem. p p p inf pe ; T T 2G+λ ; r r R ; u u r ; K K K 0 ; an the equations become: (u + 2u r ) δ p (10) (r2 K p ) r2 γ 2 (p p e) (11) K u M[α e +(1 α) u r ] (12) having γ 2 Lp V R2 e δ p inf p e 2G+λ Its aroun the aimensional value of δ that the linearization is execute an stopping at the rst orer we have: K 0 S u u 0 + u 1δ K K 0 + K 1 δ p p 0 + p 1δ Please note that we are taking into account the hypothesis of small movements an therefore we can express the exponential as the evelopment in series of mclaurin stopping at the rst orer 2
3 until a R R an from 1 to R. p(a) + p(a) (a a) K 1 + M[α u + (1 α)u r ] Now lets move on to express the conitions on the bounaries of points a R an R Rassuming we are moving linearly from point a R For example p(a ) p inf with a a u(a). We can now write our equations in orerδ 0 : c.c. ( u 0 + 2u 0 r ) 0 p 0 ) r 2 γ 2 (p 0 p e) (r 2 K 0 K0 1 + M[α u 0 + (1 α) u 0 r ] p 0 0 r 1 p 0 0 r a R u 0 + 2ν u 0 1 ν r 0 r 1 u 0 + 2ν u 0 1 ν r 0 r a R It can also be emonstrate thatu 0 0 K 0 1. Even p 0 has an analytical answer such as: p 0 p e + A r eγr + B r e γr with A an B known from bounary conitions At the orer δ 1 the equations are: ( u 1 + 2u 1 r ) p 0 c.c. p 0 ) r 2 γ 2 p 1 (r 2 K 0 K1 1 + M[α u 1 + (1 α) u 1 r ] p 1 u p 0 1 r 1 p 1 u p 0 1 r a R u 1 + 2ν u 1 1 ν r 0 r 1 u 1 + 2ν u 1 1 ν r 0 r a R The carrying that goes through the tumor in the non linear moel has value (aimensional): Q 4πr 2 q In the linear case it woul be where Q 4πr 2 (q 0 + δq 1) (13) q 1 K 1 q 0 p 0 p 0 p 1 Optimal K 1 etermination why K 1? It has been ecie to optimize the conuctivity of the hyaulic mean K for essentially two reasons: - The relationship between K an the eformation u is crucial for it's use in the poroelastic theory. It is as a matter of fact the mainly responsible for the interation between ui an pharmaceutical - Both in literature an in our moel it is still not clear how to combine K with the eformation u. This stuy oes not oubt the exponential relationship that goes on between eformation an conuctivity, but wants to verify if the exponential curve that better approximates the experimental values (strictly etermine by a mathematical metho) is qualitatively similar to the one chosen in our moel( etermine with a physical-empiric metho). The experimental ata reporte in the 1were taken from the work of McGuire et al. an for each value of p inf (37,52,69 mmhg)... is calculate Q inf meia (Q I 0.15,Q II 1.4,Q III 0.35). 3
4 - χ 0 (Q I ) π 2 a 4 ( p0 )2 2 p0 + 8Qπa - χ 1 32π 2 p0 ɛr4 + 8Q I πr 2 ɛ - χ 2 32π 2 ɛr 4 ( p0 )2 + 8Q I πr 2 ɛ p0 We can now ene the lagrangiana enition as follows: L(p 1, K 1, α, b, c) J 1 a α ( 2 p p 1 r γ 2 p 1 + ( 2 p 0 r + 2 p 0 )K ( p0 K1 )( )) b (p 1(a) + u 1 (a) p0(a) ) c (p 1 (1) + u 1 (1) p0(1) ) where α, b, c are moltiples of Lagrange. imposing: Figure 1: experimental values Lagrangian approach To make it easier from now on we will omit remembering that all measurements are aimensional From the optimization we have K ott for every value of Q inf. For this reason the size that has to be minimize has been expresse as follows: J (Q(a ) Q I ) 2 + ξ 2 ˆ 1 a K 2 1 (14) The optimization has been one for the linear equations stopping at the 1 orer, furthermore, for reasons tie to the metho of functionality... is expresse as follows Q(a ) ˆ 1 a Q(r)δ D (15) where δ D is a known function elta i Dirac. Taking into account that an remembering that 14 linearization has stoppe on orer 1 the 14 becomes : J χ 0 + where: ˆ 1 a (χ 1 p 1 + χ 2K 1 )δ D + ξ 2 ˆ 1 a K 1 (16) L 0 (17) We obtain the conitions necessary for the resolution of the problem. Therefore we have: L p α 2 α r γ2 α 8πɛ α(a) 0 α(1) 0 b α (a)) c α (1) ((4π p0 r4 +Q I r 2 )δ D ) L K 1 0 K 1 ξ [ 8πɛ ((4π 2 p 0 r 4 + Q I p0 2 r2 )δ D ) + α( 2 p 0 r + 2 p 0 ) p0 2 (α )] L α 0 2 p 1 ( p0 + 2 r K1 )( p 1 γ 2 p 1 + ( 2 p 0 r + 2 p 0 ) 0 L b 0 p 1 (a) u(a) p0(a) L c 0 p 1 (1) u 1 (1) p0(1) please note that: 2 )K Deriving from p 1 we are able to etermine the three multiplicators of lagrange α, b, c - - Deriving in respect to K 1 we obtain K 1,ott 4
5 - eriving in repect to α we obtain the irect - eriving in respect tob e c we obtain the initial conitions of the irect Discretization of the metho The equation system generate from the 15 becomes iscretize by an explicit metho at the secon orer Figure 2: rst case an 6000 points Results As state before, we face the problem separately for each of the three values of the infusion range. Given the obvious impossibility to impose a punctual conition on the function it has become necessary the introuction of a gaussiana function such as f C(θ)e θr2 (elta i Dirac approssimato)....which woul have the eect of blocking the information on all of r given ierent weights for every point. In this way, acting onθ We can choose the strength of our control. Obviously accoring to the choices mae on δ D the calculus gri must be moi- e. To make the calculations easier, the stuy will be execute with a δ D the most ample possible t.c. 1 a δ D which oesn't have to move from the unit for more than a variation of 1.5%. Finally, for a particular case we must emonstrate the convergence of the results at the varying of δ D As long as the gri was chosen accurately. In this case stuy will be shown 3 ierent cases with a δ D ever more force. 1 a δ D ; 1 a δ D ; 1 a δ D 3 1. δ D1 : If we know ecie to keep on using a gri of 6000 points even forδ D3 we woul obtain something like this: Figure 3: thir case an 6000 points It is clear that if we wish to keep on using this gri we woul have to raise the number of points on the gri. The following will show a gri of points: 5
6 In the gure below are reporte the values of the solutions not optimize an not linearize for each of the three experiments (Q I, Q II, Q III ) confronte with the values of the optimization at the change of the parameter ξ: Case 1 (Q I ) Figure 4: fourth case an points Whileδ D2 has the following form (gri of 8000) Figure 5: secon case an 8000 points The optimization is obviously strongly inuence by the freeom that is given to the control. The parameter responsible for such choice is ξ,the more it assumes small values the more the control is free to go its own way espite the energy spent to make all this possible. Figure 6: Q for ierent values of ξ Which enlarge in the point of interest shows the valiity of the optimization: 6
7 Figure 7: zoom Figure 9: zoom The black curve represents the non linear solution the others have a value of ξ respectively of 1, 0.1, 0.05, 0.04, The gure beneath, is in reference to the ierent hyaulic conuctivity of the mean K: Beneath are reporte the values taken by J at the varying of the control parameter ξ: Figure 10: J for ierent values of ξ Figure 8: K for ierent values of ξ Which enlarge in the point of interest: In the en, to verify if the chosen gri is sucient the values are reporte for one of the cases calculate above before, with the 6000 point gri an then with a 1000 point gri 7
8 Figure 11: conversetion Figure 13: The problem has now reache perfect conversion. Case 2 Q II We report in the same orer of the rst case the results obtaine. Figure 12: Figure 14: 8
9 Figure 15: Figure 17: case 3 (Q III ) The optimization in the thir case has taken the following values: Figure 16: Figure 18: 9
10 Figure 19: Convergence The graphic shows the convergence of the compute when you increase the ots on the gri: Figure 20: 10
Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More information1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a
Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..
More informationLecture 2 - First order linear PDEs and PDEs from physics
18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More information3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects
IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More information1 Boas, p. 643, problem (b)
Physics 6C Solutions to Homework Set #6 Fall Boas, p. 643, problem 3.5-3b Fin the steay-state temperature istribution in a soli cyliner of height H an raius a if the top an curve surfaces are hel at an
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationECE341 Test 2 Your Name: Tue 11/20/2018
ECE341 Test Your Name: Tue 11/0/018 Problem 1 (1 The center of a soli ielectric sphere with raius R is at the origin of the coorinate. The ielectric constant of the sphere is. The sphere is homogeneously
More informationA note on the Mooney-Rivlin material model
A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel
More informationS10.G.1. Fluid Flow Around the Brownian Particle
Rea Reichl s introuction. Tables & proofs for vector calculus formulas can be foun in the stanar textbooks G.Arfken s Mathematical Methos for Physicists an J.D.Jackson s Classical Electroynamics. S0.G..
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More information(3-3) = (Gauss s law) (3-6)
tatic Electric Fiels Electrostatics is the stuy of the effects of electric charges at rest, an the static electric fiels, which are cause by stationary electric charges. In the euctive approach, few funamental
More informationSemianalytical method of lines for solving elliptic partial dierential equations
Chemical Engineering Science 59 (2004) 781 788 wwwelseviercom/locate/ces Semianalytical metho of lines for solving elliptic partial ierential equations Venkat R Subramanian a;, Ralph E White b a Department
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationKramers Relation. Douglas H. Laurence. Department of Physical Sciences, Broward College, Davie, FL 33314
Kramers Relation Douglas H. Laurence Department of Physical Sciences, Browar College, Davie, FL 333 Introuction Kramers relation, name after the Dutch physicist Hans Kramers, is a relationship between
More informationPARALLEL-PLATE CAPACITATOR
Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationNon-Equilibrium Continuum Physics TA session #10 TA: Yohai Bar Sinai Dislocations
Non-Equilibrium Continuum Physics TA session #0 TA: Yohai Bar Sinai 0.06.206 Dislocations References There are countless books about islocations. The ones that I recommen are Theory of islocations, Hirth
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationCrack-tip stress evaluation of multi-scale Griffith crack subjected to
Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationIn the usual geometric derivation of Bragg s Law one assumes that crystalline
Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationEXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL
Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
More informationMathcad Lecture #5 In-class Worksheet Plotting and Calculus
Mathca Lecture #5 In-class Worksheet Plotting an Calculus At the en of this lecture, you shoul be able to: graph expressions, functions, an matrices of ata format graphs with titles, legens, log scales,
More informationCharge { Vortex Duality. in Double-Layered Josephson Junction Arrays
Charge { Vortex Duality in Double-Layere Josephson Junction Arrays Ya. M. Blanter a;b an Ger Schon c a Institut fur Theorie er Konensierten Materie, Universitat Karlsruhe, 76 Karlsruhe, Germany b Department
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationSIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES
SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES Jan Macek, Miloš Polášek Czech Technical University in Prague, Josef Božek Research Center Introuction Improvement of emissions from reciprocating internal
More informationAnalytical accuracy of the one dimensional heat transfer in geometry with logarithmic various surfaces
Cent. Eur. J. Eng. 4(4) 014 341-351 DOI: 10.478/s13531-013-0176-8 Central European Journal of Engineering Analytical accuracy of the one imensional heat transfer in geometry with logarithmic various surfaces
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationFluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall.
Flui Mechanics EBS 189a. Winter quarter, 4 units, CRN 52984. Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 eihmeyer Hall. Course Description: xioms of flui mechanics, flui statics,
More informationPhysics 170 Week 7, Lecture 2
Physics 170 Week 7, Lecture 2 http://www.phas.ubc.ca/ goronws/170 Physics 170 203 Week 7, Lecture 2 1 Textbook Chapter 12:Section 12.2-3 Physics 170 203 Week 7, Lecture 2 2 Learning Goals: Learn about
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationSYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d
SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)
Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationMONTE CARLO SIMULATION OF VIRTUAL COMPTON SCATTERING AT MAMI L. Van Hoorebeke c, J. Berthot b,p.y. Bertin b, V. Breton b, W.U. Boeglin g,r. Bohm, N. Degrane c, N. D'Hose a, M. Distler, J.E. Ducret a, R.
More informationTMA4195 Mathematical modelling Autumn 2012
Norwegian University of Science an Technology Department of Mathematical Sciences TMA495 Mathematical moelling Autumn 202 Solutions to exam December, 202 Dimensional matrix A: τ µ u r m - - s 0-2 - - 0
More informationTutorial Test 5 2D welding robot
Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationApproaches for Predicting Collection Efficiency of Fibrous Filters
Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State
More informationCalculus Class Notes for the Combined Calculus and Physics Course Semester I
Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average
More informationA Second Time Dimension, Hidden in Plain Sight
A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will
More informationDAMTP 000/NA04 On the semi-norm of raial basis function interpolants H.-M. Gutmann Abstract: Raial basis function interpolation has attracte a lot of
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports On the semi-norm of raial basis function interpolants H.-M. Gutmann DAMTP 000/NA04 May, 000 Department of Applie Mathematics an Theoretical Physics Silver
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationElectromagnetic surface and line sources under coordinate transformations
PHYSICAL REVIEW A 8, 3382 29 Electromagnetic surface an line sources uner coorinate transformations Steven A. Cummer,* Nathan Kuntz, an Bogan-Ioan Popa Department of Electrical an Computer Engineering
More informationCrack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion
Focuse on Mechanical Fatigue of Metals Crack onset assessment near the sharp material inclusion tip by means of moifie maximum tangential stress criterion Onřej Krepl, Jan Klusák CEITEC IPM, Institute
More informationMAE 210A FINAL EXAM SOLUTIONS
1 MAE 21A FINAL EXAM OLUTION PROBLEM 1: Dimensional analysis of the foling of paper (2 points) (a) We wish to simplify the relation between the fol length l f an the other variables: The imensional matrix
More informationHomework 2 Solutions EM, Mixture Models, PCA, Dualitys
Homewor Solutions EM, Mixture Moels, PCA, Dualitys CMU 0-75: Machine Learning Fall 05 http://www.cs.cmu.eu/~bapoczos/classes/ml075_05fall/ OUT: Oct 5, 05 DUE: Oct 9, 05, 0:0 AM An EM algorithm for a Mixture
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More informationAn analytical investigation into filmwise condensation on a horizontal tube in a porous medium with suction at the tube surface
Heat Mass Transfer (29) 45:355 361 DOI 1.17/s231-8-436-y ORIGINAL An analytical investigation into filmwise conensation on a horizontal tube in a porous meium with suction at the tube surface Tong Bou
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationDebond crack growth in fatigue along fiber in UD composite with broken fibers
Debon crack growth in atigue along iber in UD composite with broken ibers Johannes Eitzenberger an Janis arna Dept o Applie Physics an Mechanical Engineering Lulea University o Technology, Sween SE 97
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationQUANTUMMECHANICAL BEHAVIOUR IN A DETERMINISTIC MODEL. G. t Hooft
QUANTUMMECHANICAL BEHAVIOUR IN A DETERMINISTIC MODEL G. t Hooft Institute for Theoretical Physics University of Utrecht, P.O.Box 80 006 3508 TA Utrecht, the Netherlans e-mail: g.thooft@fys.ruu.nl THU-96/39
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More information