\begin{figure}[H]
  \begin{center}
    \begin{tikzpicture}
       \begin{loglogaxis}[log ticks with fixed point,
        xmin=0.1,xmax=1000000,ymin=0.1, ymax=100,
        extra y ticks={0.1, 0.2, 0.4, 0.6, 0.8, 1, 2, 4, 6, 8, 10, 20, 40, 60, 100},
        extra y tick style={log identify minor tick positions=false},
        ylabel={Drag coefficient, $C_D$},
        xlabel={Reynolds Number, $Re=\frac{Ud\rho}{\mu}$},
        grid=both, legend cell align={left},
        legend pos=north east,
        width=0.95\textwidth, height=0.7\textwidth]

        \addplot[blue, mark=o, only marks] table[x,y,col sep=comma] {csv/Cd_Schlichting.csv};
        \addlegendentry{H. Schlichting (1982)\cite{Sch17}};

        \addplot[green,domain=0.1:3] {8*pi/(x*(0.5-0.577216+ln(8/x))))};
        \addlegendentry{Oseen's approximation};

        \addplot[red, thick, only marks, mark=square] coordinates {(10, 3.46)(100, 1.46)(10000, 1.09)};
        \addlegendentry{Simulated};

        \addplot[orange, very thick, only marks] coordinates {(10000, 0.925)};
        \addlegendentry{k-$\epsilon$ model};

        \addplot[thick, only marks, mark=x] coordinates {(15.3, 2.2)(17.3, 2.25)(97.3, 1.26)(105, 1.24)};
        \addlegendentry{D. J. Tritton (1959)\cite{Tri59}};

        \addplot[thick, only marks, mark=o] coordinates {(10, 2.846)(100, 1.056)};
        \addlegendentry{S. C. R. Dennis (1970)\cite{Den70}};

        \addplot[thick, only marks, mark=diamond] coordinates {(100, 1.172)};
        \addlegendentry{A. E. Hamielec (1969)\cite{Ham69}};
  \end{loglogaxis}
    \end{tikzpicture}
    \caption{Drag coefficient vs. Reynolds number for circular cylinder}
    \label{fig:vali1}
  \end{center}
 \end{figure}