Week 7: Shock tube simulation project

Objective: To deploy the mesh and create 1D elements on the following component with given cross-section and DOF Procedure: Import the given IGS file into Hypermesh from the model import option.                                   …

Aim: Simulation of a Full-Hydrodynamic Port-Fuel Injection engine in Converge CFD. Introduction and Theory: Fuel injection is the introduction of fuel in an IC engine, most commonly automotive engines, through an injector. Fuel injection was replaced by carburetors from the 1990s onward. The main difference between carburetion…

Aim: Reynolds-averaged Navier-Stokes derivation and analysis Introduction and Theory: The RANS equations are time-averaged equations of motion for fluid flow. The idea behind the equations is Reynolds decomposition, wherein an instantaneous quantity is decomposed into its time-averaged and fluctuating quantities. The RANS…

Aim: To set up transient shock tube simulation using Converge CFD. Aim: To perform CHT simulation on an aluminum pipe with air flowing through it. Introduction and Theory:  A simple shock tube is a tube, rectangular or circular in cross-section, usually constructed of metal, in which gas at low pressure…

Aim: To perform CHT simulation on an aluminum pipe with air flowing through it. Introduction and Theory: CHT(Conjugate Heat Transfer) corresponds to the heat transfer in solids and heat transfer in fluids. In solids, conduction often dominates and in fluids, convection dominates.CHT analysis allows us to analyze…

Aim: To simulate the Prandtl-Meyer Shock flow in Converge CFD. Theory: 1. What is a Shock wave? A shock wave is a type of propagating disturbance that moves faster than the local speed of sound in the medium. Like an ordinary wave, a shock wave carries energy and can propagate through a medium but is characterized…

Aim: Transient simulation for flow through an elbow joint containing a throttle valve. Introduction and Theory: A throttle is a mechanism by which fluid flow is managed by constriction or obstruction. In a traditional IC engine, the throttle valve is part of the air-intake system, which controls the amount of air…

Aim: Steady-state simulation for flow through an elbow joint containing a throttle valve. Introduction and Theory: A throttle is a mechanism by which fluid flow is managed by constriction or obstruction. In a traditional IC engine, the throttle valve is part of the air-intake system, which controls the amount of air…

Aim: To simulate flow over a backward-facing step for different mesh sizes. Introduction and Theory: Backward-Facing Step (BFS) flow is a representative model for separation flows, which can be widely seen in aerodynamic flows (airfoil, spoiler, high attack angle process), engine flows, condensers, vehicles (cars,…

Aim - To simulate flow through a channel for 3 different mesh sizes using Converge CFD software. Introduction and Theory - Flows completely bounded by solid surfaces are called internal flows which include flows through pipes (Round cross-section), ducts (Not Round cross-section), nozzles, diffusers. The basic principles…

Aim: To simulate the Gearbox sloshing effect in Ansys Fluent by varying the lubricating fluid and immersion levels. Introduction and Theory (Part 2): Fluid Sloshing Effect: In fluid dynamics, slosh refers to the movement of liquid inside another object (which is, typically, also undergoing motion). Strictly speaking,…

Aim - To perform and study combustion simulation of a non-premixed model. Introduction and Theory - Combustion or burning is a high-temperature exothermic redox chemical reaction between a fuel(the reductant) and an oxidant, usually, atmospheric oxygen. Complete combustion is stoichiometric concerning the fuel, where there…

Aim: To conduct a parametric study on gate valve simulation by varying the opening of the gate disc. Introduction and Theory: Gate valves work by inserting a rectangular gate or wedge into the path of a moving fluid. They are operated by a threaded stem which connects the actuator (generally a hand wheel or a motor) to…

Aim: To perform steady-state analysis on Cyclone Separator using DPM(Discrete Phase model). Introduction and Theory:  Cyclonic separation is a method of removing particulates from an air, gas, or liquid stream, without the use of filters, through vortex separation. When removing particulate matter from liquid,…

Aim: To perform steady-state conjugate heat transfer analysis on graphics card by varying air velocity. Introduction and Theory: CHT(Conjugate Heat Transfer) corresponds to the heat transfer in solids and heat transfer in fluids. In solids, conduction often dominates and in fluids, convection dominates.CHT analysis…

Aim - To conduct the Rayleigh Taylor CFD simulation. Introduction  - Rayleigh Taylor Instability (RTI) is a phenomenon that occurs at the interface of fluids of different density wherein the lighter fluid pushes the heavier fluid due to the density difference. It is a dynamic process wherein the system tries to reduce…

Aim: To simulate flow through an exhaust port/manifold and analyze the heat transfer characteristics. Introduction and Theory: CHT(Conjugate Heat Transfer) corresponds to the heat transfer in solids and heat transfer in fluids. In solids, conduction often dominates and in fluids, convection dominates.CHT analysis allows…

Aim: To study the aerodynamic effects of flow over the Ahmed body and conduct a suitable grid independence test. Introduction: As the burning of fossil fuels becomes a more pressing issue, manufacturers are introducing more fuel-efficient cars to the market. One main contributor to fuel burn is the car’s aerodynamic…

Aim - Simulation of flow over the cylinder and analyze the von Karman vortex street at various Reynolds numbers. Introduction and Theory - In fluid dynamics, vortex shedding is the oscillating flow that takes place when a fluid such as water or air moves past a bluff body. It takes place at certain velocities and depends…

Given Second Order ODE:   Code 1: %Matlab program to animate a second order ODE close all; clear all; clc; %Given conditions related to ODE L=1; m=1; b=0.05; g=9.81; %Time period is from 0-20 seconds;creating a span of 500 entities for finer %animation t_span=linspace(0,20,500); %1st term indicates angular displacement;2nd…

Aim:To check the mixing effectiveness of a mixing tee geometry with different lengths and momentum ratios. momentum ratio=velocity at cold inlet(y axis in geometry)/velocity at hot inlet(x axis in geometry) The geometry was prepared on spaceclaim by extracting the fluid volume. We get the following 2 geometries: Short…

Wedge angles used - 4 degrees,3 degrees,2 degrees,1 degrees Given data: Reynolds Number=2100(Laminar flow) Density = 997 kg/m^3 Dynamic Viscosity(nu) = 8.9e-4 Pa-s Diameter of pipe =0.025 m Formulae Used: Entry Length=0.06*Re*D Average Velocity=(Re*nu)/(rho*D) Max Velocity=2*Average Velocity Pressure Difference=(32*nu*L*V(avg))/D^2…

Wedge angles used - 4 degrees,3 degrees,2 degrees,1 degrees Given data: Reynolds Number=2100(Laminar flow) Density = 997 kg/m^3 Dynamic Viscosity(nu) = 8.9e-4 Pa-s Diameter of pipe =0.025 m Formulae Used: Entry Length=0.06*Re*D Average Velocity=(Re*nu)/(rho*D) Max Velocity=2*Average Velocity Pressure Difference=(32*nu*L*V(avg))/D^2…

Flow simulation through a Backward facing step. blockMeshDict File: In this folder we specify the geometry related details such as vertices,blocks,boundaries etc. /*--------------------------------*- C++ -*----------------------------------*\\ ========= | \\\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\\\ / O peration…

Non-Conservative Form: Conservative form Equations are: The Solution terms(U),Flex terms(F) and Souce terms(J) respectively are: Primitive variables and their corresponding relations with Solution,Flex and Source terms Main code: clear all close all clc %Inputs n=31; x=linspace(0,3,n); dx=x(2)-x(1); gamma=1.4; %TimeSteps…

Here ,k1 and k2 are the CFL numbers. Code 1: %Matlab program to solve the 2D Heat Conduction problem clear all; close all; clc; %Inputs L=1; nx=20; ny=nx; x=linspace(0,1,nx); y=linspace(0,1,ny); dx=x(2)-x(1); dy=y(2)-y(1); T=298*ones(nx,ny); %Boundary Conditions T_left=400; T_top=600; T_right=800; T_bottom=900; %Initializing…

Assignment: 1.Central Difference Scheme For calculating the fourth order approximation of second order derivative we will require 5 stencils,namely(i-2,i-1,i,i+1,i+2) 2.Skewed Right sided difference Scheme For calculating the fourth order approximation of second order derivative we will require 6 stencils,namely(i,i+1,i+2,i+3,i+4,i+5)…