Air Standard Cycle in MATLAB

AIM: Finite Volume Method What is the Finite Volume Method(FVM),  Major differences between FDM & FVM  Need for interpolation schemes and flux limiters in FVM. INTRODUCTION: The Finite Volume (or control volume) Method or FVM is a numerical technique that transforms the partial differential…

OBJECTIVES: Derive both the conservation and non-conservation forms of the governing equations and solve them using MacCormack's technique. Steady-state distribution of primitive variables inside the nozzle. Time-wise variation of the primitive variables. Variation of Mass flow rate distribution inside the nozzle at different…

1. OBJECTIVE: To investigate the effect of preheating (298K to 600K) on the adiabatic flame temperature. To investigate the effect of preheating on combustion efficiency (LHV = 50MJ/kg, and furnace exit temperature = 1700K). To investigate the effect of preheating on fuel saving.   2. THEORETICAL BACKGROUND:  …

OBJECTIVES: To study the effect of equivalence ratio on the final Adiabatic Flame Temperature (in a Constant Volume chamber) and compare the results with Cantera. To study the Adiabatic Flame Temperature variation with heat loss in a Constant Pressure chamber. To calculate the Adiabatic Flame Temperature for different…

OBJECTIVE: To set-up and perform the simulation in OpemFOAM. To measure the velocity profile at 0.085 m from the inlet of the geometry. To explain why the Grading/Growth factor is essential for meshing. INTRODUCTION: OpenFOAM is a free, open-source CFD software that is a three-dimensional CFD tool, based on the finite…

GIVEN: Coefficient Matrix, A=[5 1 2 ; -3 9 4; 1 2 -7]              Solution Matrix, X=[x1;x2;x3]              Residual Matrix, B=[10 -14 33]   OBJECTIVES: To solve the system of matrices using Jacobi, Gauss-seidel, and SOR methods. To compute…

OBJECTIVES: To write a code in Matlab for : The PV diagram of an Air Standard Otto Cycle To calculate the Thermal Efficiency of the engine.   INTRODUCTION: Otto cycle is a gas power cycle that is used in the spark-ignition internal combustion engine. An Otto cycle consists of two isentropic (reversible adiabatic)…

OBJECTIVES: To write code to fit a linear and cubic polynomial for the Cp data. To plot the linear and cubic fit curves along with the raw data points. To explain the parameters used to measure the fitness characteristics for both the curves.   INTRODUCTION: Curve fitting examines the relationship between one or more…

AIM: To write a program in Matlab to simulate the Forward Kinematics of a 2R Robotic Arm.   INTRODUCTION:                          The Robotic Arm consists of 2 arms i.e. Link 1 (length l1) and Link 2 (length l2). Link 1 is fixed at base coordinates…

AIM: To calculate the Drag Force for Flow over a Bicycle using Matlab programming.   INTRODUCTION: Drag is a mechanical force. It is generated by the interaction and contact of a solid body with a fluid (liquid or gas). It is not generated by a force field, in the sense of a gravitational field or an electromagnetic…

AIM:  To derive fourth-order approximations of the second-order derivatives for Central Difference, Skewed Right-Sided Difference, and Skewed Left-Sided Difference schemes. To prove that these schemes are fourth-order accurate. To write programs in Matlab to determine the second-order derivative of the function…

AIM: To write a program in Matlab to simulate the transient behavior of a simple pendulum and to create an animation of its motion.                                                 THEORY: The Ordinary…

OBJECTIVE: To write a function that accepts the time step as an argument and solve the 1D Convection Equation. Also, to study the effect of time step size on the numerical solution. To study the CFL criteria that govern the linear convection equation and see if this criterion applies to an implicit solution. To explain…

AIM: To solve the 2D heat conduction equation by using the point iterative techniques for steady-state and transient state analysis using Matlab.   GIVEN: Top Boundary = 600 K Bottom Boundary = 900 K Left Boundary = 400 K Right Boundary = 800 K Absolute error criteria =  1e-4   OBJECTIVES: To solve the steady-state…

OBJECTIVES: To study the effect of time step on the stability of the solution for implicit and explicit schemes To discuss the CFL number and its effect on the stability of the solution. Also, to plot the unstable solution. THEORY: Courant-Friedrichs-Lewy (CFL) condition:  In order to converge explicit and semi-implicit…

AIM: To compare the convergence rate (number of iterations) of the steady-state and transient (unsteady)-state heat conduction analysis using Matlab GIVEN: Top Boundary = 600 K; Bottom Boundary = 900 K; Left Boundary = 400 K; Right Boundary = 800 K; Absolute error criteria =  1e-4 THEORY:      (This…

AIM:  To derive fourth-order approximations of the second-order derivatives for Central Difference, Skewed Right-Sided Difference, and Skewed Left-Sided Difference schemes. To prove that these schemes are fourth-order accurate. To write programs in Matlab to determine the second-order derivative of the function…