tiales. Ecuaciones en Derivadas Par 1. Introduccion. ). por su semejanza con Jas ODE, es una ‘nicht por unit relacion entre sus derivadas Una ccuacién d ecuacidn donde una cierta funcidn med parciales con respecte a tas variables independicnies cencial en derivadas parciates (PI) nita w viene de seria Si wh u(xy.z), une ccuacide diferencial en derivadhs parciale: Se denomina arden de la PDE al mas alto geaclo de derivacién pareial que aparece en la exptesion, Asi(1) es una PDE de 2 orden. mientras que 2) es na PDE de primer orcen La eouacién (1) ¢s lineal ya que wy sus derivackis aparecen sin multipticarse y no apareecn clevadas a potencias, 1a ccuncidn (2) es, cn cambio, no Hincal una funcidn de (xy) que verifies que es de ba farina uOuy) = xb vets = ¥) donde Fy g son funcianes arbitrarias, Lintonves asu verau = flig-+2-ftig! 4 Fag” a | f Bites! SY fly 2 Pg eb gt a ftegt = Big —Pogt de donde se deduce que Pero -(3)-@) n diferencial de 2° orden cn derivadas parciales cuya solucién tiene la Sc ha Hlegado a una ccus forma uGxy) © Rx yy B(x = y) Considérese otro ejemplo es la ecuacién diferencial de primer orden cuya soluci6n tiene 1a forma. u=xity)donde f es une funcién arbitraria. Otro ejemplo es el siguiente: Otro ejemplo u ox ei au Php} ay Ss ®u = 20 cots FY co gat am Pu oe ee Repiiese que segiin los resutados obtenidos existen infinitas soluciones posibles de la PDE. Pero ahora la arbitrariedad de la solucién general viene dada en términos de funciones, apareciendo tantas como el orden de la ecuacién. Desde el punto de vista de la Matematica pucde parecer mis preciso obtener en cualquier caso la solucién general, sin embargo, se van a buscar soluciones dentro del campo de la Fisica por lo que solo interesara una solucién particular concreia, Estas soluciones particulares van a satisfacer unas determinadas condiciones de contorno y de valor inicial Es decir, se va a tratar de obtener Ia solucién de una cierta PDE que verifique unas condiciones en el contorno del dominio en que esta definida (condiciones de contorno), y si ademas una variable es cl tiempo "t" las condiciones en t=O se darén como dato (condiciones inticiales)Algunos ejemplos tipicos de eeuaciones en derivadas parciales son Ecuacién de difusion: Es [a cldsica ecuacion unidimensional de difusisn del calor, de segundo orden, tineal, homogénea ¥ de cocficientes constantes. Ecuacién de onda: Ecuacién de Laplace Esta es una ecuacion bidimensional, de seguro orden, lineal, homagénea y de coeficientes constantes, describiendo poteneiaics eléetricos © gravitatorios 9 procesos de difusién cn los que se ha alcanzado un equilibrio térmico, Ecuaciin de Po Pu ae Es también una ecuacion bidimensional, de segundo orden, lineal, de coeficientes constantes, pero no homogénea. 2, FCUACIONES DE SEGUNDO ORDEN LINEALES HOMOGENEAS DE COEFICIENTES CONSTANTES ‘Una ecuaci6n de 2 orden en derivadas parciales lineal, homogenea, con coeficientes constantes tiene Ja forma (suptiesto dos variables indepenicntes): donde, a,h,b,£g,c son constantes,se puede decir que estas ecuaciones se clasifican > vliplicas, patabdlicas ¢ hiperbéticas Esto es, si aero bn e5 eliptica: z vt ta ecuacién es parabatien: ab-W oo fa eovacion es hiperbiilica asiens ecuaciones de difusidn, de ondas y de Laplace perienceen a los tipos Segim esto, Vas el Ecuacién de difusién: parabolica Ecuacién de onda: hiperbolica Ecuacidn de Laplace: elipticn fasificaciin sigue sicnde vdliha mclise cunts los coeficiettes de ta ecuacion a, h, dewey, Errestos casos fe condciin puede cantbiar ee tipo al pasar de wn Nos: Este «280 flaiciones variahles Cuadlrante «otra, Par egenipte fa eenacrin y oxy Oy = 0, ¢ hiperbélica ay es eliptica en la regién ¥ ~* > 0, parabélica a lo largo de las rectas ¥ 2 enlaregisn? ~¥® <0,CAPITULO XvItr ECUACIONES DE LA FISICA MATEMATICA § 1, TIPOS FUNDAMENTALES DE LAS ECUACIONES DE LA FISICA MATEMATICA Se llaman ecuaciones fundamentales de la fisica matematica (para el caso de unas funciones de dos variables independientes) a las siguientes ecuaciones diferenciales con derivadas parciales de segundo orden. I. Ecuacién de onda: (1) A la necesidad de analizar esta ecuacién conduce el examen de los procesos de vibraciones transversales de una cuerda, vibraciones longitudinales del vdstago, oscilaciones eléctricas en un conductor, oscilaciones torsionales de un Arbol, oscilaciones de un gas, etc. Esta ecuacién es la mas sencilla de tipo hiperbélico, II. Ecuacién de conduccién de calor o ecuacién de Fourier: ep Fu at ace @ A la necesidad de analizar esta ecuacién conduce el examen de los procesos de propagacién de calor, filtracién de liquido y gas en un medio poroso (por ejemplo, filtracién de petréleo y gas en are- niscos subterréneos), algunos problemas de la teoria de probabili- dades, etc. Esta ecuacién es la més sencilla de tipo parabélico. III. Ecuacién de Laplace: fu, ou etait306. Ecuaciones de Ia fisica matemdtica El examen de los problemas sobre campos eléctricos y magné- ticos, sobre el estado térmico estacionario, problemas de la hidro- dindmica, difusion, etc., plantea la necesidad de analizar esta ecuacién. Esta ecuacién es la més sencilla de tipo eliptico. En las ecuaciones (1), (2) y (3) la funcién desconocida u depende de dos variables. Se estudian también las ecuaciones correspondien- tes para unas funciones con mayor numero de variables, Asi, la ecuacién de onda con tres variables independientes tiene la forma: Pua f(Fu, Fu oe Geta (1) at az" ay’ la ecuacién de conduccién de calor con tres variables independientes tiene la forma: ou (24 a) re eee 2 at aa ay? (2) la ecuacién de Laplace con tres variables independiente tiene la forma: Fu, Fu tH (3) at" ay’ § 2. ECUACION DE OSCILACIONES DE UNA CUERDA. FORMULACION DEL PROBLEMA CON VALORES DE CONTORNO. ECUACIONES DE OSCILACIONES ELECTRICAS EN LOS CONDUCTORES En la fisica matematica bajo el término de cuerda se entiende un hilo flexible y eldstico. Las tensiones que surgen en la cuerda en cualquier momento, estan dirigidas por la tangente a su perfil. Supongamos que en el momento inicial la cuerda de longitud | esta dirigida a lo largo del segmento del eje Oz desde 0 hasta 1. Supon- gamos también que los extremos de la cuerda estan fijados en los puntos z = 0 y z = 1. Si desviamos la cuerda de su posicién inicial y luego la dejamos en libertad, o bien, sin desviar la cuerda, si damos en el momento inicial a sus puntos cierta velocidad, o si desviamos la cuerda y damos a sus puntos cierta velocidad entonces los puntos de la cuerda comienzan a realizar unos movimientos y se dice que la cuerda comienza a vibrar, El ptoblema consiste en la determinacién de la forma de la cuerda en cualquier momento y en la determinacién de la ley de movimiento de cada punto de la cuerda en funcién del tiempo. Analicemos pequefias desviaciones de los puntos de la cuerda a partir de la posicién inicial. Por tanto, se puede suponer que elEcuacién de oscilaciones de una cuerda 367 movimiento de los puntos de la cuerda se efecttia perpendicularmente al eje Ox y en un mismo plano. Bajo esta suposicién el proceso de oscilacién de la cuerda se describe per una funcién u (z, ft), que da la magnitud de desplazamiento de un punto de la cuerda de abscisa + en el momento ¢ (fig. 371). Puesto que examinamos pequefias desviaciones de la cuerda en el plano (z, u), supongamos que la longitud del elemento de la cuerda — M,Mz es igual a su proyeccién ‘sobre el eje Oz, 0 sea *), ~ MM, = 2, — 2;. Supongamos tambiéA que la tensién es igual en todos los puntos de la cuerda; designemos esta tension por T. Bue z Fig. 372 Analicemos el elemento de la cuerda MM" (fig. 372). En los extremos de este elemento actian las fuerzas 7 dirigidas por las tangentes a la cuerda. Supongamos que lds tangentes forman con el eje Oz los angulos p y @ + Ag. Entonces, la proyeccién sobre el eje Ou de las fuerzas que acttian sobre el elemento MM’ sera igual a Tsen (p + Ag) — Tsen p. Puesto que el angulo @ es peque- *) Esta suposicién equivale al despreciamiento de la magnitud uk en comparacién con 1. En efecto, Ps 4 Jers iy deany—7y, 1 ~ matin) Vinee f (t+ A m4368 Ecuaciones de la fisica matemética fio, podemos poner tg p = sen @, entonces tenemos: T sen (p+ Ag) — Tseng ~ HTiglot Ao) —Tegg ar | Mets ) du(z, 4] ax ex 4 Fu(x+0Az, t) Fu(z, yy = 7 are TOY, az * ae SS o 0 (luego ue 0372 Ecuaciones de la fisica matemética examinemos también el caso de 4 <0). Pues, De estas igualdades obtenemos dos ecuaciones: X"44X =0, (8) T’ + a2T = 0. (9) Las soluciones generales de estas ecuaciones serdn: (véase cap. XIII, § 21): X (x) = Acos Vax + Bsen Viz, (40) T (&=Ccosa Vit + DsenaVit, (11) donde A, B, C, D son constantes arbitrarias. Sustituyendo las expresiones X (x) y 7 (t) en la igualdad (6), obtenemos u(z, t) = (A cos Viz + Bsen Viz) (Cosa Vit + D senaVit). Escojamos ahora las constantes A y B de modo que se cumplan las condiciones (2) y (3). Puesto que 7 (t) #0 (en caso contrario u (z, t) = 0 lo que contradice a la condicién planteada), la funcién X (z) debe satisfacer a las condiciones (2) y (3), es decir, ha de ser X (0) = 0, X (J) = 0. Poniendo los valores z= 0 yz=Jen la igualdad (10), en virtud de (2) y (3) obtenemos: O=A-1+8B-0, O-= AcosVAl+ BsenVAl=0. De la primera ecuacién hallamos que A = 0. De la segunda obte- nemos: BsenVAl=0. B +0, puesto que en caso contrario seria X = 0 y u=0, lo que contradice a la condicién. Por consiguiente, debe ser: sen VAl= 0, de donde (n0 tomamos el valor x = 0, puesto qué en este caso seria X =0 y u = 0). Pues, hemos obtenido X=Bsen™e, (13)Soluctén de la ecuacién de vibraciones de una cuerda 373 Los valores hallados de % se laman propios para el problema dado con valores de contorno. Las funciones X (z), que les corresponden, se Ilaman funciones propias. Observacién. Si tomamos en vez de — i la expresién + 4 =F, entonces la ecuacién (8) adquiere la forma: X" — BX =0. La solucién general de esta ecuacion: X= Ae + Be™, La solucién en semejante forma que difiere de cero no puede satis- facer a las condiciones de contorno (2) y (3). Sabiendo VX y aprovechando la igualdad (11) podemos escribir: T(t) =Coos "t+ Dsen Te (n=4, 2, ...). (14) Para cada valor de n, por consiguiente, para cada A, ponemos las expresiones (13) y (14) en la igualdad (6) y obtenemos la solu- cién de la ecuacidn (1) que satisface a las condiciones de contorno (2) y (3). Designemos esta solucin con u, (z, t): Un (x, pasen™e(c, cos t+ Da sen 1). (15) Para cada valor de n podemos tomar sus constantes C y D, y por eso escribimos C, y D, (la constante B esta incluida en C, y D,). Puesto que la ecuacién (1) es lineal y homogénea, entonces la suma de solu- ciones también es una solucién, y por eso la funcién representada por la serie u(z, )= 3 Un (z, t) n= ulz, => (ccos 2 + Dy sent) son, (18) también sera la solucién de la ecuacién diferencial (1) que satisfa- cera a las condiciones de contorno (2) y (3). Es evidente, que la serie (16) seré la solucién de la ecuacién (1) solamente en el caso en que los coeficientes C, y D, son tales que esta serie converge y_convergen las series que se obtienen después de la derivacién doble término a término respecto a z y a t.314 Ecuaciones de ta fisica matemética La solucién (16) debe, ademés, satisfacer a las condiciones ini- ciales (4) y (5). Logremos esto, eligiendo las constantes C, y D,,_ Poniendo en la igualdad (16) t = 0, obtenemos (véase la condicién (a): (17) Si la funcién f (x) es tal que en el intervalo (0, 1) podemos desarro- Marla en la serie de Fourier (véase § 1, cap. XVII), entonces se cum- ple la condicién (17), si ponemos (18) Luego, derivamos los términos de la igualdad (16) respecto a t y po- nemos t = 0. De la condicién (5) se obtiene la igualdad Determinemos los coeficientes de Fourier de esta serie: t | ot sen “Fads o =| otesen rae (19) ann o Pues, hemos demostrado que si la serie (16), donde los coefi- cientes C, y D, estén determinados segiin las formulas (18) -y (19), permite la derivacién doble término a término, entonces, esta serie representa la‘funcién wu (z, ) que es la solucién de la ecuacién (1) y satisface a las condiciones de contorno e iniciales (2) — (5). Observacién. Al resolver el problema examinado para la ecua- cién de onda mediante otro procedimiento, se puede demostrar que la serie (16) representa la solucién, también, en el caso en que esta serie no permite la derivacién término a término. Aqui, la funcién f (x) debe ser dos veces derivable y @ (z), una vez derivable.Ecuacién de propagacién del calor en un véstago 378 § 4. ECUACION DE PROPAGACION DEL CALOR EN UN VASTAGO. PLANTEO DEL PROBLEMA CON VALORES DE CONTORNO Examinemos un véstago homogéneo de longitud 1. Supongamos que la superficio lateral del véstago no conduce el calor y que la temperatura es igual en todos los puntos de la seccién transversal del vastago. Estudiemos el proceso de difusién del calor en el vastago. [Link] eje Or de modo que un extremo del vastago coincida con el punto z = 0, y el otro, con el punto x = J (fig. 373). ( a 5% z Fig. 373 Sea u(z, t) la temperatura en la seccién de lvdstago con abscisa z en el momento t. Experimentalmente esta establecido que la velo- cidad de difusién del calor, es decir, la cantidad del calor que fluye a través de la seccién de abscisa z por la unidad de tiempo, se deter- mina por la férmula: g=—K&s, 0 donde S es el rea de la’ seccién del vastago examinado y k es el coeficiente de conductibilidad térmica *). ‘Analicemos un elemento del vdstago comprendido entre las secciones de abscisas x, y 22 (t, — 2; = Az). La cantidad del calor que pasa a través de la seccién de abscisa z,, durante el tiempo At es igual a: ou ee 2) AQ — hE] Sat (2) lo mismo tiene lugar para la seccién de abscisa 2: nou AQ FEE] Sat (3) *) La velocidad de propagacién del calor o velocidad del flujo térmico so de- termina del modo siguiente: read. q=lim 2, donde AQ es la cantidad de calor que pasa por la seccién S en el tiempo At.