HOME

HW 2 Mathematics 503, Mathematical Modeling, CSUF June 2 2007

Nasser M. Abbasi

June 15, 2014

Contents

1 Problem 1 (section 1.2,#9, page 33)
2 Problem 2 (section 1.1, #11, page 18)

1 Problem 1 (section 1.2,#9, page 33)

problem:

PIC

solution:

Let [Θ ]  be the dimension of temperature.

A small note on how to scale the dependent variable and the independent variable in the problem.

For the dependent variable, say y(t)  , if the problem does not say which problem parameters to scale against, then look for a parameter or a combination of these parameters with the appropriate dimension of that of the dependent variable such the parameter is the largest value. This way we are measuring the depndent variable relative to the largest parameter in the problem. Hence when we write ~y = ys  where s is the largest parameter in the problem with dimension of y.

For the indepenent variable (typically this will be time), we look for the larget rate of change in the problem involving this variable, and scale the time relative to this rate. Hence we write    t
~t = s  where now s is the smallest parameter in the problem with dimension of t .

Now we start with the solution of the problem.

(a) dT= q×
dt number, hence

     [  ]
      Θ-
[q ]=  T

i.e. the dimension of q is temperature over time.

Now to find the dimension of k. sincedT
dt = k(T − Tf)  , then [Θ]
 T  = [k][Θ]   hence

     [1 ]
[k]=  T-

Now to find the dimension of A . Since  − A
e T   must be a numerical dimensionless value, then A must have the same dimensions as T (which has dimension [Θ]  ), hence

[A]=  [Θ ]

(b) The constants in this problem are Tf,T0,q,k,A , and the dependent variable is T which is the temperature of the sample, and the independent variable is t which is time.

dT-     −TA
 dt = qe  − k(T − Tf)

Hence to reduce the above ODE to dimensionless form, we need to transform the dependent and the independent variables to dimensionless variables.  Start with the dependent variable T .

We write      T
T~=  s  , so we need to find s as a combination using the constants Tf,T0,q,k,A in the ODE which would have the dimension of temperature to make ~T dimensionless.

This is easy, since we are told to use Tf  , hence s = Tf  and so

     T
T~= --
    Tf
(1)

Now we need to transform the independent variable, which is time t , hence we write ~  t
t = s  , and we need to find s as a combination of the constants Tf,T0,q,k,A which has a dimension of time.

Looking above at the dimensions of these constants, we see that the following combinations are possible: {         }
  1,A,Tf,T0- .
  k q q  q Now since we want the new time scale to be large when the heat loss is large, and since the heat loss is proportional to k , hence we pick s= 1k  , this will make s small and st  will then be large. Therefor

~
t = kt
(2)

Now we convert the ODE to non dimensional using (1) and (2)

dT-=  dT-d~t
 dt   d~t dt

but from (2) we have d~t
dt = k , hence

dT    dT
---= k---
dt     d~t

and since T = T ~T
     f , then the above becomes

dT      dT~
---= Tfk---
dt       d~t

Hence the ODE dT     −A
-dt = qeT − k(T − Tf)  becomes (after replacing T by TfT~ in the RHS)

pict

Now -q-
Tfk  must be dimensionless. Verify -[ΘT]-
[Θ][1T] = 1  . OK, so the above ODE can be written as

|----------------------|
|  ~     -−A~  (     )  |
| dTd~t = βeTfT − ~T − 1   |
-----------------------|

where β = Tqk-
     f

Note  −A-
eTf~T   must be dimensionless also, since T−fA ~T  is dimensionless.

For the initial conditions, now it becomes T~=  T-
 o   Tf

2 Problem 2 (section 1.1, #11, page 18)

problem:

PIC

solution:

First draw a free body diagram. Assume the ball is moving upwards, and assume the drag force is given by        ′   2
Fd = k(y (t))   where k is constant of proportionality, and weight of ball is given by W  = mg

PIC

Now to obtain the equation of motion, apply Newtons second law F = ma , hence we obtain

   ( dy(t))2         d2y(t)
− k  -----  − mg = m ---2--
      dt              dt

Notice that we assumed the positive y direction is upwards. Hence we obtain the equation of motion

           (      )
  d2y(t)     dy(t) 2
m   dt2  + k   dt    + mg = 0

with the initial conditions

         dy-(0-)
y(0) = y0,  dt  = V

This problem has the following constants: {V,m,g,k,y}
          0 . To find the dimensions of k , since  (    )
k  dy(t) 2
   dt   must have the dimension of Newton [N]  then we write

     [ (    ) ]
         dy(t) 2    [  L-]     [L]2
[N ]=  k   dt    →   M T2 = [k] T

    [   ]
      M
[k]=   L-

and     [ L]
[g]=   T2 , and [m ]= [M ]  , and      [ L]
[V]=   T , and y0 = [L]  , hence we write

           {     [L ]             [ L ]    [M ]        }
constants →   V =  -- ,m = [M],g =  -2- ,k = --  ,y0 = [L ]
                  T                T         L

Now we can start scaling the ODE.

The dependent variable is y(t)  , and the independent variable is t . Hence we write y~= ys  , where [s]= [L]  , so using the above constants, we need to find a combination which has the dimension of length [L ]  . Clearly y
0   is one possibility. There are other combinations which give the dimension of length. We find the following: {   m V2}
 y0,k,g- .

Now we get to the hard part of these scaling problems. Which combination to choose? The problem did not give us a hint on this. I choose s= V2
    g  , therefor

~y=  y-g
    V2

Now for the time scale. We write ~t = ts  where [s]= [T]  , looking at the above constants, we need a combination with dimension [T]  , we see that s= V
   g  has the dimension of time, so we have

~   t-
t = V g

Now that we have the scaling completed, we apply them to scale the ODE. In otherwords, we need to rewrite the ODE so that instead of using{y,t} it will use {~y,~t}

dy= dyd~t
dt  d~tdt  , but d~t= -g
dt  V  , hence

dy-  g-dy-
dt = V d~t

and d2y2 = d(dy) = d(dy) d~t=  d(g-dy) g
dt   dt  dt    d~t dt  dt   d~t Vd~t V  hence

      (  )
d2y    -g 2 d2y
dt2 =  V    d~t2

Hence the scaled ODE is

pict

Now replace      2
y→  Vg-~y

pict

Lets us now verify the above is dimensionless. Looking at [    ]              [  ]           [  ][  ]
 -kV2 →  [M][-1][L]2 -1L  →  [M][-1] L22  T2
 m g      L  M   T    T2-     L  M   T    L

Hence [-kV2]→  1
 m g  (like magic)

So the dimensionless ODE is

d2y~   ( dy~)2
--2 + β  ---  + 1=  0
 d~t      dt~

where       2
β = kmVg

Now for the initial conditions in terms of the new scaled variables. Since

y(0) = y0,dy-(0-)= V
           dt

hence

       y
y~(0) = -02g
       V

and      dy(0)   dy(0)d~t  V2d~y(0)d~t
V =  -dt-=  -d~t-dt =-g-d~t-dt  , but d~t   g
dt = V  hence

dy(0)-= V2d~y(0)~-g→  V d~y(0~)
 dt   g  dtV      dt

So d~y(0)-  1dy(0)
 dt~ = V dt  , but dy(0)
 dt = V , therefor

dy~(0)
--d~t- = 1

So the equation of motion in scaled dimensionless form is

|--------------------|
|       (  )         |
| d2~y2 + β d~y 2+ 1= 0  |
--d~t-----d~t-----------

with initial conditions

|------------|
|       y    |
| ~y(0)= V02g  |
--------------

and

|----------|
| d~y(0)     |
| -d~t-= 1  |
------------