summaryrefslogtreecommitdiff
path: root/CHAPTERS/chap3/chapter3.tex
diff options
context:
space:
mode:
Diffstat (limited to 'CHAPTERS/chap3/chapter3.tex')
-rw-r--r--CHAPTERS/chap3/chapter3.tex143
1 files changed, 143 insertions, 0 deletions
diff --git a/CHAPTERS/chap3/chapter3.tex b/CHAPTERS/chap3/chapter3.tex
new file mode 100644
index 0000000..9b400e8
--- /dev/null
+++ b/CHAPTERS/chap3/chapter3.tex
@@ -0,0 +1,143 @@
+\chapter{Importing Mesh From Third Party Software in OpenFOAM}
+\thispagestyle{empty}
+\label{sec:chap3}
+\newcommand{\LocCHthreefig}{\Origin/CHAPTERS/chap3/figures}
+
+OpenFOAM can be used for creating and meshing geometrical shapes like Box, Pipe. When dealing with complex geometries like a turbine blade, aircraft,
+ship etc, we cannot use the blockMesh utility. In such cases it is always better to create the geometry and mesh in dedicated CAD and Meshing softwares
+and solve those usiing OpenFOAM. As a prerequisite it is expected the user should have knowledge about creating geometry and generating mesh in softwares
+like Gmabit, Gmsh, Salome, ICEM etc. This chapter deals with the steps involved in importing mesh files in OpenFOAM using different mesh conversion tools.
+
+\section{Geometry}
+
+We will use the above problem of Flow over a square cylinder as an example for importing mesh file in OpenFOAM. Here we have a square cylinder
+of length 1m and height 1 m. Inlet velocity is set at 1 $\frac{m}{s}$ for Reynolds number (Re) 100. The size of the domain choosen is 60 m by 40 m.
+The boundary conditions are as shown in the , Fig \ref{square} below.
+
+\begin{figure}[t]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/square.png}
+\caption{Flow over square Cylinder}
+\label{square}
+\end{figure}
+
+\section{Meshing}
+
+We have generated a hexhedral mesh for the above geometry with 40000 cells and saved the mesh file as cylmesh.msh.
+The mesh generated is as shown below, Fig \ref{mesh}
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/cylmesh.jpeg}
+\caption{Mesh}
+\label{mesh}
+
+\end{figure}
+
+\section{Importing the mesh file}
+
+In incompressibel solvers go to icoFoam and create a solver inside it by the name \textbf{cylinder}. Now go inside the cavity case and copy the
+\begin{itemize}
+\item 0
+\item system
+\end{itemize}
+
+\flushleft folder and paste it inside the cylinder folder. Please make a not here that we do not need the \textbf{constant} folder here. After this copy the
+cylmesh.msh mesh file create earlier and paste this inside this folder. Thus the our case file is now ready. Now open the command terminal and type the
+path for the cylinder folder. Now since we have a Fluent (.msh) mesh file we will use the mesh conversion command as shown below followed by the file name \\
+\centering \textbf{fluentMeshToFoam file-name.msh} \newline
+
+\flushleft In the terminal window type the above command with the file name and press enter.
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/conversion.png}
+\caption{convert}
+\label{mesh}
+\end{figure}
+
+\flushleft In case you have a 3D mesh file then you can use the command \\ \vspace{0.5cm}
+\centering {\textbf{fluent3DMeshToFoam file-name.msh}} \newline
+
+\flushleft The Fluent mesh file is converted into OpenFOAM mesh file. Now if we look back into our cylinder folder we can see that the "constant"
+folder is now generated. When we open the constant folder we will see that the transport properties file is missing. Since we had converted the
+fluent mesh file into openfoam the fluid property files were missing. Copy the transport property file from the constant folder of cavity case
+and paste this inside the constant folder of cylinder. The trasnportProperties file contains the value of fluid viscosity, we can either change it
+or keep it default. \newline
+
+\flushleft Make a note here that we do not use the \textbf{blockMesh} command here
+
+\section{Boundary Conditions}
+
+When we import the geometry in OpenFOAM we need to be very careful with the boudnary names used while creating the mesh file. Since OpenFOAM is case
+sensitive in case of any mistake with the boundary names can create an error while running the solver. To view the boundary names in the command terminal
+go to polyMesh folder inside the constant. Inside polyMesh you can see a file by the name \textbf{boundary}. Open this file in any editor of your
+choice, eg, gedit boundary, Fig \ref{boundary}.
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/boundary.png}
+\caption{Boundary file}
+\label{boundary}
+\end{figure}
+
+\flushleft The boundary names will be as shown in the domain shown above, Fig \ref{square}. In case of any error with the boundary names you can
+always refer to this boundary file. Now in your command terminal go to the 0 folder and open the pressure file. Make sure that the boundary names
+match exactly the names in the boundary file, in case of errors make the necessary changes.
+
+\section{Solver settings}
+
+In the terminal window go to the controlDict file inside system and open it in any editor of your choice. Change the endTime from 0.5 to 1.5 seconds.
+Save the file and close it, Fig \ref{cd} and come back to the cylinder folder.
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/controldict.png}
+\caption{controlDict file}
+\label{cd}
+\end{figure}
+
+\flushleft After making the necessary changes we can now run the solver. In the temrinal window type the name of the solver \textbf{icoFoam}
+and press enter. The iterations will be seen running on the terminal window. After the iterations stop we can now start with the visualization.
+
+\section{Post-Processing}
+
+Launch paraview by typing \textbf{paraFoam} in the terminal window and once it opens click on the Apply button to view the geometry, Fig \ref{geom}.
+In the active varialble control menu change from Solid Color to Velocity (U). You can now see the initial conditions for velocity, Fig \ref{vel}.
+To view the animation on the right hand top of paraview click on the play button of VCR menu. You can see the change in velocity in the paraview
+window with the passage of time, Fig \ref {vel-1}.
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/geom-paraview.png}
+\caption{Geometry in Paraview}
+\label{geom}
+\end{figure}
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/vel.png}
+\caption{Initial velocity condition}
+\label{vel}
+\end{figure}
+
+\begin{figure}[h]
+\centering
+\includegraphics[width=\lgfig]{\LocCHthreefig/vel-1.png}
+\caption{Velocity at 1 sec}
+\label{vel-1}
+\end{figure}
+
+\section{Mesh Conversion Commands}
+
+The user can also import mesh files from other meshing softwares as well. Here is a list of commands to import mesh files in OpenFOAM.
+
+\begin{itemize}
+\item ANSYS : ansysToFoam file-name
+\item IDEAS : ideasToFoam file-name
+\item CFX : cfxToFoam file-name
+\item SALOME : ideasUnvToFoam file-name
+\end{itemize}
+
+
+