Recent Changes - Search:

Menu

editer le Menu

Modélisation du transfert radiatif dans le milieu interstellaire



Équipe 3

Attention, contenu explicite et très physique...

NomResponsabilitéCourrielIdentifiant
MaudeCoordination du projetmailto:unnom@hotmail.com roylabma
RaphaëlResponsable du matérielmailto:unnom@hotmail.com brochura
JohannesResponsable des rapportsmailto:unnom@hotmail.com brustljo
GuillaumeResponsable des exposésmailto:unnom@hotmail.com hivertgu
ÉlieResponsable de la qualité du françaismailto:unnom@hotmail.com bouffael
CatherineÀ déterminermailto:unnom@hotmail.com massonct
AlexandreÀ déterminermailto:unnom@hotmail.com carbonal

Équipe de Laval:

ProfesseurLaurent Drissenldrissen@phy.ulaval.ca
ProfesseurGilles JoncasGilles.Joncas@phy.ulaval.ca
Ph.D.Alexandre Alariealexandre.alarie.1@ulaval.ca
Post-Doc coordonnateur du projetDominic Lagroisdominic.lagrois.1@ulaval.ca

Courriel à l'équipe profs/étudiants gradués


Équipe internationale d'experts:

  • Instituto de Astronomía - UNAM Mexico
Expert transfer radiatif/physique atomiqueLuc Binettelbinette@astro.unam.mx

Liens importants:

Table des matières

Session H 2013

Session A 2012


Journal de bord

Le but du journal de bord est de consigner vos notes et questions en lien avec l'évolution de vos apprentissage. Le professeur ou toute autre personne de l'équipe interpelée par votre intervention répondra directement à même le journal de bord.

Catherine?08 avril 2013, 13:26

J'ai changer et ajouter des choses dans le rapport pour que ca fonctionne avec ce quon a fait jusqua maintenant. je vous lai envoyé par mio :)

Raph?08 avril 2013, 12:01

J'ai améliorer le rapport pas mal aujourd'hui, il est toujours dans le même liens j'ai seulement updater le fichier. Vous pouvez relire question de vérifier la cohérence parce que j'ai changer pas mal de truc de place question de ne pas trop mettre l'accent sur le 3D. Il reste quelques trucs à ajouter pour notre partie. Sinon la grande partie du travail est de remanier pour qu'il ait de l'allure sur le wiki haha! Bon travail!

Raph?04 avril 2013, 11:27

J'ai mis ci-dessous le document word de la nouvelle version du rapport. Il y a beaucoup de travail à faire encore pour la version à remettre mardi.

Raph?04 avril 2013, 11:21

Maude et Guillaume?20 mars 2013, 10:58

Notre partie du compte-rendu, temporairement dans le journal de bord comme nous ne savons pas ou le remettre: Lors de la dernière session, nous nous étions chargés de trouver le programme de modélisation le plus efficace pour visualiser nos résultats. Après quelques essais non concluants avec Mappings et Vis5D, notre choix s'est arrêté sur Mayavi2, qui nous apporte une représentation adéquate. Nous avons aussi fait des recherches sur les projets similaires au nôtre pour nous assurer de la pertinence de notre recherche. Cependant, les matrices à modéliser sont beaucoup trop grandes et Mayavi n’est pas capable de les modéliser. Il a donc fallu écrire un programme qui réduit les matrices tout en gardant leur forme pour permettre la visualisation. Ce programme prend donc un bloc 3x3x3 et le réduit à un seul bloc, donc, réduit le fichier d’un facteur 27. Aussi, nous nous sommes attardés à la visualisation de notre projet final: comme nous voudrons comparer le résultat obtenu avec une image en 2D, il fallait trouver un moyen de transformer la matrice 3D en image 2D. Un programme ressortant une vue de face, de droite et de dessus d'une matrice 3D a donc été créé. Dans le futur, il nous reste à programmer l'étape de comparaison qui calculera les différences entre l’image de base et celle modélisée et modifiera les paramètres de l’analyse spatiale pour essayer de réduire ces différences en dessous d’un seuil acceptable. Il faudra aussi porter notre attention sur la création du squelette du programme principal: une fois toutes les parties individuellement complétées, il faudra un programme qui coordonne le tout.

pour maude?21 février 2013, 08:46
        open(unit=1,file='geometry.tmp',status='unknown')
           read(1,*) ncols, nlines
           read(1,*) nfiles
           do n=1,nfiles
              read(1,*) namef(n)
           enddo
	close(unit=1)
        do n=1,nfiles

c print*,namef(n)

          open(unit=2,file=namef(n),status='unknown')
             longueur=(ncols-1)+(ncols+1)*nlines
             read(2,*) (vect(ii),ii=1,longueur)

          close(unit=2)

c vmin=1000000000. c vmax=0.

          do i=1,ncols
             do j=1,nlines
                pos=i+(ncols+1)*j
                out(n,i,j)=vect(pos)
                valeur(i,j)=vect(pos)

c if (vect(pos).lt.vmin) then c vmin=vect(pos) c endif c if (vect(pos).gt.vmax) then c vmax=vect(pos) c endif

             enddo
           enddo

Guillaume?14 février 2013, 11:24

Programme de réduction de matrice, utile pour visualiser (réduire la taille du fichier)

Guillaume?14 février 2013, 11:20
program CreateMatrix

  ! Ce programme prend une matrice VTK et la réduit d'un facteur 27 pour en réduire la taille.

  implicit none

  integer :: i,j,k,m,n,o,temp,numx,numy,numz,x,y,z
  real,allocatable :: B(:,:,:)
  character nomf

16 format(a)
17 format(i6)
18 format(i9)

! Cette section lit un ficher VTK et en sort la matrice.

open(unit=3,file='reduction.in',status="old")

read(3,*) nomf

 close(unit=3)

open(unit=1,file='nomf',status="old")

! Si la matrice débute à la ligne 1 (n'a pas d'en-tête), passer cette section et définir manuellement numx,numy,numz
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*) numx
  read(1,*) numy
  read(1,*) numz
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*)
  read(1,*)


allocate(B(numx,numy,numz))


do z = 1, numz
  do y = 1, numy
    do x = 1, numx
      read(1,*) B(x,y,z)
    end do
  end do
end do


 close(unit=1)


open(unit=2,file='formationmatrice3.vtk',status="replace")


! Le bloc suivant créé l'entête du format .vtk

write(unit=2, fmt=16)'# vtk DataFile Version 2.0'
write(unit=2, fmt=16)'Modélisation de nébuleuse'
write(unit=2, fmt=16)'ASCII'
write(unit=2, fmt=16)''    
write(unit=2, fmt=16)'DATASET STRUCTURED_POINTS'
write(unit=2, fmt=16)'DIMENSIONS'
write(unit=2, fmt=18) int((numx)/3.),int((numy)/3.),int((numz)/3.)
write(unit=2, fmt=16)'ORIGIN    0.000   0.000   0.000'
write(unit=2, fmt=16)'SPACING    1.000   1.000   1.000'
write(unit=2, fmt=16)''        
write(unit=2, fmt=16,advance="no")'POINT_DATA   '
write(unit=2, fmt=18)int((numx)/3.)*int((numy)/3.)*int((numz)/3.)
write(unit=2, fmt=16)'SCALARS scalars float'
write(unit=2, fmt=16)'LOOKUP_TABLE default'

  ! Chaque itération créé une donnée de la matrice
  ! Dans le ficher VTK, la hiérarchie des axes est x --> y --> z  pour ce qui est de l'ordre

  do k = 2, (int((numx+1)/3)-1)*3+1, 3
    do j = 2, (int((numy+1)/3)-1)*3+1, 3
      do i = 2, (int((numz+1)/3)-1)*3+1, 3
        temp=0
        do m = i-1, i+1
          do n = j-1, j+1
            do o = k-1, k+1
              temp=temp+B(m,n,o)
            end do
          end do
        end do

        write(unit=2, fmt=17) temp/27
      end do
    end do
  end do

end program CreateMatrix
 

pour maude?07 février 2013, 10:41

c routine d ecriture des extrants 2d c c c c Copyright (C) 2010 Martin Aube c c This program is free software: you can redistribute it and/or modify c it under the terms of the GNU General Public License as published by c the Free Software Foundation, either version 3 of the License, or c (at your option) any later version. c c This program is distributed in the hope that it will be useful, c but WITHOUT ANY WARRANTY; without even the implied warranty of c MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the c GNU General Public License for more details. c c You should have received a copy of the GNU General Public License c along with this program. If not, see <http://www.gnu.org/licenses/>. c c Contact: martin.aube@cegepsherbrooke.qc.ca c c

       subroutine extrants2d (outfile,valeur,nom,xcell0,ycell0,pixsiz,
     + gain,offset,nbx,nby,valmax)
       integer width
       real valeur(700,1200),xcell0,ycell0,gain,offset
       integer i,j,nbx,nby,valmax
       character*12 nom
       character*40 outfile
       print*,'Ecriture du fichier de ',nom,': ',outfile
       open(unit=1,file=outfile,status='unknown')
       write(1,1000) ! Tag P2).
       write(1,1001) xcell0
       write(1,1002) ycell0
       write(1,1003) pixsiz
       write(1,1004) gain
       write(1,1005) offset
       write(1,*) nbx,nby,valmax
       do i=1,nbx ! Debut de la boucle sur toutes les cases en x.
       do j=1,nby ! Debut de la boucle sur toutes les cases en y.
       valeur(i,j)=(valeur(i,j)-offset)/gain ! Transformation des donnees avec le gain et l'offset et recherche
       enddo ! Fin de la boucle sur toutes les cases en y.
       enddo
       write(1,*) ((nint(valeur(i,j)),i=1,nbx),j=1,nby) ! Ecriture de toutes les donnees qui sont ensuite inscrites

c ! dans la matricel. Ce sont des boucles imbriquees dans la c ! fonction "write" qui couvrent tout le domaine delimite par nbx et nby. c ! L'increment de la boucle sur les ranges (latitute) est de -1 car c ! on considere que dans les fichiers, c ! le nombre en haut a gauche des donnees est a la c ! coordonnee (1,nby).

       close(1) ! Fermeture du fichier relief.pgm.

 1000 format('P2')
 1001 format('# lon0 ',F9.5)
 1002 format('# lat0 ',F9.5)
 1003 format('# pixsiz ',F9.3)
 1004 format('# gain ',F15.11)
 1005 format('# offset ',F15.11)
       return
       end

Martin?05 février 2013, 12:06

Offre d'emploi disponibles pour l'été prochain à l'astrolab du Mont-Mégantic. Je vous recommande fortement de considérer cette possibilité qui serait certainement très agréables et conviviale mais aussi aurait des impacts sur la progression de votre carrière car il s'agit d'un emploi étudiant à caractère scientifique ce qui est plutôt rare.

Si vous voulez en savoir plus ecrivez moi au martin.aube@cegepsherbrooke.qc.ca et je vous enverrai les documents.

Maude?21 décembre 2012, 14:38

Nous avons fini notre partie du rapport et fait l'abstract. Bonnes vacances!!!

Raph?20 décembre 2012, 14:58

Voilà, j'ai modifier quelques petites choses qui étaient dans la liste et fait la conclusion pour notre partie.

Cath?20 décembre 2012, 08:10

C'est bon ca,on pourrait aussi recopier les choses a changer et les effacer au fur et a mesurer quon les a corriger

Maude?19 décembre 2012, 07:03

Bonjour, je propose qu'on modifie chacun les parties du rapport qui nous correspondent. Sinon dans la partie conclusion chaque équipe pourrait écrire ce qu'on a fait et ce qu'on a à faire. Aussi, vu qu'on avait une moins grande partie à faire dans le rapport, guillaume et moi on peut s'occuper de faire l'abstract.

Martin?18 décembre 2012, 08:44

Il faudra penser à convertir nos statistiques spatiales qui viennent de la moyenne des colonnes en z en statistiques spatiales par voxel. Il faut comprendre que les statistiques que nous trouvons pour les lignes et colonnes x et y sont le résultat de la moyenne en z. Il y a donc un atténuation des variations spatiales si ces dernières sont relativement aléatoires. Autrement dis les fluctuation spatiales de voxel a voxel devraient être beaucoup plus grandes que de colonne z intégrée voisines.

Il serait intéressant d'explorer le fameux facteur de clumpiness mentionné par le prof. Joncas pour avoir une idée de cette fluctuation voxel a voxel.

Martin?17 décembre 2012, 12:13

J'ai retrouvé un programme que j'avais fait pour mon directeur de maîtrise il au siècle dernier. Le programme calculait les paramètres physique d'une région HII à partir de son spectre. Peut-être qu'on pourra récupérer certains éléments. Attach:physiqueHII.f

Maude?12 décembre 2012, 21:27

Bonjour, j'ai passé le rapport dans Antidote, même si je crois qu'il est pas nécessairement 100% fini, au moins il y aura ça de moins de fautes

alex?12 décembre 2012, 18:17

J'embarke dessus! Éli t mieux de pas être dessus comme la dernière fois :P

Maude?12 décembre 2012, 16:22

Fini de travailler (pour l'instant), vous pouvez y aller.

Cath?11 décembre 2012, 22:32

Bon, j'ai cleaner le cadre métho. Il faudrait juste que l'équipe Raphael/Élie/Johannes écrive sur la partie statistique du commentaire en rouge, parce que je ne sais pas c'est quoi. voila! pour notre partie, je vais commencer à l'expliquer plus en profondeur à la suite du cadre métho, je pense que ca serait plus clair comme ca (un résumer de tout le projet au début et plus d'explication plus loin)

Cath?11 décembre 2012, 22:05

ou 3 parties

Cath?11 décembre 2012, 22:01

Normalement, vous expliquez votre partie dans le cadre méthodologique plus en profondeur. Est-ce qu'on laisse le cadre comme il est et on rajoute 2 parties qui expliquent chaque équipe ce quon a fait en dessous?

Raph?11 décembre 2012, 21:55

Est-ce qu'on a des parties pré-établies pour travailler sur le rapport ou on y va free-for-all? Parce que ça doit être fait pour jeudi

Elie?07 décembre 2012, 14:39

Elie?07 décembre 2012, 14:38

Elie?07 décembre 2012, 14:37

Elie?07 décembre 2012, 14:37

Elie?07 décembre 2012, 14:36

Maude?30 novembre 2012, 14:42

> triFusionR := proc (a, b) local c1, c2, i; global nombreComparaisonsFusion, ListeRetour, T; if a <> b then c1 := a; c2 := trunc((1/2)*a+(1/2)*b)+1; triFusionR(a, c2-1); triFusionR(c2, b); i := a; while c1 <= trunc((1/2)*a+(1/2)*b) and c2 <= b do if ListeRetour[c1] < ListeRetour[c2] then T[i] := ListeRetour[c1]; c1 := c1+1 else T[i] := ListeRetour[c2]; c2 := c2+1 end if; i := i+1 end do; while c1 <= trunc((1/2)*a+(1/2)*b) do T[i] := ListeRetour[c1]; c1 := c1+1; i := i+1 end do; while c2 <= b do T[i] := ListeRetour[c2]; c2 := c2+1; i := i+1 end do; for i from a to b do ListeRetour[i] := T[i] end do end if end proc;

> triFusion := proc (Liste) local i; global ListeRetour, nombreComparaisonsFusion; nombreComparaisonsFusion := 0; ListeRetour := Liste; triFusionR(1, numelems(Liste)); ListeRetour end proc;

Elie?30 novembre 2012, 13:43
c programme pour trouver les coefficients de Fourier d une ligne ou colonne
c
c      subroutine fourier(datai,ndata,datao,an,bn)
c declaration des variables
      integer ndata,i,n
      real*8 datai(3000),datao(3000),an(2000),bn(2000),L,lambda(2000)
      real*8 f(3000),pi, bnres(2000), anres(2000), a0, maxb, maxa
      integer nha,nhb,Ta(2000),Tb(2000)
c      character
c initialisation des variables  

      open(unit=1,file='inputline.txt',status='unknown')
         do i=1,153
            read(1,*) datai(i)
            f(i)=0.
      enddo
      do i=154,306
        datai(i)=datai(i-153)
      ndata=153*2
      enddo





      do i=1,2000
         an(i)=0.
         bn(i)=0.
      enddo
      do i=1,3000
         datao(i)=0.
      enddo
      a0=0.
      L=real(ndata)/2.
      pi=3.14159
      maxb=0.
      maxa=0.



c main
c      print*,'Nombre d harmoniques?'
c      read*,nh
      nha=1
      nhb=1
      do i=1, ndata
       a0=a0+1./(2.*L)*datai(i)
      enddo
      do n=1,ndata/4
         lambda(n)=2.*L/real(n)
         do i=1,ndata
            an(n)=an(n)+2./(2.*L)*datai(i)*cos(real(n)*real(i-1)
     +      *pi/L)
            bn(n)=bn(n)+2./(2.*L)*datai(i)*sin(real(n)*real(i-1)
     +      *pi/L)
         enddo
      enddo
      do n=1,ndata/4
         if (abs(bn(n)).gt.(maxb)) then
            maxb=abs(bn(n))
         endif
         if (abs(an(n)).gt.(maxa)) then
            maxa=abs(an(n))
         endif
      enddo
      do n=1,ndata/4
         if (abs(bn(n)).lt.(maxb/15.)) then
            Tb(n)=0
         else
            Tb(n)=1
         endif
            bnres(nhb)=Tb(n)*bn(n)
            print*,n, nhb, lambda(n), bnres(nhb)
            nhb=nhb+1

         if (abs(an(n)).lt.(maxa/15.)) then
            Ta(n)=0
         else
            Ta(n)=1
         endif
            anres(nha)=Ta(n)*an(n)
            print*,n, nha, lambda(n), anres(nha)
            nha=nha+1

      enddo
      open (unit=1,file='spectre.txt',status='unknown')
      close(unit=1)

c     maxim=bn(1)
c      do i=1,ndata
c         if(bn(i)).

      print*, 'nhb=',nhb,'nha=',nha,'a0=',a0


c         lambda(n)=2.*L/real(n)
      open (unit=1,file='res.txt',status='unknown')
      open (unit=2,file='inp.txt',status='unknown')
      do i=1,ndata
       f(i)=a0
      do n=1,nhb-1
               f(i)=f(i)+bnres(n)*sin(real(n)
     +      *real(i-1)*pi/L)
         enddo
      do n=1,nha-1
              f(i)=f(i)+anres(n)*cos(real(n)
     +      *real(i-1)*pi/L)
         enddo
         write(1,*) i,f(i)
         write(2,*) i,datai(i)
      enddo    
      close(unit=2)
      close(unit=1)
      stop
      end
 

Martin?26 novembre 2012, 11:36

Guillaume?23 novembre 2012, 15:29
program CreateMatrix

  ! Ce programme en F95 créé une matrice d'entrée pour Mayavi2 (.vtk)

  implicit none

  integer :: i,j,k,numx,numy,numz,numxx,numyy,numzz,w
  real data(50,50,50),valmax1,valmax2,distance1,distance2,xcentre1,ycentre1,zcentre1,xcentre2,ycentre2,zcentre2

  numx=50
  numy=50  ! À titre d'exemple seulement; elles devraient être déclarées dans le programme principal
  numz=50

  valmax1=1000.
  valmax2=500.
  xcentre1=25.
  ycentre1=25.
  zcentre1=25.
  xcentre2=10.
  ycentre2=10.
  zcentre2=10.

  w=1

16 format(a)

18 format(i6)

open(unit=w,file='formationmatrice.vtk',status="replace")

  ! Le bloc suivant créé l'entête du format .vtk

write(unit=w, fmt=16)'# vtk DataFile Version 2.0'
write(unit=w, fmt=16)'Modélisation de nébuleuse'
write(unit=w, fmt=16)'ASCII'
write(unit=w, fmt=16)''    
write(unit=w, fmt=16)'DATASET STRUCTURED_POINTS'
write(unit=w, fmt=16)'DIMENSIONS'
write(unit=w, fmt=18)numx,numy,numz
write(unit=w, fmt=16)'ORIGIN    0.000   0.000   0.000'
write(unit=w, fmt=16)'SPACING    1.000   1.000   1.000'
write(unit=w, fmt=16)''        
write(unit=w, fmt=16,advance="no")'POINT_DATA   '
write(unit=w, fmt=18)(numx)*(numy)*(numz)
write(unit=w, fmt=16)'SCALARS scalars float'
write(unit=w, fmt=16)'LOOKUP_TABLE default'

  do i = 0, (numx-1)
    do j = 0, (numy-1)
      do k = 0, (numz-1)
        distance1=sqrt((xcentre1-real(i))**2.+(ycentre1-real(j))**2.+(zcentre1-real(k))**2.)
        distance2=sqrt((xcentre2-real(i))**2.+(ycentre2-real(j))**2.+(zcentre2-real(k))**2.)
        if (distance1.eq.0.) then
        distance1=1.
        end if
        if (distance2.eq.0.) then
        distance2=1.
        end if
        write(unit=w, fmt=*) valmax1/distance1+valmax2/distance2
      end do
    end do
  end do
end program CreateMatrix
 

Maude?23 novembre 2012, 15:20
real data(50,50,50),valmax1,valmax2,distance1,distance2,xcentre1,ycentre1,zcentre1,xcentre2,ycentre2,zcentre2
integer i,j,k
valmax1=1000
valamx2=500
xcentre1=25
ycentre1=25
zcentre1=25
xcentre2=10
ycentre2=10
zcentre2=10
do i=1,50
do j=1,50
do k=1,50
distance1=sqrt((xcentre1-real(i))**2.+(ycentre1-real(j))**2.+(zcentre1-real(k))**2.)
distance2=sqrt((xcentre2-real(i))**2.+(ycentre2-real(j))**2.+(zcentre2-real(k))**2.)
data(i,j,k)=valmax1/distance1+valmax2/distance2
enddo
enddo
enddo
end

Elie?16 novembre 2012, 16:00

Attach:fourier16nov.f

  • Il fonctionne*

Guillaume?10 novembre 2012, 21:42

Je sais que c'est pas mal en retard comme réponse, mais oui, on l'avait trouvée.

Martin?08 novembre 2012, 22:13

Aviez vous vu cette page pour l'installation de vis5d?

http://vis5d.sourceforge.net/doc/ch02sec2.html

Martin?06 novembre 2012, 22:40

Confirmation, vous allez faire le cours d'intégration en astrophysique la prochaine session sous ma "supervision". Le cours comporte 2h de théorie et 3h de projet. Dans votre cas le projet sera la poursuite du graphycs.

Guillaume?05 novembre 2012, 16:03

Faudrait contacter Martin avant mercredi...

alex?04 novembre 2012, 10:35

Ok en tant que deuxième véhicule! est-ce qu'il va y avoir 2 départs différent finalement? où nous partons tous en même temps?

Catherine?02 novembre 2012, 18:18

Je peux venir mercredi :)

Martin?29 octobre 2012, 18:47

Des nouvelles,

-La demande pour un portable performant qui servira à afficher les données 3D est déposée à la direction. Je ne sais pas combien de temps ça prendra.

-Le souper avec les étudiants/prof de St-Georges et le prof Drissen de U. Laval est confirmé pour le 7 nov.

-J'ai contacté via mio la prof de philo pour lui demander si elle accepte votre absence.

Guillaume?26 octobre 2012, 14:29

Laissez faire ce site-la, fallait juste pouvoir le copier sur un autre ordinateur.

Guillaume?26 octobre 2012, 14:28

Guillaume?26 octobre 2012, 13:41

Voici les info pour le cours de philo quElie Alexandre guillaume et maude manquerait : Letre humain, 340-102-MQ gr 03726 local 3-32-192 prof: Guy Francavilla (cependant, a la base, le cours est donne par Danielle Terzan mais il y a eu une echange de groupe entre les profs)

alex?18 octobre 2012, 15:16

ok je ne suis pu dessus, il me reste encore quelques affaires a faire... je vais probablement complêter se soir

alex?18 octobre 2012, 14:39

Je vais travailler sur le document du compte-rendu... donc touché s'y pas... je vais metre un autre message quand je vais avoir fini!

Catherine?17 octobre 2012, 16:43

Pour le cadre méthodologique, je vais mettre plus d'information sur la partie que je fais (car jen sais plus) que sur vos parties car je ne connais trop ce que vous allez faire... alors vous pouvez rajouter des informations à cet effet.:)

Catherine?13 octobre 2012, 08:19

Séparation du lab : maude guillaume lintro, johannes elie raphael théo, moi alexandre métho

Guillaume?12 octobre 2012, 13:11

Il fonctionne pour moi ici au CEGEP. Essaie d'updater ton Adobe Reader peut-être.

Élie?11 octobre 2012, 21:26

Dans les liens important, le lien «Article de Lagrois, Joncas, Drissen 2012 sur l'objet W4 nébuleuse du cygne» ne fonctionne pas pour moi (il me donne une espèce de longue série de codes).

Martin?10 octobre 2012, 22:57

Les images de la nébuleuse dans chacune de ces raies.

Catherine?10 octobre 2012, 15:44

Quelles données?

Martin?10 octobre 2012, 08:50

Voici la liste des données qu'on a sur la nébuleuse W4 NII 6548A, NII 6584A, SII 6716A, SII 6731A, Ha 6563A, Le bruit

Catherine?09 octobre 2012, 18:52

J'aimerais mieux faire dans la partie méthodologique vu que je n'avais rien chercher sur la théorie au début...

Maude?09 octobre 2012, 12:28

Bonjour, nous (moi, Catherine et Guillaume) avons écrit des sous-points dans la partie compte-rendu pour faire un "plan" si vous pensez a d'autres points rajoutez-les! On pourra par la suite se diviser les points pour avoir des tâches équitables :)

Catherine?07 octobre 2012, 13:39

Qui pourrait mardi midi?

Guillaume?06 octobre 2012, 22:54

Pour le rapport, il faudrait trouver un moment où on peut se rencontrer et se séparer le travail. Peut-être pendant la période de SN mardi (c'est la journée du lundi)?

Sinon, pour la compétition, c'est vrai qu'il faudrait chercher, mais je saurais vraiment pas par où commencer. Il faudrait qu'on trouve le projet du chercheur allemand que Martin a mentionné.

Catherine?05 octobre 2012, 21:24

Alors, qu'est-qu'on fait pour le rapport pour dans 2 semaines?

Et j'y pense,on n'a pas dit qui allait lire sur «la compétition» pour savoir si notre projet vallait la peine.

Catherine?05 octobre 2012, 06:41

Salut!, moi je suis encore au début de mes recherches :P mais jai trouvé ça qui est un bon résumer je trouve des sortes de nébuleuses :Un nuage interstellaire de gaz et de poussières. Les nébuleuses brillantes rayonnent la lumière du gaz qui les composent (nébuleuse à émission) ou reflétent la lumière des étoiles (nébuleuse à réflection), voir les deux. Les nébuleuses sombres sont des nuages de gaz et de poussières qui ne sont pas illuminés (et cache la lumière provenant des autres étoiles). Les Nebuleuses planétaire sont des coquilles de gaz éjectées par les étoiles en fin de vie.

Martin?02 octobre 2012, 08:55

Bonjour, ma journée à U. Laval fut très profitable pour notre projet. J'ai maintenant la piste de recherche très claire (du moins je crois). Selon notre discussion et surtout avec Luc Binette de l'U de Mexico et créateur du modèle MAPPINGS, cet outil serait trop sophistiqué dans un premier temps pour ce qu'on veut faire. Toujours selon l'équipe de Laval une bonne piste pour tout de suite serait d'utiliser les formules de Osterbrock pour des rapports de raies spectrales vs densité, métallicité, température etc. Je dois avoir ce volume dans mon bureau. Sinon je le commanderai. Selon les profs de Laval ce volume est souvent nommé le nouveau testament... Nous commenceront nos analyses avec un objet qui a été observé par le Post-Doc du projet. Il s'agit de W4 un section de la nébuleuse du voile. Il faudra garder ces images confidentielles (très très important) car elle n'ont pas encore fait l'objet de publication.

Je pense que mon idée d'assimilation de donnée recursive les a pas mal intéressé. C'est une méthode très commune en météo mais elle ne semble pas du tout utilisée en astro. Ca pourrait donc être un apport significatif de notre part.

Enfin je vais prendre contact avec un chercheur allemand de Mexico qui aurait fait quelque chose de semblable à ce qu'on veut faire. Il faudra évaluer les différences et la pertinence de faire notre projet de notre côté ou de s'associer à lui pour améliorer le sien.

Guillaume?01 octobre 2012, 13:51

D'abord, bienvenue aux deux nouveaux membres de l'équipe!

Ensuite, je crois que le résumé de Catherine est assez adéquat. Reste donc à analyser les 4 (5?) éléments de la liste et étudier leur comportement typique dans les nébuleuses.

Pour entrer un peu plus dans les détails, il faudra traduire le comportement typique des différents facteurs (température, pression) en fonction de la distance par rapport au centre de la nébuleuse en équation qu'on pourra entrer dans un "générateur de nébuleuse". On pourrait ensuite ajouter un élément de hasard à ces fonctions pour éviter d'obtenir une nébuleuse trop...parfaite.

Catherine?29 septembre 2012, 08:04

Rencontre de vendredi : commencer par une nébuleuse à un temps donné (pas de dynamique), créer une nébuleuse à l'aide d'une matrice(plien de cube) et ¸grace à des fonctions qui correspondent à un certain élément (ex: la température) remplir les cubes. Placer les axes de façon à avoir l'angle le plus simple. Faire une image 2D en additionnant le contenu de chaque cube. Comparer avec une vraie nébuleuse pour améliorer la notre. Commencer avec la nébuleuse la plus simple, mais avec des données réalistes. Une fois qu'on va être capable de faire ca, nous pouvons le mettre en model 3D grâce à des logiciels et avec le temps nous pourrons créer une nébuleuse à l'aide des observations (le contraire). Etpuis nous pourrons ajouter la dynamique (ou avant). La liste qu'on voulait faire sur les éléments à mettre dans notre model.

1. Température 2. Densité ou pression 3. poussière 4. supernova 5. diffusion (?) demander aux autres ce qu'ils observent J'ai bien compris?

Raphaël?28 septembre 2012, 15:11

http://alobel.freeshell.org/wind3d.html Codes existants sur le transfert radiatif

Maude ?28 septembre 2012, 15:04

Une idée pour perfectionner notre modèle de base (à un temps précis) une fois qu'il sera créé: comparer notre modèle pour une nébuleuse précise avec les données réelles de la nébuleuse et "soustraire" les deux "images" pour faire ressortir les différences et ainsi savoir ce qu'il nous reste à corriger/améliorer/ajouter. il faudrait ensuite répeter le processus jusqu'à obtenir une représentation parfaite.

Élie?21 décembre 2012, 14:28

Élie?26 septembre 2012, 18:40

J'ai trouvé des liens intéressant parlant du 3D: http://www.sdsc.edu/~nadeau/PhD/VisualizingStarsAndEmissionNebulas.pdf https://www.mpi-inf.mpg.de/~lintu/papers/vis05.pdf Le deuxième contient même l'explication d'un algorithme!

J'ai aussi trouvé un texte sur les pulsars au cas où: http://iopscience.iop.org/0004-637X/715/2/1248/fulltext/

Martin?25 septembre 2012, 08:25

Les comptes de Maude et Guillaume sont crées sur mammouth

Maude?24 septembre 2012, 16:12

Bonjour, il faudrait que vous ajoutiez votre code d'étudiants (les 6 premières lettres de votre nom de famille et les 2 premières de votre prénom) dans le tableau en haut, Martin en aura besoin pour nous assigner des noms d'utilisateurs pour galileo.

Maude?21 septembre 2012, 16:10

Un manuel sur les "multibody systems" et les différentes méthodes (Euler, Runge-Kutta): http://mat21.etsii.upm.es/mbs/bookPDFs

Maude?21 septembre 2012, 14:05

Guillaume?21 septembre 2012, 09:04

J'ai remarqué que plusieurs articles nommés ici sont sur le même site, alors voici l'adresse où on retrouve un grand nombre d'articles sur l'astrophysique, on pourra y jeter un coup d'oeil si nécessaire.

http://arxiv.org/list/astro-ph/recent

Johannes?20 septembre 2012, 21:39

Pour plus d'infos sur des modèles existants sur la dynamique des milieux interstellaires: http://arxiv.org/pdf/astro-ph/0012444v1.pdf

Martin?20 septembre 2012, 14:18

Il y a effectivement un lien entre la dynamique et l'émission de photons dans le MI via les ondes de choc. Principalement autour de supernovae.

Raphaël?20 septembre 2012, 13:17

J'ai aussi trouvé un document qui explique un modèle en trois dimension de la nébuleuse d'Orion. Je n'ai pas encore pris le temps de le lire, mais il me semble intéressant, parce qu'il explique comment la modélisation a été faite. http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1995ApJ...438..784W&data_type=PDF_HIGH&whole_paper=YES&type=PRINTER&filetype=.pdf

Raphaël?19 septembre 2012, 13:25

en faisant ma recherche sur les vents stellaires, j'ai trouvé un document qui pourrait être intéressant pour pas mal tout le monde pour ce qui est de la compréhension. Voici le lien: http://www.bartol.udel.edu/~owocki/preprints/encyc_hsw.pdf

Élie?13 septembre 2012, 09:44

Je me propose comme responsable de la qualité du français.

Guillaume?12 septembre 2012, 23:07

Comme le dit si bien Bob;

"Jsuis par sur de comprendre... non attend un peu. Jsuis sur de pas comprendre!

Martin?10 septembre 2012, 14:13

Un message de Laurent Drissen de l'U Laval

Salut Martin,

SpIOMM travaille entre 350 et 850 nm. Tu trouveras quelques articles décrivant l'instrument et les résultats scientifiques sur la page web du groupe: http://www.astro.phy.ulaval.ca/publications/2012.htm (même chose pour 2011 et 2010).

Lors d'une réunion de groupe, avec les étudiants, je pourrai faire une présentation sur SpIOMM.

J'imagine qu'une combinaison de données SpIOMM/FantOMM serait parfaite pour le projet, même si on a mis l'emphase sur SpIOMM dans la demande. Gilles et Dominic observent Sh-258 avec le Fabry-Perot cette semaine, et ils ont obtenu un cube SpIOMM il y a un mois.

Laurent

Martin?09 septembre 2012, 20:33

Est-ce que quelqu'un peut faire une recherche sur les modèles existants pour simuler la dynamique du milieu interstellaire?

Martin?09 septembre 2012, 20:33

Salut, j'ai commencé à penser au concept qui pourrait être à la base de notre modèle. Pour la dynamique, il faudra surement utiliser un schéma de type Runge-Kutta ou Verlet. Il faut penser à modéliser le centre de masse d'une cellule du modèle plus une distribution de vitesse autour de la vitesse du centre de masse. On pourrait séparer cette distribution de vitesses en bin et suivre chaque bin. Après un pas de calcul temporel alors on projette la répartition des particules et des vitesses sur la grille et on reprends le calcul. J'ai aussi l'impression qu'on devra modéliser deux gaz en simultané et de manière couplée. Le gaz de particules et le gaz de photons. Peut-être que ce commentaire vous semblera rédigé dans une langue étrangère mais à court terme au moins cela me servira d'aide mémoire.

Guillaume?07 septembre 2012, 20:57

Regardez la section 5.6.8, ca aide un peu à comprendre les variables r et z.

Guillaume?07 septembre 2012, 20:42

Trouvé! Je suis allé voir la section 5.6.8 et...

"For radiative transfer in circumstellar envelopes and disks spherical coordinates are particularly well suited. But for defining the structure of the circumstellar disk it will be more convenient to use cylindrical coordinates (rc, z), where z = 0 is the equatorial plane. The two coordinate systems are related via:

rc = rsin(theta) , z = rcos(theta)"

p.88 du chapitre 5

Pas trop compliqué en fait.

Guillaume?07 septembre 2012, 20:22

Il semblerait que "z" soit bien une hauteur, puisqu'on parle plus loin de la "height z". Quelle hauteur, jsais pas trop...

Je pousserais pas trop loin la recherche puisqu'on ne sait pas vraiment comment on va approcher le projet. Pour le moment, j'essaie simplement de m'introduire au sujet et de connaître les grandes lignes.

Johannes?06 septembre 2012, 19:24

Peut-être un article qui pourrait éventuellement nous aider à trouver la densité de la surface du gaz: http://arxiv.org/pdf/1201.3826v1.pdf

Johannes?06 septembre 2012, 19:16

Salut, j'essaie depuis quelque temps de comprendre la partie de l'article qui décrit comment on peut s'y prendre pour estimer la structure des disques de gaz qui forment, plus tard, des systèmes planétaires. Apparemment il faut deux types de données, la densité de la surface du gaz, Σg(r), et la température T de la surface qui coupe le disque puisque T dépend de z, dans la plupart des situations. Ils ont la fonction p(r,z) qui donne la pression du gaz, où r est le radiant, mais ils ne disent pas ce qui est z...par contre r peut aussi être en fonction de z. Donc z c'est...la composante verticale? Tsé, au bouts, à z(min) et z(max) j'imagine que le radiant diminue...: http://www.repro-tableaux.com/a/grissemann/otmargrissemann-kreisel.html ok, c'est pas la meilleure représentation...peut être que se sont des disques très étranges, j'en ai aucune idée. En tout cas, à partir d'une bonne estimation de T(r) on peut rendre le résultat plus exact, en assumant que T(r) = Tm(r), et répéter une suite d'étapes en boucles, jusqu'à ce que le résultat se stabilise. Mais l'ajout de la variable température me semble très compliqué. :s En somme ça me semble être un article intéressant pour nous, puisqu'on y explique comment trouver l'épaisseur du disque, pas comme dans la section 5.6.8... je ne voudrais quand même pas une modélisation finale ressemblant à une crêpe tsé.. Donc avez-vous compris l'article de la même façon? Savez-vous ce qui est z? Et je ne suis pas vraiment au courant de comment on mesure la densité d'un gaz situé si loin de nous.

Raphaël?03 septembre 2012, 21:43

Bon je ne suis pas vraiment habitué à ce genre de page. Quelqu'un peut me dire comment faire un post-it? .. Ce ne sera pas nécessaire j'ai trouvé :).

Martin?21 août 2012, 16:01

Pour les post-it je vous suggère d'indiquer en début de ligne la date et heure de votre commentaire et de l'effacer lorsqu'il est périmé. Ce sera très utile surtout lorsque vous serez en rédaction des rapports afin de minimiser les risques de conflits de version.

Maude Roy-Labbé?20 août 2012, 21:47

J'ai modifié la tâche "faire un test" pour qu'elle ne soit plus en retard mais bien accomplie.

Guillaume Hivert?10 août 2012, 15:18

Pas une mauvaise idée les post-it Maude, mais je pense pas qu'on va vraiment les utiliser assez pour avoir besoin d'un code de couleur.

Guillaume Hivert?09 août 2012, 17:29

Ca fonctionné avec Internet Explorer, mais pas avec Google Chrome. J'ai réussi à le faire marcher en jouant dans les cookies, merci!

Martin?09 août 2012, 15:00

Salut Guillaume, Je le fais ici sans problème. Essaie de changer de navigateur. J'ai utilisé firefox. Aussi véfifie si les cookies sont activés. Tu verras que mon test a ajouté ddd à la fin de la première ligne

Guillaume Hivert?09 août 2012, 09:16

Je ne peux toujours pas éditer les notes dans les tâches. Le site me demande un mot de passe.

Martin?08 août 2012, 21:32

Je vois que tu as bien compris l'utilité du post-it. Il permettra d'éviter certains conflits d'édition.

Martin?08 août 2012, 21:31

C'est normal car les messages sont générés comme du code wiki directement dans la page.

Martin?2 août 2012, 14:50

Le mot de passe pour la rétroaction est "boboche".

Edit - History - Print - Recent Changes - Search
Page last modified on November 02, 2013, at 01:08 am UTC