Methods for solving the radiative transfer equation with multiple scattering. Part 3: Exact methods: Discrete-ordinate and Adding.

Similar documents
Methods for solving the radiative transfer equation. Part 3: Discreteordinate. 1. Discrete-ordinate method for the case of isotropic scattering.

C.11 Bang-bang Control

1. Inclusion of the surface reflection into the radiative transfer equation.

CS 2750 Machine Learning Lecture 8. Linear regression. Supervised learning. a set of n examples

Load balancing by MPLS in differentiated services networks

THE THEORY OF AN INSCRIBABLE QUADRILATERAL AND A CIRCLE THAT FORMS PASCAL POINTS

M2S1 - EXERCISES 8: SOLUTIONS

A Note on the Almost Sure Central Limit Theorem in the Joint Version for the Maxima and Partial Sums of Certain Stationary Gaussian Sequences *

MATRIX ANALYSIS OF ANCHORED STRUCTURES

SOLVING FUZZY SOLID TRANSPORTATION PROBLEM BASED ON EXTENSION PRINCIPLE WITH INTERVAL BUDGET CONSTRAINT

ELEMENTS OF NUMBER THEORY. In the following we will use mainly integers and positive integers. - the set of integers - the set of positive integers

7.0 Equality Contraints: Lagrange Multipliers

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

Solutions for HW4. x k n+1. k! n(n + 1) (n + k 1) =.

Hamilton s principle for non-holonomic systems

Generalized Linear Regression with Regularization

UNIT 7 RANK CORRELATION

Convergence Analysis of Piecewise Polynomial Collocation Methods for System of Weakly Singular Volterra Integral Equations of The First Kind

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

are positive, and the pair A, B is controllable. The uncertainty in is introduced to model control failures.

ROOT-LOCUS ANALYSIS. Lecture 11: Root Locus Plot. Consider a general feedback control system with a variable gain K. Y ( s ) ( ) K

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

Non-degenerate Perturbation Theory

LINEARLY CONSTRAINED MINIMIZATION BY USING NEWTON S METHOD

Basic Concepts in Numerical Analysis November 6, 2017

Abstract. 1. Introduction

Some Different Perspectives on Linear Least Squares

UNIT 6 CORRELATION COEFFICIENT

A stopping criterion for Richardson s extrapolation scheme. under finite digit arithmetic.

KURODA S METHOD FOR CONSTRUCTING CONSISTENT INPUT-OUTPUT DATA SETS. Peter J. Wilcoxen. Impact Research Centre, University of Melbourne.

ENGI 4421 Propagation of Error Page 8-01

Maximum Walk Entropy Implies Walk Regularity

Functions of Random Variables

Solutions to problem set ); (, ) (

3.1 Introduction to Multinomial Logit and Probit

Band structure calculations

The Mathematics of Portfolio Theory

Stationary states of atoms and molecules

Johns Hopkins University Department of Biostatistics Math Review for Introductory Courses

Numerical Analysis Formulae Booklet

Johns Hopkins University Department of Biostatistics Math Review for Introductory Courses

General Method for Calculating Chemical Equilibrium Composition

Chapter 3. Differentiation 3.3 Differentiation Rules

= 2. Statistic - function that doesn't depend on any of the known parameters; examples:

STK4011 and STK9011 Autumn 2016

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

CS 2750 Machine Learning. Lecture 7. Linear regression. CS 2750 Machine Learning. Linear regression. is a linear combination of input components x

ρ < 1 be five real numbers. The

CH E 374 Computational Methods in Engineering Fall 2007

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming

X ε ) = 0, or equivalently, lim

PART ONE. Solutions to Exercises

A Primer on Summation Notation George H Olson, Ph. D. Doctoral Program in Educational Leadership Appalachian State University Spring 2010

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

Coherent Potential Approximation

THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 9, Number 3/2008, pp

On the Pfaff problem

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

How to break tetrahedral symmetry

Given a table of data poins of an unknown or complicated function f : we want to find a (simpler) function p s.t. px (

Chapter 4 Multiple Random Variables

Chapter 3. Differentiation 3.2 Differentiation Rules for Polynomials, Exponentials, Products and Quotients

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

Solutions to Odd-Numbered End-of-Chapter Exercises: Chapter 17

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION

Chapter 5 Properties of a Random Sample

Long blade vibration model for turbine-generator shafts torsional vibration analysis

A Family of Generalized Stirling Numbers of the First Kind

STAT 400 Homework 09 Spring 2018 Dalpiaz UIUC Due: Friday, April 6, 2:00 PM

3D Reconstruction from Image Pairs. Reconstruction from Multiple Views. Computing Scene Point from Two Matching Image Points

9 U-STATISTICS. Eh =(m!) 1 Eh(X (1),..., X (m ) ) i.i.d

Optimality of Strategies for Collapsing Expanded Random Variables In a Simple Random Sample Ed Stanek

1. Linear second-order circuits

Plate Bending Analysis by Two-dimensional Non-linear Partial Differential Equations

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

MOLECULAR VIBRATIONS

( x) min. Nonlinear optimization problem without constraints NPP: then. Global minimum of the function f(x)

The theoretical background of

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

Principles of multiple scattering in the atmosphere. Radiative transfer equation with scattering for solar radiation in a plane-parallel atmosphere.

April 1980 TR/96. Extrapolation techniques for first order hyperbolic partial differential equations. E.H. Twizell

1 Onto functions and bijections Applications to Counting

Third handout: On the Gini Index

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA

DATA DOMAIN DATA DOMAIN

Econometric Methods. Review of Estimation

Lecture 8 IEEE DCF Performance

VIII Dynamics of Systems of Particles

Lecture 07: Poles and Zeros

Lecture 14. Review for Exam 1.

Landé interval rule (assignment!) l examples

Mu Sequences/Series Solutions National Convention 2014

Linear regression (cont) Logistic regression

. The set of these sums. be a partition of [ ab, ]. Consider the sum f( x) f( x 1)

ESS Line Fitting

Example of CLT for Symmetric Laminate with Mechanical Loading

Overview. Review Superposition Solution. Review Superposition. Review x and y Swap. Review General Superposition

A Penalty Function Algorithm with Objective Parameters and Constraint Penalty Parameter for Multi-Objective Programming

Week 10 Spring Lecture 19. Estimation of Large Covariance Matrices: Upper bound Observe. is contained in the following parameter space,

Transcription:

Lecture. Methos for sovg the raatve trasfer equato wth utpe scatterg. art 3: Exact ethos: screte-orate a Ag. Obectves:. screte-orate etho for the case of sotropc scatterg.. Geerazato of the screte-orate etho for hoogeeous atosphere. 3. uerca peetato of the screte-orate etho: O 4. rcpes of varace. 5. Ag etho. equre reag: L: 6. 6.3.- 6.3.4 6.4 cchazz... Yag et a. BA: A research a teachg software too for ae-parae raatve trasfer the earth's atosphere. Buet of the Aerca Meteoroogca ocety 79-4 998. http:www.cess.ucsb.euesrgba.ht Avace reag: hoas G.E. a K. taes aatve trasfer the atosphere a ocea Chapter 8.- 8.. Lu Q. a. Weg Avace oubg Ag Metho for aatve rasfer aetary Atospheres. Joura of Atospherc ceces 63 3459 3465 6. taes K.. say W. Wscobe a K. Jayaweera uercay stabe agorth for screteorate-etho raatve trasfer utpe scatterg a ettg ayere ea." App. Opt. 7 5-59 998.. screte-orate etho for the case of sotropc scatterg. OE: A screte-orate etho has bee eveope by Charasekhar the 95s Charasekhar. aatve trasfer 96 over ubcatos. eca the raatve trasfer equato Lecture 8 for azuthay epeet ffuse testy: ' ' ' exp 4

or sotropc scatterg the scatterg phase fucto s. Hece we have ' ' exp 4 Let s appy the Gaussa quaratures to repace the tegra Eq.[.] [.] a exp [.] 4 hoogeeous part where - ters a a are the Gaussa weghts costats a are quarature ages or pots. Eq.[.] s a syste of hoogeeous ffereta equatos: outo of Eq.[.] geera souto partcuar souto where the geera souto s a souto of the hoogeeous part of the Eq.[.] eotg the geera souto of Eq.[.] ca be fou as g exp k [.3] sertg Eq.[.3] to Eq.[.] we obta g k a g [.4] We ca f g the for g L k where L s a costat to be etere. ubsttutg ths expresso for g Eq.[.4] we have a k Eq.[.5] gves soutos for -k. a k [.5] hus geera souto s where L are costats. L exp k [.6] k

he partcuar souto ca be fou as h exp [.7] 4 where h are costats. sertg Eq.[.7] to Eq.[.] we have h a h [.8] ro Eq.[.8] h s fou as h γ where γ s etere fro γ { } [.9] a Ag the geera souto Eq.[.6] a the partcuar souto Eq.[.7] we have L γ exp k exp [.] k 4 where L are costats to be etere fro the bouary cotos. H-fucto has bee trouce by Charasekhar as H... k [.] Expressg γ the H-fucto Eq.[.] becoes L H H exp k exp [.] k 4 Eq.[.] gves a spe souto for the se-fte sotropc atosphere see L:6.. H H 4 [.3] 3

4. Geerazato of the screte-orate etho for a hoogeeous atosphere. Let s coser the atosphere wth o-sotropc scatterg. We ca expa the ffuse testy the cose seres cos ϕ ϕ ϕ o we ee to sove exp 4 4 ϖ ϖ δ o he geera souto ay be wrtte as k L exp φ φ k L are coeffcets to be etere. he partcuar souto ay be wrtte as exp p Z Where Z s the foowg fucto 4 H H Z ζ ϖ he copete souto of the raatve trasfer s Z k L exp exp φ [.4] - Let s geeraze the copete souto Eq.[.4] of the raatve trasfer for the hoogeeous atosphere. he atosphere ca be ve to the hoogeeous

ayers each s characterze by a sge scatterg abeo phase fucto a optca epth. OE: f a atospherc ayer has gases aerosos aor cous oe ees to cacuate the effectve optca propertes of ths ayer. or -th ayer we ca wrte the souto usg Eq.[.4]. o spfy otatos et s coser the azutha epeet case.e. so we have L φ exp k Z exp [.5] ow we ee to atch the bouary a cotuty cotos betwee ayers. At the top of the atosphere OA: o owwar ffuse testy [.6] At the ayer s bouary: upwar a owwar testes ust be cotuous [.7] At the botto of the atosphere assug the Laerta surface: rsur [ exp ] Eqs.[.6]-[.8] prove ecessary equatos to f the ukow coeffcets. [.8] 3. uerca peetato of the screte-orate etho: O O s a OA uerca coe base o the screte-orate etho eveope by taes Wscobe et a. O s opey avaabe a has a goo user-gue. oe features: O appes to the hoogeeous othothera pae-parae atosphere. A user ay set-up ay ubers of the pae-parae ayers. 3 Each ayer ust be characterze by the effectve optca epth sge scatterg abeo a asyetry paraeter f the Heyey-Greeste phase fucto s use. 4 A user ay use ay phase fucto by provg the Legere poyoa expaso coeffcets. 5

5 A user seects a uber of streas keepg that the coputato te vares as 3. 6 A key probe s to obta a souto for fuxes for strogy forwar-peake scatterg. 7 O aows prectg the testy as a fucto of the recto a posto at ay pot the atosphere.e. ot oy at the bouares of the ayers. O s corporate to the BA raatve trasfer coe. BA ata Barbara O Atospherc aatve rasfer s a OA coputer coe: http:www.cess.ucsb.euesrgba.ht see troucto to BA 4. rcpes of varace. eca the eftos of refecto a trassso of a ayer trouce Lecture 9. f the soar fux s cet o a ayer of optca epth : ϕ ϕ ϕ Or the geera case: ϕ ϕ * ϕ t r ϕ ϕ ϕ ϕ ϕ c * ϕ ϕ ϕ ϕ ϕ c he prcpe of varace for the se-fte atosphere Abartzua 94: the ffuse refecte testy caot be chage f a ayer of fte optca epth havg the sae optca propertes as those of the orga ayer s ae see L: 6.3.. 6

he prcpes of varace for a fte atosphere Charasekhar 95: he refecte upwar testy at ay gve optca epth resuts fro the refecto of a the atteuate soar fux exp a b the owwar ffuse testy at the eve : exp * [.9] ' * * he ffusey trastte owwar testy at the eve resuts fro a the trassso of cet soar fux a b the refecto of the upwar ffuse testy above the eve : [.] ' 7

8 3 he refecte upwar testy at the top of the fte atosphere s equvaet to a the refecto of soar fux pus b the rect a ffuse trassso of the upwar ffuse testy above the eve : exp [.] 4 he ffusey trastte owwar testy at the botto of the fte atosphere s equvaet to a the trassso of the atteuate soar fux at the eve pus b the rect a ffuse trassso of the owwar ffuse testy at the eve fro above: * exp * * exp * [.] ' ' o * * ' '

5. Ag etho. Ag etho s a exact etho for sovg the raatve trasfer equato wth utpe scatterg. t uses geoetrca ray-tracg approach a the refecto a trassso of each vua atospherc ayer. trategy: kowg the refecto a trassso of two vua ayers the refecto a trassso of the cobe ayer ay be obtae by cacuatg the successve refectos a trasssos betwee these two ayers. OE: f optca epths of these two ayers are equae ths etho s referre to as the oubg-ag etho. Coser two ayers wth refecto a a tota rect pus ffuse trassso a fuctos respectvey. Let s eote the cobe refecto a tota trassso fuctos by a a cobe refecto a tota trassso fuctos betwee ayers a by a respectvey. 9

he cobe refecto fucto s...] [... [.3] OE: Eq.[.3] we use that x x he cobe tota trassso fucto s...] [... [.4] he cobe refecto fucto betwee ayers a :...] [... [.5] he cobe tota trassso fucto betwee ayers a :...] [... [.6] ro Eqs.[.3]-[.6] we f that ; ; [.7] Let s trouce sg that exp fro Eqs.[.6]-[.7] we f exp exp exp exp [.8]

exp exp exp exp exp δ [.9] hus we ca wrte a syste of teratve equatos for the coputato of ffuse trassso a refecto for the two ayers the for: Q Q Q exp exp exp exp exp [.3] OE: Eq.[.3] the prouct of two fuctos pes tegrato over the approprate age so that a utpe-scatterg cotrbutos are cue. or stace uerca proceure of the ag etho: As the startg pot oe ay cacuate the refecto a trassso fuctos of a ta ayer of very sa optca epth e.g. -8 that the sge scatterg approxato s appcabe. he usg Eq.[.3] oe coputes the refecto a trassso fuctos of the ayer of. 3 sg Eq.[.3] oe repeats the cacuatos ag the ayers ut a esrabe optca epth s acheve.