Skip to content

About

To model supersonic flow around a bullet, calculate Cd and drag forces using MATLAB, and validate these with Fluent CFD simulations, comparing analytical and numerical approaches for accuracy.

Topics

Resources

Stars

0 stars

Watchers

1 watching

Forks

Latest commit

 

History

6 Commits

Folders and files

NameName
Last commit message
Last commit date
 
 
 
 
 
 
 
 
 
 

Repository files navigation

Supersonic Flow Analysis of a Bullet

2D steady RANS simulation of supersonic flow over a bullet at Mach 1.5 and Mach 2 in ANSYS Fluent, with a mesh-convergence study and comparison against an empirical drag correlation in MATLAB.

ANSYS Fluent MATLAB Flow

Pressure Contour Static pressure at Mach 1.5 (100k mesh): bow shock at the nose, expansion over the shoulder, recompression at the base.


Objectives

  • Resolve the shock and expansion structure around a bullet at Mach 1.5 and 2.
  • Compute the drag force and drag coefficient, and check their sensitivity to the mesh.
  • Compare CFD results against a simple empirical drag estimate.

Simulation Setup

Item Setting
Geometry 2D bullet, 8 mm length (planar 2D because of ANSYS Student limits)
Solver Pressure-based, steady state
Turbulence model Standard k-ε
Fluid Air, ideal gas
Free stream Mach 1.5 and Mach 2.0
Meshes Quad-dominant, 50k / 100k / 300k cells
Outputs Pressure, velocity and Mach contours; drag force and $C_d$ monitors
Geometry 100k mesh
Bullet Geometry Mesh

Results

Mesh sensitivity (CFD)

Mach Mesh $C_d$ Drag force*
1.5 50k 0.45 290 N
1.5 300k 0.42 270 N
2.0 50k 0.49 295 N
2.0 300k 0.46 310 N

* Planar 2D model, so forces are per unit depth (1 m).

CFD vs. empirical estimate

The MATLAB scripts use the free-stream conditions (ISA sea level, $T$ = 288.15 K, $p$ = 101 325 Pa) and the empirical supersonic drag correlation

$$ C_d \approx 0.15 + \frac{0.6}{M^2}, \qquad D = \tfrac{1}{2}\rho V^2 A_{ref} C_d $$

Mach CFD $C_d$ (300k) Empirical $C_d$ CFD drag Empirical drag
1.5 0.42 0.417 270 N 266 N
2.0 0.46 0.30 310 N 340 N
Velocity at Mach 2 (100k)
Velocity Contour

Key findings

  • Mach 1.5: CFD and the empirical estimate agree within about 1 % in $C_d$.
  • Mach 2: The correlation predicts a falling $C_d$ as Mach number increases, while the CFD predicts a rising one. The gap (0.30 vs. 0.46) shows the limits of the simple correlation, and also shows that the Mach 2 CFD needs more work (see below).
  • Refining the mesh from 50k to 300k lowers $C_d$ by about 6–7 % at both Mach numbers.

Limitations and next steps

  • Residuals at Mach 2 did not converge well on any of the meshes. A density-based coupled solver, shock-adapted mesh refinement and the k-ω SST model would make these results more reliable.
  • An axisymmetric model would give 3D-equivalent drag values that can be compared directly with ballistic data.

Repository Contents

├── code/
│   ├── mach_number_1.5_calculation.m   # Empirical Cd and drag, Mach 1.5
│   └── mach_number_2_calculation.m     # Empirical Cd and drag, Mach 2
├── docs/
│   └── bullet_analysis_1.pdf           # Full report (30 pages)
└── figures/                            # Geometry, mesh and contour plots

📄 Full report (PDF): theory, all residual, contour and drag histories for the 50k, 100k and 300k meshes.

References

  • Anderson, J. D. (2016). Fundamentals of Aerodynamics, 6th ed. McGraw-Hill.
  • Liepmann, H. W., & Roshko, A. (1957). Elements of Gas Dynamics. Wiley.
  • Carlucci, D. E., & Jacobson, S. S. (2013). Ballistics: Theory and Design of Guns and Ammunition. CRC Press.

Author

Burak Yörükçü, Mechanical Engineering, Çukurova University · GitHub · burakyorukcu@outlook.com

About

To model supersonic flow around a bullet, calculate Cd and drag forces using MATLAB, and validate these with Fluent CFD simulations, comparing analytical and numerical approaches for accuracy.

Topics

Resources

Stars

0 stars

Watchers

1 watching

Forks

Releases

Packages

Contributors

Languages