1 Computational Modeling for Engineering MECN 6040Professor: Dr. Omar E. Meza Castillo Department of Mechanical Engineering

2 Ecuaciones lineales homogéneas con coeficientes constantes

3 Ecuaciones lineales homogéneas con coeficientes constantesdonde ai son constantes, an  0. Ecuación o polinomio auxiliar : Para n = 2, Si probamos y(x) = emx, obtenemos la ecuación auxiliar.

4 Las dos raíces del polinomio auxiliar son: (1) b2 – 4ac > 0: reales y distintas, m1  m2 (2) b2 – 4ac = 0: reales e iguales, m1 = m2 = -b/(2a) (3) b2 – 4ac < 0: complejas conjugadas,

5 Caso 1: Raíces reales y distintasLa solución general es Caso 2: Raíces reales repetidas

6 Caso 3: Raíces complejas conjugadas EscribimosCaso 3: Raíces complejas conjugadas Escribimos , una solución general es Usando la fórmula de Euler:

7 Como es solución general, tomando C1=C2=1 y C1 =1, C2 =-1Tenemos dos soluciones: Así, ex cos x y ex sen x son un conjunto fundamental de soluciones y la solución general es:

8 Resolver las EDs siguientes:(c)

9 Resolver Solución:

10 Resolver las ecuaciones:Para la primera ecuación : Para la segunda ecuación : Como Luego

11 Separación de VARIABLESSOLUCIÓN DE ECUACIONES DIFERENCIALES PARCIALES MÉTODO ANALÍTICO

12 Separación DE VariablesAsumiendo que u=X(x)Y(y), entonces

13 Determine las soluciones producto de:Solución Sea u=X(x)Y(y) y entonces Introducimos una constante de separación real como 2.

14 Tenemos tres casos: Caso I: 2 > 0 X” – 42X = 0, Y’ − 2Y = 0

15 Caso II: -2 < 0 X” + 42X = 0, Y’ + 2Y = 0

16 Caso III: 2 =0 X” = 0, Y’ = 0

17 Ecuación de Calor La ecuación de calor puede describirse así:(1) (2) (3)

18 Ecuación de Calor f(x) T x Distribución de temperatura a lo largo de la barra en un instante de tiempo cualquiera u=0 u=0 L

19 Ecuación de Calor Usando u(x,t) = X(x)T(t), y −2 como la constante de separación: (4) (5)

20 Ecuación de Calor Resolviendo: (6) (7) Condiciones de Frontera

21 Cuando las condiciones de frontera X(0) = X(L) = 0 se aplican a (4), estas soluciones son sólo X(x) = 0. Aplicando la primera condición a (6) se obtiene c1=0. Por tanto X(x)=c2sinx. La condición X(L) = 0 implica que X(L)=c2sinL=0. Tenemos que sinL=0 para c2  0 y =n/L, donde n=1, 2, 3, … y las soluciones correspondientes

22 La solución general de (7) esy por tanto (10) donde An = c2c3.

23 Ahora usando las condiciones iniciales u(x,0)=f(x), 0 < x < L, tenemos(11) Por el principio de superposición, la función (12) Debe cumplir (1) y (2). Si ponemos t=0, entonces

24 Se conoce como un desarrollo de semi-intervalo para f en una serie seno. Si ponemos An = bn, n= 1, 2, 3, … entonces: Llegamos a la conclusión de que la solución del PVF descrito por (1), (2) y (3) se expresa mediante la serie infinita

25 En el caso especial, donde u(x,0)=100, L = , y k = 1, entoncesDe modo que la solución es:

26 Mathcad

27 Otro caso La ecuación de calor puede describirse así: (1) (2) (3)

28 En el caso especial, donde L = , y k = 1, entoncesDe modo que la solución es:

29 Mathcad

30 Ecuación de LAPLACE Considere el siguiente problema de valores en la frontera: (1) (2) (3)

31 T2 L T1 T1 W

32

33

34

35

36

37

38

39 Mathcad