In this method, one defines
|
|
- Claude Freeman
- 5 years ago
- Views:
Transcription
1 1 Feautrier s Method for Radiative Transfer This project is to use Feautrier s method to solve the Gray Atmosphere problem for a plane-parallel atmosphere Although we have developed alternate techniques for solving this problem, the extension from the gray atmosphere situation to the general problem in which the opacities are frequency dependent becomes a formidable problem The technique developed here can be extended to the frequency-dependent case more easily A gray atmosphere lacks frequency dependence The problem is to find the intensity as a function of depth and angle; in particular, to find the outgoing intensity as a function of angle The equation of transfer is di τ, µ) µ = I τ, µ) S τ) = I τ, µ) B τ), where we can use the gray relation S = B We begin by assuming that Bτ) is known, and later implement a better alogrithm For this problem, assume it is given by B τ) = 3 4 F τ + 2 ) 3 where F is the constant) flux At the end of this project, you will compute F to determine how accurate this approximation is In the real world, you would then have to alter the above approximation to ensure that F remains constant through the atmosphere By such iterations, the full atmosphere model is then found Also, in the real world, S B and one has to make use of opacity information, which is also a function of density To determine the density, the equation of hydrostatic equilibrium has to be used However, in this project, we will neither consider this complication, nor do this iteration It is customary to divide the intensity into outgoing and ingoing streams, so that the variable µ is now a positive quantity: di τ, +µ) µ di τ, µ) µ In this method, one defines =I τ, +µ) B τ) =I τ, µ) B τ) u τ, µ) = 1 [I τ, +µ) + I τ, µ)] 2 v τ, µ) = 1 [I τ, +µ) I τ, µ)] 2
2 2 Thus, dv τ, µ) µ =u τ, µ) B τ) du τ, µ) µ =v τ, µ) We can combine them to eliminate v: µ 2d2 u τ, µ) 2 = u τ, µ) B τ) The boundary conditions on I are taken to be that of no incoming flux at the surface, and the diffusion approximation at the atmosphere s base: I, µ) = I τ, +µ) = [ db τ) B τ) + µ + µ 2d2 B τ) d 2 τ ] + τ We ignore second and higher derivatives of B The boundary conditions become du τ, µ) µ =v, µ) = u, µ) [ ] [ ] du τ, µ) db τ) [u τ, µ) + v τ, µ)] τ = u τ, µ) + µ = B τ) + µ τ 1) To solve these equations, we must discretize them Choose angles based upon Gauss-Legendre quadrature for ease, choose 3 positive and 3 negative values for µ j ) Choose optical depths ranging from to, say, 1, with a fixed interval of = 1 Thus τ i+1 τ i It would be better in general to choose ln τ for the independent variable for accuracy s sake, but this complicates the algebra and a linear grid is sufficient here We can write the system of equations for 1 < i < N) as A i,j u i 1,j + D i,j u i,j C i,j u i+1,j = E i 2) For 2 i N 1 one has, using central differencing, τ A i,j = µ 2 j / 2, D i,j = 1+2µ 2 j / 2, C i,j = µ 2 j / 2, E i = B i B τ i ) Near τ =, we expand u in a Taylor series: du τ, µ) u 2,j u 1,j d 2 u τ, µ) =u 1,j + u 1,j + ) 2 u µ 1,j B 1 j 2µ 2 j
3 3 Multiplying this by µ j /, we can identify A 1,j =, D 1,j = 1 + µ j / + 2µ j, C 1,j = µ j /, E 1 = B 1 2µ j At the other boundary, we find A N,j =µ j / 1/2, D N,j = 1/2 + µ j /, C N,j =, E N = B N 1/2 + µj / ) + B N 1 1/2 µj / ) Eq 2) is linear, and we can use a substitution scheme similar to that of the Henyey technique Assume for each j, using the shorthand u i = u i,j, By substitution, one finds u i = α i+1 + β i+1 u i+1 3) α i = E i 1 + α i 1 A i 1, D i 1 β i 1 A i 1 C β i = i 1 D i 1 β i 1 A i 1 Beginning at the surface, β 2 = C 1 /D 1 since A 1 = Now one can determine all α i s and β i s up to i = N + 1 In turn, one uses Eq 3) to obtain each u i from i = N to i = 1 After the u i s are determined, one reconstructs the flux: F τ i ) = µi τ i, µ) dµ = 2 j a j µ j [ I τi, +µ j ) I τi, µ j )] 4) =4 j a j µ j v i,j = 4 j a j µ 2 j u i,j Here, w j are the Gauss-Legendre weights, and the sums are only over postive values of j To find u i,j = du/) i,j, use the discretized relation u i,j = u i+1,j u i 1,j, 2 which is valid for 2 i N 1 For the boundaries Taylor-series expansions can be used: u 1,j = 4u 2,j 3u 1,j u 3,j 2, u N,j = 3u N,j + u N 2,j 4u N 1,j 2
4 4 For some reason, using the simpler result from the inner boundary condition that u 1,j = u 1,j/µ j produces an instability in the iterations described below How constant is F as a function of depth? Fig 1 shows that near the surface, the flux is especially poorly determined To ensure that F remains constant as a function of depth, it is necessary to alter the approximation B = 3F/4)τ + 2/3) This procedure is commonly known as Λ-iteration, since Λ is the operator that yields S Figure 1: The model flux F mod as a function of optical depth for cases of different angular resolutions, m = 1, 3, 5 Improvement with increasing m is slow In the gray atmosphere case, the constancy of the flux implies B = J, so one procedure would be to update B by using B τ) =J τ) = 1 +1 Idµ = 1 [ ) )] a 2 2 j I τ, +µj + I τ, µj 1 = a j u ) τ, µ j j j Once this value for B is used in the Feautrier scheme, new values for u i,j are determined This iteration must be repeated until convergence Unfortunately, typically thousands of iterations are necessary see Fig 2) This
5 5 happens because at large optical depth, we will always find J B no matter what, and the changes in B are exponentially small Note that the flux was not used in the computation of the correction of B Figure 2: Λ iteration for different angular resolutions Each curve shows the result of 1 successive iterations Results for m = 15) are shifted by +1-1); somewhat better convergence with increasing m is evident A better approach is the so-called Unsöld-Lucy procedure Begin with the transfer equation and its first two moments, in the gray case with no scattering, in which case S = B: µ di = I B, 1 df 4 = J B, dk = 1 4 F Since B is not the correct source function, F will not be constant Integrate the last of the above relations: K τ) = 1 F τ ) C, where C is a constant Now use the Eddington approximation J = 3K, and C can be found if we further approximate J) = F )/2, so J 3 4 F τ ) F )
6 6 The first moment equation becomes B τ) = J τ) 1 df τ) F τ ) F ) 1 4 df τ) Now, the correction Bτ) should be that which makes F constant: B τ) 3 4 F F + B τ), where F is the correct, constant, flux So, B τ) = 3 [ F τ ) F ] [F τ) F ] 1 d [F τ) F ] 4 Even though we used a number of approximations, the correction term uses only flux information, so that convergence to the proper solution is guaranteed This technique works well Fig 3) Note the expanded scale of the abcissa in Fig 3 Typically, only 5 or 6 iterations are adequate Figure 3: Unsöld-Lucy iteration Each curve represents the result of 1 successive iteration Results for m = 15) are shifted by +1-1)
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Zářivě (magneto)hydrodynamický seminář Ondřejov 20.03.2008 p. Outline 1.
More informationLecture 3 Numerical Solutions to the Transport Equation
Lecture 3 Numerical Solutions to the Transport Equation Introduction I There are numerous methods for solving the transport problem numerically. First we must recognize that we need to solve two problems:
More informationLecture 2 Solutions to the Transport Equation
Lecture 2 Solutions to the Transport Equation Equation along a ray I In general we can solve the static transfer equation along a ray in some particular direction. Since photons move in straight lines
More informationPHAS3135 The Physics of Stars
PHAS3135 The Physics of Stars Exam 2013 (Zane/Howarth) Answer ALL SIX questions from Section A, and ANY TWO questions from Section B The numbers in square brackets in the right-hand margin indicate the
More informationRadiation in climate models.
Lecture. Radiation in climate models. Objectives:. A hierarchy of the climate models.. Radiative and radiative-convective equilibrium.. Examples of simple energy balance models.. Radiation in the atmospheric
More informationModel Atmospheres. Model Atmosphere Assumptions
Model Atmospheres Problem: Construct a numerical model of the atmosphere to estimate (a) Variation of physical variables (T, P) with depth (b) Emergent spectrum in continuum and lines Compare calculated
More informationTRANSFER OF RADIATION
TRANSFER OF RADIATION Under LTE Local Thermodynamic Equilibrium) condition radiation has a Planck black body) distribution. Radiation energy density is given as U r,ν = 8πh c 3 ν 3, LTE), tr.1) e hν/kt
More information2. Stellar atmospheres: Structure
2. Stellar atmospheres: Structure 2.1. Assumptions Plane-parallel geometry Hydrostatic equilibrium, i.e. o no large-scale accelerations comparable to surface gravity o no dynamically significant mass loss
More information1. Runaway Greenhouse
1. Runaway Greenhouse Making the simple complicated is commonplace; making the complicated simple, awesomely simple, that s creativity. - Charlie Mingus (Jazz Musician) Initially we assume that the radiative
More informationAstro 305 Lecture Notes Wayne Hu
Astro 305 Lecture Notes Wayne Hu Set 1: Radiative Transfer Radiation Observables From an empiricist s point of view there are 4 observables for radiation Energy Flux Direction Color Polarization Energy
More informationEnergy transport: convection
Outline Introduction: Modern astronomy and the power of quantitative spectroscopy Basic assumptions for classic stellar atmospheres: geometry, hydrostatic equilibrium, conservation of momentum-mass-energy,
More informationModelling stellar atmospheres with full Zeeman treatment
1 / 16 Modelling stellar atmospheres with full Zeeman treatment Katharina M. Bischof, Martin J. Stift M. J. Stift s Supercomputing Group FWF project P16003 Institute f. Astronomy Vienna, Austria CP#AP
More informationSolution of the radiative transfer equation in NLTE stellar atmospheres
Solution of the radiative transfer equation in NLTE stellar atmospheres Jiří Kubát kubat@sunstel.asu.cas.cz Astronomický ústav AV ČR Ondřejov Non-LTE Line Formation for Trace Elements in Stellar Atmospheres,
More information1 3 µ2 p ν (τ ν ) = I ν. s ν (τ ν ) + (1) µ dq ν. dτ ν. and. = Q ν {( 1 µ 2) p ν (τ ν ) }, (2)
1. Plane-parallel atmospheres with polarized scattering. In the treatment by Chandrasehar (196, pp 38-5), the transfer equations are expressed in terms of I l and I r, the intensities in directions parallel
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 informationFundamental Stellar Parameters
Fundamental Stellar Parameters Radiative Transfer Specific Intensity, Radiative Flux and Stellar Luminosity Observed Flux, Emission and Absorption of Radiation Radiative Transfer Equation, Solution and
More informationPhysics Dec The Maxwell Velocity Distribution
Physics 301 7-Dec-2005 29-1 The Maxwell Velocity Distribution The beginning of chapter 14 covers some things we ve already discussed. Way back in lecture 6, we calculated the pressure for an ideal gas
More informationOpacity and Optical Depth
Opacity and Optical Depth Absorption dominated intensity change can be written as di λ = κ λ ρ I λ ds with κ λ the absorption coefficient, or opacity The initial intensity I λ 0 of a light beam will be
More informationTo remain at 0 K heat absorbed by the medium must be removed in the amount of. dq dx = q(l) q(0) dx. Q = [1 2E 3 (τ L )] σ(t T 4 2 ).
294 RADIATIVE HEAT TRANSFER 13.7 Two infinite, isothermal plates at temperatures T 1 and T 2 are separated by a cold, gray medium of optical thickness τ L = L (no scattering). (a) Calculate the radiative
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationStellar Atmospheres: Basic Processes and Equations
Stellar Atmospheres: Basic Processes and Equations Giovanni Catanzaro Abstract The content of this chapter is a very quick summary of key concepts that concern the interaction between photons created in
More informationModule 6: Implicit Runge-Kutta Methods Lecture 17: Derivation of Implicit Runge-Kutta Methods(Contd.) The Lecture Contains:
The Lecture Contains: We continue with the details about the derivation of the two stage implicit Runge- Kutta methods. A brief description of semi-explicit Runge-Kutta methods is also given. Finally,
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 informationChapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow
Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already
More informationRadiative transfer in planetary atmospheres
Chapter 9 Radiative transfer in planetary atmospheres Warning: This chapter is just a very rough introduction to planetary atmospheres. Since the discovery of planets around other stars than the sun, the
More informationMTH 3311 Test #1. order 3, linear. The highest order of derivative of y is 2. Furthermore, y and its derivatives are all raised to the
MTH 3311 Test #1 F 018 Pat Rossi Name Show CLEARLY how you arrive at your answers. 1. Classify the following according to order and linearity. If an equation is not linear, eplain why. (a) y + y y = 4
More informationDeriving the Asymptotic Telegrapher s Equation (P 1 ) Approximation for thermal Radiative transfer
Deriving the Asymptotic Telegrapher s Equation (P 1 ) Approximation for thermal Radiative transfer Shay I. Heizler Department of Physics, Bar-Ilan University, Ramat-Gan Table of Contents Introduction and
More informationarxiv: v1 [physics.comp-ph] 22 Jul 2010
Gaussian integration with rescaling of abscissas and weights arxiv:007.38v [physics.comp-ph] 22 Jul 200 A. Odrzywolek M. Smoluchowski Institute of Physics, Jagiellonian University, Cracov, Poland Abstract
More informationSection 11.5 and Problem Radiative Transfer. from. Astronomy Methods A Physical Approach to Astronomical Observations Pages , 377
Section 11.5 and Problem 11.51 Radiative Transfer from Astronomy Methods A Physical Approach to Astronomical Observations Pages 365-375, 377 Cambridge University Press 24 by Hale Bradt Hale Bradt 24 11.5
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationEXACT EQUATIONS AND INTEGRATING FACTORS
MAP- EXACT EQUATIONS AND INTEGRATING FACTORS First-order Differential Equations for Which We Can Find Eact Solutions Stu the patterns carefully. The first step of any solution is correct identification
More informationInverse Problems, Information Content and Stratospheric Aerosols
Inverse Problems, Information Content and Stratospheric Aerosols May 22, 214 Table of contents 1 Inverse Problems 2 Aerosols 3 OSIRIS 4 Information Content 5 Non-Linear Techniques 6 Results Inverse Problems
More informationAnalysis of Scattering of Radiation in a Plane-Parallel Atmosphere. Stephanie M. Carney ES 299r May 23, 2007
Analysis of Scattering of Radiation in a Plane-Parallel Atmosphere Stephanie M. Carney ES 299r May 23, 27 TABLE OF CONTENTS. INTRODUCTION... 2. DEFINITION OF PHYSICAL QUANTITIES... 3. DERIVATION OF EQUATION
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationBeer-Lambert (cont.)
The Beer-Lambert Law: Optical Depth Consider the following process: F(x) Absorbed flux df abs F(x + dx) Scattered flux df scat x x + dx The absorption or scattering of radiation by an optically active
More informationVII. Hydrodynamic theory of stellar winds
VII. Hydrodynamic theory of stellar winds observations winds exist everywhere in the HRD hydrodynamic theory needed to describe stellar atmospheres with winds Unified Model Atmospheres: - based on the
More informationLimb Darkening. Limb Darkening. Limb Darkening. Limb Darkening. Empirical Limb Darkening. Betelgeuse. At centre see hotter gas than at edges
Limb Darkening Sun Betelgeuse Limb Darkening Stars are both redder and dimmer at the edges Sun Limb Darkening Betelgeuse Limb Darkening Can also be understood in terms of temperature within the solar photosphere.
More informationSupplement: Statistical Physics
Supplement: Statistical Physics Fitting in a Box. Counting momentum states with momentum q and de Broglie wavelength λ = h q = 2π h q In a discrete volume L 3 there is a discrete set of states that satisfy
More informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationRadiative Transfer Plane-Parallel Frequency-Dependent
4 Radiative Transfer Plane-Parallel Frequency-Dependent variables I ν J ν H ν K ν in astronomy 4. Basic Equations After Fukue, J. 2, PASJ, 63, in press We here assume the followings: i) The disk is steady
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationB.V. Gudiksen. 1. Introduction. Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della
Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della À Ø Ò Ø ËÓÐ Ö ÓÖÓÒ B.V. Gudiksen Institute of Theoretical Astrophysics, University of Oslo, Norway e-mail:boris@astro.uio.no Abstract. The heating mechanism
More informationProfessor Terje Haukaas University of British Columbia, Vancouver Size Effect Models. = 1 F(τ ) (1) = ( 1 F(τ )) n (2) p s
Size Effect Models Weibull Theory This theory is based on a weaest lin approach to material strength, and essentially simplifies the problem to the calibration of two parameters m and. To understand the
More informationEART164: PLANETARY ATMOSPHERES
EART164: PLANETARY ATMOSPHERES Francis Nimmo Last Week Radiative Transfer Black body radiation, Planck function, Wien s law Absorption, emission, opacity, optical depth Intensity, flux Radiative diffusion,
More informationCHAPTER 26. Radiative Transfer
CHAPTER 26 Radiative Transfer Consider an incoming signal of specific intensity I ν,0 passing through a cloud (i.e., any gaseous region). As the radiation transits a small path length dr through the cloud,
More information42. Change of Variables: The Jacobian
. Change of Variables: The Jacobian It is common to change the variable(s) of integration, the main goal being to rewrite a complicated integrand into a simpler equivalent form. However, in doing so, the
More informationPreface to the Second Edition. Preface to the First Edition
Contents Preface to the Second Edition Preface to the First Edition iii v 1 Introduction 1 1.1 Relevance for Climate and Weather........... 1 1.1.1 Solar Radiation.................. 2 1.1.2 Thermal Infrared
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationChapter 2: Numeric, Cell, and Structure Arrays
Chapter 2: Numeric, Cell, and Structure Arrays Topics Covered: Vectors Definition Addition Multiplication Scalar, Dot, Cross Matrices Row, Column, Square Transpose Addition Multiplication Scalar-Matrix,
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationNumerical Methods. Equations and Partial Fractions. Jaesung Lee
Numerical Methods Equations and Partial Fractions Jaesung Lee Solving linear equations Solving linear equations Introduction Many problems in engineering reduce to the solution of an equation or a set
More informationCMFGEN. Associate every atomic species with a given identification number
CMFGEN Associate every atomic species with a given identification number For each species, set: atomic mass atomic number solar abundance species abbreviations species identification Read in MODEL_SPEC
More informationNumerical Methods for Engineers, Second edition: Chapter 1 Errata
Numerical Methods for Engineers, Second edition: Chapter 1 Errata 1. p.2 first line, remove the Free Software Foundation at 2. p.2 sixth line of the first proper paragraph, fe95.res should be replaced
More informationStellar Atmospheres. University of Denver, Department of Physics and Astronomy. Physics 2052 Stellar Physics, Winter 2008.
Stellar Atmospheres University of Denver, Department of Physics and Astronomy Physics 2052 Stellar Physics, Winter 2008 By Kathy Geise Introduction A star does not have a solid surface, so the definition
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationNext we shall remove the only unphysical approximation we made in this CASE, which is that we ignored delayed neutrons.
71 Several observations are in order: 1. Because the terms in the [ ] are defined in terms of n, which depends on t, they will in general depend on t. (But not on position or energy, because position and
More informationStars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11)
Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11) Kenneth Wood, Room 316 kw25@st-andrews.ac.uk http://www-star.st-and.ac.uk/~kw25 What is a Stellar Atmosphere? Transition from dense
More informationOutline. Today we will learn what is thermal radiation
Thermal Radiation & Outline Today we will learn what is thermal radiation Laws Laws of of themodynamics themodynamics Radiative Radiative Diffusion Diffusion Equation Equation Thermal Thermal Equilibrium
More informationNon-LTE models for synthetic spectra of Type Ia supernovae
A&A 569, A62 (24) DOI:.5/4-636/2322265 c ESO 24 Astronomy & Astrophysics Non-LTE models for synthetic spectra of Type Ia supernovae IV. A modified Feautrier scheme for opacity-sampled pseudo-continua at
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationMTF071 Computational Fluid Dynamics of Turbulent
CHALMERS TEKNISKA HÖGSKOLA Termo- och Fluiddynamik 42 96 Göteborg MTF07 Computational Fluid Dynamics of Turbulent Flow http://wwwtfdchalmersse/gr-kurs/mtf07 Task K2 Lars Davidson 2003-02-03 In Task K you
More informationStellar atmospheres: an overview
Stellar atmospheres: an overview Core M = 2x10 33 g R = 7x10 10 cm 50 M o 20 R o L = 4x10 33 erg/s 10 6 L o 10 4 (PN) 10 6 (HII) 10 12 (QSO) L o Photosphere Envelope Chromosphere/Corona R = 200 km ~ 3x10
More informationRadiative Equilibrium Models. Solar radiation reflected by the earth back to space. Solar radiation absorbed by the earth
I. The arth as a Whole (Atmosphere and Surface Treated as One Layer) Longwave infrared (LWIR) radiation earth to space by the earth back to space Incoming solar radiation Top of the Solar radiation absorbed
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationTHREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS
THREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS Carsten Beckert 1. Introduction To calculate the neutron transport in a reactor, it is often
More informationSYDE 112, LECTURE 7: Integration by Parts
SYDE 112, LECTURE 7: Integration by Parts 1 Integration By Parts Consider trying to take the integral of xe x dx. We could try to find a substitution but would quickly grow frustrated there is no substitution
More informationAtmospheric Sciences 321. Science of Climate. Lecture 6: Radiation Transfer
Atmospheric Sciences 321 Science of Climate Lecture 6: Radiation Transfer Community Business Check the assignments Moving on to Chapter 3 of book HW #2 due next Wednesday Brief quiz at the end of class
More informationRecap Lecture + Thomson Scattering. Thermal radiation Blackbody radiation Bremsstrahlung radiation
Recap Lecture + Thomson Scattering Thermal radiation Blackbody radiation Bremsstrahlung radiation LECTURE 1: Constancy of Brightness in Free Space We use now energy conservation: de=i ν 1 da1 d Ω1 dt d
More informationThe Solar Neutrino Problem. There are 6 major and 2 minor neutrino producing reactions in the sun. The major reactions are
The Solar Neutrino Problem There are 6 major and 2 minor neutrino producing reactions in the sun. The major reactions are 1 H + 1 H 2 H + e + + ν e (PP I) 7 Be + e 7 Li + ν e + γ (PP II) 8 B 8 Be + e +
More informationNIMEQ: MHD Equilibrium Solver for NIMROD
NIMEQ: MHD Equilibrium Solver for NIMOD E.C.Howell, C..Sovinec University of Wisconsin-Madison 5 th Annual Meeting of Division of Plasma Physics Dallas, Texas, Nov. 17-Nov. 1,8 1 Abstract A Grad-Shafranov
More information(Refer Slide Time: 01:11) So, in this particular module series we will focus on finite difference methods.
Computational Electromagnetics and Applications Professor Krish Sankaran Indian Institute of Technology, Bombay Lecture 01 Finite Difference Methods - 1 Good Morning! So, welcome to the new lecture series
More informationComputational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras
Computational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras Module No. # 05 Lecture No. # 24 Gauss-Jordan method L U decomposition method
More informationAccretion disks. AGN-7:HR-2007 p. 1. AGN-7:HR-2007 p. 2
Accretion disks AGN-7:HR-2007 p. 1 AGN-7:HR-2007 p. 2 1 Quantitative overview Gas orbits in nearly circular fashion Each gas element has a small inward motion due to viscous torques, resulting in an outward
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationOn the General Solution of Initial Value Problems of Ordinary Differential Equations Using the Method of Iterated Integrals. 1
On the General Solution of Initial Value Problems of Ordinary Differential Equations Using the Method of Iterated Integrals. 1 Ahsan Amin ahsanamin2999@gmail.com This First version January 2017 First Preliminary
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More information(a) x cos 3x dx We apply integration by parts. Take u = x, so that dv = cos 3x dx, v = 1 sin 3x, du = dx. Thus
Math 128 Midterm Examination 2 October 21, 28 Name 6 problems, 112 (oops) points. Instructions: Show all work partial credit will be given, and Answers without work are worth credit without points. You
More informationPhysics 607 Final Exam
Physics 607 Final Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationEFFECTS OF THE TEMPERATURE DISCRETIZATION ON NUMERICAL METHODS FOR THERMAL RADIATION TRANSPORT
International Conference on Mathematics, Computational Methods & Reactor Physics M&C 2009 Saratoga Springs, New York, May 3 7, 2009, on CD-ROM, American Nuclear Society, LaGrange Park, IL 2009 EFFECTS
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationCS 542G: Robustifying Newton, Constraints, Nonlinear Least Squares
CS 542G: Robustifying Newton, Constraints, Nonlinear Least Squares Robert Bridson October 29, 2008 1 Hessian Problems in Newton Last time we fixed one of plain Newton s problems by introducing line search
More informationTopics ASTR 3730: Fall 2003
Topics Qualitative questions: might cover any of the main topics (since 2nd midterm: star formation, extrasolar planets, supernovae / neutron stars, black holes). Quantitative questions: worthwhile to
More informationarxiv: v2 [astro-ph.ep] 7 Nov 2014
Astronomy & Astrophysics manuscript no. Parmentier24c c ESO 24 November, 24 A non-grey analytical model for irradiated atmospheres. II: Analytical vs. numerical solutions Vivien Parmentier,2,3, Tristan
More information5. Solving the Bellman Equation
5. Solving the Bellman Equation In the next two lectures, we will look at several methods to solve Bellman s Equation (BE) for the stochastic shortest path problem: Value Iteration, Policy Iteration and
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationNotes on arithmetic. 1. Representation in base B
Notes on arithmetic The Babylonians that is to say, the people that inhabited what is now southern Iraq for reasons not entirely clear to us, ued base 60 in scientific calculation. This offers us an excuse
More informationRegression. Oscar García
Regression Oscar García Regression methods are fundamental in Forest Mensuration For a more concise and general presentation, we shall first review some matrix concepts 1 Matrices An order n m matrix is
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationA high accuracy approximation for half - space problems with anisotropic scattering
IOSR Journal of Applied Physics (IOSR-JAP e-iss: 78-486. Volume 4, Issue (May. - Jun., PP 7-46 A high accuracy approximation for half - space problems with anisotropic scattering M. H. Haggag,, H.F. Mashali
More information