Transient dynamics and structure of optimal excitations in thermocapillary spreading: Precursor film model

Linearized modal stability theory has shown that the thermocapillary spreading of a liquid film on a homogeneous, completely wetting surface can produce a rivulet instability at the advancing front due to formation of a capillary ridge. Mechanisms that drain fluid from the ridge can stabilize the flow against rivulet formation. Numerical predictions from this analysis for the film speed, shape, and most unstable wavelength agree remarkably well with experimental measurements even though the linearized disturbance operator is non-normal, which allows transient growth of perturbations. Our previous studies using a more generalized nonmodal stability analysis for contact lines models describing partially wetting liquids (i.e., either boundary slip or van der Waals interactions) have shown that the transient amplification is not sufficient to affect the predictions of eigenvalue analysis. In this work we complete examination of the various contact line models by studying the influence of an infinite and flat precursor film, which is the most commonly employed contact line model for completely wetting films. The maximum amplification of arbitrary disturbances and the optimal initial excitations that elicit the maximum growth over a specified time, which quantify the sensitivity of the film to perturbations of different structure, are presented. While the modal results for the three different contact line models are essentially indistinguishable, the transient dynamics and maximum possible amplification differ, which suggests different transient dynamics for completely and partially wetting films. These differences are explained by the structure of the computed optimal excitations, which provides further basis for understanding the agreement between experiment and predictions of conventional modal analysis

Heat transfer and second law analyses of forced convection in a channel partially filled by porous media and featuring internal heat sources

This paper provides a comprehensive study on the heat transfer and entropy generation rates in a channel partially filled with a porous medium and under constant wall heat flux. The porous inserts are attached to the walls of the channel and the system features internal heat sources due to exothermic or endothermic physical or physicochemical processes. Darcy-Brinkman model is used for modelling the transport of momentum and an analytical study on the basis of LTNE (local thermal non-equilibrium) condition is conducted. Further analysis through considering the simplifying, LTE (local thermal equilibrium) model is also presented. Analytical solutions are, first, developed for the velocity and temperature fields. These are subsequently incorporated into the fundamental equations of entropy generation and both local and total entropy generation rates are investigated for a number of cases. It is argued that, comparing with LTE, the LTNE approach yields more accurate results on the temperature distribution within the system and therefore reveals more realistic Nusselt number and entropy generation rates. In keeping with the previous investigations, bifurcation phenomena are observed in the temperature field and rates of entropy generation. It is, further, demonstrated that partial filling of the channel leads to a substantial reduction of the total entropy generation. The results also show that the exothermicity or endothermicity characteristics of the system have significant impacts on the temperature fields, Nusselt number and entropy generation rates

Numerical and experimental investigation of a vertical LiBr falling film absorber considering wave regimes and in presence of mist flow

The absorber represents the most critical component in absorption systems and one of the key issues. In this component complex heat and mass transfer phenomena during the absorption process takes place simultaneously. For this reason the development of mathematical models validated against experimental data always constitutes useful tools for the design and improvement of falling film absorbers. A testing device has been designed and built to reproduce absorption phenomena in vertical LiBr-H2O falling film absorbers with the primary objective to obtain experimental data. On the other hand, a mathematical model of falling film absorption of H2O vapour in LiBr aqueous solutions has been implemented. Wave regime is considered by including and solving the Free Surface Deflection Equation. The numerical results are validated using the experimental data. During the development of this work, the authors have paid careful attention to the verification of experimental data. Such verification consists of performing energy and mass balances in the fluid film side. Important discrepancies were found in our experimental data. Therefore, an extensive study was carried out in order to find the source of such errors. The conclusion is that there is a drag of LiBr solution in the water vapour which increases with the Re number. This mist flow cannot be measured experimentally, but can be evaluated in an indirect way. The mathematical models have been adapted in order to consider the influence of mist flow. On the other hand, in the literature there are not many experimental works related to falling film absorbers which expose enough information to verify the reliability of their experimental data.Peer ReviewedPostprint (author's final draft

New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid

Seismic imaging of the mantle has revealed large and small scale heterogeneities in the lower mantle; specifically structures known as large low shear velocity provinces (LLSVP) below Africa and the South Pacific. Most interpretations propose that the heterogeneities are compositional in nature, differing in composition from the overlying mantle, an interpretation that would be consistent with chemical geodynamic models. Numerical modeling of persistent compositional interfaces presents challenges, even to state-of-the-art numerical methodology. For example, some numerical algorithms for advecting the compositional interface cannot maintain a sharp compositional boundary as the fluid migrates and distorts with time dependent fingering due to the numerical diffusion that has been added in order to maintain the upper and lower bounds on the composition variable and the stability of the advection method. In this work we present two new algorithms for maintaining a sharper computational boundary than the advection methods that are currently openly available to the computational mantle convection community; namely, a Discontinuous Galerkin method with a Bound Preserving limiter and a Volume-of-Fluid interface tracking algorithm. We compare these two new methods with two approaches commonly used for modeling the advection of two distinct, thermally driven, compositional fields in mantle convection problems; namely, an approach based on a high-order accurate finite element method advection algorithm that employs an artificial viscosity technique to maintain the upper and lower bounds on the composition variable as well as the stability of the advection algorithm and the advection of particles that carry a scalar quantity representing the location of each compositional field. All four of these algorithms are implemented in the open source FEM code ASPECT

Optimal boundary control with critical penalization for a PDE model of fluid-solid interactions

We study the finite-horizon optimal control problem with quadratic functionals for an established fluid-structure interaction model. The coupled PDE system under investigation comprises a parabolic (the fluid) and a hyperbolic (the solid) dynamics; the coupling occurs at the interface between the regions occupied by the fluid and the solid. We establish several trace regularity results for the fluid component of the system, which are then applied to show well-posedness of the Differential Riccati Equations arising in the optimization problem. This yields the feedback synthesis of the unique optimal control, under a very weak constraint on the observation operator; in particular, the present analysis allows general functionals, such as the integral of the natural energy of the physical system. Furthermore, this work confirms that the theory developed in Acquistapace et al. [Adv. Differential Equations, 2005] -- crucially utilized here -- encompasses widely differing PDE problems, from thermoelastic systems to models of acoustic-structure and, now, fluid-structure interactions.Comment: 22 pages, submitted; v2: misprints corrected, a remark added in section

A Method for Geometry Optimization in a Simple Model of Two-Dimensional Heat Transfer

This investigation is motivated by the problem of optimal design of cooling elements in modern battery systems. We consider a simple model of two-dimensional steady-state heat conduction described by elliptic partial differential equations and involving a one-dimensional cooling element represented by a contour on which interface boundary conditions are specified. The problem consists in finding an optimal shape of the cooling element which will ensure that the solution in a given region is close (in the least squares sense) to some prescribed target distribution. We formulate this problem as PDE-constrained optimization and the locally optimal contour shapes are found using a gradient-based descent algorithm in which the Sobolev shape gradients are obtained using methods of the shape-differential calculus. The main novelty of this work is an accurate and efficient approach to the evaluation of the shape gradients based on a boundary-integral formulation which exploits certain analytical properties of the solution and does not require grids adapted to the contour. This approach is thoroughly validated and optimization results obtained in different test problems exhibit nontrivial shapes of the computed optimal contours.Comment: Accepted for publication in "SIAM Journal on Scientific Computing" (31 pages, 9 figures

Thermophysical Phenomena in Metal Additive Manufacturing by Selective Laser Melting: Fundamentals, Modeling, Simulation and Experimentation

Among the many additive manufacturing (AM) processes for metallic materials, selective laser melting (SLM) is arguably the most versatile in terms of its potential to realize complex geometries along with tailored microstructure. However, the complexity of the SLM process, and the need for predictive relation of powder and process parameters to the part properties, demands further development of computational and experimental methods. This review addresses the fundamental physical phenomena of SLM, with a special emphasis on the associated thermal behavior. Simulation and experimental methods are discussed according to three primary categories. First, macroscopic approaches aim to answer questions at the component level and consider for example the determination of residual stresses or dimensional distortion effects prevalent in SLM. Second, mesoscopic approaches focus on the detection of defects such as excessive surface roughness, residual porosity or inclusions that occur at the mesoscopic length scale of individual powder particles. Third, microscopic approaches investigate the metallurgical microstructure evolution resulting from the high temperature gradients and extreme heating and cooling rates induced by the SLM process. Consideration of physical phenomena on all of these three length scales is mandatory to establish the understanding needed to realize high part quality in many applications, and to fully exploit the potential of SLM and related metal AM processes

Rayleigh-B\'enard convection with a melting boundary

We study the evolution of a melting front between the solid and liquid phases of a pure incompressible material where fluid motions are driven by unstable temperature gradients. In a plane layer geometry, this can be seen as classical Rayleigh-B\'enard convection where the upper solid boundary is allowed to melt due to the heat flux brought by the fluid underneath. This free-boundary problem is studied numerically in two dimensions using a phase-field approach, classically used to study the melting and solidification of alloys, which we dynamically couple with the Navier-Stokes equations in the Boussinesq approximation. The advantage of this approach is that it requires only moderate modifications of classical numerical methods. We focus on the case where the solid is initially nearly isothermal, so that the evolution of the topography is related to the inhomogeneous heat flux from thermal convection, and does not depend on the conduction problem in the solid. From a very thin stable layer of fluid, convection cells appears as the depth -- and therefore the effective Rayleigh number of the layer increases. The continuous melting of the solid leads to dynamical transitions between different convection cell sizes and topography amplitudes. The Nusselt number can be larger than its value for a planar upper boundary, due to the feedback of the topography on the flow, which can stabilize large-scale laminar convection cells.Comment: 36 pages, 16 figure