Velocity Dependent Acceleration within Drag Force

In summary: You can solve it by using a couple of techniques, but I'm not really an expert in that area. So, you might want to consult a more experienced person for that. In summary, integrating both sides of the equation won't help you solve the problem, and you might want to consult someone more experienced in solving differential equations.
  • #1
1
0
For fun, I wanted to crate a program that dealt with 2d projectile motion w/ wind resistance. Specifically, where the object started from rest at a given height. I looked up several equations that dealt with wind resistance. Specifically, drag force. Drag force was defined as Fd = CdPAV^2/2. So, I thought acceleration would be a = CdPAV^2/2*m. Dealing with just vertical acceleration the equation would read a = g - Fd/m. I thought that I could double integrate both the vertical and horizontal accelerations and find, parametrically, a general position v. position equation, for that object. The problem that was immediately apparent was that Fd was dependent on velocity. Since this is acceleration equation, how do I integrate velocity FOR velocity. I figure I am doing something wrong. Could someone explain what?
 
Physics news on Phys.org
  • #2
Note this post only describes the math for the y-axis. (The x-axis movement would not involve g).

This can be integrated directly to determine velocity versus time, but apparently not position versus time. Wiki article with the answers, but not the derivation:

http://en.wikipedia.org/wiki/Drag_(physics)#Velocity_of_a_falling_object

http://en.wikipedia.org/wiki/Terminal_velocity#Derivation_for_terminal_velocity

derivation in this pdf file (see section 3, falling object with Newton drag):

http://math.kennesaw.edu/~sellerme/sfehtml/classes/math2202/fallingbodies.pdf

You can use numerical integration using the equation for acceleration versus velocity, to calculate velocity and position over a series of small time steps, or use the equation linked to above to calculate positions using the equation for velocity versus time. Example using acceleration versus velocity:

Acceleration = F(v) = g - ( Cd ρ A v2) / (2 m)

For each time step, you calculate

new_velocity = old_velocity + time_step x average_acceleration

One way to do this is to use trapezoidal method

average_accleration = 1/2 (old_acceleration + new_acceleration)

new_velocity = old_velocity + time_step x 1/2 (old_acceleration + new_acceleration)

where old_acceleration is the acceleration at the start of a time step and new_acceleration is the acceleration at the end of a time step. This requires you be able to calculate new_acceleration directly or to be able to predict it via an iterative method (see below)

You can then calculate position once you calculate the new velocity, using the same method:

new_position = old_position + time_step x 1/2 x (old_velocity + new_velocity)

An iterative corrector method will improve the results. In the algorithm shown below, an, vn, and pn, are successive "guesses" that should converge quickly. F(...) calculates the acceleration based on vn(t), and Δt is the elapsed time per step. You may want to do 6 to 8 interations instead of the 4 shown in this example. The first step is essentially Euler (since a1(t) is set = F(v0(t)) (= a(t-1)), the remaining steps are trapezoidal. Even though each step of this algorithm will converge to a specific set of values, the algorithm is based on trapezoidal rule, a linear approximation, so you need to use small time steps (Δt).

F(v) = g - ( Cd ρ A v2) / (2 m)

v0(t) = v(t-1)
p0(t) = p(t-1)

a1(t) = F(v0(t))
v1(t) = v(t-1) + 1/2 (a(t-1) + a1(t)) Δt
p1(t) = p(t-1) + 1/2 (v(t-1) + v1(t)) Δt

a2 = F(v1(t))
v2(t) = v(t-1) + 1/2 (a(t-1) + a2(t)) Δt
p2(t) = p(t-1) + 1/2 (v(t-1) + v2(t)) Δt

a3 = F(v2(t))
v3(t) = v(t-1) + 1/2 (a(t-1) + a3(t)) Δt
p3(t) = p(t-1) + 1/2 (v(t-1) + v3(t)) Δt

a4 = F(v3(t))
v4(t) = v(t-1) + 1/2 (a(t-1) + a4(t)) Δt
p4(t) = p(t-1) + 1/2 (v(t-1) + v4(t)) Δt

...

v(t) = vn(t)
p(t) = pn(t)
a(t) = F(pn(t))

time += Δt
t += 1

This is a predictor-corrector type algorithm:

http://en.wikipedia.org/wiki/Predictor-corrector_method#Euler_trapezoidal_example

Since you have an equation for velocity versus time v = V(t), you can use it directly and just use trapezoidal rule (time steps still need to be small):

t = 0
time = 0
v(t) = v(time)
p(t) = p(time)

loop:

t += 1
time += Δt
v(t) = V(time)
p(t) = p(t-1) + 1/2 (v(t-1) + v(t)) Δt
 
Last edited:
  • #3
Welcome to PF,

I agree that it looks nasty: right now you've got what's called an integral equation. The velocity depends on its own integral in some way. My advice would be NOT to take that step of integrating both sides of the equation. Instead, leave it as it is. Then you have a differential equation -- the velocity is related to its derivative. That's something that I, at least, have a better idea of how one might go about solving. For instance, in the x-direction where drag is the only force acting, you have:

[tex] F_x = -\frac{1}{2}\rho A C_d v_x^2 [/tex]

[tex] ma_x = -\frac{1}{2}\rho A C_d v_x^2 [/tex]

You can write this as

[tex] m\frac{dv_x}{dt} = -\frac{1}{2}\rho A C_d v_x^2 [/tex]

This differential equation is separable.
 
  • #4
rcgldr said:
v(t) = vn(t)
p(t) = pn(t)
a(t) = F(pn(t))

time += Δt
t += 1

a(t) is a function of velocity, so that third line should be:

a(t) = F(vn(t))
 
  • #5


I would like to commend you for your curiosity and initiative in creating a program that deals with 2D projectile motion with wind resistance. It is a fascinating and complex topic to explore.

Regarding your question about integrating velocity for velocity in the acceleration equation, I would like to offer some insights. The equation you have mentioned, a = CdPAV^2/2*m, is a simplified version of the drag force equation. In reality, drag force is not only dependent on velocity, but also on other factors such as the shape and size of the object, the density of the fluid, and the viscosity of the fluid.

In order to accurately model the motion of an object with wind resistance, you would need to take into account the changing drag force as the object moves through the fluid. This can be done by using differential equations and numerical methods to calculate the velocity and position of the object at different time intervals.

In addition, the integration of velocity for velocity is not a straightforward process and requires advanced mathematical techniques. Therefore, it may be more practical to use numerical methods to solve for the position and velocity of the object at different time intervals, rather than trying to derive a general position v. position equation.

I hope this helps to clarify your doubts. Keep exploring and learning about projectile motion and wind resistance – it is a fascinating and important topic in science and engineering.
 

1. What is velocity-dependent acceleration within drag force?

Velocity-dependent acceleration within drag force is a phenomenon in which the acceleration of an object is affected by the velocity at which it is moving through a fluid. This is due to the presence of a drag force, which is a resistance force that opposes the motion of the object.

2. How does drag force affect the acceleration of an object?

Drag force acts in the opposite direction of the object's motion and increases as the object's velocity increases. This means that as the object moves faster, the drag force increases and therefore, the acceleration decreases. This is known as velocity-dependent acceleration within drag force.

3. Why does drag force increase with velocity?

Drag force is caused by the interaction between the object and the fluid it is moving through. As the object's velocity increases, it pushes more fluid out of its way, resulting in a larger force in the opposite direction. This is why drag force increases with velocity.

4. How is velocity-dependent acceleration within drag force calculated?

The formula for velocity-dependent acceleration within drag force is a = (F - mg)/m, where F is the drag force, m is the mass of the object, and g is the acceleration due to gravity. This formula takes into account the decrease in acceleration due to the opposing drag force.

5. How does air resistance affect velocity-dependent acceleration?

Air resistance, which is a type of drag force, can significantly affect velocity-dependent acceleration. In cases where an object is moving through a fluid with a high air resistance, the drag force can be significant enough to cancel out the object's acceleration, resulting in a constant velocity or even a decrease in velocity over time.

Suggested for: Velocity Dependent Acceleration within Drag Force

Back
Top