Computational Modeling for Engineering MECN 6040

1. Computational Modeling for EngineeringMECN 6040 Professor: Dr. Omar E. Meza Castillo omeza@bayamon.inter.edu http://facultad.bayamon.inter.edu/omeza Department of Mechanical Engineering

2. Finite difference Methods for Hyperbolic problems Wave Equation

3. Partial Differential Equations • Our study of PDE’s is focused on linear, second-order equations • The following general form will be evaluated for B2 - 4AC (Variables – x and y/t)

4. Partial Differential Equations B2-4AC Category Example < 0 Elliptic Laplace equation (steady state with 2 spatial dimensions) = 0 Parabolic Heat conduction equation (time variable with one spatial dimension) >0 Hyperbolic Wave equation (time-variable with one spatial dimension)

5. WAVE Equation Mechanical vibrations of a guitar string, or in the membrane of a drum, or a cantilever beam are governed by a partial differential equation, called the Wave Equation. To derive the wave equation we consider an elastic string vibrating in a plane, as the string on a guitar. Assume u(x,t) is the displacement of the string away from its equilibrium position u=0.

6. WAVE Equation u u(x,t) x x x+Dx T(x+Dx,t) b a T(x,t)

7. WAVE Equation T(x+Dx,t) b Newton’s Law: F=ma Vertical: where ρ is the mass density function of the string Horizontal: As we assume there is only vertical motion of the string a T(x,t)

8. WAVE Equation Let Then, So, Fvert = ma implies, from the previous slide, that Or,

9. WAVE Equation T(x+Dx,t) b Now, tan β is just the slope of the string at x+Δx and tan a is the slope at x. Thus, So, a T(x,t)

10. WAVE Equation Now, take the limit as Δx -> 0, to get This the Wave Equation in one dimension (x)

11. Wave EquationBoundary CONDITIONS Since the wave equation is second order in space and time, we need two boundary conditions for each variable. We will assume the string x parameter varies from 0 to L. Typical set-up is to give boundary conditions at the ends of the string in time: u(0,t) = 0 u(1,t) = 0 and space initial conditions for the displacement and velocity: u(x,0) = f(x) ut(x,0)=g(x)

12. Wave EquationBoundary CONDITIONS

13. One DimensionalWave Equation • We will solve this equation for x and t values on a grid in x-t space: Approximate solution uij=u(xi, tj) at grid points

14. Finite Difference Approximation • To approximate the solution we use the finite difference approximation for the two derivatives in the wave equation. We use the centered-difference formula for the space and time derivatives:

15. Finite Difference Approximation • Then the wave equation (utt =c2uxx ) can be approximated as Or, Let p = (ck/h) Solving for ui,j+1 we get:

16. One Dimensional Wave Equation • Note that this formula uses info from the j-th and (j-1)st time steps to approximate the (j+1)-st time step:

17. One Dimensional Wave Equation • The solution is known for t=0 (j=0) but the formula uses values for j-1 , which are unknown. So, we use the initial condition ut(x,0)=g(x) and the centered difference approximation for ut to get • At the first time step, we then get

18. One Dimensional Wave Equation • Simplifying, we get • So,