Source: http://docplayer.es/328770-Un-nuevo-enfoque-para-el-tratamiento-de-los-terminos-difusivos-de-la-ecuacion-de-conveccion-difusion-en-el-metodo-de-galerkin-discontinuo.html
Timestamp: 2017-02-25 05:39:59+00:00

Document:
Un nuevo enfoque para el tratamiento de los términos difusivos de la ecuación de convección-difusión en el método de Galerkin discontinuo - PDF
Download "Un nuevo enfoque para el tratamiento de los términos difusivos de la ecuación de convección-difusión en el método de Galerkin discontinuo"
Ernesto Rivas Poblete
1 REVISTA INTERNACIONAL MÉTODOS NUMÉRICOS EN INGENIERÍA Rev. Int. Mét. Numér. Ing. Un nuevo enfoque para el tratamiento de los términos difusivos de la ecuación de convección-difusión en el método de Galerkin discontinuo H. Gómez, I. Colominas *, F. Navarrina, M. Casteleiro Grupo de Métodos Numéricos en Ingeniería, GMNI Departamento de Métodos Matemáticos y de Representación E.T.S. de Ingenieros de Caminos, Canales y Puertos Universidad de A Coruña Campus de Elviña, A Coruña, ESPAÑA RESUMEN Los modelos utilizados en la práctica para problemas de convección-difusión están habitualmente basados en la ley de Fick. En ciertas aplicaciones el uso de la ley de Fick lineal puede ser adecuado, a pesar de predecir velocidades de propagación infinitas. Sin embargo, en muchas otras aplicaciones es necesario el uso de modelos con velocidad finita de propagación [1]. A lo largo de la historia se han propuesto varios modelos que intentan solventar este inconveniente y, sin duda, el modelo lineal más aceptado entre la comunidad científica es el propuesto por vez primera por Cattaneo en 1948 [2, 3]. El modelo de Cattaneo fue desarrollado para problemas de difusión pura. Sin embargo, en aplicaciones prácticas de ingeniería el transporte se produce conjuntamente por convección y por difusión. Por ello los autores han propuesto recientemente una generalización del modelo de Cattaneo aplicable a problemas convectivo-difusivos que da lugar a un modelo de ecuaciones hiperbólicas [4]. Por otra parte, en el marco de los métodos de Galerkin discontinuo, la discretización de los términos difusivos es compleja y requiere la utilización de métodos híbridos. Sin embargo, el método de Galerkin discontinuo es más adecuado para la discretización del modelo hiperbólico de convección-difusión propuesto por los autores. Por ello, en este artículo se presenta una formulación de Galerkin discontinuo para el modelo hiperbólico de conveccióndifusión y se compara con el modelo parabólico discretizado mediante métodos de Galerkin discontinuo. SUMMARY The models that are being used to solve convection-diffusion problems in engineering are normally based on Fick s law. In some applications the use of linear Fick s law is accurate enough in spite of predicting an infinite speed of propagation. However, there are many other applications where a model that predicts a finite speed of propagation is mandatory [1]. A number of models have been proposed so far to eliminate the infinite speed paradox. The linear model that has been most commonly accepted by scientist and engineers is the one proposed initially by Cattaneo in 1948 [2, 3]. Cattaneo s law was originally proposed for pure-diffusive problems. However, in the context of most engineering applications it is important to consider both, diffusion and convection, phenomena since neither of these transport processes can be ignored. For this reason, the authors have recently proposed a generalization of Cattaneo s law that can be used in convective-diffusive problems. This constitutive equation, together with the mass conservation equation, leads to a hyperbolic model for convection-diffusion [4]. Correspondencia a: E.T.S. de Ingenieros de Caminos, Canales y Puertos, Universidad de A Coruña, Campus de Elviña, A Coruña, ESPAÑA. (2 2 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO On the other hand, in the framework of the discontinuous Galerkin methods, the discretization of the diffusive terms is by no means trivial and the use of hybrid methods is standard. However, the discontinuos Galerkin method is more appropriate for the discretization of the hyperbolic model. For this reason, we present herein a discontinuous Galerkin scheme for the hyperbolic convection-diffusion model. The numerical results are compared with those of the parabolic model solved by discontinuous Galerkin methods for parabolic problems. PALABRAS CLAVE: infinita Convección-Difusión, Ley de Cattaneo, Hiperbólico, Galerkin discontinuo, Velocidad 1. INTRODUCCIÓN El método de Galerkin discontinuo (Discontinuous Galerkin Method, DGM) es normalmente atribuido a Reed y Hill [5] y fue desarrollado para resolver la ecuación del transporte de neutrones. Aunque la relevancia del método fue advertida inmediatamente, la actual popularidad de los esquemas de Galerkin discontinuo se debe en gran medida a los trabajos de Cockburn y Shu. Desde finales de los años 80 y durante la década de los 90 publicaron (con la colaboración de otros autores) una serie de artículos relativos a lo que se conoce como método de Runge-utta/Galerkin discontinuo [6, 7, 8, 9, 10]. El gran desarrollo de estos métodos en los últimos años se debe a que poseen la flexibilidad de un método de elementos finitos continuo (es posible modelar geometrías complejas con mallas no estructuradas), una gran estabilidad, la capacidad de proporcionar soluciones con orden de convergencia arbitrario en mallas muy distorsionadas y la posibilidad de paralelizar los cálculos obteniendo eficiencias muy altas. Es cierto que, desde un punto de vista teórico, los métodos de volúmenes finitos son capaces de proporcionar soluciones de alto orden sin aumentar el número de grados de libertad del problema. Sin embargo, en la práctica, cuando se utilizan mallas no estructuradas estos métodos acaban convirtiéndose en métodos cuadráticos o en el mejor de los casos en métodos de orden entre 2 y 3. La tendencia actual en dinámica de fluidos computacional es la de utilizar métodos de orden muy alto (al menos en las zonas en las que la solución es suficientemente regular) para así beneficiarse de la convergencia espectral que se produce al aumentar el orden del polinomio interpolador [11]. Estas son algunas de las causas del gran desarrollo de los métodos de Galerkin discontinuo en los últimos años. La característica principal del método de Galerkin discontinuo es que la solución numérica puede ser discontinua en las caras de los elementos quedando éstos acoplados mediante los flujos numéricos característicos de los métodos de volúmenes finitos. Del mismo modo que las funciones interpolantes, las funciones de test son también discontinuas por lo que se puede definir una función de test que tome valor unitario dentro del elemento y nulo en el resto del dominio asegurando de este modo que la discretización obtenida es conservativa. En consecuencia, si el esquema numérico es convergente, convergerá hacia una solución débil del problema en virtud del teorema de Lax [12]. Al contrario que los métodos de elementos finitos continuos, el método de Galerkin discontinuo fue desarrollado en primer lugar para ecuaciones hiperbólicas. Esta técnica es muy adecuada para la discretización de problemas hiperbólicos, pero su extensión a problemas elípticos o parabólicos no es en absoluto trivial. Desde hace algunos años, se ha producido una intensa investigación en la extensión de los métodos de Galerkin discontinuo para ecuaciones parabólicas y elípticas con el fin de utilizarlos para la discretización de los términos viscosos de las ecuaciones de Navier-Stokes. Probablemente el primer intento satisfactorio de aplicar una técnica de Galerkin discontinuo a la discretización de3 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 3 los términos viscosos de las ecuaciones de Navier-Stokes fue llevado a cabo por Bassi y Rebay en 1997 [13]. A este estudio le siguió el artículo de Cockburn y Shu [14] en el que se introdujo el método de Galerkin discontinuo local (LDG) que es una generalización del método de Bassi-Rebay (BR). En general, en cualquier método de Galerkin discontinuo, para discretizar una ecuación de segundo orden es necesario descomponerla en un sistema equivalente de primer orden, de modo que la función incógnita y sus derivadas son tratadas como variables independientes. Los resultados obtenidos mediante estos métodos son buenos pero el coste computacional aumenta de forma considerable debido al aumento del número de incógnitas. Por ello, en los últimos años se han publicado una gran cantidad de artículos en los que se aborda el tratamiento de los términos difusivos [15, 16, 17, 18, 19]. Sin embargo, son necesarias todavía mejoras considerables en la aplicación de este método a problemas elípticos y parabólicos. El modelo hiperbólico de convección-difusión propuesto por los autores en [4] y basado en la ley de Cattaneo se establece en términos de un sistema acoplado de 2 ecuaciones diferenciales en las que la concentración y el flujo de contaminante son las incógnitas. Este sistema de ecuaciones es totalmente hiperbólico, por lo que su discretización por el método de Galerkin discontinuo es muy eficiente, al no ser necesario discretizar ningún término de segundo orden. Este hecho constituye una ventaja del modelo de convección-difusión hiperbólico frente al modelo parabólico. En este artículo se comparan ambos modelos analizando su idoneidad para ser discretizados en el espacio por el método de Galerkin discontinuo. El resto del artículo está organizado del siguiente modo: en la sección 2 se revisa la formulación parabólica del problema de convección-difusión (se incluye también un término reactivo para generalizar la exposición). En la sección 3 se presenta y revisa el modelo hiperbólico propuesto por los autores en [4]. Las secciones 4 y 5 están dedicadas a la discretización por el método de Galerkin discontinuo de los modelos de convección-difusión parabólico e hiperbólico respectivamente. Estas secciones se encuentran complementadas con la sección 6, en la que se analizan diferentes aspectos computacionales. En la sección 7 se presentan y discuten varios ejemplos numéricos utilizando ambas formulaciones. La sección 8 está dedicada al análisis de un ejemplo práctico en ingeniería. El ejemplo consiste en el estudio de la evolución de un vertido accidental en el puerto de A Coruña. En este caso se estudia, además, la influencia en la solución de un término reactivo lineal. 2. FORMULACIÓN CLÁSICA DEL PROBLEMA DE CONVECCIÓN-DIFUSIÓN-REACCIÓN Para obtener las ecuaciones que gobiernan el problema se asumirá que se pueden desacoplar las ecuaciones hidrodinámicas de las ecuaciones del transporte, y así, calcular primero el campo de velocidades que se utiliza luego para resolver las ecuaciones que proporcionan los niveles de contaminación. Por tanto, supondremos conocido el campo de velocidades, teniendo presente que debe verificar las ecuaciones de la hidrodinámica. Además, se asumirá que el fluido es incompresible y que no existen fuentes ni sumideros, por lo que el campo de velocidades debe ser solenoidal. Teniendo en cuenta las consideraciones anteriores y suponiendo válida la ley de Fick, las ecuaciones que gobiernan el problema de convección-difusión-reacción son: u t + a x (u) + x (q) = r(u) (1.1) q = D x (u) (1.2)4 4 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO donde (1.1) es la ecuación de conservación de masa de contaminante y (1.2) es la ecuación de Fick (ecuación constitutiva). En las ecuaciones anteriores u es la concentración de contaminante, r(u) es un término reactivo, a es el vector velocidad, q es el flujo difusivo de contaminante por unidad de densidad de fluido y D es el tensor de difusividad (se asume que es definido positivo). El sistema de ecuaciones (1) se puede desacoplar ya que es posible introducir (1.2) en (1.1) y resolver la ecuación escalar u t + a x (u) x (D x (u)) = r(u) (2) Por último, introduciremos la hipótesis de que el medio es homogéneo e isótropo (por tanto D = DI). Teniendo presente lo anterior, consideremos el transporte por convección-difusión-reacción en un dominio Ω R 2 con un contorno regular. Denotamos Γ a este contorno y n a su normal exterior unitaria. Supongamos que el contorno está compuesto por una parte Γ D no vacía en la que se prescribe el valor de u (condición tipo Dirichlet) y una parte complementaria Γ N en la que se prescribe el flujo (condición tipo Neumann). Además, se supone conocida la distribución inicial de la cantidad transportada u. Con estas consideraciones podemos plantear el problema de valores iniciales y de contorno como sigue: dado un campo de velocidades solenoidal a, dada la difusividad D y dadas las convenientes condiciones iniciales y de contorno, hallar u: Ω [0, T ] R tal que u t + a x (u) D x (u) = r(u) en Ω (0, T ) (3.1) u(x, 0) = u 0 (x) en Ω (3.2) u = u D en Γ D (0, T ] (3.3) D x (u) n = h en Γ N (0, T ] (3.4) El problema (3) es un problema bien planteado en el sentido de Hadamard, pero se puede probar (ver [20] ó [21] para una demostración) que predice velocidades de propagación infinitas Ecuaciones fundamentales 3. UN MODELO HIPERBÓLICO PARA EL PROBLEMA DE CONVECCIÓN-DIFUSIÓN-REACCIÓN El modelo más comúnmente aceptado para describir la difusión lineal evitando la paradoja del transporte de masa a velocidad infinita es el propuesto por Cattaneo. La ley de Cattaneo se puede escribir como: q + τ q t = D x (u) (4) donde τ se denomina tensor de relajación y sus componentes tienen dimensiones de tiempo. La ley de Cattaneo fue propuesta para problemas de difusión pura y no se puede utilizar en su forma original para problemas de convección-difusión. Sin embargo, en las aplicaciones en ingeniería el transporte se produce conjuntamente por convección y por difusión. Por esta razón, los autores han propuesto recientemente la siguiente generalización de la ley de Cattaneo (ver [22] o la referencia posterior [23]) ( ) q q + τ t + x (q) a = D x (u) (5)5 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 5 Por tanto, la formulación completa del problema del transporte por convección-difusión-reacción en un medio incompresible y utilizando un modelo de transporte con velocidad de propagación finita es la siguiente: u t + a x (u) + x (q) = r(u) (6.1) ( ) q q + τ t + x (q) a = D x (u) (6.2) donde (6.1) es la ecuación de conservación de masa de contaminante y (6.2) es la ecuación de Cattaneo con término convectivo. En estas ecuaciones se observa que la formulación propuesta es más general que la formulación basada en la ley de Fick ya que tomando τ = 0 en (6.2) se recuperan las ecuaciones clásicas. Como se demostrará a continuación, el sistema (6) es totalmente hiperbólico. Además, en el caso de difusión pura (a = 0) y bajo las hipótesis de homogeneidad e isotropía (τ = τi, D = DI) el sistema (6) se puede escribir como una ecuación hiperbólica con carácter dispersivo (relación de dispersión no lineal) donde la velocidad máxima de propagación de las perturbaciones viene dada por [20]: c = D/τ (7) 3.2. Estudio del modelo propuesto como un sistema de ecuaciones en forma conservativa Si consideramos que el dominio es bidimensional y usamos la notación q = (q 1, q 2 ) T, a = (a 1, a 2 ) T, el sistema (6) se puede escribir (bajo la hipótesis de que el medio es homogéneo e isótropo y, por tanto, τ = τi, D = DI) del siguiente modo: donde U = u τq 1 τq 2 U t + x (F ) = S (8) ; F = ua 1 + q 1 τq 1 a 1 + Du ua 2 + q 2 τq 1 a 2 ; S = τq 2 a 1 τq 2 a 2 + Du r(u) q 1 q 2 (9) Para establecer las propiedades básicas de la ecuación (8) es necesario escribirla en forma no conservativa. De este modo, si denominamos F i a la columna i-ésima de la matriz F, definimos las matrices jacobianas A 1, A 2 como A 1 = U (F 1 ) = a 1 1/τ 0 D a a 1 ; A 2 = U (F 2 ) = a 2 0 1/τ 0 a 2 0 D 0 a 2 e introducimos la hipermatriz A T = (A 1, A 2 ), la ecuación (8) se puede reescribir como (10) U t + (A x )U = S (11) Si ahora definimos κ = (κ 1, κ 2 ) T como un vector arbitrario con la única restricción de que tenga norma Euclídea unitaria, se puede demostrar que el sistema (11) es totalmente hiperbólico si la ecuación det(ωi A κ) = 0 (12)6 6 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO tiene tres soluciones reales diferentes ω para valores arbitrariamente prescritos de κ [24]. Es fácil probar que las soluciones de (12) son ω 1 = a κ (13a) ω 2 = a κ c (13b) ω 3 = a κ + c (13c) Y, por tanto, (8) es un sistema totalmente hiperbólico. Además, las soluciones de (12) son los autovalores de la matriz A κ que se conoce normalmente como matriz de proyección. Llegados a este punto es sencillo calcular los autovectores de A κ. Una vez hecho esto, podemos definir la matriz cuadrada C como la matriz cuyas columnas son los autovectores de la matriz de proyección, esto es, C = κ 2 Dκ 1 Dκ 1 0 c c κ 1 Dκ 2 Dκ 2 (14) En consecuencia, sea Λ la matriz diagonal, tal que los elementos de la diagonal son los autovalores ω 1, ω 2, ω 3, entonces, se verifica la siguiente relación A κ = CΛC 1 (15) Es importante notar que, aunque las matrices A 1 y A 2 son diagonalizables, el sistema (8) no lo es, ya que A 1 y A 2 son diagonalizables en diferentes bases. En consecuencia, los cuasi-invariantes de Riemann [24] no se pueden definir en 2D. En este punto, es útil definir el número adimensional H = a c que es el homólogo del número de Mach [25] en problemas de flujo compresible o del número de Froude [26] en flujos con superficie libre. La definición (16) permite clasificar el flujo como sigue: H < 1 Flujo subcrítico H > 1 Flujo supercrítico H = 1 Flujo crítico Con la definición anterior, en flujo supercrítico no es posible el transporte de contaminante hacia aguas arriba porque el modo de propagación difusivo avanza con una velocidad menor que el modo de propagación convectivo. La imposición de las condiciones de contorno al sistema (8) no es trivial y una elección inadecuada puede provocar que la solución no exista o que no sea única. Una situación muy frecuente es tener que imponer una condición de contorno en una pared sólida fija. La condición de contorno en este caso es q n = 0, siendo n la normal exterior unitaria en cada punto del contorno. Sin embargo, los dominios computacionales están normalmente compuestos también de contornos de entrada de flujo y contornos de salida de flujo. En lo que sigue utilizaremos la notación Γ in para el contorno de entrada de flujo (aquel que verifica a n < 0) y Γ out para el contorno de salida de flujo (aquel que verifica a n > 0). En estos tipos de contornos la imposición de las condiciones de contorno es más complicada. Es conocido que para problemas hiperbólicos escalares sólo se debe imponer una condición en el contorno de entrada de flujo (ver, por ejemplo [24]), y para sistemas de ecuaciones hiperbólicos sólo se deben (16)7 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 7 prescribir las ondas que entran por el contorno. En consecuencia, es necesario un análisis de Riemann en la dirección de la normal exterior. Procediendo de este modo, se encuentran tres velocidades de propagación diferentes: ω 1 = a n (17a) ω 2 = a n c (17b) ω 3 = a n + c (17c) Para obtener un problema bien planteado se deben imponer en un determinado contorno tantas componentes de U como velocidades negativas (ondas entrantes) haya en (17). Por tanto, en un contorno de entrada de flujo supercrítico se deben prescribir todas las componentes de U mientras que en un contorno de salida de flujo supercrítico no se debe imponer ninguna. Por contra, en un contorno de salida de flujo subcrítico, sólo se debe prescribir una componente de U, mientras que en un contorno de entrada de flujo subcrítico se deben imponer dos componentes de U. Independientemente de las condiciones de contorno se debe dar una condición inicial para (8). Por tanto, el problema bidimensional del transporte por convección-difusión-reacción formulado con la ley de Cattaneo, se puede plantear como sigue: dados D, τ > 0, dado un campo de velocidades solenoidal a y dadas las convenientes condiciones iniciales y de contorno, hallar U : Ω [0, T ] R 3 tal que U t + x (F ) = S en Ω (0, T ) (18a) U(x, 0) = U 0 (x) en Ω (18b) U in = U in D en Γ D (0, T ] (18c) F in n = F N en Γ N (0, T ] (18d) siendo F y S los vectores definidos en (9), U in las componentes entrantes de U y F in las componentes entrantes de la matriz de flujo. 4. EL MÉTODO DE GALERIN DISCONTINUO PARA EL MODELO PARABÓLICO DE CONVECCIÓN-DIFUSIÓN-REACCIÓN Como se ha dicho anteriormente, para aplicar el método de Galerkin discontinuo a la ecuación de convección-difusión-reacción es necesario obtener la forma débil a partir del sistema (1) (en general para aplicar el método de Galerkin discontinuo a una ecuación de segundo orden es necesario escribirla como un sistema mayor de primer orden). Sea un subconjunto de Ω con contorno regular. Si denotamos n a su normal unitaria exterior, la forma débil de (1) se escribe como: v u t dx = w qdx = ua x (v)dx + q x (v)dx + vr(u)dx vua n dx vq n dx (19.1) ud x (w)dx udw n dx (19.2)8 8 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO A continuación procedemos a utilizar el método de Galerkin discontinuo. Sea Ω h = {} una partición de Ω. Definimos los espacios rotos V(Ω h ) y Σ(Ω h ) como los espacios de funciones cuya restricción a cada elemento pertenece a los espacios de Sobolev H 1 () y ( H 1 () ) 2. De este modo V = {v L 2 (Ω) : v H 1 (), Ω h } (20a) Σ = {w ( L 2 (Ω) ) 2 : w ( H 1 () ) 2, Ωh } (20b) A continuación introducimos los subespacios de elementos finitos V h V y Σ h Σ V h = {v L 2 (Ω) : v P m (), Ω h } (21a) Σ h = {w ( L 2 (Ω) ) 2 : w (P m ()) 2, Ω h } (21b) En la expresión anterior P m () es el espacio de polinomios de grado menor o igual que m. Nótese que los espacios discretos V h y Σ h no implican continuidad C 0 de la solución y que, en general, m puede ser diferente en cada elemento. Teniendo presente lo anterior, el problema se reduce a obtener u h V h, q h Σ h tales que: v u h t dx = w q h dx = u h a x (v)dx + q h x (v)dx + n dx e e vh u,e c e u h D x (w)dx e e e vr(u h )dx vh q,e d n dx v V h, Ω h (22.1) h u,e d Dw n w Σ h, Ω h (22.2) donde las sumas recorren las caras e del elemento y los flujos numéricos h u,e c, h q,e d, hu,e d son aproximaciones de au e, q e y u e respectivamente. La elección del flujo numérico convectivo h u,e c es trivial, ya simplemente se toma una de las funciones upwind estándar del método de volúmenes finitos. Por contra, la elección de los flujos numéricos difusivos (h q,e d y h u,e d ) es crucial para la estabilidad y la precisión del esquema, así como para los dominios de dependencia del método que determinarán la eficiencia de una eventual paralelización. Para definir los flujos numéricos difusivos seguiremos a Arnold et al [27]. Sea de nuevo e una cara compartida por los elementos 1 y 2 y sean n 1 y n 2 las normales exteriores a 1 y 2 respectivamente. Sea v una función definida en 1 2 posiblemente discontinua sobre e. Definimos v i = (v i ) e, i = 1, 2. Para una función escalar v se define y para una función vectorial w {v} = 1 2 (v 1 + v 2 ), [[v]] = v 1 n 1 + v 2 n 2 (23) {w} = 1 2 (w 1 + w 2 ), [w]] = w 1 n 1 + w 2 n 2 (24) Con la notación anterior los flujos numéricos se escriben como: h u,e d = {u h } + γ e [[u h ]] (25) h q,e d = {q h } α e ([[u h ]) + β e [[q h ] (26)9 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 9 En las expresiones anteriores β e y γ e son funciones vectoriales definidas en e. A menudo son constantes y, de hecho, para varios métodos ambas son nulas. En el caso más sencillo el término α e ([[u h ]]) se calcula como α e ([[u h ]]) = η e [[u h ]] (27) siendo η e una constante. En este artículo se analizarán el método de Bassi-Rebay y el método de Galerkin discontinuo local. El primero de los métodos se obtiene tomando γ e = β e = 0 y η e = 0. Para obtener el método de Galerkin discontinuo local se eligen γ e = β e con β e = 1 2 (S 12n 1 + S 21 n 2 ) (28) donde S 12, S 21 {0, 1} y verifican la relación S 12 + S 21 = 1. Para la constante η e se toma un valor positivo. Los ejemplos que se presentan en este artículo se han resuelto tomando el valor η e = EL MÉTODO DE GALERIN DISCONTINUO PARA EL MODELO HIPERBÓLICO DE CONVECCIÓN-DIFUSIÓN-REACCIÓN Para obtener la forma débil de (8) procederemos del mismo modo que en la sección 4. Si es un subconjunto de Ω con contorno regular y normal exterior unitaria n, la forma débil de (8) resulta: w U t dx F x (w)dx + wf ndx = w Sdx (29) Consideremos de nuevo la partición Ω h y sea W el espacio roto W = {w ( L 2 (Ω) ) 3 : w ( H 1 () ) 3, Ωh } (30) Denominamos W h a su espacio discreto asociado W h = {w ( L 2 (Ω) ) 3 : w (P m ()) 3, Ω h } (31) El problema se reduce, por tanto, a obtener U h W h, tal que w U h t dx F h x (w)dx+ e e w H U,e C dx = w S h dx w W h, Ω h (32) un flujo numérico convectivo upwind. Para determinar el flujo numérico H U,e C siendo H U,e C supongamos que la cara e es compartida por los elementos 1 y 2 donde 1 =. Con esta notación el flujo numérico se puede escribir como H U,e C = Ae +U e 1 + A e U e 2 (33) donde las matrices A e + y A e vienen definidas por donde Λ e ± = a n± a n A e ±C e = C e Λ e ± (34) a n c± a n c La matriz C ha sido definida en (14) para una dirección arbitraria κ. a n+c± a n+c 2 (35)10 10 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO 6.1. Integración numérica 6. OTROS ASPECTOS COMPUTACIONALES La integración numérica se lleva a cabo en todos los casos mediante cuadraturas de Gauss- Legendre. En todos los ejemplos se utilizarán elementos Lagrangianos bicuadráticos, por lo que para la integración numérica se emplearán 3 puntos de integración en cada dirección. Para las integrales de línea se utilizan también 3 puntos para cada cara Integración temporal Tanto para el modelo parabólico como para el modelo hiperbólico, ensamblando todas las contribuciones elementales, se obtiene un sistema de ecuaciones diferenciales ordinarias que gobierna la evolución temporal de las ecuaciones discretas. Este sistema se puede escribir en ambos casos como M du dt = R(U) (36) donde M es la matriz de masas global, U es el vector global de grados de libertad y R(U) es el vector residuo. Como las funciones de forma sólo toman valores distintos de cero dentro del elemento correspondiente, la matriz de masas M es diagonal por bloques. Debido a esta estructura tan particular de la matriz de masas, podría calcularse la matriz M 1 en el primer paso de tiempo para posteriormente resolver los sistemas de ecuaciones lineales pertinentes mediante una multiplicación (este procedimiento es bastante habitual en la bibliografía; ver por ejemplo [28]). Sin embargo, en este trabajo se ha optado por realizar la factorización de Cholesky de la matriz en el primer paso de tiempo y luego realizar en cada iteración posterior las correspondientes substituciones hacia delante o hacia atrás. Dada la estructura de la matriz, la factorización de M se obtiene ensamblando la factorización de cada uno de los bloques. Por tanto, esta operación tiene un coste computacional muy bajo. El sistema semidiscreto (36) se puede integrar utilizando el método TVD-Runge-utta de tercer orden propuesto por Shu y Osher [29]. Dadas las variables en el paso previo n, U n, las variables U n+1 se calculan de forma explícita en tres pasos U (1) = U n + tl(u n ) U (2) = 3 4 U n U (1) tl(u (1) ) U n+1 = 1 3 U n U (2) tl(u (2) ) (37a) En las ecuaciones anteriores el operador L( ) representa la derivada temporal del sistema de ecuaciones diferenciales ordinarias (36), de modo que L(U) = M 1 R(U). El método anterior es linealmente estable para números de Courant menores o iguales que 1/(2p + 1) donde p es el orden de los polinomios interpolantes Condiciones de contorno Los métodos de Galerkin discontinuo permiten la imposición de las condiciones de contorno tipo Dirichlet en forma débil. En la bibliografía es muy habitual que las condiciones de contorno se impongan de forma débil tanto en los contornos de salida de flujo como en los contornos de entrada de flujo. Muchos autores consideran la imposición de las condiciones de contorno en forma débil más11 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 11 Figura 1. Convección oblicua a la malla. Planteamiento del problema. D = 10 6, a = (cos α, sin α). Para el modelo hiperbólico se toma τ = 1. eficiente en los contornos de salida de flujo [30], ya que de este modo se evitan las oscilaciones que se obtienen al intentar resolver numéricamente una capa límite. Sin embargo, la imposición en forma débil de las condiciones Dirichlet en los contornos de entrada de flujo no ofrece ventajas claras [31]. En este artículo todas las condiciones de contorno se imponen de forma fuerte ya que no se aprecian inestabilidades significativas en ninguno de los contornos. 7. EJEMPLOS NUMÉRICOS. COMPARACIÓN DE LOS MODELOS PARABÓLICO E HIPERBÓLICO A continuación se resolverán varios de los problemas clásicos de convección-difusión utilizando tanto el modelo hiperbólico como el modelo parabólico. En ninguno de los ejemplos de esta sección se considerarán términos reactivos Convección oblicua a la malla El planteamiento del problema se presenta en la figura 1. El problema se resuelve en el dominio cuadrado [0, 1] [0, 1]. Sobre este dominio se define una malla regular elementos Lagrangianos bicuadráticos, por lo que la longitud de las caras de los elementos es h = En consecuencia, el número de Peclet elemental es P e = a h/(2d) = Se resolverá el problema para α = π/4 (ver figura 1). Para el modelo hiperbólico se deben prescribir todas la incógnitas en el contorno de entrada de flujo, mientras que no se debe imponer ninguna condición en el contorno de salida de flujo. En el contorno de entrada de flujo se ha impuesto la condición q = 0 junto con la condición discontinua de concentraciones que se encuentra dibujada en la figura 1. Para el modelo parabólico es necesario imponer una condición de contorno en todos los puntos de la frontera del dominio. En el contorno de salida de flujo se ha impuesto la condición q n = 0, mientras que en el contorno de entrada de flujo se prescribe la misma condición discontinua de concentraciones que se ha impuesto para el modelo hiperbólico. Las soluciones estacionarias utilizando los métodos de Bassi-Rebay, Galerkin discontinuo12 12 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO Figura 2. Convección oblicua a la malla. Solución estacionaria. Bassi-Rebay (izquierda), LDG (centro) y modelo hiperbólico (derecha). D = 10 6, a = (cos α, sin α), α = π/4. Para el modelo hiperbólico se toma τ = 1 Figura 3. Convección de un sinusoide en una rotación pura. Planteamiento del problema. D = 10 6, a = ( x 2, x 1 ). Para el modelo hiperbólico se toma τ = 1. local y el modelo hiperbólico se encuentran en la figura 2. Como se puede observar, las soluciones numéricas son prácticamente iguales. Esto se debe a la condición de contorno que se ha impuesto para el modelo parabólico en el contorno de salida de flujo. Si se impusiese para el modelo parabólico una condición tipo Dirichlet en el contorno de salida de flujo la solución numérica del modelo hiperbólico sería claramente más precisa que las soluciones del modelo parabólico Convección de un sinusoide en una rotación pura En esta sección se resuelve un problema de flujo circular alrededor del centro del dominio cuadrado [ 1, 1] [ 1, 1]. Este es un test clásico de precisión de un método numérico (ver, por ejemplo, [32]). El planteamiento del problema se encuentra en la figura 3. El campo de velocidades viene dado por la expresión a = ( x 2, x 1 ). La difusividad es D = 10 6 y para el modelo hiperbólico se toma una tiempo de relajación τ = 1. Por tanto, el problema está claramente dominado por la convección. Utilizaremos una malla regular de elementos Lagrangianos bicuadráticos, por lo que el número13 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 13 Figura 4. Convección de un sinusoide en una rotación pura. Solución estacionaria. Bassi-Rebay (izquierda), LDG (centro) y modelo hiperbólico (derecha). D = 10 6, a = ( x 2, x 1). Para el modelo hiperbólico se toma τ = 1. Figura 5. Convección de un sinusoide en una rotación pura. Solución estacionaria en x 2 = 0. Bassi-Rebay (rojo), LDG (azul) y modelo hiperbólico (verde). D = 10 6, a = ( x 2, x 1 ). Para el modelo hiperbólico se toma τ = 1. de Peclet elemental es P e En cuanto a las condiciones de contorno, para el modelo parabólico se impone la condición u = 0 en todo el contorno y en el segmento x 2 = 0, 0 x 1 1 se impone u(x 1, 0) = sen(πx 1 ). Para el modelo hiperbólico el flujo es supercrítico en todos los contornos del problema. En los contornos de entrada de flujo se imponen condiciones homogéneas para todas las variables y en los contornos de salida de flujo no se impone ninguna variable. Además, en el segmento x 2 = 0, 0 x 1 1 se impone u(x 1, 0) = sen(πx 1 ) tal y como se puede ver en la figura 3. La solución exacta es esencialmente una convección pura de esta condición de contorno a lo largo de las líneas de corriente. En la figura 4 se representan las soluciones estacionarias utilizando el modelo parabólico y los métodos de Bassi- Rebay y de Galerkin discontinuo local y el modelo hiperbólico y el método de Galerkin discontinuo propuesto. En la figura 5 se dibujan las soluciones en x 2 = 0. Todas las soluciones son muy precisas y no existe difusión apreciable fuera de las líneas de corriente.14 14 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO Figura 6. Simulación de un vertido accidental en el puerto de A Coruña. Fotografía digital de la zona. 8. UN EJEMPLO DE APLICACIÓN PRÁCTICA: VERTIDO DE UN CONTAMINANTE EN EL PUERTO DE A CORUÑA En esta sección se simulará la evolución de un vertido accidental en el puerto de A Coruña [33]. El dominio del problema será la zona del espejo de agua del puerto, según se puede apreciar en las figuras 6 y 7. Se simulará la evolución de un vertido en la zona de los pantalanes petroleros. El objetivo de este ejemplo es demostrar que la metodología propuesta se puede aplicar a un problema real de ingeniería. Por ello, no se ha considerado necesario realizar una estimación de los parámetros del problema, lo que hubiese requerido la realización de ensayos experimentales y mediciones de campo. Del mismo modo, el campo de velocidades no se ha calculado sino que se ha supuesto un campo de velocidades solenoidal que sigue la dirección predominante de los vientos en la zona como se puede observar en el esquema de la figura 8 (izquierda). Para limitar el dominio del problema se ha definido un contorno en mar abierto mediante un arco que va desde el extremo del dique de abrigo hasta el extremo del muelle de la Dársena de Oza. La malla del problema se presenta en la figura 8. Como se puede apreciar en esta figura se han introducido algunas modificaciones respecto al dominio real representado en la figura 7. Estos cambios se han realizado con el fin de mejorar las características geométricas de la malla. Se han modificado los siguientes elementos: Los espigones situados entre la Dársena de la Marina y la Antedársena Los pantalanes del muelle de petroleros15 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 15 A Coruña A Coruña Figura 7. Simulación de un vertido accidental en el puerto de A Coruña. Esquema en planta del puerto. Realmente, la supresión de estos últimos elementos no constituye una simplificación del dominio ya que, al ser estas estructuras diáfanas, permiten el paso de agua, por lo que a efectos del flujo es como si no existiesen. También en la figura 8 se diferencian los tres tipos de contorno que existen en este problema: el contorno tipo pared se ha representado en color verde; el contorno donde se produce el vertido se ha representado en color rojo; y el contorno de mar abierto se ha dibujado en azul. En el contorno tipo pared se impone como condición de contorno que el flujo normal de contaminante sea nulo. En el contorno donde se produce el vertido se impone un flujo normal de valor q n = Por último, en el contorno de mar abierto se impone la condición de tipo mixto q n = cu, donde c es la celeridad de la onda difusiva. Asimismo, la figura 8 muestra el campo de velocidades utilizado para este problema. El flujo está caracterizado por un valor máximo del número H (ver la definición en (16)) que viene dado por H máx = , por lo que estamos en condiciones de flujo subcrítico en todo el dominio. El tiempo se mide mediante la variable adimensional t = t/τ Transporte por convección-difusión En primer lugar se ha resuelto el problema sin tener en cuenta ningún término reactivo. De este modo se pretende ver los efectos del vertido en el caso más desfavorable posible. En las figuras 9 10 se muestra la condición inicial de concentraciones y las soluciones de concentraciones para t = 25, t = 50, t = 75, t = 100, t = 150, t = 200 y t = 400.16 16 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO Figura 8. Simulación de un vertido accidental en el puerto de A Coruña. Campo de velocidades y tipos de contornos (izquierda) y malla de elementos finitos (derecha). En la parte izquierda de la figura el contorno tipo pared se ha coloreado en verde; el contorno donde se produce el vertido en rojo; y el contorno de mar abierto se encuentra dibujado en azul. La malla de cálculo consta de 2023 elementos Lagrangianos bicuadráticos y ha sido generada con el código GEN4U [34] Transporte por convección-difusión-reacción En esta sección se analiza la influencia en la solución de un término reactivo lineal del tipo r(u) = λu, λ > 0. La inclusión de este tipo de términos reactivos es muy importante desde el punto de vista práctico, ya que uno de los métodos paliativos que se utilizan en situaciones de vertidos marinos consiste en forzar reacciones químicas para que el contaminante se extienda lo menos posible. La intensidad de la reacción con respecto a la convección viene dada por el número de Damköhler, esto es, D a = λh (38) a donde h es una longitud característica de la malla [35]. El número de Damköhler mínimo para el problema que se ha resuelto es D a En la figura 11 se comparan las soluciones estacionarias con término reactivo y sin término reactivo. Se puede observar que los efectos del vertido son mucho menores cuando se considera el término reactivo. 9. CONCLUSIONES Y LÍNEAS FUTURAS DE DESARROLLO En este artículo se ha presentado un nuevo enfoque para el tratamiento de los términos difusivos de la ecuación de convección-difusión en el método de Galerkin discontinuo. El procedimiento está basado en la utilización de un modelo hiperbólico para el problema de convección-difusión que ha sido recientemente propuesto por los autores. Este modelo tiene, con respecto al modelo parabólico, la ventaja de que predice una velocidad finita de propagación para el contaminante. Además, el modelo propuesto puede ser discretizado de forma muy eficiente mediante el método de Galerkin discontinuo17 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 17 Figura 9. Simulación de un vertido accidental en el puerto de A Coruña. Se presentan (de izquierda a derecha y de arriba a abajo) la condición inicial de concentraciones y las soluciones de concentración para los tiempos adimensionales t = 25, t = 50 y t = 75. ya que se puede escribir como un sistema de ecuaciones en forma conservativa. Desde el punto de vista numérico, la metodología de cálculo presentada (modelo matemático hiperbólico junto con el método de Galerkin discontinuo propuesto) compite (tanto en coste computacional como en calidad de resultados) con el empleo del modelo parabólico junto con la utilización de los modelos de Galerkin discontinuo más avanzados, como son los métodos de Bassi- Rebay y Galerkin discontinuo local. AGRADECIMIENTOS18 18 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO Figura 10. Simulación de un vertido accidental en el puerto de A Coruña. Se presentan (de izquierda a derecha y de arriba a abajo) las soluciones de concentración para los tiempos adimensionales t = 100, t = 150, t = 200 y t = 400. Este trabajo ha sido parcialmente financiado por la Subdirección Xeral de I+D de la Xunta de Galicia (Proyectos PGIDIT06TAM11801PR, PGIDT03PXIC18001PN y PGIDT05PXIC18002PN), el Ministerio de Educación y Ciencia (Proyectos DPI y DPI ), la Universidad de A Coruña y la Fundación de la Ingeniería Civil de Galicia. REFERENCIAS 1. Compte A, The generalized Cattaneo equation for the description of anomalous transport processes, Journal of Physics A: Mathematical and General, 30 (1997) Cattaneo MC, Sulla conduzione de calore, Atti del Semin. Mat. e Fis. Univ. Modena, 3:83 101, 1948.19 TRATAMIENTO DE TÉRMINOS DIFUSIVOS EN EL MÉTODO DE GALERIN DISCONTINUO 19 Figura 11. Comparación de las soluciones estacionarias con y sin término reactivo 3. Cattaneo MC, Sur une forme de l equation de la chaleur éliminant le paradoxe d une propagation instantaneé, Comptes Rendus de L Academie des Sciences: Series I-Mathematics, 247: , Gómez H, Colominas I, Navarrina F, Casteleiro M, A finite element formulation for a convection-diffusion equation based on Cattaneo s law, Computer Methods in Applied Mechanics and Engineering, 196: , Reed WH, Hill TR, Triangular mesh methods for the neutron transport equation, Technical Report LA-UR , Los Alamos Scientific Laboratory, Cockburn B, arniadakis GE, Shu C-W, (Eds.), Discontinuous Galerkin methods. Theory, computation and applications, Springer, Cockburn B, Lin S-Y, Shu C-W, TVB Runge-utta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems, Journal of Computational Physics, 84:90 113, Cockburn B, Shu C-W, TVB Runge-utta local projection discontinuous Galerkin finite element method for scalar conservation laws II: general framework, Mathematics of Computation, 52: , Cockburn B, Shu C-W, The Runge-utta discontinuous Galerkin finite element method for conservation laws V: multidimensional systems, Journal of Computational Physics, 141: , Cockburn B, Shu C-W, Runge-utta discontinuous Galerkin methods for convection dominated flows, Journal of Scientific Computing, 16(3): , Barth T, An introduction to upwind finite volume and finite element methods: some unifying themes, VI Lecture Series 2005 on Higher Order Discretization Methods in CFD, Lax, PD, Hyperbolic systems of conservation laws and the mathematical theory of shock waves, SIAM, Philadelphia, Bassi F, Rebay S, A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations, Journal of Computational Physics, 131: , Cockburn B, Shu C-W, The local discontinuous Galerkin method for time-dependent convection-diffusion systems, SIAM Journal of Numerical Analysis, 35(6): , Baumann CE, Oden JT, A discontinuous hp finite element method for convection-diffusion problems, Computer Methods in Applied Mechanics and Engineering, 175: , Baumann CE, Oden JT, A discontinuous hp finite element method for the Euler and the Navier-Stokes equations, International Journal for Numerical Methods in Engineering, 31:79 95, Castillo P, Performance of discontinuous Galerkin methods for elliptic PDEs, SIAM Journal of Scientific Computing, 24(2): , Castillo P, Cockburn B, Perugia I, Schötzau D, An a priori error analysis of the local discontinuous Galerkin method for elliptic problems, SIAM Journal on Numerical Analysis, 38(5): , Sherwin SJ, irby RM, Peiró J, Taylor RL, Zienkiewicz OC, On 2D elliptic discontinuous Galerkin methods, International Journal for Numerical Methods in Engineering, 65: , Gómez H, Una formulación hiperbólica para el problema de convección-difusión en mecánica de fluidos computacional, Tesis Doctoral, Universidad de A Coruña, Antontsev S, Díaz JI, Shmarev S, Energy methods for free boundary problems. Applications to nonlinear PDEs and Fluid20 20 H. GÓMEZ, I. COLOMINAS, F. NAVARRINA, M. CASTELEIRO Mechanics, Series Progress in Nonlinear Differential Equations and Their Applications, No. 48, Birkäuser, Gómez H, Una nueva formulación para el problema del transporte por convección-difusión, Proyecto Técnico, Universidad de A Coruña, Christov CI, Jordan PM, Heat conduction paradox involving second-sound propagation in moving media, Physical Review Letters, 94(15): , Courant R, Hilbert D, Methods of mathematical physics. Vol II., John Wiley & Sons, Courant R, Friedrichs O, Supersonic flow and shock waves, Springer Verlag, Whitham GB, Linear and nonlinear waves, John Wiley & Sons, Arnold DN, Brezzi F, Cockburn B, Marini LD, Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems, SIAM Journal Numerical Analysis, 39: , Luo H, Baum JD, Löhner R, A p-multigrid discontinuous Galerkin method for the Euler equations on unstructured grids, Journal of Computational Physics, 211: , Shu C-W, Osher S, Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of Computational Physics, 77: , Bazilevs Y, Hughes TJR, Weak imposition of Dirichlet boundary conditions in fluid mechanics, Computers & Fluids, in press, Hughes TJR, Scovazzi G, Bochev PB, Buffa A, A multiscale discontinuous Galerkin method with the computational structure of a continuous Galerkin method, Computer Methods in Applied Mechanics and Engineering, 195: , Brooks A, Hughes TJR, Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations, Computer Methods in Applied Mechanics and Engineering, 32(1 3): , Figueroa CA, Colominas I, Mosqueira G, Navarrina F, Casteleiro M, A stabilized finite element approach for advectivediffusive transport problems, Proceedings of the XX Iberian Latin-American Congress on Computational Methods in Engineering (CDROM), eds. P.M. Pimenta, R.Brasil, E.Almeida, São Paulo, Brasil, Sarrate J, Huerta A, Efficient unstructured quadrilateral mesh generation, International Journal for Numerical Methods in Engineering, 49: , Hauke G, A simple subgrid scale stabilized method for the advection-diffusion-reaction equation, Computer Methods in Applied Mechanics and Engineering, 191: , 2002. Mostrar más
Congresso de Métodos Numéricos em Engenharia 2015 Lisboa, 29 de Junho a 2 de Julho 2015 c APMTAC, Portugal 2015 SUBESCALAS DE ORDEN REDUCIDO PARA MODELOS POD EN PROBLEMAS DE DINÁMICA DE FLUIDOS Joan Baiges,1,2, Más detalles La solución de algunas EDO de Riccati
Revista digital Matemática, Educación e Internet (http://tecdigital.tec.ac.cr/revistamatematica/). Vol 15, No 2. Marzo Agosto 2015. ISSN 1659-0643 La solución de algunas EDO de Riccati José Alfredo Jiménez Más detalles EL ANÁLISIS POR ELEMENTOS FINITOS: UNA METODOLOGÍA MUY RECIENTE EN ECONOMÍA
EL ANÁLISIS POR ELEMENTOS FINITOS: UNA METODOLOGÍA MUY RECIENTE EN ECONOMÍA Miguel Escribano Ródenas Departamento de Economía Financiera y Contabilidad I Escuela Universitaria de Estudios Empresariales Más detalles MÉTODO DE ELEMENTOS FINITOS EN MALLAS NO ESTRUCTURADAS APLICADO A FLUJO COMPRESIBLE
MÉTODO DE ELEMENTOS FINITOS EN MALLAS NO ESTRUCTURADAS APLICADO A FLUJO COMPRESIBLE Weht German, Giovacchini Juan Pablo, Sacco Carlos, D Errico Mario Dpto. Mec. Aer., Facultad de Ingeniería, Instituto Más detalles Tópicos de implementación en paralelo
Apéndice C Esquema PGP. Implementación en un entorno paralelo. Para poder abordar eficazmente la simulación computacional de grandes problemas estructurales, es necesario utilizar estrategias numéricas Más detalles Tema 2. Espacios Vectoriales. 2.1. Introducción
Tema 2 Espacios Vectoriales 2.1. Introducción Estamos habituados en diferentes cursos a trabajar con el concepto de vector. Concretamente sabemos que un vector es un segmento orientado caracterizado por Más detalles MODELO DE CONTAMINACIÓN DEL AIRE
ENFOQUTE. : 62-73 Copyright 200 Universidad Tecnológica Equinoccial ISSN: 390-6542 MODELO DE CONTMINCIÓN DEL IRE Iván Naula RESUMEN El presente documento estudia un modelo matemático de contaminación del Más detalles ESTUDIO DE FLUIDOS NO NEWTONIANOS CON LOS MÉTODOS DE VOLÚMENES Y ELEMENTOS FINITOS
REVISTA FACULTAD DE INGENIERÍA, U.T.A. (CHILE), VOL 10, 00, pp. 3-34 ESTUDIO DE FLUIDOS NO NEWTONIANOS CON LOS MÉTODOS DE VOLÚMENES Y ELEMENTOS FINITOS Guillermo Sánchez M. 1 César Vial R. Nelson Moraga Más detalles Una Aplicación del Cálculo Matricial a un Problema de Ingeniería
Divulgaciones Matemáticas Vol. 9 No. 221, pp. 197 25 Una Aplicación del Cálculo Matricial a un Problema de Ingeniería An Application of Matrix Calculus to an Engineering Problem P. R. Almeida Benítez Dpto. Más detalles I.- CONCEPTOS BÁSICOS DEL MÉTODO DE LOS ELEMENTOS FINITOS
I.- CONCEPTOS BÁSICOS DEL MÉTODO DE LOS ELEMENTOS FINITOS 1.1.- Introducción Muchos problemas de importancia práctica que frecuentemente aparecen en ingeniería, resultan de una complejidad matemática tal Más detalles Mecánica de Fluidos Computacional
Mecánica de Fluidos Computacional Patricio Bohórquez Escuela Politécnica Superior de Jaén Ingeniería Industrial Universidad de Jaén 22/Febrero/2011 Mecánica de Fluidos Computacional p.1 Curso de Mecánica Más detalles Agradecimientos AGRADECIMIENTOS
Agradecimientos AGRADECIMIENTOS Sirvan estas breves palabras como muestra de mi sentido y sincero agradecimiento hacia todas aquellas personas que han colaborado en la realización de esta tesina. Gracias Más detalles O.C. Zienkiewicz, R.L. Taylor. El Método de los Elementos Finitos. Vols 1 y 2. CIMNE-Mc Graw Hill, 1994.
Índice de la teoría 1. Presentación. Estas lecciones sólo pretenden ser una introducción que sirva para orientar sobre los conceptos, para un estudio más amplio se recomienda leer alguna publicación especializada, Más detalles Computational simulation of fluid flow by using element based finite volume method on unstructured grids
Rev. Téc. Ing. Univ. Zulia. Vol. 31, Nº 1, 13-25, 2008 Computational simulation of fluid flow by using element based finite volume method on unstructured grids Carlos D. Araujo, José A. Rincón, Gilberto Más detalles Espacios vectoriales con producto interno
Capítulo 8 Espacios vectoriales con producto interno En este capítulo, se generalizarán las nociones geométricas de distancia y perpendicularidad, conocidas en R y en R 3, a otros espacios vectoriales. Más detalles Flujo incompresible. Mario Storti. Centro Internacional de Métodos Numéricos en Ingeniería - CIMEC INTEC, (CONICET-UNL), Santa Fe, Argentina
Flujo incompresible Mario Storti Centro Internacional de Métodos Numéricos en Ingeniería - CIMEC INTEC, (CONICET-UNL), Santa Fe, Argentina mstorti@intec.unl.edu.ar http://www.cimec.org.ar/mstorti September Más detalles Resolución numérica de un modelo de frontera libre para el crecimiento tumoral
XX Congreso de Ecuaciones Diferenciales y Aplicaciones X Congreso de Matemática Aplicada Sevilla, 24-28 septiembre 2007 (pp. 1 8) Resolución numérica de un modelo de frontera libre para el crecimiento Más detalles Tema 1 Generalidades sobre Ecuaciones Diferenciales Ordinarias (E.D.O.)
Tema 1 Generalidades sobre Ecuaciones Diferenciales Ordinarias (E.D.O.) 1.1 Definiciones Se llama ecuación diferencial a toda ecuación que contiene las derivadas de una o más variables dependientes respecto Más detalles ASADES Avances en Energías Renovables y Medio Ambiente Vol. 8, Nº 2, 2004. Impreso en la Argentina. ISSN 0329-5184
ASADES Avances en Energías Renovables y Medio Ambiente Vol. 8, Nº 2, 2004. Impreso en la Argentina. ISSN 0329-5184 ANALISIS BIDIMENSIONAL DEL EFECTO DEL VIENTO SOBRE UN GRUPO DE VIVIENDAS N. Salvo, V. Más detalles INGENIERÍA INDUSTRIAL
PROGRAMA DE ASIGNATURA ACTIVIDAD CURRICULAR: ANÁLISIS NUMÉRICO Y CÁLCULO AVANZADO Código: 032535 Área: Matemática Aplicada Bloque: Ciencias Básicas Nivel: 3º Tipo: Obligatoria Modalidad: Anual Carga Horaria Más detalles Ejemplo 1.2 En el capitulo anterior se demostró que el conjunto. V = IR 2 = {(x, y) : x, y IR}
Subespacios Capítulo 1 Definición 1.1 Subespacio Sea H un subconjunto no vacio de un espacio vectorial V K. Si H es un espacio vectorial sobre K bajo las operaciones de suma y multiplicación por escalar Más detalles Computación Científica en Paralelo
Computación Científica en Paralelo Métodos de Descomposición de Dominio Luis Miguel de la Cruz luiggix@gmail.com www.dci.dgsca.unam.mx/lmcs Unidad de Investigación en Cómputo Aplicado DGSCA-UNAM. Posgrado Más detalles GENERACIÓN DE MALLAS DE HEXAEDROS EN GEOMETRÍAS MÚLTIPLEMENTE CONEXAS MEDIANTE SUBMAPPING
7th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 08) Vall de Núria, 9 a 11 de enero de 2008 c LaCàN, www.lacan-upc.es GENERACIÓN DE MALLAS DE HEXAEDROS EN GEOMETRÍAS MÚLTIPLEMENTE Más detalles Nombre de la asignatura: METODOS NUMERICOS. Carrera : Ingeniería Mecánica. Clave de la asignatura: ACB- 9311 Clave local:
Nombre de la asignatura: METODOS NUMERICOS Carrera : Ingeniería Mecánica Clave de la asignatura: ACB- 9 Clave local: Horas teoría horas practicas créditos: -0-8.- UBICACIÓN DE LA ASIGNATURA A) RELACIÓN Más detalles ECUACIONES DIFERENCIALES Y DIFERENCIAS FINITAS DIFFERENTIAL EQUATIONS AND FINITE DIFFERENCE
ECUACIONES DIFERENCIALES Y DIFERENCIAS FINITAS DIFFERENTIAL EQUATIONS AND FINITE DIFFERENCE Arturo Pérez París: Universidad de Alcalá. Madrid (España). arturo.perez@mixmail.com Julio Gutiérrez Muñoz: Universidad Más detalles PLAN DE ASIGNATURA. Presentación
PLAN DE ASIGNATURA Presentación Nombre de la asignatura: Técnicas de Modelación y Simulación Curso Académico: 2012-2013 Departamentos: o Ingeniería Mecánica. Área de Ingeniería Térmica y de Fluidos Área Más detalles 1 v 1 v 2. = u 1v 1 + u 2 v 2 +... u n v n. v n. y v = u u = u 2 1 + u2 2 + + u2 n.
Ortogonalidad Producto interior Longitud y ortogonalidad Definición Sean u y v vectores de R n Se define el producto escalar o producto interior) de u y v como u v = u T v = u, u,, u n ) Ejemplo Calcular Más detalles MÉTODO DE ELEMENTOS FINITOS APLICADO A FLUJO COMPRESIBLE CON GAS EN EQUILIBRIO
MÉTODO DE ELEMENTOS FINITOS APLICADO A FLUJO COMPRESIBLE CON GAS EN EQUILIBRIO Germán Weht, Juan P. Giovacchini, Carlos G. Sacco, Mario A. D Errico Dpto. Mec. Aer., Facultad de Ingeniería, Instituto Universitario Más detalles TEMA 4: CALCULO NUMERICO DE AUTOVALORES
Lino Alvarez - Aurea Martinez METODOS NUMERICOS TEMA 4: CALCULO NUMERICO DE AUTOVALORES 1 INTRODUCCION La determinación de autovalores y autovectores de una matriz cuadrada A de orden n es un problema Más detalles MÉTODOS DIRECTOS PARA LA RESOLUCIÓN DE ECUACIONES ALGEBRAICAS LINEALES
CAPÍTULO 4 EJERCICIOS RESUELTOS: MÉTODOS DIRECTOS PARA LA RESOLUCIÓN DE ECUACIONES ALGEBRAICAS LINEALES Ejercicios resueltos 1 1. Determine el número de operaciones aritméticas necesarias para calcular Más detalles Dominios ficticios y conjuntos de nivel para problemas de frontera libre. Aplicación al crecimiento tumoral
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, 21-25 septiembre 2009 (pp. 1 8) Dominios ficticios y conjuntos de nivel para problemas de frontera Más detalles Solución numérica de la ecuación de difusión-advección-reacción con un esquema de separación de operadores
Miscelánea Matemática 39 (4) 6 SMM Solución numérica de la ecuación de difusión-advección-reacción con un esquema de separación de operadores David Parra Guevara Centro de Ciencias de la Atmósfera Universidad Más detalles Apuntes sobre algunos teoremas fundamentales de análisis complejo, con 20 ejemplos resueltos (2007-08)
Variable Compleja I (3 o de Matemáticas) Apuntes sobre algunos teoremas fundamentales de análisis complejo, con ejemplos resueltos (7-8) En estos apuntes, consideraremos las funciones anaĺıticas (holomorfas) Más detalles Fascículo 2. Álgebra Lineal. Cursos de grado. Gabriela Jeronimo Juan Sabia Susana Tesauri. Universidad de Buenos Aires
Fascículo 2 Cursos de grado ISSN 1851-1317 Gabriela Jeronimo Juan Sabia Susana Tesauri Álgebra Lineal Departamento de Matemática Facultad de Ciencias Exactas y Naturales Universidad de Buenos Aires 2008 Más detalles De la invariancia de las ecuaciones proporcionadas por la mecánica analítica, se resolverá cualquier problema mecánico del mismo modo.
Programa de la asignatura Curso: 2006 / 2007 AMPLIACIÓN DE FÍSICA (3170) PROFESORADO Profesor/es: LUIS ANDRÉS VEGA GONZÁLEZ - correo-e: verga@ubu.es FICHA TÉCNICA Titulación: INGENIERÍA DE CAMINOS, CANALES Más detalles Clasificación de métricas.
Clasificación de métricas. 1. El problema de clasificación. Como bien sabemos, el par formado por una métrica T 2 (esto es, un tensor 2-covariante simétrico) sobre un espacio vectorial E, (E, T 2 ), constituye Más detalles Valores propios y vectores propios
Capítulo 6 Valores propios y vectores propios En este capítulo investigaremos qué propiedades son intrínsecas a una matriz, o su aplicación lineal asociada. Como veremos, el hecho de que existen muchas Más detalles ESCUELA TÉCNICA SUPERIOR DE INEGENIEROS AERONÁUTICOS UNIVERSIDAD POLITÉCNICA DE MADRID. Juan Manuel Tizón Pulido
SEMINARIO DE SIMULACIÓN NUMÉRICA EN SISTEMAS DE PROPULSIÓN Juan Manuel Tizón Pulido DEFINICIÓN y UTILIDAD Computational Fluid Dynamics (CFD) engloba las técnicas que se encargan de predecir el comportamiento Más detalles Uso de un software como propuesta metodológica para la enseñanza de ecuaciones diferenciales
Sociedad Colombiana de Matemáticas XV Congreso Nacional de Matemáticas 2005 Apuntes Lecturas Matemáticas Volumen Especial (2006), páginas 361 369 Uso de un software como propuesta metodológica para la Más detalles 1. Producto escalar, métrica y norma asociada
1. asociada Consideramos el espacio vectorial R n sobre el cuerpo R; escribimos los vectores o puntos de R n, indistintamente, como x = (x 1,..., x n ) = n x i e i i=1 donde e i son los vectores de la Más detalles 4. Ecuaciones integrales de la transferencia de calor.
Departamento de Ingeniería Química 76.47 Fenómenos De Transporte -76.03 Operaciones I PROGRAMA ANALÍTICO 1. Nociones fundamentales de la mecánica de fluidos. Concepto de medio continuo. El fluido como Más detalles CÁLCULO PARA LA INGENIERÍA 1
CÁLCULO PARA LA INGENIERÍA 1 PROBLEMAS RESUELTOS Tema 3 Derivación de funciones de varias variables 3.1 Derivadas y diferenciales de funciones de varias variables! 1. Derivadas parciales de primer orden.! Más detalles Tema 3 Resolución de Sistemas de Ecuaciones Lineales
Tema Resolución de Sistemas de Ecuaciones Lineales Índice Introducción 2 Método de Gauss 2 Resolución de sistemas triangulares 22 Triangulación por el método de Gauss 2 Variante Gauss-Jordan 24 Comentarios Más detalles Aplicaciones lineales
Capítulo 4 Aplicaciones lineales 4.1. Introduccción a las aplicaciones lineales En el capítulo anterior encontramos la aplicación de coordenadas x [x] B que asignaba, dada una base del espacio vectorial, Más detalles Descomposición de dominios
Descomposición de dominios Miguel Vargas 27/10/10 1/29 Contenido Contenido Solución de ecuaciones diferenciales con descomposición de dominios Dominios sin traslape, complemento de Schur Método alternante Más detalles CURSO DE MÉTODOS DE LA FÍSICA MATEMÁTICA ANÁLISIS FUNCIONAL OPERADORES NO ACOTADOS
CURSO DE MÉTODOS DE LA FÍSICA MATEMÁTICA ANÁLISIS FUNCIONAL H. FALOMIR DEPARTAMENTO DE FÍSICA FACULTAD DE CIENCIAS EXACTAS - UNLP OPERADORES NO ACOTADOS 1. Extensiones de operadores lineales Sea A un operador Más detalles 1. Examen 21/Junio/1994. Para la inversión de una matriz cuadrada A de orden n n, cuya inversa existe, se ha definido la siguiente iteración
CAPÍTULO 5 EJERCICIOS RESUELTOS: MÉTODOS ITERATIVOS PARA ECUACIONES LINEALES Ejercicios resueltos 1 1. Examen 21/Junio/1994. Para la inversión de una matriz cuadrada A de orden n n cuya inversa existe Más detalles Vectores y Valores Propios
Capítulo 11 Vectores y Valores Propios Las ideas de vector y valor propio constituyen conceptos centrales del álgebra lineal y resultan una valiosa herramienta en la solución de numerosos problemas de Más detalles Homogeneización de problemas (elasto)hidrodinámicos en lubricación
XX Congreso de Ecuaciones Diferenciales y Aplicaciones X Congreso de Matemática Aplicada Sevilla, 4-8 septiembre 7 pp. 8 Homogeneización de problemas elastohidrodinámicos en lubricación G. Bayada, S. Martin, Más detalles DIPLOMADO DINÁMICA DE FLUIDOS COMPUTACIONAL CFD
DIPLOMADO DINÁMICA DE FLUIDOS COMPUTACIONAL CFD 2015 INTRODUCCIÓN El Diplomado en Dinámica de Fluidos Computacional (CFD, por sus siglas en inglés Computational Fluid Dynamics), constituye uno de los pilares Más detalles Análisis comparativo del flujo y modelado de la transferencia de calor para una mezcla de gas y partículas en secadores ciclónicos por aspersión
CURSO CERO Departamento de Matemáticas Profesor: Raúl Martín Martín Sesiones 18 y 19 de Septiembre Capítulo 1 La demostración matemática Demostración por inducción El razonamiento por inducción es una Más detalles Nota 1. Los determinantes de orden superior a 3 se calculan aplicando las siguientes propiedades:
Capítulo 1 DETERMINANTES Definición 1 (Matriz traspuesta) Llamaremos matriz traspuesta de A = (a ij ) a la matriz A t = (a ji ); es decir la matriz que consiste en poner las filas de A como columnas Definición Más detalles MATEMÁTICAS EMPRESARIALES II:
MATEMÁTICAS EMPRESARIALES II: FUNCIÓN REAL DE VARIAS VARIABLES ÓPTIMOS DE UNA FUNCIÓN ESCALAR MATERIAL DIDÁCTICO DE SOPORTE González-Vila Puchades, Laura Ortí Celma, Francesc J. Sáez Madrid, José B. Departament Más detalles Repaso de matrices, determinantes y sistemas de ecuaciones lineales
Tema 1 Repaso de matrices, determinantes y sistemas de ecuaciones lineales Comenzamos este primer tema con un problema de motivación. Problema: El aire puro está compuesto esencialmente por un 78 por ciento Más detalles FLUJOS TURBULENTOS BIFASICOS 3D EN LLENADO DE MOLDES DE FUNDICION Y EN IMPACTO DE OLAS EN ESTRUCTURAS Nelson Moraga, Carlos Garrido y Daniel Garrido
FLUJOS TURBULENTOS BIFASICOS 3D EN LLENADO DE MOLDES DE FUNDICION Y EN IMPACTO DE OLAS EN ESTRUCTURAS Nelson Moraga, Carlos Garrido y Daniel Garrido Departamento de Ingeniería Mecánica Universidad de La Más detalles Computational simulation of fluid flow by using element based finite volume method on unstructured grids
Rev. Téc. Ing. Univ. Zulia. Vol. 3, Nº, 3-2, 28 Computational simulation of fluid flow by using element based finite volume method on unstructured grids Carlos D. Araujo, José A. Rincón, Gilberto I. Materano Más detalles Grado en Ingeniería Energética
CURSO 2014/2015 Grado en Ingeniería Energética DATOS DE LA ASIGNATURA Nombre: Matemáticas II Denominación en inglés: Applied Mathematics Código: Carácter: 606711105 Básico Créditos: Totales: Teóricos: Más detalles Capitán de fragata ingeniero AGUSTÍN E. GONZÁLEZ MORALES ÁLGEBRA PARA INGENIEROS
Capítulo 4 Ortogonalidad y Series de Fourier El adjetivo ortogonal proviene del griego orthos (recto) y gonia (ángulo). Este denota entonces la perpendicularidad entre dos elementos: dos calles que se Más detalles Panorama del curso Métodos Numéricos I
Espacios Vectoriales Departamento de Matemáticas, CCIR/ITESM 4 de enero de 2 Índice 3.. Objetivos................................................ 3.2. Motivación............................................... Más detalles Espacios generados, dependencia lineal y bases
Espacios generados dependencia lineal y bases Departamento de Matemáticas CCIR/ITESM 14 de enero de 2011 Índice 14.1. Introducción............................................... 1 14.2. Espacio Generado............................................ Más detalles Tema 3. Espacios vectoriales
Tema 3. Espacios vectoriales Estructura del tema. Definición y propiedades. Ejemplos. Dependencia e independencia lineal. Conceptos de base y dimensión. Coordenadas Subespacios vectoriales. 0.1. Definición Más detalles APLICACIÓN DE LOS FILTROS DE KALMAN A SISTEMAS DE CONTROL. José Gregorio Díaz (1); Ana María Mejías (2); Francisco Arteaga (3)
APLICACIÓN DE LOS FILTROS DE KALMAN A SISTEMAS DE CONTROL José Gregorio Díaz (1); Ana María Mejías (2); Francisco Arteaga (3) (1) Universidad de Carabobo, Escuela de Ing. Eléctrica. Dpto. de Sistemas y Más detalles 6 CONCLUSIONES Y LÍNEAS FUTURAS DE TRABAJO
6 CONCLUSIONES Y LÍNEAS FUTURAS DE TRABAJO Sirva este capítulo como resumen de las conclusiones que han ido apareciendo durante el desarrollo de este trabajo. Si bien a lo largo de los distintos capítulos Más detalles ALGEBRA LINEAL. Héctor Jairo Martínez R. Ana María Sanabria R.
ALGEBRA LINEAL Héctor Jairo Martínez R. Ana María Sanabria R. SEGUNDO SEMESTRE 8 Índice general. SISTEMAS DE ECUACIONES LINEALES.. Introducción................................................ Conceptos Más detalles Matrices equivalentes. El método de Gauss
Matrices equivalentes. El método de Gauss Dada una matriz A cualquiera decimos que B es equivalente a A si podemos transformar A en B mediante una combinación de las siguientes operaciones: Multiplicar Más detalles Departamento Aeronáutica
Departamento Aeronáutica 3 Encuentro Nacional de Computación de Alto Rendimiento para Aplicaciones Científicas - WHPC 2014 Córdoba, 5 y 6 de agosto Grupo de Trabajo Grupo de Trabajo Este grupo de trabajo Más detalles PROBLEMA 1. 1. [1.5 puntos] Obtener la ecuación de la recta tangente en el punto ( 2, 1) a la curva dada implícitamente por y 3 +3y 2 = x 4 3x 2.
PROBLEMA. ESCUELA UNIVERSITARIA POLITÉCNICA DE SEVILLA Ingeniería Técnica en Diseño Industrial Fundamentos Matemáticos de la Ingeniería Soluciones correspondientes a los problemas del Primer Parcial 7/8. Más detalles suarez@us.es, madelgado@us.es, cristianm@us.es
Estudio teórico de algunos modelos de terapias anti-angiogénicas Antonio Suárez, Manuel Delgado, Cristian Morales-Rodrigo, Dpto. de Ecuaciones Diferenciales y Análisis Numérico, Universidad de Sevilla, Más detalles Universidad de Costa Rica Escuela de Matemática ALGEBRA LINEAL. x x1 n. θ y. 1 n x1 n ȳ1 n. Carlos Arce S. William Castillo E. Jorge González V.
Universidad de Costa Rica Escuela de Matemática ALGEBRA LINEAL x x x1 n θ y y ȳ1 n 1 n x1 n ȳ1 n Carlos Arce S. William Castillo E. Jorge González V. 2003 Algebra Lineal Carlos Arce S., William Castillo Más detalles ERNESTO ARANDA. 1. Introducción 1 1.1. Ejemplo... 2
EL MÉTODO DE LOS ELEMENTOS FINITOS DESDE UNA PERSPECTIVA MATEMÁTICA ERNESTO ARANDA Índice. Introducción.. Ejemplo....................................... 2 2. Formulación variacional de problemas de contorno Más detalles CONDICIONES DE CONTORNO EN EL METODO DE DIFERENCIAS FINITAS CON MALLAS IRREGULARES
Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería. Vol. 4, 1, 31-40( 1988) CONDICIONES DE CONTORNO EN EL METODO DE DIFERENCIAS FINITAS CON MALLAS IRREGULARES BIBIANA M. LUCCIONI* Más detalles Estrategia de Refinamientos Automáticos hp Integrada con un Resolvedor de Dos Mallas: Aplicaciones al Electromagnetismo. David Pardo Zubiaur
Estrategia de Refinamientos Automáticos hp Integrada con un Resolvedor de Dos Mallas: Aplicaciones al Electromagnetismo dzubiaur@yahoo.es Supervisores: L. Demkowicz, C. Torres-Verdin, L. Tabarovski Otros Más detalles Departamento de Matemática Aplicada FUNDAMENTOS DE MATEMATICAS. Ingeniería Química (Curso 2005-06) Álgebra Lineal Práctica 3
1. Matrices en Matlab Departamento de Matemática Aplicada FUNDAMENTOS DE MATEMATICAS. Ingeniería Química (Curso 2005-06) Álgebra Lineal Práctica 3 Para introducir una matriz en Matlab se procede de la Más detalles Capitán de fragata ingeniero AGUSTÍN E. GONZÁLEZ MORALES. ÁLGEBRA PARA INGENIEROS (Solucionario)
Capitán de fragata ingeniero AGUSTÍN E. GONZÁLEZ MORALES ÁLGEBRA PARA INGENIEROS (Solucionario) 2 Í N D I C E CAPÍTULO : MATRICES, DETERMINANTES Y SISTEMAS DE ECUACIONES LINEALES CAPÍTULO 2: ESPACIOS VECTORIALES Más detalles INGENIERÍA MECÁNICA E INDUSTRIAL MATERIALES Y MANUFACTURA INGENIERÍA MECÁNICA División Departamento Licenciatura
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO FACULTAD DE INGENIERÍA PROGRAMA DE ESTUDIO MANUFACTURA III 10 10 Asignatura Clave Semestre Créditos INGENIERÍA MECÁNICA E INDUSTRIAL MATERIALES Y MANUFACTURA INGENIERÍA Más detalles Universidad Autónoma de Yucatán Facultad de Matemáticas
Universidad Autónoma de Yucatán Facultad de Matemáticas Ecuaciones Diferenciales Licenciatura en Ciencias de la Computación Cuarto semestre DESCRIPCIÓN: Durante este curso se presentarán técnicas analíticas, Más detalles Aplicaciones Lineales
Aplicaciones Lineales Concepto de aplicación lineal T : V W Definición: Si V y W son espacios vectoriales con los mismos escalares (por ejemplo, ambos espacios vectoriales reales o ambos espacios vectoriales Más detalles Espacio afín. Transformaciones afines y movimientos
Capítulo Espacio afín. Transformaciones afines y movimientos. Espacio afín y espacio afín métrico Definición. El espacio afín (tridimensional) está constituido por los siguientes elementos. El espacio Más detalles Universidad Politécnica de Madrid
Universidad Politécnica de Madrid Escuela Técnica Superior de Ingenieros de Minas Departamento de Matemática Aplicada y Métodos Informáticos x r 1 x 1 x 2 = x * d 1 * d 1 f(x)=cte. r Resolución de sistemas Más detalles Funciones lineales, diádicas y tensores en la mecánica newtoniana
Funciones lineales, diádicas y tensores en la mecánica newtoniana En este documento se introducen los conceptos de diádicas y tensores en un espacio vectorial euclídeo tridimensional, con objeto de dotar Más detalles SEMINARIO DE SIMULACIÓN PROPULSIÓN
SEMINARIO DE SIMULACIÓN NUMÉRICA EN SISTEMAS DE PROPULSIÓN Juan Manuel Tizón Pulido jmtizon@aero.upm.es http://webserver.dmt.upm.es/zope/dmt/members/jmtizon/libre-eleccion-1 Motores Cohete: Capítulo 2 Más detalles Introducción a la difusión en el tratamiento de imágenes
2 Introducción a la difusión en el tratamiento de imágenes Guillermo Gallego La idea que motiva el estudio de la difusión anisótropa en el tratamiento de imágenes es la búsqueda de métodos de suavizamiento Más detalles Cinemática Inversa del Robot. CI-2657 Robótica M.Sc. Kryscia Ramírez Benavides
Cinemática Inversa del Robot M.Sc. Kryscia Ramírez Benavides Introducción Resuelve la configuración que debe adoptar el robot para una posición y orientación del extremo conocidas. 2 Introducción (cont.) Más detalles GEOMETRÍA DEL ESPACIO EUCLÍDEO
CAPÍTULO I. GEOMETRÍA DEL ESPACIO EUCLÍDEO SECCIONES 1. Vectores. Operaciones con vectores. 2. Rectas y planos en R 3. 3. Curvas y superficies en R 3. 4. Nociones de topología métrica. 1 1. VECTORES. OPERACIONES Más detalles NUMERO DE CONDICION Y DETERMINANTE DE UNA MATRIZ Jorge Lemagne Pérez, Facultad de Matemática y Computación, Universidad de La Habana
REVISTA INVESTIGACION OPERACIONAL Vol. 2, No., 2000 NUMERO DE CONDICION Y DETERMINANTE DE UNA MATRIZ Jorge Lemagne Pérez, Facultad de Matemática y Computación, Universidad de La Habana RESUMEN En la resolución Más detalles Una nota sobre las estimas a priori del error para formulaciones mixtas aumentadas
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, 1-5 septiembre 009 (pp. 1 8) Una nota sobre las estimas a priori del error para formulaciones mixtas Más detalles UNIVERSIDAD LIBRE PROGRAMAS ANALÍTICOS PROGRAMA ACADÉMICO: SISTEMAS JORNADA: UNICA
PROGRAMAS ANALÍTICOS SECCIONAL: PEREIRA FACULTAD: INGENIERIAS PROGRAMA ACADÉMICO: SISTEMAS JORNADA: UNICA 1.- INFORMACION DEL DOCENTE: E-mail institucional: ingenierias@unilibrepereira.edu.co E-mail personal: Más detalles Métodos de Elementos Finitos y Aplicaciones
Métodos de Elementos Finitos y Aplicaciones Teoremas de Extensión para espacios de Sobolev Formulación del problema Dado un dominio abierto Ω R N y el espacio de Sobolev W 1,p (Ω) correspondiente, se busca Más detalles Ecuación de Bellman y Aplicaciones a la Optimización de Trayectorias
Resumen Implementación de Esquemas Numéricos para la Ecuación de Bellman y Aplicaciones a la Optimización de Trayectorias Departemento de Ingeniería Matemática Universidad de Chile 7 de enero de 2005 Resumen Más detalles Tema III. Capítulo 2. Sistemas generadores. Sistemas libres. Bases.
Tema III Capítulo 2 Sistemas generadores Sistemas libres Bases Álgebra Lineal I Departamento de Métodos Matemáticos y de Representación UDC 2 Sistemas generadores Sistemas libres Bases 1 Combinación lineal Más detalles Máster Universitario en Ingeniería de Caminos, Canales y Puertos Introducción al Análisis Numérico
Máster Universitario en Ingeniería de Caminos, Canales y Puertos Introducción al Análisis Numérico Departamento de Matemática Aplicada Universidad Granada Introducción El Cálculo o Análisis Numérico es Más detalles Matias Avila Salinas Fecha del documento: 22/11/2013 v 1.3.0 d3a5aa318ef63e709c908d2b3a54231e
Matias Avila Salinas Fecha del documento: 22/11/2013 v 1.3.0 d3a5aa318ef63e709c908d2b3a54231e Este fichero electrónico (PDF) contiene incrustada la tecnología CVN (CVN-XML). La tecnología CVN de este fichero Más detalles Factorizaciones de Cholesky, matrices definidas. semidefinidas positivas.
Factorizaciones de Cholesky, matrices definidas y semidefinidas positivas Héctor Manuel Mora Escobar Universidad Central, Bogotá hectormora@yahoo.com Junio de 2011 1 Introducción Este documento presenta, Más detalles Análisis III. Joaquín M. Ortega Aramburu
Análisis III Joaquín M. Ortega Aramburu Septiembre de 1999 Actualizado en julio de 2001 2 Índice General 1 Continuidad en el espacio euclídeo 5 1.1 El espacio euclídeo R n............................... Más detalles Las matrices tienen un número cada vez mas creciente de aplicaciones en la solución de problemas en Ciencia y Tecnología.
Aplicaciones de las Matrices a la Solución de Problemas de Redes Eléctricas Resumen Se muestra como obtener, sistemas de ecuaciones lineales que permitan calcular intensidades de corrientes en los ramales Más detalles 2017 © DocPlayer.es Política de privacidad | Condiciones del servicio | Feedback

References: Resolución 
 Resolución 
 RESOLUCIÓN 
 RESOLUCIÓN 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 resolución