Turbulent Transfer Research Articles

Two-point statistics in non-homogeneous rotating turbulence

Time-resolved two-dimensional two-component particle image velocimetry measurements with high spatial resolution are carried out in a water tank agitated by four blades rotating at constant speed. Different blade geometries and rotation speeds are tested for the purpose of modifying turbulent flow conditions. In all cases where no baffles are used to break the rotation, the Zeman length is an order of magnitude smaller than the Taylor length. Compared with the cases with baffles which break the rotation, in the unbaffled cases turbulence production and/or mean advection are significant and the turbulence nonlinearity is dramatically reduced for the horizontal (i.e. normal to the axis of rotation) two-point turbulence fluctuating velocities. This nonlinearity reduction is manifest not only in the interscale turbulent energy transfer but also in the interspace turbulent energy transfer, which nearly vanishes. However, the nonlinearity is not reduced for the vertical two-point turbulence fluctuating velocities: the corresponding interscale turbulent transfer rate is in fact intensified, and its dependence on the two-point separation distance, as well as that of the corresponding interspace turbulent transfer rate, which does not vanish, is significantly modified. Even though non-homogeneities are very different for different blades and rotation speeds in the unbaffled cases, the horizontal fluctuating velocity's second-order structure function collapses with scalings which resemble predictions for homogeneous turbulence subject to strong rotation. The vertical fluctuating velocity's second-order structure function does not collapse for different blade geometries by neither these nor the Kolmogorov predictions.

THE IMPACT OF SOLAR RADIATION ON THE TEMPERATURE REGIME OF A ROOM IN WINTER

The article is devoted to the computational studies of the air-thermal state of office premises taking into account the effect of solar radiation coming through window openings. The study was conducted for a room with two windows using two heating devices installed under them. The air-temperature regime of premises in winter, characterized by the highest thermal energy consumption for heat supply, is considered. The study is based on the solution of a three-dimensional nonlinear heat transfer problem described by a system of equations of turbulent momentum and energy transfer. The k-e turbulence model is used to close this system. The results of numarical modeling of the physical situation under study are presented. The research data on the features of the air-thermal state of the premises under solar radiation conditions are given. The results of a comparative analysis of the air-temperature conditions of office premises corresponding to the solution of the specified heat exchange problems in the presence and absence of solar radiation are presented and the effects of solar radiation on the structure of the air flow and the thermal state of the premises are established. It is shown that in the presence of solar radiation, the air flow picture and the character of the temperature fields in the premises change significantly. In particular, in these conditions, the increase in the average temperature of the premises for the studied period is 2.5 °С. The possibility of a certain reduction of the load on the heating system in the presence of solar radiation is noted. Bibl. 17, Fig. 5.

Numerical investigation of turbulent mass transfer processes in turbulent fluidized bed by computational mass transfer

Turbulent fluidized bed possesses a distinct advantage over bubbling fluidized bed in high solids contact efficiency and thus exerts great potential in applications to many industrial processes. Simulation for fluidization of fluid catalytic cracking (FCC) particles and the catalytic reaction of ozone decomposition in turbulent fluidized bed is conducted using the Eulerian–Eulerian approach, where the recently developed two-equation turbulent (TET) model is introduced to describe the turbulent mass diffusion. The energy minimization multi-scale (EMMS) drag model and the kinetic theory of granular flow (KTGF) are adopted to describe gas–particles interaction and particle–particle interaction respectively. The TET model features the rigorous closure for the turbulent mass transfer equations and thus enables more reliable simulation. With this model, distributions of ozone concentration and gas–particles two-phase velocity as well as volume fraction are obtained and compared against experimental data. The average absolute relative deviation for the simulated ozone concentration is 9.67% which confirms the validity of the proposed model. Moreover, it is found that the transition velocity from bubbling fluidization to turbulent fluidization for FCC particles is about 0.5 m·s–1 which is consistent with experimental observation.

Polymer-doped two-dimensional turbulent flow to study the transition from Newtonian turbulence to elastic instability

Detecting the flow regimes of Newtonian turbulence (NT), elasto-inertial filament (EIF), elasto-inertial turbulence (EIT), and maximum drag reduction (MDR) of polymer solution and their transition has been a hot topic in the last decade. We attempted to detect NT, EIF, EIT, and MDR by visualizing vortex shedding downstream of an array of cylinders that was inserted perpendicular to polymer-doped two-dimensional (2D) flow. Since polymers are stretched at the cylinders, the consequent vortex shedding is affected by viscoelasticity. The flow regimes are characterized based on Weissenberg (Wi) and Reynolds numbers (Re) with the relaxation time of the polymeric solution estimated from capillary-thinning experiments. The flow regimes are observed for different molecular weights of polyethylene oxide and polyacrylamide in solution and are categorized as either vortex type 1, type 2, and type 3 on a Re–Wi map based on flow visualization using particle image velocimetry. In addition, turbulent statistics of these flow regimes are calculated to more fully quantify these flow regimes. We found that vortex types from 1 to 3 have a similarity to NT, EIF, EIT, and MDR. In addition, characteristic turbulent energy transfer without an increase in turbulent energy production was found in the flow regimes of vortex types 2 and 3 of each polymer solution. Our results suggest intriguing parallels between pipe, jet, and 2D turbulent flow for drag-reducing polymeric solutions.

Contemporary and Conventional Passive Methods of Intensifying Convective Heat Transfer—A Review

The ever-increasing demand for effective heat dissipation and temperature control in industrial and everyday applications highlights a critical research problem. The need for development is not only in terms of providing thermal comfort to humans but also forms the basis for the efficient operation of machines and equipment. Cooling of industrial machinery and household electronic equipment is a crucial element in any manufacturing process, and the planning and design of appropriate cooling systems continues to be an integral part of the machine design and construction process. Manufacturers aim to maximize performance while minimizing size and weight. This article reviews widely used passive methods to enhance heat transfer, focusing on their effectiveness in improving convective heat transfer. The techniques examined include surface modifications and advanced materials like foamed metals and nanostructured coatings, which influence turbulence and heat transfer coefficients. The key findings demonstrate that surface roughness, perforated fins, and twisted tapes enhance fluid mixing but may increase flow resistance. The review underscores the significance of these passive methods in optimizing cooling system efficiency across various applications. Despite the variety of techniques available, many areas, especially those involving laser beam modifications, remain underexplored, indicating a need for further research in this field.

Characterization of turbulent multiphase particle mass transfer in gas-liquid flows

Turbulent mass transfer is crucial in the design and operation of industrial multiphase gas-liquid flow reactors. Despite its ubiquity, due to the complexity of multiphase turbulent flows, there is a lack of computational models and experimental data on mass transfer fluxes at the local scale for non-ideal bubbles. We develop a new measuring technique to visualize a well-known multiphase transfer problem: the local retention of water-soluble aerosols in large rising cap bubbles. The hydrodynamics of the bubble and the liquid phase electrolyte concentration of the water-soluble aerosols are resolved with high spatial and temporal resolution via conductivity measurements using a high-resolution Wire-mesh sensor of 1mm pitch with an acquisition frequency of 3200Hz. We obtain a representation of prototypical cap bubbles of 8.25mm equivalent spherical volume bubble diameter, and the associated wakes containing the dissolved aerosol material based on reconstruction from the available large ensemble of bubbles. These high-fidelity grade data can be directly used to develop and validate particle mass transfer models for heterogeneous multiphase flows. At the multiphase flow regime in this work, the particle mass transfer is governed by inertial effects in the gas phase and by turbulent diffusivity and interface-related oscillations in the liquid phase. We formulate a physical model of the phenomenon that considers turbulence's stochastic nature via the Reynolds decomposition of the investigated variables. Our model supports the interpretation of the measured data. It also enables the calculation of the particle mass transfer flux directly from the concentration field in the mass transfer wake. The predictions and assumptions of our model are validated against the experimental measurements. This work will be extended to resolve the transient behavior of the aerosol retention, and to calculate the associated particle mass transfer coefficients.

A 1+ Mechanism Model for Predicting the Mixed-Oil Concentration in Multiproduct Pipelines

Summary Petroleum products are frequently transported successively through the same multiproduct pipeline. Due to turbulent and convective diffusion mass transfer, two adjacent oils will mix with each other, forming a mixed-oil segment. Accurate and rapid prediction of mixed-oil concentration is crucial for the precise management of mixed-oil segments. Conventional 1D modeling methods exhibit shortcomings in accurately representing the asymmetric distribution characteristics of mixed-oil concentration curves, and high-dimensional models are not practically applied due to their prohibitive computational time costs. Building on the 1D model framework, this paper proposes a “1+” mechanism model by considering the convective mass transfer behavior between the turbulent core region and the laminar boundary layer, and new governing equations and corresponding numerical solution methods are also introduced. Simulation experiments affirm the ability of the new model to characterize the asymmetric distribution features of mixed-oil concentration curves, along with its high computational efficiency in engineering applications. This is demonstrated by the computational time of approximately 30 seconds for simulating a pipeline of 300 km in length (Δx = 10 m, Δt = 1 second, CPU: i5-12500H, RAM: 16 GB). When applied to pipelines in industrial scenarios, the new model is shown to accurately predict the distribution of mixed-oil concentration curves. The research findings are significantly beneficial in assisting field personnel to gain advanced insights into the mixed-oil concentration distribution at the station, enabling timely and well-informed strategies for handling mixed-oil segment, thereby enhancing the operational efficiency of multiproduct pipelines.

Numerical investigation of natural convection from a horizontal heat sink with an array of rectangular fins

This research focuses on the numerical analysis of a rectangular fin array with a horizontal heat sink to examine its performance in natural convection heat transfer. The study uses air as the primary working medium. Steady-state 3D modeling is carried out using the finite volume method (FVM), allowing velocity and temperature distributions to be evaluated by solving equations relating to mass conservation, fluid dynamics, turbulence, and heat transfer. This research details temperature and velocity variations as well as fluid trajectories. It also explores how the intensification of heat flow, varying from 361 W/m2 to 2527 W/m2, influences the cooling efficiency of the fins. It is observed that the increase in heat flux reduces the formation of vortices and intensifies natural convection by rising the temperature gradient between the radiator and the surrounding. The temperature increases at all parts of the heat sink fin array as heat flux increases. Besides, it is noticed that the maximum temperature is generated at the mid-region of the base plate heat sink and gradually dissipates to the surroundings through the bodies of the fins. Moreover, the study indicates that the average convection heat transfer coefficient (hav) and the average Nusselt number (Nus) increase by 86 % and 84 %, respectively, when the heat flux is increased from 361 W/m2 to 2527 W/m2. These observations are corroborated by a comparison with existing experimental data, showing an agreement of less than 9 %. The digital model is validated by comparing pre-existing data on radiators with horizontal rectangular louvers, revealing minimal deviations of less than 2 %.

Modeling of turbulence, chemistry, and heat transfer of a rocket nozzle

Calculating the composition of the combustion products serves as the starting point for the high-area ratio (HAR) rocket nozzle. The HAR rocket nozzle finds application in next-generation hypersonic vehicles and missiles. The gas dynamics, gas radiation, and turbulence-chemistry interactions (TCI) play a crucial role in simulations of hypersonic turbulent reactive flows through a rocket nozzle. The interaction of turbulence and chemistry leads to random species mixing and enhanced heat transfer. The k-omega shear stress transport (SST) model is employed to calculate the average flow fields. The eddy dissipation concept (EDC) with reduced chemistry is applied to capture the TCI in the nozzle. The spectral properties of H2O are calculated using the full spectrum k-distribution (FSK) model in conjunction with the first-order spherical harmonics method. The mass fractions of combustion products obtained from Particle Swarm Optimization (PSO) and Monte Carlo (MC) methods agree well with the Chemical Equilibrium with Applications (CEA) code. The flow fields obtained with the current solver are validated with the published experimental data. A significant rise and fall of 7–8% have been observed in the centerline flow properties when not including the TCI models in the flow solver. The centerline flow fields obtained with different discretization schemes are in good agreement. The radiative heat flux from the H2O species is more dominant than the convective.

Renormalization group model of turbulent bioconvection of gyrotactic microorganisms

The problems connected with processes of turbulent bioconvection are considered. These processes are typical for technologies of microbiological and other industries. The concept of using renormalization groups in the theory of turbulence is shown. A step-by-step algorithm for applying the renormalization methodology to deriving the equation for mass transfer in bioconvective processes, which includes renormalization of the turbulent diffusion coefficient, is presented. The dependence for the turbulent Schmidt number was obtained. It has been shown that turbulent Schmidt number is function of six parameters: molecular Schmidt number, turbulent viscosity, molecular viscosity, parameter of geometric shape of microorganisms, average velocity of motion of microorganisms relative to the fluid, orientation parameter. The turbulent Schmidt number enables closing the diffusion equation. The effect of the geometric shape of microorganisms on the development of turbulence in a bioconvective flow in the presence of gyrotactic microorganisms is analyzed. Analysis shows that microorganisms with spherical symmetry make a small contribution to the formation of turbulent mass transfer processes. The results of the analysis are in a good agreement with data from studies of the nonlinear hydrodynamic stability of bioconvective processes based on the Lorentz model.

A turbulent mass diffusivity model for analyzing the mixing characteristics in an impinging stream-rotating packed bed

In this study, the fluid flow and mixing process in an impinging stream-rotating packed bed (IS-RPB) is simulated by using a new three-dimensional computational fluid dynamics model. Specifically, the gas–liquid flow is simulated by the Euler–Euler model, the hydrodynamics of the reactor is predicted by the RNG k–ε method, and the high gravity environment is simulated by the sliding mesh model. The turbulent mass transfer process is characterized by the concentration variance and its dissipation rate εc formulations, and therefore the turbulent mass diffusivity can be directly obtained. The simulated segregation index Xs are in agreement with our previous experimental results. The simulated results reveal that the fringe effect of IS can be offset by the end effect at the inner radius of RPB, so the investigation for the coupling mechanism between IS and RPB is critical to intensify the mixing process in IS-RPB.

Uniturbulence statistics and analysis of factors influencing the energy spectrum

This paper investigates the dynamics of unidirectionally propagating surface Alfvén waves, employing magnetohydrodynamic numerical simulations and statistical methodologies. The primary goal of this work is to enhance our understanding of the nonlinear self-cascade of surface Alfvén waves, which we term as uniturbulence, by unraveling the complex relationships between various length scales and their interplay with turbulent energy transfer mechanisms. To achieve this, we extensively analyze the phenomenon of uniturbulence using methods such as power spectrum analysis, radially averaged Fourier transform, and kurtosis. We employ these techniques to investigate the spatiotemporal distributions of kinetic and magnetic energy in uniturbulent flows. We also reveal the crucial role of the density contrast's variations and the role of Yaglom's law in characterizing energy transfer mechanisms. Our findings reveal that the inertial range of the perpendicular kinetic energy and magnetic energy along the z-axis depicts a progressive change in slope values, ultimately approaching the often-observed values of −5/3 and −3/2, respectively. Furthermore, our kurtosis analysis highlights the non-Gaussian behavior of the flow field at different length scales and over time, offering a perspective on uniturbulence dynamics. The correlations observed among diverse statistical approaches emphasize the complex interplay between different length scales in the context of uniturbulence. Our findings contribute to understanding this phenomenon, establishing a basis for future investigations to clarify the connections regulating these turbulent dynamics.

Numerical and experimental analysis of a cross-finned solar receiver for parabolic dish collectors

A parabolic dish solar collector (PDC) is used to concentrate solar radiation on a receiver to produce hot water or steam for various thermal applications. The solar receiver of the PDC is one of the significant components to improve the overall efficiency of the PDC. The present study is a numerical and experimental investigation of the thermal performance of a solar receiver with vertical fins with and without cross fins. The fin orientation was investigated using 0.056 kg/s, 0.083 kg/s, and 0.111 kg/s water flow rates. COMSOL Multiphysics® was used to simulate the receiver to determine the receiver’s temperature behavior. The thermal performance parameters such as heat gain, heat loss, convective heat transfer coefficient, and thermal, and exergy efficiency were analyzed under various flow conditions. The cross fins effectively prolong the water residence, increasing turbulence, and heat transfer and reducing heat loss to the surroundings. The receiver with cross fins at higher flow rate conditions shows higher thermal performance than lower flow rate and without fins conditions. The solar radiation during the outdoor testing is in a range of 300–950 W/m2 and the ambient temperature varies from 30-37 °C. The peak average thermal and exergy efficiency obtained for the cross-finned receiver at 0.111 kg/s is 64.4 %, and 6.3 %, respectively. Cross-finned receiver operating at 0.111 kg/s reduced CO2 mitigation rate, estimated at approximately 0.29 tons of CO2 per year, along with a lower carbon credit gain at $0.47. Thus, the proposed cross-fin solar receiver demonstrates improved thermoeconomic performance than conventional receiver designs, offering a more cost-effective design with minimal investment.

Numerical study of structural optimization for improving gas-liquid flow and microalgae carbon fixation rate in air-lift photobioreactors

In this paper, the gas-liquid flow pattern and CO2 dissolution transport process of microalgae solution in an air-lift photobioreactor (PBR) are investigated by numerical simulations. In order to promote the gas-liquid mass transfer efficiency to enhance the carbon fixation, we propose different optimization schemes for the structure. Firstly, a two-phase flow model and a mass transfer model are employed to investigate the fluid flow characteristics. The flow state is evaluated by velocity in light direction, turbulent kinetic energy (TKE), gas hold-up and CO2 mass transfer coefficient. The results show that the gas flow rate of 0.2 V/V/min is the optimum value for the inlet. Secondly, CO2 mass transfer and carbon fixation rates are analyzed for different inlet CO2 volume fractions. With the increase of CO2 volume fraction, the carbon fixation rate first increases and then decreases, the highest fixation rate reaching 2.21 × 10−4 mol m−3 s−1 at 6 %. Finally, based on numerical simulations, three optimized structures with different spacing of baffles are proposed and evaluated the performance in comparison with the original. The results show that the gas distribution and CO2 mass transfer are significantly improved, with the one-sided baffle structure is the best. In summary, this study provides a new and useful reference to evaluate the gas-liquid mixing and improve the microalgae carbon fixation rate.

Numerical investigation into turbulent drag reduction via the application of pufferfish spine-inspired cone microstructures in Suboff models

The effective reduction of seawater drag is pivotal in enhancing the speed and minimizing the energy consumption of submarines, which has significant implications in the fields of energy and defense. Surface bionics has emerged as one of the leading techniques for drag reduction. Current research primarily focuses on replicating the groove-like structures observed on shark skins and the flexible properties of dolphin skins. However, the application of cone microstructures on submarine surfaces remains relatively underexplored. In this study, a novel arrangement of bionic drag-reducing microstructures is employed to modify the turbulence structure surrounding the submarine by incorporating bionic cone microstructures at both the front and rear ends of the submarine. Numerical simulations were performed using the SST k-ω turbulence model to evaluate the impact of these frontal microstructures on drag reduction under varying Reynolds numbers, spacings, and positions, as well as the tail microstructures’ effect at different Reynolds numbers, heights, and circumferential separation angles. The findings reveal that positioning microstructures at the submarine’s head increases the drag reduction rate proportionally with the distance from the apex, displaying an inverse relationship between spacing and drag reduction rate. Conversely, an increase in cone separation angle at the tail leads to a decrease in the overall drag reduction rate. At the same time, an inverse proportionality is observed between cone height and drag reduction rate. This suggests that cone microstructures play a dual role: mitigating friction drag greatly and augmenting pressure drag, thereby achieving overall drag reduction. Moreover, these cone microstructures disrupt eddy currents within the boundary layer surrounding the submarine, restraining the propagation of turbulent momentum transfer in both the head and tail regions. This research not only pioneers a novel drag reduction strategy for underwater vehicles but also sparks new avenues for their optimized surface design.

Development of a modified thermal phase change model based on the interfacial area concentration determined by an interface reconstruction algorithm

In order to enhance the physical basis of the mass transfer model in numerical simulation of gas-liquid direct contact condensation and improve the fidelity, a set of methods and models are developed. Based on the concept of the piecewise-linear interface calculation (PLIC), an interface reconstruction algorithm is proposed to calculate the interfacial area concentration in three-dimensional non-orthogonal structured grids. The actual interface in the cell is replaced by a plane, whose position is determined iteratively to satisfy the local volume fraction and orientation of the actual interface. The area of the polygonal cut by the plane is taken as the interfacial area in the cell. The proposed algorithm is validated by calculating the surface area of a series of basic geometries. Obvious superiority of the proposed algorithm over the traditional models (gradient model and algebraic model) is seen. Then, the reconstructed area is used to determine the interfacial area concentration in the thermal phase change model. The process is implemented by ANSYS Fluent user-defined functions. The Stefan problem and steam bubble condensation experiment are used to test the performance of the developed methods and models. The simulation results agree well with the theoretical and experimental results with the quantitative relative errors no more than 4.3 % and 15 %, respectively. The oxygen jets condensation experiments are also used as test cases. The primary dominant frequencies predicted by the simulations are 163.01 Hz and 131.65 Hz, which agree well with the 139.72 Hz and 125.75 Hz in the experiments. The relative errors are only 17 % and 4.7 %. The unstable flow pattern of oxygen jet in liquid oxygen crossflow is explained based on the dynamic characteristics of turbulence, heat and mass transfer. The periodic variation of mass transfer rate is considered as the cause of condensation oscillation. The method for calculating the interfacial area proposed in this paper can serve as an optional and practical improvement to the numerical simulation of direct contact condensation.

Turbulent momentum and kinetic energy transfer of channel flow over three-dimensional wavy walls

Turbulent momentum and kinetic energy transfer of channel flow over three-dimensional wavy walls

Spectral energy transfer analysis of a forced homogeneous isotropic turbulence using triple decomposition of velocity gradient tensor

Triple decomposition of the velocity gradient tensor is exploited to investigate the role of the simple-shear, normal-strain, and rigid-body rotation, in the turbulent kinetic energy transfer of a forced, high-Reynolds-number homogeneous isotropic turbulent flow. Splitting the total energy flux into three elementary partial fluxes shows that the rigid-body rotation has a secondary contribution to turbulent kinetic energy transfer, while within 80% of the total energy flux is produced through the simple-shear and normal-strain processes. Energy is dominantly extracted through the shearing process from larger motion scales and injected into smaller scales by the normal-straining process. In addition, each partial energy flux contains a strong conservative part and a weak non-conservative part. Besides the total kinetic energy transfer, all partial energy fluxes also show scale-invariance behaviour over the inertial subrange. Statistically, total energy flux shows a strong correlation with the simple-shear and normal-strain processes. However, the relatively poor correlation between these two dominant processes, implies a fairly low efficiency for turbulent kinetic energy transfer. These results can provide new insights into sub-grid scale turbulence modelling based on the physical mechanism of simple-shear and normal-strain processes.

Physics-informed neural network integrate with unclosed mechanism model for turbulent mass transfer

Turbulent mass transfer is widespread in chemical engineering processes including separation, reaction, and others. Traditionally, turbulent mass transfer processes can be simulated by computational fluid dynamics (CFD) methods. However, CFD methods require adequate turbulent models for the fluid flow and mass (and heat) transfer, which are normally difficult to develop for reliable solutions. Moreover, the CFD simulation is computationally intensive and hard to use in the cases like process optimization or analysis where the simulation needs to be called repeatedly. In this paper, physics-informed neural network (PINN) is developed and trained by unclosed mechanism model and sparse observation data of turbulent mass transfer process. The PINN method shows stronger generalization ability in solving the velocity, pressure and concentration fields in a turbulent mass transfer process than traditional deep neural network (DNN) method and turbulent Schmidt model. Under different boundary conditions, the PINN developed in the present paper can instantly predict the concentration distributions with sufficient accuracy and be used for inverse computation for estimating the turbulent viscosity and mass diffusivity as output results. The PINN is also capable of handling data noise by adjusting parameters, suggesting its potential in integrating experimental data.

Modeling snowpack dynamics and surface energy budget in boreal and subarctic peatlands and forests

Abstract. The snowpack has a major influence on the land surface energy budget. Accurate simulation of the snowpack energy and radiation budget is challenging due to, e.g., effects of vegetation and topography, as well as limitations in the theoretical understanding of turbulent transfer in the stable boundary layer. Studies that evaluate snow, hydrology and land surface models against detailed observations of all surface energy balance components at high latitudes are scarce. In this study, we compared different configurations of the SURFEX land surface model against surface energy flux, snow depth and soil temperature observations from four eddy-covariance stations in Finland. The sites cover two different climate and snow conditions, representing the southern and northern subarctic zones, as well as the contrasting forest and peatland ecosystems typical for the boreal landscape. We tested different turbulent flux parameterizations implemented in the Crocus snowpack model. In addition, we examined common alternative approaches to conceptualize soil and vegetation, and we assessed their performance in simulating surface energy fluxes, snow conditions and soil thermal regime. Our results show that a stability correction function that increases the turbulent exchange under stable atmospheric conditions is imperative to simulate sensible heat fluxes over the peatland snowpacks and that realistic peat soil texture (soil organic content) parameterization greatly improves the soil temperature simulations. For accurate simulations of surface energy fluxes, snow and soil conditions in forests, an explicit vegetation representation is necessary. Moreover, we demonstrate the high sensitivity of surface fluxes to a poorly documented parameter involved in snow cover fraction computation. Although we focused on models within the SURFEX platform, the results have broader implications for choosing suitable turbulent flux parameterization and model structures depending on the potential use cases for high-latitude land surface modeling.