Changeset 3632 in Sophya for trunk/Cosmo/SimLSS

Timestamp:

May 25, 2009, 1:57:27 PM (16 years ago)

Author:

cmv

Message:

add lobe N-S pour dipole avec regroupement analogique, cmv 25/05/2009

Location:

trunk/Cosmo/SimLSS

Files:

• Makefile (modified) (3 diffs)
• cmvhshsns.cc (added)
• geneutils.cc (modified) (1 diff)
• geneutils.h (modified) (1 diff)

trunk/Cosmo/SimLSS/Makefile

@@ -21 +21 @@
        $(EXE)cmvtstsch $(EXE)cmvtstblack $(EXE)cmvtvarspec $(EXE)cmvdefsurv \
        $(EXE)cmvtintfun $(EXE)cmvconcherr $(EXE)cmvtinterp $(EXE)cmvtstagn \
-       $(EXE)cmvschdist
+       $(EXE)cmvschdist $(EXE)cmvhshsns
  
 PROGSOBJ = \
@@ -28 +28 @@
           $(OBJ)cmvtstsch.o $(OBJ)cmvtstblack.o $(OBJ)cmvtvarspec.o $(OBJ)cmvdefsurv.o \
           $(OBJ)cmvtintfun.o $(OBJ)cmvtinterp.o \
-          $(OBJ)cmvconcherr.o $(OBJ)cmvtluc.o $(OBJ)cmvtstagn.o $(OBJ)cmvschdist.o
+          $(OBJ)cmvconcherr.o $(OBJ)cmvtluc.o $(OBJ)cmvtstagn.o $(OBJ)cmvschdist.o \
+          $(OBJ)cmvhshsns.o
  
 LIBROBJ = \
@@ -272 +273 @@
 
 ##############################################################################
+cmvhshsns: $(EXE)cmvhshsns
+        echo $@ " done"
+$(EXE)cmvhshsns: $(OBJ)cmvhshsns.o $(LIBR)
+        $(CXXLINK) $(CXXREP) -o $@ $(OBJ)cmvhshsns.o $(MYLIB) 
+$(OBJ)cmvhshsns.o: cmvhshsns.cc
+        $(CXXCOMPILE) $(CXXREP) -I$(MYEXTINC) -o $@ cmvhshsns.cc
+
+##############################################################################

trunk/Cosmo/SimLSS/geneutils.cc

@@ -361 +361 @@
 
 //-------------------------------------------------------------------
+/*=========================================================
+Antenne de longueur totale d (lue au milieu cad 2 brins de d/2)
+see Jackson p401  (antenne selon Oz, centre z=0 en O)
+k = w/C = 2 Pi Nu / C = 2 Pi / Lambda
+- densite de courant: J = I sin(kd/2 - k|z|)*delta(x)*delta(y)
+                        (selon Oz) et pour |z|<d/2
+                        I current peak value if kd>=Pi
+- current at gap (center): I0 = I sin(kd/2)
+#
+Z=sqrt(Mu/Epsilon) is the impedance of the vacum (376.7 Ohms)
+  cad pour le vecteur de pointing: |S| = |E0|^2/Z
+- Time average radiated power by solid angle unit:
+ t is theta the polar angle, t=0 along antenna, Pi/2 perpendicular
+ dP/dOmega = Z*I^2/(8Pi^2) * (cos(kd/2*cos(t)) - cos(kd/2))^2 / sin(t)^2
+           if we define  L = d/lambda ,  kd/2= Pi * L
+ dP/dOmega = Z*I^2/(8Pi^2) * (cos(Pi*L*cos(t)) - cos(Pi*L))^2 / sin(t)^2
+#
+- pour t -> 0
+          dP/dOmega = Z*I^2/(8Pi^2) * ( 1/4 * (Pi*L*sin(Pi*L)*t)^2 )
+- pour t -> Pi
+          dP/dOmega = Z*I^2/(8Pi^2) * ( 1/4 * (Pi*L*sin(Pi*L)*(t-Pi))^2 )
+#
+- Polarisation in plane containing antenna and the direction of propagation
+#
+- pour kd << 1  (L->0)  -> approximation du dipole:
+          dP/dOmega = Z*I^2/(8Pi^2) * ( (Pi*L)^4/4 * sin(t)^2 )
+- pour kd=Pi  (L=1/2)  -> half-wave length:
+          dP/dOmega = Z*I^2/(8Pi^2) * cos(Pi/2*cos(t))^2/sin(t)^2
+- pour kd=2*Pi  (L=1)  -> full-wave length:
+          dP/dOmega = Z*I^2/(8Pi^2) * 4*cos(Pi/2*cos(t))^4/sin(t)^2
+#
+- Le programme retourne Val tel que: dP/dOmega = Z*I^2/(8Pi^2) * Val
+    L<0 (arg[0]<0) alors Val est normalise tel que le maximum = 1
+#
+La carte antenne fournie par Peterson:
+         d=8cm, distance entre 2 centres consecutifs 10cm
+#
+# ---- AntCentFed et AntDipole
+# Input:
+#   t = angle entre la direction de la radiation et le fil de l'antenne (radian)
+#   L = d/Lambda  (d=longueur totale de l'antenne)
+# Output:
+#   (8Pi^2) dP/dOmega / (Z*I^2)
+========================================================= */
+
+double AntCentFed(double L, double trad)
+{
+ double pil = M_PI*L;
+ double st = sin(trad), ct = cos(trad);
+ // L'antenne standing-wave
+ double v = 0.;
+ if(fabs(st)<1e-3) {
+   if(ct>0.) {  // theta ~= 0
+     v = pil*sin(pil)*trad/2.;  v *= v;
+   } else {  // theta ~= 180
+     v = -pil*sin(pil)*(trad-M_PI)/2.;  v *= v;
+   }
+ } else {
+   v = (cos(pil*ct) - cos(pil))/st;  v *= v;
+ }
+
+ return v;
+}
+
+double AntDipole(double L, double trad)
+{
+ double pil = M_PI*L;
+ double st = sin(trad);
+ // L'approximation du dipole
+ double vd = pil*pil*st/2.;  vd *= vd;
+ return vd;
+}
+
+//-------------------------------------------------------------------
 
 static unsigned short IntegrateFunc_GlOrder = 0;

trunk/Cosmo/SimLSS/geneutils.h

@@ -94 +94 @@
 double SinNXsX_Sqr(double x,unsigned long N,bool app=false);
 
+double AntCentFed(double L, double trad);
+double AntDipole(double L, double trad);
+
 //----------------------------------------------------
 double IntegrateFunc(GenericFunc& func,double xmin,double xmax