Multifluid flows with weak and strong discontinuous interfaces using an elemental enriched space

In a previous paper, the authors presented an elemental enriched space to be used in a finite-element framework (EFEM) capable of reproducing kinks and jumps in an unknown function using a fixed mesh in which the jumps and kinks do not coincide with the interelement boundaries. In this previous publication, only scalar transport problems were solved (thermal problems). In the present work, these ideas are generalized to vectorial unknowns, in particular, the incompressible Navier-Stokes equations for multifluid flows presenting internal moving interfaces. The advantage of the EFEM compared with global enrichment is the significant reduction in computing time when the internal interface is moving. In the EFEM, the matrix to be solved at each time step has not only the same amount of degrees of freedom (DOFs) but also the same connectivity between the DOFs. This frozen matrix graph enormously improves the efficiency of the solver. Another characteristic of the elemental enriched space presented here is that it allows a linear variation of the jump, thus improving the convergence rate, compared with other enriched spaces that have a constant variation of the jump. Furthermore, the implementation in any existing finite-element code is extremely easy with the version presented here because the new shape functions are based on the usual finite-element method shape functions for triangles or tetrahedrals, and once the internal DOFs are statically condensed, the resulting elements have exactly the same number of unknowns as the nonenriched finite elements.Peer ReviewedPreprin

General treatment of essential boundary conditions in reduced order models for non-linear problems

Inhomogeneous essential boundary conditions must be carefully treated in the formulation of Reduced Order Models (ROMs) for non-linear problems. In order to investigate this issue, two methods are analysed: one in which the boundary conditions are imposed in an strong way, and a second one in which a weak imposition of boundary conditions is made. The ideas presented in this work apply to the big realm of a posteriori ROMs. Nevertheless, an a posteriori hyper-reduction method is specifically considered in order to deal with the cost associated to the non-linearity of the problems. Applications to nonlinear transient heat conduction problems with temperature dependent thermophysical properties and time dependent essential boundary conditions are studied. However, the strategies introduced in this work are of general application.Fil: Cosimo, Alejandro. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Cardona, Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Idelsohn, Sergio Rodolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina. Institució Catalana de Recerca i Estudis Avancats; España. International Center for Numerical Methods in Engineering; Españ

Modelación numérica de deslizamientos de ladera en embalses mediante el método de partículas y elementos finitos (PFEM)

Se presentan resultados obtenidos mediante la aplicación del MétododePartículasyElementosFinitos (PFEM) en la simulación dedeslizamientosdeladeraenembalses. Es un fenómeno complejo, por la interacción entre el material deslizado, la masa de agua del embalse, yel conjunto formado por el vaso y la presa. PFEM es un esquema numérico original con el que se ha afrontado con éxito la resolución de problemas de interacción fluido-estructura. Combina un enfoque Lagrangiano con la resolución de las ecuaciones deelementosfinitosmediante la generación de una malla, que se actualiza en cada paso de tiempo. Se presentan resultados de casos de validación en los que se han reproducido ensayos en laboratorio existentes en la bibliografía técnica. Se muestran también otros cálculos más complejos, sobre la cartografía a escala real de un embalse, donde se aprecia el fenómeno de generación de la ola, su propagación por el embalse y la afección a la presa. Por último, se ha modelado el deslizamiento ocurrido en 1958 en la bahía de Lituya (Alaska), enel que la caída de 90 millones de toneladas de roca produjo una ola que alcanzó una sobreelevación máxima de 524 m en la ladera opuesta. Los resultados permiten afirmar que PFEM puede ser una herramienta útil enel análisis de riesgos frente a este tipo de fenómenos, ofreciendo una buena aproximación de las afecciones potenciale

On the use of deep, shallow or flat shell finite elements for the analysis of thin shell structures

This paper is confined to the study of thin shells. The aim is to summarize the different theories used and to examine the assumptions upon which each of them is based. The intention is to show when it is more suitable to use a particular approximation and to indicate the errors it introduces. Beginning with the general deep shell theory, some simplifications are introduced to obtain the shallow shell theories. The special implications of this theory for the finite element method are also examined. Finally the particular case of flat elements is discussed

A Fully Coupled Formulation For Incompressible Fluid-Elastic Structure-Interactions

We present a general formulation for analysis of fluid flows with structural interactions using the particle finite element method (PFEM). The fluids are fully coupled to the structures that can undergo highly non-linear response due to large deformations. The key feature of the PFEM is the use of an updated Lagrangian description to model the motion of nodes(particles) in both the fluid and the structure domains. A mesh connects the nodes defining the discretized domains where the governing equations, expressed in an integral form are solved as in the standard FEM. A fractional step scheme for the transient coupled fluid-structure solution is described. Examples of application of the PFEM method to solve a number of fluidstructure interaction problems including free-fluid-surfaces, breaking waves and fluid particle separation may be easily solved with this formulation are presented.Fil: Marti, Julio Marcelo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Desarrollo Tecnológico para la Industria Química. Universidad Nacional del Litoral. Instituto de Desarrollo Tecnológico para la Industria Química; ArgentinaFil: Idelsohn, Sergio Rodolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Desarrollo Tecnológico para la Industria Química. Universidad Nacional del Litoral. Instituto de Desarrollo Tecnológico para la Industria Química; ArgentinaFil: Limache, Alejandro Cesar. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Desarrollo Tecnológico para la Industria Química. Universidad Nacional del Litoral. Instituto de Desarrollo Tecnológico para la Industria Química; Argentin

Modelling the vertical UL 94 test: competition and collaboration between melt dripping, gasification and combustion

An experimental and numerical investigation of the effect of bisphenol A bis(diphenyl phosphate) (BDP) and polytetrafluoroethylene (PTFE) on the fire behaviour of bisphenol A polycarbonate/acrylonitrile butadiene styrene (PC/ABS) in the vertical UL 94 scenario is presented. Four PC/ABS blends were discussed, which satisfy different UL 94 classifi cations d ue to the  competing ef fects of gasifica ti on, charring, flame inhibition and melt flow/dripping. For numerical investigation, the particle finite element method (PFEM) is used. Its capability to model the complex fire behaviour of polymers in the UL 94 is analysed. The materials’ properties are characterised, in particular the additives impact on the dripping behaviour during thermal exposure. BDP is an efficie nt p lasticiser; adding PTFE p reve nts dripping  by causing a flo w limit. P FEM simulation s reproduce the dripping and burning behaviour, in particular the competition between gasification and dripping. The thermal impact of both the burner and the flame is approximated taking into account flame inhibition, charring and effective heat of combustion. PFEM is a promising numerical tool for the investigation of the fire behaviour of polymers, particularly when large deformations are involved. Not only the principal phenomena but also the different UL 94 classi fi cations and t he exti nc tion times are well predicted

Consistent pressure Laplacian stabilization for incompressible continua via higher‐order finite calculus

We present a stabilized numerical formulation for incompressible continua based on a higher‐order Finite Calculus (FIC) approach and the finite element method. The focus of the paper is on the derivation of a stabilized form for the mass balance (incompressibility) equation. The simpler form of the momentum equations neglecting the non‐linear convective terms, which is typical for incompressible solids, Stokes flows and Lagrangian flows is used for the sake of clarity. The discretized stabilized mass balance equation adds to the standard divergence of velocity term a pressure Laplacian and an additional boundary term. The boundary term is relevant for the accuracy of the numerical solution, especially for free surface flow problems. The Laplacian and boundary stabilization terms are multiplied by non‐linear parameters that have an extremely simple expression in terms of element sizes, the pressure and the discrete residuals of the incompressibility equation and the momentum equations, thus ensuring the consistency of the method. The stabilized formulation allows solving the incompressible problem iteratively using an equal‐order interpolation for the velocities (or displacements) and the pressure, which are the only unknowns. The use of additional pressure gradient projection variables, typical of many stabilized methods, is unnecessary. The formulation is particularly useful for heterogeneous incompressible materials with discontinuous material properties, as it allows computing all the stabilization matrices at the element level. Details of the finite element formulation are given. The good behaviour of the new pressure Laplacian stabilization (PLS) technique is shown in simple but demonstrative examples of application. A very accurate solution was obtained in all cases in 2–3 iterations. Copyright © 2010 John Wiley & Sons, Ltd

Towards a New Turbulence Model Based on Lagrangian Flows, Reduced Order Models, and Global-Local Approximations

Although the Navier-Stokes equations represent correctly both, laminar fluid flows as well as turbulent ones, the current power of the computers does not allow solving the latter without making empirical approximations that make the results are only predictable within the same margins from which the empirical approximation was realized. Taking into account that the vast majority of fluid flows that must be simulated by the industries are indeed of a turbulent nature, this makes it worthwhile to continue improving the models so that they fit more and more with the physicsmathematic equations. The current project in which we are working is in that sense. Based on previous work, we are proposing a new turbulence model that fits more with the physic and can be solved in one of the current computer. Always within the multi-scale models, in which all the turbulence model are based, the new idea consists in: a) treating the macro-scale with lagrangian particles, which convect and diffuse the turbulence; b) solving the micro-scale as a problem of unstable laminar fluid flow subjected to velocity gradients; c) applying a Reduction Order Model (ROM) to the micro-scale; d) going from the microscale to the macro-scale as in a Global-Local Model with a ROM.Publicado en: Mecánica Computacional, vol. XXXV, no 1.Facultad de Ingenierí