I Share
i want to share my ideas
Monday, August 25, 2014
Jagamantha Kutumbam naadi by Sirivennela Sitaramasatrytelugu song meaning
Wednesday, December 19, 2012
Magical Night at a Beach
ఆవేదన
Wednesday, January 25, 2012
ఆలోచన
నే వింటూనే .. ఓ ఆలోచన
నే చూస్తున్నది వింటున్నది .. ఓ ఆలోచన
నా స్పర్శయిన .. ఓ ఆలోచన
నా శ్వాసయిన .. ఓ ఆలోచన
నా శ్వాసను స్పర్శించినది ఓ ఆలోచన
నా స్నేహము కోరి
నా ముందు వాలి
నా లోన చేరి
నా లోన దాగెను .. ఈ ఆలోచన.
అనువనువుకు చేరి
అనుకువగా మారి
మురిపించెను నన్నే
మరిపించెను నన్నే .. ఈ ఆలోచన.
నా స్వార్ధముగా నను సంధించి
నా స్వ అర్ధమును బంధించి
నా రథమును చేపట్టి
సారథి అయి సంచరించి
సర్వము తానే అనిపించేలచేసేను
సత్యము తానై కనిపించేలా .. ఈ ఆలోచన.
అంధకారములో మునిగేలా
ఆవేశముతో కదిలేలా
ఆనందము తోలగేల
అలజడులు సృష్టించెను .. ఈ ఆలోచన.
కలలాగ చేరినా
కళలాగ మారినా
కలకాలము ఉండునా .. ఈ ఆలోచన ..
ఆలోచనలకు ఆది ఒక ఆలోచనయే
ఆలోచనలకు అంతము ఒక ఆలోచనయే
... రాజేష్ కాట్రెడ్డి
Saturday, January 7, 2012
Direct Numerical Simulation of the Unsteady Navier-Stokes Equations
Problem Statement:
To develop a solver to numerically solve the unsteady Navier-Stokes equation using Direct numerical simulation for a homogeneous isotropic turbulent shear flow.
Theoretical background:
Unsteady Navier-Stokes equation for incompressible flow is given by following equations. The equations shown are non-dimensional form of the NS equations.
There are several methods available to solve the equations numerically. Finite difference method is adopted to solve the simultaneous system of equations 2d simulations can be performed using stream function method but obviously not valid for a 3d simulation. Hence direct numerical simulation of the momentum equations is adopted and hence a poisson equation has to solved simultaneously for pressure distribution to have divergence free flow.
Several methods are established so far in literature to solve these equations like MAC, SIMPLE, SIMPLEC and so on. But the most well know method is the Marker And Cell(MAC) approach.
Staggered Grid:
The staggered grid configuration is a shown in the figure. The pressure values are stored in the center of the cell and the velocity values are stored along the surface of the cell. This way of storing the values helps in eliminating the formation of checkered pressure distribution of the flow field.
Marker and Cell method:
The finite difference form of equation is solved using Marker and Cell method in the present assignment. This formulation gives explicit equations for the velocity in time from the momentum equations. The pressure field is also solved along with these equations as a poisson equation using continuity. The formulation and the pressure estimation are extensively done by T.Sundarajan in Computational fluid flow text book. The derivations have minor correction in notations.
Initial and Boundary Conditions:
2D Lid driven cavity:
The top wall is the moving boundary and has the boundary
Conditions of U= Uwall and V = 0
And all other walls have U = 0, and V= 0 as the boundary conditions,
3D Lid driven cavity:
The top wall is the moving boundary and has the boundary
Conditions of U= Uwall and V = 0
And all other walls have U = 0, and V= 0 as the boundary conditions
Including the front and the back side wall.
3D Lid driven cavity with top and bottom walls Moving:
The top wall is the moving boundary and has the boundary
Conditions of U= Uwall and V = 0
The bottom wall is the moving boundary and has the boundary
Conditions of U = - Uwall and V = 0
And all other walls have U = 0, and V= 0 as the boundary conditions
DNS simulations: 3D
3D Lid driven cavity with top and bottom walls Moving:
The top wall is the moving boundary and has the boundary, Conditions of U= Uwall and V = 0
The bottom wall is the moving boundary and has the boundary,Conditions of U = - Uwall and V = 0
And all other walls have symmetry B.C U wall = U inside, and V wall= V inside.
Results and Discussion:
2D Lid Driven Cavity:
The Results obtained for the cases of Re 100 and Re 400 are shown below. The velocity fields evaluated agrees with the results obtained by Ghia et al. The Velocity contour plot shown in the figure below shows the recirculation zone and the overall rotation of the fluid along with the boundary layer. The velocities are compared on the lines passing through the mid sections.
The 2D code converges for a under relaxation factor of 0.7. The flow is also divergence free.
Re 100:
Re 400:
3D Lid Driven Cavity:
A similar attempt is made to solve flow in a 3d cavity with moving top wall. The velocity contour plots are as shown below are on a plane passing through the mid-point of the cavity. The flow is observed to be similar to the pattern observed in 2d case. An attempt is also made to capture the flow structure when both the top and bottom walls are moving for a 3D cavity. The code was able to capture the flow structures. The comparison data was not available for the 3D case but the contour plots agree with the flow structures obtained in the literature.
Re 100: Top Wall moving from left to right
Re 100: Top Wall moving from right to left , bottom Wall moving from left to right
DNS of Homogeneous and Isotropic flow:
An attempt was made to simulate the DNS results of the homogeneous isotropic flow for a Reynolds number 1000. The velocity and pressure contours at the mid-plane are shown at a time instant. The results show the formation of various structures initial over time of 0.5sec. The results later observed to diverge and the code is unstable with the given input parameters of dx, dy, dz and dt. Hence an attempt has to made to reduce the grid spacing and time step in order to obtain converged values.
\
Velocity contours of a plane close to the symmetric boundary show the development of the roll ups.
REFRERENCES:
1)V.Ghia, K.N.Ghia and C.T.Shin, "High-Re solutions for incompressible flow using the Navier-Stokes equations and a multi-grid method.”
2)T. Sundararajan,K. Muralidhar,”Computational fluid flow and Heat Transfer”, IIT Kanpur series texts.