Smooth and Robust Solutions for Dirichlet Boundary Control of Fluid-Solid Conjugate Heat Transfer Problems

This work offers new computational methods for the optimal control of the conjugate heat transfer (CHT) problem in thermal science. Conjugate heat transfer has many important industrial applications, such as heat exchange processes in power plants and cooling in electronic packaging industry, and has been a staple of computational methods in thermal science for many years. This work considers the Dirichlet boundary control of fluid-solid CHT problems. The CHT system falls into the category of multi-physics problems. Its domain typically consists of two parts, namely, a solid region subject to thermal heating or cooling and a conjugate fluid region responsible for thermal convection transport. These two different physical systems are strongly coupled through the thermal boundary condition at the fluid-solid interface. The objective in the CHT boundary control problem is to select optimally the fluid inflow profile that minimizes an objective function that involves the sum of the mismatch between the temperature distribution in the system and a prescribed temperature profile and the cost of the control. This objective is realized by minimizing a nonlinear objective function of the boundary control and the fluid temperature variables, subject to partial differential equations (PDE) constraints governed by the coupled heat diffusion equation in the solid region and mass, momentum and energy conservation equations in the fluid region. Although CHT has received extensive attention as a forward problem, the optimal Dirichlet velocity boundary control for the coupled CHT process to our knowledge is only very sparsely studied analytically or computationally in the literature [131]. Therefore, in Part I, we describe the formulation of the optimal control problem and introduce the building blocks for the finite element modeling of the CHT problem, namely, the diffusion equation for the solid temperature, the convection-diffusion equation for the fluid temperature, the incompressible viscous Navier-Stokes equations for the fluid velocity and pressure, and the model verification of CHT simulations. In Part II, we provide theoretical analysis to explain the nonsmoothness issue which has been observed in this study and in Dirichlet boundary control of Navier-Stokes flows by other scientists. Based on these findings, we use either explicit or implicit numerical smoothing to resolve the nonsmoothness issue. Moreover, we use the numerical continuation on regularization parameters to alleviate the difficulty of locating the global minimum in one shot for highly nonlinear optimization problems even when the initial guess is far from optimal. Two suites of numerical experiments have been provided to demonstrate the feasibility, effectiveness and robustness of the optimization scheme. In Part III, we demonstrate the strategy of achieving parallel scalable algorithms for CHT models in Simulations of Reactor Thermal Hydraulics. Our motivation originates from the observation that parallel processing is necessary for optimal control problems of very large scale, when the simulation of the underlying physics (or PDE constraints) involves millions or billions of degrees of freedom. To achieve the overall scalability of optimal control problems governed by PDE constraints, scalable components that resolve the PDE constraints and their adjoints are the key. In this Part, first we provide the strategy of designing parallel scalable solvers for each building blocks of the CHT modeling, namely, for the discrete diffusive operator, the discrete convection-diffusion operator, and the discrete Navier-Stokes operator. Second, we demonstrate a pair of effective, robust, parallel, and scalable solvers built with collaborators for simulations of reactor thermal hydraulics. Finally, in the the section of future work, we outline the roadmap of parallel and scalable solutions for Dirichlet boundary control of fluid-solid conjugate heat transfer processes

Development of a conjugate heat transfer solver

The current research study presents a numerical approach in modelling the conjugate heat transfer system of the gas-turbine rotating discs-cavities. The work was undertaken to understand such phenomena and, more specifically, to numerically investigate the thermal interactions in rotating discs-cavities. The developed solver is capable of dealing with complex heat transfer problems, such as unsteady three-dimensional compressible rotating-flows. The development was based on integrating an inhouse computational fluid dynamics code (SURF) with a heat conduction solver internally. Method of interpolation using mapped area was also introduced for treating non-matching meshes at interface, which plays an effective role in exchanging boundary data. This thesis also documents the development of a numerical finite volume cell-vertex hybrid edgebased heat conduction code by the author using FORTRAN. The heat conduction solver was developed and validated to deal with three dimensional solid-domains using unstructured elements. The validation process was carried out on several test cases for investigating the temperature distribution. The test results were presented to show good agreement with the analytical, experimental and other commercial numerical solutions where they exist

Koopman analysis of the long-term evolution in a turbulent convection cell

We analyse the long-time evolution of the three-dimensional flow in a closed cubic turbulent Rayleigh-B\'{e}nard convection cell via a Koopman eigenfunction analysis. A data-driven basis derived from diffusion kernels known in machine learning is employed here to represent a regularized generator of the unitary Koopman group in the sense of a Galerkin approximation. The resulting Koopman eigenfunctions can be grouped into subsets in accordance with the discrete symmetries in a cubic box. In particular, a projection of the velocity field onto the first group of eigenfunctions reveals the four stable large-scale circulation (LSC) states in the convection cell. We recapture the preferential circulation rolls in diagonal corners and the short-term switching through roll states parallel to the side faces which have also been seen in other simulations and experiments. The diagonal macroscopic flow states can last as long as a thousand convective free-fall time units. In addition, we find that specific pairs of Koopman eigenfunctions in the secondary subset obey enhanced oscillatory fluctuations for particular stable diagonal states of the LSC. The corresponding velocity field structures, such as corner vortices and swirls in the midplane, are also discussed via spatiotemporal reconstructions.Comment: 32 pages, 9 figures, article in press at Journal of Fluid Mechanic

Fourier spectral methods for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is computationally demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reactiondiffusion equations. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is show-cased by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models,together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator

Steklov Spectral Geometry for Extrinsic Shape Analysis

We propose using the Dirichlet-to-Neumann operator as an extrinsic alternative to the Laplacian for spectral geometry processing and shape analysis. Intrinsic approaches, usually based on the Laplace-Beltrami operator, cannot capture the spatial embedding of a shape up to rigid motion, and many previous extrinsic methods lack theoretical justification. Instead, we consider the Steklov eigenvalue problem, computing the spectrum of the Dirichlet-to-Neumann operator of a surface bounding a volume. A remarkable property of this operator is that it completely encodes volumetric geometry. We use the boundary element method (BEM) to discretize the operator, accelerated by hierarchical numerical schemes and preconditioning; this pipeline allows us to solve eigenvalue and linear problems on large-scale meshes despite the density of the Dirichlet-to-Neumann discretization. We further demonstrate that our operators naturally fit into existing frameworks for geometry processing, making a shift from intrinsic to extrinsic geometry as simple as substituting the Laplace-Beltrami operator with the Dirichlet-to-Neumann operator.Comment: Additional experiments adde

Convection in the Melt

A physical problem involving the melting/freezing of a phase-change material (PCM) is the applied setting of this research. The development of models that couple the partial differential equations for energy transport and fluid motion with phases of differing densities is a primary goal of the research. In Chapter 2, a general framework is developed for the formulation of conservation laws that admit interfaces. A notion of weak solution is developed and its relation with classical and other weak formulations is discussed. Conditions that hold across various kinds of interfaces are also developed. The formulation is examined for the conservation of mass, momentum and energy in Chapter 3. In Chapter 4, a numerical method for the solution of conservation law equations is given. The method uses a Crank-Nicolson time discretization and solves the implicit equations with a Newton/Approximate Factorization technique. The method captures interfaces and is consistent with the control volume weak formulations of Chapter 2. The numerical solution converges to the distributional solution of the conservation law. In Chapter 5, three applications of the theory are developed and numerical computations are presented. First, a one dimensional problem is studied involving conservation of mass. momentum and energy in a phase-change material with a liquid density larger than that of the solid. The second application is a suction problem in two dimensions. The bulk movement of a liquid and void are simulated with and without the effects of surface tension. The third application is to a three-dimensional simulation of the heating of a cylindrical canister of PCM in 1-g and 0-g. For this simulation the Marangoni stress is the important driving force on the flow

Spectral methods for CFD

One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched