Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary layer analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Exploring Causal Relationships and Adjustment Timescales of Aerosol-Cloud Interactions in Geostationary Satellite Observations and CAM6 Using Wavelet Phase Coherence Analysis

We present for the first time within the cloud physics context, the application of wavelet phase coherence analysis to disentangle counteracting physical processes associated with the lead-lag phase difference between cloud-proxy liquid water path (LWP) and aerosol-proxy cloud droplet number concentration ( N d ) in an Eulerian framework using satellite-based observations and climate model outputs. This approach allows us to identify the causality and dominant adjustment timescales governing the correlation between LWP and N d . Satellite observations indicate a more prevalent positive correlation between daytime LWP and N d regardless of whether LWP leads or lags N d . The positive cloud water response, associated with precipitation processes, typically occurs within 1 hr, while the negative response resulting from entrainment drying, usually takes 2–4 hr. CAM6 displays excessively rapid negative responses along with overly strong negative cloud water response and insufficient positive response, leading to a more negative correlation between LWP and N d compared to observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Observations of Offshore Low‐Level Jets Off the U.S. East Coast Reveal Systematic Biases in ERA5 and HRRR

Low-level jets (LLJs)—wind speed maxima typically occurring a few hundred meters above the surface—are common off the U.S. East Coast and influence many atmospheric processes with societal importance, including cloud formation, aviation safety, and search-and-rescue. However, their vertical structure and frequency remain poorly quantified due to limited offshore observations. This study presents new scanning Doppler LiDAR and infrared spectroradiometer data from the 2024 summer deployment of an offshore barge during the Wind Forecast Improvement Project 3. These coupled wind and temperature profiles provide unprecedented resolution to assess LLJ behavior and model performance. LLJs occurred in over 21% of observed profiles, with a weak diurnal preference for nighttime and early morning hours and maximum winds typically near 300 m. Both ERA5 and High-Resolution Rapid Refresh analysis underestimate jet wind speeds and misrepresent the boundary layer thermal structure. These results highlight persistent model biases and the critical need for high-resolution offshore observations.

17 WIND ENERGY↗
DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗
DOE OSTI2025

Impact of inland boundary conditions on seawater intrusion modeling in stratified coastal aquifers under sea-level rise

This study systematically examines the impact of inland boundary conditions on seawater intrusion modeling in general stratified coastal aquifers under sea-level rise (SLR). We derived analytical solutions for seawater wedge toe positions in stratified aquifers with general head boundary (GHB) conditions and compiled solutions for constant head boundary (CHB), constant flux boundary (CFB), and GHB conditions in both homogeneous and stratified aquifers. Additionally, we introduced the concept of equivalent boundary conditions and developed equations to ensure consistent steady-state toe positions across different boundary types. Unlike previous studies focusing on specific layering configurations, we conducted a comprehensive sensitivity analysis of seawater intrusion in general stratified aquifers using transmissivity centroid elevation (TCE). Our results show that SLR generally exacerbates seawater intrusion, with unconfined aquifers exhibiting higher sensitivity across all boundary conditions. Among boundary types, CHB conditions were the most sensitive to SLR, CFB the least, and GHB showed an intermediate response depending on its distance from the model boundary. Furthermore, seawater intrusion increases with the TCE of the stratified aquifer under SLR, suggesting that high-permeability or preferential flow layers near the surface amplify the impact of SLR. These findings emphasize the importance of selecting appropriate boundary conditions when evaluating SLR-driven seawater intrusion in stratified coastal aquifers.

Boundary condition↗
DOE OSTI2024

Stress localization investigation of additively manufactured GRCop-42 thin-wall structure

A full-field crystal plasticity (CP) framework is presented for the GRCop-42 alloy to study microscopic mechanical behavior and local stress heterogeneities. The microstructures of additively manufactured (AM) materials are often unique relative to conventionally processed materials, and the local thermal histories drive these differences during the build process. These thermal histories depend on the process parameters (laser power, scan speed, and scan strategy) and the part geometry. Prior research has shown that the mechanical properties of thin-walled structures can vary significantly with wall thickness due to changes in the thermal boundary conditions during manufacturing. It is, therefore, desirable to perform CP simulations based on the phenomenological constitutive model to predict the local mechanical responses induced by microstructural heterogeneities. This work generates representative microstructures based on experimentally collected grain information (i.e., texture) for grain scale stress analysis, and the material constitutive parameters are calibrated using the experimental mechanical testing data. Here, we specifically investigated the effect of crystallographic texture and grain morphologies on the size-dependent mechanical properties of AM GRCop-42. The selection of appropriate material properties for implementing an effective free surface boundary condition and the influence of adjacent buffer layers are also discussed. Analysis of local field results reveals a strong correlation between stress localization and the initial grain orientation. However, no significant relationship between the misorientation of the individual adjacent grains and the average misorientation is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectral distortions to momentum and scalar exchanges by non-turbulent motion and patchy landscape variability

Modifications to the spectra of turbulent velocity and scalars and co-spectra of vertical fluxes of momentum and scalars due to patchy landscape heterogeneity and non-stationarity are explored for a Mediterranean ecosystem. About 9 months of high frequency measurements of the three velocity components, water vapor concentration, carbon dioxide concentration, and air temperature were analyzed for different seasons (spring/summer) and prevalent wind directions (southeast/northwest). The two wind directions sampled a contrast of clumped and patchy landscape comprised of olive trees (southeast) and wall bounded flow disturbed by the presence of few upwind trees (northwest). The measured spectra and co-spectra were also compared to theoretical scaling forms from stationary, planar homogeneous flow, in the absence of subsidence as derived from the Kansas experiment. To assess the role of low frequency non-turbulent motion on the spectral and co-spectral content, a 5-min Fourier cutoff was introduced and the analysis was limited to near-neutral conditions where the boundary layer depth is shallow compared to its unstable counterpart. It was shown that the velocity statistics were not appreciably impacted by the low-frequency motion causing non-stationarity. Moreover, the turbulent scalar fluxes were also shown not to be significantly impacted by such low frequency motion. The scalar variances were impacted, especially the water vapor variance and its concomitant spectral shape. When the non-turbulent motion was filtered, the scalar spectra at low wavenumbers followed expectations from the so-called attached eddy hypothesis (i.e. exhibited a $k^{-1}_x$ scaling with $k_x$ defining the longitudinal wavenumber) applicable for near-neutral conditions. For momentum co-spectra, the canonical shapes from the Kansas experiment appear to describe well the measurements here and in both dominant directions and seasons with some adjustment to the integral time scales based on wind direction. For the scalar co-spectra, deviations from the Kansas experiment were prevalent. The most noticeable and surprising deviations were their slow decay with increased sampling frequency at inertial subrange scales. This slow decay was shown not to contribute appreciably to the overall scalar fluxes. At those fine scales, predictions from local isotropy were expected to hold. The scalar co-spectral deviations from local isotropy were then discussed using a simplified co-spectral budget model where scalar–scalar co-spectra naturally emerged and the interplay between landscape heterogeneity and a scale-dependent pressure-scalar de-correlation time was postulated. It is also envisaged that the findings here offer a preliminary template for analyzing eddy-covariance data in situations that deviate from ideal conditions, especially regarding low-frequency modulations of scalar spectra and vertical scalar flux co-spectra.

Canopy turbulence↗
DOE OSTI2025

Characterizing Seasonal Variation of the Atmospheric Mixing Layer Height Using Machine Learning Approaches

As machine learning becomes more integrated into atmospheric science, XGBoost has gained popularity for its ability to assess the relative contributions of influencing factors in the atmospheric boundary layer height. To examine how these factors vary across seasons, a seasonal analysis is necessary. However, dividing data by season reduces the sample size, which can affect result reliability and complicate factor comparisons. To address these challenges, this study replaces default parameters with grid search optimization and incorporates cross-validation to mitigate dataset limitations. Using XGBoost with four years of data from the atmospheric radiation measurement (ARM) (Southern Great Plains (SGP) C1 site, cross-validation stabilizes correlation coefficient fluctuations from 0.3 to within 0.1. With optimized parameters, the R value can reach 0.81. Analysis of the C1 site reveals that the relative importance of different factors changes across seasons. Lower tropospheric stability (LTS, ~0.53) is the dominant factor at C1 throughout the year. However, during DJF, latent heat flux (LHF, 0.44) surpasses LTS (0.22). In SON, LTS (0.58) becomes more influential than LHF (0.18). Further comparisons among the four long-term SGP sites (C1, E32, E37, and E39) show seasonal variations in relative importance. Notably, during JJA, the differences in the relative importance of the three factors across all sites are lower than in other seasons. This suggests that boundary layer development in the summer is not dominated by a single factor, reflecting a more intricate process likely influenced by seasonal conditions such as enhanced convective activity, higher temperatures, and humidity, which collectively contribute to a balanced distribution of parameter impacts. Furthermore, the relative importance of LTS gradually increases from morning to noon, indicating that LTS becomes more significant as the boundary layer approaches its maximum height. Consequently, the LTS in the early morning in autumn exhibits greater relative importance compared to other seasons. This reflects a faster development of the mixing layer height (MLH) in autumn, suggesting that it is easier to retrieve the MLH from the previous day during this period. The findings enhance understanding of boundary layer evolution and contribute to improved boundary layer parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characterizing the effect of hypersonic boundary layer turbulence on antenna performance: A computational approach

The degradation of antenna performance during hypersonic re-entry is a well known phenomenon that can lead to complete radio blackout. Recent additions to the Empire code establish it as a tool for the study and analysis of the problem. Coupling to the Sandia Parallel Aerodynamics and Reentry Code (SPARC) enables the electromagnetic analysis of realistic re-entry plasma profiles. The geometric flexibility afforded by both Empire and SPARC allow the consideration of arbitrary vehicle and antenna configurations. We have used this tool to study antenna performance during re-entry when the boundary layer becomes turbulent. A concise description of line-of-sight transmissions, which employs advanced statistical methods, was developed. New insights into the low altitude reflectometer readings of RAM-C2 are offered. Techniques for the reconstruction of the re-entry plasma profile from reflectometer data were explored.

42 ENGINEERING↗
DOE OSTI2026

Observed Land Surface Influence on Atmospheric Heat and Moisture Profiles During Interstorms

Land-atmospheric (L-A) feedbacks have historically been studied using models whose structure and parameterizations influence outcomes and insights. The representation of L-A feedbacks based on observations alone remains an ongoing challenge for understanding boundary layer development and precipitation. To address this gap, we use ground-based passive remote sensing and in-situ observations to present an analysis of the atmosphere during 103 interstorm soil moisture drydown events spanning nine warm seasons (2016–2024) in the U.S. Southern Great Plains region. By separating events based on local L-A coupling signals and characterizing the profiles of atmospheric heat and moisture to surface energy flux behavior, we investigate the physical mechanisms linking land surface processes to boundary layer development. We find that during interstorm drydowns, the atmospheric column follows a consistent pattern: moisture increases within the boundary layer, peaks near its top, and declines rapidly above, while warming occurs through the depth. Drydowns that shift toward evaporation produce stronger and deeper thermodynamic responses than cases dominated by sensible heating, which are weaker and shallower. Additionally, moisture is accumulated faster within the boundary layer during shorter drydowns, with longer drydowns representing slower, moisture-limited growth. Drydowns with wetter initial soil moisture will sustain stronger moistening within and above the boundary layer, accelerating buoyancy growth and convective potential toward the next storm. These results provide observational evidence linking surface flux evolution to boundary layer thermodynamics and offer a process-level benchmark for evaluating coupled L-A representations in models and demonstrating the influence of soil moisture on short-term weather forecasting skill.

Zhang, M. S. [Massachusetts Inst. of Technology (M↗
DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗
DOE OSTI2024

New insight into tuning magnetic phases of RMn6Sn6 kagome metals

Abstract Predicting magnetic ordering in kagome compounds offers the possibility of harnessing topological or flat-band physical properties through tuning of the magnetism. Here, we examine the magnetic interactions and phases of ErMn 6 Sn 6 which belongs to a family of R Mn 6 Sn 6 , R = Sc, Y, Gd–Lu, compounds with magnetic kagome Mn layers, triangular R layers, and signatures of topological properties. Using results from single-crystal neutron diffraction and mean-field analysis, we find that ErMn 6 Sn 6 sits close to the critical boundary separating the spiral-magnetic and ferrimagnetic ordered states typical for non-magnetic versus magnetic R layers, respectively. Finding interlayer magnetic interactions and easy-plane Mn magnetic anisotropy consistent with other members of the family, we predict the existence of a number of temperature and field dependent collinear, noncollinear, and noncoplanar magnetic phases. We show that thermal fluctuations of the Er magnetic moment, which act to weaken the Mn-Er interlayer magnetic interaction and quench the Er magnetic anisotropy, dictate magnetic phase stability. Our results provide a starting point and outline a multitude of possibilities for studying the behavior of Dirac fermions in R Mn 6 Sn 6 compounds with control of the Mn spin orientation and real-space spin chirality.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2026

Turbulence and Cloudiness in Cumulus-Topped Marine Boundary Layers

Turbulence and cloudiness in cumulus-topped marine boundary layers are studied using data from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site. The analysis includes eight periods of nonprecipitating shallow cumulus clouds, spanning 141 h and encompassing 603 individual clouds. On average, the cumulus had bases at 558 m, were 99 m thick, had a chord length of 500 m, and exhibited an hourly base-layer cloudiness of 12%. Changes in cloud fraction were primarily driven by variations in cloud number not chord length. High-resolution Doppler lidar and cloud radar observations were used to estimate updraft and downdraft mass fluxes in the cumulus base layer under both clear and cloudy regions. Cloudy updrafts contributed only ∼23% of the total updraft mass flux at cloud base, indicating that most upward mass flux originated from clear-air eddies. Updraft strength, rather than updraft fraction, was found to predominantly control both clear and cloudy updraft mass fluxes. Cloud-base cloudiness strongly correlated with cloudy updraft mass flux but showed negligible correlation with clear-air updraft mass flux. Clear-air downdraft mass flux exhibited a strong relationship with the ratio of surface buoyancy to inversion strength. Mesoscale analysis revealed that moist patches had greater cloudiness, updraft mass flux, and vertical velocity variance compared to dry patches. Additionally, mesoscale base-layer cloud fraction was highly and significantly correlated with cloudy updraft mass flux in both dry and moist environments. Results presented herein have implications for cumulus parameterization development along with climatological and model evaluation studies conducted at the ENA site.

cumulus clouds↗
DOE OSTI2025

Wind and Weather Variability within the Californian Offshore Wind Energy Areas

Weather variability over the Northeast Pacific (NEP) region and its influence on wind resources within the Californian offshore wind energy areas (WEAs) at Humboldt and Morro Bay are characterized using 20-years reanalysis model and satellite data. The hub-height (180 m) winds at both locations are predominantly northwesterly driven by the NEP high pressure system, with strong coastal gradients in surface pressure, fluxes, planetary boundary layer (PBL) depths and cloudiness. These sharp coastal gradients and strong annual cycles of temperature and moisture advections pose potential challenges in accurately modeling the local wind resource. Hub-height wind speeds and power capacity factors significantly vary for different regimes of PBL depths, surface fluxes and rain area fractions. This highlights the importance of studying the physical mechanisms driving these weather regimes, hence our analysis of how large-scale NEP weather variability drives the local meteorology at the WEAs. Furthermore, at both WEAs, PBL tops and cloud boundaries intersect the rotor layer (80-280 m) more than 30% and 20% of the time, respectively. While PBL depths significantly modulates hub-height winds and power, cloud boundaries do not have a similar impact, likely due to reanalysis errors in simulating cloud boundaries accurately. These findings underscore the challenges in deploying tall wind turbines in shallow cloudy boundary layers, where the interaction between clouds, precipitation, and atmospheric layers can impact turbine efficiency. As turbines grow taller and are deployed in more complex meteorological conditions, understanding these interactions is crucial for improving wind power forecasting and optimizing energy production in coastal regions.

17 WIND ENERGY↗
DOE OSTI2024

Diurnal evolution of non-precipitating marine stratocumuli in a large-eddy simulation ensemble

Abstract. We explore the cloud system evolution of non-precipitating marine stratocumuli with a focus on the impacts of the diurnal cycle and free-tropospheric (FT) humidity based on an ensemble of 244 large-eddy simulations generated by perturbing initial thermodynamic profiles and aerosol conditions. Cases are categorized based on their degree of decoupling and the cloud liquid water path (LWPc, based on model columns with cloud optical depths greater than one). A budget analysis method is proposed to analyze the evolution of cloud water in both coupled and decoupled boundary layers. More coupled clouds start with a relatively low LWPc and cloud fraction (fc) but experience the least decrease in LWPc and fc during the daytime. More decoupled clouds undergo greater daytime reduction in LWPc and fc, especially those with higher LWPc at sunrise because they suffer from faster weakening of net radiative cooling. During the nighttime, a positive correlation between FT humidity and the LWPc emerges, consistent with higher FT humidity reducing both radiative cooling and the humidity jump, both of which reduce entrainment and increase LWPc. The LWPc is more likely to decrease during the nighttime for a larger LWPc and greater inversion base height (zi), conditions under which entrainment dominates as turbulence develops. In the morning, the rate of the LWPc reduction depends on the LWPc at sunrise, zi, and the degree of decoupling, with distinct contributions from subsidence and radiation.

Chen, Yao-Sheng (ORCID:0000000208354132)↗
DOE OSTI2026

Aircraft-derived particle fluxes distinguish entrainment zone and decoupled layer nucleation in marine boundary layers

The vertical distribution of freshly nucleated aerosol particles in the marine boundary layer remains poorly constrained, limiting our ability to represent new particle formation in climate models. Here we characterize 3–10 nm particle events, termed small particle events (SPEs), by deriving their vertical turbulent fluxes from aircraft measurements during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. To overcome stationarity limitations of traditional eddy covariance methods, we applied continuous wavelet transform analysis to data collected during June–July 2017 and January–February 2018 flights over the Azores. Our flux-based analysis revealed two distinct SPE scenarios with different vertical structures and spatial extents. The first featured nucleation in the entrainment zone, where free tropospheric air entrains into the boundary layer. The second showed nucleation in the decoupled layer, a stratified region between the well-mixed surface layer and cloud-topped upper boundary layer. In both cases, convergence of air masses from different layers diluted preexisting aerosol surface area to very low levels, creating conditions favorable for nucleation and generating strong downward particle fluxes. SPEs occurred in 15 % of flights, challenging prevailing theoretical expectations that new particle formation should rarely occur in marine boundary layers due to high condensation and coagulation sink capacity of sea spray aerosols. Aircraft-derived particle fluxes provide first observational constraints on the vertical location and source strength of likely nucleation regions in the remote marine boundary layer, improving aerosol source representations in climate models and reducing uncertainties in aerosol-cloud interactions.

Rasheeda Satheesh, Ajmal [North Carolina State Uni↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗