Transient heat conduction problems using B.I.E.M

A method, using boundary elements, is presented as a solution to plane transient heat conduction. The proposed method considers the governing equation to be a Helmholtz's equation and solves the problem of time variation using step by step integration. A numerical procedure is developed and its effectiveness verified. Several examples are provided and their results compared with the theoretical ones.

FTCS finite difference scheme GPGPU parallel computing for the heat conduction equation = Programación en paralelo GPGPU del método en diferencias finitas FTCS para la ecuación del calor

En el presente artículo se muestran las ventajas de la programación en paralelo resolviendo numéricamente la ecuación del calor en dos dimensiones a través del método de diferencias finitas explícito centrado en el espacio FTCS. De las conclusiones de este trabajo se pone de manifiesto la importancia de la programación en paralelo para tratar problemas grandes, en los que se requiere un elevado número de cálculos, para los cuales la programación secuencial resulta impracticable por el elevado tiempo de ejecución. En la primera sección se describe brevemente los conceptos básicos de programación en paralelo. Seguidamente se resume el método de diferencias finitas explícito centrado en el espacio FTCS aplicado a la ecuación parabólica del calor. Seguidamente se describe el problema de condiciones de contorno y valores iniciales específico al que se va a aplicar el método de diferencias finitas FTCS, proporcionando pseudocódigos de una implementación secuencial y dos implementaciones en paralelo. Finalmente tras la discusión de los resultados se presentan algunas conclusiones. In this paper the advantages of parallel computing are shown by solving the heat conduction equation in two dimensions with the forward in time central in space (FTCS) finite difference method. Two different levels of parallelization are consider and compared with traditional serial procedures. We show in this work the importance of parallel computing when dealing with large problems that are impractical or impossible to solve them with a serial computing procedure. In the first section a summary of parallel computing approach is presented. Subsequently, the forward in time central in space (FTCS) finite difference method for the heat conduction equation is outline, describing how the heat flow equation is derived in two dimensions and the particularities of the finite difference numerical technique considered. Then, a specific initial boundary value problem is solved by the FTCS finite difference method and serial and parallel pseudo codes are provided. Finally after results are discussed some conclusions are presented.

A new boundary condition solved with B.I.E.M.

Among the classical operators of mathematical physics the Laplacian plays an important role due to the number of different situations that can be modelled by it. Because of this a great effort has been made by mathematicians as well as by engineers to master its properties till the point that nearly everything has been said about them from a qualitative viewpoint. Quantitative results have also been obtained through the use of the new numerical techniques sustained by the computer. Finite element methods and boundary techniques have been successfully applied to engineering problems as can be seen in the technical literature (for instance [ l ] , [2], [3] . Boundary techniques are especially advantageous in those cases in which the main interest is concentrated on what is happening at the boundary. This situation is very usual in potential problems due to the properties of harmonic functions. In this paper we intend to show how a boundary condition different from the classical, but physically sound, is introduced without any violence in the discretization frame of the Boundary Integral Equation Method. The idea will be developed in the context of heat conduction in axisymmetric problems but it is hoped that its extension to other situations is straightforward. After the presentation of the method several examples will show the capabilities of modelling a physical problem.

Simulation of Temperatures and Stresses during Construction of a RCC Dam

Se describen las variaciones de temperaturas y de tensiones durante la construcción de presas de hormigón compactado. The curing of concrete is an exothermic process. The heat of hydration generated induces temperature increases in the concrete, which will disappear in the long term by heat conduction in the concrete mass and thermal exchanges with the environment. The problem is of particularly interest for large concrete masses, as is the case of dams, because the time involved in the heat diffusion process grows with the square of the dimensions and a hotter dam interior implies the possibility of cracking the exposed surfaces of the dam. The Cuira dam, currently being built in Venezuela using roller compacted concrete, is a 134 m high, arch-gravity dam. In support of the design, different strategies were analysed, including various combinations of cooling of the water and the aggregate in order to achieve acceptable results. The calculations were conducted with Abaqus, taking into account all the necessary mechanical and thermal characteristics, as well as the relevant non-linearities. The analyses led to the conclusion that no cooling was required, even taking into account the stress state imposed by an early and rapid filling of the reservoir.

Self-similar motion of laser half-space plasmas. I. Deflagration regime

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density 0,produced by the (anomalous) absorption of a laser pulse of irradiation = (j>0f/T(0< (< T) at the critical density nc(«c/«0=eheat conduction and ion-electron energy exchange and retains three dimensionless numbers: e, Zt (ion charge number), and a = (9/c/4m,) (T/C 2n l/4>oKe)213, where k, m, are Boltzmann's constant and the ion mass, and Ke X (electron temperature)5'2 = heat conductivity. If a >e- 4 ' 3 , a deflagration wave separates an isentropic compression with a shock bounding the undisturbed plasma, and an isentropic expansion flow to the vacuum. The structures of these three regions are completely determined; in particular, the Chapman-Jouguet condition is proved and the density behind the deflagration is found. The deflagration-compression thickness ratio is large (small) for a^e- 5 ' 3(a>e- 5 ' 3 ) . The compression to expansion ratio for both energy and thickness is 0(e"2). For Z,- large, a deflagration exists even if a~e~413. Condition a>e~4'3 may be applied to pulses that are not linear.

Self-similar motion of laser half-space plasmas. II. Thermal waves and intermediate regimes

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density n,produced by the (anomalous) absorption of a laser pulse of irradiation

heat conduction and ion-electron energy exchange, involves three dimensionless numbers: e = nc/n0 assumed small, Z, (ion charge number), and a parameter a€~4'3, a qualitative discussion of how plasma behavior changes with a, is given.

Transition from isentropic to isothermal expansion in laser produced plasma

The transition that the expansion flow of laser-produced plasmas experiences when one moves from long, low intensity pulses (temperature vanishing at the isentropic plasma-vacuum front,lying at finite distance) to short, intense ones (non-zero, uniform temperature at the plasma-vacuum front, lying at infinity) is studied. For plznar geometry and lqge ion number Z, the transition occurs for dq5/dt=0.14(27/8)k712Z’1zn$/m4f, 12nK,,; mi, and K are laser intensity, critical density,ion mass, and Spitzer’s heat conduction coefficient. This result remains valid for finite Zit,h ough the numerical factor in d$/dt is different. Shorter wavelength lasers and higher 4 plasmas allow faster rising pulses below transition.

Electron cooling in probe collection from magnetized plasmas with anomalous transport

The electron-retarding range of the current-voltage characteristic of a flat Langmuir probe perpendicular to a strong magnetic field in a fully ionized plasma is analysed allowing for anomalous (Bohm) cross-field transport and temperature changes in the collection process. With probe size and ion thermal gyroradius comparable, and smaller than the electron mean free path, there is an outer quasineutral region with ion viscosity determinant in allowing nonambipolar parallel and cross flow. A potential overshoot lying either at the base or inside the quasineutral region both makes ions follow Boltzmann's law at negative bias and extends the electron-retarding range to probe bias e(j)p ~ +2Too. Electron heating and cooling occur roughly at positive and negative bias, with a re-minimum around efa ~ - 2 7 ^ ; far from the probe heat conduction cools and heats electrons at and radially away from the probe axis, respectively. The potential overshoot with no thermal effects would reduce the electron current Ie, making the In Ie versus 4>p graph downwards-concave,but cooling further reduces Ie substantially, and may tilt the slope upwards past the temperature minimum. The domain of strict validity of our analysis is narrow in case of low ion mass (deuterium), breaking down with the ion Boltzmann law.

Mejora y optimización de un código de simulación hidrodinámico con transporte de radiación para el estudio de sistemas en condiciones de alta densidad de energía

Se presentan las mejoras introducidas en un código de transporte de radiación acoplada a la hidrodinámica llamado ARWEN para el estudio de sistemas en el rango de física de alta densidad de energía (High Energy Density Physics). Los desarrollos introducidos se basan en las siguientes áreas: ít>,~ Ecuaciones de estado: se desarrolla una nueva metodología mediante la cual es posible ajustar los resultados de un modelo simple de ecuaciones de estado como QEOS a datos experimentales y resultados de AIMD. Esta metodología tiene carácter general para poder ser aplicada en gran cantidad de materuales de interés y amplia la flexibilidad de ajuste de los métodos de los que ha partido como base este trabajo. En segundo lugar, se ha desarrollado una librería para la gestión de tablas de datos de ecuaciones de estado que también incluye la gestión de tablas con datos de opacidades y de ionización. Esta nueva librería extiende las capacidades de la anterior al tener llamadas más específicas que aceleran los cálculos, y posibilidad de uso de varias tablas a la vez. Solver de difusión: se ha desarrollado un nuevo paquete para resolver la ecuación de difusión que se aplicará a la conducción de calor dentro del plasma. El método anterior no podía ser ejecutado en paralelo y producía resultados dependientes de la resolución de la malla, mientras que este método es paralelizable y además obtiene una solución con mejor convergencia, lo que supone una solución que no depende del refinamiento del mallado. Revisión del paquete de radiación: en primer lugar se ha realizado una revisión de la implementación del modelo de radiación descubriendo varios errores que han sido depurados. También se ha incluido la nueva librería de gestión de tablas de opacidades que permiten la obtención de las propiedades ópticas del plasma en multigrupos de energía. Por otra parte se ha extendido el cálculo de los coeficientes de transporte al esquema multimaterial que ha introducido David Portillo García en el paquete hidrodinámico del código de simulación. Por último se ha revisado el esquema de resolución del transporte y se ha modificado para hacerlo paralelizable. • Se ha implementado un paquete de trazado de rayos para deposición láser que extiende la utilidad del anterior al ser en 3D y poder utilizar entonces diferentes configuraciones. • Una vez realizadas todas estas tareas se ha aplicado el código ARWEN al estudio de la astrofísica de laboratorio simulando los experimentos llevados a cabo en la instalación PALS por Chantal Stehlé acerca de ondas de choque radiativas. Se han comparado los resultados experimentales frente a las predicciones del código ARWEN obteniéndose una gran concordancia en la velocidad de la onda de choque generada y en las dimensiones del precursor. El código de simulación sobre el que se ha trabajado, junto con los desarrollos aportados por otros investigadores durante la realización de esta tesis, ha permitido participar en colaboraciones con laboratorios de Francia o Japón y se han producido resultados científicos publicados basados en el trabajo descrito en esta tesis. ABSTRACT Improvements in radiation hydrodynamic code ARWEN for the study of systems in the range of physics high energy density (High Energy Density Physics) are presented. The developments introduced are based on the following áreas: • Equations of state: a new methodology was developed to adjust the results of a simple Equation of State model like QEOS to experimental data and results of AIMD. This methodology can be applied to a large amount of materials and it increases the flexibility and range of the previous methods used as basis for this work. Also a new computer library has been developed to manage data tables of thermodynamic properties as well as includes the management of opacity and ionization data tables. This new library extends the capabilities of the previous one with more specific routines, and the possibility of using múltiple tables for several materials. • Diffusion solver: a new package has been developed to solve the diffusion equation applied to the heat conduction of the plasma. The previous method is not parallelizable and it produced mesh dependent results, while this new package can be executed in parallel and achieves a more converged solution that does not depend on the refinement of the mesh. • Radiation package: the check of the radiation model rose several bugs in the implementation that had been removed. The new computer library for EOS managing includes capabilities to store opacity tables for multigroups of energy. Moreover the transport coefficients calculations have been extended for the new multimaterial hydrodynamic package developed by David Portillo García. Also the solving methodology for the transport equation has been modified to make the code run in parallel. • A new ray tracing package has been introduced to extend the previous one to 3D. Once all these tasks has been implemented, the ARWEN code has been applied to study laboratory astrophysics systems. Simulations have been done in order to reproduce the results of the experiments carried out in PALS facility by Chantal Stehlé in radiative shock production. Experimental results are in cióse agreement to the ARWEN estimations of the speed of the shock wave and the length of the precursor. The simulation code used in this thesis, including the work done in ARWEN by other colleagues at the time of this research, allowed the collaboration with other research institution in France and Japan and some of the results presented in this thesis have been published in scientific journals.

Sheath vaporization of a monodisperse fuel-spray jet

The group vaporization of a monodisperse fuel-spray jet discharging into a hot coflowing gaseous stream is investigated for steady flow by numerical and asymptotic methods with a two-continua formulation used for the description of the gas and liquid phases. The jet is assumed to be slender and laminar, as occurs when the Reynolds number is moderately large, so that the boundary-layer form of the conservation equations can be employed in the analysis. Two dimensionless parameters are found to control the flow structure, namely the spray dilution parameter 1, defined as the mass of liquid fuel per unit mass of gas in the spray stream, and the group vaporization parameter e, defined as the ratio of the characteristic time of spray evolution due to droplet vaporization to the characteristic diffusion time across the jet. It is observed that, for the small values of e often encountered in applications, vaporization occurs only in a thin layer separating the spray from the outer droplet-free stream. This regime of sheath vaporization, which is controlled by heat conduction, is amenable to a simplified asymptotic description, independent of ε,in which the location of the vaporization layer is determined numerically as a free boundary in a parabolic problem involving matching of the separate solutions in the external streams, with appropriate jump conditions obtained from analysis of the quasi-steady vaporization front. Separate consideration of dilute and dense sprays, corresponding, respectively, to the asymptotic limits λ<<1 and λ>>1, enables simplified descriptions to be obtained for the different flow variables, including explicit analytic expressions for the spray penetration distance.

Desarrollo de un método numérico para flujos multimaterial en 2D en ARWEN

En esta tesis se presenta un método numérico para resolver las ecuaciones de Euler para flujos multimaterial en malla euleriana. Este solver se ha acoplado en el código hidrodinámico en dos dimensiones con transporte de radiación desarrollado en el Instituto de Fusión Nuclear de la UPM bajo la dirección del profesor Pedro Velarde, ARWEN. Los objetivos de este trabajo son: Desarrollo e implementación de un método de Godunov unsplit de alto orden multimaterial en 2D para malla euleriana en geometría cartesiana y geometría cilíndrica. Se presenta una extensión del trabajo realizado por Miller y Puckett (36) a una formulación unsplit. Además, se ha prestado especial atención al acoplamiento con el transporte de radiación y la conducción de calor. El método presentado se ha probado en una gran cantidad de problemas. Aplicación del código multimaterial al estudio de experimentos reales: • Simulación de una propuesta de experimento de laboratorio para reproducir la etapa de arrancamiento de material de la interacción entre el gas proveniente de la explosión de una supernova y la estrella secundaria en un escenario degenarado (SD). • Formación de jets en el laboratorio producidos por la colisión de dos plasmas. ABSTRACT We present a solver for the Euler equations for multimaterial flows in eulerian mesh. This solver has been coupled in the 2D AMR radiation transport code developed at Instituto de Fusión Nuclear (UPM) under the direction of professor Pedro Velarde, ARWEN. The main goals of this thesis are: Development and implementation of an 2D unsplit high-order Godunov method for multimaterial flows in eulerian mesh for cartesian and axialsimetry geometry. We present an extension of the work of Miller and Puckett (36) to an unsplit formulation. Also, we have paid special attention to the coupling with radiation transport and heat conduction. The method has been tested in a wide variety of problems. Application of the multimaterial solver to the study of real experiments: • Simulation of a proposal of a laboratory experiment aimed to reproducing the stripping stage of the interaction between the gas ejected during a supernova explosion and the secondary star in the Single Degenerate scenario. • Experiments of plasma jets in the laboratory obtained by the collission of two hot plasmas.

Effects of radiation losses on very lean methane/air flames propagating upward in a vertical tube

The stationary upward propagation of a very lean methane/air flame in a long vertical tube open at the bottom and closed at the top is simulated numerically using a single overall chemical reaction to model combustion and assuming an optically thin gas and a transparent or non-reflecting tube wall to approximately account for radiation losses from CO2CO2 and H2OH2O. Buoyancy plays a dominant role in the propagation of these flames and causes a large region of low velocity of the burnt gas relative to the flame to appear below the flame front when the equivalence ratio is decreased. The size of this region scales with the radius of the tube, and its presence enhances the effect of radiation losses, which would be otherwise negligible for a standard flammability tube, given the small concentration of radiating species. Heat conduction is found to be important in the low velocity region and to lead to a conduction flux from the flame to the burnt gas that causes extinction at the flame tip for a value of the equivalence ratio near the flammability limit experimentally measured in the standard tube. The effect of radiation losses decreases with the radius of the tube. Numerical results and order-of-magnitude estimates show that, in the absence of radiation, a very lean flame front fails to propagate only after recirculation of the burnt gas extends to its reaction region and drastically changes its structure. This condition is not realized for the standard flammability tube, but it seems to account for the flammability limit measured in a tube of about half the radius of the standard tube.

Nonlocal electron heat-flux

Electron thermal conduction in a not quite collisional unmagnetlzed plasma is analysed. The failure of classical results for temperature scale-length up to 100 times larger than thermal mean-free-path for electron scattering, and large ion-charge number Z , is discussed. Recent results from a nonlocal model of conduction at large Z are reviewed. Closed form expressions for Braginskii's coefficients a ,/3 , y for Z =0(1) are derived. An extension of the nonlocal model for Z =0(1) is discussed.

Laminar counterflow parallel-plate heat exchangers: Exact and approximate solutions

Multilayered, counterflow, parallel-plate heat exchangers are analyzed numerically and theoretically. The analysis, carried out for constant property fluids, considers a hydrodynamically developed laminar flow and neglects longitudinal conduction both in the fluid and in the plates. The solution for the temperature field involves eigenfunction expansions that can be solved in terms of Whittaker functions using standard symbolic algebra packages, leading to analytical expressions that provide the eigenvalues numerically. It is seen that the approximate solution obtained by retaining the first two modes in the eigenfunction expansion provides an accurate representation for the temperature away from the entrance regions, specially for long heat exchangers, thereby enabling simplified expressions for the wall and bulk temperatures, local heat-transfer rate, overall heat-transfer coefficient, and outlet bulk temperatures. The agreement between the numerical and theoretical results suggests the possibility of using the analytical solutions presented herein as benchmark problems for computational heat-transfer codes.