Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DRAG”

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

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

At least 55 records · Page 3

Design of a Preliminary Family of Airfoils for High Reynolds Number Wind Turbine Applications

For the past 30 years, offshore wind turbines exhibited a continual pattern of growth that is expected to continue as the industry pushes for higher efficiency. Current designs for the next generation of wind turbines are so large that the chordwise Reynolds number of the blades is well beyond the design range of existing open-source airfoil families. This paper presents a preliminary family of new airfoils designed specifically for the needs of these next-generation offshore turbines, ranging from 21% thick to 30% thick with operating Reynolds numbers between 12 million and 18 million. These airfoils are intended to be alternative to the FFA airfoils that are commonly used on reference turbines such as the IEA 15MW and 22MW designs. In this work, airfoil performance metrics and design targets are developed, the design process is outlined, an optimization scheme is presented, and finally the airfoils and their simulated performance are compared to existing baselines. Lift to drag ratios in a clean condition were improved by up to 49.3% from the baseline FFA airfoil, and rough condition lift to drag ratio was improved by up to 9.3%. It is estimated that the cumulative improvements provided by this airfoil family would result in an approximately 1% increase of Annual Expected Power (AEP) for the 22 MW turbine compared to the current baseline.

airfoils↗

Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Astronomical simulations↗

OC7 phase I: CFD investigation of viscous forces on rectangular members of semisubmersibles

In Phase I of the Offshore Code Comparison Collaboration 7 (OC7) project, we aim to address key questions related to the modeling of viscous hydrodynamic loads on floating platforms. Here, this paper presents results from a collaborative effort within Phase I of the OC7 project. We investigate viscous forces acting on the rectangular members of semisubmersible platforms using high-fidelity computational fluid dynamics (CFD) simulations. The parametric CFD analyses include two-and-a-half-dimensional (2.5D) simulations of pontoon sections of two different semisubmersible platform designs (the VolturnUS-S and the INO WINDMOOR) and full three-dimensional (3D) structure simulations of the VolturnUS-S platform. We extract both the added mass and drag coefficients that are typically used in engineering-level hydrodynamic models from the 2.5D forced oscillation simulations with various Keulegan-Carpenter numbers and compare them to the coefficients from the 3D forced oscillation simulations. We observe consistent variation of the coefficients along the pontoon legs due to 3D platform and free-surface effects. With the breadth of the simulations and component-level detailed analyses, this work establishes a useful database and several empirical formulas for selecting hydrodynamic coefficients when modeling semisubmersible platforms in engineering tools.

13 HYDRO ENERGY↗

Identification and Exploration of a Series of SARS-Cov-2 M Pro Cyano-Based Inhibitors Revealing Ortho-Substitution Effects within the P3 Biphenyl Group

Starting from a simple scaffold hopping exercise based on our previous exploration of cysteine protease inhibitors against legumain, compound 6a was identified as a starting point for the development of a SARS-CoV-2 main protease (M Pro ) inhibitor. Compound 6a displayed submicromolar biochemical potency in the ultrasensitive assay developed by Drag and coworkers. Through an iterative structure−activity relationship campaign, we discovered an unexpected improvement in both biochemical and cellular potency through the incorporation of an ortho substituent within the P3 benzamide. X-ray crystallography revealed that incorporation of the ortho substituent caused a subtle but important binding enhancement of the P1 glutamate group within the M Pro S1 pocket. While incorporation of the ortho substituent improved the potency, the off-target selectivity against a panel of cysteine proteases and cell activity remained suboptimal. Further scanning of the P2 core revealed that incorporation of the 3.1.0 proline could address these issues and afford compound 22e, a highly potent and cellularly active M Pro inhibitor.

COVID-19 infection↗

Shape effects on the local dynamics of suspensions of spheroidal particles

The effect of shape on the dynamics of suspensions of non-spherical heavy particles is examined by fully resolved numerical simulations of oblate and prolate spheroids, as well as spheres, for a density ratio of ten, volume fractions ranging from 0.5% to 5%, and Reynolds numbers between 20 and 30. The dynamics is determined both by the interactions of the particles with the fluid as well as by collisions, with the number and importance of collisions increasing with volume fractions. A single isolated oblate or prolate spheroid falling under gravity generally falls broadside on, for the governing parameters examined here, and at low-volume fractions, the majority of particles in a suspension fall that way. At higher-volume fractions, the orientation is more random. The slip velocity decreases as the volume fraction increases for all shapes, as expected, but the effect of the shape is much less than seen for a single particle. Furthermore, this seems to be due to two effects. For all volume fractions, the most deformed particles cluster more than spheres and less deformed particles, which increases their slip velocity. As the concentration increases, the increased particle interactions also causes more particles to fall short side-on, which reduces the frontal area and the resulting drag, increasing the slip velocity. This second effect is, of course, absent for spherical particles.

42 ENGINEERING↗

TEAMER: Pitching Foil Crossflow Turbine Efficiency Data

This dataset documents the efficiency testing of a pitching foil crossflow turbine, conducted at the University of New Hampshire's (UNH) Chase Ocean Engineering Laboratory tow tank facility. The tests explored various pitch phases and amplitudes, ranging from 0 to 18 degrees, across different flow speeds and turbine RPMs. Specifically, testing was performed at pitch amplitudes of 0,3,6,9,12,15 and 18 degrees. The results indicate that turbine efficiency improves with increased pitch amplitude, reaching optimal performance at 12 degrees. Modifications were made to the test frame by UNH to enhance the stiffness of the support structure. The testing specifically mapped the efficiency of a single quadrant of the Ocean Renewable Power Company (ORPC) Autonomous Turbine Generator Unit (ATGU), covering the full range of allowable pitching amplitudes and phases. The dataset provides both raw and processed data, including drag and torque measurements, and contains Python scripts used for data processing and visualization. This research was supported by funding from TEAMER RFTS 1 (Request for Technical Support).

16 TIDAL AND WAVE POWER↗

Wireless Frequency‐Multiplexed Acoustic Array‐Based Acoustofluidics

Abstract Acoustofluidics has shown great potential in enabling on‐chip technologies for driving liquid flows and manipulating particles and cells for engineering, chemical, and biomedical applications. To introduce on‐demand liquid sample processing and micro/nano‐object manipulation functions to wearable and embeddable electronics, wireless acoustofluidic chips are highly desired. This paper presents wireless acoustofluidic chips to generate acoustic waves carrying sufficient energy and achieve key acoustofluidic functions, including arranging particles and cells, generating fluid streaming, and enriching in‐droplet particles. To enable these functions, the wireless acoustofluidic chips leverage mechanisms, including inductive coupling‐based wireless power transfer (WPT), frequency multiplexing‐based control of multiple acoustic waves, and the resultant acoustic radiation and drag forces. For validation, the wirelessly generated acoustic waves are measured using laser vibrometry when different materials (e.g., bone, tissue, and hand) are inserted between the WPT transmitter and receiver. Moreover, the wireless acoustofluidic chips successfully arrange nanoparticles into different patterns, align cells into parallel pearl chains, generate streaming, and enrich in‐droplet microparticles. This research is anticipated to facilitate the development of embeddable wireless on‐chip flow generators, wearable sensors with liquid sample processing functions, and implantable devices with flow generation and acoustic stimulation abilities for engineering, veterinary, and biomedical applications.

Li, Jiali↗

Designing Spin‐Crossover Systems to Enhance Thermopower and Thermoelectric Figure‐of‐Merit in Paramagnetic Materials

Thermoelectric materials, capable of converting temperature gradients into electrical power, have been traditionally limited by a trade‐off between thermopower and electrical conductivity. This study introduces a novel, broadly applicable approach that enhances both the spin‐driven thermopower and the thermoelectric figure‐of‐merit (zT) without compromising electrical conductivity, using temperature‐driven spin crossover. Our approach, supported by both theoretical and experimental evidence, is demonstrated through a case study of chromium doped‐manganese telluride, but is not confined to this material and can be extended to other magnetic materials. By introducing dopants to create a high crystal field and exploiting the entropy changes associated with temperature‐driven spin crossover, we achieved a significant increase in thermopower, by approximately 136 μV K −1 , representing more than a 200% enhancement at elevated temperatures within the paramagnetic domain. Our exploration of the bipolar semiconducting nature of these materials reveals that suppressing bipolar magnon/paramagnon‐drag thermopower is key to understanding and utilizing spin crossover‐driven thermopower. These findings, validated by inelastic neutron scattering, X‐ray photoemission spectroscopy, thermal transport, and energy conversion measurements, shed light on crucial material design parameters. We provide a comprehensive framework that analyzes the interplay between spin entropy, hopping transport, and magnon/paramagnon lifetimes, paving the way for the development of high‐performance spin‐driven thermoelectric materials.

magnons↗

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗

Leveraging interpolation models and error bounds for verifiable scientific machine learning

Effective verification and validation techniques for modern scientific machine learning workflows are challenging to devise. Statistical methods are abundant and easily deployed, but often rely on speculative assumptions about the data and methods involved. Error bounds for classical interpolation techniques can provide mathematically rigorous estimates of accuracy, but often are difficult or impractical to determine computationally. Here, in this work, we present a best-of-both-worlds approach to verifiable scientific machine learning by demonstrating that (1) multiple standard interpolation techniques have informative error bounds that can be computed or estimated efficiently; (2) comparative performance among distinct interpolants can aid in validation goals; (3) deploying interpolation methods on latent spaces generated by deep learning techniques enables some interpretability for black-box models. We present a detailed case study of our approach for predicting lift-drag ratios from airfoil images. Code developed for this work is available in a public Github repository.

97 MATHEMATICS AND COMPUTING↗

Urban wash-off of tire wear particles

Tire wear particles (TWPs) are an important class of microplastics due to their toxicity and abundance. Because most TWPs are generated on impervious road surfaces, urban wash-off is the critical first phase of waterborne transport from their zone of production to stormwater drainage. However, little is known about the driving factors behind their mobilization. In this study, we use a rainfall simulator to investigate how surface roughness, rainfall intensity, and surface slope affect wash-off behaviors of TWPs. We also analyze how the size and shape of mobilized TWPs change over the course of simulated storm events. We found that low surface roughness, high rainfall intensity (most significant factor), and low slope result in the most rapid conveyance of TWP load. On average, large particles (>1000 µm) travelled faster than small particles (<125 µm). Particle shape explained a very small amount of variance in TWP wash-off velocity but was found to be more important under higher surface roughness conditions. In addition to wash-off velocity, we found similar conditions controlled the percent mobilization of TWPs. Low surface roughness and high rainfall intensity resulting in higher TWP wash-off rates is consistent with mineral sediment wash-off behavior. Conversely, low surface slope and large particle size leading to faster conveyance is directly opposed to mineral sediment wash-off. Our findings suggest drag-dominated flow and that sufficient runoff depth is the most important parameter governing TWP wash-off. These findings are important first steps to understanding wash-off behaviors of TWPs and informing future modeling efforts and mitigation strategies.

13 HYDRO ENERGY↗

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Understanding the effect of minor alloying elements on helium bubble formation in ferritic-martensitic steels

Ferritic-martensitic steels are promising structural materials for advanced nuclear reactors. To minimize long-term radioactivity, reduced-activation ferritic-martensitic steels have been developed by substituting high-activation elements like Ni and Mo with low-activation elements such as W. However, the impact of these alloying modifications on helium bubble formation, which plays a key role in material swelling, remains unclear. Here, in this study, we compared helium bubble formation in ferritic-martensitic steel T91 and reduced-activation ferritic-martensitic steel F82H. Both materials were irradiated with sequential 100 keV, 150 keV, and 200 keV helium ions to a dose of 0.5 dpa and a helium concentration of 9,000 appm at 500°C. The helium bubbles in F82H exhibited a larger average size and a lower density than those in T91, suggesting differences in minor alloying elements may influence the bubble growth. Here, to investigate the effects of these alloying elements, we characterized radiation-induced segregation near bubbles and grain boundaries. Prominent Ni-Mn-Si enriched clusters were found near bubbles in T91, while only Mn-Si enriched clusters were found near bubbles in F82H. In addition, the obvious Cr enrichment near grain boundaries was absent around bubbles in both steels. The different segregation trends among elements revealed the variations in element diffusion mechanisms and the different sink biases between bubbles and grain boundaries. Cr enrichment near grain boundaries is mostly driven by interstitial-mediated diffusion. However, since bubble growth relies on net vacancy flux, vacancy-mediated diffusion plays a dominant role in controlling element segregation near bubbles. Therefore, Cr enrichment was not found near bubbles. Because of preferential vacancy-drag diffusion for Ni, Si and Mn, these elements were enriched near bubbles. Due to the strong binding energies of vacancies with these solute atoms, the vacancy diffusivity can be reduced near these solutes. Therefore, the more prominent Ni-Si-Mn clustered near helium bubbles in T91 lead to stronger suppression of helium bubble growth compared to F82H.

36 MATERIALS SCIENCE↗

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗