[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (5,321)

Search Parameters:
Keywords = convection

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
12 pages, 543 KiB  
Article
Assessment of Safety and Efficacy of Expanded Hemodialysis with Medium Cut-Off Dialyzer Compared to Haemodiafiltration
by Matteo Marcello, Marco Simonini, Anna Lorenzin, Valentina Corradi, Grazia Maria Virzì, Carlotta Caprara, Alessandra Brendolan, Claudia Benedetti, Paolo Lentini, Monica Zanella and Claudio Ronco
J. Clin. Med. 2025, 14(6), 1798; https://doi.org/10.3390/jcm14061798 - 7 Mar 2025
Viewed by 50
Abstract
Background: Removal of large uraemic toxins is still a challenge. Haemodiafiltration (HDF) has produced some results, although large convective volume, optimal vascular access to increase the blood flow rate and strict water quality management are required. Medium cut-off, high-retention-onset membranes have been recently [...] Read more.
Background: Removal of large uraemic toxins is still a challenge. Haemodiafiltration (HDF) has produced some results, although large convective volume, optimal vascular access to increase the blood flow rate and strict water quality management are required. Medium cut-off, high-retention-onset membranes have been recently developed, introducing the concept therapy called expanded haemodialysis (HDx). Furthermore, vitamin E-coated membrane has potential beneficial effects on inflammation and oxidative stress. Methods: A prospective longitudinal multicentre study was conducted for 3 months among 24 chronic haemodialysis patients. Patients were randomly assigned into either HDF with high-flux membrane or HDx with Theranova or ViE-X membrane. The primary goal was to assess albumin loss among the three types of dialyzers. Secondary goals included assessment of depurative efficacy for uraemic toxins and clinical outcomes. Results: Mean albumin loss was significantly higher in patients undergoing HDx with Theranova membrane, without any difference in serum albumin concentration among the three groups. Instantaneous clearance of small and middle molecules was significantly higher in patients undergoing HDF, but we did not find differences in removal ratio and Kt/V. Reduction in the erythropoietin resistance index was observed in patients treated with ViE-X membrane due to their lower dialysis vintage. Conclusions: The higher albumin loss during HDx has no effects on pre-dialysis serum albumin. HDx with Theranova in the presence of lower session length, lower Qb, lower convective dose, and lower instantaneous clearance reached the same dialysis efficacy compared to HDF. Full article
(This article belongs to the Special Issue New Insights into Peritoneal Dialysis and Hemodialysis: 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Serum albumin at different timepoints. Timepoint 0 refers to the first observation, and subsequent timepoints (1–3) refer to month from the first observation.</p>
Full article ">
21 pages, 10727 KiB  
Article
Co-Combustion of Coal and Biomass: Heating Surface Slagging and Flue Gases
by Andrey Zhuikov, Dmitrii Glushkov, Andrey Pleshko, Irina Grishina and Stanislav Chicherin
Fire 2025, 8(3), 106; https://doi.org/10.3390/fire8030106 - 7 Mar 2025
Viewed by 238
Abstract
An experimental study was carried out on the ignition and combustion processes of particles (100–200 µm in size) of coals of different degrees of metamorphism and biomass, as well as mixtures based on them, under conditions of conductive and convective heating, which correspond [...] Read more.
An experimental study was carried out on the ignition and combustion processes of particles (100–200 µm in size) of coals of different degrees of metamorphism and biomass, as well as mixtures based on them, under conditions of conductive and convective heating, which correspond to the conditions of fuel ignition in boiler furnaces at grates and flaring combustion. The biomass contents in the composition of the coal-based fuel mixtures were 10, 20, and 30 wt.%. Under the conductive (at 700–1000 °C) and convective (at 500–800 °C) heating of fuel particles, ignition delay times were determined using a hardware–software complex for the high-speed video registration of fast processes. The ignition delay times were found to vary from 1 to 12.2 s for conductive heating and from 0.01 to 0.19 s for convective heating. The addition of 10–30 wt.% biomass to coals reduced the ignition delay times of fuel mixtures by up to 70%. An analysis of the flue gas composition during the combustion of solid fuels allowed us to establish the concentrations of the main anthropogenic emissions. The use of biomass as an additive (from 10 to 230 wt.%) to coal reduced NOx and SOx emissions by 19–42% and 24–39%, respectively. The propensity of fuels to cause slagging depending on their component composition was established. The use of up to 30 wt.% of biomass in the mixture composition did not affect the increase in the tendency to cause slagging on heating surfaces in the boiler furnace and did not pose a threat to the layer agglomeration during the layer combustion of the mixtures. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the experimental setup for fuel combustion under conditions of conductive heating.</p>
Full article ">Figure 2
<p>A schematic diagram of the experimental setup for fuel combustion under convective heating conditions.</p>
Full article ">Figure 3
<p>A schematic diagram of the experimental setup for studying the composition of flue gases during the combustion of fuels.</p>
Full article ">Figure 4
<p>SEM images of investigated fuels: (<b>a</b>) lignite particles; (<b>b</b>) bituminous coal particles; (<b>c</b>) biomass particles.</p>
Full article ">Figure 5
<p>Video frames of videograms of ignition and combustion of lignite particles (composition Cb) under conditions of conductive heating at <span class="html-italic">T</span><sub>g</sub> = 800 °C.</p>
Full article ">Figure 6
<p>Video frames of videograms of ignition and combustion of bituminous coal particles (composition Ch) under conditions of conductive heating at <span class="html-italic">T</span><sub>g</sub> = 800 °C.</p>
Full article ">Figure 7
<p>Video frames of videograms of ignition and combustion of biomass particles (composition B) under conditions of conductive heating at <span class="html-italic">T<sub>g</sub></span> = 800 °C.</p>
Full article ">Figure 8
<p>Dependences of ignition delay times of coal, biomass, and their mixtures on oxidizer temperature under conductive heating conditions: (<b>a</b>) brown coal (Cb), biomass (B), and their mixtures (CbB-1, CbB-2, and CbB-3); (<b>b</b>) bituminous coal (Ch), biomass (B), and their mixtures (ChB-1, ChB-2, and ChB-3).</p>
Full article ">Figure 9
<p>Video frames of videograms of ignition and combustion of lignite particles (composition Cb) under conditions of movement in flow of heated air at <span class="html-italic">T</span><sub>a</sub> = 700 °C (Δ<span class="html-italic">t</span> = 0.01 s): (<b>a</b>) <span class="html-italic">t</span><sub>d</sub> = 0.052 s; (<b>b</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + ∆<span class="html-italic">t</span>; (<b>c</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + 2∆<span class="html-italic">t</span>; (<b>d</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + 3∆<span class="html-italic">t</span>.</p>
Full article ">Figure 10
<p>Video frames of videograms of ignition and combustion of bituminous coal particles (composition Ch) under conditions of motion in flow of heated air at <span class="html-italic">T<sub>a</sub></span> = 700 °C (Δ<span class="html-italic">t</span> = 0.01 s): (<b>a</b>) <span class="html-italic">t</span><sub>d</sub> = 0.095 s; (<b>b</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + ∆<span class="html-italic">t</span>; (<b>c</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + 2∆<span class="html-italic">t</span>; (<b>d</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + 3∆<span class="html-italic">t</span>.</p>
Full article ">Figure 11
<p>Video frames of videograms of ignition and combustion of biomass particles (composition B) under conditions of motion in flow of heated air at <span class="html-italic">T<sub>a</sub></span> = 700 °C (Δ<span class="html-italic">t</span> = 0.01 s): (<b>a</b>) <span class="html-italic">t</span><sub>d</sub> = 0.030 s; (<b>b</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + ∆<span class="html-italic">t</span>; (<b>c</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + 2∆<span class="html-italic">t</span>; (<b>d</b>) <span class="html-italic">t</span> = <span class="html-italic">t</span><sub>d</sub> + 3∆<span class="html-italic">t</span>.</p>
Full article ">Figure 12
<p>Dependences of ignition delay times of coal, biomass, and their mixtures on oxidant temperature under convective heating conditions: (<b>a</b>) brown coal (Cb), biomass (B), and their mixtures (CbB-1, CbB-2, and CbB-3); (<b>b</b>) bituminous coal (Ch), biomass (B), and their mixtures (ChB-1, ChB-2, and ChB-3).</p>
Full article ">Figure 13
<p>Composition of flue gases during combustion of lignite (Cb), biomass (B), and their mixtures (CbB-1, CbB-2, and CbB-3) at 800 °C: (<b>a</b>) CO and CO<sub>2</sub>; (<b>b</b>) NO<sub>x</sub> and SO<sub>x</sub>.</p>
Full article ">Figure 14
<p>Composition of flue gases during combustion of bituminous coal (Ch), biomass (B), and their mixtures (CbB-1, CbB-2, and CbB-3) at 800 °C: (<b>a</b>) CO and CO<sub>2</sub>; (<b>b</b>) NO<sub>x</sub> and SO<sub>x</sub>.</p>
Full article ">Figure 15
<p>Boiler furnace for layer combustion of solid fuel mixtures: (<b>a</b>) with one fuel bunker; (<b>b</b>) with two fuel bunkers.</p>
Full article ">Figure 16
<p>Boiler furnace for flaring combustion of solid fuel mixture: (<b>a</b>) with single common fuel bunker; (<b>b</b>) with two fuel bunkers.</p>
Full article ">Figure 17
<p>Flare combustion boiler furnace equipped with down draft nozzle.</p>
Full article ">
10 pages, 426 KiB  
Review
The Blue Supergiant Problem and the Main-Sequence Width
by Jorick S. Vink and Rene D. Oudmaijer
Galaxies 2025, 13(2), 19; https://doi.org/10.3390/galaxies13020019 - 6 Mar 2025
Viewed by 147
Abstract
Using Gaia DR3 we derive new distances and luminosities for a sample of Galactic B supergiants which were thought to be post main-sequence (MS) objects from their HR diagram location beyond the terminal-age MS (TAMS). When applying the newer Gaia distances in addition [...] Read more.
Using Gaia DR3 we derive new distances and luminosities for a sample of Galactic B supergiants which were thought to be post main-sequence (MS) objects from their HR diagram location beyond the terminal-age MS (TAMS). When applying the newer Gaia distances in addition to enhanced amounts of core-boundary mixing, aka convective overshooting, we show that these Galactic B supergiants are likely enclosed within the MS band, indicating an evolutionary stage of steady core hydrogen burning. We discuss the importance of considering enhanced overshooting and how vectors in the mass-luminosity plane (ML-plane) can be used to disentangle the effects of wind mass loss from interior mixing. We finish with the key message that any proposed solution to the BSG problem should consider not only an explanation for the sheer number of B supergiants inside the Hertzsprung gap, but should at the same time also account for the steep drop in rotation rates identified at spectral type B1—corresponding to an effective temperature of ∼21 kK, and for which two distinct families of solutions have been proposed. Full article
(This article belongs to the Special Issue Circumstellar Matter in Hot Star Systems)
Show Figures

Figure 1

Figure 1
<p>HR diagram comparison of Galactic B supergiants with solar metallicity [<a href="#B34-galaxies-13-00019" class="html-bibr">34</a>] MESA model tracks with initial masses from 8–60 <math display="inline"><semantics> <msub> <mi mathvariant="normal">M</mi> <mo>⊙</mo> </msub> </semantics></math>, in steps 8, 12, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60 <math display="inline"><semantics> <msub> <mi mathvariant="normal">M</mi> <mo>⊙</mo> </msub> </semantics></math>. The figure on the left-hand side is for a small amount of overshooting with <math display="inline"><semantics> <msub> <mi>α</mi> <mi>ov</mi> </msub> </semantics></math>= 0.1, while the figure on the right-hand side is for <math display="inline"><semantics> <msub> <mi>α</mi> <mi>ov</mi> </msub> </semantics></math>= 0.5. The dashed lines on both the left and right-hand side denote the TAMS location. The blue stars represent the observed B supergiant sample from Crowther et al. [<a href="#B22-galaxies-13-00019" class="html-bibr">22</a>] with our updated luminosities utilising Gaia DR3 distances.</p>
Full article ">Figure 2
<p>An illustration of the mass-luminosity plane, showing a typical evolutionary track that begins at the ZAMS, indicated by the red dot, and progresses along the black arrow towards the TAMS. The dotted vector indicates how factors such as increased rotation and/or convective overshooting can extend the M-L vector. The curved dashed line represents the gradient where mass-loss rates influence this M-L vector. The red solid area is prohibited, as set by the mass-luminosity relationship. Adapted from Higgins &amp; Vink [<a href="#B34-galaxies-13-00019" class="html-bibr">34</a>].</p>
Full article ">Figure 3
<p>Luminosities (Left-hand side) and rotational velocities (right-hand side) versus effective temperature from the VLT Flames survey of massive stars for evolutionary masses over 15 <math display="inline"><semantics> <msub> <mi mathvariant="normal">M</mi> <mo>⊙</mo> </msub> </semantics></math>. Luminosity classes are represented by blue plus signs (for luminosity classes II–V) and red stars (for luminosity class I). The evolutionary tracks for LMC metallicity (50% solar), including the predicted bi-stability jump, are shown in grey, with initial rotation velocities of 250 km/s for five masses: 15, 20, 30, 40, and 60 <math display="inline"><semantics> <msub> <mi mathvariant="normal">M</mi> <mo>⊙</mo> </msub> </semantics></math>. Note that the critical mass for bi-stability braking is approximately 35 <math display="inline"><semantics> <msub> <mi mathvariant="normal">M</mi> <mo>⊙</mo> </msub> </semantics></math> at LMC metallicity in these specific models computed with the BONN stellar evolution code for a step overshooting value of <math display="inline"><semantics> <msub> <mi>α</mi> <mi>ov</mi> </msub> </semantics></math> = 0.5. The rotational velocity tracks can be compared to the angular momentum conservation case, depicted as grey dotted background lines, showing that bi-stability braking can be steeper than would be attributable to angular momentum conservation. The black dots along the evolutionary tracks correspond to <math display="inline"><semantics> <msup> <mn>10</mn> <mn>5</mn> </msup> </semantics></math> year time-steps. Adapted from Vink et al. [<a href="#B20-galaxies-13-00019" class="html-bibr">20</a>].</p>
Full article ">
24 pages, 6847 KiB  
Article
Comparing Reflectivity from Space-Based and Ground-Based Radars During Detection of Rainbands in Two Tropical Cyclones
by Corene J. Matyas, Stephanie E. Zick and Kimberly M. Wood
Atmosphere 2025, 16(3), 307; https://doi.org/10.3390/atmos16030307 - 6 Mar 2025
Viewed by 78
Abstract
With varying tangential winds and combinations of stratiform and convective clouds, tropical cyclones (TCs) can be difficult to accurately portray when mosaicking data from ground-based radars. This study utilizes the Dual-frequency Precipitation Radar (DPR) from the Global Precipitation Measurement Mission (GPM) satellite to [...] Read more.
With varying tangential winds and combinations of stratiform and convective clouds, tropical cyclones (TCs) can be difficult to accurately portray when mosaicking data from ground-based radars. This study utilizes the Dual-frequency Precipitation Radar (DPR) from the Global Precipitation Measurement Mission (GPM) satellite to evaluate reflectivity obtained using four sampling methods of Weather Surveillance Radar 1988-Doppler data, including ground radars (GRs) in the GPM ground validation network and three mosaics, specifically the Multi-Radar/Multi-Sensor System plus two we created by retaining the maximum value in each grid cell (MAX) and using a distance-weighted function (DW). We analyzed Hurricane Laura (2020), with a strong gradient in tangential winds, and Tropical Storm Isaias (2020), where more stratiform precipitation was present. Differences between DPR and GR reflectivity were larger compared to previous studies that did not focus on TCs. Retaining the maximum value produced higher values than other sampling methods, and these values were closest to DPR. However, some MAX values were too high when DPR time offsets were greater than 120 s. The MAX method produces a more consistent match to DPR than the other mosaics when reflectivity is <35 dBZ. However, even MAX values are 3–4 dBZ lower than DPR in higher-reflectivity regions where gradients are stronger and features change quickly. The DW and MRMS mosaics produced values that were similar to one another but lower than DPR and MAX values. Full article
Show Figures

Figure 1

Figure 1
<p>Reflectivity values produced during the MAX sampling method (maximum value is retained) (shading) and circles denoting locations where data are available for Dual-Frequency Precipitation Radar (DPR) and Ground Radar (GR) in the Global Precipitation Measurement Mission Ground Validation System Validation Network (GPM GVS VN) database for (<b>a</b>) Laura and (<b>b</b>) Isaias.</p>
Full article ">Figure 2
<p>Hurricane Laura reflectivity data from sampling methods (<b>a</b>) DPR, (<b>b</b>) MAX, (<b>c</b>) GR, (<b>d</b>) MRMS, and (<b>e</b>) DW, and the difference of each sampling method subtracted from DPR (<b>f</b>) MAX, (<b>g</b>) GR, (<b>h</b>) MRMS, and (<b>i</b>) DW.</p>
Full article ">Figure 3
<p>Tropical Storm Isaias reflectivity data from sampling methods (<b>a</b>) DPR, (<b>b</b>) MAX, (<b>c</b>) GR, (<b>d</b>) MRMS, and (<b>e</b>) DW, and the difference of each sampling method subtracted from DPR (<b>f</b>) MAX, (<b>g</b>) GR, (<b>h</b>) MRMS, and (<b>i</b>) DW.</p>
Full article ">Figure 4
<p>Flowchart showing the procedure for processing reflectivity from a single radar and then placing those data on the TC-centric grid. The “identify final value” step varies by approach, and the same approach is used for all single-station grids and the TC-centered grid: it either computes the distance-weighted average for the DW method or retains the maximum value for the MAX method.</p>
Full article ">Figure 5
<p>Distribution of reflectivity at each primary radar from each sampling method. Circles and asterisks denote minimum and maximum values that are outliers, whiskers denote minimum and maximum values that are not outliers, and boxes show lower and upper quartiles and median reflectivity.</p>
Full article ">Figure 6
<p>For each primary radar, each WSR sampling method is subtracted from DPR, and these statistics are portrayed as follows: (<b>a</b>) mean error, (<b>b</b>) standard deviation of the mean error, and (<b>c</b>) mean absolute error. (<b>d</b>) Correlation coefficients calculated between DPR and each WSR sampling method.</p>
Full article ">Figure 7
<p>The percent that each radar in the MAX mosaic contributed to the points within 100 km of each primary radar. The locations of the contributing or secondary radars are visible in <a href="#atmosphere-16-00307-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 8
<p>Map of difference in reflectivity between DPR and MAX mosaic created with radars (<b>a</b>) KHGX, KLCH, KPOE (as in <a href="#atmosphere-16-00307-f002" class="html-fig">Figure 2</a>f), and (<b>b</b>) KLCH and KPOE only.</p>
Full article ">Figure 9
<p>Difference in reflectivity between DPR and each WSR sampling method (DPR-WSR) with values binned according to DPR reflectivity. Boxplot bars, whiskers, asterisks and circles are as in <a href="#atmosphere-16-00307-f005" class="html-fig">Figure 5</a>.</p>
Full article ">
23 pages, 17262 KiB  
Review
Research Progress on Solar Supergranulation: Observations, Theories, and Numerical Simulations
by Chong Huang and Rui Wang
Universe 2025, 11(3), 87; https://doi.org/10.3390/universe11030087 - 6 Mar 2025
Viewed by 75
Abstract
Solar supergranulation is a large-scale convective structure on the solar surface, whose formation mechanism and dynamical properties are closely related to key physical processes such as solar magnetic field evolution, coronal heating, and solar wind acceleration. This paper reviews recent research progress on [...] Read more.
Solar supergranulation is a large-scale convective structure on the solar surface, whose formation mechanism and dynamical properties are closely related to key physical processes such as solar magnetic field evolution, coronal heating, and solar wind acceleration. This paper reviews recent research progress on solar supergranulation, focusing on the latest achievements in high-resolution observations, theoretical models, and numerical simulations. By analyzing the flow field structure, magnetic field distribution, and their relationship with the solar activity cycle, the crucial role of supergranulation in solar physics is revealed. Studies indicate that supergranulation is not only a crucial component of the solar convection zone but also drives coronal heating and solar wind acceleration through mechanisms such as magnetic reconnection and Alfvén wave propagation. Furthermore, the interaction between supergranulation and larger-scale convective patterns (e.g., giant cells) provides new insights into the dynamics of the solar interior. Despite significant progress in recent years, the formation mechanism and dynamical nature of supergranulation remain unresolved. Future research should combine high-resolution observations, theoretical modeling, and numerical simulations to further elucidate the complex dynamical processes and the central role of supergranulation in solar physics. Full article
Show Figures

Figure 1

Figure 1
<p>Horizontal velocity field of supergranulation derived from granule tracking using the CALAS camera at the Pic-du-Midi Observatory. Bottom right: CALAS field of view on the Sun. The large rectangle indicates the full image size, the smaller rectangle marks the region where velocities were computed, and the small square represents the field of view of the Hinode/SOT instrument for granule tracking. Top: Supergranulation velocity field with divergence contours overlaid. Arrows indicate the direction and magnitude of the horizontal velocity field obtained from granule tracking. Bottom left: Zoomed-in view of granulation, illustrating the relative scales of granules and supergranules [<a href="#B4-universe-11-00087" class="html-bibr">4</a>].</p>
Full article ">Figure 2
<p>Irregular distribution of magnetic fields at supergranulation boundaries [<a href="#B5-universe-11-00087" class="html-bibr">5</a>].</p>
Full article ">Figure 3
<p>Supergranulation boundaries are overlaid on Doppler velocity maps for mid-latitude (left) and polar (right) regions. The top images (<b>A</b>) were generated using a moments technique applied to unbinned data, with off-limb solar observations serving as the “zero” velocity reference. The bottom images (<b>B</b>) were produced via Gaussian fitting of binned data (to enhance the signal-to-noise ratio) using the “zero velocity” wavelength reference. Blue denotes blueshifts or outflows, while red indicates redshifts or downflows. Each panel has a spatial scale of approximately 540<sup>″</sup> × 300<sup>″</sup>, and the Doppler velocity maps illustrate flow patterns on the solar surface [<a href="#B6-universe-11-00087" class="html-bibr">6</a>].</p>
Full article ">Figure 4
<p>Influence of the Coriolis force on supergranulation flows. (<b>a</b>) Vertical vorticity (curl) averaged across areas with positive horizontal divergence (〈curl<math display="inline"><semantics> <msub> <mrow> <mo>〉</mo> </mrow> <mo>+</mo> </msub> </semantics></math>, depicted by the blue curve) and negative horizontal divergence (〈curl<math display="inline"><semantics> <msub> <mrow> <mo>〉</mo> </mrow> <mo>−</mo> </msub> </semantics></math>, shown by the red curve), plotted as functions of sin(<math display="inline"><semantics> <mi>λ</mi> </semantics></math>)<math display="inline"><semantics> <mo>Ω</mo> </semantics></math>(<math display="inline"><semantics> <mi>λ</mi> </semantics></math>)/<math display="inline"><semantics> <msub> <mo>Ω</mo> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> </semantics></math>, where <math display="inline"><semantics> <mi>λ</mi> </semantics></math> represents the heliographic latitude and <math display="inline"><semantics> <mo>Ω</mo> </semantics></math>/<math display="inline"><semantics> <msub> <mo>Ω</mo> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> </semantics></math> is the local surface angular velocity relative to the equator. A vorticity value of 10<sup>−6</sup> s<sup>−1</sup> corresponds to an angular velocity of 2.5° per day or a typical tangential velocity of 10 m s<sup>−1</sup>. (<b>b</b>) Horizontally averaged 〈curl div〉 as a function of sin(<math display="inline"><semantics> <mi>λ</mi> </semantics></math>)<math display="inline"><semantics> <mo>Ω</mo> </semantics></math>(<math display="inline"><semantics> <mi>λ</mi> </semantics></math>)/<math display="inline"><semantics> <msub> <mo>Ω</mo> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> </semantics></math> [<a href="#B31-universe-11-00087" class="html-bibr">31</a>].</p>
Full article ">Figure 5
<p>Variation in the horizontal and vertical field scales of supergranulation with latitude [<a href="#B55-universe-11-00087" class="html-bibr">55</a>].</p>
Full article ">Figure 6
<p>Snapshots of plume positions from n-body advective-interaction models with differing parameter values. The rate of new granular plume creation (<math display="inline"><semantics> <msub> <mi mathvariant="normal">N</mi> <mi>n</mi> </msub> </semantics></math>) increases from left to right, and the advection profile width (<math display="inline"><semantics> <msub> <mi>σ</mi> <mi>f</mi> </msub> </semantics></math>) increases from bottom to top. In each sub-figure, the x and y axes represent spatial dimensions in units of simulation grid points (pixels) [<a href="#B61-universe-11-00087" class="html-bibr">61</a>].</p>
Full article ">Figure 7
<p>Flow velocity maps of solar supergranules. This sequence of Mercator projection maps illustrates the longitudinal and latitudinal flow velocities of solar supergranules observed from May to August 2010. The maps reveal large-scale flow structures, with prograde and southward velocities shown in red and retrograde and northward velocities in blue. The data indicate the existence of giant convection cells that play a crucial role in maintaining the Sun’s differential rotation [<a href="#B63-universe-11-00087" class="html-bibr">63</a>].</p>
Full article ">Figure 8
<p>Flow structure in a high-resolution simulation of solar convection. The radial velocity is shown near (<b>a</b>) the top and (<b>b</b>) deeper within the convection zone (r = 0.98 R and r = 0.92 R, respectively). Yellow/white tones denote upflow, and blue/black tones denote downflow. Each spherical horizontal surface is shown in a Mollweide projection, with dashed lines indicating latitudes of 0°, ±30°, and ±60°, and longitudes of 0° and ±90°. The color table is linear, and peak velocity amplitudes are indicated by the color bar [<a href="#B64-universe-11-00087" class="html-bibr">64</a>].</p>
Full article ">Figure 9
<p>Maps of Intensity and flow fields at the surface of the solar convection zone, derived from local MHD simulations. The figure presents a comparison between the kinematic growth phase (panels (<b>a</b>–<b>c</b>)) and the saturated phase (panels (<b>d</b>–<b>f</b>)) of the simulation run [<a href="#B19-universe-11-00087" class="html-bibr">19</a>].</p>
Full article ">
23 pages, 29026 KiB  
Article
Urban Impacts on Convective Squall Lines over Chicago in the Warm Season—Part I: Observations of Multi-Scale Convective Evolution
by Michael L. Kaplan, S. M. Shajedul Karim and Yuh-Lang Lin
Atmosphere 2025, 16(3), 306; https://doi.org/10.3390/atmos16030306 - 6 Mar 2025
Viewed by 169
Abstract
In this study, our aim is to diagnose how two quasi-linear convective systems (QLCS) are organized so one can determine the possible role of the city of Chicago, IL, USA, in modifying convective precipitation systems. In this Part I of a two-part study, [...] Read more.
In this study, our aim is to diagnose how two quasi-linear convective systems (QLCS) are organized so one can determine the possible role of the city of Chicago, IL, USA, in modifying convective precipitation systems. In this Part I of a two-part study, we employ large-scale analyses, radiosonde soundings, surface observations, and Doppler radar data to diagnose the precursor atmospheric circulations that organize the evolution of two mesoscale convective systems and compare those circulations to radar and precipitation. Several multi-scale processes are found that organize and modify convection over the Chicago metroplex. Two sequential quasi-linear convective systems (QLCS #1 and #2) were organized that propagated over Chicago, IL, USA, during an eight-hour period on 5–6 July 2018. The first squall line (QLCS #1) built from the southwest to the northeast while strengthening as it propagated over the city, and the second (QLCS #2) propagated southeastwards and weakened as it passed over the city in association with a polar cold front. The weak upper-level divergence associated with a diffluent flow poleward of an expansive ridge built over and strengthened a low-level trough and confluence zone, triggering QLCS #1. Convective downdrafts from QLCS #1 produced a cold pool that interacted with multiple confluent low-level jets surrounding and focused on the metroplex urban heat island, thus advecting the convection poleward over the metroplex. The heaviest precipitation occurred just south-southeast of Midway Airport, Chicago. Subsequently, a polar cold front propagated into the metroplex, which triggered QLCS #2. However, the descending air above it under the polar jet and residual cold pool from QLCS #1 rapidly dissipated the cold frontal convection. This represents a case study where very active convection built over the metroplex and was likely modified by it, as evidenced in numerical simulations to be described in Part II. Full article
(This article belongs to the Section Meteorology)
Show Figures

Figure 1

Figure 1
<p>Study area focuses on Chicago, IL, USA, and its adjacent states. For analysis, focus is mainly on Chicago metropolitan region, annotated by a green rectangle box. Analysis of meso-γ scale observations included National Weather Service (NWS) stations, namely Romeoville, IL (KLOT), Midway Airport (KMDW), O’Hare Airport (KORD), Wheeling, IL (KPWK), DuPage County, IL (KDPA), Rockford Airport (KRFD), Lincoln, IL (KILX), and Gary, IN (KGYY). Shaded red color indicates urban area. Black dots indicate major adjacent urban cities.</p>
Full article ">Figure 2
<p>RUC one-hour forecasts of 300 hPa streamlines (black solid arrows), 400–250 hPa isobaric potential vorticity (IPV; color fill, K hPa<sup>−1</sup> s<sup>−1</sup>), and IPV advection (positive blue solid, negative red dashed) valid at (<b>a</b>) 00007/4, (<b>b</b>) 12007/4, (<b>c</b>) 00007/5, and (<b>d</b>) 12007/5. Red arrow in (<b>d</b>) indicates the area of split streamlines south of Chicago, IL, USA.</p>
Full article ">Figure 3
<p>RUC one-hour forecasts of total column precipitable water (black solid &lt; 1 inch, green fill &gt; 1 inch), 850 hPa moisture transport (red vectors, color fill &gt; 10 g kg<sup>−1</sup> s<sup>−1</sup>), height (black solid, m), and equivalent potential temperature (green dashed, K), respectively valid at (<b>a</b>,<b>b</b>) 00007/4, (<b>c</b>,<b>d</b>) 12007/4, (<b>e</b>,<b>f</b>) 00007/5, and (<b>g</b>,<b>h</b>) 12007/5. Red oval-shaped area in (<b>e</b>,<b>f</b>) highlights eastward surge of moisture.</p>
Full article ">Figure 4
<p>Zoomed in area of RUC analyses of (<b>a</b>,<b>c</b>) 3-h (0900–1200 UTC 5 July 2018) mixed layer CAPE, MUCAPE change, and SBCAPE change (thin red line for smaller values and thick red line for larger values in J kg<sup>−1</sup>), CIN (convective inhibition = blue fill in J kg<sup>−1</sup>) and in (<b>b</b>) MUCAPE change, LCL height (parcel saturation level = green hatched in m AGL) and (<b>d</b>) Lincoln, Illinois (KILX) sounding valid at 12007/5. Star indicates location of Chicago, IL, USA, and square indicates location of Lincoln, IL, USA. Green oval-shaped area in (<b>c</b>) highlights area of significant SBCAPE &gt; 2000 J kg<sup>−1</sup> southwest of Chicago, IL, USA.</p>
Full article ">Figure 5
<p>Zoomed in area of RUC analysis of 250 hPa velocity divergence (purple solid, 10<sup>−5</sup> s<sup>−1</sup>), 850 hPa velocity convergence (black or red dashed, 10<sup>−5</sup> s<sup>−1</sup>), and vertically differential velocity divergence (green or blue fill, 10<sup>−5</sup> s<sup>−1</sup>) valid at (<b>a</b>) 12007/5, (<b>b</b>) 14007/5, (<b>c</b>) 16007/5, (<b>d</b>) 18007/5, (<b>e</b>) 19007/5, and (<b>f</b>) 20007/5; and 300 hPa ageostrophic wind barbs (kt) and 700–500 hPa omega (magenta solid upward/red dashed downward, μbs<sup>−1</sup>) valid at (<b>g</b>) 12007/5, (<b>h</b>) 14007/5, (<b>i</b>) 16007/5, (<b>j</b>) 18007/5, (<b>k</b>) 19007/5, and (<b>l</b>) 20007/5.</p>
Full article ">Figure 6
<p>Zoomed in area of RUC analysis of (<b>a</b>) 250 hPa wind barbs (kt) and velocity tensor magnitude (red solid, 10<sup>−4</sup> s<sup>−1</sup>) and (<b>b</b>) 300 hPa wind barbs (kt), height (black solid, m), isotachs (color fill &gt; 60 kt) and divergence (magenta solid, 10<sup>−4</sup> s<sup>−1</sup>) valid at 19007/5. Star indicates the location of Chicago, IL, USA. Blue circle in (<b>a</b>) highlights area of wind barb speed variation just southwest of Chicago, IL, USA.</p>
Full article ">Figure 7
<p>Zoomed in area of RUC analysis of (<b>a</b>,<b>b</b>) 925 hPa wind barbs (kt), height (black solid, m), isotachs (color fill &gt; 60 kt), temperature (red dashed, °C), and dewpoint temperature (green solid, fill &gt; 14 °C) valid at 00007/6 and 02007/6. (<b>c</b>,<b>d</b>) 700 hPa wind barbs (kt), height (solid black in m), mean relative humidity (green fill &gt; 70%), and temperature (dotted red line in °C) valid at 00007/6 and 02007/6. (<b>e</b>,<b>f</b>) and 300 hPa ageostrophic wind barbs (kt), isotach (blue fill &gt; 60 kt), and 700–500 hPa omega (magenta solid upward/red dashed downward, μbs<sup>−1</sup>) valid at 00007/6 and 02007/6. Star indicates location of Chicago, IL, USA.</p>
Full article ">Figure 8
<p>Mean sea level pressure (MSLP) tendency [hPa (3 h<sup>−1</sup>)] observations valid at (<b>a</b>) 15007/5, (<b>b</b>) 18007/5, and (<b>c</b>) 21007/5. Surface streamline analyses are valid at (<b>d</b>) 15007/5, (<b>e</b>) 18007/5, and (<b>f</b>) 21007/5. L and H, just south-southwest of Chicago, signify meso-low and meso-high transitions, respectively, consistent with evolving surface confluence (conf) and difluence (diff) from northwestern Indiana to Chicago. Star indicates location of Chicago, IL, USA.</p>
Full article ">Figure 9
<p>Surface temperature (°C) observations valid at (<b>a</b>) 17007/5, (<b>b</b>) 18007/5, (<b>c</b>) 19007/5, and (<b>d</b>) 20007/5. Star indicates location of Chicago, IL, USA.</p>
Full article ">Figure 10
<p>Zoomed in area of RUC analysis of the lowest 100 hPa three-hour mean mixing ratio change (color fill, g kg<sup>−1</sup> [3 h]<sup>−1</sup>) and surface wind barbs (kt) valid at (<b>a</b>) 12007/5 and (<b>b</b>) 14007/5. (<b>c</b>) Pressure (MSLP, black solid, hPa) and MSLP tendency (red solid positive, dashed blue negative, hPa [2 h]<sup>−1</sup>) valid at 16007/5. (<b>d</b>) Surface wind barbs (kt), surface vorticity (color fill, 10<sup>−5</sup> s<sup>−1</sup>), and surface divergence (red solid positive, blue dashed negative, 10<sup>−5</sup> s<sup>−1</sup>) valid at 16007/5. RUC analysis of three-hour SBCAPE change (red solid positive, blue dashed negative, J kg<sup>−1</sup>) and SBCIN (blue fill, J kg<sup>−1</sup>) valid at (<b>e</b>) 12007/5, (<b>f</b>) 16007/5, (<b>g</b>) 18007/5, and (<b>h</b>) 22007/5. Star indicates location of Chicago, IL, USA. Red oval shape in (<b>d</b>) highlights area of maximum mixing ratio change just southeast of Chicago, IL, USA. Black contour in (<b>c</b>) just southwest of Chicago indicates the trough.</p>
Full article ">Figure 11
<p>Observed Plymouth State Weather Center’s surface meteograms at (<b>a</b>) Romeoville, Illinois (KLOT), (<b>b</b>) Midway Airport (KMDW), (<b>c</b>) O’Hare Airport (KORD), (<b>d</b>) Wheeling, Illinois (KPWK), (<b>e</b>) DuPage County, Illinois (KDPA), and (<b>f</b>) Gary, Indiana (KGYY) for period 0000–2300 UTC 5 July 2018. Temperature drop (sky oval shape) and pressure rise (red oval shape) at KMDW (<b>b</b>) and KGYY (<b>f</b>) between 1900 and 2000 UTC.</p>
Full article ">Figure 12
<p>NOAA NEXRAD Doppler composite reflectivity (dBz) valid at (<b>a</b>) 1700 UTC, (<b>b</b>) 1800 UTC, (<b>c</b>) 1810 UTC, (<b>d</b>) 1820 UTC, (<b>e</b>) 1840 UTC, (<b>f</b>) 1845 UTC, (<b>g</b>) 1855 UTC, (<b>h</b>) 1905 UTC, (<b>i</b>) 1930 UTC, (<b>j</b>) 2000 UTC, (<b>k</b>) 2100 UTC, (<b>l</b>) 2200 UTC, (<b>m</b>) 2300 UTC, (<b>n</b>) 2330 UTC 5 July 2018, and (<b>o</b>) 0000 UTC 6 July 2018. Purple color annotation indicates evolution of QLCS #1, and red color annotation indicates evolution of QLCS#2.</p>
Full article ">Figure 13
<p>NWS Chicago/Romeoville (KLOT) one-hour Doppler-derived rainfall totals (in) valid at approximately (<b>a</b>) 18307/5, (<b>b</b>) 19007/5, (<b>c</b>) 20007/5, (<b>d</b>) 21007/5, (<b>e</b>) 22007/5, (<b>f</b>) 23007/5, (<b>g</b>) 00007/6 at KLOT. For this event, precipitation scale for hourly rates in inches roughly ranges from 0.25–0.50 light blue, 0.50–0.75 green, 0.75–1.25 magenta, and 1.25–1.75 dark blue.</p>
Full article ">Figure 13 Cont.
<p>NWS Chicago/Romeoville (KLOT) one-hour Doppler-derived rainfall totals (in) valid at approximately (<b>a</b>) 18307/5, (<b>b</b>) 19007/5, (<b>c</b>) 20007/5, (<b>d</b>) 21007/5, (<b>e</b>) 22007/5, (<b>f</b>) 23007/5, (<b>g</b>) 00007/6 at KLOT. For this event, precipitation scale for hourly rates in inches roughly ranges from 0.25–0.50 light blue, 0.50–0.75 green, 0.75–1.25 magenta, and 1.25–1.75 dark blue.</p>
Full article ">
22 pages, 3903 KiB  
Article
Ultrasound-Assisted Extraction, Characterization, and Antioxidant Activities of the Polysaccharides from Fermented Astragalus membranaceus
by Jingyan Zhang, Zijing Liang, Kang Zhang, Xi Tang, Lei Wang, Xueyan Gu, Huub F. J. Savelkoul and Jianxi Li
Molecules 2025, 30(5), 1159; https://doi.org/10.3390/molecules30051159 - 4 Mar 2025
Viewed by 274
Abstract
This study aimed to optimize the ultrasound-assisted extraction (UAE) of polysaccharides from fermented Astragalus membranaceus (FAPS) and to investigate the physicochemical properties and antioxidant activities of the extracted polysaccharides. Using a combination of single-factor experiments and response surface methodology based on a Box–Behnken [...] Read more.
This study aimed to optimize the ultrasound-assisted extraction (UAE) of polysaccharides from fermented Astragalus membranaceus (FAPS) and to investigate the physicochemical properties and antioxidant activities of the extracted polysaccharides. Using a combination of single-factor experiments and response surface methodology based on a Box–Behnken design, we improved the extraction of crude FAPS without deproteinization. Under optimal conditions (50 °C, 60 min, 8 mL/g, 480 W), the yield of crude FAPS obtained by UAE (7.35% ± 0.08) exceeded the yield from convectional hot water extraction (6.95% ± 0.24). After protein removal, the FAPS was subjected to comprehensive chemical analyses, including HPLC, HPGPC, FT-IR, UV spectroscopy, and a Congo red assay. The results showed that FAPS had a significantly higher carbohydrate content compared to the non-fermented group (95.38% ± 6.20% vs. 90.938% ± 3.80%), while the protein content was significantly lower than that of the non-fermented Astragalus polysaccharides (APS) group (1.26% ± 0.34% vs. 6.76% ± 0.87%). In addition, FAPS had a higher average molecular weight and a lower Mw/Mn ratio compared to APS. The primary monosaccharides in FAPS were identified as Glc, Ara, Gal and GalA, with a molar ratio of 379.72:13.26:7.75:6.78, and FAPS lacked a triple helix structure. In vitro, antioxidant assays showed that FAPS possessed superior antioxidant properties compared to APS. These results emphasize the significant potential of FAPS as an antioxidant, possibly superior to that of APS. The results of this study suggest that fermentation and UAE offer promising applications for the development and utilization of Astragalus membranaceus for human and animal health. Full article
Show Figures

Figure 1

Figure 1
<p>Flow chart of UAE and antioxidant activity analysis of FAPS.</p>
Full article ">Figure 2
<p>Effects of different extraction temperatures (<b>A</b>), extraction times (<b>B</b>), ratios of water to material (<b>C</b>) and extraction powers (<b>D</b>) on the yield of crude FAPS (CFAPS).</p>
Full article ">Figure 3
<p>The three dimensional response surface plots and two dimensional contour plots show the interaction effects between water to material ratio and extraction time (<b>A</b>,<b>D</b>), extraction power and extraction time (<b>B</b>,<b>E</b>), extraction power and the ratio of material on the yield of CFAPS (<b>C</b>,<b>F</b>).</p>
Full article ">Figure 4
<p>The HPLC chromatograms of monosaccharide standards (<b>A</b>), APS (<b>B</b>) and FAPS (<b>C</b>).</p>
Full article ">Figure 5
<p>The HPGPC spectra of APS (<b>A</b>) and FAPS (<b>B</b>).</p>
Full article ">Figure 6
<p>Molecular weight distribution of FAPS and APS.</p>
Full article ">Figure 7
<p>FT-IR spectrum of APS (<b>A</b>) and FAPS (<b>B</b>).</p>
Full article ">Figure 8
<p>UV spectrum of FAPS and APS (<b>A</b>). Changes in absorption wavelength maximum of mixture of Congo red, FAPS and APS at various concentrations of NaOH (<b>B</b>).</p>
Full article ">Figure 9
<p>Scavenging effects of APS and FAPS at different concentrations on DPPH radical scavenging assay (<b>A</b>), hydroxyl radical scavenging assay (<b>B</b>), ABTS radical scavenging assay (<b>C</b>) and ferric reducing power assay (<b>D</b>).</p>
Full article ">
19 pages, 2274 KiB  
Article
Construction and Application of a Coupled Temperature and Pressure Model for CO2 Injection Wells Considering Gas Composition
by Hang Lai, Peng Chen, Lingang Lv and Song Lu
Energies 2025, 18(5), 1238; https://doi.org/10.3390/en18051238 - 3 Mar 2025
Viewed by 278
Abstract
Accurate prediction of the temperature and pressure fields in carbon dioxide (CO2) injection wells is critical for enhancing oil recovery efficiency and ensuring safe carbon sequestration. At present, the prediction model generally assumes that CO2 is pure and does not [...] Read more.
Accurate prediction of the temperature and pressure fields in carbon dioxide (CO2) injection wells is critical for enhancing oil recovery efficiency and ensuring safe carbon sequestration. At present, the prediction model generally assumes that CO2 is pure and does not consider the influence of impurities in CO2 components. This study takes into account the common impurities, such as air and various alkanes in CO2, and uses Refprop 9.0 software to calculate the physical parameters of the mixture. A comprehensive coupling model was developed to account for axial heat conduction, convective heat transfer, frictional heat generation, the soup coke effect, pressure work, and gas composition. The model was solved iteratively using numerical methods. We validated the accuracy of the calculated results by comparing our model with the Ramey model using measured injection well data. Compared with the measured bottom hole temperature and pressure data, the error percentage of our model to predict the bottom hole temperature and pressure is less than 1%, while the error percentage of Ramey model to predict the bottom hole temperature and pressure is 5.15% and 1.33%, respectively. Our model has higher bottom hole temperature and pressure prediction accuracy than the Ramey model. In addition, we use the model to simulate the influence of different injection parameters on wellbore temperature and pressure and consider the influence of different gas components. Each injection parameter uses three components. Based on the temperature and pressure data calculated by the model simulation, the phase state of CO2 was analyzed. The results show that the impurities in CO2 have a great influence on the predicted wellbore pressure, critical temperature, and critical pressure. In the process of CO2 injection, increasing the injection pressure can significantly increase the bottom hole pressure, and changing the injection rate can adjust the bottom hole temperature. The research provides valuable insights for CO2 sequestration and enhanced oil recovery (EOR). Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the heat transfer model.</p>
Full article ">Figure 2
<p>Model solution workflow diagram.</p>
Full article ">Figure 3
<p>The measured wellbore temperature and pressure data and the comparison of the prediction results of the two models (our model and Ramey model).</p>
Full article ">Figure 4
<p>The measured wellbore temperature and pressure data and the comparison of the predicted results of our model using two components (100% CO<sub>2</sub> and 90.5% CO<sub>2</sub> + 9.5% N<sub>2</sub>).</p>
Full article ">Figure 5
<p>The temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at different injection temperatures: (<b>a</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection temperature of 20 °C; (<b>b</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection temperature of 30 °C; (<b>c</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection temperature of 40 °C.</p>
Full article ">Figure 6
<p>The temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at different injection pressures: (<b>a</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection pressure of 30 MPa; (<b>b</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection pressure of 40 MPa; (<b>c</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection pressure of 50 MPa.</p>
Full article ">Figure 7
<p>The temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at different injection rates: (<b>a</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection rate of 60 t/d; (<b>b</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection rate of 80 t/d; (<b>c</b>) the temperature distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection rate of 100 t/d.</p>
Full article ">Figure 8
<p>The pressure distribution of the fluid in the tubing under three kinds of CO<sub>2</sub> components at different injection temperatures: (<b>a</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection temperature of 20 °C; (<b>b</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection temperature of 30 °C; (<b>c</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection temperature of 40 °C.</p>
Full article ">Figure 9
<p>The pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at different injection pressures: (<b>a</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection pressure of 30 MPa; (<b>b</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection pressure of 40 MPa; (<b>c</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection pressure of 50 MPa.</p>
Full article ">Figure 10
<p>The pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at different injection rates: (<b>a</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection rate of 60 t/d; (<b>b</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection rate of 80 t/d; (<b>c</b>) the pressure distribution of the fluid in the tubing with three kinds of CO<sub>2</sub> components at an injection rate of 100 t/d.</p>
Full article ">
17 pages, 5812 KiB  
Article
Significance of Cloud Microphysics and Cumulus Parameterization Schemes in Simulating an Extreme Flood-Producing Precipitation Event in the Central Himalaya
by Ujjwal Tiwari and Andrew B. G. Bush
Atmosphere 2025, 16(3), 298; https://doi.org/10.3390/atmos16030298 - 3 Mar 2025
Viewed by 164
Abstract
Between 11 and 14 August 2017, the southern belt of the central Himalaya experienced extreme precipitation, with some stations recording more than 500 mm of accumulated rainfall, which resulted in widespread, devastating flooding. Precipitation was concentrated over the sub-Himalaya, and the established forecasting [...] Read more.
Between 11 and 14 August 2017, the southern belt of the central Himalaya experienced extreme precipitation, with some stations recording more than 500 mm of accumulated rainfall, which resulted in widespread, devastating flooding. Precipitation was concentrated over the sub-Himalaya, and the established forecasting systems failed to predict the event. In this study, we evaluate the performance of six cloud microphysics schemes in the Weather Research and Forecasting (WRF) model forced with the advanced ERA5 dataset. We also examine the importance of the cumulus scheme in WRF at 3 km horizontal grid spacing in highly convective events like this. Six WRF simulations, each with one of the six different microphysics schemes with the Kain–Fritsch cumulus scheme turned off, all fail to reproduce the spatial variability of accumulated precipitation during this devastating flood-producing precipitation event. In contrast, the simulations exhibit greatly improved performance with the cumulus scheme turned on. In this study, the cumulus scheme helps to initiate convection, after which grid-scale precipitation becomes dominant. Amongst the different simulations, the WRF simulation using the Morrison microphysics scheme with the cumulus turned on displayed the best performance, with the smallest normalized root mean square error (NRMSE) of 0.25 and percentage bias (PBIAS) of −6.99%. The analysis of cloud microphysics using the two best-performing simulations reveals that the event is strongly convective, and it is essential to keep the cumulus scheme on for such convective events and capture all the precipitation characteristics showing that in regions of extreme topography, the cumulus scheme is still necessary even down to the grid spacing of at least 3 km. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>Total accumulated precipitation during 11–16 August 2017 (<b>a</b>) recorded over 95 surface weather stations across Nepal and (<b>b</b>) measured by the Global Precipitation Measurement (GPM) satellite.</p>
Full article ">Figure 2
<p>WRF domains: (<b>a</b>) outer domain (D1, 15 km) and (<b>b</b>) inner domain (D2, 3 km). Topography is shaded in both domains. B1, B2, and B3 in (<b>b</b>) are the regions of intense precipitation where the hourly precipitation variability and cloud microphysics analyses are performed.</p>
Full article ">Figure 3
<p>Spatial distribution of accumulated precipitation over 11–16 August 2017. (<b>a</b>–<b>f</b>) WRF simulated precipitation accumulation without the cumulus scheme using (<b>a</b>) Lin, (<b>b</b>) Milbrandt, (<b>c</b>) Morrison, (<b>d</b>) Thompson, (<b>e</b>) WDM6, (<b>f</b>) WSM6 microphysics schemes; (<b>g</b>) represents station observations and (<b>h</b>) represents GPM data.</p>
Full article ">Figure 4
<p>Same as <a href="#atmosphere-16-00298-f003" class="html-fig">Figure 3</a>, but with the cumulus scheme turned on.</p>
Full article ">Figure 5
<p>Modelled versus observed precipitation, along with respective values of the compared statistical parameters.</p>
Full article ">Figure 6
<p>Comparison of 3-hourly accumulated precipitation (mm) between GPM satellite observations and WRF simulations using various microphysics schemes from 10 to 14 August 2017 over (<b>a</b>) B1, (<b>b</b>) B2, and (<b>c</b>) B3.</p>
Full article ">Figure 7
<p>Time series of spatial averaged column density hydrometeors from Morrison_cu_sim and WSM6_cu_sim. The solid line represents the region B1, the broken line represents B2, and the dotted line represents B3.</p>
Full article ">Figure 8
<p>Vertical profile of spatially averaged mixing ratio of hydrometeors. The solid line represents the region B1, the broken line represents B2, and the dotted line represents B3.</p>
Full article ">Figure 9
<p>Vertically averaged vertical wind speed (w, m/s) at 00:00, 12 August 2017 (local time) plotted from (<b>a</b>) Morrison_cu_sim and (<b>b</b>) Morrison_sim.</p>
Full article ">Figure 10
<p>Hourly grid scale, cumulus, and total accumulated precipitation simulated by Morrison_cu_sim averaged over (<b>a</b>) B1, (<b>b</b>) B2, and (<b>c</b>) B3.</p>
Full article ">Figure 11
<p>(<b>a</b>) Mean Sea level pressure anomaly (shaded) and 850 hPa wind anomalies; (<b>b</b>) 500 hPa geopotential height (shaded) and wind anomalies, and (<b>c</b>) 250 hPa geopotential height and wind anomalies averaged over 11–14 August 2017.</p>
Full article ">Figure 12
<p>Skew-T Log-P diagram from the Patna station, India for (<b>a</b>) 11 August 2017 and (<b>b</b>) 12 August 2017.</p>
Full article ">
27 pages, 12001 KiB  
Article
Numerical Simulation of Convective Heat Transfer in Gyroid, Diamond, and Primitive Microstructures Using Water as the Working Fluid
by Jie Zhang and Xiaoqing Yang
Energies 2025, 18(5), 1230; https://doi.org/10.3390/en18051230 - 3 Mar 2025
Viewed by 210
Abstract
With the continuous increase in the thermal power of electronic devices, air cooling is becoming increasingly challenging in terms of meeting heat dissipation requirements. Liquid cooling media have a higher specific heat capacity and better heat dissipation effect, making it a more efficient [...] Read more.
With the continuous increase in the thermal power of electronic devices, air cooling is becoming increasingly challenging in terms of meeting heat dissipation requirements. Liquid cooling media have a higher specific heat capacity and better heat dissipation effect, making it a more efficient cooling method. In order to improve the heat dissipation effect of liquid cooling, a TPMS structure with a larger specific surface area, which implicit function parameters can control, can be arranged in a shape manner and it is easy to expand the structural design. It has excellent potential for application in the field of heat dissipation. At present, research is still in its initial stage and lacks comparative studies on liquid cooled convective heat transfer of TPMS structures G (Gyroid), D (Diamond), and P (Primitive). This paper investigates the heat transfer performance and pressure drop characteristics of a sheet-like microstructure composed of classic TPMS structures, G (Gyroid), D (Diamond), and P (Primitive), with a single crystal cell length of 2π (mm), a cell number of 1 × 1 × 5, and a microstructure size of 2π (mm) × 2π (mm) × 22π (mm) using a constant temperature surface model. By analyzing the outlet temperature tout, structural pressure p, average convective heat transfer coefficient h0, Nusselt number Nu, and average wall friction factor f of the microstructure within the speed range of 0.01–0.11 m/s and constant temperature surface temperature is 100 °C, the heat transfer capacity D > G > P and pressure drop D > G > P were obtained (the difference in pressure drop between G and P is very small, less than 20 Pa, which can be considered consistent). When flow velocity is 0.01 m/s, the maximum temperature difference at the outlet of the four structures reached 17.14 °C, and the maximum difference in wall friction factor f reached 103.264, with a relative change of 646%. When flow velocity is 0.11 m/s, the maximum pressure difference among the four structures reached 8461.84 Pa, and the maximum difference in h0 reached 7513 W/(m2·K), with a relative change of 63.36%; the maximum difference between Nu reached 76.32, with a relative change of 62.09%. This paper explains the reasons for the above conclusions by analyzing the proportion of solid area on the constant temperature surface of the structure, the porosity of the structure, and the characteristics of streamlines in the microstructure. Full article
(This article belongs to the Section J1: Heat and Mass Transfer)
Show Figures

Figure 1

Figure 1
<p>TPMS-G structure with cycles of (2, 2, 7).</p>
Full article ">Figure 2
<p>TPMS to be repaired in SpaceClaim.</p>
Full article ">Figure 3
<p>TPMS initially presented in COMSOL.</p>
Full article ">Figure 4
<p>Imported TPMS Defects.</p>
Full article ">Figure 5
<p>Embedded rectangular prism in TPMS.</p>
Full article ">Figure 6
<p>TPMS-G structure with Boolean post cycle numbers of (1, 1, 5).</p>
Full article ">Figure 7
<p>Physical Model of TPMS-G Microstructure.</p>
Full article ">Figure 8
<p>Single cells of four TPMS structures.</p>
Full article ">Figure 9
<p>Four types of TPMS imported into COMSOL.</p>
Full article ">Figure 10
<p>Microstructure physical models of four types of TPMS.</p>
Full article ">Figure 11
<p>Mesh division of different refinement levels on the entrance surface of the microstructure TPMS.</p>
Full article ">Figure 12
<p>Outlet temperature under different grids.</p>
Full article ">Figure 13
<p>Structural pressure under different grids.</p>
Full article ">Figure 14
<p>Relative outlet temperature under different grids.</p>
Full article ">Figure 15
<p>Relative pressure under different grids.</p>
Full article ">Figure 16
<p>Flow velocity in the yz section of four TPMS microstructures.</p>
Full article ">Figure 17
<p>Velocity cloud map through the yz plane of microstructure geometry center.</p>
Full article ">Figure 18
<p>Velocity cloud map of microstructure TPMS inlet surface.</p>
Full article ">Figure 19
<p>Flow velocity in the yz section of TPMS in four microstructures.</p>
Full article ">Figure 20
<p>Structure of TPMS along the mainstream direction.</p>
Full article ">Figure 21
<p>Partial streamline of TPMS.</p>
Full article ">Figure 22
<p>Partial streamline of TPMS and outlet section of microstructure.</p>
Full article ">Figure 23
<p>Streamline of TPMS-G, D1, P, D2 outlet.</p>
Full article ">Figure 24
<p>Velocity cloud map through microstructure geometric center xy plane.</p>
Full article ">Figure 25
<p>Temperature cloud map of xy plane through microstructure geometry center.</p>
Full article ">Figure 26
<p>Temperature cloud map through the yz plane of microstructure geometry center.</p>
Full article ">Figure 27
<p>Outlet temperatures of four microstructures at different flow velocities.</p>
Full article ">Figure 28
<p>Four microstructure pressures at different flow velocities.</p>
Full article ">Figure 29
<p>Average convective heat transfer coefficients <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>h</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math> of four microstructures at different flow velocities.</p>
Full article ">Figure 30
<p>Nusselt numbers <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>u</mi> </mrow> </semantics></math> of four microstructures at different flow velocities.</p>
Full article ">Figure 31
<p>Average wall friction factors <math display="inline"><semantics> <mrow> <mi>f</mi> </mrow> </semantics></math> of four microstructures at different flow velocities.</p>
Full article ">
18 pages, 5664 KiB  
Article
Magnetohydrodynamic Blood-Carbon Nanotube Flow and Heat Transfer Control via Carbon Nanotube Geometry and Nanofluid Properties for Hyperthermia Treatment
by Nickolas D. Polychronopoulos, Evangelos Karvelas, Lefteris Benos, Thanasis D. Papathanasiou and Ioannis Sarris
Computation 2025, 13(3), 62; https://doi.org/10.3390/computation13030062 - 3 Mar 2025
Viewed by 231
Abstract
Hyperthermia is a promising medical treatment that uses controlled heat to target and destroy cancer cells while minimizing damage to the surrounding healthy tissue. Unlike conventional methods, it offers reduced risks of infection and shorter recovery periods. This study focuses on the integration [...] Read more.
Hyperthermia is a promising medical treatment that uses controlled heat to target and destroy cancer cells while minimizing damage to the surrounding healthy tissue. Unlike conventional methods, it offers reduced risks of infection and shorter recovery periods. This study focuses on the integration of carbon nanotubes (CNTs) within the blood to enable precise heat transfer to tumors. The central idea is that by adjusting the concentration, shape, and size of CNTs, as well as the strength of an external magnetic field, heat transfer can be controlled for targeted treatment. A theoretical model is developed to analyze laminar natural convection within a simplified rectangular porous enclosure resembling a tumor, considering the composition of blood, and the geometric characteristics of CNTs, including the interfacial nanolayer thickness. Using an asymptotic expansion method, ordinary differential equations for mass, momentum, and energy balances are derived and solved. Results show that increasing CNT concentration decelerates fluid flow and reduces heat transfer efficiency, while elongated CNTs and thicker nanolayers enhance conduction over convection, to the detriment of heat transfer. Finally, increased tissue permeability—characteristic of cancerous tumors—significantly impacts heat transfer. In conclusion, although the model simplifies real tumor geometries and treatment conditions, it provides valuable theoretical insights into hyperthermia and nanofluid applications for cancer therapy. Full article
(This article belongs to the Special Issue Post-Modern Computational Fluid Dynamics)
Show Figures

Figure 1

Figure 1
<p>Flow configuration and boundary conditions. Vertical and horizontal walls are electrically insulated (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>V</mi> </mrow> <mrow> <mi>E</mi> </mrow> </msub> </mrow> </semantics></math> is the electrical potential). Vertical walls are isothermal and horizontal walls are adiabatic.</p>
Full article ">Figure 2
<p>Effect of the different CNT geometrical factors and volume fraction on <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>k</mi> </mrow> <mrow> <mi>b</mi> <mi>i</mi> <mi>o</mi> </mrow> </msub> </mrow> </semantics></math>: (<b>a</b>) diameter, <math display="inline"><semantics> <mrow> <mi>d</mi> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>γ</mi> </mrow> </semantics></math> parameter, (<b>c</b>) length, <math display="inline"><semantics> <mrow> <mi>S</mi> </mrow> </semantics></math>, (<b>d</b>) nanolayer thickness, <math display="inline"><semantics> <mrow> <mi>t</mi> </mrow> </semantics></math>, and (<b>e</b>) volume fraction, <math display="inline"><semantics> <mrow> <mi>φ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 2 Cont.
<p>Effect of the different CNT geometrical factors and volume fraction on <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>k</mi> </mrow> <mrow> <mi>b</mi> <mi>i</mi> <mi>o</mi> </mrow> </msub> </mrow> </semantics></math>: (<b>a</b>) diameter, <math display="inline"><semantics> <mrow> <mi>d</mi> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>γ</mi> </mrow> </semantics></math> parameter, (<b>c</b>) length, <math display="inline"><semantics> <mrow> <mi>S</mi> </mrow> </semantics></math>, (<b>d</b>) nanolayer thickness, <math display="inline"><semantics> <mrow> <mi>t</mi> </mrow> </semantics></math>, and (<b>e</b>) volume fraction, <math display="inline"><semantics> <mrow> <mi>φ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Results for different CNT diameter values, <math display="inline"><semantics> <mrow> <mi>d</mi> </mrow> </semantics></math>: (<b>a</b>) dimensionless temperature profile, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) streamfunction, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ψ</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) dimensionless vertical velocity,<math display="inline"><semantics> <mrow> <mo> </mo> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>ο</mi> </mrow> </msub> </mrow> </semantics></math>, and (<b>d</b>) average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3 Cont.
<p>Results for different CNT diameter values, <math display="inline"><semantics> <mrow> <mi>d</mi> </mrow> </semantics></math>: (<b>a</b>) dimensionless temperature profile, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) streamfunction, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ψ</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) dimensionless vertical velocity,<math display="inline"><semantics> <mrow> <mo> </mo> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>ο</mi> </mrow> </msub> </mrow> </semantics></math>, and (<b>d</b>) average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Results for different <math display="inline"><semantics> <mrow> <mi>γ</mi> </mrow> </semantics></math> values in the core region: (<b>a</b>) dimensionless temperature profile, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) streamfunction, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ψ</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) dimensionless vertical velocity,<math display="inline"><semantics> <mrow> <mo> </mo> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>ο</mi> </mrow> </msub> </mrow> </semantics></math>, and (<b>d</b>) average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Results for different CNT length (<math display="inline"><semantics> <mrow> <mi>S</mi> </mrow> </semantics></math>) in the core region: (<b>a</b>) dimensionless temperature profile, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) streamfunction, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ψ</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) dimensionless vertical velocity,<math display="inline"><semantics> <mrow> <mo> </mo> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>ο</mi> </mrow> </msub> </mrow> </semantics></math>, and (<b>d</b>) average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Results for different CNT interfacial nanolayer thickness (<math display="inline"><semantics> <mrow> <mi>t</mi> </mrow> </semantics></math>) values in the core region: (<b>a</b>) dimensionless temperature profile, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) streamfunction, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ψ</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) dimensionless vertical velocity,<math display="inline"><semantics> <mrow> <mo> </mo> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>ο</mi> </mrow> </msub> </mrow> </semantics></math>, and (<b>d</b>) average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Results for different CNT volume fraction <math display="inline"><semantics> <mrow> <mi>φ</mi> </mrow> </semantics></math> values in the core region: (<b>a</b>) dimensionless temperature profile, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>b</b>) streamfunction, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ψ</mi> </mrow> <mrow> <mn>0</mn> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>, (<b>c</b>) dimensionless vertical velocity,<math display="inline"><semantics> <mrow> <mo> </mo> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>ο</mi> </mrow> </msub> </mrow> </semantics></math>, and (<b>d</b>) average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>Dependence of average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> </mrow> </msub> <mo> </mo> </mrow> </semantics></math> on the Hartmann number <math display="inline"><semantics> <mrow> <mi>H</mi> <mi>a</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Dependence of average Nusselt number <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>N</mi> <mi>u</mi> </mrow> <mrow> <mi>a</mi> <mi>v</mi> <mo>.</mo> </mrow> </msub> <mo> </mo> </mrow> </semantics></math> on the Darcy number <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>a</mi> </mrow> </semantics></math> for different values of the <math display="inline"><semantics> <mrow> <mi>γ</mi> </mrow> </semantics></math> parameter.</p>
Full article ">
23 pages, 17900 KiB  
Article
Unveiling the Impact of Microfractures on Longitudinal Dispersion Coefficients in Porous Media
by Muyuan Wang, Keliu Wu, Qingyuan Zhu and Jiawei Ye
Processes 2025, 13(3), 722; https://doi.org/10.3390/pr13030722 - 2 Mar 2025
Viewed by 295
Abstract
Longitudinal dispersion coefficient is a key parameter governing solute transport in porous media, with significant implications for various industrial processes. However, the impact of microfractures on the longitudinal dispersion coefficient remains insufficiently understood. In this study, pore-scale direct numerical simulations are performed to [...] Read more.
Longitudinal dispersion coefficient is a key parameter governing solute transport in porous media, with significant implications for various industrial processes. However, the impact of microfractures on the longitudinal dispersion coefficient remains insufficiently understood. In this study, pore-scale direct numerical simulations are performed to analyze solute transport in microfractured porous media during unstable miscible displacement. Spatiotemporal concentration profiles were fitted to the analytical solution of the convection–dispersion equation to quantify the longitudinal dispersion coefficient across different microfracture configurations. The results indicate that the longitudinal dispersion coefficient is highly sensitive to microfracture characteristics. Specifically, an increased projection length of microfractures in the flow direction and a reduced lateral projection length enhance longitudinal dispersion at the outlet. When Peclet number ≥1, the longitudinal dispersion coefficient follows a three-stage variation pattern along the flow direction, with microfracture connectivity and orientation dominating its scale sensitivity. Furthermore, both diffusion-dominated and mixed advective-diffusion regimes are observed. In diffusion-dominated regimes, significant channeling alters the applicability of traditional scaling laws, with the relationship between longitudinal dispersion coefficient and porosity holding only when the Peclet number is below 0.07. These results provide a comprehensive scale-up framework for CO2 miscible flooding in unconventional reservoirs and CO2 storage in saline aquifers, offering valuable insights for the numerical modeling of heterogeneous reservoir development. Full article
(This article belongs to the Section Energy Systems)
Show Figures

Figure 1

Figure 1
<p>Granular porous medium before microfracture incorporation. (Note: The computational domain (<span class="html-italic">L</span> × <span class="html-italic">W</span>) dimensions are 1000 μm × 50 μm; the white regions represent solid particles and the blue regions represent the pore space).</p>
Full article ">Figure 2
<p>SSD and PSD in porous medium.</p>
Full article ">Figure 3
<p>Porous media domains with embedded microfractures. Microfractures in M-PM1 to M-PM4 are parallel to the flow direction but vary in length and density. (<b>a</b>) M-PM1. (<b>b</b>) M-PM2. (<b>c</b>) M-PM3. (<b>d</b>) M-PM4. (<b>e</b>) M-PM5. (<b>f</b>) M-PM6. (<b>g</b>) M-PM7. (<b>h</b>) M-PM8. (Note: In M-PM1 to M-PM4, microfractures are aligned parallel to the flow direction but differ in length and density. M-PM5 and M-PM6 contain microfractures oriented at 90° and 45° to the flow direction, respectively, while M-PM8 exhibits enhanced microfracture connectivity compared to M-PM7).</p>
Full article ">Figure 4
<p>Outlet concentration curves for mesh resolutions during CO<sub>2</sub> miscible displacements in M-PM5 at <span class="html-italic">Pe</span> = 100 (<b>a</b>) and local mesh details in Mesh 3 (<b>b</b>).</p>
Full article ">Figure 5
<p>Comparison of simulated and analytical outlet concentration data in M-PM1 at <span class="html-italic">Pe</span> = 0.01.</p>
Full article ">Figure 6
<p>Comparison of numerical simulation and analytical solutions (Berkowitz and Zhou [<a href="#B48-processes-13-00722" class="html-bibr">48</a>]; Wang et al. [<a href="#B49-processes-13-00722" class="html-bibr">49</a>]) for <span class="html-italic">D<sub>L</sub></span> and <span class="html-italic">Pe</span> in stable miscible flow between two parallel flat plates.</p>
Full article ">Figure 7
<p>Variation in dimensionless <span class="html-italic">D<sub>L</sub></span> (<span class="html-italic">D<sub>L</sub></span>/<span class="html-italic">D</span>) along flow direction from M-PM1 to M-PM8 under different <span class="html-italic">Pe</span>. (<b>a</b>) <span class="html-italic">Pe</span> = 0.01. (<b>b</b>) <span class="html-italic">Pe</span> = 0.1. (<b>c</b>) <span class="html-italic">Pe</span> = 1. (<b>d</b>) <span class="html-italic">Pe</span> = 10. (<b>e</b>) <span class="html-italic">Pe</span> = 50. (<b>f</b>) <span class="html-italic">Pe</span> = 100.</p>
Full article ">Figure 7 Cont.
<p>Variation in dimensionless <span class="html-italic">D<sub>L</sub></span> (<span class="html-italic">D<sub>L</sub></span>/<span class="html-italic">D</span>) along flow direction from M-PM1 to M-PM8 under different <span class="html-italic">Pe</span>. (<b>a</b>) <span class="html-italic">Pe</span> = 0.01. (<b>b</b>) <span class="html-italic">Pe</span> = 0.1. (<b>c</b>) <span class="html-italic">Pe</span> = 1. (<b>d</b>) <span class="html-italic">Pe</span> = 10. (<b>e</b>) <span class="html-italic">Pe</span> = 50. (<b>f</b>) <span class="html-italic">Pe</span> = 100.</p>
Full article ">Figure 8
<p>Impact of strong flow direction heterogeneity on CO<sub>2</sub> concentration front propagation. (<b>a</b>) M-PM1. (<b>b</b>) M-PM2. (<b>c</b>) M-PM8. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM1, M-PM2, and M-PM3 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 0.1).</p>
Full article ">Figure 9
<p>Impact of single horizontal microfracture length on CO<sub>2</sub> concentration front propagation. (<b>a</b>) M-PM1, <span class="html-italic">Pe</span> = 1. (<b>b</b>) M-PM2, <span class="html-italic">Pe</span> = 1. (<b>c</b>) M-PM1, Pe = 10. (<b>d</b>) M-PM2, <span class="html-italic">Pe</span> = 10. (<b>e</b>) M-PM1, <span class="html-italic">Pe</span> = 100. (<b>f</b>) M-PM2, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM1 and M-PM2 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 1, 10, 100).</p>
Full article ">Figure 10
<p>Impact of microfracture parallel density on CO<sub>2</sub> concentration front propagation. (<b>a</b>) M-PM3, <span class="html-italic">Pe</span> = 1. (<b>b</b>) M-PM4, <span class="html-italic">Pe</span> = 1. (<b>c</b>) M-PM3, <span class="html-italic">Pe</span> = 10. (<b>d</b>) M-PM4, <span class="html-italic">Pe</span> = 10. (<b>e</b>) M-PM3, <span class="html-italic">Pe</span> = 100. (<b>f</b>) M-PM4, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM3 and M-PM4 at <span class="html-italic">t<sub>D</sub></span> = 0.8 and <span class="html-italic">Pe</span> = 1, 10, 100).</p>
Full article ">Figure 11
<p>Dimensionless CO<sub>2</sub> concentration profiles at <span class="html-italic">x<sub>D</sub></span> = 1 for M-PM1 to M-PM4 as a function of <span class="html-italic">t<sub>D</sub></span> at <span class="html-italic">Pe</span> = 1.</p>
Full article ">Figure 12
<p>Impact of microfracture orientation on CO<sub>2</sub> concentration front propagation at high <span class="html-italic">Pe</span>. (<b>a</b>) M-PM3, <span class="html-italic">Pe</span> = 50. (<b>b</b>) M-PM3, <span class="html-italic">Pe</span> = 100. (<b>c</b>) M-PM5, <span class="html-italic">Pe</span> = 50. (<b>d</b>) M-PM5, <span class="html-italic">Pe</span> = 100. (<b>e</b>) M-PM6, <span class="html-italic">Pe</span> = 50. (<b>f</b>) M-PM6, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM3, M-PM5, and M-PM6 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 50, 100).</p>
Full article ">Figure 13
<p>Impact of microfracture connectivity on CO<sub>2</sub> concentration front propagation at high <span class="html-italic">Pe</span>. (<b>a</b>) M-PM7, <span class="html-italic">Pe</span> = 100. (<b>b</b>) M-PM7, <span class="html-italic">Pe</span> = 50. (<b>c</b>) M-PM8, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> for M-PM7 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 50, 100; <span class="html-italic">c<sub>D</sub></span> distribution for M-PM8 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 100).</p>
Full article ">Figure 14
<p>Perturbation of CO<sub>2</sub> concentration front in Stage 1 (<b>a</b>–<b>d</b>), Stage 2 (<b>e</b>), and Stage 3 (<b>f</b>).</p>
Full article ">Figure 14 Cont.
<p>Perturbation of CO<sub>2</sub> concentration front in Stage 1 (<b>a</b>–<b>d</b>), Stage 2 (<b>e</b>), and Stage 3 (<b>f</b>).</p>
Full article ">Figure 15
<p>Impact of microfracture orientation and connectivity on CO<sub>2</sub> concentration distribution along domain length. (<b>a</b>) M-PM1. (<b>b</b>) M-PM6. (<b>c</b>) M-PM7. (<b>d</b>) M-PM8.</p>
Full article ">Figure 16
<p>Identification of dispersion regimes in unstable displacement process.</p>
Full article ">
19 pages, 4329 KiB  
Article
Experimental Analysis of Heat and Flow Characteristics on Inclined and Multiple Impingement Jet Heat Transfer Using Optimized Heat Sink
by Altug Karabey and Dogan Yorulmaz
Appl. Sci. 2025, 15(5), 2657; https://doi.org/10.3390/app15052657 - 1 Mar 2025
Viewed by 312
Abstract
Thermal management at a high heat flux is crucial for electronic devices, and jet impingement cooling is a promising solution. The heat transfer properties of a rectangular-finned heat sink are investigated under angled and multi-impingement jet configurations in this study. Experiments were conducted [...] Read more.
Thermal management at a high heat flux is crucial for electronic devices, and jet impingement cooling is a promising solution. The heat transfer properties of a rectangular-finned heat sink are investigated under angled and multi-impingement jet configurations in this study. Experiments were conducted with three different nozzle diameters, three different heat sink angles, three dimensionless nozzle-to-heat sink distance ratios, and five different velocity values. As a result, the obtained data are presented as Nu-Re graphs, and the impacts of the parameters on heat transfer (HT) are analyzed. It is concluded that the Nusselt number increases with the increasing nozzle diameter and Reynolds number, whereas it decreases with increasing distance between the nozzle and the heat sink. When comparing the angle values under an identical flow velocity, nozzle diameter, and dimensionless h/d distance experimental conditions, it was found that the Nusselt numbers were very close to each other. Under constant heat flux and for all investigated angles, the highest Nusselt number for the rectangular-finned inclined heat sink was observed at a 10° heat sink inclination, a nozzle diameter of D = 40 mm, a dimensionless distance of h/d = 6, and a flow velocity of 9 m/s. This study deepens the understanding of the heat transfer mechanism of impinging jets and provides an efficient method framework for practical applications. Full article
Show Figures

Figure 1

Figure 1
<p>Test setup. (1) Velocity regulator, (2) fan, (3) valve, (4) regulator, (5) variac, (6) mechanism for the lifting jack and heat sink angle, (7) impingement multi-jet nozzles, (8) angularly adjustable finned heat sink, (9) thermocouples, (10) datalogger, and (11) computer.</p>
Full article ">Figure 2
<p>The general characteristics of the heat sink.</p>
Full article ">Figure 3
<p>Comparison of experimental results for inclined heat sink (α = 10°) at constant heat flux (200 W) and h/d = 6.</p>
Full article ">Figure 4
<p>Comparison of experimental results for inclined heat sink (α = 10°) at constant heat flux (200 W) and h/d = 7.</p>
Full article ">Figure 5
<p>Comparison of experimental results for inclined heat sink (α = 10°) at constant heat flux (200 W) and h/d = 8.</p>
Full article ">Figure 6
<p>Comparison of experimental results for inclined heat sink (α = 20°) at constant heat flux (200 W) and h/d = 6.</p>
Full article ">Figure 7
<p>Comparison of experimental results for inclined heat sink (α = 20°) at constant heat flux (200 W) and h/d = 7.</p>
Full article ">Figure 8
<p>Comparison of experimental results for inclined heat sink (α = 20°) at constant heat flux (200 W) and h/d = 8.</p>
Full article ">Figure 9
<p>Comparison of experimental results for inclined heat sink (α = 30°) at constant heat flux (200 W) and h/d = 6.</p>
Full article ">Figure 10
<p>Comparison of experimental results for inclined heat sink (α = 30°) at constant heat flux (200 W) and h/d = 7.</p>
Full article ">Figure 11
<p>Comparison of experimental results for inclined heat sink (α = 30°) at constant heat flux (200 W) and h/d = 8.</p>
Full article ">
15 pages, 4352 KiB  
Article
Unraveling Mass Transfer and Reaction Processes in CVD-Grown MoS2 Films: A Multiphysical Field Coupling Study
by Zhen Yang, Jinwei Lin, Qing Zhang, Yutian Liu, Shujun Han, Yanbin Zhou, Shuo Chen, Shenlong Zhong, Xianli Su, Qingjie Zhang and Xinfeng Tang
Appl. Sci. 2025, 15(5), 2627; https://doi.org/10.3390/app15052627 - 28 Feb 2025
Viewed by 233
Abstract
The two-dimensional semiconductor material MoS2, grown via chemical vapor deposition, has shown significant potential to surpass silicon in advanced electronic technologies. However, the mass transfer and chemical reaction processes critical to the nucleation and growth of MoS2 grains remain poorly [...] Read more.
The two-dimensional semiconductor material MoS2, grown via chemical vapor deposition, has shown significant potential to surpass silicon in advanced electronic technologies. However, the mass transfer and chemical reaction processes critical to the nucleation and growth of MoS2 grains remain poorly understood. In this study, we conducted an in-depth investigation into the mass transfer and chemical reaction processes during the chemical vapor deposition of MoS2, employing a novel multi-physics coupling model that integrates flow fields, temperature fields, mass transfer, and chemical reactions. Our findings reveal that the intermediate product Mo3O9S4 not only fails to participate directly in MoS2 film growth but also hinders the diffusion of MoS6, limiting the growth process. We demonstrate that increasing the growth temperature accelerates the diffusion rate of MoS6, mitigates the adverse effects of Mo3O9S4, and promotes the layered growth of MoS2 films. Additionally, lowering the growth pressure enhances the convective diffusion of reactants, accelerating grain growth. This research significantly advances our understanding of the mass transport and reaction processes in MoS2 film growth and provides critical insights for optimizing chemical vapor deposition systems. Full article
Show Figures

Figure 1

Figure 1
<p>Finite element geometric model of dual temperature zone CVD.</p>
Full article ">Figure 2
<p>Optical micrographs of the sample surface grown at (<b>a</b>) 1033 K, (<b>b</b>) 1053 K, and (<b>c</b>) 1073 K under a constant pressure of 70 Torr, along with the corresponding (<b>d</b>) Raman spectra and (<b>e</b>) PL spectra.</p>
Full article ">Figure 3
<p>Radial cross-section of Mo<sub>3</sub>O<sub>9</sub> concentration distribution at the center of the substrate when the pressure is 70 Torr and the temperature in growth region is set to (<b>a</b>,<b>b</b>) 1033 K, (<b>c</b>,<b>d</b>) 1053 K, and (<b>e</b>,<b>f</b>) 1073 K.</p>
Full article ">Figure 4
<p>Radial cross-section of the concentration distribution of Mo<sub>3</sub>O<sub>9</sub>S<sub>4</sub> at the center of the substrate when the temperature in growth region is set to (<b>a</b>) 1033 K, (<b>b</b>) 1053 K, and (<b>c</b>) 1073 K at a pressure of 70 Torr.</p>
Full article ">Figure 5
<p>Concentration distribution of MoS<sub>6</sub> when the temperature in growth region is set to (<b>a</b>–<b>c</b>) 1033 K, (<b>d</b>–<b>f</b>) 1053 K, and (<b>g</b>–<b>i</b>) 1073 K at a pressure of 70 Torr, where (<b>b</b>,<b>e</b>,<b>h</b>) are the radial cross-sectional views at the center of the substrate and (<b>c</b>,<b>f</b>,<b>i</b>) are the axial cross-sectional views at the center of the substrate.</p>
Full article ">Figure 6
<p>Optical microscopic images of the sample grown at a growth temperature of 1053 K under (<b>a</b>) 10 Torr, (<b>b</b>) 40 Torr, and (<b>c</b>) 100 Torr, along with the corresponding (<b>d</b>) Raman spectra and (<b>e</b>) PL spectra.</p>
Full article ">Figure 7
<p>Velocity distribution diagrams of the system at a pressure of (<b>a</b>–<b>c</b>) 10 Torr, (<b>d</b>–<b>f</b>) 40 Torr, and (<b>g</b>–<b>i</b>) 100 Torr when the temperature is set to 1053 K in growth region, where (<b>b</b>,<b>e</b>,<b>h</b>) are radial cross-section diagrams at the center of the substrate and (<b>c</b>,<b>f</b>,<b>i</b>) are corresponding local magnified diagrams.</p>
Full article ">Figure 8
<p>Radial cross-sectional view of the Mo<sub>3</sub>O<sub>9</sub> concentration distribution at the center of the substrate at a set temperature of 1053 K in growth region, with a pressure of (<b>a</b>,<b>b</b>) 10 Torr, (<b>c</b>,<b>d</b>) 40 Torr, and (<b>e</b>,<b>f</b>) 100 Torr.</p>
Full article ">Figure 9
<p>Radial cross-sectional view of Mo<sub>3</sub>O<sub>9</sub>S<sub>4</sub> concentration distribution at the center of the substrate in growth region, with a temperature set at 1053 K and pressures of (<b>a</b>) 10 Torr, (<b>b</b>) 40 Torr, and (<b>c</b>) 100 Torr.</p>
Full article ">Figure 10
<p>Concentration distribution of MoS<sub>6</sub> when the temperature is set to 1053 K in a growth region, with the gas pressure at (<b>a</b>–<b>c</b>) 10 Torr, (<b>d</b>–<b>f</b>) 40 Torr, and (<b>g</b>–<b>i</b>) 100 Torr, where (<b>b</b>,<b>e</b>,<b>h</b>) are the radial cross-sectional views at the center of the substrate and (<b>c</b>,<b>f</b>,<b>i</b>) are the axial cross-sectional views at the center of the substrate.</p>
Full article ">
17 pages, 4259 KiB  
Article
Analyzing an Extreme Rainfall Event in Himachal Pradesh, India, to Contribute to Sustainable Development
by Nitin Lohan, Sushil Kumar, Vivek Singh, Raj Pritam Gupta and Gaurav Tiwari
Sustainability 2025, 17(5), 2115; https://doi.org/10.3390/su17052115 - 28 Feb 2025
Viewed by 307
Abstract
In the Himalayan regions of complex terrains, such as Himachal Pradesh, the occurrence of extreme rainfall events (EREs) has been increasing, triggering landslides and flash floods. Investigating the dynamics and precipitation characteristics and improving the prediction of such events are crucial and could [...] Read more.
In the Himalayan regions of complex terrains, such as Himachal Pradesh, the occurrence of extreme rainfall events (EREs) has been increasing, triggering landslides and flash floods. Investigating the dynamics and precipitation characteristics and improving the prediction of such events are crucial and could play a vital role in contributing to sustainable development in the region. This study employs a high-resolution numerical weather prediction framework, the weather research and forecasting (WRF) model, to deeply investigate an ERE which occurred between 8 July and 13 July 2023. This ERE caused catastrophic floods in the Mandi and Kullu districts of Himachal Pradesh. The WRF model was configured with nested domains of 12 km and 4 km horizontal grid resolutions, and the results were compared with global high-resolution precipitation products and the fifth-generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis dataset. The selected case study was amplified by the synoptic scale features associated with the position and intensity of the monsoon trough, including mesoscale processes like orographic lifting. The presence of a western disturbance and the heavy moisture transported from the Arabian Sea and the Bay of Bengal both intensified this event. The model has effectively captured the spatial distribution and large-scale dynamics of the phenomenon, demonstrating the importance of high-resolution numerical modeling in accurately simulating localized EREs. Statistical evaluation revealed that the WRF model overestimated extreme rainfall intensity, with the root mean square error reaching 17.33 mm, particularly during the convective peak phase. The findings shed light on the value of high-resolution modeling in capturing localized EREs and offer suggestions for enhancing disaster management and flood forecasting. Full article
(This article belongs to the Section Air, Climate Change and Sustainability)
Show Figures

Figure 1

Figure 1
<p>The daily evolution of infrared brightness temperature (Unit: Kelvin) was derived from the INSAT-3DR satellite product. Panel figures (<b>a</b>–<b>f</b>) are plotted from 8 to 13 July 2023, respectively.</p>
Full article ">Figure 2
<p>The plotted area of the figure demonstrates the dimensions of the outer domain (D01). A rectangular box indicates the dimensions of the inner domain (D02), along with the topography of the study domains.</p>
Full article ">Figure 3
<p>(<b>a</b>) Climatological mean rainfall distribution (mm/day) for the six days (8 July to 13 July) over the 40 years (1984 to 2023); (<b>b</b>) Rainfall anomaly for the period from 8 July to 13 July for 2023.</p>
Full article ">Figure 4
<p>Spatiotemporal distribution of daily rainfall (mm) valid for six days (8 to 13 July 2023) from the IMD gridded data (<b>top row</b>), ERA5 (<b>second row</b>), MSWEP data (<b>third row</b>), and the WRF model’s inner domain simulation (<b>bottom row</b>).</p>
Full article ">Figure 5
<p>The Equitable Threat Score (ETS) for simulated rainfall (inner domain) validated against the MSWEP product at various threshold values from 8 July to 13 July 2023.</p>
Full article ">Figure 6
<p>Vertically integrated moisture transport (VIMT; kg.m<sup>−1</sup>.s<sup>−1</sup>) for all six days from the ERA5 data. The contours are presenting the VIMT and vectors denote the flow of moisture transport.</p>
Full article ">Figure 7
<p>Vertically integrated moisture transport (VIMT; kg.m<sup>−1</sup>.s<sup>−1</sup>) for all six days from the WRF model simulation.</p>
Full article ">Figure 8
<p>Area-averaged pressure vs. time vertical distribution of relative humidity (%) from (<b>a</b>) ERA5 and (<b>b</b>) WRF simulation for the inner domain.</p>
Full article ">Figure 9
<p>700 hPa daily geopotential height (m) and wind flow (m/s) from ERA5 (<b>first</b> and <b>second</b> rows) and WRF model’s outer domain simulation (<b>third</b> and <b>fourth</b> rows) valid for 8–13 July 2023.</p>
Full article ">Figure 10
<p>Extreme rainfall events disaster preparedness block diagram.</p>
Full article ">
Back to TopTop