Finite element approximation of the viscoelastic flow problem: a non-residual based stabilized formulation

In this paper, a three-field finite element stabilized formulation for the incompressible viscoelastic fluid flow problem is tested numerically. Starting from a residual based formulation, a non-residual based one is designed, the benefits of which are highlighted in this work. Both formulations allow one to deal with the convective nature of the problem and to use equal interpolation for the problem unknowns View the MathML sources-u-p (deviatoric stress, velocity and pressure). Additionally, some results from the numerical analysis of the formulation are stated. Numerical examples are presented to show the robustness of the method, which include the classical 4: 1 planar contraction problem and the flow over a confined cylinder case, as well as a two-fluid formulation for the planar jet buckling problem.Peer ReviewedPostprint (author's final draft

An adaptive fixed-mesh ALE method for free surface flows

In this work we present a Fixed-Mesh ALE method for the numerical simulation of free surface flows capable of using an adaptive finite element mesh covering a background domain. This mesh is successively refined and unrefined at each time step in order to focus the computational effort on the spatial regions where it is required. Some of the main ingredients of the formulation are the use of an Arbitrary-Lagrangian–Eulerian formulation for computing temporal derivatives, the use of stabilization terms for stabilizing convection, stabilizing the lack of compatibility between velocity and pressure interpolation spaces, and stabilizing the ill-conditioning introduced by the cuts on the background finite element mesh, and the coupling of the algorithm with an adaptive mesh refinement procedure suitable for running on distributed memory environments. Algorithmic steps for the projection between meshes are presented together with the algebraic fractional step approach used for improving the condition number of the linear systems to be solved. The method is tested in several numerical examples. The expected convergence rates both in space and time are observed. Smooth solution fields for both velocity and pressure are obtained (as a result of the contribution of the stabilization terms). Finally, a good agreement between the numerical results and the reference experimental data is obtained.Postprint (published version

Palabras para un amigo. El triste día.

Cali, Enero 18 de 2014 Todo empezó con una llamada a las 7:49 de la mañana. Yo tenía la resaca del día anterior y no quise contestar. Por lo tanto, para que no me molestaran más, puse mi celular en silencio. No quería ser interrumpido de ese sueño inútil que no servía para sacarme los litros de alcohol que había ingresado a mi cuerpo la noche anterior. Y como una señal divina, hubo un sol que penetró las maderas de mi persiana, un calor que ni el aire podía opacar. Mi garganta se secó tanto que salí de mi cuarto a beber un poco de agua. Recuerdo mucho que ese día, no tomé un vaso sino que agarré la jarra pues estaba solo en mi casa y nadie me regañaría. Cuando volví a mi cuarto, vi que tenía treinta y siete llamadas perdidas. De inmediato, me asuste. Le escribí a cada uno de los que me habían llamado preguntándoles qué pasó. Garzón fue el primero en responder. Recuerdo muy bien que me dijo “ya te llamo”. A continuación, vino una serie de palpitaciones seguidas por el suspenso de la llamada, dejé que el teléfono sonara una vez pues no aguantaba las ganas de saber qué era lo que había ocurrido

Dynamic term-by-term stabilized finite element formulation using orthogonal subgrid-scales for the incompressible Navier–Stokes problem

In this paper, we propose and analyze the stability and the dissipative structure of a new dynamic term-by-term stabilized finite element formulation for the Navier–Stokes problem that can be viewed as a variational multiscale (VMS) method under some assumptions. The essential point of the formulation is the time dependent nature of the subscales and, contrary to residual-based formulations, the introduction of two velocity subscale components. They represent the components of the convective and the pressure gradient terms, respectively, of the momentum equation that cannot be captured by the finite element mesh. A key idea of the proposed method is that the convective subscale is close to a solenoidal field and the pressure gradient subscale is close to a potential field. The method ensures stability in anisotropic space–time discretizations, which is proved using numerical analysis for a linearized problem and demonstrated in classical numerical tests. The work includes a detailed description of the proposed formulation and several numerical examples that serve to justify our claims.Peer ReviewedPostprint (author's final draft

Stabilized finite element formulations for the three-field viscoelastic fluid flow problem

The Finite Element Method (FEM) is a powerful numerical tool, that permits the resolution of problems defined by partial differential equations, very often employed to deal with the numerical simulation of multiphysics problems. In this work, we use it to approximate numerically the viscoelastic fluid flow problem, which involves the resolution of the standard Navier-Stokes equations for velocity and pressure, and another tensorial reactive-convective constitutive equation for the elastic part of the stress, that describes the viscoelastic nature of the fluid. The three-field (velocity-pressure-stress) mixed formulation of the incompressible Navier-Stokes problem, either in the elastic and in the non-elastic case, can lead to two different types of numerical instabilities. The first is associated with the incompressibility and loss of stability of the stress field, and the second with the dominant convection. The first type of instabilities can be overcome by choosing an interpolation for the unknowns that satisfies the two inf-sup conditions that restrict the mixed problem, whereas the dominant convection requires a stabilized formulation in any case. In this work, different stabilized schemes of the Sub-Grid-Scale (SGS) type are proposed to solve the three-field problem, first for quasi Newtonian fluids and then for solving the viscoelastic case. The proposed methods allow one to use equal interpolation for the problem unknowns and to stabilize dominant convective terms both in the momentum and in the constitutive equation. Starting from a residual based formulation used in the quasi-Newtonian case, a non-residual based formulation is proposed in the viscoelastic case which is shown to have superior behavior when there are numerical or geometrical singularities. The stabilized finite element formulations presented in the work yield a global stable solution, however, if the solution presents very high gradients, local oscillations may still remain. In order to alleviate these local instabilities, a general discontinuity-capturing technique for the elastic stress is also proposed. The monolithic resolution of the three-field viscoelastic problem could be extremely expensive computationally, particularly, in the threedimensional case with ten degrees of freedom per node. A fractional step approach motivated in the classical pressure segregation algorithms used in the two-field Navier-Stokes problem is presented in the work.The algorithms designed allow one the resolution of the system of equations that define the problem in a fully decoupled manner, reducing in this way the CPU time and memory requirements with respect to the monolithic case. The numerical simulation of moving interfaces involved in two-fluid flow problems is an important topic in many industrial processes and physical situations. If we solve the problem using a fixed mesh approach, when the interface between both fluids cuts an element, the discontinuity in the material properties leads to discontinuities in the gradients of the unknowns which cannot be captured using a standard finite element interpolation. The method presented in this work features a local enrichment for the pressure unknowns which allows one to capture pressure gradient discontinuities in fluids presenting different density values. The stability and convergence of the non-residual formulation used for viscoelastic fluids is analyzed in the last part of the work, for a linearized stationary case of the Oseen type and for the semi-discrete time dependent non-linear case. In both cases, it is shown that the formulation is stable and optimally convergent under suitable regularity assumptions.El Método de los Elementos Finitos (MEF) es una herramienta numérica de gran alcance, que permite la resolución de problemas definidos por ecuaciones diferenciales parciales, comúnmente utilizado para llevar a cabo simulaciones numéricas de problemas de multifísica. En este trabajo, se utiliza para aproximar numéricamente el problema del flujo de fluidos viscoelásticos, el cual requiere la resolución de las ecuaciones básicas de Navier-Stokes y otra ecuación adicional constitutiva tensorial de tipo reactiva-convectiva, que describe la naturaleza viscoelástica del fluido. La formulación mixta de tres campos (velocidad-presión-tensión) del problema de Navier-Stokes, tanto en el caso elástico como en el no-elástico, puede conducir a dos tipos de inestabilidades numéricas. El primer grupo, se asocia con la incompresibilidad del fluido y la pérdida de estabilidad del campo de tensiones, y el segundo con la convección dominante. El primer tipo de inestabilidades, se puede solucionar eligiendo un tipo de interpolación entre las incógnitas que satisfaga las dos condiciones inf-sup que restringen el problema mixto, mientras que la convección dominante requiere del uso de formulaciones estabilizadas en cualquier caso. En el trabajo, se proponen diferentes esquemas estabilizados del tipo SGS (Sub-Grid-Scales) para resolver el problema de tres campos, primero para fluidos del tipo cuasi-newtonianos y luego para resolver el caso viscoelástico. Los métodos estabilizados propuestos permiten el uso de igual interpolación entre las incógnitas del problema y estabilizan la convección dominante, tanto en la ecuación de momento como en la ecuación constitutiva. Comenzando desde una formulación de tipo residual usada en el caso cuasi-newtoniano, se propone una formulación no-residual para el caso viscoelástico que muestra un comportamiento superior en presencia de singularidades numéricas y geométricas. En general, una formulación estabilizada produce una solución estable global, sin embargo, si la solución presenta gradientes elevados, oscilaciones locales se pueden mantener. Con el objetivo de aliviar este tipo de inestabilidades locales, se propone adicionalmente una técnica general de captura de discontinuidades para la tensión elástica. La resolución monolítica del problema de tres campos viscoelástico puede llegar a ser extremadamente costosa computacionalmente, sobre todo, en el caso tridimensional donde se tienen diez grados de libertad por nodo. Un enfoque de paso fraccionado motivado en los algorítmos clásicos de segregación de la presión usados en el caso del problema de dos campos de Navier-Stokes, se presenta en el trabajo, el cual permite la resolución del sistema de ecuaciones que definen el problema de una manera completamente desacoplada, lo que reduce los tiempos de cálculo y los requerimientos de memoria, respecto al caso monolítico. La simulación numérica de interfaces móviles que envuelve los problemas de dos fluidos, es un tópico importante en un gran número de procesos industriales y situaciones físicas. Si se resuelve el problema utilizando un enfoque de mallas fijas, cuando la interfaz que separa los dos fluidos corta un elemento, la discontinuidad en las propiedades materiales da lugar a discontinuidades en los gradientes de las incógnitas que no pueden ser capturados utilizando una formulación estándar de interpolación. Un enriquecimiento local para la presión se presenta en el trabajo, el cual permite la captura de gradientes discontinuos en la presión, asociados a fluidos de diferentes densidades. La estabilidad y la convergencia de la formulación no-residual utilizada para fluidos viscoelásticos es analizada en la última parte del trabajo, para un caso linealizado estacionario del tipo Oseen y para un problema transitorio no-lineal semi-discreto. En ambos casos, se logra mostrar que la formulación es estable y de convergencia óptima bajo supuestos de regularidad adecuados.Postprint (published version

Numerical analysis of a stabilized finite element approximation for the three-field linearized viscoelastic fluid problem using arbitrary interpolations

The original publication is available at www.esaimm2an.org.In this paper we present the numerical analysis of a three-field stabilized finite element formulation recently proposed to approximate viscoelastic flows. The three-field viscoelastic fluid flow problem may suffer from two types of numerical instabilities: on the one hand we have the two inf-sup conditions related to the mixed nature problem and, on the other, the convective nature of the momentum and constitutive equations may produce global and local oscillations in the numerical approximation. Both can be overcome by resorting from the standard Galerkin method to a stabilized formulation. The one presented here is based on the subgrid scale concept, in which unresolvable scales of the continuous solution are approximately accounted for. In particular, the approach developed herein is based on the decomposition into their finite element component and a subscale, which is approximated properly to yield a stable formulation. The analyzed problem corresponds to a linearized version of the Navier-Stokes/Oldroyd-B case where the advection velocity of the momentum equation and the non-linear terms in the constitutive equation are treated using a fixed point strategy for the velocity and the velocity gradient. The proposed method permits the resolution of the problem using arbitrary interpolations for all the unknowns. We describe some important ingredients related to the design of the formulation and present the results of its numerical analysis. It is shown that the formulation is stable and optimally convergent for small Weissenberg numbers, independently of the interpolation used.Peer ReviewedPostprint (author's final draft

Sensitivity analysis in calculus of variations. Some applications

This paper deals with the problem of sensitivity analysis in calculus of variations. A perturbation technique is applied to derive the boundary value problem and the system of equations that allow us to obtain the partial derivatives (sensitivities) of the objective function value and the primal and dual optimal solutions with respect to all parameters. Two examples of applications, a simple mathematical problem and a slope stability analysis problem, are used to illustrate the proposed method

La condizione in diritto romano: tra negozio e rapporto giuridico.

Il lavoro si propone di analizzare alcuni brani, prevalentemente tratti dal Digesto giustinianeo, che si sono occupati della definizione della disciplina della clausola condizionale e del negozio condizionato. La riflessione che ne risulta riguarda il ruolo svolto dalla clausola condizionale nell'evoluzione del diritto privato, da diritto legato all'atto giuridico, a diritto costruito sugli assetti di interessi concretamente perseguiti dalle parti

Cumplimiento de los estándares de calidad ambiental y su impacto en la liquidez de una empresa de servicios ambientales, 2016-2019

La presente investigación tuvo como objetivo determinar el impacto del cumplimiento de los estándares de calidad ambiental en la liquidez de una empresa de servicios ambientales. Fue una investigación aplicada de diseño no experimental transeccional-causal y la muestra estuvo conformada por una empresa del sector de servicios ambientales, los instrumentos fueron la recolección de datos como los estados financieros, el análisis de datos estuvo conformado por el cálculo de ratios de liquidez y un modelo de regresión. Dentro de los principales resultados se encontró que la empresa en estudio sigue presentando problemas de liquidez. Siendo la capacidad de solvencia a sus obligaciones a corto plazo, a razón de 0,98 veces, podemos concluir que la empresa corre riesgos a suspender pagos, apenas puede cubrir sus obligaciones, aplicando un análisis de regresión se encontró que el p-valor para mi variable dependiente es mayor que el 5%, por los cual rechazamos la hipótesis alterna y aceptamos la hipótesis nula que niega el efecto del cumplimiento de los estándares ambientales. Se concluyo que los estándares de calidad ambiental tienen un impacto sobre la liquidez de la empresa en un 4,5%, aunque es bajo pero tiene una leve tendencia positiva en los últimos años