Academic Projects

FLUID DYNAMICS

CC 299. High-Resolution Numerical Methods for CFD. (Instituto Tecnológico de Aeronáutica - ITA)

Projects.

- Project 4. Own Two-Dimensional Solver - Euler Equations in Curvilinear Coordinate System. [Report, C code included]

Roe scheme, with entropy fix (1st Order)
Explicit Harten TVD scheme (2nd Order)

Results for:

Shock Reflection Test Case Double Mach Reflection

Gun Tunnel - Mach 4 Implosion Symmetry Test Case

Application of this solver for a real life application: Landing Gear Certification - Tire Burst. Effective economy of ~10 destructive burst tests.

- Project 3. Shock Tube Problem - Numerical Fluxes Part III. [Report, C code included]

Explicit Harten TVD scheme (1st Order)
Explicit Harten TVD scheme (2nd Order)
- Implicit Harten TVD scheme (1st Order)
- Implicit Harten TVD scheme (2nd Order)

- Project 2. Shock Tube Problem - Numerical Fluxes Part II. [Report, C code included]

Explicit Steger-Warming Flux Vector Splitting (FVS) Method (1st Order, 2nd Order)
Explicit Van-Leer non-MUSCL scheme (1st Order, 2nd Order)
- Explicit Liou's FVS Method (AUSM+)
Roe's Method (including Entropy Fix)
- Implicit Steger-Warming FVS (1st Order, 2nd Order)

- Project 1. Shock Tube Problem - Numerical Fluxes Part I. [Report, C code included]

Centered Scheme - Second Order Linear dissipation model
Centered Scheme - Fourth Order Linear dissipation model
Centered Scheme - Non-Linear dissipation model (Pulliam)
MacCormack method
Lax-Wendroff method
Implicit Beam-Warming method

CC 297. Fundamentals of Computational Fluid Dynamics. (Instituto Tecnológico de Aeronáutica - ITA)

Projects.

- Project 3. Full-Potential Solver 2-D, structured mesh [Report, C code included (soon)], [Ref. Holst, 1983]

Full Potential flow over airfoil configurations.
α sequence for convergence acceleration of AF2
Mesh Refinement
Elliptic vs Parabolic mesh comparison
- Lifting airfoil
Shock capturing

Pressure Contours Pressure Contours

over biconvex airfoil over cylinder

- Project 2. Parabolic and Elliptic mesh generation. [Report, C code included]

O-grid.
Biconvex Airfoil, NACA 0012, Circular Cylinder
ADI, AF2 relaxation methods
Periodic Boundary Conditions (Periodic Thomas Algorithm)
Control Functions (P, Q) to control point distribution
(To do) Adaptive control, Smooth out P and Q, obtain airfoil given coordinate file.

ξ and η line control Points (ξ , η) control

- Project 1. Small Perturbation Theory. [Report, C code included]

Laplace's Equation.
Iterative Methods (Jacobi, Gauss-Seidel, SOR, Line-Gauss-Seidel, SLOR, ADI, AF2)
Biconvex Airfoil
Perturbation theory

MAE 471. Computational Fluid Dynamics. (Arizona State University - ASU)

Projects.

- Project 2. MacCormack Method - Quasi 1-D Laval Nozzle [MATLAB code] (Really detailed assignment)

Predictor Corrector Method.

- Project 3. Momentum equations, Viscous Burgers Equations 2-D [MATLAB code]

Forward in Time, Centered in Space (FTCS)

- Final Project. Fractional Step Method on a staggered grid - Lid Driven Cavity [MATLAB code, C file] (Translating to C, not finished yet)

- Forward in Time, Centered in Space scheme
- Projection of Lagrange multiplier (''Pressure") onto solenoidal (Divergence free) space.

(Some problems concerning convergence for higher Re, which is expected, but not that much)

Homeworks.

- [HW 3]. Oscillating Plate Problem. [MATLAB Code]

FTCS explicit
Backwards in Time, Centered in Space implict (Laasonen scheme)

[HW 4]. 1-D Convection Problem. [Solution, MATLAB code]

- Upwind Method
- Lax-Wendroff Method
- BTCS
- MacCormack Method
- Von-Neumann stability analysis

[HW 5]. 1-D Burgers Equation.

First Order TVD scheme (FCT scheme). [MATLAB code]
Accuracy demonstration. [MATLAB code]
Second Order TVD scheme (MinMod limiter). [MATLAB code]

APM 524. Spectral Methods. (Arizona State University - ASU)

[HW 3] .

- 1D Burgers Equation - Fast Fourier Transform (FFT) [MATLAB code]

- Linear Advection Diffusion of a vortex blob - RK4 for first 2 time steps, Adams-Bashforth third order time step and FFT for spacial derivatives. [MATLAB code]

[HW 4]. Soon

MAE 362. High Speed Aerodynamics. (Arizona State University - ASU)

- Final Project. Murmann and Cole Difference Scheme - Numerical Solutions for Transonic Flow Report [MATLAB code]

Perturbation Theory
Prandtl-Glauert Equation
Poisson System (Gauss Seidel)
Pretty cool mesh definition.

EN 2217. Viscous Fluid Dynamics. (Universidade Federal do ABC - UFABC)

[Final Project] - Polhausen Method around a 2-D cylinder for boundary layer calculations. [MATLAB Code] (Portuguese Version Only)

2D Flow around cylinder - Potential flow.
Code diverges for \theta = pi/2, which is expected. However, it can't reach separation point and, hence, a more complete method (Thwaites, p.e.) should be applied.

STRUCTURES

EN2212. Vibration Analysis. (Universidade Federal do ABC - UFABC)

4-Stage Building. 4 Degrees of freedom problem. [MATLAB Code]

Modal response due to a harmonic excitation (VIV problem)
Modal response due to a non-periodic excitation

EN 3213. Theory of Elasticity. (Universidade Federal do ABC - UFABC)

Finite Element Method for the linear elasticity of a beam - Energy Method. [MATLAB Code] (Adapted from Professor Reyolando Brasil code, UFABC; the method part - which is the most important - is not my code! )