OUTPU T 29 OUTPU T 30 OUTPU T 31 OUTPU T 32 OUTPU T 25 OUTPU T 26 OUTPU T 27 OUTPU T 28 OUTPU T 21 OUTPU T 22 OUTPU T 23 OUTPU T 24 OUTPU T 17 OUTPU T 18 OUTPU T 19 OUTPU T 20 OUTPU T 13 OUTPU T 14 OUTPU T 15 OUTPU T 16 OUTPU T 9 OUTPU T 10 OUTPU T 11 OUTPU T 12 OUTPU T 5 OUTPU T 6 OUTPU T 7 OUTPU T 8 OUTPU T 1 OUTPU T 2 OUTPU T 3 OUTPU T 4 29 30 31 32 25 26 27 28 21 22 23 24 17 18 19 20 13 14 15 16 9

A User’s Guide

The phase-field method has become an important and extremely versatile technique for simulating microstructure evolution at the mesoscale. Thanks to the diffuse-interface approach, it allows us to study the evolution of arbitrary complex grain morphologies without any presumption on their shape or mutual distribution. It is also straightforward to account for different thermodynamic driving forces for microstructure evolution, such as bulk and interfacial energy, elastic energy and electric or magnetic energy, and the effect of different transport processes, such as mass diffusion, heat conduction and convection. The purpose of the paper is to give an introduction to the phase-field modeling technique. The concept of diffuse interfaces, the phase-field variables, the thermodynamic driving force for microstructure evolution and the kinetic phase-field equations are introduced. Furthermore, common techniques for parameter determination and numerical solution of the equations are discussed. To show the variety in phase-field models, different model formulations are exploited, depending on which is most common or most illustrative.

/pdf/an-introduction-to-phase-field-modeling-of-microstructure-2cm9coegv9.pdf

An introduction to phase-field modeling of microstructure evolution

The growth of a vapor bubble in a superheated liquid is controlled by three factors: the inertia of the liquid, the surface tension, and the vapor pressure. As the bubble grows, evaporation takes place at the bubble boundary, and the temperature and vapor pressure in the bubble are thereby decreased. The heat inflow requirement of evaporation, however, depends on the rate of bubble growth, so that the dynamic problem is linked with a heat diffusion problem. Since the heat diffusion problem has been solved, a quantitative formulation of the dynamic problem can be given. A solution for the radius of the vapor bubble as a function of time is obtained which is valid for sufficiently large radius. This asymptotic solution covers the range of physical interest since the radius at which it becomes valid is near the lower limit of experimental observation. It shows the strong effect of heat diffusion on the rate of bubble growth. Comparison of the predicted radius‐time behavior is made with experimental observations in superheated water, and very good agreement is found.

The growth of vapor bubbles in superheated liquids. report no. 26-6

Preface. 1 Introduction. 1.1 Importance of Heat Transfer. 1.2 Heat Transfer Modes. 1.3 The Laws of Heat Transfer. 1.4 Formulation of Heat Transfer Problems. 1.4.1 Heat transfer from a plate exposed to solar heat flux. 1.4.2 Incandescent lamp. 1.4.3 Systems with a relative motion and internal heat generation. 1.5 Heat Conduction Equation. 1.6 Boundary and Initial Conditions. 1.7 Solution Methodology. 1.8 Summary. 1.9 Exercise. Bibliography. 2 Some Basic Discrete Systems. 2.1 Introduction. 2.2 Steady State Problems. 2.2.1 Heat flow in a composite slab. 2.2.2 Fluid flow network. 2.2.3 Heat transfer in heat sinks (combined conduction-convection). 2.2.4 Analysis of a heat exchanger. 2.3 Transient Heat Transfer Problem (Propagation Problem). 2.4 Summary. 2.5 Exercise. Bibliography. 3 The Finite Elemen t Method. 3.1 Introduction. 3.2 Elements and Shape Functions. 3.2.1 One-dimensional linear element. 3.2.2 One-dimensional quadratic element. 3.2.3 Two-dimensional linear triangular elements. 3.2.4 Area coordinates. 3.2.5 Quadratic triangular elements. 3.2.6 Two-dimensional quadrilateral elements. 3.2.7 Isoparametric elements. 3.2.8 Three-dimensional elements. 3.3 Formulation (Element Characteristics). 3.3.1 Ritz method (Heat balance integral method-Goodman's method). 3.3.2 Rayleigh-Ritz method (Variational method). 3.3.3 The method of weighted residuals. 3.3.4 Galerkin finite element method. 3.4 Formulation for the Heat Conduction Equation. 3.4.1 Variational approach. 3.4.2 The Galerkin method. 3.5 Requirements for Interpolation Functions. 3.6 Summary. 3.7 Exercise. Bibliography. 4 Steady State Heat Conduction in One Dimension. 4.1 Introduction. 4.2 Plane Walls. 4.2.1 Homogeneous wall. 4.2.2 Composite wall. 4.2.3 Finite element discretization. 4.2.4 Wall with varying cross-sectional area. 4.2.5 Plane wall with a heat source: solution by linear elements. 4.2.6 Plane wall with a heat source: solution by quadratic elements. 4.2.7 Plane wall with a heat source: solution by modified quadratic equations (static condensation). 4.3 Radial Heat Flow in a Cylinder. 4.3.1 Cylinder with heat source. 4.4 Conduction-Convection Systems. 4.5 Summary. 4.6 Exercise. Bibliography. 5 Steady State Heat Conduction in Multi-dimensions. 5.1 Introduction. 5.2 Two-dimensional Plane Problems. 5.2.1 Triangular elements. 5.3 Rectangular Elements. 5.4 Plate with Variable Thickness. 5.5 Three-dimensional Problems. 5.6 Axisymmetric Problems. 5.6.1 Galerkin's method for linear triangular axisymmetric elements. 5.7 Summary. 5.8 Exercise. Bibliography. 6 Transient Heat Conduction Analysis. 6.1 Introduction. 6.2 Lumped Heat Capacity System. 6.3 Numerical Solution. 6.3.1 Transient governing equations and boundary and initial conditions. 6.3.2 The Galerkin method. 6.4 One-dimensional Transient State Problem. 6.4.1 Time discretization using the Finite Difference Method (FDM). 6.4.2 Time discretization using the Finite Element Method (FEM). 6.5 Stability. 6.6 Multi-dimensional Transient Heat Conduction. 6.7 Phase Change Problems-Solidification and Melting. 6.7.1 The governing equations. 6.7.2 Enthalpy formulation. 6.8 Inverse Heat Conduction Problems. 6.8.1 One-dimensional heat conduction. 6.9 Summary. 6.10 Exercise. Bibliography. 7 Convection Heat Transfer 173 7.1 Introduction. 7.1.1 Types of fluid-motion-assisted heat transport. 7.2 Navier-Stokes Equations. 7.2.1 Conservation of mass or continuity equation. 7.2.2 Conservation of momentum. 7.2.3 Energy equation. 7.3 Non-dimensional Form of the Governing Equations. 7.3.1 Forced convection. 7.3.2 Natural convection (Buoyancy-driven convection). 7.3.3 Mixed convection. 7.4 The Transient Convection-diffusion Problem. 7.4.1 Finite element solution to convection-diffusion equation. 7.4.2 Extension to multi-dimensions. 7.5 Stability Conditions. 7.6 Characteristic-based Split (CBS) Scheme. 7.6.1 Spatial discretization. 7.6.2 Time-step calculation. 7.6.3 Boundary and initial conditions. 7.6.4 Steady and transient solution methods. 7.7 Artificial Compressibility Scheme. 7.8 Nusselt Number, Drag and Stream Function. 7.8.1 Nusselt number. 7.8.2 Drag calculation. 7.8.3 Stream function. 7.9 Mesh Convergence. 7.10 Laminar Isothermal Flow. 7.10.1 Geometry, boundary and initial conditions. 7.10.2 Solution. 7.11 Laminar Non-isothermal Flow. 7.11.1 Forced convection heat transfer. 7.11.2 Buoyancy-driven convection heat transfer. 7.11.3 Mixed convection heat transfer. 7.12 Introduction to Turbulent Flow. 7.12.1 Solution procedure and result. 7.13 Extension to Axisymmetric Problems. 7.14 Summary. 7.15 Exercise. Bibliography. 8 Convection in Porous Media. 8.1 Introduction. 8.2 Generalized Porous Medium Flow Approach. 8.2.1 Non-dimensional scales. 8.2.2 Limiting cases. 8.3 Discretization Procedure. 8.3.1 Temporal discretization. 8.3.2 Spatial discretization. 8.3.3 Semi- and quasi-implicit forms. 8.4 Non-isothermal Flows. 8.5 Forced Convection. 8.6 Natural Convection. 8.6.1 Constant porosity medium. 8.7 Summary. 8.8 Exercise. Bibliography. 9 Some Examples of Fluid Flow and Heat Transfer Problems. 9.1 Introduction. 9.2 Isothermal Flow Problems. 9.2.1 Steady state problems. 9.2.2 Transient flow. 9.3 Non-isothermal Benchmark Flow Problem. 9.3.1 Backward-facing step. 9.4 Thermal Conduction in an Electronic Package. 9.5 Forced Convection Heat Transfer From Heat Sources. 9.6 Summary. 9.7 Exercise. Bibliography. 10 Implementation of Computer Code. 10.1 Introduction. 10.2 Preprocessing. 10.2.1 Mesh generation. 10.2.2 Linear triangular element data. 10.2.3 Element size calculation. 10.2.4 Shape functions and their derivatives. 10.2.5 Boundary normal calculation. 10.2.6 Mass matrix and mass lumping. 10.2.7 Implicit pressure or heat conduction matrix. 10.3 Main Unit. 10.3.1 Time-step calculation. 10.3.2 Element loop and assembly. 10.3.3 Updating solution. 10.3.4 Boundary conditions. 10.3.5 Monitoring steady state. 10.4 Postprocessing. 10.4.1 Interpolation of data. 10.5 Summary. Bibliography. A Green's Lemma. B Integration Formulae. B.1 Linear Triangles. B.2 Linear Tetrahedron. C Finite Element Assembly Procedure. D Simplified Form of the Navier-Stokes Equations. Index.

Fundamentals of the Finite Element Method for Heat and Fluid Flow

The fifth edition of "Numerical Methods for Engineers" continues its tradition of excellence. Instructors love this text because it is a comprehensive text that is easy to teach from. Students love it because it is written for them--with great pedagogy and clear explanations and examples throughout. The text features a broad array of applications, including all engineering disciplines. The revision retains the successful pedagogy of the prior editions. Chapra and Canale's unique approach opens each part of the text with sections called Motivation, Mathematical Background, and Orientation, preparing the student for what is to come in a motivating and engaging manner. Each part closes with an Epilogue containing sections called Trade-Offs, Important Relationships and Formulas, and Advanced Methods and Additional References. Much more than a summary, the Epilogue deepens understanding of what has been learned and provides a peek into more advanced methods. Approximately 80% of the end-of-chapter problems are revised or new to this edition. The expanded breadth of engineering disciplines covered is especially evident in the problems, which now cover such areas as biotechnology and biomedical engineering. Users will find use of software packages, specifically MATLAB and Excel with VBA. This includes material on developing MATLAB m-files and VBA macros.

Numerical methods for engineers

Numerical study of coupled slosh modes in a 3D vessel subjected to multi-directional excitations

A dissipative particle dynamics simulation of a pair of red blood cells in flow through a symmetric and an asymmetric bifurcated microchannel

Purpose




Better membrane oxygenators need to be developed to enable efficient gas exchange between venous blood and air.




Design/methodology/approach




Optimal design and analysis of such devices are achieved through mathematical modeling tools such as computational fluid dynamics (CFD). In this study, a control volume-based one-dimensional (1D) sub-channel analysis code is developed to analyze the gas exchange between the hollow fiber bundle and the venous blood. DIANA computer code, which is popular with the thermal hydraulic analysis of sub-channels in nuclear reactors, was suitably modified to solve the conservation equations for the blood oxygenators. The gas exchange between the tube-side fluid and the shell-side venous blood is modeled by solving mass, momentum and species conservation equations.




Findings




Simulations using sub-channel analysis are performed for the first time. As the DIANA-based approach is well known in rod bundle heat transfer, it is applied to membrane oxygenators. After detailed validations, the artificial membrane oxygenator is analyzed for different bundle sizes (L/W) and bundle porosity (epsilon) values, and oxygen saturation levels are predicted along the bundle. The present sub-channel analysis is found to be reasonably accurate and computationally efficient when compared to conventional CFD calculations.




Research limitations/implications




This approach is promising and has far-reaching ramifications to connect and extend a well-known rod bundle heat transfer algorithm to a membrane oxygenator community. As a variety of devices need to be analyzed, simplified approaches will be attractive. Although the 1D nature of the simulations facilitates handling complexity, it cannot easily compete with expensive and detailed CFD calculations.




Practical implications




This work has high practical value and impacts the design community directly. Detailed numerical simulations can be validated and benchmarked for future membrane oxygenator designs.




Social implications




Future membrane oxygenators can be designed and analyzed easily and efficiently.




Originality/value




The DIANA algorithm is popularly used in sub-channel analysis codes in rod bundle heat transfer. This efficient approach is being implemented into membrane oxygenator community for the first time.

Numerical investigation of membrane oxygenation using sub-channel analysis

Clinical observations indicate that the shape and tortuosity of the carotid siphon are some of the contributing factors to the initiation and growth of an aneurysm. The present study explores the validity of this observation by performing systematic numerical simulations. Computational fluid dynamics (CFD) based calculations are performed to compare and contrast four different types of patient-specific carotid siphons, viz., C-, S-, U-, and helical shape, to investigate the hemodynamic influences on flow features, secondary flow patterns, and helicity. Fewer curved regions and the presence of local acute curvature were found to result in higher velocity magnitude, leading to giant sidewall aneurysms in the distal end of this curvature. In contrast, a larger number of curved regions in the parent vessel resulted in disturbed flow and reduced maximum streamwise velocity. When the velocity is lower, smaller aneurysms are observed at the bifurcation carina. The influence of siphon tortuosity, which is exemplified through the Dean number and linked to secondary flows, causes higher helicity when the vessel is more tortuous. It is hypothesized that a highly tortuous vessel protects the further growth of an aneurysm. This is in contrast to a less tortuous vessel with single acute curvature and prone to further expansile behavior of an aneurysm.

Influence of carotid tortuosity on the hemodynamics in cerebral aneurysms

Understanding the detailed dynamics of bubble ebullition cycle is central to effectively exploit coolant phase change in subcooled flow boiling. In the present study, the bubble growth rate and the bubble departure mechanism in subcooled flow boiling conditions are investigated. The bubble growth rate is estimated by employing an energy balance model, which ensures that the applied wall heat flux contribution to components, such as (i) microlayer evaporation, (ii) conduction through superheated layer region, and (iii) condensation heat transfer, are accounted. The bubble departure diameter and the type of departure (i.e., sliding or lift-off) are predicted using a simple force balance model on the vapor bubble. The implemented models are thoroughly validated against the benchmark experimental data. Furthermore, the influence of operating conditions, viz., pressure (1–3  bar), mass flux (200–1000  kgm−2s−1), heat flux (200–500  kWm−2), and the degree of subcooling (20–30 K) on the bubble dynamics, is investigated for subcooled flow boiling conditions. Based on the parameters considered, a flow regime map is developed to identify the bubble departure type and diameter for a given set of operating conditions. It was noticed that the bubble departure diameter was maximum at low pressure (1 bar), low mass flux (200 kgm−2s−1), high heat flux (500 kWm−2), and low degree of subcooling (20 °C). For all the values of pressure and degree of subcooling, the sliding mode of departure was noticed at low and high mass flux values, whereas bubble departs by lift-off for moderate values of mass flux.

B.S.V. Patnaik

Papers

Numerical study of coupled slosh modes in a 3D vessel subjected to multi-directional excitations

A dissipative particle dynamics simulation of a pair of red blood cells in flow through a symmetric and an asymmetric bifurcated microchannel

Numerical investigation of membrane oxygenation using sub-channel analysis

Influence of carotid tortuosity on the hemodynamics in cerebral aneurysms

Bubble growth and departure behavior in subcooled flow boiling regime