This document describes a conjugate heat transfer analysis of an electronics cooling system using OpenFOAM. It outlines the objectives to develop a CFD model for CHT analysis and validate it with experiments. The methodology section describes the governing equations solved for fluid and solid regions as well as the interface coupling. A simple circuit board cooling case is modeled and tested. Additionally, a server cooling case is proposed with details on geometry, meshing, boundary conditions and results showing temperature distributions.
Thermal Simulations of an Electronic System using Ansys IcepakIJERA Editor
Present electronics industry component sizes are efficiently reducing to meet the product requirement with
compact size with greater performance in compact size products resulting in different problems from thermal
prospective to bring product better performance electrically and mechanically.
In this paper we will study how to overcome the thermal problem for a product which includes components
reliability and PCB performance by using CFD thermal simulation tool (Ansys Icepak).
Numerical methods in Transient-heat-conductiontmuliya
This file contains slides on Numerical methods in Transient heat conduction.
The slides were prepared while teaching Heat Transfer course to the M.Tech. students in Mechanical Engineering Dept. of St. Joseph Engineering College, Vamanjoor, Mangalore, India, during Sept. – Dec. 2010.
Contents: Finite difference eqns. by energy balance – Explicit and Implicit methods – 1-D transient conduction in a plane wall – stability criterion – Problems - 2-D transient heat conduction – Finite diff. eqns. for interior nodes – Explicit and Implicit methods - stability criterion – difference eqns for different boundary conditions – Accuracy considerations – discretization error and round–off error - Problems
MB3 is the block in which we will study Mass Balances in Multiple Phases.
Phases include: vapor-liquid equilibrium, liquid-liquid systems and solid-liquid systems
Extraction, Humidification, Distillation, Absorption, Crystallization are some examples of the processes involving these type of processes.
SECTION 1
-Single Phase Theory (Phase-Diagrams)
--Phase Diagram
--Phase Change
--Boiling point (T,P); Sublimation point (T,P); Freezing point (T,P)
--L-S, L-G, G-S Equilibrium lines
-Vapor Pressure (Antoine and Clapyeron Eqn)
--Vapor pressure definition
--Volatility
--Latent heat of vaporization
--Vapor Pressure estimation/calculation
-Gibbs Phase Rule
--DOF for systems (intensive properties)
SECTION 2
-Gas-Liquid Systems (1 condensable)
--Only one substance condenses
--Evaporation, Drying, Humidity
--Examples: Air-Humidity
-Gas-Liquid Systems (multi-condensable)
--Raoult Law
--Henry Law
--Dew/Bubble points
--Txy & Pxy Diagrams
--Examples: Ethanol-Water
SECTION 3
-Solid-Liquid Systems
--Crystallization
--Solubility and Saturation
--Hydrated Salts
-Liquid-Liquid Systems
--Miscibility and Extraction
�
Problems and Exercises in my webpage
www.ChemicalEngineeringGuy.com
Suscribe to my channel:
www.youtube.com/ChemEngineeringGuy
Visit my Facebook Page:
www.facebook.com/Chemical.Engineering.Guy
e-Mail me:
chemical.engineering.guy@gmail.com
This document summarizes simulation results for windshield defrost modeling using the heatTransfer/chtMultiRegionFoam solver in OpenFOAM. The simulation was run on a Core i7 950 processor with 24GB memory under the Mint17 operating system. The model contains exterior, ice, and cabin regions with defined patches and initial conditions. The blockMesh divides the thickness direction into 45 cells concentrated at the ice surface. Boundary conditions and initial conditions are defined for the cabin, exterior, and ice regions.
O documento discute o ciclo termodinâmico de Brayton, proposto por George Brayton em 1870. O ciclo é usado em turbinas a gás e consiste em quatro processos: compressão isentrópica, aquecimento a pressão constante, expansão isentrópica e resfriamento a pressão constante. O documento também aborda a eficiência térmica do ciclo e formas de melhorá-la, como o uso de regeneração de calor.
LMS Imagine.Lab Amesim/STAR-CCM+ co-simulation: solid stress analysis of a ge...Siemens PLM Software
This presentation introduces the co-simulation capabilities between LMS Imagine.Lab Amesim and Star-CCM+, illustrated in a gerotor pump example.
Please find more information on our website:
siemens.com/plm/simcenter
1) The document discusses heat transfer modeling in FLUENT including the energy equation, wall boundary conditions, conjugate heat transfer, thin walls, natural convection, and radiation models.
2) It provides details on the various terms in the energy equation including viscous dissipation, species diffusion, chemical reactions, and interphase transfer.
3) Examples are presented on conjugate heat transfer of an electronic chip on a circuit board using both a meshed wall and thin wall approach.
Alfa Laval owns the global trademarks for MISSIONTM and Alfa Laval. The document discusses Alfa Laval's thermal fluid heating systems, which use circulating oil to transfer heat from fuel-fired heaters to applications like fuel and cargo tanks. Alfa Laval is a global leader in supplying complete thermal fluid systems and has over 30 years of experience in the technology. Their systems offer advantages like high temperature capability, low maintenance, and global after-sales support.
This document discusses heat exchangers, including their types, performance parameters, and design methodologies. It introduces the log mean temperature difference method for relating heat transfer rate to inlet/outlet temperatures. It also describes the effectiveness-NTU method, where effectiveness is defined as the ratio of actual to maximum possible heat transfer, and NTU is the number of transfer units. Sample problems demonstrate the use of these methods to determine required surface areas, heat transfer rates, and outlet temperatures for given heat exchanger configurations and operating conditions.
Electric Process Heaters
0 INTRODUCTION/PURPOSE
1 SCOPE
2 FIELD OF APPLICATION
3 DEFINITIONS
4 ADVANTAGES OF ELECTRIC HEATERS
4.1 Safety
4.2 Environment
4.3 Location of Equipment
4.4 Low Temperature Applications
4.5 Cross Contamination
4.6 Control
5 DISADVANTAGES OF ELECTRIC HEATERS
6 POTENTIAL APPLICATIONS FOR ELECTRIC
PROCESS HEATERS
7 GENERAL DESIGN AND OPERATING CONSIDERATIONS
8 TYPES OF PROCESS ELECTRIC HEATERS
8.1 Pipeline Immersion Heaters
8.2 Tank Heaters and Boilers
8.3 Indirect (Fluid Bath) Heaters
8.4 Radiant Furnaces
8.5 Induction Heaters
8.6 Hot Block Heaters
9 CONTROL
10 REFERENCES
FIGURES
1 ELECTRIC HEAT EXCHANGER CONSTRUCTION
2 SHEATHED HEATING ELEMENTS
This document outlines Antonio Di Rienzo's PhD defense presentation on mesoscopic numerical methods for reactive flows. The presentation covers:
1. The Link-wise Artificial Compressibility Method (LW-ACM), which solves the incompressible Navier-Stokes equations using kinetic mock-up models in a similar way to the Lattice Boltzmann Method (LBM).
2. Using LBM to solve the Radiative Transfer Equation by updating intensity according to lattice velocities.
3. A new LBM model for simulating reactive flows that addresses current limitations in modeling combustion processes.
1d 2d heat transfer lumped capacitance model giesler chart, forced and free convection heat transfer, radiation, heat exchangers, boiling and condensation for mechanical engineering students
Lecture 1_Introduction to Process Modeling.pdfmudassar60
This document provides an overview of a Process Systems Engineering (PSE) course taught by Dr. Xiang Zhang. The course covers process modeling, including steady-state and dynamic modeling of lumped and distributed parameter systems. It describes the main steps of modeling as defining the problem, identifying controlling factors, analyzing data, developing mathematical models, solving models, and verifying solutions. An example of continuous stirred-tank reactor modeling is discussed. The document outlines the course content, lecturers, and examination details.
Université Badji-Mokhtar Annaba, Déprtement de Métallurgie et Génie des Matériaux
Cours de thermodynamique et cinétique chimique
Licence, Master métallurgie et génie des matériaux
1. The document discusses heat exchangers, including their types, components, heat transfer processes, and calculations.
2. It covers counterflow and parallel flow heat exchangers, log mean temperature difference, effectiveness-NTU method, and selection of heat exchangers.
3. Examples are provided to demonstrate calculations of heat exchanger length, surface area, rate of heat transfer, and the impact of fouling.
LCU14-101: Coresight Overview
---------------------------------------------------
Speaker: Mathieu Poirier
Date: September 15, 2014
---------------------------------------------------
Coresight is the name given to a set of IP blocks providing hardware assisted tracing for ARM based SoCs. This presentation will give an introduction to the technology, how it works and offer a glimpse of the capabilities it offers. More specifically we will go over the components that are part of the architecture and how they are used. Next will be presented the framework Linaro is working on in an effort to provide consolidation and standardization of interfaces to the coresight subsystem. We will conclude with a status of our current upstreaming efforts and how we see the coming months unfolding.
---------------------------------------------------
★ Resources ★
Zerista: http://lcu14.zerista.com/event/member/137703
Google Event: https://plus.google.com/u/0/events/cvb85kqv10dsc4k3e0hcvbr6i58
Presentation: http://www.slideshare.net/linaroorg/lcu14-101-coresight-overview
Video: https://www.youtube.com/watch?v=NzKPd3FByxI&list=UUIVqQKxCyQLJS6xvSmfndLA
Etherpad: http://pad.linaro.org/p/lcu14-101
---------------------------------------------------
★ Event Details ★
Linaro Connect USA - #LCU14
September 15-19th, 2014
Hyatt Regency San Francisco Airport
---------------------------------------------------
http://www.linaro.org
http://connect.linaro.org
A first order hyperbolic framework for large strain computational computation...Jibran Haider
An explicit Total Lagrangian momentum-strains mixed formulation in the form of a system of first order hyperbolic conservation laws, has recently been published to overcome the shortcomings posed by the traditional second order displacement based formulation when using linear tetrahedral elements.
The formulation, where the linear momentum and the deformation gradient are treated as unknown variables, has been implemented within the cell centred finite volume environment in OpenFOAM. The numerical solutions have performed extremely well in bending dominated nearly incompressible scenarios without the appearance of any spurious pressure modes, yielding an equal order of convergence for velocities and stresses.
To have more insight into my research, please visit my website:
http://jibranhaider.weebly.com/
The document provides programming tips for OpenFOAM, including:
1) How to get a patch's ID from its name using findPatchID().
2) How to calculate the sum of a field over a specified patch using gSum().
3) How to access boundary values of a variable on a patch and the cells adjacent to the patch using boundaryField() and faceCells().
4) How to read cell zone definitions and access cell labels using cellZones() and findZoneID().
It also demonstrates using logical operators on boolLists, appending to DynamicLists, and removing duplicate values from a list.
THERMAL ANALYSIS OF HEAT SINK (VARIABLE SHIELD PROFILE) USED IN ELECTRONIC CO...IAEME Publication
This paper deals with the comparative study of heat sink having fins with shield of various profiles namely Trapezoidal with curved and plane inclination and slope shield as heat sinks are the commonly used devices for enhancing heat transfer in electronic components. For the purpose of study heat sink is modeled by using the optimal geometric parameter such as fin height, fin thickness, base height, fin pitch as 48 mm, 1.6 mm, 8mm, 4mm and after that simulation is done at
heat load of 75W and with a air velocity of 4.7171m/s and air inlet temperature is taken as 295 K. The simulation is carried out with a commercial package provided by fluent incorporation.
Cfd and conjugate heat transfer analysis of heat sinks with different fin geo...eSAT Journals
This document discusses a computational fluid dynamics (CFD) and conjugate heat transfer analysis of different fin geometries for heat sinks used in electronics cooling. Five fin geometries - zigzag, fluted, slanted mirror, custom pin fin, and staggered array - were analyzed under different heat loads and air velocity. The results show that the slanted mirror geometry provided the best thermal performance with the lowest thermal resistance and highest heat transfer coefficient, while maintaining a relatively low pressure drop. CFD simulations using ANSYS Fluent were conducted to analyze fluid flow, heat transfer, temperature distribution, and thermal performance of the different heat sink designs.
This slide is describing how to set up the OpenFOAM simulations including rotating geometries.
The SRF (Single Rotating Frame) is covered and MRF (Multiple Reference Frame).will be covered in it.
Large strain computational solid dynamics: An upwind cell centred Finite Volu...Jibran Haider
Presented our research at the 12th World Congress on Computational Mechanics (WCCM) and 6th Asia Pacific Congress on Computational Mechanics (APCOM) at the COEX Convention Center in Seoul, Korea.
Large strain solid dynamics in OpenFOAMJibran Haider
The document describes a numerical methodology for simulating large strain solid dynamics using OpenFOAM. It proposes using a total Lagrangian formulation and first-order conservation laws similar to computational fluid dynamics to model solid mechanics problems involving large deformations. A cell-centered finite volume method is used for spatial discretization along with Riemann solvers and linear reconstruction to capture fluxes. A two-stage Runge-Kutta scheme is employed for time integration. Results are presented demonstrating the method's ability to handle problems involving mesh convergence, enhanced reconstruction, highly nonlinear behavior, plasticity, contact, unstructured meshes, and complex geometries.
OpenFOAM for beginners: Hands-on trainingJibran Haider
NOTE: A NEW VERSION OF THIS PRESENTATION IS AVAILABLE ON MY RESEARCH GATE ACCOUNT
(https://www.researchgate.net/publication/327594760_OpenFOAM_course_for_beginners_Hands-on_training)
OpenFOAM is an open source Computational Fluid Dynamics software package based on C++ programming language within the context of cell centred Finite Volume Method. It is developed primarily by OpenCFD Ltd and distributed by OpenCFD Ltd and the OpenFOAM Foundation.
Disclaimer:
This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software and owner of the OPENFOAM® and OpenCFD® trade marks. "
The document discusses heat sink design for cooling electronic systems. It compares different design parameters for heat sinks such as fin material, shape, layout, working fluid, and forcing mechanism. Two solutions are proposed to achieve heat fluxes of 2 kW/m^2 and 200 kW/m^2. Sample calculations are shown for a solution using aluminum fins cooled by forced air to achieve a heat flux of 209.377 kW/m^2.
The document summarizes heat sink design optimization work for a server. Experimental testing and CFD simulation were used to analyze different heat sink designs. Mini partitions and a vapor chamber were added, lowering the PCH temperature by up to 6%. The final design was a cross-cut fin heat sink with both modifications. The CFD model results matched experimental testing within 5%, validating the simulation approach.
Water cooled minichannel heat sinks for microprocessor cooling: Effect of fin...Danial Sohail
Heat sink with different fin spacing mounted on a microprocessor were tested for their heat removing capabilities by varying coolant flow rates over them
This document provides a guide to the tutorials available in OpenFOAM. It lists the main solvers in OpenFOAM along with the tutorials available for each solver. For each tutorial, it provides a short description of what is done in the tutorial and includes screenshots of results when relevant. The goal is to help users quickly select relevant tutorials to learn the basics of using different solvers in OpenFOAM. It focuses on describing the tutorials at a procedural level rather than going into the underlying physics and equations.
This document describes the design of a machine to directly measure contact line friction. It begins with background on contact angles and traditional methods for measuring contact angle hysteresis and pinning forces, such as the captive needle and tilting plate methods. It then outlines the design of the new measurement machine, which uses a linear motion stage to drag a droplet across a substrate while measuring the pinning force with a high precision force sensor. Key aspects of the design include the motion system, force sensor, vibration isolation, sample holder, enclosure box, and various adjustments for positioning and observation. The machine aims to directly measure pinning forces with improved accuracy compared to existing indirect methods.
Enhancement of open-source computational-fluid-dynamics code OpenFOAM(R) for Detached-Eddy Simulation as presented on the OSCIC'10 in Munich 2010 by Dr. Ulf Bunge.
Suse Studio: "How to create a live openSUSE image with OpenFOAM® and CFD tools"Baltasar Ortega
Una descripción de Suse Studio, además de una magnífica explicación de su utilización de la mano de Alberto Passalacqua.
"How to create a live
openSUSE image with OpenFOAM® and CFD tools"
The document summarizes research on using composite thin films embedded with microencapsulated phase change materials (PCMs) for passive thermal management of electronic devices. Through modeling, three effective configurations were identified that reduced temperature peaks and cooling times compared to a plain film. The first two used a 450μm thick film with 40% and 60% PCM, reducing temperature by 9% and 13% respectively. The third used a 400μm thick film with 60% PCM, reducing temperature by 11% and cooling time by 19%. Larger thicknesses and PCM fractions generally improved thermal performance but required optimizing the phase change temperature. Future work could refine the model and experimentally validate the results.
The document summarizes an internship focused on optimizing lubrication and cooling systems in cold rolling mills. The primary aims of the internship are to improve heat transfer modeling in the thermal simulator and identify actions to increase cooling efficiency. Work is underway to integrate a new rheology model into the simulator and analyze results from tests of 13 different roll cooling configurations. The intern will next monitor pilot tests comparing 6 configurations to determine an optimized cooling setup.
This document summarizes the development of a high-power lithium target for accelerator-based boron neutron capture therapy (BNCT). Key points:
- A water-cooled conical target is being developed to accept a 50 kW proton beam and produce neutrons via the 7Li(p,n)7Be reaction for BNCT applications.
- Computational fluid dynamics modeling was used to design the target with 20 helical water channels to keep the lithium surface below 150°C with a water flow of 80 kg/min.
- An initial prototype target was fabricated and underwent preliminary hydraulic testing matching CFD predictions. Further electron beam thermal testing is planned at Sandia National Laboratories to validate the
Performance prediction of a thermal system using Artificial Neural NetworksIJERD Editor
This document summarizes a study on using artificial neural networks (ANNs) to predict the performance of a condenser system and assess fouling over time. Experiments were conducted on an industrial condenser to collect temperature and flow rate data. An ANN model was developed and trained to predict the overall heat transfer coefficient of the clean condenser system based on the input parameters. The model was then used to calculate the fouling factor by comparing the predicted clean performance to the actual performance measured over time, indicating degradation due to fouling on the heat transfer surfaces. The developed system provides a method to monitor condenser performance and identify when cleaning is needed to improve efficiency.
This document summarizes a student project analyzing the effectiveness of different fin designs on a CPU heatsink. It includes the design of a base heatsink model with 30 fins and details simulations varying the number of fins, fin height, and rotation speed. The results show that increasing the number of fins from 30 to 50 lowers the maximum surface temperature by 5.4°C, while further increases do not impact temperature as much.
This document discusses modeling heat transfer in a tube heat exchanger using analytical and numerical methods. It will use MATLAB and FLUENT software to model the heat transfer process. Governing equations for the heat transfer and fluid flow are presented. Initial and boundary conditions are defined to solve the equations numerically using a 4th order Runge-Kutta method in MATLAB. The goal is to investigate the results from analytical, numerical and simulation approaches to optimize the heat exchanger efficiency.
This thesis examines approximate dynamic programming methods for optimizing residential water heating systems. It formulates the water heating problem as a discrete-time, average cost periodic Markov decision process (MDP) to minimize operating costs and customer discomfort over an infinite horizon. The methodology section describes solving the MDP using finite-horizon dynamic programming and average cost dynamic programming for periodic MDPs via relative value iteration. Results from numerical simulations will evaluate different control policies, including set-point methods, temperature aggregation, and usage history aggregation. Extensions to solar water heating and automated demand response are also discussed.
This document discusses methods for assessing the energy performance of heat exchangers over time. It describes calculating the overall heat transfer coefficient U to determine if fouling or other issues have reduced efficiency. The procedure involves monitoring operating parameters, calculating thermal properties, and determining U by measuring the heat duty, surface area, and log mean temperature difference. An example application to a liquid-liquid exchanger is provided, comparing test data to design specifications to identify potential fouling issues.
This document summarizes a study that uses computational fluid dynamics (CFD) to analyze the flow and heat transfer characteristics of different elliptical pin fin heat sink configurations. The study aims to examine the effects of elliptical pin fin design parameters on heat sink performance. CFD simulations are conducted using ANSYS Fluent software to compute the thermal resistance, heat transfer coefficient, and surface Nusselt number of plate fin and elliptical pin fin heat sinks. Results show that the elliptical pin fin heat sink has better heat transfer performance than the plate fin heat sink. Simulation results for thermal resistance and pressure drop are validated against experimental data from previous studies.
This document reports on analyses and experimentation conducted on a Quadrafire Classic Bay 1200 pellet stove to increase its efficiency and reduce emissions. The researchers have developed a theoretical model that characterizes the stove's heating output as a continuous function of distance along the heat exchanger, accounting for both convection and radiation heat transfer. Their experiments involved upgrading instrumentation to directly measure airflow rates and temperatures, validating assumptions in the previous model. Key findings include radiation significantly influencing heat output and an air-rich firepot reaction that reduces combustion efficiency. The updated model and experiments provide a foundation for redesigning the stove.
This document provides an overview of Kern's method for designing shell-and-tube heat exchangers. It begins with objectives and an introduction to Kern's method. It then outlines the design procedure algorithm and provides an example application. The example involves designing an exchanger to sub-cool methanol condensate using brackish water as the coolant. The document walks through each step of the Kern's method design process for this example, including calculating properties, determining duties, selecting tube/shell parameters, and estimating heat transfer coefficients.
High Pressure Die Casting Cooling calculation with application of ThermodynamicsIRJET Journal
This document discusses calculations for cooling channel design in high pressure die casting tools. It begins with an overview of high pressure die casting and importance of cooling. Thermodynamic principles of heat transfer via conduction, convection and radiation are explained. Equations for heat transfer rate via different modes are provided. An example calculation is presented for a die casting tool producing an aluminum contactor housing. Total heat input, weight of molten aluminum, and heat transfer coefficients are calculated. Based on the calculations, cooling channel length and depth are determined as 299.53mm and 15.23mm for the moving side, and 258.74mm and 22.56mm for the fixed side. This resolves prior soldering issues.
Understanding and predicting CO2 properties - Presentation by Richard Graham in the Effects of Impurities on CO2 Properties session at the UKCCSRC Cardiff Biannual Meeting 10-11 September 2014
IRJET- Experimental Evaluation of Shell & Tube Heat Exchanger with P – Toluid...IRJET Journal
The document describes an experimental evaluation of a shell and tube heat exchanger using p-Toluidine as a phase change material (PCM). P-Toluidine has a melting temperature of 44°C and was filled in the tubes of the heat exchanger. Hot water was circulated through the shell side while cold water was passed through the tubes. The temperatures of the hot and cold water inlet and outlets were measured over time. Results showed that the effectiveness of the heat exchanger, or the maximum heat transfer, increased from 38% initially to 54% as the hot water temperature increased to 60°C and melted the PCM. The heat exchanger performance was found to improve with the PCM melting and storing thermal
IRJET- Experimental Evaluation of Shell & Tube Heat Exchanger with P – Toluid...IRJET Journal
The document describes an experimental study of a shell and tube heat exchanger using p-Toluidine as the phase change material (PCM). P-Toluidine has a melting temperature of 44°C and was selected as the PCM due to its suitable thermo-physical properties. Water was used as the working fluid flowing through the tubes. Experiments were conducted to evaluate the temperature differences in the shell and tube heat exchanger with variations in the mass flow rate. The results showed that the effectiveness was higher when the PCM was fully melted compared to during the melting process. Equations related to heat exchanger effectiveness, maximum possible heat transfer, and the number of transfer units were also presented.
This document summarizes research on modeling heat transfer and bandwidth determination for shape memory alloy (SMA) actuators used in robotics applications. Key points include:
- A heat transfer model was developed using energy balance equations to analyze bandwidth and efficiency of SMA actuators.
- Bandwidth increases with higher air speed/lower temperature and smaller wire diameter, as convective heat transfer increases and heating/cooling times decrease.
- Transformation bandwidth, which accounts for phase changes, was determined and increases with smaller wire diameter and higher air speed/lower temperature.
- An equation was developed to calculate efficiency of SMA actuators based on material properties, heat transfer, latent heat of transformation, and temperature
Photovoltaic thermal (PV/T) collectors with nanofluids and nano-Phase Change ...Ali Al-Waeli
The presentation is derived from my PhD viva presentation which focuses on the topic of Photovoltaic thermal (PV/T) collectors with nanofluids and nano-Phase Change Material.
Presented by: Dr. Ali Hussein A. Alwaeli
Photovoltaic thermal (PV/T) collectors with nanofluids and nano-Phase Change ...
Avinash_PPT
1. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Conjugate Heat Transfer for Electronic Cooling
using OpenFOAM
Avinash Valiba Gorde(CMS-1310)
Internal Guide
Dr. Sukratu Barve
Dr. Mihir Arjunwadkar
Centre for Modelling and Simulation
External Guide
Dr. Vikas Kumar
Mr. Mohan Labade
Centre for Development of Advanced Computing Pune,INDIA
June 20, 2015
2. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Outline
1 Need of Project
2 Objectives
3 Methodology
4 CHT Case in OpenFOAM
5 Server Case
6 Preprocessing
7 Results
8 Conclusion
9 Future Scope
10 References
3. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Need of Project
Considering following aspect,
A very limited conjugate heat transfer case studies have been
analyzed using OpenFOAM.
4. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Objectives
To develop a CFD model for conjugate heat transfer analysis
of an electronics cooling system
Validation of CFD model with experimental results
5. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Methodology
Conjugate Heat Transfer
In CHT computational domain is divided into fluid and solid
regions. The NaviereStokes equations and the energy equation in
the fluid regions are solved first. Then the heat transfer equation
in the solid regions is solved.
Figure : Conjugate Heat Transfer
6. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Steady state heat transfer
Steady state conditions the temperature within the system does
not change with time.
Steady State Equation for Solid
− · (κ T) = Sh (1)
Unsteady state heat transfer
Unsteady state conditions the temperature within the system
changes with time.
Unsteady State Equation for Solid
ρCp
∂T
∂t
= · (κ T) + Sh (2)
7. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Solver
chtMutiRegionSimpleFoam Solver
It is steady state solver for conjugate heat transfer. It is
combination of heatConductionFoam and buoyantFoam for
conjugate heat transfer
Governing Equation Used by OpenFoam
chtMultiregionSimpleFoam Solver
Continuity Equation
∂ρ
∂t
+
−→
· (ρ−→u ) = 0 (3)
8. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Momentum Equation
The sum of all forces acting on a body is equal to the rate of
change of momentum of it
1) Body forces = gravitational, centrifugal
2)Surface forces acting on the surface of the body e.g. pressure or
viscous forces
∂ (ρu)
∂t
+ · (ρuu) = − p + · τ − ρg (4)
Energy Equation
∂(ρE)
∂t
+ · (ρuE) = · (α e) − · [u (p − ρgr)] + Sh (5)
α =
µCp
κ is the thermal diffusivity.
Sh is the thermal source term.
e is the specific internal energy.
E is the specific total energy of the gas defined as:
E = e + u2
2
9. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Equation in solid region
∂ (ρh)
∂t
= · (α h) + Sh (6)
α =
µCp
κ is the thermal diffusivity.
Sh is the heat source term.
h is the sensilbe enthalpy.
10. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Coupling at the solid-fluid interface
Tw =
ks
δxs
Ts
kf
δxf
+ ks
δxs
+ (1 − k) Tf +
Tref
δxf
=
kf
δxf
Tf + ks
δxs
Ts
kf
δxf
+ ks
δxs
(7)
Tf and Ts are the temperatures of the fluid and solid regions
kf and ks are the heat conductivities of the fluid and solid
Tref is a reference temperature
δxf and δxs are the distances between the cell centers
Figure : Interface Temp
11. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
CHT Case in OpenFOAM
Geometry For Simple Circuit Board Cooling
Figure : Simple cht case
12. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Meshing For Simple Circuit Board Cooling
Figure : Simple cht case
13. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Implementation of Heat Source
Heat source
In src/fvOptions/lnInclude
const dimensionSet ds =
rho.dimensions()*fld.dimensions()/dimTime*dimVolume;
ds = ρ * h/time * volume
ds = kg/m3 * J/(kg * S) * m3
ds = W (watt)
tmp <fvMatrix<Type > >tmtx(new fvMatrix < Type >(fld, ds));
h (200 0);
14. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Boundary Conditions
Solid Region
Initial Temperature=300K.
Wall is zeroGradient.
Heat Source 200 Watt (4 IC’s).
Interface between solid(IC), fluid(AIR) Mixed B.C.
Fluid Region
Velocity v=0.3m/s.
p = atmospheric.
Initial Temperature (T) = 299K.
Inlet=fixed Value Outlet=inletOutlet for Velocity &
Temperature.
Walls maintained at symmetry, zeroGradient.
Interfaces between solid and fluid = Mixed B.C.
15. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Temperature contour at mid plane
Figure : Temperature variation in solid along xaxis
16. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
300
320
340
360
380
400
420
0 0.2 0.4 0.6 0.8 1 1.2 1.4
TemperatureinKelvin
Distance along X-axis
'100_watt.csv' u 8:2
'200_watt.csv' u 8:2
Figure : Temperature Distribution
17. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Server Case
18. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Server Case
Geometry of Server
This is Intel Server Board S2600GZ/GL having 2 SOCKETS, 24
RAMS, PCB, 2SMPS
19. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Generalized geometry of Server
Figure : Generalized model of server
20. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Sink of Socket
High Thermal conductivity, larger surface area
Low thermal resistance
Material Aluminium, Copper
Air can be easily flow through it
Figure : Sink of Socket
21. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Meshing
For server create 3 blocks specifying coordinates of each point.
Dimensions are L=0.46m Zaxis, W=0.41m Xaxis, H=0.09m Yaxis
it is structured mesh
Figure : blockMesh for Server
22. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
topoSet
topoSet utility is use to create different domain inside the mesh
block.
Figure : blockMesh for Server
23. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Turbulance Model
RNG κ − turbulence model
Useful for rapidly strains flows, Heat and mass transfer
Low Reynold number effect
Standard κ − model for high reynolds number
Turbulent Kinetic Energy (k)
∂ (ρk)
∂t
+
∂ (ρkui )
∂xi
=
∂
∂xj
αk µeff +
∂k
∂xj
+ Gk + Gb − ρ − Ym + Sk (8)
Dissipation rate ( )
∂ (ρ )
∂t
+
∂ (ρ ui )
∂xi
=
∂
∂xj
αk µeff +
∂
∂xj
+ C1
k
(Gk + C3 Gb) − C2ρ
2
k
− R + S (9)
The coefficient updated each iteration
C
∗
1 = c1 −
η 1 − η
η0
1 + βη3
(10)
η =
k
τij τij (11)
24. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Table : Constant values for turbulence model
C1 C2 C3 Cµ η β σk σ
1.42 1.68 -0.33 0.085 4.38 0.012 0.7194 0.7194
25. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Boundary Conditions
Solid Region
Initial Temperature=299K.
Walls zeroGradient.
Heat Source :- 103 Watt (2 sockets), 4 Watt (on 1, 4, 9, 12,
13, 14, 16, 21, 24 ).
Interface between solid(SOCKETs, RAMs), fluid(AIR) Mixed
B.C.
Fluid Region
Velocity v=2.1m/s.
p = atmospheric.
Initial Temperature (T) = 299K.
Inlet=fixed Value Outlet=inletOutlet for Velocity &
Temperature.
Walls maintained at symmetry, zeroGradient.
Interfaces between solid and fluid = Mixed B.C.
26. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Solid Thermophysical Properties
Table : Solid thermophysical properties
Solid name κ(W/mK) Cp(J/kgK) rho(kg/m3) Material
SOCKET 1 & 2 124 702 2325 Silicon
SINK 1 & 2 220 903 2700 Aluminium
RAM 124 702 2325 Silicon
SMPS 1 & 2 53 450 8000 Steel
BOX 1 & 2 53 450 8000 Steel
FLAPPER 53 450 8000 Steel
BOARD 0.02 1800 1500 FR4
27. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Table : Fluid(AIR) thermophysical properties
Parameter value
Cp (J/kg K) 1000
µ 1.8e−05
Pr 0.7
28. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Results
Temperature Contours
Figure : Temperature distrubution in Server
29. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Temperature along Heat Sink
coordinates are (0.1 0.02 0.28) and (0.1 0.02 0.360)
298
300
302
304
306
308
310
312
314
316
318
320
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08
TemperatureinKelvin
distace along z-axis
SinkTemp betn fins
Figure : Temp. distrubution of fin along z-direction
30. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Temperature between fin passage
coordinates are (0.06 0.02 0.332) & (0.19 0.02 0.0332)
298
300
302
304
306
308
310
312
314
316
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Temp.inKelvin
distace along x-axis
HeatsinkTemp
Figure : Temp. distrubution of fin along x-direction
31. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Temperature Variation of RAMs & SOCKETs
coordinates are (0.095 0.0075 0) & (0.09 0.0075 0.46)
295
300
305
310
315
320
325
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Temp.inKelvin
distace along z-axis
SOCKET & RAM Temp
Figure : Temperature distribution along z-direction
32. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Comparision with experimental results
Comparision of Sockets & Rams Temperature obtained from cfd
model & experiment
Table : Comparision of Socket & Ram Temp(K)
Region OpenFOAM Result Experimental Result
SOCKET 319.838 318
RAM 301.43 303
33. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Error in Results
The % error in Sockets & Rams Temperature :
%error =
Resultexp − Resultcfd
Resultcfd
∗ 100 (12)
Sockets Temperature % error
%error =
318 − 319.838
319.838
∗ 100 = 0.57% (13)
Rams Temperature % error
%error =
303 − 301.43
301.43
∗ 100 = 0.52% (14)
34. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Conclusion
Results obtained from cfd model are good, as compared to
experimental results.
This CFD results gives three dimensional temperature
distribution in server components such as RAMs, SOCKETs,
BOARD
The % error in Sockets temperature is 0.57 % and Rams
temperature is 0.52 %
35. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Future Scope
To implement Fan boundary condition for mass flow rate
To implement micro thickness baffle for sink region
36. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
References
ANSYS.
Fluent 6.3 users guide.
www.sharcnet.ca/Software/Fluent6/html/ug/node479.htm.
Yunus A. Cengel.
Heat Transfer A Practical Approach.
McGraw-Hill, forth edition, 2003.
Yongsheng Lian Chaolei Zang.
Conjugate heat transfer analysis using simplified household refrigerator model.
Technical report, Department of Mechanical Engineering University of Lousville USA, June 2014.
Porterie B. Consalvi J.L., Pizzo Y.
Numerical analysis of heating process in upward flame spread over thick pmma slab.
Fire Safty J, V43(N5):351–362, 7/08.
John D.Anderson.
Computational Fluid Dynamics Basic with application.
Mc GrawHill, 2000.
Wilcox D.C.
Turbulence Modeling for CFD.
DCW Industries, 3 edition, 2006.
D.E.Dwyer.
Defining ventilation boundary condition for green house climate model.
Master’s thesis, Delft University of Technology, August 2014.
OpenFOAM Foundation.
Openfoam programmer guide.
www.openfoam.org.
OpenFOAM Foundation.
Openfoam user guide.
www.openfoam.org.
37. Need of Project Objectives Methodology CHT Case in OpenFOAM Server Case Preprocessing Results Conclusion Future Sco
Thank You !