Grid Sensitivity Study Research Articles

Numerical approach for the simulation of flow-induced noise around a structure with complex geometry: High-speed train bogie in a cavity

The bogie region is a significant source of aerodynamic noise on high-speed trains. Owing to its complex geometry and flow field, numerical simulations using computational fluid dynamics are especially challenging. The main challenge is to achieve a grid with adequate resolution, especially in the boundary layer, while ensuring computational affordability. This challenge is addressed here by employing a hybrid grid, integrating a structured hexahedral mesh near solid surfaces with an unstructured polyhedral mesh in the remaining volume. To limit the number of cells in the boundary layer region, the delayed detached eddy simulation method is adopted. Additionally, to achieve a further reduction in the cell count, the Reynolds number of the model is decreased by scaling down the model size and lowering the inflow speed. The hybrid grid generation and numerical settings are guided by validated simulations of flow over cylinders. A grid sensitivity study, conducted with a simplified half-width bogie model, reveals the meshing requirements for the full-width model. Aerodynamic results highlight the rear section of the cavity and bogie as primary noise sources, emphasizing the critical role of the detached shear layer from upstream components. Time-resolved surface pressure data are input into the Ffowcs Williams–Hawkings equation for far-field noise calculation. The results indicate that the sound energy is concentrated below 200 Hz in the full-scale model, with the cavity contributing more than the bogie. This study provides a practical numerical approach for simulating a structure with complex geometry, offering insights for realistic model simulations.

Numerical simulations of temperature loads on multilayer Laue lenses

We present numerical simulations of the heat loads on novel diffractive X-ray optics, known as multilayer Laue lenses, exposed to high-intensity X-ray beams produced by an X-ray free electron laser (XFEL). These lenses can be used to focus XFEL beams to nanometer spots. The temperature distribution within the lens and temperature evolution as a function of incident pulse frequencies were calculated for two different lens geometries and several material pairs and material ratios of the MLLs. Simulations considered the special pulse structure of European XFEL with X-rays being delivered in pulse trains. After defining the geometric model, computational grid, material properties, and boundary conditions, a grid sensitivity study was carried out. We solved the transient heat energy transport equation in solids for mixed boundary conditions. The results of these simulations will help select materials and lens geometry for future XFEL experiments.

Thermochemical non-equilibrium hypersonic flow over a rectangular cavity embedded on a compression ramp

This paper reports a systematic computational investigation that elucidates the fundamental thermochemical non-equilibrium physics that occurs when air at Mach number of 11 encounters a rectangular cavity of aspect ratio L/D = 2.0 embedded on a 25° compression ramp. The mechanistic details of this highly complex flow phenomenon are obtained by solving the compressible form of the Navier–Stokes equations in two dimensions using a finite-volume open-source library. Chemical and thermal non-equilibrium processes are treated using a five-species, 12-reaction chemical kinetics, and a two-temperature model, respectively. Following a detailed validation and grid sensitivity study, two simulations are conducted, one with isothermal boundary conditions and the other with conjugate heat transfer (CHT) to identify the effect of energy transmission to the material on surface heat flux. Fast Fourier transforms and near-wall velocity profiles inside and in the neighborhood of the cavity are used to identify primary oscillatory modes and shear layer dynamics. Two new descriptive states defined as “states I and II,” representative of the minimum and maximum deflection of the shear layer, are used to discuss the dynamical behaviors in the cavity, including the separation region before the cavity, trailing edge effects, frequency analysis of probe data collected at several key locations, and the effect of CHT on surface heat flux. It is found that the flow features at the cavity's center strongly influence the separation upstream of the cavity, and the transrotational temperature near the cavity's trailing edge is strongly correlated with the oscillations of the shear layer.

Thermal Hydraulics Simulation of a Water Spray System for a Cooling Fluid Catalytic Cracking (FCC) Regenerator

Olefins are crucial building blocks for petrochemical industry, serving as raw materials for the production of various products such as plastics, synthetic fibers, detergents, solvents, and other chemicals. In FCC, heavy petroleum feedstocks are injected into a catalytic cracking unit, where they are mixed with a catalyst. The catalyst aids in breaking down the large hydrocarbon molecules into smaller fragments, including olefins like Propylene and Ethylene. These polymerization reactions occur at high temperatures. They demand that heat removal occurs as quickly as possible in order to control the reactor temperature and to avoid “hot spots” in the Regenerator or localized oxidation reactions (and to avoid creep rupture of the regenerator steel cladding). The cooling of the regenerator cladding surface can be achieved by impinging water droplets (spray), ejected from a spray nozzle. Spray cooling can provide uniform cooling and handle high heat fluxes in both a single phase and two phases. This research provides a thermal hydraulic design of regenerator spray cooling systems. In the framework of this research, Fire Dynamics Simulator (FDS) software was applied in order to simulate the temperature field and the water vapor mass fraction. A COMSOL Multiphysics finite element code was used in order to calculate the temperature field inside the regenerator cladding. The calculated surface temperatures and heat transfer convective coefficient, obtained using FDS software, were validated successfully against COMSOL numerical results and previous results in the literature. The numerical simulations were carried out for two cases. The first case was carried out at a distance of 0.5 m, and the second case was carried out at a distance of 0.2 m. A grid sensitivity study was carried out on the FDS model. Numerical integrations were carried out over time in order to calculate the average temperatures. The difference between these four average temperatures, calculated by applying different grids, is less than 7.4%. The calculated surface temperatures and heat transfer convective coefficient were validated successfully against COMSOL numerical results and previous research. It was shown that the calculated temperatures decrease in the second case. The water spray system managed to cool the steel wall more effectively as the water spray system approaches the steel cladding.

Development of an Interpolated RANS-LES solver for wall-bounded turbulent fluid flows

Scale-resolving hybrid RANS-LES models are being increasingly used to numerically simulate wall-bounded turbulent flows since they are computationally efficient compared to LES (Large Eddy Simulation) while also being more accurate compared to RANS (Reynolds-Averaged Navier–Stokes) models. Existing hybrid RANS-LES models typically either suffer from modeled stress depletion or they can introduce a sharp jump in the mean Reynolds stress across the RANS-LES zone. We present a novel interpolated RANS-LES (IRL) solver implemented in OpenFOAM, in which the RANS equations (for eddy-viscosity) and LES equations (for filtered velocity) are evolved simultaneously on the same computational grid. Based on the distance from the no-slip wall, the flow domain is demarcated into a near-wall, free-stream, and a “hybrid” region, sandwiched between the near-wall and free-stream regions. Spalart Allmaras (S-A) model has been used to evolve the RANS eddy viscosity, while the Wall Adapting Local Eddy Viscosity (WALE) model is used for evolving the filtered LES velocity. Using a novel scheme based on solving a partial differential equation, the turbulent eddy-viscosity is interpolated between the RANS eddy viscosity in the near wall region and the effective eddy viscosity from the resolved LES fluctuations in the free-stream region. The mean subgrid stress in the LES momentum equation is then corrected in the near-wall and hybrid regions using the interpolated turbulent eddy-viscosity. IRL’s grid resolution requirements are the same as DES’s (Detached Eddy Simulation). Grid convergence and sensitivity studies have been carried out for two canonical flow geometries. For turbulent channel flow simulations, the IRL solver does not display the typical mean stress depletion observed in DES while also predicting resolved Reynolds stresses accurately in the free-stream region. For simulations of flow over backward facing step, the IRL solver can predict the friction factor and coefficient of pressure better than that DES, especially after the reattachment point.

Analysis of shift in discharge coefficient for contaminated multihole orifice flow meter

The aim of this paper is to provide a comprehensive numerical study of influence of contamination and multi-hole orifice (MHO) flow meter geometry on shift in discharge coefficient. In this research, authors analyzed steady, three-dimensional, turbulent flow of dry air through MHO flow meter with 3 different β parameters. Shift in discharge coefficient was calculated for 29 different combinations of contamination parameters for 7 contamination angles, 3 different Reynolds numbers and 3 β parameters. Numerical method that was used was finite volume method with SIMPLE algorithm and standard k-ε turbulence model. Grid sensitivity study was performed on four systematically refined numerical grids for MHO with contamination for all values of Reynolds number, all values of β parameters and 3 values of contamination angle. Obtained results were grid independent. Numerical results for MHO without contamination were compared with the experimental results found in the literature. It was found that contamination angle has most influence on shift in discharge coefficient. Also it was found that the contamination has influence on the change of pressure drop values, which directly affects the change of other parameters. Pressure drop and singular pressure loss coefficient of the orifice with contamination are smaller compared to the values for a pure orifice, whereby the measurement accuracy was reduced. It was found that the multi-hole orifice meter was less sensitive to the pressure drop changes due to the increasing of contamination angle in regards to the single-hole orifice meters.

Computational Study of a Swirl Stabilized Flame in a Gas Turbine Combustor

Size of recirculation zone and temperature field is examined in a swirl stabilized aero gas turbine combustor using computational fluid dynamics. A sector of 22.5° of an annular combustor is modeled for the study. Unstructured tetrahedral meshes comprising 1.2x106 elements are employed in the model where the governing equations are solved using CFD flow solver CFX. Grid sensitivity study and code validation through limited test data has also been carried out. The effect of pressure and fuel-air ratio on the size of recirculation zone and apparent size of flame zone has been studied at steady state conditions. The heat release rate due to combustion is found to affect the length of recirculation zone as well as the velocity field inside the combustion liner.

Impact of eave and roof pitch on cross ventilation for an isolated building with sawtooth roof

An eave refers to an extension attached to the building roof to protect the interior space from direct solar radiation and improve the performance on cross ventilation. In this study, the impact of eave inclination angle and roof pitch of an isolated sawtooth roof building on cross ventilation were investigated. The eave configurations at either windward or leeward openings were included. 3D steady Reynolds-Averaged Navier-Stokes (RANS) equation in combination with the Shear-Stress Transport model (SST k-ω model) was used for the Computational Fluid Dynamics (CFD) simulations. Grid sensitivity study was carried out and the performance of cross ventilation was evaluated based on the non-dimensional velocity magnitude, spatial distribution of pressure coefficient as well as the ventilation rate of the building. For the simulation model with 55° roof pitch, it is observed that a region with high velocity magnitude formed on top of the leeward eave due to the higher roof pitch and presence of the leeward eave. Results also indicated that the building model with 90° leeward eave and 55° roof pitch has the highest increment in ventilation rate which is 7.16%. On the other hand, the building model with 90° windward eave has the highest pressure coefficient because more blockage of airflow is caused by a steeper roof as the roof pitch of the building increases. Furthermore, the building model with 90° leeward eave shows a larger region with negative pressure at the leeward façade indicating higher airflow leaving the leeward opening. Therefore, the airflow behavior and characteristic are both dependent on the roof pitch and eave inclination angle for a naturally ventilated building.

Numerical Analysis of NACA 6409 and Eppler 423 Airfoils

The present study aims to numerically examine two commonly used airfoils that provide high lift at low Reynolds numbers, using various turbulence models with the commercial code ANSYS Fluent 19. Eppler 423 airfoil is widely used in wind turbine blades, wings of gliders and R/C airplanes designed for sport flying and competitions thanks to high lift capability. NACA 6409 airfoil is seen in free-flight model planes. In literature review, it has been found that k-ω SST, k-kl-ω and 𝛾-𝑅𝑒𝜃 transition models have a good prediction in Cl-Cd estimation of drag bucket region, transition flow regime, and laminar separation point. In addition to RANS models, laminar model was also investigated. In order to ensure the accuracy and precision of numerical investigations; the grid sensitivity study and near wall model were used. The method of decomposition of flow direction vector procedure has been used to provide the advantage of simulating all angle of attacks using only a single mesh. The lift coefficient curves and polar diagrams were formed from the angle of attack of -4 degree up to the stall. Experimental data, XFoil analysis data and numerical analysis data were evaluated in terms of best glide ratio-(Cl/Cd)max and minimum sink-(Cl1.5/Cd)max criteria. As expected, at high angle of attack, the transition ramp on the upper surface moves forward and the pressure gradient became more adverse. The results showed that fine mesh gives better results than other mesh sizes and transitional γ-Reθ SST model was successful in modelling the transition event than other turbulence models.

Large-eddy simulation of a planar offset-jet with heat transfer: The effects of ventilation

Several industrial and engineering applications employ ventilated offset-jets such as flow separation control devices, upper surface blowing used in short take-off and landing aircraft, drying processes, and fuel injection systems. The design and efficient operation of these devices relies on controlling the turbulence levels, heat transfer rates, and the locations of flow reattachment points. Therefore, an understanding of the flow and thermal characteristics of ventilated offset-jets helps in the design and operation of these devices. The present study aims to address this by performing large-eddy simulations (LES) of a planar turbulent offset-jet for different velocities of the ventilated-jet. To understand and quantify the effects of ventilation on the offset-jet, a range of velocity-ratios of 0, 0.1, 0.14, and 0.18 are considered for a primary offset-jet Reynolds number of 14,000. First, a grid-sensitivity study is performed to establish the convergence of the solver and the LES index of quality of resolution is obtained to ascertain the quality of the mesh. Thereafter, the mean and second-order statistics of the streamwise component of velocity are compared with reference data from the literature to validate the numerical solver. The effects of ventilation on the offset-jet are studied using the decay of streamwise velocity, jet-spread, the evolution of pressure coefficient, friction coefficient, the mean Nusselt number, and the unsteady characteristics. As the velocity-ratio increases, the suction pressure in the initial region of the domain decreases, which causes the jet to spread away from the bottom wall. Further the peak magnitudes of the mean Nusselt number and coefficient of pressure decrease exponentially with an increase in the velocity-ratio. The peak magnitude of the mean Nusselt number for the unventilated-jet case – the offset-jet without the ventilated-jet – is the largest and its value for the ventilated-jet cases decreases as the velocity-ratio is increased. It is concluded that in the initial region, the ventilated-jet has a profound effect on the offset-jet flow and as the flow develops into a wall-jet the effects of ventilation vanish.

Progress in the validation of rotor aerodynamic codes using field data

Abstract. Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

Numerical simulations of high Prandtl number liquid jets impinging on a flat plate

In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226<Re<2850, 77<Pr<161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.

Experimental and numerical investigation of beer foam

The frothy foam on top of beer is produced by bubbles of gas, predominantly carbon dioxide, rising to the surface. The chemical components that produce the head are wort protein, yeast, and hop residue. This involves a large number of chemical components and numerous physical interactions. The beer foam and, especially, its stability is an essential quality characteristic of a beer. A consumer defines a beer's head by its stability, quantity, lacing (glass adhesion or cling), whiteness, creaminess, and strength. To achieve an accurate prediction of beer foam formation and collapse is challenging because complex numerical models are required to account for these nonlinear beer foam effects. To analyze a new design of a beer bottom-to-top pouring system for the startup company Einstein 1, we first set up the experimental tests of this pouring system. Afterward, we performed the associated repeatability studies to achieve stable working conditions. To study beer foam formation and its collapse, we employed a multiphase Reynolds-averaged Navier–Stokes solver that considered two inter-penetrating continua, which allowed accounting for multi-component phases and mass and heat transfer between these phases. We numerically and experimentally investigated beer foam patterns, beer heights, beer/foam ratio, foam height, foam stability, and foam volume fractions. We performed grid sensitivity studies and validated the numerical solver by comparing results against model test measurements. The results indicated that that the higher the temperature of the beer and the higher the tap pressure is, the greater the foam development and the associated foam height are but not necessarily the foam stability.

Aerodynamic Performance Enhancement of an Axisymmetric Deflector Applied to Savonius Wind Turbine Using Novel Transient 3D CFD Simulation Techniques

Many recent studies show that the performance of Savonius turbines can be considerably increased by using wind deflectors. Axisymmetric deflectors are particularly interesting; they concentrate the wind flow in all directions. This study aims to aerodynamically optimize the truncated cone deflector shape through transient 3D CFD simulations using sliding mesh techniques. To reduce the mesh size and thus the simulation time, symmetrical boundary conditions were applied to rotating body faces. A mesh grid sensitivity study was conducted to define the optimum mesh size. Additionally, hybrid numerical approaches combining coupled and SIMPLE solvers were particularly influential in reducing computational time. Concave- and convex-arced-shaped faces deflectors were compared to the original truncated cone deflector, showing an increase in the performance for the convex type and a decrease for the concave one. Then, eight cases involving convex spline shape deflectors were simulated. All these deflectors had an equal volume to the original truncated cone deflector. One of the cases showed a 20% average increase in the performance over the original deflector. This result shows the importance of the geometrical shapes in the design of axisymmetric deflectors.

Numerical assessment of the effects of a porous medium in transient sloshing-induced impulsive hydrodynamics

A compressible VoF solver was employed to predict sloshing-induced flow phenomena with and without accounting for the effect of a porous medium situated in various zones inside partially filled tanks. Grid sensitivity studies ensured obtaining grid-independent solutions. The solver, validated against experimental measurements, was employed to systematically scrutinize the influence of positions and shapes of different porous zones as well as their number on sloshing-induced free surface elevations and associated impulsive impact pressures. Results indicated that 1) a porous medium situated at the bottom of a tank increased impulsive impact pressures; 2) sufficiently high thin single porous zones at the tank's bottom and sides reduced impulsive impact pressures; 3) placing several porous zones decreased free surface elevations and the associated impulsive impact pressure; 4) although small porous zones at the sides of a tank dampened sloshing phenomena, their position was more effective than their number; and 5) the effectiveness of porous zones to dampen sloshing was similar to that of solid baffles.

Large-Eddy Simulations of Wind-Driven Cross Ventilation, Part1: Validation and Sensitivity Study

Natural ventilation is gaining popularity in response to an increasing demand for a sustainable and healthy built environment, but the design of a naturally ventilated building can be challenging due to the inherent variability in the operating conditions that determine the natural ventilation flow. Large-eddy simulations (LES) have significant potential as an analysis method for natural ventilation flow, since they can provide an accurate prediction of turbulent flow at any location in the computational domain. However, the simulations can be computationally expensive, and few validation and sensitivity studies with respect to simulation parameters such as grid resolution and boundary conditions have been reported. The objectives of this study are to validate LES of wind-driven cross-ventilation and to quantify the sensitivity of the solution to the grid resolution and the inflow boundary conditions. We perform LES for an isolated building with two openings, using three different grid resolutions and two different inflow conditions with varying turbulence intensities. Predictions of the ventilation rate are compared to a reference wind-tunnel experiment available from literature, and we also quantify the age of air and ventilation efficiency. For the cross-ventilation case modeled in this paper, the prediction of the mean ventilation flow rate is very robust, showing negligible sensitivity to the grid resolution or the inflow characteristics with the maximum error of 1.9 and 1.3% for each sensitivity study. However, a sufficiently fine grid resolution is required to obtain accurate predictions of the detailed flow pattern and the age of air as they show comparably larger errors of 10 and 20% in the grid sensitivity study, and the standard deviation of the instantaneous ventilation rate is affected by the turbulence in the inflow condition with showing 44% difference.

Investigations of transient sloshing induced impulsive hydrodynamics

Sloshing-induced impulsive pressures may be extremely high and of short rise time duration. A compressible Volume of Fluid (VoF) method was adopted to study these pressures in a partially filled three-dimensional tank with various oblique tank ceiling angles. Comparative available experimental model test measurements of free surface profiles and velocity and pressure time histories validated this numerical method. To obtain a sufficiently fine grid to yield converging predictions, grid sensitivity studies were performed. The influence of the oblique tank ceiling angle on sloshing-induced impulsive pressures was examined. Results demonstrated that the VoF method obtained accurate predictions of sloshing phenomena. The oblique tank ceiling angle dampened impulsive pressure and, therefore, this angle had to be accounted for in the design of the tank.

Study on aerodynamic features of rod thrust vector control for physical applications

Mechanical thrust vector control is a classical and important branch in the vectoring control field, offering an extremely reliable control effect. In this article, a simple technology using a cylindrical rod has been numerically investigated to achieve jet controls for three-dimensional conical axisymmetric nozzles. Complex flow phenomena caused by the cylindrical rod on a flat plate and in a converging–diverging nozzle are elucidated with the purpose of a profound understanding of this technique for physical applications. Published experimental data are used to validate the dependability of current CFD results. A grid sensitivity study is carried through and analyzed. The result section discusses the impacts of three factors on performance, involving the rod penetration height, rod location, and nozzle pressure ratio. Significant vectoring performance variations and flow topologies descriptions are illuminated in full detail. When the rod penetration height changes, this technique has an effective control range, namely H/Rt ≤ 0.4. In this effective control range, the vectoring angle and efficiency increase and the thrust coefficient decreases with a deeper rod insertion. As the rod location moves downstream towards the nozzle exit, the vectoring angle increases and the thrust coefficient decays. Moreover, the direction of jet deflection remarkably varies for diverse rod locations. While the nozzle pressure ratio increases, the vectoring angle initially increases to reach the maximum level and then decays slightly. Meanwhile, the thrust coefficient continuously increases.

Hydrodynamic characteristic investigations of a collaborative spoiled rudder under open water conditions

This study presents the hydrodynamic performance estimation of a collaborative spoiled rudder (CSR) - a new concept of ship rudder under open water conditions. The grid sensitivity study is carried out and the numerical method used in this work is validated by comparing with experimental data. The hydrodynamic performances of different CSR configurations at various rudder angles are simulated. The result shows that the existence of the spoiler has positive effect on the rudder lift and the effect gets more obvious with the increase of spoiler open angle. The spoiler can be controlled easily owing to the small spoiler stock torque. The optimal matching results show that, under the optimal matching conditions, the lift to drag ratios of the CSR are generally higher than those of the CR for the same target lift, and the improvement of the CSR in the rudder effectiveness gets more obvious with the increase of target lift. Furthermore, the CSR having the optimal matching rudder angle and spoiler open angle is much easier to steer compared with the CR with the same lift being generated. Overall, properly adjusting the rudder angle and the spoiler open angle of the CSR can significantly improve the rudder effectiveness.

Large-eddy simulation of a planar offset wall-jet with heat transfer: Characterization, turbulent kinetic energy and Reynolds shear stress budgets

Planar wall-attaching offset jets can be found in a variety of engineering applications. In the present work, a planar turbulent offset wall-jet, at a Reynolds number of 7500 and an offset-ratio of 5, with heat transfer is numerically studied using large-eddy simulation (LES). First, a grid sensitivity study is performed and the quality of the mesh used is assessed using LES index of quality of resolution and an appropriate mesh is selected. Next, the solver is validated for mean streamwise velocity, streamwise Reynolds normal stress, and the evolution of coefficient of pressure on the wall, using reference experimental data from the literature. Thereafter, the unsteady flow features of the jet, maximum velocity decay and jet-spread, statistical features of flow and temperature, and turbulent kinetic energy (TKE) budgets and in-plane component of Reynolds shear stress are studied. The streamwise evolution of Nusselt number on the wall shows three distinct peaks and the location of these peaks correlates with a change-of-sign of the wall skin-friction coefficient. The domain can be categorized into three distinct regions namely, the recirculation, impingement and wall-jet regions. Distinct flow and thermal characteristics are observed depending on the region of interest. The peak magnitudes of Reynolds normal and shear stresses increase as one moves from the recirculation region to the impingement region and decreases thereafter attaining a lowest value in the wall-jet region. The heat transfer from the wall is observed to be most effective in the recirculation region. The production of TKE is dominant in all the three regions of the flow, whereas the Reynolds shear stress budgets indicate that the production and velocity–pressure-gradient terms are dominant and balance each other in all the three regions of the flow.