Tesis para optar al grado de Magíster en Ciencias de la Ingeniería, Mención Matemáticas Aplicadas / Memoria para optar al título de Ingeniera Civil Matemática / Consideramos un volumen $V\subset \R^3$ generado al rotar alrededor del eje $Z$ un dominio $\Omega \subset \R^2$ acotado y suave que vive en el plano $XZ$. En este trabajo se construye una solución del flujo de mapa armónico del volumen $V$ a la esfera $S^2$ que revienta en tiempo finito, el problema es \begin{eqnarray*} v_t &=& \Delta v + |\nabla v |^2 v \text{ in } V \times (0,T)\\ v &=& v_{\partial V} \text{ in } \partial V \times (0,T)\\ v(\cdot , 0) &=&v_0 \text{ in } V, \end{eqnarray*} donde $v: V \times [0,T) \to S^2$, $v_0 : \overline{V} \to S^2$ es suave y $v_{\partial V}=\left. v_0\right|_{\partial V} : \partial V \to S^2$. Dado un punto $q \in \Omega$ de define la circunferencia $c(q)$ generada al rotar el punto $q$ alrededor del eje Z. Se encuentran datos iniciales y de frontera tales que la solución $v$ revienta exactamente en la curva $c(q)$ en un tiempo finito pequeño. La construcción de la solución se hace reduciendo el problema a 2 dimensiones y usando el método de Dávila, Del Pino y Wei \cite{dav} que transforma el problema en un sistema de \textit{inner-outer gluing} que separa el efecto principal de la ecuación cerca y lejos de la singularidad. Se obtiene una solución cuyo orden principal cerca de la singularidad tiene el perfil de un mapa armónico 1-corrotacional escalado. En la introducción se recuerdan la ecuación de flujo de mapa armónico y su origen, se establece el problema y la reducción a 2 dimensiones. En el primer capítulo se enuncian resultados útiles de topología y análisis funcional, y propiedades probadas en \cite{dav} para los mapas armónicos 1-corrotacionales y el operador linealizado en torno a ellos. En el segundo capítulo se obtiene un ansatz de la solución y se usa el método de Dávila, Del Pino y Wei \cite{dav} para reducir el problema a resolver un sistema de \textit{inner-outer gluing} que después se resuelve usando punto fijo. En el capítulo cuatro se obtienen las hipótesis para el punto fijo mediante estimaciones a priori obtenidas dividiendo el sistema en tres problemas principales: el problema interior, el problema exterior y el problema de los parámetros. En la parte final se concluye con algunas observaciones sobre este trabajo y posibles trabajos futuros en torno a el. / Este trabajo ha sido parcialmente financiado por el proyecto Fondecyt 1150066 y el Centro de Modelamiento Matemático, Proyecto Basal PFB 03 / Fondecyt 1150066 y CMM - Conicyt PIA AFB170001

Tesis para optar al grado de Magíster en Ciencias de la Ingeniería, Mención Matemáticas Aplicadas / Memoria para optar al título de Ingeniero Civil Matemático / En este trabajo se presenta y estudia un modelo para el campo de temperaturas en la cara superior de una pila de lixiviación en presencia de evaporación e irrigada de forma puntual.El modelo consiste en un sistema de ecuaciones en derivadas parciales deducido del balance de los flujos de calor y masa en la cara superior de la pila. En estas ecuaciones las incógnitas son el flujo de masa en la superficie de la cara superior de la pila y la temperatura de la pila.Se muestran distintos resultados acerca de la existencia de soluciones para estas ecuaciones evaluando tres casos: Modelo radial: La cara superior es un circulo. Modelo geometría general sin evaporación: La cara superior es una superficie plana con borde suave y evaporación nula. Modelo geometría general con evaporación: La cara superior es una superficie plana con borde suave, tomando en cuenta evaporación no nula. En cada uno de los modelos se estudian las condiciones bajo las cuales los sistemas de ecuaciones en derivadas parciales tienen solución y además que propiedades cumplen estas soluciones. En el modelo radial se demuestra la existencia y unicidad del campo bajo ciertas condi-ciones. También obtiene la existencia y unicidad de soluciones en el modelo con geometría general sin evaporación asumiendo que el riego no es puntual, en caso contrario solo se puede asegurar la existencia. En el modelo con geometría general y evaporación no nula, solo se logra demostrar la existencia de soluciones en un sentido débil. En todos los modelos se obtienen cotas que permiten entender el comportamiento del campo de temperaturas en función de las variables del problema. / CMM - Conicyt PIA AFB170001

The present work focuses on Nodal Discontinuous Galerkin Method applied to the one-dimensional linear advection equation, which approximates the global solution, partitioning its domain into elements. In each element the local solution is approximated by using interpolation in such a way that the total numerical solution is a direct sum of those approximations (polynomials). This method aims at reaching a high order through a simple implementation. This model is studied by Hesthaven and Warburton [16], with the particularity of Joining the best of the Finite Volumes Method and the best of Finit Element Method . First, the main results are revised in detail concerning the Jacobi orthogonal polynomials; more precisely, its generation formula and other results which help implementing the method. Concepts regarding interpolation and best approximation are studied. Furthermore, some notions about Sobolev space interpolation is revised. Secondly, theoretical aspects of the method are explained in detail , as well as its functioning. Thirdly, both the two method consistency theorems (better approximation and interpolation), proposed by Canuto and Quarteroni [4], and error behavior theorem based on Hesthaven and Warburton [16] are explained in detail. Finally, the consistency theorem referred to the interpolation is veri ed numerically through the usage of the Python language as well as the error behavior. It is worth mentioning that, from our numerical results, we propose a new bound for the consistency (relation 4.2 (4.2)), whose demonstration will remain for a future investigation. / El presente trabajo consiste en el estudio del método numérico Galerkin Discontinuo Nodal aplicado a la ecuación de advección lineal unidimensional, el cual aproxima la solución global, particionando su dominio en elementos. En cada elemento se aproxima la solución local usando interpolación; de tal manera que la solución numérica total es una suma directa de dichas aproximaciones (polinomios). El método busca alcanzar un alto orden mediante una implementación sencilla. Este modelo es estudiado por Hesthaven y Warburton[16], con la particularidad de Fusionar lo mejor del método de Volúmenes Finitos con lo mejor del método de Elementos Finitos . Primero se revisan en detalle los principales resultados sobre los polinomios ortogonales de Jacobi; más precisamente, su fórmula de generación y otros resultados que ayudan en la implementación del método. Se estudian los conceptos de interpolación y mejor aproximación. Además, se revisan algunas nociones de interpolación de espacios de Sobolev. Segundo, se detallan aspectos teóricos del método, así como su funcionamiento. Tercero, se brinda en detalle tanto la demostración de los dos teoremas de consistencia del método (mejor aproximación e interpolación) propuestos en Canuto y Quarteroni[4] como el comportamiento del error basado en Hesthaven y Warburton [16] . Finalmente, se veri ca numéricamente, mediante el uso del lenguaje Python, el teorema de consistencia referido a interpolación, así como el comportamiento del error. Se propone una nueva cota para el consistencia (relación (4.2)) basados en los resultados numéricos, cuya demostración quedará para una futura investigación. / Tesis

Doctor en Ciencias de la Ingeniería, Mención Modelación Matemática / En este trabajo se aborda el problema de encontrar soluciones regulares para algunas EDPs elípticas e hipo-elípticas no lineales y estudiar sus propiedades cualitativas. En una primera etapa, se considera la ecuación $$ -\Delta u = \lambda e^u, $$ $\lambda > 0$, en un dominio exterior con condición de Dirichlet nula. Un esquema de reducción finito-dimensional permite encontrar infinitas soluciones regulares cuando $\lambda$ es suficientemente pequeño. En la segunda parte se estudia la existencia de soluciones de la ecuación no local $$ (-\Delta)^s u = u^{p \pm \epsilon}, u > 0, $$ en un dominio acotado y suave, con condición de Dirichlet nula; donde $s > 0$ y $p:=(N+2s)/(N-2s) \pm \epsilon$ es cercano al exponente crítico ($\epsilon > 0$ pequeño). Para hallar soluciones, se utiliza un esquema de reducción finito-dimensional en espacios de funciones adecuados, donde el término principal de la función reducida se expresa a partir de las funciones de Green y de Robin del dominio. La existencia de soluciones dependerá de la existencia de puntos críticos de este término principal y de una condición de no degeneración. Por último, se considera un problema no local en el grupo de Heisenberg $H$. En particular, se buscan propiedades de rigidez para soluciones estables de $$ (-\Delta_H)^s v = f(v) en H, $$ $s \in (0,1)$. Como paso fundamental, se prueba una desigualdad del tipo Poincaré en conexión con un problema elíptico degenerado en $R^4_+$. Esta desigualdad se usará en un procedimiento de extensión para dar un criterio bajo el cual los conjuntos de nivel de las soluciones del problema anterior son superficies mínimas en $H$, es decir, tienen $H$-curvatura media nula.

En general, encontrar una solución analítica de una ecuación diferencial parcial no es fácil, y más aún cuando ésta ecuación es no lineal. Debido a esto, surgieron varios métodos numéricos para encontrar una solución aproximada a la deseada. Los métodos numéricos más conocidos son: • Métodos de Diferencias Finitas que tuvo su gran auge en la década de 1950. • Métodos de Elementos Finitos que tuvo su gran auge en la década de 1960. • Métodos Espectrales que tuvo su gran auge en la década de 1970. Mientras que los métodos de diferencias finitas dan soluciones aproximadas en los puntos de la malla computacional elegida, los métodos de elementos finitos dan aproximaciones polinomiales continuas o continuas por partes en regiones poligonales (generalmente triangulares en dos dimensiones), mientras que los métodos espectrales brindan soluciones aproximadas en la forma de polinomios sobre todo su dominio. Los métodos espectrales son una clase de discretización espacial para ecuaciones diferenciales. Las componentes claves para su formulación son las funciones base (llamadas también funciones de aproximación o expansión) y las funciones de prueba. Las funciones base se usan para dar una representación aproximada de la solución. Las funciones de prueba se usan para asegurar que la ecuación diferencial y quizás algunas condiciones de frontera se cumplan tanto como sea posible por la serie truncada de expansión. Esto se consigue minimizando, con respecto a una norma adecuada, el residuo producido por el uso de la expansión truncada en lugar de la solución exacta. Los métodos espectrales tienen un amplio uso en diferentes áreas como: teoría cuántica ([31], [36]) basado en la ecuación Schrödinger que proporciona la descripción teórica de numerosos sistemas en química y física; teoría cinética basada en la ecuación de Boltzmann ([27], [32]) o en la ecuación de Fokker-Planck ([5], [45]); problemas en mecánica de fluidos ([4], [20], [42]). También hay importantes aplicaciones en el escape átomos de la atmósfera del planeta ([14], [51]) como la pérdida de carga de partículas de la tierra ([33], [43]) y del sol [11]. El presente trabajo pretende contribuir en sentar los fundamentos sobre métodos espectrales, para que sean aplicados en futuras investigaciones más elaboradas, así como brindar los códigos de implementación (en Matlab), los cuales raramente se encuentran en forma explícita en la literatura. Este trabajo está organizado de la siguiente manera: el Capítulo 1 abarca las propiedades más importantes de los polinomios ortogonales; en particular, los polinomios de Chebyshev, los cuales son adecuados para representar funciones de dominio finito y sus relaciones de recurrencia asociadas. Además, se presenta un breve repaso de las fórmulas de cuadratura gaussiana. En el Capítulo 2, se presenta en forma detallada los métodos espectrales polinomiales, útiles para problemas con condiciones de frontera no periódicas. Presentamos los métodos de Galerkin, Tau y de Colocación. En el Capítulo 3 se da ejemplos de la implementación numérica de la ecuación del calor usando los métodos de diferencias finitas y los métodos espectrales, usando los polinomios de Chebyshev. Además, se brindan los detalles necesarios para implementar la ecuación de Burger usando los métodos espectrales. / Tesis

Doctor en Ciencias de la Ingeniería, Mención Modelación Matemática / Esta tesis está dedicada al estudio del sistema de Boussinesq estacionario: \begin{subequations}\label{sum_sp:eqn_Boussinesq} \begin{equation} -\nu \Delta\vu +(\vu\cdot\nabla)\vu+\nabla \pi=\theta\vg \text{\quad en $\Omega$,}\qquad \div\;\vu=0 \text{\quad en $\Omega$,} \end{equation} \begin{equation} -\kappa \Delta\theta +\vu\cdot\nabla\theta=h \text{\quad en $\Omega$,} \end{equation} \end{subequations} donde $\Omega\subset\R{3}$ es un conjunto abierto, acotado y conexo; $\vu$, $\pi$ y $\theta$ representan el campo de velocidades, la presión y la temperatura del fluido, respectivamente, siendo éstas las incógnitas del sistema; $\nu>0$ es la viscosidad cinemática del fluido, $\kappa>0$ es la difusividad térmica del fluido, $\vg$ es la aceleración de la gravedad y $h$ es una fuente de calor aplicada al fluido. El objetivo de esta tesis es el estudio de la teoría $L^p$ para el sistema de Boussinesq estacionario considerando dos diferentes tipos de condiciones de frontera del campo de velocidades. En efecto, en una primera etapa, se considerará la condición de frontera de Dirichlet no homogéneo \begin{equation}\label{sum_sp:cond_Dirichlet_velocity} \vu=\vub\text{\quad sobre\quad}\Gamma, \end{equation} donde $\Gamma$ denota la frontera del dominio; mientras que en una segunda etapa, el campo de velocidades tendrá prescrito la condición de frontera de Navier no homogéneo \begin{equation}\label{sum_sp:cond_Navier_velocity} \vu\cdot\vn=0,\quad 2\left[\DT(\vu)\vn\right]_{\vt}+\alpha\;\vu_{\vt}=\bm{a},\text{\quad sobre\quad}\Gamma, \end{equation} donde $\DT(\vu)=\frac{1}{2}\left(\nabla\vu+(\nabla\vu)^T\right)$ es el tensor de deformación asociado con el campo de velocidades $\vu$, $\vn$ es el vector normal unitario exterior, $\vt$ es el correspondiente vector unitario tangente, $\alpha$ y $\vNb$ son una función de fricción y un campo vectorial tangencial definidas ambas sobre la frontera. Además, la condición de frontera para la temperatura será, en las dos primeras partes, la condición de frontera de Dirichlet no homogéneo \begin{equation}\label{sum_sp:cond_Dirichlet_temperature} \theta=\thb\text{\quad sobre\quad}\Gamma. \end{equation} Así, en primer lugar, estudiamos la existencia y unicidad de la solución débil para el problema \eqref{sum_sp:eqn_Boussinesq}, \eqref{sum_sp:cond_Dirichlet_velocity} y \eqref{sum_sp:cond_Dirichlet_temperature} en el caso hilbertiano. Además, la existencia de soluciones generalizadas para $p\geq\frac{3}{2}$ y soluciones fuertes para $1<p<\infty$ es probada. También, se estudiará la existencia y unicidad de la solución muy débil. Vale la pena señalar que debido a que la condición de Dirichlet no homogénea es considerada para la velocidad, el hecho de que la frontera del dominio pueda ser no conexa juega un papel importante, ya que aparece de manera explícita en las hipótesis de algunos de los principales resultados. Por otro lado, en la segunda etapa de la tesis, se estudiará la existencia de soluciones débiles en el caso de Hilbert, así como la existencia de soluciones generalizadas para $p>2$ y soluciones fuertes para $p\geq\frac{6}{5}$ para el problema \eqref{sum_sp:eqn_Boussinesq}, \eqref{sum_sp:cond_Navier_velocity} y \eqref{sum_sp:cond_Dirichlet_temperature}. Tenga en cuenta que la suposición hecha anteriormente acerca de la no conexidad de la frontera no aparecerá aquí debido a la restricción de impermeabilidad en la frontera. Finalmente, en la última parte de esta tesis, estudiamos la teoría $L^p$ para las ecuaciones de Stokes con la condición de Navier \eqref{sum_sp:cond_Navier_velocity}. Más precisamente, nos ocuparemos de la regularidad $W^{1,p}$ para $p\geq2$ y la regularidad $W^{2,p}$ para $p\geq\frac{6}{5}$.

Resuelve y analiza el cálculo de la solución numérica de la ecuación de advección - difusión mediante el esquema de diferencias finitas. Presenta criterios de estabilidad, los cuales garantizan la estabilidad del esquema planteado para resolver la ecuacion de adveccion - difusión longitudinal; los criterios obtenidos deben garantizar la estabilidad y convergencia, como función de los números de Couran y Péclet, con todo esto, se tiene a disposición métodos sencillos que son numéricamente estables y convergentes, por lo que no se considera necesario recurrir a métodos más complicados para resolver la ecuación de advección - difusión para el caso unidimensional. Además, realiza la simulación computacional de la solución numérica de la ecuación de advección - difusión longitudinal utilizando el software Matlab y el lenguaje de programación Python. / Trabajo de suficiencia profesional

Doctor en Ciencias de la Ingeniería, Mención Modelación Matemática / Este trabajo se concentra principalmente en estudiar el método de reducción de Lyapunov-Schmidt y sus aplicaciones al estudio de existencia de soluciones a problemas semilineales elípticos. En particular, utilizamos exitosamente este método para estudiar la ecuación de Allen-Cahn \Delta u + u(1-u^2)=0, en R^N en diferentes contextos. La geometría de los conjuntos de nivel de soluciones enteras de esta ecuación, presenta una estructura variada y compleja. En particular, esta ecuación esta presente en la famosa conjetura de E. De Giorgi, la cual afirma que si la dimensión del espacio es tal que 2\leq N\leq 8, las soluciones acotadas de esta ecuación que son monótonas en una dirección, tienen por conjuntos nivel a una familia de hiperplanos paralelos entre si, es decir, la solución depende solo de una variable. Gran progreso se ha alcanzado en la demostración de esta conjetura durante las \'ultimas décadas. La monotonía de las soluciones esta relacionada con sus propiedades de estabilidad. En el programa de entender el conjunto de soluciones enteras de esta ecuación, es interesante estudiar soluciones que tienes índice de Morse finito, de las cuales para nuestro conocimiento, pocos ejemplos se conocen hasta ahora. En la primera parte de esta investigación, utilizamos el método de reducción, en esencia no variacional, para construir una familia de soluciones acotadas axialmente simétricas a la ecuación de Allen-Cahn en R3, con la propiedad de tener múltiples transiciones sobre una dilatación grande de una catenoide. De nuestro desarrollo, se evidencia contundentemente que estas soluciones tienen indice de Morse grande a medida que la catenoide se vuelve más y más dilatada. Motivados por este descubrimiento y utilizando el mismo método, continuamos este trabajo construyendo una nueva familia de soluciones axialmente simétricas a la ecuación de Allen-Cahn en R3, cuyo conjunto nodal consiste en dos componentes conexas que provienen del grafo y su reflexión respecto al eje z, de una solución suave y radialmente simétrica de la ecuación de Liouville en R2. De igual forma, encontramos fuerte evidencia para afirmar que el índice de Morse de esta familia de soluciones es finito. Luego, presentamos el estudio de la ecuación no homogénea de Allen-Cahn en R2, en la cual presentamos otra aplicación del método reducción construyendo, bajo ciertas condiciones geométricas, una familia de soluciones cuyos conjuntos nodales, fuera de una bola grande de R2, tienen dos componentes conexas que son asintóticamente semirrectas no paralelas entre si. Finalmente, y en contraste, consideramos el contexto variacional presentando resultados de existencia de múltiples soluciones para un sistema elíptico de ecuaciones con un acoplamiento simétrico. La aplicación del método de reducción variacional, permite luego aplicar de forma clásica el teorema de paso de montaña simétrico. La importancia del método de reducción, en este caso, radica en que las propiedades de simetría del sistema de ecuaciones, las cuales provienen de la forma del sistema, en lugar de las no linealidades, son heredadas por ecuación reducida.

Estudia la existencia y unicidad de la solución débil para una ecuación de evolución de segundo orden, presentados en dos casos, semi lineal y lineal, obteniendo regularización de la solución débil para el caso semi lineal y la dependencia continua sobre los datos iniciales para el caso lineal. Para se utiliza el método de Faedo-Galerkin y la igualdad de la Energía, esta última está basada del libro Problemas Aux Limites non Homogéneas et Applications, volumen 1, de Lions J. y Magenes E. Finalmente, se abordan algunas aplicaciones a las ecuaciones diferenciales parciales. / Tesis

Doctor en Ciencias de la Ingeniería, Mención Modelación Matemática / Esta tesis esta dividida en seis partes. La primera parte está dedicada a probar propiedades de Hadamard y teoremas del tipo de Liouville para soluciones viscosas de ecuaciones diferenciales parciales elípticas completamente no lineales con término gradiente \begin{equation}\label{eq06-10-13 1} \mathcal{M}^{-}(|x|,D^2u)+\sigma(|x|)|Du|+f(x,u)\leq 0,\quad \ x\in\Omega, \end{equation} donde $\Omega=\mathbb{R}^N$ o un dominio exterior, las funciones $\sigma:[0,\infty)\to\mathbb{R}$ y $f:\Omega\times (0,\infty)\to (0,\infty)$ son continuas las cuales satisfacen algunas condiciones extras. En la segunda parte se estudia la existencia de soluciones que explotan en la frontera para ecuaciones elípticas fraccionarias semilineales \begin{equation}\label{eq06-10-13 2} \arraycolsep=1pt \begin{array}{lll} (-\Delta)^{\alpha} u(x)+|u|^{p-1}u(x)=h(x),\quad & x\in\Omega,\\[2mm] \phantom{ (-\Delta)^{\alpha} u(x)+|u|^{p-1}} u(x)=0,\quad & x\in\bar\Omega^c,\\[2mm] \phantom{ (-\Delta)^{\alpha} \ } \lim_{x\in\Omega, x\to\partial\Omega}u(x)=+\infty, \end{array} \end{equation} donde $p>1$, $\Omega$ es un dominio abierto acotado $C^2$ de $\mathbb{R}^N(N\geq2)$, el operador $(-\Delta)^{\alpha}$ con $\alpha\in(0,1)$ es el Laplaciano fraccionario y $h:\Omega\to\R$ es una función continua la cual satisface algunas condiciones extras. Por otra parte, analizamos la unicidad y el comportamiento asimptótico de soluciones al problema (\ref{eq06-10-13 2}). El objetivo principal de la tercera parte es investigar soluciones positivas para ecuaciones elípticas fraccionarias \begin{equation}\label{eq06-10-13 3} \arraycolsep=1pt \begin{array}{lll} (-\Delta)^{\alpha} u(x)+|u|^{p-1}u(x)=0,\quad & x\in\Omega\setminus\mathcal{C},\\[2mm] \phantom{ (-\Delta)^{\alpha} u(x)+|u|^{p-1}} u(x)=0,\quad & x\in\Omega^c,\\[2mm] \phantom{ (-\Delta) \ } \lim_{x\in\Omega\setminus\mathcal{C}, \ x\to\mathcal{C}}u(x)=+\infty, \end{array} \end{equation} donde $p>1$ y $\Omega$ es un dominio abierto acotado $C^2$ de $\mathbb{R}^N(N\geq2)$, $\mathcal{C}\subset \Omega$ es el frontera de dominio $G$ que es $C^2$ y satisface $\bar G\subset\Omega$. Consideramos la existencia de soluciones positivas para el problema (\ref{eq06-10-13 3}). Mas aún, analizamos la unicidad, el comportamiento asimptótico y la no existencia al problema (\ref{eq06-10-13 3}). En la cuarta parte, estudiamos la existencia de soluciones débiles de (F) $ (-\Delta)^\alpha u+g(u)=\nu $ en un dominio $\Omega$ abierto acotado $C^2$ de $\R^N (N\ge2)$ el cual se desvanece en $\Omega^c$, donde $\alpha\in(0,1)$, $\nu$ es una medida de Radon y $g$ es una función no decreciente satisfaciendo algunas hipótesis extras. Cuando $g$ satisface una condición de integrabilidad subcrítica, probamos la existencia y unicidad de una solución débil para el problema (F) para cualquier medida. En el caso donde $\nu$ es una masa de Dirac, caracterizamos el comportamiento asimptótico de soluciones a (F). Asimismo, cuando $g(r)=|r|^{k-1}r$ con $k$ supercrítico, mostramos que una condición de absoluta continuidad de la medida con respecto a alguna capacidad de Bessel es una condición necesaria y suficiente para que (F) sea resuelta. El propósito de la quinta parte es investigar soluciones singulares débiles y fuertes de ecuaciones elípticas fraccionarias semilineales. Sean $p\in(0,\frac{N}{N-2\alpha})$, $\alpha\in(0,1)$, $k>0$ y $\Omega\subset \R^N(N\geq2)$ un dominio abierto acotado $C^2$ conteniendo a $0$ y $\delta_0$ la masa de Dirac en $0$, estudiamos que la solución débil de $(E)_k$ $ (-\Delta)^\alpha u+u^p=k\delta_0 $ en $\Omega$ la cual se desvanece en $\Omega^c$ es una solución débil singular de $(E^*)$ $ (-\Delta)^\alpha u+u^p=0 $ en $\Omega\setminus\{0\}$ con el mismo dato externo. Por otra parte, estudiamos el límite de soluciones débiles de $(E)_k$ cuando $k\to\infty$. Para $p\in(0, 1+\frac{2\alpha}{N}]$, el límite es infinito en $\Omega$. Para $p\in(1+\frac{2\alpha}N,\frac{N}{N-2\alpha})$, el límite es una solución fuertemente singular de $(E^*)$. Finalmente, en la sexta parte estudiamos la ecuación elíptica fraccionaria semilineal (E1) $(-\Delta)^\alpha u+\epsilon g(|\nabla u|)=\nu $ en un dominio $\Omega$ abierto acotado $C^2$ de $\R^N (N\ge2)$, el cual se desvanece en $\Omega^c$, donde $\epsilon=\pm1$, $\alpha\in(1/2,1)$, $\nu$ es una medida de Radon y $g:\R_+\mapsto\R_+$ es una funci\'on continua. Probamos la existencia de soluciones débiles para el problema (E1) cuando $g$ es subcrítico. Además, el comportamiento asimptótico y la unicidad de soluciones son descritas cuando $\epsilon=1$, $\nu$ es una masa de Dirac y $g(s)=s^p$ con $p\in(0,\frac)$.