commit cfbe0542199d3c3f6d0f339842006262f5cf44e6
parent 0d69cec6cdb15dcbd48b65a3d0693608091f9c2b
Author: Ivan Gankevich <i.gankevich@spbu.ru>
Date: Sat, 14 Mar 2020 20:46:56 +0300
Roll/velocity plot.
Diffstat:
5 files changed, 52 insertions(+), 31 deletions(-)
diff --git a/Makefile b/Makefile
@@ -16,6 +16,7 @@ all: build/$(SLIDES).pdf
build/$(NAME).pdf: build/inkscape/law-of-reflection.eps
build/$(NAME).pdf: build/gnuplot/verification.eps
+build/$(NAME).pdf: build/gnuplot/velocity.eps
build/$(NAME).pdf: main.tex
build/$(NAME).pdf:
@echo " LATEX $<"
diff --git a/gnuplot/roll b/gnuplot/roll
@@ -1,6 +1,6 @@
-5, 0.246736
-10, 0.246744
-20, 0.251822
-30, 0.278555
-40, 0.362693
-50, 0.410339
+#5 0.246736
+10 0.246744
+20 0.251822
+30 0.278555
+40 0.362693
+50 0.410339
diff --git a/gnuplot/velocity b/gnuplot/velocity
@@ -1,6 +1,6 @@
-5, 0.078010
-10, 0.088081
-20, 0.121901
-30, 0.225737
-40, 0.446216
-50, 0.587744
+#5 0.078010
+10 0.088081
+20 0.121901
+30 0.225737
+40 0.446216
+50 0.587744
diff --git a/gnuplot/velocity.gnuplot b/gnuplot/velocity.gnuplot
@@ -1,9 +1,23 @@
set terminal svg dynamic size 1080/4*1.9, 1080/4 round enhanced font 'Liberation Serif, 12'
set output 'build/gnuplot/velocity.svg'
-s = 180.0/3.14
-static_roll = -0.247885
+s = 180.0/pi
+#static_roll = 0.247885
+static_roll = 0.246736
-#set border 1+2
-#set multiplot layout 1,2
+unset key
+set grid
+set border 1+2 back
+set xtics nomirror out
+set ytics nomirror out
+set xrange [0:*]
+set yrange [0:*]
+load 'gnuplot/rdbu.pal'
-plot 'gnuplot/velocity' using 1:2 with lines
+set multiplot layout 1,2
+
+set xlabel 'Wind velocity, m/s'
+set ylabel 'Roll angle, °' offset 1.5,0
+plot 'gnuplot/roll' using 1:(($2-static_roll)*s) with lines smooth csplines ls 1 lw 2
+
+set ylabel 'Ship velocity, m/s'
+plot 'gnuplot/velocity' using 1:2 with lines smooth csplines ls 1 lw 2
diff --git a/main.tex b/main.tex
@@ -302,20 +302,6 @@ vector in cartesian coordinates). Then the solution is written as
and reduces to general form of the solution for potential flow around a
cylinder given in~\eqref{eq-solution-cylinder}.
-\subsection{Ship roll angle and velocity}
-
-\subsection{Computational performance analysis}
-
-\section{Discussion}
-
-Solution on the boundary \eqref{eq-solution-on-the-boundary} provides simple
-explanation of areas with the highest and lowest pressure for potential flow
-around a cylinder. At left-most and right-most points on the cylinder boundary
-velocity is nought because incident and reflected particle velocities have
-opposite directions and cancel each other out. At top-most and bottom-most points
-incident and reflected particle velocities have the same direction and
-total velocity is two times larger than the velocity of the flow.
-
In order to be compatible with the surface of any object, solution near the
boundary \eqref{eq-solution-near-the-boundary} uses different term \(s\) than
the solution for potential flow around a cylinder which makes reflected
@@ -333,6 +319,26 @@ between the two is less than 1\%.
solution~\eqref{eq-solution-near-the-boundary}.\label{fig-verification}}
\end{figure}
+\subsection{Ship roll angle and velocity}
+
+\begin{figure}
+ \centering
+ \includegraphics{build/gnuplot/velocity.eps}
+ \caption{\label{fig-roll-velocity}}
+\end{figure}
+
+\subsection{Computational performance analysis}
+
+\section{Discussion}
+
+Solution on the boundary \eqref{eq-solution-on-the-boundary} provides simple
+explanation of areas with the highest and lowest pressure for potential flow
+around a cylinder. At left-most and right-most points on the cylinder boundary
+velocity is nought because incident and reflected particle velocities have
+opposite directions and cancel each other out. At top-most and bottom-most points
+incident and reflected particle velocities have the same direction and
+total velocity is two times larger than the velocity of the flow.
+
\section{Conclusion}
