Menu

All Courses / undefined

All Courses / undefined / undefined

I am going to write a program in MATLAB for solving Quasi one dimensional nozzle flows using MacCormack's Technique. I will be using Initial and boundary conditions which are provided in Jhon D Anderson's textbook. I am posting the solutions of non-conservative form and working on the conservative part which I will be…

Rohan Shah

updated on 03 Jul 2018

Project Details

Loading...

Leave a comment

Thanks for choosing to leave a comment. Please keep in mind that all the comments are moderated as per our comment policy, and your email will not be published for privacy reasons. Please leave a personal & meaningful conversation.

PleaseÂ loginÂ to add a comment

Other comments...

No comments yet!

Be the first to add a comment

Be the first to add a comment

Read more Projects by Rohan Shah (31)

Planck's Radiation Law | Quantum Physics

Objective:

I was going through Max Planck's radiation law for a couple of days and the way that Planck had derived his beautiful equations just by assuming that the energy is not continuous rather it is an integral multiple of h (Planck constant) just amazed me. I couldn't resist my temptation after realizing that these formulas…

25 Aug 2019 12:45 PM IST

Planck's Radiation Law | Quantum Physics

Objective:

I was going through Max Planck's radiation law for a couple of days and the way that Planck had derived his beautiful equations just by assuming that the energy is not continuous rather it is an integral multiple of h (Planck constant) just amazed me. I couldn't resist my temptation after realizing that these formulas…

24 Aug 2019 07:55 AM IST

Determining minimum cushion pressure to break the thickness of ice field | Newton - Raphson Method

Objective:

I am going to write a code in python for determining the minimum cushion pressure required to break the given thickness of ice by implementing the Newton-Raphson technique. I will also be running a setup to find the optimum relaxation factor for the thickness of 0.6 feet. Equations:The equation that was derived by Muller…

15 Oct 2018 07:57 PM IST

Constraint Minimization Problem

Objective:

I am going to minimize a non-linear function by using Lagrange Multipliers. Problem statement:Minimize the function `5-(x-2)^2 -2(y-1)^2` subject to `x + 4y - 3 = 0`. Solution: `L = f - lambdag`, where ` f = 5-(x-2)^2 -2(y-1)^2` and `g = x + 4y - 3` Thus, I solve `L = 5-(x-2)^2…

06 Aug 2018 01:40 AM IST

Numerical Solution of Quasi 1D Nozzle Flows: MacCormack's Technique

Objective:

I am going to write a program in MATLAB for solving Quasi one dimensional nozzle flows using MacCormack's Technique implementing both of the conservative and non-conservative forms of governing equations and perform the grid dependence test. Note, I will use the set of initial conditions provided in Anderson's…

20 Jul 2018 10:22 AM IST

Numerical Solution of Quasi 1D Nozzle Flows: MacCormack's Technique

Objective:

I am going to write a program in MATLAB for solving Quasi one dimensional nozzle flows using MacCormack's Technique. I will be using Initial and boundary conditions which are provided in Jhon D Anderson's textbook. I am posting the solutions of non-conservative form and working on the conservative part which I will be…

03 Jul 2018 06:19 PM IST

A 2R Robotic Arm simulator

Objective:

I am going to write a program in python to simulate a 2R Robotic Arm. I've used online GIF maker to stitch all of my figures and generate the GIF. Note that, I have just used it to stitch the figures, everything else is done by the program. All the figures (200) have been uploaded here. Program: # Simulating a 2R Robotic…

03 Jul 2018 02:51 PM IST

Numerical Solution of Quasi 1D Nozzle Flows: MacCormack's Technique

Objective:

I am going to write a program in MATLAB for solving Quasi one dimensional nozzle flows using MacCormack's Technique implementing both of the conservative and non-conservative forms of governing equations and perform the grid dependence test. Note: I will use the set of initial conditions provided in Anderson's textbook…

03 Jul 2018 02:09 PM IST

Data Visualizer Tool

Objective:

I am going to write a program in python to parse information from the CONVERGE output file (IC engine simulation results) and use that data for basic performance calculations for engine and generating appropriate plots. Aim:1. Plot the data for any columns and save it by giving an appropriate name and labels. …

19 Jun 2018 02:17 AM IST

Demonstrating the Curve Fitting

Objective:

I am going to write a program in python to demonstrate the curve fit on the data set provided by increasing the degree of polynomials namely first degree, third degree and fourth degree. Program: # Program to demonstrate Curve Fit by increasing the degree of the polynomials import matplotlib.pyplot as plt import numpy…

14 Jun 2018 08:56 PM IST

Otto cycle - Calculating the efficiency of the cycle and generating P - V Diagram

Objective:

I am going to write a program in python to simulate the Otto cycle. I will be calculating the thermal efficiency of the engine and be generating the P - V diagram. Problem Setup and geometric parameters:1. Gamma is set to `1.4`.2. The pressure at state point `1` is set to `101325` bars.3. The temperature…

14 Jun 2018 09:07 AM IST

Calculating Drag Force over a Bicycle

Objective:

I am going to write a program in python which will calculate the drag force over a bicycle and plot it against velocity and drag coefficients. Problem Setup:1. Velocities to be used are `2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12` and `13`2. Drag Coefficients to be used are `0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9`…

14 Jun 2018 08:18 AM IST

Simulating the Ideal Rankine Cycle

Objective:

I am going to write a MATLAB program to simulate the Ideal Rankine Cycle. I will be using standard Xsteam library to get the steam table data for water. Problem setup: 1. The Pressure at the Turbine Inlet is set to `44` bars. 2. The Temperature at the Turbine inlet is set to `538` °C.3.…

14 Jun 2018 01:18 AM IST

Units of quantities

Objective:

I am going to write down the expressions and units for the quantities below:1. Kinematic viscosity2. Dynamic viscosity3. Peclet number4. Nusselt number5. Momentum flux6. Fourier's law of heat conduction7. Fick's law of diffusion8. Wall shear stress9. Gradient of Temperature field10. Divergence of velocity field …

25 May 2018 07:32 AM IST

Basic Math

Objective:

I am going to explain the meaning of divergence, gradient, curl and rank of a matrix along with examples. Vector differential operator:The vector differential operator `grad`, pronounced as 'del' or 'nabla', is defined in three dimensions to be`nabla = partial/(partialx)hati + partial/(partialy)hatj + partial/(partialz)hatk`…

25 May 2018 07:32 AM IST

Coding Basics

Objective:

Introduction Programming statements: %1st statement A=[1,2,3,4,5] %2nd statement B=[1;2;3;4;5] %3rd statement C = A*B' %4th statement C = B*A %5th statement A(1:end) %6th statement A(1:2) %7th statement B(1:2,3) %8th statement B(length(B)-2) Here I will be explaining what does each of the above programming statement mean…

25 May 2018 07:32 AM IST

Row and column vector operations in MATLAB

Objective:

How do we perform row and column vector operations: Code for the above operations: %Row vector operation x = [1 2 3 4 5] %Column vector operation y = [6;7;8;9;10]

25 May 2018 07:31 AM IST

Implicit form - transient case

Objective:

Crank Nicolson method implementation `(delT)/(dt) = alpha((del^2T)/(delx^2) + (del^2T)/(dely^2))` `(T_(i,j)^(n+1) - T_(i,j)^n)/(Deltat) = alpha(1/2(T_(i+1,j)^(n+1) - 2T_(i,j)^(n+1) + T_(i-1,j)^(n+1) + T_(i,j+1)^n - 2T_(i,j)^n + T_(i,j-1)^n))/(Deltax^2) + alpha(1/2(T_(i+1,j)^(n+1) - 2T_(i,j)^(n+1) + T_(i-1,j)^(n+1) + T_(i,j+1)^n…

24 May 2018 02:42 AM IST

Calculating drag force over a bicycle.

Objective:

I am going to calculate drag force and plot it vs velocity and drag coefficient using MATLAB Coding for drag force vs Velocity: As `F prop v^2` where `F` is drag force and `v` is velocity, the plot will have a quadratic curve. Code: %Calculating the drag force over the bicycle and plotting drag force vs velocity…

18 May 2018 09:32 AM IST

Spectral radius

Objective:

I am going to code in MATLAB for solving linear system `Ax = b` using iterative methods namely, Jacobi, Gauss-Seidel, and Successive over-relaxation. I will be calculating the eigenvalues and spectral radius of iteration matrix. The aim of this project is to notice the behavior of spectral radius as the diagonal…

13 May 2018 04:33 AM IST

Solving 2nd order ode for simple pendulum and simulating it in MATLAB

Objective:

I am going to solve a second order ODE for the simple pendulum and created a simulation of that pendulum in MATLAB. How does the function program to solve the pendulum ODE look like? Code for calling the function: %Calling a function - pendulum ode solver function [dtheta_dt] = pendulum_ode_solver(t, theta) g =…

05 May 2018 10:43 AM IST

Parsing the NASA's thermodynamic data file

Objective:

I am going to write a MATLAB program to parse the NASA's thermodynamic data file namely 'therm.txt'. I will then calculate the thermodynamic properties of all the 53 species and plot them against the temperature. Function filesThe series of functions used in the driver script is as follows:1. 'extract_data.m': %…

05 May 2018 06:08 AM IST

Solve a 2D steady state heat conduction equation explicitly using point iterative techniques.

Objective:

I am going to write a program in Matlab to solve a two-dimensional steady-state equation using point iterative techniques namely, Jacobi, Gauss-Seidel, and Successive Over-relaxation methods. Problem Setup:1. The domain is a unit square.2. The number of points along the x-direction is equal to the number…

29 Apr 2018 04:11 PM IST

Effect of time step size on the numerical solution.

Objective:

I am going to write a function in MATLAB that accepts time step as an argument and solves the same problem as mentioned here and compare the effect of time step on the numerical solution. Problem Setup - Addition of assumptions to the problem setup discussed above:1. Grid points set to `80`.2. Time…

29 Apr 2018 12:17 AM IST

Numerical solution to the wave equation - Explicit Method.

Objective:

I am going to write a program in MATLAB which will compare initial and final velocity profile for 1D Linear convection for different value of grid points. 1D linear convection is described as follows: …

25 Apr 2018 08:08 PM IST

Deriving 4th order approximation of the 2nd ordered derivative using CDS, FDS and BDS.

Objective:

I am going to derive fourth order accurate approximations of second order derivative using three differencing schemes namely, central differencing, skewed right sided differencing, and skewed left sided differencing. I will then be writing a program in MATLAB which will compare these schemes. DERIVING SCHEMES: …

19 Apr 2018 06:26 AM IST

Stoichiometric Combustion Calculator for Alkane, Alkene and Alkyne fuels.

Objective:

I am going to write a program which will save the plot generated using MATLAB for comparing the effect of a number of moles of carbon on stoichiometric coefficient (ar). I will also describe how does this trend change with the fuel type (alkanes, alkenes, alkynes) by looking at the plot generated. Below…

19 Apr 2018 06:19 AM IST

Otto cycle - Calculating the efficiency of the cycle and generating P-V plot.

Objective:

In this project, I have written a program in MATLAB which calculates the efficiency of the Otto cycle and creates a PV diagram for it. How does the function code look like? Calling a function code which will solve the piston kinematics equation to calculate the volume trace: %Calling a function [V] function…

19 Apr 2018 06:11 AM IST

Solving the Stalagmite mathematical model.

Objective:

In this program, I have used MATLAB to code for getting an optimum solution (maximum value in this case) of stalagmite function using the genetic algorithm. The genetic algorithm is an optimization technique based on the evolution theory which was proposed by Charles Darwin and it is used for solving optimization…

19 Apr 2018 06:10 AM IST

Discritization Basics - Generating bar graph displaying errors.

Objective:

I am going to write a program in MATLAB which will plot a bar graph showing the global truncation error for the 1st, 2nd and 4th order approximation of the first derivative at a step size of 0.1. Analytical function and its derivative:Analytical function: `f(x) = sin(x)/x^3`, Analytical derivative: `f'(x) = (x^3cos(x)…

19 Apr 2018 06:09 AM IST

Discretization basics - Looking at the effect of dx on error.

Objective:

I am going to write a program in MATLAB which will show the trend of how global truncation error for 1st, 2nd and 4th order approximations of the first derivative changes when we change the step size. Analytical function and its analytical derivative: Analytical function: `f(x) = sin(x)/x^3`, Analytical derivative:…

19 Apr 2018 06:08 AM IST

Showing 1 of 31 projects

Try our top engineering courses, projects & workshops today!Book a FREE Demo