Processing math: 71%

Amateur Astronomie

Benutzer-Werkzeuge

Webseiten-Werkzeuge

Zuletzt angesehen:
demo_limbdarkening

Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen angezeigt.

Link zu dieser Vergleichsansicht

Beide Seiten der vorigen RevisionVorhergehende Überarbeitung
Nächste Überarbeitung		Vorhergehende Überarbeitung
demo_limbdarkening [2022/03/29 09:13] – torsten.roehldemo_limbdarkening [2022/04/04 17:44] (aktuell) – torsten.roehl
Zeile 1: Zeile 1:
-====== Lichtkurve ohne Randverdunkelung ====== +====== Lichtkurve mit Randverdunkelung (engl. limb darkening) ====== 
-// Um die Effekte der Randverdunkelung (engl. limb darkening) besser verstehen zu können, wird zuerst die Lichtkurve __ohne__ Randverdunkelung in Abhängigkeit der Planetengrö0e und des Impaktparameters **b** untersucht.+// Unter Randverdunkelung (engl. limb darkening) versteht man, das die Helligkeit eines Sterns nicht gleichmäßig über die Fläche verteilt ist. An den Rändern ist der Stern dunkler (kühler) in der Mitte heller (heißer). 
 +Direkt beobachten kann man das allerdings nur an der Sonne, bei den anderen Sternen kann man nur durch physikalische Modelle und Analysen des Sternenlichtes (Photometrie, Spektroskopie) indirekt auf die jeweilige  
 +Randverdunkelung schließen. Da aber Sterne nun mal nicht gleichmäßig hell über ihre Oberfläche leuchten, bedeutet das, das dies bei der Auswertung der Lichtkurve eines Exoplaneten berücksichtigt werden muss. Es zeigt sich weiterhin, dass diese Effekt auch von der beobachteten Wellenlängen (also des verwendeten Filters) abhängig ist. Die Modellierung der Lichtkurve (z.b. mit AstroimageJ) berücksichtigt die Randverdunkelung durch die Annahme eines Quadratischen Modells. Dieses dieses Modell (ANALYTIC LIGHT CURVES FOR PLANETARY TRANSIT SEARCHES [[https://ui.adsabs.harvard.edu/abs/2002ApJ...580L.171M/abstract| Mandel  and Algol (2002)]]) liegt der Animation zur Grunde. 
 +// 
 + 
 + 
 +<sortable> 
 +| {{ :sun_limb.png?120 | }} **Randverdunkelung** (Foto ohne Sonnenflecken und Venus) | {{ :sun_orginal.jpg?120 |}} Foto der Sonne im sichtbaren Licht (Venustransit 2012) https://en.wikipedia.org/wiki/Limb_darkening   |  {{ :sun_uniform.png?120 |}} __Ohne__ Randverdunkelung (uniforme Helligkeit)| 
 +</sortable> 
 + 
 +// 
 +Man kann aber auch umgekehrt vorgehen, und mit Hilfe der experimentell gewonnen Lichtkurven die Physik des Zentralgestirns studieren. 
 +Beispielsweise lassen sich  die physikalische Prozesse die den theoretischen  Modellen zugrunde liegen untersuchen, indem  die gewonnen realen Daten der Lichtkurven mit den theoretischen Modellen verglichen werden.
 // //
 <WRAP center  box 100%> <WRAP center  box 100%>
Zeile 7: Zeile 19:
 </WRAP> </WRAP>
  
-===== Lichkurve ohne Randverdunkelung  ===== +===== Lichkurve mit Randverdunkelung  ===== 
-Bei der Lichtkurve ohne Randverdunkelung, gehen wir von eine gleichmäßigen (uniformen) Beleuchtung, sowohl des Himmelshintergrund als auch des Sterns aus. Die Animationen berechnet die exakten Schnittfläche zwischen Exoplaneten und Stern, d.h. wir beziehen uns auf eine 2-dimensionale Geometrie und vernachlässigen z.B., dass Stern und Planet Kugelförmig sind.  Die Skalierung der Lichtkurve ist relativ, um möglichst viele Details darstellen zu können.+ 
 <PRELOAD> <PRELOAD>
 +
 ##########################################################################################################  ########################################################################################################## 
  /lib/scripts/jqplot/jquery.jqplot.js  /lib/scripts/jqplot/jquery.jqplot.js
Zeile 24: Zeile 38:
 # #
 /lib/scripts/src/header.js  /lib/scripts/src/header.js 
-/lib/scripts/src/pointCharge.js+/lib/scripts/src/mandelAlgol.js  
 ############################################################################################################ ############################################################################################################
 </PRELOAD> </PRELOAD>
Zeile 35: Zeile 50:
 var bText = true; var bText = true;
 var bImpact = true; var bImpact = true;
 +var bPlanet = true;
 +
 var xmin   = -7.5; var xmin   = -7.5;
 var xmax   =  7.5; var xmax   =  7.5;
 // slider // slider
-var planet = 0.1;+var planet = 0.12;
 var impact = 0.0;  var impact = 0.0; 
  
-var unit = 12.5; + 
-function onInit() {  + 
-writeMainHeader();+function onInit() {  // first Animation context 
 + 
 + 
  
 /////////////// BOARD  /////////////// BOARD 
Zeile 60: Zeile 80:
  
 //   //  
- +pathBild = 'http://wiki.linux-astronomie.de/lib/exe/fetch.php?media=sun_limb.png' 
 +brd.create('image', [pathBild,          [ -3,2], [6,6]] ,      {visible: true, fixed: true} );
  
 //////////// STAR  //////////// STAR 
Zeile 70: Zeile 90:
             size:75,             size:75,
             fixed: true,             fixed: true,
 +            visible: false,
             layer:8,             layer:8,
             fillColor: 'yellow',                       fillColor: 'yellow',          
Zeile 96: Zeile 117:
  
 function showImpact0(){  function showImpact0(){ 
-   if (bText && bImpact && impact <= 0.3) return true;+   if (bText && bImpact && impact <= 0.4) return true; 
 +   return false; 
 +} 
 +function showImpact04(){  
 +   if (bText && bImpact && impact > 0.4) return true;
    return false;    return false;
 } }
Zeile 107: Zeile 132:
  
  
 +brd.create('text',[-5.5, function() {return 5-(5-impactPosition())/2;},"b"],{ visible: function() { return showImpact04();},needsRegularUpdate:true});
 +
 +
 +brd.create('arrow', [ [-6,function() {return 5-(5-impactPosition())/2;}],[-6,5]] ,
 +{visible: function() {return showImpact04();},strokeColor:'#000'});
 +
 +brd.create('arrow', [ [-6,function() {return 5-(5-impactPosition())/2;} ],[-6,function(){ return impactPosition()}]] 
 +,{visible: function() {return showImpact04();},strokeColor:'#000', stokeWidth: 1});
 +
 +/////////// R_p/ R_Star
 +
 +function showPlanet(){ 
 +   if (bText && bPlanet) return true;
 +   return false;
 +}
 +
 +function planetText(){ return 'R_p/R_* ≈ ' + planet.toFixed(2) + ' ' ;   }
 +function planetText1(){ 
 +  
 +if (planet <= 0.01 ){
 +    return 'R_p/R_* ≈ R_{ERDE} ' ; 
 +}
 +if (planet > 0.01 && planet < 0.10 ){
 +    return 'R_p/R_* > R_{ERDE} ' ; 
 +}
 +if (planet == 0.10 ){
 +    return 'R_p/R_* ≈ R_{JUPITER} ' ; 
 +}
 +if (planet > 0.10 && planet != 0.12){
 +    return 'R_p/R_* > R_{JUPITER} ' ; 
 +}
 +if (planet == 0.12 ){
 +    return 'R_p/R_* ≈ R_{HD2099458b} ' ; 
 +}
 +}
 +
 +
 +
 +brd.create('text',[5,-5.5, function(){ return planetText();} ],{ visible: function() { return showPlanet();},needsRegularUpdate:true});
 +
 +
 +
 +brd.create('text',[5,-7, function(){ return planetText1();} ],{ visible: function() { return showPlanet();},needsRegularUpdate:true});
  
 // //
Zeile 131: Zeile 199:
 ////////// //////////
  
-function areaOfIntersection(x0, y0, r0, x1, y1, r1) 
-{ 
-    var rr0 = r0 * r0; 
-    var rr1 = r1 * r1; 
-    var d = Math.sqrt((x1 - x0) * (x1 - x0) + (y1 - y0) * (y1 - y0)); 
  
-    // Circles do not overlap +// var uniformDip = brd.create('functiongraph', 
-    if (d > r1 + r0) +//                       [function(x){ 
-    { +// 
-    return 0; +// 
-    }+//  x0 =0; y0 = 0; r0 = 75/2; 
 +// 
 +// 
 +//  x1 = 12.5*x; 
 +//  r1 = r0*planet; 
 +// 
 +//  scale = 300; 
 +//  shift = 3; 
 +//  
 +// // var ppp = hallo(planet); 
 +//  res = 160*uniformLD(0.35*x,impact,planet)-2;   
 +// // console.log("x=",0.2*x,"y=",res) 
 +// 
 +//  return res; 
 +// 
 +//       }, xmin, xmax ] ,{strokeWidth:3,visible:false}); 
 +////////////////////////////////////////////////////// 
 +var dipNone = brd.create('functiongraph', [function(x){
  
-    // Circle1 is completely inside circle0 +  scale = 300; 
-    else if (d <= Math.abs(r0 - r1) && r0 >= r1) +  shift = 2; 
-    { +  g1=0; 
-    // Return area of circle1 +  g2=0; 
-    return Math.PI * rr1; +  res = scale*quadraticLD(0.35*x,impact,g1,g2,planet)-shift-scale;   
-    }+ /// console.log("x=", x, "res=",res); 
 +  return res;
  
-    // Circle0 is completely inside circle1 +       }, xmin, xmax ] ,{strokeWidth:1.2,strokeColor:"#000000"}); 
-    else if (d <= Math.abs(r0 - r1) && r0 < r1) +////////////////////////////////////////////////////// 
-    { +var dipU = brd.create('functiongraph',  [function(x){
-    // Return area of circle0 +
-    return Math.PI * rr0; +
-    }+
  
-    // Circles partially overlap +  scale = 250; 
-    else +  shift = 2; 
-    { +  g1=0.67462231; 
-    var phi = (Math.acos((rr0 + (d * d) - rr1) / (2 * r0 * d))) * 2; +  g2= 0.16275864;  
-    var theta = (Math.acos((rr1 + (d * d) - rr0) / (2 * r1 * d))) * 2; +  res = scale*quadraticLD(0.35*x,impact,g1,g2,planet)-shift-scale;   
-    var area1 = 0.5 * theta * rr1 - 0.5 * rr1 * Math.sin(theta); +  
-    var area2 = 0.5 * phi * rr0 - 0.5 * rr0 * Math.sin(phi);+  return res;
  
-    // Return area of intersection +       }, xmin, xmax ] ,{strokeWidth:1.2,strokeColor:"#BF00FF"}); 
-    return area1 + area2; +////////////////////////////////////////////////////// 
-    } +////////////////////////////////////////////////////// 
-}+var dipB = brd.create('functiongraph',   [function(x){
  
 +  scale = 250;
 +  shift = 2;
 +  g1=0.57939128;
 +  g2=0.21110543; 
 +  res = scale*quadraticLD(0.35*x,impact,g1,g2,planet)-shift-scale;  
 +  
 +  return res;
  
 +       }, xmin, xmax ] ,{strokeWidth:1.2,strokeColor:"#0080FF"});
 +//////////////////////////////////////////////////////
 +//////////////////////////////////////////////////////
 +var dipV = brd.create('functiongraph',   [function(x){
  
-var uniformDip = brd.create('functiongraph', +  scale = 250; 
-                       [function(x){+  shift = 2; 
 +  g1= 0.40433600; 
 +  g2=0.29138968;  
 +  res = scale*quadraticLD(0.35*x,impact,g1,g2,planet)-shift-scale;  
  
 +  return res;
  
-  x0 =0; y0 = 0; r0 = 75/2;+}, xmin, xmax ] ,{strokeWidth:1.2, strokeColor:"#04B404",visible: true}); 
 +////////////////////////////////////////////////////// 
 +////////////////////////////////////////////////////// 
 +var dipR = brd.create('functiongraph',  [function(x){
  
-  y1 = impact*37;  +  scale = 250; 
-  x1 = 12.5*x; +  shift = 2; 
-  r1 = r0*planet; +  g1=0.31786107; 
- +  g2= 0.30489961;  
-  scale =-0.04; +  res = scale*quadraticLD(0.35*x,impact,g1,g2,planet)-shift-scale;   
-  shift = -2; + 
-  res = scale*areaOfIntersection(x0, y0, r0, x1, y1, r1)+shift;+
   return res;   return res;
  
-       }, xmin, xmax ] ,{strokeWidth:3});+}, xmin, xmax ] ,{strokeWidth:1.2, strokeColor:"#ff0000",visible: true}); 
 +////////////////////////////////////////////////////// 
 +////////////////////////////////////////////////////// 
 +var dipI = brd.create('functiongraph',   [function(x){
  
 +  scale = 250;
 +  shift = 2;
 +  g1= 0.24515439;
 +  g2= 0.29645038; 
 +  res = scale*quadraticLD(0.35*x,impact,g1,g2,planet)-shift-scale;  
 +  
 +  return res;
  
 +}, xmin, xmax ] ,{strokeWidth:1.2,strokeColor:"#61210B", visible: true});
 +//////////////////////////////////////////////////////
  
 +onlInit();  // second Animation context
  
  
- 
- 
- 
-// 
-// x parameter 
-// p Start/Endpunkte 
-// 
-function line2Points(x,p1X,p1Y,p2X,p2Y){ 
-    m = (p2Y-p1Y) / (p2X-p1X); 
-    y = m* (x-p1X) + p1Y;   
-    return y; 
-} 
- 
-    
-////////////////////////////////////////////////////// 
 } // END OF BLOCK } // END OF BLOCK
  
Zeile 228: Zeile 321:
       brd.fullUpdate();                brd.fullUpdate();         
 }   }  
 +function onShowImpactChanged(value) {
 +  if (value.checked ){
 +              bImpact = true;             
 +            }else { bImpact = false ;} 
 +            brd.fullUpdate();   
 +           
 +                   
 +}
  
  
-function onShowPlanetChanged() {           +function onShowPlanetChanged(value) {     
-            brd.fullUpdate();             + if (value.checked ){ 
-} +              bPlanet = true;              
-function onShowImpactChanged() { +            }else { bPlanet = false ;}  
-            +            brd.fullUpdate();    
-            brd.fullUpdate();            +       
 +                   
 } }
  
 function onPlanetChanged() { function onPlanetChanged() {
             planet = parseFloat($("#planetSlider").val());                           planet = parseFloat($("#planetSlider").val());              
-            $("#planetOut").text(planet.toFixed(2));+            $("#planetOut").text(planet.toFixed(2));      
             brd.fullUpdate();                         brd.fullUpdate();            
 } }
Zeile 257: Zeile 359:
  
 <html> <html>
 +
 <body onload="onInit()"> <body onload="onInit()">
 <center> <center>
 <section>   <div id="jbox" class="jxgbox" style="width:500px; height:500px;"></div> <section>   <div id="jbox" class="jxgbox" style="width:500px; height:500px;"></div>
-</center>  Idealisierte Lichtkurve zur Demostration der Transitphasen. Bei der primären Phase läuft der Exoplanet vor dem Zentralgestirn und dem Beobachter vorbei, bei der sekundären Phase wird der Exoplanet vom Zentralgestirn verdeckt. Obwohl beide Phasen einen Lichtabfall erzeugen, stellt bereits die primäre Phase die Grenze des derzeit machbaren in der Amateurastronomie dar. +</center> Lichtkurve mit Randverdunkelung nach dem Model von Mandel und Algol am Beispiel des Exoplaneten <i>HD 209458 b</i>. Die schwarze Kurve ist ohne Randverdunkelung, die farbigen Kurven entsprechen Kurven mit den Filtern U,B,V,R und I. Die grundsätzliche  Veränderung der Lichtkurve kann in Abhängigkeit der Planetengröße (Rp/R∗) und des Impaktparameters <b>b</b> studiert werden.  Für <i>HD 209458 b</i> sollte <b>b</b> kleine Werte annehmen  und Rp/R∗ 0.12 in der Nähe von 0.12 sein. 
- +<br><br><b>Einstellungen</b><br>   Planet (rel. Größe):  
-<br>Einstellungen<br><br>   Planet (Größe):  +<span id="planetOut" style="display:none"></span>&nbsp;&nbsp;&nbsp;&nbsp; 
-<span id="planetOut" style="display:true"></span>&nbsp;&nbsp;&nbsp;&nbsp; +<input id="planetSlider" type="range" min="0.02" max="0.14" value="0.12"
-<input id="planetSlider" type="range" min="0.01" max="0.19" value="0.1"+
                step="0.01" style="width:250px" oninput="onPlanetChanged()">  <br><br>   Impaktparameter b:   <span id="bOut" style="display:true"></span>&nbsp;&nbsp;&nbsp;&nbsp;                step="0.01" style="width:250px" oninput="onPlanetChanged()">  <br><br>   Impaktparameter b:   <span id="bOut" style="display:true"></span>&nbsp;&nbsp;&nbsp;&nbsp;
-            <input id="bSlider" type="range" min="0.0" max="1.5" value="0.0"+            <input id="bSlider" type="range" min="0.0" max="1.3" value="0.0"
                 step="0.01" style="width:250px" oninput="onImpactChanged()">  <br>                 step="0.01" style="width:250px" oninput="onImpactChanged()">  <br>
                                            
Zeile 277: Zeile 379:
 </div> </div>
  
-</section>   +</section>   
  
 </body> </body>
Zeile 298: Zeile 401:
 </div> </div>
 <div> <div>
-todo+Der Impaktparamter b, hat einen wesentlichen Einfluss auf die Form der Lichtkurve. Der Parameter ist hier  
 +auf den Sternenradius skaliert, so dass der Planet bei einen Wert von b=1 einen 
 +partiellen Transit vollzieht.  
  </div> <br>  </div> <br>
  
Zeile 308: Zeile 414:
 </div> </div>
 <div> <div>
-todo +Die Größenverhältnisse (Rp/R∗) zwischen Exoplanet und Zentralgestirn (engl. Hoststar), verändern 
- </div> </html>+die Form der Lichtkurve ebenfalls. Die Amateurastronomie muss sich (mit derzeitigen Amateurequipment) auf die großen Planeten konzentrieren, d.h. auf das Sonnensystem bezogen gilt: R_J/R_* &asymp; 0.1, da Jupiter der größte Planet im Sonnensystem ist. 
 +</div> 
 +<div><br>Randverdunkelung</div> 
 +<div>Für die Randverdunkelung wurden die Koeffizienten für den Exoplaneten <i>HD 209458 b</i> 
 +in das Modell übernommen. 
 +</div></html> 
 +<sortable> 
 +|Filter|γ_1|γ_2| 
 +|<html><p style="color:#BF00FF"> <b>U</b></p></html>|0.67462231|0.16275864| 
 +|<html><p style="color:#0080FF"> <b>B</b></p></html>|0.57939128|0.21110543| 
 +|<html><p style="color:#04B404"> <b>V</b></p></html>|0.40433600|0.29138968| 
 +|<html><p style="color:#ff0000"> <b>R</b></p></html>|0.31786107|0.30489961|  
 +|<html><p style="color:#61210B"> <b>I</b></p></html>|0.24515439|0.29645038| 
 +</sortable> 
 + 
 +[[https://astroutils.astronomy.osu.edu/exofast/limbdark.shtml| Koeffizienten von Exofast]] 
 +[[https://astroutils.astronomy.osu.edu/exofast/| Exofast]] 
 ++++ ++++
  
-<html> hallo </html> 
  
demo_limbdarkening.1648545187.txt.gz · Zuletzt geändert: 2022/03/29 09:13 von torsten.roehl