2 The wave equation in the finite element formalism
|
|
- William Turner
- 5 years ago
- Views:
Transcription
1 1 Introduction In flat space-time, the constant speed wave equation ( 2 2 t )φ = ρ (1) is simple enough, but when we let space-time curve, the equation becomes very complex. This is especially the case for the vector wave equation, which I will not have the time to consider in this project, but also, to a lesser degree, for the scalar equation. However, space can always be approximated as flat locally, and so as long as one looks at a small enough element of space, it is sufficient to solve the flat-space equations in that region. The goal if this project is to see if this approach combined with the finite element method satisfactory solves the scalar wave equation for the Schwarzschild space-time. 2 The wave equation in the finite element formalism Using the Galerkin method to minimize the residual for each element, we get Ω e d x j (φ j ( 2 N j )N i φ j N j N i ) = ρn i d x (2) Ω e M φ + Kφ =b (3) with M i j = Ω e d xn i N j, K i j = Ω e d x( N i ) ( N j ) and b i = Ω e d xρn i. The surface term has been set equal to zero here, since I will enforce zero boundary conditions here. Since I will be solving the wave equation the open space around a black hole, there shouldn t really be any boundary - space should extent infinitly. The simulation, however, has to be finite, and so there must be a boundary with appropriate conditions. Simply enforcing zero derivative and zero value at the border will not give correct behavior; in this case the field will be reflected with opposite phase from the border. This is clearly wrong: the border is just a random line in space, after all. What is needed is non-reflecting boundary conditions. These are very hard to implement correctly, but a common approximation is to extend the area considered slightly, and introduce a damping term (a term proportional to φ) in the new area only. If the damping isn t too dramatic (sudden damping would be akin to a wall, and cause reflection), a wave headed out of the area will be gradually damped to negligibility before being reflected by the boundary. 1
2 The analytical equation with a simple damping term becomes 2 tφ =ρ + 2 φ β φ (4) which leads to the discrete equation M φ + B φ + Kφ =b (5) with B i j = Ω e d xβn i N j. We can now use finite difference in time, and get ( t 2 b =M φ l+1 2φ l +φ l 1) + t ( 2 B φ l+1 φ l 1) + t 2 Kφ l φ l+1 =2φ l φ l 1 + M + t 2 B }{{} A 1 (6) ( ) [B + tk]φ l + Bφ l 1 + tb t Using lumped M and B, A becomes diagonal, and so this should be efficient to solve. (7) 3 Test of non-reflection boundary conditions I modified the fem/wave1 program to include a source and damping term. The damping variable beta is a scalar field set to inside the main area of the simulation, and increasing linearly to a given value at the end of the extended area. To make any reflections more apparent, I let the source exist only a t =, generating a wave that, ideally, would be sharp and circular, with both outside and inside the circle, propagating at constant speed from the source. I generated a homogenous circular grid t1.grid with 1 elements and 1 border nodes using the command makegrid --batch -m PreproStdGeom -g DISK r= degrees=36 \ -p e=elmt3n2d nel=1 resol=1 --casename nreftest" solving it for 4 time steps of size.1 with the input file nreftest.i. It turns out that the initial wave propagates outwards while dispersing, leaving a trail of concentric, weaker, circular waves. The effect of the damping can be measured by measuring the amplitude of the reflected waves as they converge on the center. Varying the damping parameter by orders of magnitude around 1, I found β max = 1 to be ideal, with a reflected 1 wave of amplitude about 1 of the original wave. Figure 1 to 4 illustrate this. Be ware that the scale on these figures is not the same. 2
3 u at time= e 17 Figur 1: Flat-space wave equation with damping outside r = 1. t = 1. u at time= Figur 2: Flat-space wave equation with damping outside r = 1. t = 11. 3
4 u at time= e e Figur 3: Flat-space wave equation with damping outside r = 1. t = 1. u at time= e Figur 4: Flat-space wave equation with damping outside r = 1. t = 1. 4
5 u at time= e Figur 5: Pulse after 4 timeunits in big grid with dense center but coarse otherwise. 3.1 Other ways of simulating infinite space Another way of simulating infinite space (which is the point of nonreflecting boundary conditions), is to make a huge but coarse grid, and gradually refine it towards the middle, where the interesting area is. I have not studied this alternative closely, but a test with a coarse grid of originally 1 elements with radius 1, refined at radii 25, 2 and shows that there is quite a bit of reflection even with this model. Figure 5 shows how a pulse starting in the center at t = looks like at t = 4. It seems like more of the wave has been reflected than what is still continuing outwards. I will stick with damping in the following. 4 Introducing space-time curvature Until now, the finite element simulation has assumed that space is flat, that is, that every element has time running at the same speed, and space being stretched equally. In the Schwarzschild metric, this is not the case. 5
6 The metric (for 2D space with rectangular coordinates) is given by g = 1 µ 2 r 2 µr 2 x 2 + µy 2 xy R r, r = x 2 + y 2, xy R r µx 2 + y 2 µ = 1 R r (8) This metric can be seperated into a time and space part, ( g tt = 1 R ) r γ = 1 ( x 2 + µy 2 xy R ) µr 2 r xy R r µx 2 + y 2 γ 1 = 1 ( µx 2 + y 2 xy R ) r 2 r xy R r x 2 + µy 2 (9) (1) (11) The flat-space method used earlier should still be valid at the element level, but a correction is needed to stitch these correctly together. As seen from the global Schwarzschild system, the locally flat coordinate system of each element is compressed or streteched corresponding to a linear coordinate transformation dx i = Λ i i dxi, with Λ 1i i Λ 1 j j γ {}}{ i j = γ i j, where marked coordinates belong the locally flat system. The last equation is, in matrix form, Λ 2 = γ, so it seems like one needs to find a sqare root of the matrix γ (can be done as an eigenvalue problem). This, however, will turn out to be unneccesary. The translation between these global and local coordinates is entierly similar to that already present in the finite element formalism, and we can reuse results from here. Every volume element must be inflated by the determinant of the transformation matrix, so d x Λ d x = γ d x, and every derivate must be modified according to i = Λ 1i i i. In our case, s only occur in the combination ( N i ) ( N j ) =Λ 1 k k k N i Λ 1kl l N j (12) =Λ 1Tl k Λ 1 k k k N i l N j (13) =γ 1l k k N i l N j (14) Finally, the time speed must be taken into consideration. This is a simple as letting dt = µdt. Since B occurs everywhere with a factor t and K and b with t 2, time speed is accounted for by including respectively µ and µ in their integrals. This is all taken care of in the functionintegrands. δ i j 6
7 u at time= e 5 Figur 6: R = 3, particle at 5, damping from 1, t = 1. 5 Results To see how waves propagate in the Schwarzschild geometry, I made a new run run with almost the same settings as the last run, but now with a Schwarzschild radius of 3 and a particle at a radius of 5, which is still, for the time being, only present at t =. The input file for this run is hole1.i, and figure 6 to 11 display the result. This simulation displays many of the aspects of what would really happen: The wave propagates only slowly in the direction of the black hole, since time runs slower there, and space is compressed there. Instead, the wave front winds itself around the horizon, meets itself on the other side of the hole, and continues winding up while diluting. In the end, it seems like a relic field is trapped at the horizen, while the field is practically elsewhere, which also seems reasonable. Closer inspection, however, reveals a problem. The relic field turns out to be inside the horizon, where the field isn t supposed to be able to go in finite time, and where the time speed in this simulation is set to be. This must be due to elements that cross the horizon. Some of the integration points of an element lie outside the horizon, then the basis functions in that element will get non-zero contributions. In that case, this problem should be avoidable by ensuring that element boundaries line up with the horizon. If the DISK_WITH_HOLE grid generator supported full circles, this would have been easy to achieve, while at the same time avoiding wasting time 7
8 u at time= Figur 7: R = 3, particle at 5, damping from 1, t = 6. u at time= Figur 8: R = 3, particle at 5, damping from 1, t = 11. 8
9 u at time= Figur 9: R = 3, particle at 5, damping from 1, t = 16. u at time= Figur 1: R = 3, particle at 5, damping from 1, t = 21. 9
10 u at time= Figur 11: R = 3, particle at 5, damping from 1, t = 31. on elements inside the horizon. Lacking the means to line up the element boundaries with the black hole horizon, an alternative could be to look at the horizon as a boundary, and force the field to be there, even though it is in the middle of the grid, and not a real boundary. One could also increase the grid resolution around the horizon, so there is a buffer of elements with very low time speed between the wave and the horizon, which could prevent the wave from reaching the horizon during the simulation. This would most easily be achieved with an adaptive grid with refinements based on time speed, which is what I originally implemented, only to discover that grid refinement is an extension of diffpack, which didn t seem to be available. Instead, I will use makegrid to make the grid extra dense in this area. I ran new simulation with doubly refined grid within r = 3.2 and forced value of inside the field, but otherwise identical settings (hole2.i). Figure 12 to 14 show the results. Increasing the density of the grid around the horizon helped. It is now more obvious that the field is compressed as it enters an area with slower time speed (much like the density of cars on a road increases when the speed limit decreases) and more compressed space. Figure 12 illustrates this nicely. The next figures, however, reveal that forcing the field to be zero inside the horizon leads to reflection. I will still use this forcing in the following section, to avoid the accumulation of the field inside horizon. 1
11 u at time= Figur 12: R = 3, extra dense within 3.2, forced zero inside hole, damping from 1, t = 6. u at time= Figur 13: R = 3, extra dense within 3.2, forced zero inside hole, damping from 1, t =
12 u at time= Figur 14: R = 3, extra dense within 3.2, forced zero inside hole, damping from 1, t = Other ways to avoid the horizon problem One could also try to avoid getting nonzero values inside the black hole by using a coordinate system where only the outside of the hole exist, for example by letting r a = r R, and introducing coordinates x = a cosθ, ȳ = a sinθ. The program includes an option for doing this, but the results are unphysical, probably due to horizon being replaced by a single point, and I will not discuss this further here. 5.2 The field of a stationary particle outside a black hole By letting the source stay inside the simulation region during the whole simulation, instead of just initially, I can calculate the field of a stationary particle outside a black hole. The input for this simulation is in hole3.i. The field took about 1 timeunits to stabilize and figure illustrates the final, stationary field. The same situation without zero inside the horizon can be seen in figure 16. It is unclear which of these is most correct. The first stabilizes relatively quickly, but does not display a strong field near the horizon. The latter has this, but also has artifacts inside the horizon, which prevent the solution from stabilizing completely. 12
13 u at time= Figur : The field of a stationary particle at r = 5 outside a black hole with R = 3 with forcing inside the horizon. 5.3 The field of an orbiting source outside a black hole Using the same parameters as in the last simulation, but letting the source orbit, results in a field that is surprisingly similar, even for high rotation speeds, as figure 17 illustrates. I am not sure if this is physical or not (a moving particle needs the vector wave equation in the case of electromagnetism in any case). 6 Conclusion These simulations have demonstrated that it is possible to get reasonablelooking results for the scalar wave equation in curved space-time by using this simple method of solving the flat-space equation locally, and applying a transformation to get the global solution. Any involvement of the horrible and numerous Christoffel symbols and their derivatives can thus be avoided. There were some difficulties with this method, particularily with the boundary conditions, both the ones at big and small radius, and it is not which of the approaches used for handling the horizon is closest to an exact solution. 13
14 u at time= Figur 16: The field of a stationary particle at r = 5 outside a black hole with R = 3 without forcing inside the horizon. 14
15 u at time= Figur 17: The field of an orbiting particle at r = 5 outside a black hole with R = 3 with forcing inside the horizon.
Quick Questions. 1. Two charges of +1 µc each are separated by 1 cm. What is the force between them?
92 3.10 Quick Questions 3.10 Quick Questions 1. Two charges of +1 µc each are separated by 1 cm. What is the force between them? 0.89 N 90 N 173 N 15 N 2. The electric field inside an isolated conductor
More informationCosmology and particle physics
Cosmology and particle physics Lecture notes Timm Wrase Lecture 9 Inflation - part I Having discussed the thermal history of our universe and in particular its evolution at times larger than 10 14 seconds
More informationSTEP Support Programme. Mechanics STEP Questions
STEP Support Programme Mechanics STEP Questions This is a selection of mainly STEP I questions with a couple of STEP II questions at the end. STEP I and STEP II papers follow the same specification, the
More informationCosmological Issues. Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab
Cosmological Issues 1 Radiation dominated Universe Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab ρ 0 0 0 T ab = 0 p 0 0 0 0 p 0 (1) 0 0 0 p We do not often
More informationAstronomy 421. Lecture 24: Black Holes
Astronomy 421 Lecture 24: Black Holes 1 Outline General Relativity Equivalence Principle and its Consequences The Schwarzschild Metric The Kerr Metric for rotating black holes Black holes Black hole candidates
More informationEM Waves in Media. What happens when an EM wave travels through our model material?
EM Waves in Media We can model a material as made of atoms which have a charged electron bound to a nucleus by a spring. We model the nuclei as being fixed to a grid (because they are heavy, they don t
More informationCalculus II. Calculus II tends to be a very difficult course for many students. There are many reasons for this.
Preface Here are my online notes for my Calculus II course that I teach here at Lamar University. Despite the fact that these are my class notes they should be accessible to anyone wanting to learn Calculus
More informationMetric-based classifiers. Nuno Vasconcelos UCSD
Metric-based classifiers Nuno Vasconcelos UCSD Statistical learning goal: given a function f. y f and a collection of eample data-points, learn what the function f. is. this is called training. two major
More informationBonus Section II: Solving Trigonometric Equations
Fry Texas A&M University Math 150 Spring 2017 Bonus Section II 260 Bonus Section II: Solving Trigonometric Equations (In your text this section is found hiding at the end of 9.6) For what values of x does
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationUniversity of Maryland Department of Physics. Spring 2009 Final Exam 20. May (175 points) Post grades on web? (Initial, please) Yes No
University of Maryland Department of Physics Physics 122 20. May 2009 (175 points) Post grades on web? (Initial, please) Yes No (If you agree, I will post your grades and your detailed scores for each
More informationElectro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras
Electro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras Lecture - 7 Gauss s Law Good morning. Today, I want to discuss two or three
More informationPhys101 First Major-111 Zero Version Monday, October 17, 2011 Page: 1
Monday, October 17, 011 Page: 1 Q1. 1 b The speed-time relation of a moving particle is given by: v = at +, where v is the speed, t t + c is the time and a, b, c are constants. The dimensional formulae
More informationBlack Hole Astrophysics Chapters 7.5. All figures extracted from online sources of from the textbook.
Black Hole Astrophysics Chapters 7.5 All figures extracted from online sources of from the textbook. Recap the Schwarzschild metric Sch means that this metric is describing a Schwarzschild Black Hole.
More informationFinal Physics of Schwarzschild
Physics 4 Lecture 32 Final Physics of Schwarzschild Lecture 32 Physics 4 Classical Mechanics II November 6th, 27 We have studied a lot of properties of the Schwarzschild metric I want to finish with a
More informationCosmic Bubble Collisions
Outline Background Expanding Universe: Einstein s Eqn with FRW metric Inflationary Cosmology: model with scalar field QFTà Bubble nucleationà Bubble collisions Bubble Collisions in Single Field Theory
More informationKinematics vs. Dynamics
NEWTON'S LAWS Kinematics vs. Dynamics Kinematics describe motion paths mathematically No description of the matter that may travel along motion path Dynamics prediction of motion path(s) of matter Requires
More informationLab I. 2D Motion. 1 Introduction. 2 Theory. 2.1 scalars and vectors LAB I. 2D MOTION 15
LAB I. 2D MOTION 15 Lab I 2D Motion 1 Introduction In this lab we will examine simple two-dimensional motion without acceleration. Motion in two dimensions can often be broken up into two separate one-dimensional
More informationSums of Squares (FNS 195-S) Fall 2014
Sums of Squares (FNS 195-S) Fall 014 Record of What We Did Drew Armstrong Vectors When we tried to apply Cartesian coordinates in 3 dimensions we ran into some difficulty tryiing to describe lines and
More informationNotes on General Relativity Linearized Gravity and Gravitational waves
Notes on General Relativity Linearized Gravity and Gravitational waves August Geelmuyden Universitetet i Oslo I. Perturbation theory Solving the Einstein equation for the spacetime metric is tremendously
More informationChapter S3 Spacetime and Gravity Pearson Education, Inc.
Chapter S3 Spacetime and Gravity What are the major ideas of general relativity? Spacetime Special relativity showed that space and time are not absolute. Instead, they are inextricably linked in a four-dimensional
More informationTOPIC V BLACK HOLES IN STRING THEORY
TOPIC V BLACK HOLES IN STRING THEORY Lecture notes Making black holes How should we make a black hole in string theory? A black hole forms when a large amount of mass is collected together. In classical
More informationAnalytical Physics 1B Lecture 5: Physical Pendulums and Introduction to Mechanical Waves
Analytical Physics 1B Lecture 5: Physical Pendulums and Introduction to Mechanical Waves Sang-Wook Cheong Friday, February 16 th, 2017 Two Exam 1 Questions with errors Correct answer: L = r X p = (2000
More informationLab I. 2D Motion. 1 Introduction. 2 Theory. 2.1 scalars and vectors LAB I. 2D MOTION 15
LAB I. 2D MOTION 15 Lab I 2D Motion 1 Introduction In this lab we will examine simple two-dimensional motion without acceleration. Motion in two dimensions can often be broken up into two separate one-dimensional
More informationDiscrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 27
CS 70 Discrete Mathematics for CS Spring 007 Luca Trevisan Lecture 7 Infinity and Countability Consider a function f that maps elements of a set A (called the domain of f ) to elements of set B (called
More information. (70.1) r r. / r. Substituting, we have the following equation for f:
7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical
More information4.4: Optimization. Problem 2 Find the radius of a cylindrical container with a volume of 2π m 3 that minimizes the surface area.
4.4: Optimization Problem 1 Suppose you want to maximize a continuous function on a closed interval, but you find that it only has one local extremum on the interval which happens to be a local minimum.
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationLAB Exercise #4 - Answers The Traction Vector and Stress Tensor. Introduction. Format of lab. Preparation reading
LAB Exercise #4 - Answers The Traction Vector and Stress Tensor Due: Thursday, 26 February 2009 (Special Thanks to D.D. Pollard who pioneered this exercise in 1991) Introduction Stress concentrations in
More informationBLACKHOLE WORMHOLE THEORY
BLACKHOLE WORMHOLE THEORY By - ASHU PRAKASH Black hole, a name which has infinite knowledge to define, but very difficult to define. What is a black hole? In general, a black hole is a gravitationally
More informationSurvey of Astrophysics A110
Black Holes Goals: Understand Special Relativity General Relativity How do we observe black holes. Black Holes A consequence of gravity Massive neutron (>3M ) cannot be supported by degenerate neutron
More informationASTR2050 Spring In this class we will cover: Hints: Escape Velocity. Relativity and the Equivalence Principle Visualization of Curved Spacetime
ASTR2050 Spring 2005 Lecture 11am 8 March 2005 In this class we will cover: Hints: Escape Velocity Black Holes Relativity and the Equivalence Principle Visualization of Curved Spacetime 1 Escape Velocity
More informationXIII. The Very Early Universe and Inflation. ASTR378 Cosmology : XIII. The Very Early Universe and Inflation 171
XIII. The Very Early Universe and Inflation ASTR378 Cosmology : XIII. The Very Early Universe and Inflation 171 Problems with the Big Bang The Flatness Problem The Horizon Problem The Monopole (Relic Particle)
More informationFRW cosmology: an application of Einstein s equations to universe. 1. The metric of a FRW cosmology is given by (without proof)
FRW cosmology: an application of Einstein s equations to universe 1. The metric of a FRW cosmology is given by (without proof) [ ] dr = d(ct) R(t) 1 kr + r (dθ + sin θdφ ),. For generalized coordinates
More informationVolume in n Dimensions
Volume in n Dimensions MA 305 Kurt Bryan Introduction You ve seen that if we have two vectors v and w in two dimensions then the area spanned by these vectors can be computed as v w = v 1 w 2 v 2 w 1 (where
More informationAAPT UNITED STATES PHYSICS TEAM AIP 2008
8 F = ma Exam AAPT UNITED STATES PHYSICS TEAM AIP 8 8 F = ma Contest 5 QUESTIONS - 75 MINUTES INSTRUCTIONS DO NOT OPEN THIS TEST UNTIL YOU ARE TOLD TO BEGIN Use g = N/kg throughout this contest. You may
More informationLecture 14: Cosmological Principles
Lecture 14: Cosmological Principles The basic Cosmological Principles The geometry of the Universe The scale factor R and curvature constant k Comoving coordinates Einstein s initial solutions 3/28/11
More informationFlat Geometry. Spherical Geometry
The Geometry of the Universe What does the constant k in the Friedmann equation really mean? In this lecture we will follow Chapter 4 of Liddle to show that it has close connections with the geometry of
More informationLast Update: March 1 2, 201 0
M ath 2 0 1 E S 1 W inter 2 0 1 0 Last Update: March 1 2, 201 0 S eries S olutions of Differential Equations Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections
More informationCURRENT-CARRYING CONDUCTORS / MOVING CHARGES / CHARGED PARTICLES IN CIRCULAR ORBITS
PHYSICS A2 UNIT 4 SECTION 4: MAGNETIC FIELDS CURRENT-CARRYING CONDUCTORS / MOVING CHARGES / CHARGED PARTICLES IN CIRCULAR ORBITS # Questions MAGNETIC FLUX DENSITY 1 What is a magnetic field? A region in
More informationAsk class: what is the Minkowski spacetime in spherical coordinates? ds 2 = dt 2 +dr 2 +r 2 (dθ 2 +sin 2 θdφ 2 ). (1)
1 Tensor manipulations One final thing to learn about tensor manipulation is that the metric tensor is what allows you to raise and lower indices. That is, for example, v α = g αβ v β, where again we use
More informationComputational Problem: Keplerian Orbits
Computational Problem: Keplerian Orbits April 10, 2006 1 Part 1 1.1 Problem For the case of an infinite central mass and an orbiting test mass, integrate a circular orbit and an eccentric orbit. Carry
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
More information(x x 0 ) 2 + (y y 0 ) 2 = ε 2, (2.11)
2.2 Limits and continuity In order to introduce the concepts of limit and continuity for functions of more than one variable we need first to generalise the concept of neighbourhood of a point from R to
More information1 Introduction. 1.1 The Problem. 1.2 Historical Perspective
1 Introduction 1.1 The Problem Consider a spherical surface with radius. Suppose that the surface has been cut by four (4) right circular cylinders with parallel axes. Consider a plane that passes by the
More informationSEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS
3 Underwater propagation 3. Ray acoustics 3.. Relevant mathematics We first consider a plane wave as depicted in figure. As shown in the figure wave fronts are planes. The arrow perpendicular to the wave
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More information(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf =
PHY 302 K. Solutions for problem set #9. Non-textbook problem #1: (a) Rotation frequency of 1 Hz means one revolution per second, or 60 revolutions per minute (RPM). The pre-lp vinyl disks rotated at 78
More informationLecture 1 General relativity and cosmology. Kerson Huang MIT & IAS, NTU
A Superfluid Universe Lecture 1 General relativity and cosmology Kerson Huang MIT & IAS, NTU Lecture 1. General relativity and cosmology Mathematics and physics Big bang Dark energy Dark matter Robertson-Walker
More information1 Multiplicity of the ideal gas
Reading assignment. Schroeder, section.6. 1 Multiplicity of the ideal gas Our evaluation of the numbers of microstates corresponding to each macrostate of the two-state paramagnet and the Einstein model
More informationAtom Model and Relativity
Atom Model and Relativity Kimmo Rouvari September 8, 203 Abstract What is the theoretical explanation for fine structure? What is the mechanism behind relativity? These questions have bothered numerous
More informationLecture XIX: Particle motion exterior to a spherical star
Lecture XIX: Particle motion exterior to a spherical star Christopher M. Hirata Caltech M/C 350-7, Pasadena CA 95, USA Dated: January 8, 0 I. OVERVIEW Our next objective is to consider the motion of test
More informationExamples and Explanations for a Flat Earth 16/10/2006 Andy Roff
Why do we believe the Earth is a sphere? You ve been told this by your teachers and parents since you were small, but have you really got any proof? Can you prove by sound reasoning and factual evidence
More information2.1 Basics of the Relativistic Cosmology: Global Geometry and the Dynamics of the Universe Part I
1 2.1 Basics of the Relativistic Cosmology: Global Geometry and the Dynamics of the Universe Part I 2 Special Relativity (1905) A fundamental change in viewing the physical space and time, now unified
More informationOutline. Black Holes. Schwartzchild radius River Model of a Black Hole Light in orbit Tidal forces
Outline Black Holes Schwartzchild radius River Model of a Black Hole Light in orbit Tidal forces Black Holes Black Holes What happens as the star shrinks / its mass increases? How much can spacetime be
More informationSpan and Linear Independence
Span and Linear Independence It is common to confuse span and linear independence, because although they are different concepts, they are related. To see their relationship, let s revisit the previous
More informationLecture 17 - The Secrets we have Swept Under the Rug
1.0 0.5 0.0-0.5-0.5 0.0 0.5 1.0 Lecture 17 - The Secrets we have Swept Under the Rug A Puzzle... What makes 3D Special? Example (1D charge distribution) A stick with uniform charge density λ lies between
More informationCircular Motion.
1 Circular Motion www.njctl.org 2 Topics of Uniform Circular Motion (UCM) Kinematics of UCM Click on the topic to go to that section Period, Frequency, and Rotational Velocity Dynamics of UCM Vertical
More informationCharged objects in Conducting Fluids
Charged objects in Conducting Fluids Net charge in a sphere of radius λ D is approximately zero. λ D 2 = ε 0κ k B T c 0 e 2 Z 2 k B T k c = 1 / 4πε 0 c 0 Z e κ Thermal energy (Joules) Coulomb constant
More informationUniversity of Maryland Department of Physics
Spring 2002 University of Maryland Department of Physics Laura Lising Physics 122 April 24, 2003 Exam #2 Solutions Multiple choice questions. Just the answer counts for these. (8 points each) 1) Suppose
More informationMath 20C Homework 2 Partial Solutions
Math 2C Homework 2 Partial Solutions Problem 1 (12.4.14). Calculate (j k) (j + k). Solution. The basic properties of the cross product are found in Theorem 2 of Section 12.4. From these properties, we
More informationMA 262, Fall 2017, Final Version 01(Green)
INSTRUCTIONS MA 262, Fall 2017, Final Version 01(Green) (1) Switch off your phone upon entering the exam room. (2) Do not open the exam booklet until you are instructed to do so. (3) Before you open the
More informationGCE AS and A Level. Physics A. AS exams 2009 onwards A2 exams 2010 onwards. Unit 4A: Approved specimen question paper. Version 1.3
GE S and Level Physics S exams 009 onwards exams 010 onwards Unit 4: pproved specimen question paper Version 1.3 General ertificate of Education 010 dvanced Examination abc version 1.3 PHYSIS Unit 4: Fields
More informationMr Green sees the shorter, straight, green path and Mr. Red sees the longer, curved, red path.
Mr Green sees the shorter, straight, green path and Mr. Red sees the longer, curved, red path. In an accelerated frame, time runs slow compared to a non-accelerated frame. The Equivalence Principle tells
More informationComputer Problems for Methods of Solving Ordinary Differential Equations
Computer Problems for Methods of Solving Ordinary Differential Equations 1. Have a computer make a phase portrait for the system dx/dt = x + y, dy/dt = 2y. Clearly indicate critical points and separatrices.
More informationLecture 4: Constructing the Integers, Rationals and Reals
Math/CS 20: Intro. to Math Professor: Padraic Bartlett Lecture 4: Constructing the Integers, Rationals and Reals Week 5 UCSB 204 The Integers Normally, using the natural numbers, you can easily define
More informationName Solutions Linear Algebra; Test 3. Throughout the test simplify all answers except where stated otherwise.
Name Solutions Linear Algebra; Test 3 Throughout the test simplify all answers except where stated otherwise. 1) Find the following: (10 points) ( ) Or note that so the rows are linearly independent, so
More informationcarroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general
http://pancake.uchicago.edu/ carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationElectromagnetic Forces on Parallel Current-
Page 1 of 5 Tutorial Models : Electromagnetic Forces on Parallel Current-Carrying Wires Electromagnetic Forces on Parallel Current- Carrying Wires Introduction One ampere is defined as the constant current
More informationSolving Differential Equations on 2-D Geometries with Matlab
Solving Differential Equations on 2-D Geometries with Matlab Joshua Wall Drexel University Philadelphia, PA 19104 (Dated: April 28, 2014) I. INTRODUCTION Here we introduce the reader to solving partial
More informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Essential physics for game developers Introduction The primary issues Let s move virtual objects Kinematics: description
More informationPHY 111L Activity 9 Moments of Inertia
PHY 111L Activity 9 Moments of Inertia Name: Section: ID #: Date: Lab Partners: TA initials: Objectives 1. Introduce moment of inertia for different objects 2. Understand the moment of inertia apparatus
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.01T Fall Term 2004 Experiment 06: Work, Energy and the Harmonic Oscillator Purpose of the Experiment: In this experiment you allow a cart
More informationPRELIMINARY EXAMINATION Department of Physics University of Florida Part A, January, 2012, 09:00 12:00. Instructions
Student ID Number: PRELIMINARY EXAMINATION Part A, January, 2012, 09:00 12:00 Instructions 1. You may use a calculator and CRC Math tables or equivalent. No other tables or aids are allowed or required.
More informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
More informationBREAKING SPAGHETTI. Yulia Peregud. Find the conditions under which dry spaghetti falling on a hard floor does not break.
BREAKING SPAGHETTI Yulia Peregud Belarusian State University Lyceum, Belarus Formulation of the problem: Find the conditions under which dry spaghetti falling on a hard floor does not break. Introduction:
More information0 questions at random and keep in order
Page 1 of 10 This chapter has 54 questions. Scroll down to see and select individual questions or narrow the list using the checkboxes below. 0 questions at random and keep in order s - (46) - (6) Fill
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationThe Cosmological Observables of Pre-Inflationary Bubble Collisions
The Cosmological Observables of Pre-Inflationary Bubble Collisions Spencer Chang (NYU) w/ M. Kleban, T. Levi 0712.2261 [hep-th] Also at Youtube, search for When Worlds Collide Trailer Cosmology Wealth
More informationSmooth Bang Theory: The Universe Steadily Gains Mass With Rate c 3 /G
Smooth Bang Theory: The Universe Steadily Gains Mass With Rate c 3 /G Abstract Mikhail Vlasov Email: vlasovm@hotmail.com October 31, 2015 According to spectroscopic data of light coming from distant galaxies
More informationAAPT UNITED STATES PHYSICS TEAM AIP 2017
2017 F = ma Exam 1 AAPT UNITED STATES PHYSICS TEAM AIP 2017 2017 F = ma Contest 25 QUESTIONS - 75 MINUTES INSTRUCTIONS DO NOT OPEN THIS TEST UNTIL YOU ARE TOLD TO BEGIN Use g = 10 N/kg throughout this
More informationElectricity & Magnetism Qualifier
Electricity & Magnetism Qualifier For each problem state what system of units you are using. 1. Imagine that a spherical balloon is being filled with a charged gas in such a way that the rate of charge
More informationThe Restricted 3-Body Problem
The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two
More informationPhysics 523, General Relativity Homework 7 Due Wednesday, 6 th December 2006
Physics 53, General elativity Homework 7 Due Wednesday, 6 th December 006 Jacob Lewis Bourjaily Problem Consider a gyroscope moving in circular orbit of radius about a static, spherically-symmetric planet
More informationQ: What does angular momentum mean? What is its equation for calculation?
Ch 10: Conservation of Angular Momentum Reeeeecap. Q: What does angular velocity mean? What is its symbol? A: The rate of change of angular displacement.. Q: What does angular acceleration mean? What is
More informationVibrations of string. Henna Tahvanainen. November 8, ELEC-E5610 Acoustics and the Physics of Sound, Lecture 4
Vibrations of string EEC-E5610 Acoustics and the Physics of Sound, ecture 4 Henna Tahvanainen Department of Signal Processing and Acoustics Aalto University School of Electrical Engineering November 8,
More informationv lim a t = d v dt a n = v2 R curvature
PHY 02 K. Solutions for Problem set # 6. Textbook problem 5.27: The acceleration vector a of the particle has two components, the tangential acceleration a t d v dt v lim t 0 t (1) parallel to the velocity
More informationSection 5-7 : Green's Theorem
Section 5-7 : Green's Theorem In this section we are going to investigate the relationship between certain kinds of line integrals (on closed paths) and double integrals. Let s start off with a simple
More informationQuest Chapter 09. Eliminate the obviously wrong answers. Consider what is changing: speed, velocity, some part of velocity? Choose carefully.
1 A dragster maintains a speedometer reading of 100 km/h and passes through a curve with a constant radius. Which statement is true? 1. The dragster rounded the curve at a changing speed of 100 km/h. 2.
More informationLecture 21 Reminder/Introduction to Wave Optics
Lecture 1 Reminder/Introduction to Wave Optics Program: 1. Maxwell s Equations.. Magnetic induction and electric displacement. 3. Origins of the electric permittivity and magnetic permeability. 4. Wave
More informationhas a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity.
http://preposterousuniverse.com/grnotes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been framed
More information8 Eigenvectors and the Anisotropic Multivariate Gaussian Distribution
Eigenvectors and the Anisotropic Multivariate Gaussian Distribution Eigenvectors and the Anisotropic Multivariate Gaussian Distribution EIGENVECTORS [I don t know if you were properly taught about eigenvectors
More information8. The Expanding Universe, Revisited
8. The Expanding Universe, Revisited A1143: History of the Universe, Autumn 2012 Now that we have learned something about Einstein s theory of gravity, we are ready to revisit what we have learned about
More informationEinstein Toolkit Workshop. Joshua Faber Apr
Einstein Toolkit Workshop Joshua Faber Apr 05 2012 Outline Space, time, and special relativity The metric tensor and geometry Curvature Geodesics Einstein s equations The Stress-energy tensor 3+1 formalisms
More informationPhysics 351, Spring 2015, Homework #6. Due at start of class, Friday, February 27, 2015
Physics 351, Spring 2015, Homework #6. Due at start of class, Friday, February 27, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationcarroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general
http://pancake.uchicago.edu/ carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been
More informationGeneral Relativity and Cosmology. The End of Absolute Space Cosmological Principle Black Holes CBMR and Big Bang
General Relativity and Cosmology The End of Absolute Space Cosmological Principle Black Holes CBMR and Big Bang The End of Absolute Space (AS) Special Relativity (SR) abolished AS only for the special
More informationComputer Science 324 Computer Architecture Mount Holyoke College Fall Topic Notes: Digital Logic
Computer Science 324 Computer Architecture Mount Holyoke College Fall 2007 Topic Notes: Digital Logic Our goal for the next few weeks is to paint a a reasonably complete picture of how we can go from transistor
More information