Introduction

Objectives
  • Formulate mathematical models based on scientific principles for physical systems.

  • Utilize numerical methods for solution generalization applicable on digital computers.

  • Understand conservation laws in engineering models; recognize steady-state vs dynamic solutions.

  • Explore various numerical methods.

Mathematical Model Overview
  • Mathematical models express the essential features of a physical system.

  • Represented by relationships involving dependent variables, independent variables, parameters, and forcing functions.

Components of Model Functions
  • Dependent variable: Reflects the behavior/state of the system.

  • Independent variables: Dimensions like time and space for system behavior analysis.

  • Parameters: Constants defining system properties.

  • Forcing functions: External influences affecting the system.

Standard Form of Model Functions
  • Generic form for a first-order differential equation:
    dydt=f(t,y,p,F)\frac{dy}{dt} = f(t, y, p, F)

  • yy: Dependent variable (e.g., temperature, velocity)

  • tt: Independent variable (e.g., time, position)

  • pp: Parameters (e.g., system constants)

  • FF: Forcing functions (e.g., external influences)

Analytical vs Numerical Solutions
  • Analytical solutions: Involve exact formula derivations for the state of the system, often expressed as y(t)y(t).

  • Numerical solutions: Use iterative calculations, such as discrete time-stepping, to approximate system behavior.

Euler's Method
  • A numerical technique for approximating solutions of differential equations by iterating calculations in finite intervals.

  • Standard Formula:
    v<em>i+1=v</em>i+(slope×step)v<em>{i+1} = v</em>i + (\text{slope} \times \text{step})

  • Iterative Process: Calculate the next value by adding the product of the estimated slope and the step size to the current value.

Case Studies in Standard Form
  • Newton's Law of Cooling:
    dTdt=k(TTa)\frac{dT}{dt} = -k(T - T_a)

  • Calculates the rate of temperature change relative to ambient temperature TaT_a.

  • Evaporation of Liquid Droplet:
    dVdt=kA\frac{dV}{dt} = -k A

  • Models volume reduction as proportional to the surface area AA of the droplet.

  • Fluid in Storage Tank:
    dhdt=Q<em>inQ</em>outA\frac{dh}{dt} = \frac{Q<em>{in} - Q</em>{out}}{A}

  • Solves for changes in liquid depth hh based on inflow Q<em>inQ<em>{in} and outflow Q</em>outQ</em>{out} rates relative to cross-sectional area AA.