提交: 86a4b8d573167ecfa08b73a06cded33d5aae654f - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: 86a4b8d573167ecfa08b73a06cded33d5aae654f (tree)
• 时间: 2014-08-07 01:05:18
• 作者: Lorenzo Isella <lorenzo.isella@gmai...>
• Commiter: Lorenzo Isella

Log Message

I also calculate the mobility radius in the free molecular regime.

更改概述

• modified: R-codes/projected-gamma.R

差异

diff -r 6a49efe3f11f -r 86a4b8d57316 R-codes/projected-gamma.R
--- a/R-codes/projected-gamma.R	Tue Aug 05 20:58:21 2014 +0200
+++ b/R-codes/projected-gamma.R	Wed Aug 06 18:05:18 2014 +0200
@@ -19,7 +19,7 @@
                                      size = 12, hjust = 0.5,
                                      vjust = 0.5, face = 'bold'),
                 axis.title.x = element_text(size = 20),
-                axis.title.y = element_text(size = 20, angle=90, vjust=.4),
+                axis.title.y = element_text(size = 20, angle=90, vjust=1),
                 axis.text.x = element_text(size=15, colour="black", vjust=1),
                 axis.text.y = element_text(size=15, colour="black", hjust=1),
                 legend.text = element_text(size = 14, vjust=0.4),
@@ -46,11 +46,23 @@
 return(res)
 }
 
+
+Rm_over_R <- function(k,f_gamma){
+
+res <- sqrt(k-(k-1)*f_gamma)
+
+}
+
+
 ##################################################################
-glength <- 100
+glength <- 101
+
+k_list <- c(2,4,8, 10, 15, 20, 30, 50 ,100)
+
 
 N <- 100000
 
+############################
 
 gamma_list <- seq(0,2,length=glength)
 
@@ -79,6 +91,11 @@
 
 names(data) <- c("area", "gamma")
 
+write.table(data,"projected-area-sintering.csv",
+            col.names=TRUE,row.names=FALSE, sep=",")
+
+
+
 gpl <- ggplot(data, aes(x=gamma, y=area
                        )) +
 
@@ -118,8 +135,82 @@
 
 tools::texi2dvi(fn,pdf=T)
 
+############################################
+#Now some calculations of the mobility radius
 
 
+Rm <- c()
+gl <- c()
+kl <- c()
+
+f_gamma <- data$area
+gamma <- data$gamma
+leng <- length(gamma)
+
+for (k in k_list){
+
+
+Rm <-c(Rm, Rm_over_R(k,f_gamma))
+gl <- c(gl, gamma)
+kl <- c(kl, rep(k, leng))
+
+}
+
+data3 <- cbind(Rm,gl,kl)
+
+data3 <- as.data.frame(data3)
+
+names(data3) <- c("Rm", "gamma", "k")
+
+
+
+lbls <- levels(data3$k)
+
+
+
+gpl <- ggplot(data3, aes(x=gamma, y=Rm,
+                        ## colour=variable,
+                        ## ## shape=d
+                        ## ,linetype=variable
+                       )) +
+
+facet_wrap( ~ k, nrow = 3, scales = "free_y"  )+
+
+scale_x_continuous(limits=c(0,2), breaks=seq(0,2)## ,
+                   ## labels=c("$0$", "$\\pi/4$", "$\\pi/2$", "$3/4\\pi$", "$\\pi$
+## ")
+                   ## ,expand = c(0.01,0.01)
+                   )+
+
+
+    
+my_ggplot_theme(c(0.8, 0.1))+
+
+## geom_point(size=4)+
+geom_line(size=1, color="black") + 
+## geom_point(size=3.6) +
+
+    ## theme(legend.text=element_text("ddd"))+
+
+## scale_shape_manual("$d$", breaks=lbls, labels=lbls, values=c(1,2,3,4,5)) +
+## scale_colour_manual("", breaks=lbls, labels=c("$R_g/R_0$ [numerical]", "$R_g/R_0$ [asymptotic]"), values=c(1,2)) +
+## scale_linetype_manual("", breaks=lbls, labels=c("$R_g/R_0$ [numerical]", "$R_g/R_0$ [asymptotic]"), values=c(1,2))+
+
+   ## theme(    axis.text.y = element_text(size=15, colour="black", vjust=1))+
+
+ xlab("$\\gamma$")+
+ ylab("$R_m(\\infty)/R$")
+
+
+fn <-'sintering_Rm_free_molecular.tex' 
+
+tikz(fn, standAlone = TRUE, width=8,height=5)
+
+
+print(gpl)
+dev.off()
+
+tools::texi2dvi(fn,pdf=T)