POV-Ray : Newsgroups : povray.advanced-users : Air resistance Server Time
26 Nov 2024 04:41:58 EST (-0500)
  Air resistance (Message 1 to 10 of 13)  
Goto Latest 10 Messages Next 3 Messages >>>
From: Rune
Subject: Air resistance
Date: 15 Feb 2001 16:08:22
Message: <3a8c4546@news.povray.org>
How can I simulate air resistance to the movement of an object? If the
object is moving with the velocity V, what will happen to the velocity (or
speed) over time? In which way will it decrease? Linearly, exponentially, or
something different?

Just to check if my general understanding of the terms is correct: I'd say
that gravity increases the speed of an object linearly (assuming the object
is falling downwards). Is that correct?

Rune
--
\ Include files, tutorials, 3D images, raytracing jokes,
/ The POV Desktop Theme, and The POV-Ray Logo Contest can
\ all be found at http://rsj.mobilixnet.dk (updated January 28)
/ Also visit http://www.povrayusers.org


Post a reply to this message

From: Greg M  Johnson
Subject: Re: Air resistance
Date: 15 Feb 2001 16:29:54
Message: <3A8C48E9.1434CBAD@my-dejanews.com>
Rune wrote:

> Just to check if my general understanding of the terms is correct: I'd say
> that gravity increases the speed of an object linearly (assuming the object
> is falling downwards). Is that correct?

Say that
    dP=V*dT   (change in position is velocity times time increment)
    dV=a*dT     (change in velocity is acceleration times time increment.

So for every frame,
    P[t+1]=P[t]+V*constant
    V[t+1]=V[t]+a*constant

In my systems,
    #declare velocity=velocity-0.1*y;
is how I usually do gravity.


As for air drag, I see a formulat at
http://web.syr.edu/~smdemar/rocketdrag.html  which I would convert to

#declare velocity=velocity-dragconstant*velocity^2;

where dragconstant=0.005;   //for starters


Post a reply to this message

From: Christoph Hormann
Subject: Re: Air resistance
Date: 15 Feb 2001 17:07:25
Message: <3A8C531E.AE4DC9D5@gmx.de>
Rune wrote:
> 
> How can I simulate air resistance to the movement of an object? If the
> object is moving with the velocity V, what will happen to the velocity (or
> speed) over time? In which way will it decrease? Linearly, exponentially, or
> something different?
> 

It's usually much more complicated depending on the object's geometry for
example.  You can quite well approximate it with a quadratic formula for
the acceleration/force depending on the speed.  

> Just to check if my general understanding of the terms is correct: I'd say
> that gravity increases the speed of an object linearly (assuming the object
> is falling downwards). Is that correct?
> 

Yes, The gravitational force is m*g where g is the acceleration.  

You might check:

http://www.lerc.nasa.gov/WWW/K-12/airplane/bga.html
http://dmoz.org/Science/Technology/Aerospace/Aeronautics/Aerodynamics/

Christoph

-- 
Christoph Hormann <chr### [at] gmxde>
IsoWood include, radiosity tutorial, TransSkin and other 
things on: http://www.schunter.etc.tu-bs.de/~chris/


Post a reply to this message

From: Francois Labreque
Subject: Re: Air resistance
Date: 15 Feb 2001 20:13:03
Message: <3A8C7E00.664C3EA2@videotron.ca>
Rune wrote:
> 
> How can I simulate air resistance to the movement of an object? If the
> object is moving with the velocity V, what will happen to the velocity (or
> speed) over time? In which way will it decrease? Linearly, exponentially, or
> something different?

Well, a few years ago, some english chap name Newton came up with the
silly idea that a force apllied to an object resulted in a acceleration
proportional to its mass, of f=ma for short.

the force caused by air resistence is computed by the following formula:

F = (Cd)A(V^2) in the direction opposing the movement.

Where:
- Cd is the drag coefficient ( 1 for a box, 0.5 for a sphere and around
0.30 for most modern cars)
- A is the cross-sectional area of the object (area "facing" the air
flow)
- V is the speed of the object in that direction.

Provided that the shape and mass of your object dont change, you can see
that the force, hence the decelleration, will decrease at a rate which
is proportional to the square of the current speed.

> 
> Just to check if my general understanding of the terms is correct: I'd say
> that gravity increases the speed of an object linearly (assuming the object
> is falling downwards). Is that correct?

For all intents and purposes, let's say that the force of gravitation is
constant (the mass and size of your object are negligeable compared to
the mass and size of the Earth) and that the mass of your object is also
constant, hence the acceleration will be constant.  On Earth, it is 9.81
m/s^2 (or 32.2 ft/s^2 for the metric-impaired)

So yes, the downward velocity of your falling object will increase
linearly.

-- 
Francois Labreque |   //\\    Wear an ASCII ribbon!
    flabreque     |  ||  ||   
        @         |   \\//    Support the campain
   videotron.ca        \\     against HTML e-mail
                      //\\    and news!


Post a reply to this message

From: Christoph Hormann
Subject: Re: Air resistance
Date: 16 Feb 2001 02:57:53
Message: <3A8CDD82.2AFF0279@gmx.de>
Francois Labreque wrote:
> 
> Well, a few years ago, some english chap name Newton came up with the
> silly idea that a force apllied to an object resulted in a acceleration
> proportional to its mass, of f=ma for short.
> 
> the force caused by air resistence is computed by the following formula:
> 
> F = (Cd)A(V^2) in the direction opposing the movement.
> 
> Where:
> - Cd is the drag coefficient ( 1 for a box, 0.5 for a sphere and around
> 0.30 for most modern cars)
> - A is the cross-sectional area of the object (area "facing" the air
> flow)
> - V is the speed of the object in that direction.
> 
> Provided that the shape and mass of your object dont change, you can see
> that the force, hence the decelleration, will decrease at a rate which
> is proportional to the square of the current speed.
> 

You should mention, that Cd is in fact not a constant but is a highly
complicated function of the reynolds number and thereby also depends on
the speed.  Values given for cars etc. are usually only valid for a
limited range of speed.  

BTW with some shapes (like the wings of a plane) it is even more
difficult, because the force no more only goes in direction opposite to
the air flow.

Christoph

-- 
Christoph Hormann <chr### [at] gmxde>
IsoWood include, radiosity tutorial, TransSkin and other 
things on: http://www.schunter.etc.tu-bs.de/~chris/


Post a reply to this message

From: Warp
Subject: Re: Air resistance
Date: 16 Feb 2001 09:38:01
Message: <3a8d3b48@news.povray.org>
You have to take into account that air resistance can cause uneven forces
in different parts of the object (for example if the density of the object is
not constant along its volumen), which can cause the object to rotate.

-- 
char*i="b[7FK@`3NB6>B:b3O6>:B:b3O6><`3:;8:6f733:>::b?7B>:>^B>C73;S1";
main(_,c,m){for(m=32;c=*i++-49;c&m?puts(""):m)for(_=(
c/4)&7;putchar(m),_--?m:(_=(1<<(c&3))-1,(m^=3)&3););}    /*- Warp -*/


Post a reply to this message

From: Greg M  Johnson
Subject: Re: Air resistance
Date: 16 Feb 2001 09:46:35
Message: <3A8D3BD9.DADAC8CF@my-dejanews.com>
"Greg M. Johnson" wrote:

> In my systems,
>     #declare velocity=velocity-0.1*y;
> is how I usually do gravity.

IMHO I think my advice is the best here for frame-by-frame pov coding of gravity.

> #declare velocity=velocity-dragconstant*velocity^2;
> where dragconstant=0.005;   //for starters

Oops!  an error above.

//   for a frame-by-frame system, experiment with small numbers for drag; toss
textbook values in dustbin.
#declare dragconstant=0.005;

//   Wind is some varying vector.
#declare wind=foo;

//   Wind blowing at craft means double effect; blowing with craft negates
effect.
#declare netdirection=wind-velocity;

//   I forgot that velocity was a vector in last version I shared.
#declare
velocity=velocity+vnormalize(netdirection)*dragconstant*vlength(netdirection)^2;


Post a reply to this message

From: Francois Labreque
Subject: Re: Air resistance
Date: 16 Feb 2001 10:07:31
Message: <3A8D418F.C73827C3@videotron.ca>
Warp wrote:
> 
>   You have to take into account that air resistance can cause uneven forces
> in different parts of the object (for example if the density of the object is
> not constant along its volumen), which can cause the object to rotate.
> 

It's got nothing to do with wind resistence per se, but with the uneven
density of the object.  the Force IS constant over the area, but since
the density isn't, the decelleration won't be, so you have to divide you
surface into smaller subsurfaces and divide the force applied to each
individual subsurface to take into account the change of density.

Of course, if your shape starts to rotate, its cross-sectional area will
also change, hence the Cd will also chance and the wind resistence will
also change.

Welcome to the fun world of finite element analysis.


> --
> char*i="b[7FK@`3NB6>B:b3O6>:B:b3O6><`3:;8:6f733:>::b?7B>:>^B>C73;S1";
> main(_,c,m){for(m=32;c=*i++-49;c&m?puts(""):m)for(_=(
> c/4)&7;putchar(m),_--?m:(_=(1<<(c&3))-1,(m^=3)&3););}    /*- Warp -*/

-- 
Francois Labreque |   //\\    Wear an ASCII ribbon!
    flabreque     |  ||  ||   
        @         |   \\//    Support the campain
   videotron.ca        \\     against HTML e-mail
                      //\\    and news!


Post a reply to this message

From: Rune
Subject: Re: Air resistance
Date: 16 Feb 2001 10:59:03
Message: <3a8d4e47@news.povray.org>
Thank you all! Your replies were all helpful.

I needed just a simple approximation, nothing complicated. And what I need
to simulate it on is not really objects, but different kinds of particles
such as smoke, water drops, sparks and other things a particle system could
be used for.

I suppose that doesn't change anything; the basic drag formula should still
work the same, right? It's just a matter of finding the right drag
coefficient, drag constant or whatever you want to call it.

Rune
--
\ Include files, tutorials, 3D images, raytracing jokes,
/ The POV Desktop Theme, and The POV-Ray Logo Contest can
\ all be found at http://rsj.mobilixnet.dk (updated January 28)
/ Also visit http://www.povrayusers.org


Post a reply to this message

From: Warp
Subject: Re: Air resistance
Date: 16 Feb 2001 11:16:53
Message: <3a8d5273@news.povray.org>
Francois Labreque <fla### [at] videotronca> wrote:
: It's got nothing to do with wind resistence per se, but with the uneven
: density of the object.

  The density of the object can be constant thorough its volume, but the
shape of the object can affect the air resistance so that it makes the
object rotate.
  A windmill can have a constant density. It rotates. :)

-- 
char*i="b[7FK@`3NB6>B:b3O6>:B:b3O6><`3:;8:6f733:>::b?7B>:>^B>C73;S1";
main(_,c,m){for(m=32;c=*i++-49;c&m?puts(""):m)for(_=(
c/4)&7;putchar(m),_--?m:(_=(1<<(c&3))-1,(m^=3)&3););}    /*- Warp -*/


Post a reply to this message

Goto Latest 10 Messages Next 3 Messages >>>

Copyright 2003-2023 Persistence of Vision Raytracer Pty. Ltd.