Advanced Tutorial: Natural Convection of a LED Spotlight

This tutorial shows how to perform a natural convection simulation on a LED spotlight, aiming to satisfy maximum temperature requirements and optimize LED performance.

This tutorial teaches how to:

• Set up and run a natural convection simulation using the Conjugate heat transfer v2.0 solver;
• Assign boundary conditions, material, and other properties to the simulation;
• Mesh with the SimScale standard meshing algorithm.

We are following the typical SimScale workflow:

1. Prepare the CAD model for the simulation.
2. Set up the simulation.
3. Create the mesh.
4. Run the simulation and analyze the results.

1. Prepare the CAD Model and Select the Analysis Type

As a first step, please click on the button below. It will import the tutorial project directly into your Workbench.

The following picture demonstrates what should be visible after importing the tutorial project:

Note that the project contains two geometries. The first one consists of the complete LED spotlight geometry, and the second one consists of a 90 degrees slice.

1.1. CAD Edit

The initial CAD geometry only contains the LED spotlight parts. Before starting to work on the setup, we have to create the external flow volume, and also imprint the geometry.

The external flow region should be large enough to prevent the boundary conditions from interfering with the results. In Figure 5, you will find a good rule of thumb for the flow volume dimensions:

• Top: 6 times the height of the LED spotlight
• Bottom: 4 times the height of the geometry
• Sides: 8 times the width of the geometry

To perform the necessary operations, we create a draft of the existing CAD model in the same workbench window, by clicking on the ‘Edit a copy’ option as shown below. You’ll then be able to perform several CAD operations to the draft within the same user interface.

Once the draft is created, we can proceed to create an External Flow Region. Based on the rule of thumb presented in figure 5, we can set the minimum and maximum coordinates:

1. Create an ‘External’ flow volume operation
2. Set the flow volume dimensions to the following:
   • X min value: ‘0’ meters
   • Y min value: ‘-0.3’ meters
   • Z min value: ‘0’ meters
   • X max value: ‘0.25’ meters
   • Y max value: ‘0.4’ meters
   • Z max value: ‘0.25’ meters
   • Hit ‘Apply’ to run the operation.

Another important step is to run an Imprint operation. In our LED geometry, a series of interfaces between solid/solid and solid/flow regions are present. The imprint operation improves the automatic detection of such interfaces.

Therefore, please follow the steps below:

Once you finish running the imprint operation, please make sure to click on ‘Save’. This will save the latest version of your draft as a new model under the Geometries list.

1.3. Create the Simulation

Once the new geometry is imported to the Workbench, it will be named Copy of LED Spotlight – Quarter Geometry. In case you would like, you can always change the name of the new geometry to something more representative. After saving the name change, we are ready to ‘Create a Simulation’.

The analysis type choice widget opens up. From the options, pick ‘Conjugate heat transfer v2.0’ and ‘Create the Simulation’:

As a result of a strongly coupled energy equation, the Conjugate heat transfer v2.0 analysis type provides much faster convergence rates, in comparison to the conventional Conjugate heat transfer (CHT).

After creating a simulation, the simulation tree will be visible in the left-hand side panel. It’s necessary to set up all entries to be able to run the simulation. At this point, SimScale will also automatically detect all Contacts within the geometry. In our case, there are a total of 64 interfaces.

The global simulation settings are shown in Figure 11. Since it’s a natural convection simulation, the velocities in the domain are small, and a ‘Laminar’ turbulence model is appropriate.

2. Simulation Setup

In the following sections, we will set up the physics of the simulation.

2.1. Model

In the Model tab, it’s possible to define the direction and magnitude of gravity:

Please define gravity as ‘-9.81’ $m/s^2$ in the y-direction.

2.2. Defining the Materials

In conjugate heat transfer simulations, it’s necessary to define both solid and fluid materials. Let’s first go through the fluid domain.

2.2.1 Fluid Materials

To add a fluid material, please click on the ‘+ button’ next to Fluid:

The fluid material library pops up. From the list, choose ‘Air‘ and hit ‘Apply’:

The geometry contains two air regions: one internal to the LED, and the external flow region. You can use the right-hand side panel to quickly select the volumes.

2.2.2 Solid Materials

Let’s now add the solid materials. First, click on the ‘+ button’ next to Solids, as in figure 12. This time, a library of solid materials appear.

From the list, please choose ‘PVC’, which is a common material for circuit boards:

The materials library contains common properties for PVC. In this tutorial, PVC is used in a circuit board, so we will adjust the material properties to account for the circuits. Figure 17 shows the changes:

• Set $(\kappa )$ Thermal conductivity to ‘2.5’ $\frac{W}{m.K}$.
• The new Specific heat is ‘1004’ $\frac{J}{kg.K}$.
• Set $(\rho )$ Density to ‘1900’ $\frac{kg}{m^3}$.
• Lastly, assign the Board volume.

All the remaining solid materials maintain the default settings. Find below a list of materials and the respective assignments.

• Aluminium: Core
• Brass: Connector Base
• Copper: Connector Pin
• Glass: Glass
• Silicon: 0.25 W Chip, 0.5 W Chip, and 1 W Chip

2.3. Initial Conditions

In CFD simulations, it’s a good practice to initialize the parameters close to the expected solution. With this approach, it’s possible to improve the convergence rate of a simulation, achieving the final result faster.

In this tutorial, we will change the initialization of temperature and velocity.

2.3.1 Temperature

The LED geometry contains three chips, which are dissipating heat. From a design perspective, it’s important to keep the chip temperatures under control, otherwise, the lifetime of the LED, along with the LED performance, becomes shorter.

Usual temperatures on the chips are around 100 $°C$, thus, initializing the LED parts at 75 $°C$ is a good initial guess. Figure 18 shows the initial steps for initialization with subdomains.

1. Add a new subdomain by clicking on the ‘+ button’ next to Subdomains.
2. Input the initialization value, in this case, ’75’ $°C$.

We want to assign all the LED parts to this subdomain. The quickest way to do that is by using an Invert visible assignment function. Please proceed as below:

1. Select the Flow region volume.
2. By right-clicking in the viewer, a window with options appears.
3. Click on Invert visible assignment. This way, the 9 volumes of interest are assigned.

2.3.2 Velocity

As a result of the high temperatures on the LED geometry, a natural convection plume develops. Therefore, we strongly recommend initializing the velocity field around the LED, which enhances the convergence rate of the simulation.

In this tutorial, we will initialize the velocity field using a cylinder geometry primitive, which is a good option for the shape of the LED geometry. Figure 21 shows the initial steps:

1. Add a new subdomain by clicking on the ‘+ button’ next to Subdomains.
2. Using the orientation cube as a reference, initialize the velocity in the y-direction with ‘0.05’ $m/s$.
3. Click on the ‘+ button’ next to Geometry primitives and select a ‘Cylinder’.

Now, a window opens up, where the user can define the cylinder orientation and dimensions. Using a cylinder slightly wider than the LED geometry yields good results. Therefore, please input the values from Figure 22:

After saving the cylinder primitive, it is automatically assigned to the velocity subdomain.

2.4. Boundary Conditions

For boundary condition definition, we will use Figure 23 as reference:

2.4.1 Natural Convection Inlet/Outlet

The natural convection inlet/outlet boundary condition is exclusive to SimScale. This boundary condition allows fluid to go in and out of the domain, therefore it’s a good option when it’s not clear what the flow behavior will be. For more details on the mathematical implementation, see this dedicated documentation page.

To create a new boundary condition, click on the ‘+ button’ next to Boundary conditions, and select the desired type from the drop-down menu.

Using figure 23 as a reference, the top, bottom, and both side faces will receive a natural convection inlet/outlet boundary condition.

2.4.2 Symmetry

Please create a second boundary condition, as in figure 24. This time, however, select a Symmetry condition:

2.5. Advanced Concepts

Additionally, the geometry contains three chips, which will be defined as power sources. Figure 28 contains further details:

In SimScale, there are two types of Power sources:

• Absolute power source: the user defines the heat flux in $W$ or $\frac{btu}{s}$
• Specific power source: the definition of the flux is in $W/m^3$ or $\frac{btu}{s.i{n}^3}$

Hence, the Absolute power source works best in our case, as we already know the exact heat flux in each chip.

Please proceed to set up the first power source, as in Figure 29:

1. Create a new power source by clicking on the ‘+ button’ next to Power sources.
2. Select an ‘Absolute power source’ and set the heat flux to 1 watt.
3. Assign the 1 W Chip volume to this power source, which can be done using the right-hand side geometry panel in the Workbench.

Please create two new power sources, for the 0.5 W Chip and the 0.25 W Chip volumes. Follow the same steps outlined above, remembering to adjust the Heat flux accordingly.

2.6. Numerics & Simulation Control

The default settings in the Numerics tab work well for most simulations. In this tutorial, no changes are required.

In the Simulation control tab, it’s possible to define a series of parameters that control the simulation process, including the number of iterations to perform, and the maximum runtime for the simulation run. Figure 30 highlights the changes in the simulation control tab:

• Define an End time of ‘750’ iterations and a Write interval of ‘750’ iterations. With these settings, only the result set from the last iteration will be written, which is a common practice for steady-state simulations.
• Change the Maximum runtime of the simulation to ‘30000’ seconds.

2.7. Results Control

Setting up result controls is a crucial step in a simulation. Since they help us to assess convergence, it’s important to set meaningful result controls for our parameters of interest.

For a LED simulation, some parameters of interest are, for example, the temperature on the chips. Therefore, we will set Area average result controls on them. To do that, please proceed as in Figure 31:

1. Create a new Area average control by clicking on the ‘+ button’ next to Surface data.
2. Define a Write interval of ‘5’ iterations. This way, the algorithm will plot the area averages every 5 iterations, which is good enough for a steady-state simulation.
3. To make the face selection easier, please hide the Board volume, by clicking on the ‘eye icon’ next to it.
4. For the first result control, please select the top face of the 1-watt chip.

Afterward, following the same procedure, please create two additional area average controls, for the top faces of the 0.5 W Chip and the 0.25 W Chip.

3. Mesh

For the meshing operation, we recommend using the standard algorithm, which is a good choice in general, as it is quite automated and delivers good results for most geometries.

After clicking on Mesh in the simulation tree, a setup window opens. The main default settings are good for an initial study. Under Advanced settings, please set the Gap refinement factor to ‘3’.

With these mesh settings, we already have a good mesh that allows us to run a successful simulation.

In the box below, we will show optional steps to configure a volume custom sizing for the mesh, which allows us to capture the convection plume more accurately. If you would like to proceed without volume custom sizing, please skip to section 4.

4. Start the Simulation

To start a new simulation, please click on the ‘+ button’ next to Simulation Runs.

This way, the mesh will be generated, and afterward, the simulation run will start automatically. While the simulation results are being calculated, you can already have a look at the intermediate results in the post-processor. They are being updated in real-time! You will also find a link to the finished project at the end of the tutorial.

The simulation run takes from 1 to 2 hours to finish. At this point, you can access the post-processing environment by clicking on ‘Solution Fields’ or ‘Post-process results’:

5. Post-Processing

5.1 Result Controls

Before heading over to the post-processing environment, let’s inspect the data obtained with the result controls. The average chip temperatures are available under the ‘Area averages’ tab:

The plot shows how the chip temperature is evolving as the algorithm performs iterations. In a converged simulation, the temperatures will stabilize and no longer change between iterations. As a result of the strongly coupled formulation from the CHT v2.0 algorithm, the temperature converges quickly, in just over 300 iterations.

For more notes on assessing convergence in CFD simulations, please refer to this article.

5.2 Temperature Contours

For electronics cooling applications, it is useful to visually evaluate the temperature on and around the chips. After accessing the post-processor by clicking on ‘Solution Fields‘, please proceed as below:

1. Set the Coloring under Parts Color to ‘Temperature’. Now the temperature is shown in the viewer.
2. To make the visualization better, we can also adjust the render mode. Both Surfaces with wireframe and Translucent surfaces with wireframe work nicely for this purpose.
3. On the right-hand side panel, click on the ‘eye icon’ next to Board. This way, we will be able to see the temperatures on the chips clearly.

Now that we have a free view of the chips, we can also use the Inspect point feature to dynamically inspect points of interest. The image below shows the steps:

1. Select the ‘Inspect point’ feature on the top bar
2. After making sure that the Board volume is hidden, you will be able to inspect the temperature directly on the chips. Adjust the camera angle and hide other volumes if necessary

5.3 Particle Traces and Animations

After accessing the post-processor by clicking on ‘Solution Fields‘, you can use several post-processing filters to further analyze the results.

For example, the Particle Trace filter tracks the particles as they move through the domain. This filter is very useful to capture the convection plume that forms due to the convection effects. Please proceed as below:

1. Create a new filter by clicking on the ‘Particle Trace’ button on the top ribbon
2. Make sure that the Pick Position button is on. This way, you can define the seeds for the traces
3. By selecting a face under the LED, you will be able to see the convection plume forming

It is possible to create even more particle trace filters, so we can see how the air moves from multiple different positions. For example, some air is drawn from the sides of our flow region, due to the convection plume that is generated. Create a second particle trace filter, this time, selecting one of the side walls:

Particle trace filters are even more powerful when animated. Now that we have a couple of particle traces, let’s create an ‘Animation’ filter for them. Remember to change the Animation type to ‘Particle Trace’ before starting the animation:

After pressing the play button (step 3 in the image above), we can see how the particles travel through the domain in the animation:

With the current configuration, we are getting temperatures around 375 $K$ on our chips. With the insights obtained from this initial simulation, we have more information to optimize the following iterations and achieve lower temperatures.

Have a look at our post-processing guide to learn how to use other filters from the post-processor for this simulation.

Last updated: June 30th, 2025