A diffusion-split method to deal with thermal shocks using standard linear tetrahedral finite elements

International audienceThe thermal analysis using linear standard tetrahedral finite elements may be affected by spurious local extrema in the regions affected by thermal shocks, in such a severe way to directly discourage the use of these elements. The present work proposes a slight modification to the discrete heat equation in order to obtain a system matrix in M-matrix form, which assures an oscillation-free solution. The performance of this method is evaluated by means of test case with analytical solution, as well as an industrial application, for which a well-behaved numerical solution is available

A computational multi-objective optimization method to improve energy efficiency and thermal comfort in dwellings

In the last years, multi-objective optimization techniques became into one of the main challenges of the building energy efficiency area. The objective of this paper is to develop and validate a computational code for multi-objective building performance optimization by linking an evolutionary algorithm and a building simulation software in a powerful cluster. A sophisticated version of the multi-objective Non-dominated Sorting Genetic Algorithm-II (NSGA-II) was implemented in Python code to determine the optimal building design, which allows working with categorical and discrete variables, and the objectives were evaluated using the building energy simulation software EnergyPlus. NSGA-II was implemented to run in a high-performance cluster for the parallel computing of the fitness of each population (set of possible designs). In this work, the strengths of the proposed method were demonstrated by its application to the optimal design of a typical single-family house, located in the Argentine Littoral region. This house has some rooms conditioned only by natural ventilation, and other rooms with natural ventilation supplemented by mechanical air-conditioning (hybrid ventilation). The most influential design variables like roof types, external and internal wall types, solar orientation, solar absorptance, size, type, and windows shading of this house among others were studied in two complex cases of 108 and 1016 possibilities to obtain the best trade-off (Pareto front) between heating and cooling performance. Finally, a decision-making method was applied to select one configuration of the Pareto front. Optimal simulation results for the study cases indicated that is possible to improve up to 95% the thermal comfort in naturally ventilated rooms and up to 82% energy performance in air-conditioned rooms of the building with respect to the original configuration by using a design that takes simultaneous advantage of passive strategies like thermal inertia and natural ventilation. The methodology was proved to give a robust and powerful tool to design efficient dwellings reducing the optimization time from almost 12 days to 4.4 h.Fil: Bre, Facundo. 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. Universidad Tecnológica Nacional. Facultad Regional Concepción del Uruguay; ArgentinaFil: Fachinotti, Victor Daniel. 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; Argentin

A two-phase two-dimensional finite element thermomechanics and macrosegregation model of mushy zone. Application to continuous casting

International audienceThe main lines of a coupled thermomechanical - solute transport model are first summarized. Macroscopic conservation equations for mass, momentum, energy and solute are obtained by a spatial averaging method. The mechanical model is a "sponge-like" one: assuming a semi-solid saturated mushy zone, the solid phase is macroscopically modeled as a compressible viscoplastic continuum, while the liquid phase flow obeys Darcy's law. Regarding solute transport, the study is limited to a binary alloy for which the solidification path is not given a priori but results from a microsegregation model (here the lever rule). A validation check of the correct implementation of this coupled model is achieved by comparing with an analytical solution in the case of a free compression of a saturated semi-solid medium. Application to the study of the solidification during secondary cooling in steel continuous casting is considered

Linear tetrahedral finite elements for thermal shock problems

International audiencePurpose - The paper seeks to present an original method for the numerical treatment of thermal shocks in non-linear heat transfer finite element analysis. Design/methodology/approach - The 3D finite element thermal analysis using linear standard tetrahedral elements may be affected by spurious local extrema in the regions affected by thermal shocks, in such a severe ways to directly discourage the use of these elements. This is especially true in the case of solidification problems, in which melted alloys at very high temperature contact low diffusive mould materials. The present work proposes a slight modification to the discrete heat equation in order to obtain a system matrix in M-matrix form, which ensures an oscillation-free solution. Findings - The proposed "diffusion-split" method consists basically of using a modified conductivity matrix. It allows for solutions based on linear tetrahedral elements. The performance of the method is evaluated by means of a test case with analytical solution, as well as an industrial application, for which a well-behaved numerical solution is available. Originality/value - The proposed method should be helpful for computational engineers and software developers in the field of heat transfer analysis. It can be implemented in most existing finite element codes with minimal effort

Computational design of thermo-mechanical metadevices using topology optimization

The present work has been conducted in order to introduce a novel approach for the design of mechanical devices conceived to manipulate the displacements field in linear elastic materials subjected to thermal gradients. Such an approach involves the solution of a topology optimization problem where the objective function defines the error in achieving a prescribed displacement field, and the mechanical device consists of two macroscopically distinguishable isotropic candidate materials. The material distribution is defined as a continuous function by following the solid isotropic microstructure (or material) with penalization (SIMP) method. The so-designed devices are easy to manufacture, since the design variables dictate the candidate materials distribution. Based on such an approach it is not necessary to devise further ways to simultaneously mimicking several thermal and mechanical effective properties, as required by coordinates transformation-based metamaterial design methods. Although the candidate materials are isotropic, the mechanical device behaves as a metamaterial allowing the desired manipulation of the displacements field. As an example, this topology optimization-based approach is applied to the design of an elastostatic cloaking device subjected to thermal gradients, considering also thermo-dependent mechanical properties.Fil: Álvarez Hostos, Juan Carlos. 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. Universidad Central de Venezuela; VenezuelaFil: Fachinotti, Victor Daniel. 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: Peralta, Ignacio. 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. Universidad Tecnológica Nacional; Argentin

Simulación Interactiva de Dinámica de Fluidos con Transferencia de Calor mediante Métodos de Partículas

Este trabajo esta dedicado a mostrar los primeros resultados obtenidos en simulación de problemas acoplados de dinámica de fluidos y de transferencia de calor usando el método de partículas llamado Smoothed Particle Hydrodynamics (SPH). La técnica empleada consiste en la solución simultánea de las ecuaciones de dinámica de fluidos y de transferencia de calor en formulación Lagrangiana usando una discretización tipo SPH desarrollada en el CIMEC. Aquí se presentan las primeras validaciones del modelo y los primeros ejemplos de su aplicación. En el trabajo se podrá apreciar como influyen los fenómenos de advección en la transferencia de calor. Los ejemplos y el código han sido corridos en una plataforma de simulación también desarrollada en el CIMEC. La plataforma permite cambiar interactivamente propiedades físicas del fluido, condiciones de contorno como el movimiento de paredes, o su temperatura, todo ello interactivamente mientras transcurre la simulación. El desarrollo permitirá la solución y simulación interactiva de variados e interesantes problemas de convección natural y transferencia de calor.Fil: Limache, Alejandro Cesar. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Santa Fe. Instituto de Desarrollo Tecnológico Para la Industria Química (i); ArgentinaFil: Rojas Fredini, Pablo Sebastián. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Santa Fe. Instituto de Desarrollo Tecnológico Para la Industria Química (i); ArgentinaFil: Fachinotti, Victor Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Santa Fe. Instituto de Desarrollo Tecnológico Para la Industria Química (i); Argentin

A 3D-fem model solving thermomechanics and macrosegregation in binary alloys solidification

International audienceThis paper introduces a three-dimensional numerical model for the coupled solution of momentum, energy and solute conservation equations, for binary alloys solidification. The spatial discretisation is carried out using linear tetrahedral finite elements, particularly those of P1+/P1 type for the velocity-pressure resolution of momentum equation. The liquid flow in the mushy zone is assumed to be governed by the Darcy's law. Thermal and buoyancy forces are taken into account by means of the Boussinesq's model. Microsegregation obeys the lever rule. The resulting solute transport equation is solved by the SUPG method. Coupling strategy between momentum, energy and solute equations is discussed and two applications are studied

Application of the arbitrary Eulerian Lagrangian finite element formulation to the thermomechanical simulation of casting processes, with focus on pipe shrinkage prediction

International audienceThe Arbitrary Lagrangian-Eulerian formulation (ALE) has become an indispensable component of finite element thermomechanical computations of casting processes. As it is an intermediate formulation between the Lagrangian formulation (material convected mesh) and the Eulerian one (fixed mesh), it allows the simultaneous computation of important phenomena: Deformation and stresses affecting solidified regions, yielding the computation of air gap evolution at part/mold interfaces. In such regions, the formulation is essentially Lagrangian. Thermosolutal convection flow in the non solidified regions; here the ALE formulation tends to a pure Eulerian one (stationary mesh). Free surface evolution at top of risers, leading to the prediction of pipe defects (macroshrinkage). In this case the ALE formulation allows the follow up of the free surface. After a brief reminder of the constitutive equations to be used in thermomechanical modeling of solidification, the mechanical equations are presented and their resolution in the context of FEM-ALE. We insist on the transport analysis, a key-point of ALE, and present a validation of the original scheme that is used here. Finally, we focus on the prediction of pipe shrinkage formation and show two industrial examples

Problema térmico inverso de conducción de calor aplicado al tratamiento térmico de austemperado usando optimización basada en simulaciones

La etapa de enfriamiento desde la temperatura de austenización hasta la temperatura ambiente es la etapa más crítica del tratamiento térmico de austemperado de una pieza de fundición nodular. La estimación tanto espacial como temporal del flujo de calor a través la superficie externa de la pieza durante dicha etapa, a fin de que la evolución de la temperatura en puntos interiores de la pieza siga unacurva temperatura-tiempo de austemperado prestablecida, constituye un problema térmico inverso de conducción de calor (IHCP, por sus siglas en inglés). En el presente trabajo, se propone la resolución de dicho problema inverso mediante el algoritmo de punto interior conocido como IPOPT. La función objetivo, que es el error en la consecución de la temperatura deseada a lo largo del tratamiento en una seriede puntos, requiere la solución de la ecuación de calor en estado transitorio durante todo el tratamiento, obtenida mediante el Método de Elementos Finitos (MEF). Las variables de diseño del problema definen la evolución de las condiciones de borde enfriamiento superficial a lo largo del tratamiento. En particular, se adoptó como variables de diseño la magnitud del flujo de calor impuesto en distintas porciones de la frontera, a distintos instantes de tiempo. De la aplicación del modelo a un caso de estudio representativo del austemperado de una pieza de fundición nodular, se demuestra que se puede puede determinar con precisión las condiciones de enfriamiento óptimas en ciertas porciones de la pieza, pero el resto de la pieza puede sufrir historias de temperaturas que perjudiquen sus propiedades mecánicas. A partir de esas observaciones, se determinan finalmente las limitaciones del tratamiento de austemperado de fundición nodular por enfriamiento superficial.Publicado en: Mecánica Computacional vol. XXXV no.30Facultad de Ingenierí

Optimization of Multilayered Walls for Building Envelopes Including PCM-Based Composites

This work proposes a numerical procedure to simulate and optimize the thermal response of a multilayered wallboard system for building envelopes, where each layer can be possibly made of Phase Change Materials (PCM)-based composites to take advantage of their Thermal-Energy Storage (TES) capacity. The simulation step consists in solving the transient heat conduction equation across the whole wallboard using the enthalpy-based finite element method. The weather is described in detail by the Typical Meteorological Year (TMY) of the building location. Taking the TMY as well as the wall azimuth as inputs, EnergyPlusTM is used to define the convective boundary conditions at the external surface of the wall. For each layer, the material is chosen from a predefined vade mecum, including several PCM-based composites developed at the Institut für Werkstoffe im Bauwesen of TU Darmstadt together with standard insulating materials (i.e., EPS or Rockwool). Finally, the optimization step consists in using genetic algorithms to determine the stacking sequence of materials across the wallboard to minimize the undesired heat loads. The current simulation-based optimization procedure is applied to the design of envelopes for minimal undesired heat losses and gains in two locations with considerably different weather conditions, viz. Sauce Viejo in Argentina and Frankfurt in Germany. In general, for each location and all the considered orientations (north, east, south and west), optimal results consist of EPS walls containing a thin layer made of the PCM-based composite with highest TES capacity, placed near the middle of the wall and closer to the internal surface