To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
The aim of this paper is to characterise the internal structures and ice-flow history of representative valley glaciers in Svalbard and infer from them dynamic changes over centennial timescales. Three polythermal and one cold valley glacier are investigated using field- and laboratory-based techniques and remote sensing. Structures along flow-unit boundaries indicate that ice-flow configuration in three of the glaciers has remained stable spanning the residence time of the ice. Deformation of a flow-unit boundary in the fourth reveals an ice-flow instability, albeit one that has been maintained since its most recent advance. Macro-crystallographic, sedimentological and isotopic analyses indicate that basal ice is elevated to the glacier surface, as shown by entrained sediments and enrichment in heavy isotopes. In narrow zones of enhanced cumulative strain, new ice facies are generated through dynamic recrystallisation. The surface density of longitudinal foliation is shown to represent the relative magnitude of cumulative strain. Geometric similarities between flow-unit boundaries in Svalbard valley glaciers and larger scale longitudinal surface structures in ice sheets suggest that deformation mechanisms are common to both.
We explore the interaction of natural convection and mechanical ventilation in a room where fresh air is supplied at low level and stale air is extracted at high level. Turbulent buoyant plumes rising from heat sources interact with this upward airflow and establish a steady-state stratification with a warm upper layer above a layer of the cold supply air. Adapting the volume balance model used in natural ventilation (Linden et al., J. Fluid Mech., vol. 212, 1990, pp. 309–335) leads to the prediction that the upper layer will vent from the room when the ventilation volume flux exceeds the volume flux in the plumes at the ceiling. However, our new laboratory experiments establish that there is still a stable two-layer stratification beyond this point of critical ventilation. Motivated by our observations, we propose that the kinetic energy flux supplied by the plume leads to turbulent mixing in the upper layer. We propose a new model of this mixing which is consistent with our experiments in both the over- and under-ventilated regimes. This has important implications for air recirculation in buildings with large ventilation flows, particularly hospital operating theatres and clean rooms.
Liquid flowing down a fibre readily destabilises into a train of beads, commonly called a bead-on-fibre pattern. Bead formation results from capillary-driven instability and gives rise to patterns with constant velocity and time-invariant bead frequency $f$ whenever the instability is absolute. In this study, we develop a scaling law for $f$ that relates the Strouhal number $St$ and capillary number $Ca$ for Ostwaldian power-law liquids with Newtonian liquids recovered as a limiting case. We validate our proposed scaling law by comparing it with prior experimental data and new experimental data using xanthan gum solutions to produce a low capillary number $Ca \leq 10^{-2}$ regime. The experimental data encompasses both Ostwaldian and Newtonian flow, as well as symmetric and asymmetric patterns, and we find the data collapses along the predicted trend across seven orders of magnitude in $Ca$. Our proposed scaling law is a powerful tool for studying and applying bead-on-fibre flows where $f$ is critical, such as heat and mass transfer systems.
A calcium-silicate xenolith (no. 11) from the ignimbrite of the Upper Chegem Caldera in Kabardino-Balkaria, Russia, has revealed a diverse mineral assemblage with As- and B-bearing phases from the apatite supergroup such as the svabite and johnbaumite–hydroxylellestadite series, in addition to cahnite and datolite. Three distinct zones of variable arsenic content have been investigated. Notably, the outermost altered zone adjacent to the ignimbrite hosts the highest concentration of arsenic and arsenate minerals. A detailed structural analysis using Raman spectroscopy was carried out to investigate the distribution of boron and arsenic in tetrahedral coordination. This has provided the basis for describing a solid-solution system between hydroxylellestadite, svabite and johnbaumite and can be used as a novel technique for identifying apatite-supergroup minerals. One aim of the analysis was to elucidate the origin of various elements and content levels, particularly in relation to the distance from the xenolith–ignimbrite contact. The presence of boron and arsenic, probably derived from ignimbrites, highlights the important role of volcanic rocks as potential contributors of these elements in mineral formation processes.
Three metamict thorium orthosilicate samples from the syenite pegmatites of the Larvic Plutonic Complex, Norway, were thoroughly examined using Raman spectroscopy, electron probe microanalyses (EPMA), electron back-scatter diffraction (EBSD) and differential scanning calorimetry (DSC). Their thermal evolution upon heating was investigated using in situ powder X-ray diffraction (HTXRD) in the range of 25–1200°C. One of the samples is a colour-zoned metamict thorium silicate with a preserved tetragonal shape. The zonation is due to the increasing hydration and element distribution. The EBSD indicates that the ratio of huttonite to thorite after the crystallisation significantly varies from zone to zone within the same sample. The crystallisation of thorite starts in the range of 420–480°C (lower than reported previously for mineral samples), while the emergence of huttonite peaks in HTXRD patterns occurs at 870–930°C. In contrast to huttonite, no thorite crystallisation peak is observed in the DSC curve. A wide temperature range is observed where both thorite and huttonite can coexist. Several fluorite-type phases form upon heating. Thorianite exists in the range of 810–1140°C. After the cooling, except for huttonite and thorite, the minor crystallised phases vary and may be represented by Ca–Th oxides and rhombohedral CaUO4.
The thermal expansion of the crystalline huttonite and thorite was determined as $\overline\alpha$V = 20.66 × 10–6 deg–1 for huttonite and $\overline\alpha$V = 12.54 × 10–6 deg–1 for thorite in the temperature range 25–1200°C. These findings contribute to a more in-depth understanding of the behaviour of thorium orthosilicates with complex compositions, both metamict and crystalline, at elevated temperatures. They have potential applications in mineralogy, nuclear chemistry and high-level waste management.
This work presents the results of an investigation of an assemblage of secondary Sc-minerals from the intraplutonic metaluminous pegmatite Kožichovice II, Třebíč Pluton, Czech Republic. The assemblage was formed by hydrothermally-induced dissolution of primary Sc-enriched (≈1.6 wt.% Sc2O3) columbite-(Mn) followed by in situ reprecipitation of volumetrically dominant fersmite (≈0.16 wt.% Sc2O3) and minor nioboheftetjernite (ScNbO4). Subsequent hydrothermal processes resulted in the formation of fluorcalciomicrolite + Sc-minerals (thortveitite + kristiansenite) + titanite. The mass balance calculations (based on EPMA-derived mineral compositions, mineral proportions obtained from TIMA automated mineralogy and textural observations) revealed that the amount of Sc released from the replaced mass of columbite-(Mn) is significantly higher than the amount of Sc incorporated in the mass of the secondary minerals. This indicates that part of the Sc was mobilised and released to the host rocks (pegmatite and granite). The secondary mineral assemblages indicate elevated Ca activity in the alteration fluids. Other occurrences of Sc-minerals in pegmatites (Baveno Pluton and Heftetjern pegmatite) show remarkable similarities in the paragenetic position of Sc-minerals (late hydrothermal/replacement minerals), including the high activity of Ca in fluids during their formation. The high activity of Ca in fluids during the metasomatic replacement of Sc-enriched precursors causes the formation of the volumetrically dominant Sc incompatible phases, followed by a local supersaturation of Sc resulting in the crystallisation of secondary Sc-minerals.
As a natural clay mineral, halloysite (Hal) possesses a distinctive nanotubular morphology and surface reactivity. Hal calcined at 750°C (Hal750°C; 0.0, 1.0, 2.0, 4.0, 6.0, 8.0 wt.%) was used to replace ground granulated blast furnace slag (GGBFS; 50.0, 49.5, 49.0, 48.0, 47.0, 46.0 wt.%) and fly ash (FA; 50.0, 49.5, 49.0, 48.0, 47.0, 46.0 wt.%) for the preparation of geopolymer in this study. The effects of the replacement ratio of Hal750°C on setting time, compressive strength, flexural strength, chemical composition and microstructure of the geopolymer were investigated. The results indicated that Hal750°C did not significantly alter the setting time. The active SiO2 and Al2O3 generated from Hal750°C participated in the geopolymerization, forming additional geopolymer gel phases (calcium (aluminate) silica hydrate and sodium aluminosilicate hydrate), improving the 28 day compressive strength of the geopolymers. When the amount of Hal750°C was 2.0 wt.%, the 28 day compressive strength of the ternary (GGBFS-FA-Hal750°C-based) geopolymer was 72.9 MPa, 34.8% higher than that of the geopolymer without the addition of Hal750°C. The special nanotubular morphology of residual Hal750°C mainly acted like reinforcing fibres, supplementing the flexural strength of the geopolymer. However, excessive Hal750°C addition (>4.0 wt.%) reduced compressive and flexural strength values due to the low degrees of geopolymerization and the porous microstructure in the ternary geopolymer. These findings demonstrate that the appropriate addition of Hal750°C improved the compressive strength of geopolymers prepared using GGBFS/FA, which provides essential data for future research and supports the utilization of low-value Hal-containing clays in geopolymer preparation.
The eddy-viscosity model, as initially used to model the mean Reynolds stress, has been widely used in the linear analysis of turbulence recently by direct extension. In this study, the mechanism of the eddy viscosity in improving the prediction of fluctuation structures with linear analysis is clarified by investigating the statistical properties of forcing, eddy-viscosity term and their correlations. From the direct numerical simulation (DNS) results of turbulent channel flows with $Re_{\tau }=186$–$2003$, the spatial correlation of forcing is partially cancelled due to its interaction with eddy-viscosity terms. The stochastic forcing after excluding the eddy-viscosity term is nearly uncorrelated spatially, which better matches the condition where the resolvent modes are consistent with the spectral proper orthogonal decomposition (SPOD) modes from DNS. With the above findings, an optimisation framework is developed by minimising the spatial correlations of the stochastic forcing. The optimised eddy-viscosity profiles nearly overlap with the mean-quantity-based model in the near-wall region, but have different maximum values. Compared with the mean-quantity-based model, the optimised results enhance the consistency between the resolvent and DNS results significantly. Based on the optimised results, a new modelling framework is further abstracted, leaving only one to-be-modelled parameter of the maximum value of the eddy-viscosity profile. This parameter follows distinctive rules with spanwise flow scales, based on which a simplified predictive model is constructed. The resolvent modes predicted by the new model exhibit high consistency with SPOD modes, which are essentially comparable to the optimised results for wide ranges of streamwise and spanwise scales.
Parisite-(Nd) (IMA2024-013), ideally CaNd2(CO3)3F2, as the Nd-dominant analogue of parisite-(Ce), occurs in dolomitic marble in the Bayan Obo Fe–Nb–REE deposit, Inner Mongolia, China. It is associated with calcite, aegirine, magnetite, hematite, fluorite, riebeckite, bastnäsite-(Ce), baryte, aeschynite-(Ce), aeschynite-(Nd), monazite and parisite-(Ce). Parisite-(Nd) occurs as subhedral to anhedral irregular grains from 0.02 mm to 1 mm. Parisite-(Nd) is transparent, yellowish-brown colour, with pale yellow streak and displays vitreous to resinous lustre. Cleavage is distinct on pseudo-{001}; fracture is uneven, or conchoidal. The Mohs hardness is 4 to 5, and it is brittle. The calculated density of parisite-(Nd) is 4.357 g/cm3. Parisite-(Nd) is pseudo-uniaxial (+), ω = (1.679) and ε = (1.754). The empirical formula is (Ca0.945Fe0.058Sr0.015Ba0.007)Σ1.025(Nd0.967Ce0.529La0.191Pr0.137Gd0.070Sm0.029Th0.022Y0.016Nb0.011Ho0.003)Σ1.975(CO3)3F1.893OH0.023. The Raman spectra of parisite-(Nd) show strong and sharp peaks at 1113, 1090, 825, 635 and 1608 cm–1 and moderate to weak bands centred at 255, 392, 739, 924, 1183, 1228, 1296, 1640, 2247, 2924 and 3065 cm–1. Powder X-ray diffraction and TEM studies give the following results: monoclinic, space group: Cc (# 9), a = 12.3283(13) Å, b = 7.1185(4) Å, c = 28.4633(37) Å, β = 98.529(14)°, V = 2470.28(42) Å3 and Z = 12.
The relevance of education and outreach (E&O) activities about the Antarctic Treaty has been recognized at the Antarctic Treaty Consultative Meetings (ATCM) and at the Committee for Environmental Protection (CEP). This study examines the key topics and the target audiences detailed in papers submitted to the ATCM on E&O. Since the Antarctic Treaty entered into force in 1961, a total of 216 ATCM papers on E&O have been produced. The number of papers has increased substantially since the mid-1990s. ‘Science’ (76.9%) and ‘Wildlife/Biodiversity/Environment’ (75.5%) were the most addressed topics in these papers, while the ‘Public’ (81.0%) and those attending ‘Schools’ (69.0%) are the main target audiences. ‘Science’ in ATCM papers increased ~120-fold from 1961–1997 to 2015–2023, while ATCM papers discussing engagement with the ‘Public’ increased ~40-fold during the same period. ‘Climate change’ was first mentioned in 2006, and the number of papers per year increased fourfold by 2015–2023. This study shows the increasing interest in E&O through time, addressing key topics to relevant audiences related to the Antarctic region. From an educational perspective, attention should be paid to emerging topics (e.g. equity, diversity and inclusion), and the engagement of early-career professionals and educators should be made a priority.
The presence of dispersed-phase droplets can result in a notable increase in a system's drag. However, our understanding of the mechanism underlying this phenomenon remains limited. In this study, we use three-dimensional direct numerical simulations with a modified multi-marker volume-of-fluid method to investigate liquid–liquid two-phase turbulence in a Taylor–Couette geometry. The dispersed phase has the same density and viscosity as the continuous phase. The Reynolds number $Re\equiv r_i\omega _i d/\nu$ is fixed at 5200, the volume fraction of the dispersed phase is up to $40\,\%$, and the Weber number $We\equiv \rho u^2_\tau d/\sigma$ is approximately 8. It is found that the increase in the system's drag originates from the contribution of interfacial tension. Specifically, droplets experience significant deformation and stretching in the streamwise direction due to shear near the inner cylinder. Consequently, the rear end of the droplets lags behind the fore head. This causes opposing interfacial tension effects on the fore head and rear end of the droplets. For the fore head of the droplets, the effect of interfacial tension appears to act against the flow direction. For the rear end, the effect appears to act in the flow direction. The increase in the system's drag is attributed primarily to the effect of interfacial tension on the fore head of the droplets which leads to the hindering effect of the droplets on the surrounding continuous phase. This hindering effect disrupts the formation of high-speed streaks, favouring the formation of low-speed ones, which are generally associated with higher viscous stress and drag of the system. This study provides new insights into the mechanism of drag enhancement reported in our previous experiments.
In this paper, aerobreakup in the stagnation region of high-Mach-number flow over a bluff body is studied with experiments and computations. Water drops of diameter 0.51–2.30 mm were acoustically levitated at sea level along the flight path of a rectangular $100\ {\rm mm} \times 150\ {\rm mm}$ rail-gun launched projectile. This enabled the study of aerobreakup at high Mach (3.03–5.12), post-shock Mach (1.5–1.9), Weber $(5 \times 10^4\unicode{x2013}4 \times 10^5)$ and Reynolds $(6 \times 10^4\unicode{x2013}3 \times 10^5)$ numbers. High-speed backlit shadowgraphy is used to record the flow structure. Computations are made for two cases, and it was found that the drop behaviour is not dominated by viscous or surface-tension effects and can be adequately captured by treating the gas as calorically perfect with the ratio of specific heats set to 1.3 to account for thermochemical effects. To assess drop surface stability at early breakup times, results from Newton's inclination method are used to determine the flow along the drop surface and input to a linear-stability analysis; from this, it was found that viscosity and surface tension can be neglected. Moreover, the acceleration term dominates the shear term at the stagnation point, a point accentuated as a drop flattens; this relation inverts closer to the drop equator. Linear-stability analysis was insufficient for modelling late-stage aerobreakup because the predicted wavelengths were too small and the expected aerobreakup times were non-physically short. To address this discrepancy, a nonlinear instability model with constant-rate growth is used that treats the accelerated drop surface as analogous to bubbles rising through a liquid; agreement with computations is good.
Modern fluvial sediments provide important information about source-to-sink process and regional tectono-magmatic events in the source area, but many factors, e.g., chemical weathering, sedimentary cycles and source-rock types, can interfere with the establishment of the source-sink system. The Lalin River (LR) and the Jilin Songhua River (JSR) are two important tributaries of the Songhua River in the Songnen Plain in NE China. They have similar flow direction, topography and identical climate backgrounds, but have notably different parent-rock types in the headwater, which provides an opportunity to explore the influencing factors of river sediment composition. To this end, the point bar sediments in the two rivers were sampled for an analysis of geochemistry (including element and Sr-Nd isotopic ratios), heavy mineral and detrital zircon U-Pb dating. The results are indicative of the fact that the two rivers have the similar geochemical composition (e.g., elements and Sr isotopes) as well as chemical weathering (CIA = 51.41–57.60, CIW = 59.68–66.11, PIA = 51.95–60.23, WIP = 56.00–65.47, Rb/Sr = 0.38–0.42) and recycling (SiO2/Al2O3 = 5.79 and 5.03, ICV = 1.0 and 1.2, CIA/WIP = 0.81–1.03) characteristics, showing a major control of climate on the low-level weathering and recycling of the river sediments. However, there are significant differences in the detrital zircon U-Pb age (a significant Mesozoic age peak for the LR but an additional Precambrian peak for the JSR), Nd isotope ratio (−6.2812–8.5830 and −8.1149–10.2411 for the LR and the JSR, respectively) and to a certain extent heavy mineral composition (e.g., for the < 63 μm fraction, a dominance of hornblende and magnetite in the LR, but haematite-limonite in the JSR) in the two river sediments, indicating that source rocks largely control the composition of the river sediments. Some of the major tectono-magmatic events (e.g., crustal growth and cratonisation of the North China Craton, closure of the Paleo-Asian Ocean, subduction and rollback of the Paleo-Pacific plate) occurring in the eastern Songnen Plain are well documented in the JSR sediments but not in the LR, the difference of which is largely regulated by the source rocks in the source area.
Different two-dimensional structural units of layered silicate minerals have different chemical and reaction properties. Sulfuric acid solution mineral-leaching systems with pH of 2.0, 4.0 and 6.0 were constructed to investigate the differential dissolution properties of lizardite (1:1 type), chlorite and talc (2:1 type minerals) and the chemical kinetic mechanism of the mineral–water interface reaction. The results showed that the dissolution efficiency of Mg in lizardite is higher than that of chlorite and talc in acidic environments (pH of 2.0, 4.0 and 6.0). The dissolution efficiency of Mg in chlorite is greater than that of talc for acidic environments when pH is 2.0 and 4.0, but chlorite and talc have nearly identical Mg dissolution efficiencies at a pH of 6.0. This phenomenon is related to the defect site on the tetrahedral sheet of chlorite and is controlled by the change of the dissolution efficiency of Al. The dissolution rates of Mg and Si in lizardite, chlorite and talc decreased with the increase of reaction time in the acidic medium for pH = 2.0, 4.0 and 6.0, and there are two linear dissolution trends at different pH values. The dissolution efficiencies of Mg and Si in lizardite, chlorite and talc were simulated and predicted by establishing a logistic model. It was found that the maximum dissolution efficiency of 2:1 type minerals chlorite and talc are only 4.72% and 1.58%, which is much lower than that of 1:1 type lizardite. This research on the reaction mechanism and dissolution kinetics of lizardite, chlorite and talc not only helps to deepen the understanding of the mineral–water interface interaction, but also reveals the different rules for Mg, Si and Al dissolution in different types of trioctahedral mineral–water interface reactions, and provides a crystal chemical basis for the ion migration and action mechanism of minerals.
The study area Sonapahar is an integral part of Shillong Meghalaya Gneissic Complex (SMGC), which is located in the Northeastern part of India. This complex mainly comprises metamorphic formations spanning from Upper Amphibolite to Ultra-high temperature granulite, interspersed with various igneous intrusions. In this study, particular attention is directed towards unravelling the metamorphic history of Mg-Al granulite. For the very first time, we establish the pressure–temperature (P-T) trajectory of the Mg-Al granulite from Sonapahar, SMGC. Employing conventional thermobarometry along with winTWQ analysis, the inferred metamorphic conditions for this granulite reveal temperatures exceeding 900°C and pressures of approximately >8 kbar. These conditions firmly indicate the presence of ultra-high-temperature metamorphism. By utilizing the Perple_X software in the MnO-Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O-TiO2-Fe2O3 compositional system, we construct a P-T pseudosection. This gives a clockwise P-T path, signifying an episode of cooling (+ minor decompression). Such a pattern also suggests rapid cooling of the tectonically-thickened crust. Concurrently, a geochemical exploration of trace and rare earth elements in the rocks offers further insights. These investigations give an idea about the protolith, having a clay-to-sandstone in nature. Additionally, chemical data from monazite within the studied rock provide a weighted mean age of 682 Ma for the peak metamorphic stage. This age aligns with the global Pan-African orogenic events. The biotite K-Ar isotopic geochronology from the symplectite position provides decompression history or cooling age of 442 Ma. This age corresponds to a period after the last peak metamorphic phase that occurred during the Pan-African thermal event.
Turbulent mixing is a pivotal phenomenon in fusion research with profound implications for energy gain. A Reynolds-averaged Navier–Stokes model capable of predicting realistic mixing transition processes is of significant importance for fusion applications, yet such a model remains elusive. This work addresses the limitations of prevalent global transition criteria, proposing a new idea to quantify local transition characteristics based on the mixing state, recognizing its direct relevance to fusion reaction rates. We delve into the description and analysis of the spatiotemporal evolution of the mixing state and its interplay with the transition process. Then, a local transition indicator is developed and compared with conventional global criteria using the large-eddy simulation (LES) of Rayleigh–Taylor turbulent mixing. Building upon this foundation, we introduce a novel eddy viscosity model based on the local transition indicator. A posterior assessment using LES data validates that it significantly outperforms standard gradient transport models during the transition stage. Consequently, we integrate this new eddy viscosity model with the Besnard–Harlow–Rauenzahn model to construct a comprehensive transition model, which demonstrates reasonably good performance in comparison with LES results. This work paves the way for future research in developing advanced modelling strategies that can effectively address the complexities of transitional flows in fusion engineering applications.
Peatlands, covering approximately one-third of global wetlands, provide various ecological functions but are highly vulnerable to climate change, with their changes in space and time requiring monitoring. The sub-Antarctic Prince Edward Islands (PEIs) are a key conservation area for South Africa, as well as for the preservation of terrestrial ecosystems in the region. Peatlands (mires) found here are threatened by climate change, yet their distribution factors are poorly understood. This study attempted to predict mire distribution on the PEIs using species distribution models (SDMs) employing multiple regression-based and machine-learning models. The random forest model performed best. Key influencing factors were the Normalized Difference Water Index and slope, with low annual mean temperature, with low annual mean temperature, precipitation seasonality and distance from the coast being less influential. Despite moderate predictive ability, the model could only identify general areas of mires, not specific ones. Therefore, this study showed limited support for the use of SDMs in predicting mire distributions on the sub-Antarctic PEIs. It is recommended to refine the criteria used to select environmental factors and enhance the geospatial resolution of the data to improve the predictive accuracy of the models.
Developing large-eddy simulation (LES) wall models for separated flows is challenging. We propose to leverage the significance of separated flow data, for which existing theories are not applicable, and the existing knowledge of wall-bounded flows (such as the law of the wall) along with embedded learning to address this issue. The proposed so-called features-embedded-learning (FEL) wall model comprises two submodels: one for predicting the wall shear stress and another for calculating the eddy viscosity at the first off-wall grid nodes. We train the former using the wall-resolved LES (WRLES) data of the periodic hill flow and the law of the wall. For the latter, we propose a modified mixing length model, with the model coefficient trained using the ensemble Kalman method. The proposed FEL model is assessed using the separated flows with different flow configurations, grid resolutions and Reynolds numbers. Overall good a posteriori performance is observed for predicting the statistics of the recirculation bubble, wall stresses and turbulence characteristics. The statistics of the modelled subgrid-scale (SGS) stresses at the first off-wall grids are compared with those calculated using the WRLES data. The comparison shows that the amplitude and distribution of the SGS stresses and energy transfer obtained using the proposed model agree better with the reference data when compared with the conventional SGS model.
Convective flows near a cylindrical solid inclusion in a fluid-saturated porous medium heated from above are studied using the Darcy–Boussinesq equations. The impermeability condition is imposed on the inclusion surface. Two-dimensional convective flows uniform in the direction of the cylinder axis are considered. An Oseen-like approximation is implemented, however, different from the conventional Oseen approximation, quasilinearization is applied to the nonlinear terms in the energy equation and not in the momentum equation. It is shown that, when the thermal conductivity of the inclusion is higher than that of the fluid, the convective flow at a distance from the inclusion larger than the inclusion size takes the form of horizontal vortices directed away from the inclusion. In the case of low thermal conductivity of the inclusion, the direction of convective circulation is opposite: in the horizontal plane passing through the cylinder axis, the fluid is leaking to the inclusion.