Algorithm for acceleration of projectile undergoing squared velocity drag?

Click For Summary
SUMMARY

The discussion centers on simulating the motion of a projectile experiencing squared velocity drag using Java. The user presents their algorithm for calculating acceleration, defined as a = g - kv², and attempts to derive the change in acceleration over time. They express uncertainty about the correctness of their approach and reference a Wikipedia article on free fall with air resistance, which provides a direct solution for constant air density. The user seeks validation of their mathematical derivation and algorithmic implementation.

PREREQUISITES
  • Understanding of classical mechanics, specifically projectile motion.
  • Familiarity with Java programming language and its syntax.
  • Knowledge of calculus, particularly differentiation and integration.
  • Concept of drag force and its mathematical representation in physics.
NEXT STEPS
  • Study the derivation of motion equations under drag forces, focusing on squared velocity drag.
  • Implement numerical methods for solving differential equations in Java.
  • Explore the use of Java libraries for physics simulations, such as JBox2D.
  • Review advanced topics in fluid dynamics related to drag forces and their effects on projectile motion.
USEFUL FOR

Students and professionals in physics, software developers working on simulation projects, and anyone interested in computational modeling of motion under drag forces.

Levis2
Messages
43
Reaction score
0
I am trying to write a java program to simulate the motion of a projectile undergoing a drag proportional to the velocity squared, but i am having some issues writing the acceleration part. This is my attempt, not sure if its right though;

a=g-kv^2
da/dt=-k dv^2/dt
since dv^2/dt=2vdv/dt;
da/dt=-2vkdv/dt
da=-2vkdv
and
a=dv/dt
dv=adt inserting;
da=-2vkadt
trying to incorporate some indexes; (bear with me, first prog. algorithm :)
a_(n+1)-a_n=-2vka_n Δt
a_(n+1)=a_n-2vka_n Δt

Is this all wrong ?
 
Physics news on Phys.org
wiki already worked out a direct solution for free fall with air resistance (for constant air density which matches your original equation):

wiki_falling_with_air_resistance.htm

Wiki doesn't show the math. To get there, you start with

a = dv/dt = g - k v2

dv / (g - k v2) = dt
 

Similar threads

  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 7 ·
Replies
7
Views
7K
  • · Replies 5 ·
Replies
5
Views
5K
Replies
3
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 18 ·
Replies
18
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
Replies
4
Views
4K