Richardson number Research Papers - Academia.edu (original) (raw)

2025

1 Département de Génie Mécanique, Université Mentouri Constantine, Route de Ain El Bey, 25000, Constantine, Algérie 2 Département Génie Mécanique, Université Oum El Bouaghi, Algérie 3 Département des Sciences et Techniques, Faculté des Sciences et de la Technologie, Centre Universitaire Abdelhafid Boussouf – Mila, Algérie 4 Département de Génie des Transport, Université Mentouri Constantine, Route de Ain El Bey, 25000, Constantine, Algérie 5 Département Génie Mécanique, Ecole National Polytechnique de Constantine ENPC, Constantine, Algérie

2025, Periodica Mathematica Hungarica

The order of accuracy of any convergent time integration method for systems of differential equations can be increased by using the sequence acceleration method known as Richardson extrapolation, as well as its variants (classical Richardson extrapolation and multiple Richardson extrapolation). The original (classical) version of Richardson extrapolation consists in taking a linear combination of numerical solutions obtained by two different time-steps with time-step sizes h and h/2 by the same numerical method. Multiple Richardson extrapolation is a generalization of this procedure, where the extrapolation is applied to the combination of some underlying numerical method and the classical Richardson extrapolation. This procedure increases the accuracy order of the underlying method from p to p + 2, and with each repetition, the order is further increased by one. In this paper we investigate the convergence of multiple Richardson extrapolation in the case where the underlying numerical method is an explicit Runge-Kutta method, and the computational efficiency is also checked.

2025, Journal of Geophysical Research: Oceans

A vertical gradient of suspended sediment concentration often exists in estuaries, particularly within the bottom boundary layer where sediment erosion and deposition take place. This results in a vertical density gradient and hence modification of the flow. However, this important effect has often been ignored in past sediment studies. Because of this and because of other empirical assumptions, exisiting erosion models cannot be used as predictive tools. This paper employs a simplified second‐order closure model to simulate the effect of sediment‐induced stratification on bottom boundary layer dynamics, and particularly the erosion process. Numerical models, which employ the concept of Richardson number dependent eddy viscosity, have been developed in the past for stratified flows. These models require a large number of data for tuning parameters of the eddy viscosity formula. The simplified second‐order closure model used here consists of the dynamic equations of motion for mean v...

2025, Astronomy & Astrophysics

In this paper we present two new results. The first concerns the proper identification of the critical Richardson number Ri(cr) above which there is no longer turbulent mixing. Thus far, all studies have assumed that: Ri(cr) = Ri (cr) =... more

2025, Dynamics of Atmospheres and Oceans

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

2025, Journal of Physical Oceanography

Three autonomous profiling Electromagnetic Autonomous Profiling Explorer (EM-APEX) floats were air deployed one day in advance of the passage of Hurricane Frances (2004) as part of the Coupled Boundary Layer Air–Sea Transfer (CBLAST)-High field experiment. The floats were deliberately deployed at locations on the hurricane track, 55 km to the right of the track, and 110 km to the right of the track. These floats provided profile measurements between 30 and 200 m of in situ temperature, salinity, and horizontal velocity every half hour during the hurricane passage and for several weeks afterward. Some aspects of the observed response were similar at the three locations—the dominance of near-inertial horizontal currents and the phase of these currents—whereas other aspects were different. The largest-amplitude inertial currents were observed at the 55-km site, where SST cooled the most, by about 2.2°C, as the surface mixed layer deepened by about 80 m. Based on the time–depth evolutio...

2025, MATEC Web of Conferences

In the present work, laminar mixed convection of a Newtonian fluid around a hot obstacle in a square cavity with moving vertical walls is studied numerically. The objective of this study is to analyze the effect of the Richardson number (0 ≼ Ri ≼ 10) and Reynolds number (50 ≼ Re ≼ 500) on both hydrodynamic and thermal characteristics around a hot obstacle in the enclosure. The analysis of the obtained results shows that the heat transfer is enhanced for high values of Richardson and Reynolds numbers.

2025, SSRN Electronic Journal

In the present work, laminar mixed convection of TiO 2 -Water nanofluid around a hot obstacle in a square cavity with moving vertical walls is studied numerically. The objective of this study is to analyze the effect of the Richardson number (0 ≼Ri≼ 10), Reynolds number (50 ≼ Re ≼ 500) and the nanoparticles volume fraction (0% ≼ ≼ 5%) on both hydrodynamic and thermal characteristics around a hot obstacle in the enclosure. The analysis of the obtained results shows that the heat transfer is enhanced for high values of Richardson and Reynolds numbers. In addition, the volume fraction of nanoparticles has a significant effect on the heat transfer within the cavity.

2025, Journal of Thermal Analysis and Calorimetry

The present numerical study consists of a mixed convection heat transfer, taking place within a ventilated cubic cavity crossed by an alumina-water nanofluid. The ventilation is assured by three openings, practiced on the walls of the cavity. So, the cold nanofluid enters by an opening placed at the top of the left vertical wall and exits by two opening: one is placed at the bottom of the right vertical wall, and the second opening occupied several possibility of location along the four walls of the cavity. All the enclosure's walls are maintained at a same temperature, higher than that of the entering fluid, except of the side walls which are adiabatic. The governing equations are discretized using the finite volume method and the SIMPLER algorithm to treat the coupling velocity-pressure. The method line by line is used to solve iteratively the algebraic equations. The simulations were carried for different values of the Richardson number (0 ^Ri ^10), the Hartmann number (0 ^Ha ^100) and a solid volume fractions u = 4% using Koo-Kleinstreuer-Lee model for the evaluation of effective thermal conductivity and dynamic viscosity of nanofluid. The results are presented in the form of streamlines, isotherms, velocity vectors and average Nusselt number.

2025

The present work relates a numerical study of mixed convection around a horizontal cylinder heated and equipped with fins and placed in a cylindrical cavity. In addition, the length of the cylinder is assumed to be infinite. Moreover, the heat transfer by radiation from the hot surface is assumed to be negligible, buoyancy effects are also considered, with the Boussinesq approximation. The forced circulation is induced by the inner cylinder in rotation with an angular velocity (ω), with its axis at the center of the annular space. Searches are made for various combinations of dimensionless numbers; Number of Reynolds (Re), number of Grashof (Gr). A finite volume scheme is used to solve the equations governing our system, which are continuity, bidimensional and energy, by the SIMPLE algorithm. The results found indicates that the heat exchange on the outer and inner cylinder characterized by the number of Nusselt decreases by increasing the Reynolds number and for low Reynolds number...

2025, International Journal of Heat and Mass Transfer

In this work, mixed convection in a Bingham plastic fluid from a heated hemi-sphere has been investigated numerically in the so-called aiding buoyancy configuration. The results reported herein, spanning the following ranges of conditions: Reynolds number, 1 6 Re 6 100, Prandtl number, 1 6 Pr 6 100, Richardson number, 0 6 Ri 6 10 and Bingham number, 0 6 Bn 6 100, are used to describe the momentum and heat transfer characteristics of a hemi-sphere oriented with its curved surface towards the oncoming fluid stream. In particular, flow characteristics are analyzed in terms of the streamlines and the morphology of fluid-like/solid-like regions and the corresponding heat transfer characteristics are examined via isotherm contours, local Nusselt number and average Nusselt number values to delineate the role of each of these parameters. The flow separates itself from the hemi-sphere at a much lower Reynolds numbers as compared to that for a sphere due to the sudden loss of the rear surface. The average value of the Nusselt number is influenced by four non-dimensional parameters, namely, Reynolds number, Prandtl number, Bingham number and Richardson number. Finally, the present results have been consolidated using the Colburn j-factor for predicting the value of the average Nusselt number in a new application.

2025, arXiv (Cornell University)

A Finite-Volume based POD-Galerkin reduced order modeling strategy for steady-state Reynolds averaged Navier-Stokes (RANS) simulation is extended for low-Prandtl number flow. The reduced order model is based on a full order model for which the effects of buoyancy on the flow and heat transfer are characterized by varying the Richardson number. The Reynolds stresses are computed with a linear eddy viscosity model. A single gradient diffusion hypothesis, together with a local correlation for the evaluation of the turbulent Prandtl number, is used to model the turbulent heat fluxes. The contribution of the eddy viscosity and turbulent thermal diffusivity fields are considered in the reduced order model with an interpolation based data-driven method. The reduced order model is tested for buoyancy-aided turbulent liquid sodium flow over a vertical backward-facing step with a uniform heat flux applied on the wall downstream of the step. The wall heat flux is incorporated with a Neumann boundary condition in both the full order model and the reduced order model. The velocity and temperature profiles predicted with the reduced order model for the same and new Richardson numbers inside the range of parameter values are in good agreement with the RANS simulations. Also, the local Stanton number and skin friction distribution at the heated wall are qualitatively well captured. Finally, the reduced order simulations, performed on a single core, are about 10 5 times faster than the RANS simulations that are performed on eight cores.

2025, CFD Letters

A numerical study of double-diffusive mixed convection within a horizontal rotating annulus has been investigated. The outer cylinder is fixed but the inner cylinder is considered to rotate in clockwise and anti-clockwise directions to introduce the forced convection effect. In addition, the solutal and thermal buoyancy forces are sustained by maintaining the inner and outer cylinder at uniform temperatures and concentrations but their values for the inner are higher than the outer. The flow is considered laminar regime under steady state conditions. The transport equations for the continuity, momentum, energy and mass transfer are solved using the finite volume technique. The considered domains in this investigation are: -15 ≤ N ≤ 15, 0.01 ≤ Ri ≤ 100 and 0.01 ≤ Pr ≤ 100. While the thermal Grashof number, Lewis number and the radius ratio are kept constant at values equal to10 4 , 1 and 2 respectively. The effect of the selected parameters on the local and average Nusselt and Sherwood numbers are presented and studied. Finally, this investigation concerned with selection the best direction of the inner cylinder rotation to enhance both heat and mass transfer. A comparison was made with the published results and a good agreement was found.

2025, WIT transactions on modelling and simulation

In this work mixed convection in a horizontal channel with the lower wall heated at a uniform heat flux is studied experimentally. The experiments are performed in air. The analysis is accomplished for several heat fluxes and forced air velocities. The Reynolds numbers investigated are between 10 and 500, these being in the laminar regime. The Richardson number Ri=Gr/Re 2 holds values in the range 0.1-100. Flow visualization is performed to detect the flow patterns into the channel. The wall temperature distribution is given as a function of the Gr/Re 2 values.

2025, Advances in Mechanical Engineering

Convective heat transfer can be enhanced passively by changing flow geometry and boundary conditions or by improving the thermal conductivity of the working fluid, for example, introducing suspended small solid nanoparticles. In this paper, a numerical investigation on laminar mixed convection in a water-Al 2 O 3 -based nanofluid, flowing in a triangular cross-sectioned duct, is presented. The duct walls are assumed at uniform temperature, and the single-phase model has been employed in order to analyze the nanofluid behaviour. The hydraulic diameter is equal to 0.01 m. A fluid flow with different values of Richardson number and nanoparticle volume fractions has been considered. Results show the increase of average convective heat transfer coefficient and Nusselt number for increasing values of Richardson number and particle concentration. However, also wall shear stress and required pumping power profiles grow significantly.

2025, Heat Transfer Engineering

Laminar mixed convection in a two-dimensional symmetrically and partially heated vertical channel is investigated. The heaters are located on both walls and uniform temperature is applied on the heated sections. The number of heaters is considered as 1, 4, 8, and 10. Aluminum oxide/water nanofluid is considered as working fluid and the inlet velocity is uniform. The continuity, momentum and energy equations with appropriate boundary conditions are solved in dimensionless form, numerically. The study is performed for Richardson number of 0.01 and 10, Reynolds number of 100 and 500, and nanofluid volume fraction of 0% and 5%. Based on the obtained velocity and temperature distributions, the local and mean Nusselt number is calculated and plotted for different cases. The variation of the mean Nusselt number with the number of the heated portions is also discussed. It is found that the addition of nanoparticles into the base fluid increases mean Nusselt number but the rate of increase depends on Reynolds, Richardson numbers and number of heated portions. It is possible to increase mean Nusselt number 138% by increasing Reynolds number from 100 to 500, Richardson number from 0.01 to 10 and number of heated portions from 1 to 10 when volume fraction value is 5%.

2025, Al-Qadisiyah Journal for Engineering Sciences

In this work numerical investigation of mixed convection by non-concentric different positions of rotating circular cylinder in square enclosure is investigated. The assumption of using six rotating cylinders at constant temperature Th and square closed enclosure with temperature Tc for all walls are maintained. Continuity, Momentum and energy governing equations are solved by using finite volume technique using Ansys-Fluent 16 commercial code with two dimension and steady state case. Richardson number (Ri=Gr/Re 2 ) varying over the wide range of 0.1, 1, 10,100, 1000 and ∞. This research investigate the effect of changing the rotating circular cylinder positions along the vertical centerline on the fluid flow and heat transfer inside the enclosure. Six different positions are tested. The phenomenon inside the square enclosure is mathematically analyzed. A streamlines, isothermal patterns, local Nusselt numbers, and average Nusselt numbers are calculated. The analytical results give that, the Richardson numbers is playing the dominant role on temperature distributions and flow patterns inside the enclosure. The average Nusselt number variation demonstrates the best location of rotating cylinder which gives best heat transfer value. The bottom left corner (P4) is the best location for hot rotating cylinders. While the rotating in middle right (P6) is the inferior position.

2025, CFD Letters

Understanding mixed convection in engineering applications such as heat exchangers, electronics cooling devices, and solar energy collectors have urged researchers to investigate this phenomenon deeper. This study investigates the fluid flow and heat transfer pattern in a two-dimensional (2D) rectangular cavity with sinusoidal heating on the moving top lid numerically. The bottom wall is kept cool while the vertical walls are insulated. The effect of Hartmann number, 𝐻𝑎 on the thermal characteristics and fluid flow are analyzed for Richardson number, 𝑅𝑖 = 1 which indicate mixed convection dominated regime. The governing equations are solved numerically using a SIMPLE algorithm with the finite volume method. The numerical results are displayed in streamlines and isotherms plots. The value of the Nusselt number indicating the heat transfer rate is also discussed. It is found that 𝐻𝑎 has a significant effect on the heat transfer process and fluid flow. It can be seen clearly when the value of 𝐻𝑎 = 30, the rate of heat transfer dropped significantly on the cold wall. Generally, the heat transfer rate decreases with the increase of 𝐻𝑎.

2025, International Journal of Mechanical Sciences

Ching-Chang Cho , Mixed convection heat transfer and entropy generation of Cu-water nanofluid in wavy-wall lid-driven cavity in presence of inclined magnetic field, International

2025, Journal of Geophysical Research

In November-December 1984 we carried out an intensive 12-day upper ocean sampling program on the equator at 140øW as part of the Tropic Heat Experiment. From our observations we constructed hourly averaged profiles of temperature, salinity, at, turbulent kinetic energy dissipation rate, and horizontal velocity. These data were used to examine the correspondence between hydrographic and velocity fields and to compare the measured turbulent dissipations with the calculated Richardson numbers. We found that the core of the Equatorial Undercurrent tracked a density surface (a t = 25.25) on times as short as 1 hour. The variability in both hydrographic and velocity fields was greatest at the semidiurnal frequency. The supertidal energy was not significantly different from the Garrett-Munk mid-latitude level once latitudinal scaling was removed from the Garrett-Munk model parameters. Horizontal velocity spectra were found to be contaminated by displacement of the background shear. Turbulent dissipation was dominated by a dirunal cycle, with high values of dissipation occurring at night above the undercurrent core. Shear and buoyancy frequency, calculated over 12-m vertical scales, were observed to track each other above the core and were dominated by a diurnal period above 40 m and by a semidiurnal period below 40 m. When shear and buoyancy frequency were combined to form a Richardson number, neither diurnal nor semidiurnal cycles were present. Above the undercurrent core, the Richardson numbers were uniformly small (0.3 to 0.6). 1. INTRODUCTION Studies have suggested that space and time variability in the equatorial ocean may be quite different from variability at mid-latitudes. Internal wave spectra from the equatorial Indian Ocean indicate more energy at frequencies above the tidal than is predicted by the Garrett-Munk (GM) universal internal wave spectrum [Eriksen, 1980]. Other deviations from GM [Garrett and Munk, 1972, 1975] include less spatial coherence [Wunsch and Webb, 1979; Eriksen, 1980] and an excess of horizontal kinetic energy over potential energy at subtidal frequencies [Eriksen, 1980]. One obvious difference is that at the equator there is no Coriolis parameter to limit the frequency bandwidth of internal waves. From their fine-structure measurements, Toole and Hayes [1984] found enhanced shear and strain variance and a greater proportion of low values of the Richardson number on the equator. McPhaden [1985] noted anomalously high fine-scale temperature and density variance confined within 1 ø of the equator. Microstructure measurements indicated a peak in turbulent kinetic energy dissipation rate within 1 ø of the equator [Crawford, 1982], but more recent results [Moumet al., 1986b-I have shown this peak to be an artifact of the sampling limitations of the previous study. In November-December 1984 we carried out an intensive 12-day sampling program near the equator at 140øW as part of the Tropic Heat experiment. Vertical profiles of temperature, conductivity, and small-scale shears were obtained every 10 min (on average) with the rapid sampling vertical Copyright 1986 by the American Geophysical Union. Paper number 6C0326. 0148-0227/86/006C-0326505.00 profiler (RSVP). Vertical profiles of horizontal currents were obtained every 30 s with a shipboard acoustic Doppler current profiler (ADCP). From these data we constructed hourly averaged profiles of temperature T, salinity S, at, kinetic energy dissipation rate e, and horizontal velocity. Using these series of profiles we made the following observations and comparisons: 1. We examined the correspondence between hydrographic and velocity fields. We found that the core of the Equatorial Undercurrent (EUC) tracked a density surface (a t = 25.25) on times as short as 1 hour. The core tracked the salinity maximum almost as well. 2. We compared the high-frequency variability with that observed at other latitudes in terms of the Garrett-Munk model. We found that (1) the time variability in both hydrographic and velocity fields was greatest at the semidiurnal frequency, (2) the supertidal energy was not significantly different from the Garrett-Munk mid-latitude level once the latitudinal scaling was removed from the Garrett-Munk model parameters, and (3) the horizontal velocity spectra were contaminated by displacement of the background shear. 3. We compared calculated Richardson numbers with the observed values of turbulent kinetic energy dissipation rate and found that (1) the dissipation was dominated by a diurnal cycle, with high values occurring at night above the undercurrent core, (2) shear and buoyancy frequency calculated over 12-m scales tracked each other extremely well, with a diurnal cycle above 40 m and a semidiurnal cycle below 40 m and (3) the Richardson number was uniformly small (between 0.3 and 0.6) above the undercurrent core. For the 12-day period, the Richardson number displayed neither a diurnal nor a semidiurnal periodicity. (In longer times series, periods were observed when the Richardson number cycled diurnally near the surface.) 12,887 12,888

2025, Environmental Fluid Mechanics

The measurements taken during the Vertical Transport and Mixing Experiment (VTMX, October, 2000) on a northeastern slope of Salt Lake Valley, Utah, were used to calculate the statistics of velocity fluctuations in a katabatic gravity current in the absence of synoptic forcing. The data from ultrasonic anemometer-thermometers placed at elevations 4.5 and 13.9 m were used. The contributions of small-scale turbulence and waves were isolated by applying a high-pass digital (Elliptical) filter, whereupon the filtered quantities were identified as small-scale turbulence and the rest as internal gravity waves. Internal waves were found to play a role not only at canonical large gradient Richardson numbers (Ri g > 1), but sometimes at smaller values (0.1 < Ri g < 1), in contrast to typical observations in flatterrain stable boundary layers. This may be attributed, at least partly, to (critical) internal waves on the slope, identified by Princevac et al. , which degenerate into turbulence and help maintain an active internal wave field. The applicability of both Monin-Obukhov (MO) similarity theory and local scaling to filtered and unfiltered data was tested by analyzing rms velocity fluctuations as a function of the stability parameter z/L, where L is the Obukhov length and z the height above the ground. For weaker stabilities, z/L < 1, the MO similarity and local scaling were valid for both filtered and unfiltered data. Conversely, when z/L > 1, the use of both scaling types is questionable, although filtered data showed a tendency to follow local scaling. A relationship between z/L and Ri g was identified. Eddy diffusivities of momentum K M and heat K H were dependent on wave activities, notably when Ri g > 1. The

2025, Journal of Sea Research

A conceptual model of the physical behavior of a shallow (6 m deep) micro-tidal estuary (Alfacs Bay) is proposed, based on the interpretation of a field data set, and subsequently tested against the results of three-dimensional hydrodynamic simulations. At seasonal timescales, the buoyancy associated with freshwater inflows dominates the tidal forcing, yielding a strongly stratified two-layered system, with the surface and the bottom layers flowing in opposite directions (classical estuarine circulation). Wind controls the physical behavior of the bay at shorter (days to weeks) timescales. Three scenarios or states have been defined, depending on the strength of stratification relative to the predominant direction and magnitude of the wind forcing, parameterized through the Wedderburn number, W. For weak winds (scenario 1), with W ≫ 1, mixing occurs as a consequence of stirring and convective cooling, and the mixed layer deepens slowly. For strong winds, with W ≪ 1/2, mixing is fast and is largely driven by shear at the pycnocline. Two scenarios are further identified for W ≪ 1/2 depending on the directionality of the winds: one for persistent NW winds (scenario 2) and another for diurnal SW winds (scenario 3). In scenario 2, the water is pushed laterally, overturning the stratification and generating transverse density gradients. In scenario 3, the estuarine circulation is weakened and even reversed, yielding strong longitudinal density gradients. The bay waters relax quickly (within 10-18 h) back to the original state, after the wind ceases, as a result of the horizontal density gradients developed under wind forcing. Bay-sea exchange rates are shown to decrease significantly in the low W scenarios, the magnitudes of these changes being largely dependent on wind direction. These scenarios have important implications for the ecology of the bay, including the occurrence of phenomena such as harmful algal blooms.

2025, ZAMP Journal of Applied Mathematics and Physics

2025, Atmospheric Science Letters

The critical Richardson number, Ric, is used in studies of stably stratified turbulence as a measure of flow laminarization. The accepted range of Ric is between 0.2 and 1. A growing body of experimental and observational data indicates, however, that turbulence survives for Ri ≫ 1. This result is supported by a new spectral theory of turbulence that accounts for strong anisotropy and waves. The anisotropization results in the enhanced horizontal mixing of both momentum and scalar. Internal wave contribution preserves vertical momentum mixing above its molecular level. In the absence of laminarization, Ric becomes devoid of its conventional meaning. Copyright © 2007 Royal Meteorological Society

2025, Physical science international journal

The mixed convection of heat transfer and fluid flow in a -lid drivencubical cavity filled with air is investigated numerically in this study. The computational procedure is based on the finite volume method and a full multigrid acceleration solver. The top wall of the cavity is maintained at a constant high -temperature Th, and it can move with a constant velocity U0. The bottom wall is immobile and maintained at a cold temperature Tc. While, the remaining boundary parts of the cavity are motionless and kept thermally insulated. Several numerical simulations were conducted to investigate mixed convection heat transfer in a sliding cubical cavity for a range of Reynolds numbers from 1000 to 5000 and Richardson numbers from 0.001 to 10. The influence of mixed convection parameters, Reynolds number, Richardson number, and heat transfer rate on the flow behavior was analyzed through parametric studies. The results include flow and heat transfer characteristics, iso-surfaces, and streamlines for the entire range of Richardson numbers and Reynolds numbers investigated. The study shows that as Reynolds number is increased beyond a critical value, the flow becomes unstable and bifurcates.

2025, Physical science international journal

2024, Hydrology Research

Currents and vertical mixing characteristics were investigated on the basis of time series of current meter and temperature data from a summer-stratified period in Lake Ontario. The experimental set up consisted of seven current meters distributed in one vertical line from 12 m below the surface to 1 m above the lake bottom at a total depth of 143 m. The period considered for the analysis was from June to September, 1991. The currents showed pronounced oscillations with two significant kinetic energy peaks, one at about 17 hours due to inertial motions, and one at 10 days, probably due to meteorological forcing. The current shear in the hypolimnion was strong enough to overcome stability and generate turbulence (Richardson numbers below 0.25) and there was probably turbulence enough available to keep the matter (almost neutral buoyant particles) in the whole Nepheloid bottom layer in suspension. In the thermocline region the turbulence was mainly damped (Richardson numbers above I ) , but some events with lower Richardson numbers were also calculated indicating increased mixing during these events. By analysing filtered and unfiltered current meter data it was found that the shear-generated turbulence in the hypolimnion was mainly due to the meteorologically forced currents. In the thermocline region, however, the vertical shcar associated with the inertial oscillation had a greater impact on the mixing.

2024, Thermal Science

In the present study, numerical simulation of magnetohydrodynamic (MHD) mixed convection heat transfer and fluid flow has been analyzed in a lid-driven enclosure provided with a constant flux heater. Governing equations were solved via differential quadrature (DQ) method. Moving wall of the enclosure has constant temperature and speed. The calculations were performed for different Richardson number ranging from 0.1 to 10, constant heat flux heater length from 0.2 to 0.8, location of heater center from 0.1 to 0.9, Hartmann number from 0 to 100 and aspect ratio from 0.5 to 2. Two different magnetic field directions were tested as vertical and horizontal. It was found that results of DQ method show good agreement with the results of literature. The magnetic field was more effective when it applied horizontally than that of vertical way. In both direction of magnetic field, it reduced the flow strength and heat transfer. Thus, it can be used as an important control parameter for heat an...

2024, Boundary-Layer Meteorology

It is shown that the observationally determined roughness relation z,, = au$g in which g is the acceleration of gravity, u* is the friction velocity in air, and a = 0.0185 (Wu, 1982) for the wind profile over the sea surface relative to the surface current, is consistent with the existence of a Richardson Number criterion at the air-sea interface in which the critical Richardson Number, Ri, = 1, such that all the shear energy is converted into potential energy.

2024, WIT Transactions on the Built Environment

The preliminary results of a statistical analysis of the atmospheric stability effects on the growth of sea waves are given in this paper. In unstable conditions (air colder than water) an increase of vertical mixing is experienced due to turbulence which makes the wind profile more uniform with height (deviating it from the logarithmic law). This circumstance increases the wind stress on the sea surface and thus gives rise to higher wave energy. In order to quantify the wave energy increase, the analysis of wave, wind and temperature data was performed, on the basis of the SWAN (Sea WAve monitoring Network) data set regarding both the Tyrrhenian and Adriatic Sea. The statistical analysis was first performed in terms of relative deviation of non-dimensional wave energy as a function of the Bulk Richardson Number in terms of the wind velocity U 10. Then another stability parameter, the Monin-Obukhov length, was related to the friction velocity u *. This approach needed the implementation of a numerical procedure to estimate the set of Monin-Obukhov parameters, which account for atmospheric stability.

2024, Journal Paper

This study investigates the cooling of a central processing unit (CPU) using a nano-encapsulated phase change material (NEPCM)-water mixture in a trapezoidal cavity with rotating cylinders and baffles. A numerical model based on the finite element method (FEM) is employed to solve the governing equations. The system is subjected to a sinusoidal temperature profile from the CPU and a constant magnetic field. Key parameters examined include Reynolds number (Re: 10-100), Richardson number (Ri: 0.1-10), Hartmann number (Ha: 5-80), NEPCM volume fraction (ϕ: 0.015-0.035), Lewis number (Le: 0.1-10), buoyancy ratio (Nz: 1-5), NEPCM fusion temperature (θ f : 0.1-0.9), and Stefan number (Ste: 0.1-0.9). Results show that increasing Re and Ri significantly enhances heat and mass transfer, with the average Nusselt number (Nu av) increasing by up to 80.5 % and average Sherwood number (Sh av) by up to 147.9 %. The magnetic field suppresses convection, reducing Nu av by 12.7 % and Sh av by 39.5 % as Ha increases. Increasing ϕ improves heat transfer (Nu av up by 32.5 %) with minimal effect on mass transfer. Le strongly influences mass transfer, with Sh av increasing by 284.6 % as Le increases. The NEPCM fusion temperature exhibits a non-monotonic effect on Nu av , with an optimal value at θ f = 0.5. In conclusion, the study reveals complex interactions between parameters, with Re, Ri, and Le having the most significant impacts on system performance. These findings provide valuable insights for optimizing CPU cooling systems using NEPCM-water mixtures and magnetohydrodynamic (MHD) effects.

2024, Energy and Thermofluids Engineering

Fluid flow in the lid driven cavity is a well-known phenomenon in the realm of fluid flow and heat transfer. Different cavity shapes such as square, circular, trapezoidal, and hexagonal have already been studied. However, most of the parametric analyses have previously been done using one specific cavity shape design. The present study intends to generalize the cavity design for a larger range of applications starting from the selection of cavity form due to the enhancement of heat transfer inside that selected cavity for different boundary conditions. Numerical studies using a finite element solver have been carried out to investigate the physics of fluid flow and heat transmission in a cavity where a wall of the enclosure is moving. In this research, the average Nusselt number for six distinct shaped lid-driven cavities is studied (e.g., square, rhombus, circular, rounded rectangular, trapezoidal and hexagonal). The hexagonal-shaped cavity has the largest Nusselt number of all of ...

2024, Quarterly Journal of the Royal Meteorological Society

Local similarity theory is suggested based on the Brunt-Väisälä frequency and the dissipation rate of turbulent kinetic energy instead the turbulent fluxes used in the traditional Monin-Obukhov similarity theory. Based on dimensional analysis (Pi theorem), it is shown that any properly scaled statistics of the small-scale turbulence are universal functions of a stability parameter defined as the ratio of a reference height z and the Dougherty-Ozmidov length scale which in the limit of z-less stratification is linearly proportional to the Obukhov length scale. Measurements of atmospheric turbulence made at five levels on a 20-m tower over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to examine the behaviour of different similarity functions in the stable boundary layer. It is found that in the framework of this approach the non-dimensional turbulent viscosity is equal to the gradient Richardson number whereas the non-dimensional turbulent thermal diffusivity is equal to the flux Richardson number. These results are a consequence of the approximate local balance between production of turbulence by the mean flow shear and viscous dissipation. The turbulence framework based on the Brunt-Väisälä frequency and the dissipation rate of turbulent kinetic energy may have practical advantages for estimating turbulence when the fluxes are not directly available.

2024, Journal of Advanced Research Design

The chemical reaction of heat and mass transfer in an unsteady magnetohydrodynamics (MHD) mixed convection stagnation point flow towards a vertical stretching plate embedded in porous medium was studied. Richardson number and concentration were taken into consideration in this study. By selecting suitable similarity variables, the non-linear partial differential equations (PDEs) are transformed into a set of ordinary differential equations (ODEs), and then solved numerically using the boundary value problem solver (bvp4c) in MATLAB software. The solution is dependent on the governing parameters including Prandtl number, reaction rate parameter, Schmidt number, thermal Richardson number and buoyancy ratio parameter. The numerical results are depicted graphically in velocity, temperature and concentration profiles. The findings have shown that as the chemical reaction parameter increases, the temperature profile and Sherwood number also increases.

2024, Applied and Computational Mechanics

In the present study, fluid flow, heat transfer, and entropy generation for mixed convection inside a water-filled square cavity were investigated numerically. The sidewalls of the cavity, which move upwards, are kept at low-temperature Tc while only a part in the center of the bottom wall is kept at high-temperature Th and the remaining parts are kept adiabatic. The governing equations, in stream function-vorticity form, are discretized and solved using the finite difference method. Particular attention was paid to the influence of the Prandtl numbers of 5.534, 3.045 and 2, corresponding respectively to the water temperatures of 303.15 K, 333.15 K and 363.15 K. The numerical results are presented in the form of streamlines, isotherms, and entropy generation contours for different values of the Richardson numbers at an arbitrary Reynolds number Re=10 2. Besides this, the evolution of the average Nusselt number and the average entropy generation is also reported. The obtained results show interesting behaviors of the flow and thermal fields, which mainly involve stable symmetric and non-symmetric steady-state solutions, as well as unsteady regimes, depending on specific values of the Richardson and Prandtl numbers. It is additionally observed that the average Nusselt number increases and the average entropy generation decreases when both the Richardson and Prandtl numbers increase.

2024, Springer eBooks

We present airplane measurements of the stably strattied nocturnal boundary layer obtained during the Severe Environmental Storms and Mesoscale Experiment (SESAME) in 1979. The cases presented here were obtained over rolling terrain in central Oklahama, with a mean slope of about 0.003. The results are in general agreement with previous modeling and observational studies for the mean and turbulence structure of the nocturnal boundary layer, with the exception that the eddy diffusivity of heat, and consequently the flux Richardson number are less than expected. DONALD H. LENSCHOW ET AL.

2024, International Journal of Numerical Methods for Heat & Fluid Flow

Purpose The purpose of this paper is to investigate the mixed convection of a two-phase water–aluminum oxide nanofluid in a cavity under a uniform magnetic field. Design/methodology/approach The upper wall of the cavity is cold and the lower wall is warm. The effects of different values of Richardson number, Hartmann number, cavitation length and solid nanoparticles concentration on the flow and temperature field and heat transfer rate were evaluated. In this paper, the heat flux was assumed to be constant of 10 (W/m2) and the Reynolds number was assumed to be constant of 300 and the Hartmann number and the volume fraction of solid nanoparticles varied from 0 to 60 and 0 to 0.06, respectively. The Richardson number was considered to be 0.1, 1 and 5. Aspect ratios were 1, 1.5 and 2. Findings Comparison of the results of this paper with the results of the numerical and experimental studies of other researchers showed a good correlation. The results were presented in the form of veloci...

2024, International Journal of Numerical Methods for Heat & Fluid Flow

Purpose-The purpose of this paper is to investigate the mixed convection of a two-phase water-aluminum oxide nanofluid in a cavity under a uniform magnetic field. Design/methodology/approach-The upper wall of the cavity is cold and the lower wall is warm. The effects of different values of Richardson number, Hartmann number, cavitation length and solid nanoparticles concentration on the flow and temperature field and heat transfer rate were evaluated. In this paper, the heat flux was assumed to be constant of 10 (W/m 2) and the Reynolds number was assumed to be constant of 300 and the Hartmann number and the volume fraction of solid nanoparticles varied from 0 to 60 and 0 to 0.06, respectively. The Richardson number was considered to be 0.1, 1 and 5. Aspect ratios were 1, 1.5 and 2. Findings-Comparison of the results of this paper with the results of the numerical and experimental studies of other researchers showed a good correlation. The results were presented in the form of velocity and temperature profiles, stream and isotherm lines and Nusselt numbers. The results showed that by increasing the Hartmann number, the heat transfer rate decreases. An increase from 0 to 20 in Hartmann number results in a 20 per cent decrease in Nusselt numbers, and by increasing the Hartmann number from 20 to 40, a 16 per cent decrease is observed in Nusselt number. Accordingly, it is inferred that by increasing the Hartmann number, the reduction in the Nusselt number is decreased. As the Richardson number increased, the heat transfer rate and, consequently, the Nusselt number increased. Therefore, an increase in the Richardson number results in an increase of the Nusselt number, that is, an increase in Richardson number from 0.1 to 1 and from 1 to 5 results in 37 and 47 per cent increase in Nusselt number, respectively. Originality/value-Even though there have been numerous investigations conducted on convection in cavities under various configurations and boundary conditions, relatively few studies are conducted for the case of nanofluid mixed convection in square lid-driven cavity under the effect of magnetic field using twophase model.

2024, Mağallaẗ al-qādisiyyaẗ li-l-ʻulūm al-handasiyyaẗ

Numerical investigation of mixed convective in a vented square cavity with fin. The horizontal walls are adiabatic, while the left and right walls are hot (ℎ) and cold () temperatures, respectively. The fluid inlet to the cavity from the lower left open area(), and exits from the upper right open area (). In this study, a finite element scheme is employed. The analysis is done for specific Prandtl number (= 7), Reynolds number (50 ≤ ≤ 200), fin length (0.2 ≤ ≤ 0.6), Richardson number (0.1 ≤ ≤ 1), and the location of the fin (0.2 ≤ ℎ ≤ 0.6). The finding indicates that the increases when high the location of the fin increases at the maximum height of this fin location is estimated to be 17% due to an increase in the fluid flow area on the hot wall caused by rising convective. The highest heat transfer occurs when the fin length equals 0.6 at the location(ℎ = 0.2).

2024

A series of laboratory experiments has been performed to investigate the flow characteristics of 2-dimensional density currents induced by selective withdrawal, which is commonly suggested as a measure for removal of high turbid water from reservoirs. Saltwater has been used to simulate the density stratification over depth and PIV(Particel Image Velocimetry) for observing the velocity structure. Experimental conditions have been established according to Richardson number, which is the dimensionless number that expresses the ratio of potential to kinetic energy. From the experiments, the patterns of longitudinal decay of centerline axial velocity induced by the withdrawal have been distinguished from other experimental cases. The rate of longitudinal decay increase as the Richardson number increases. The variations of volumetric and momentum flux along the longitudinal axis have also shown to be dependent on Richardson number.

2024

Atmospheric conditions during two weather related aircraft incidents in Iceland and elements of the climatology of windstorms

2024, Geophysical Research Letters

Most geophysical flows encompass turbulence and internal and/or Rossby waves. We demonstrate that these two different classes of waves cause remarkably similar anomalies in the turbulent transport. While all scales in both types of flows contribute to the momentum diffusion, the vertical (diapycnal) scalar diffusion in stratified flows and lateral diffusion in β‐plane turbulence can be carried out only by turbulent eddies whose size is smaller than the thresholds of turbulence anisotropization. Beyond these thresholds, both flows become dominated by waves that provide no contribution to the scalar diffusion. Stably stratified flows exhibit enhanced isopycnal diffusion of both momentum and scalar. These results shed new light on the Osborn mixing model, diapycnal and isopycnal viscosity and diffusivity, absence of the critical gradient Richardson number, and large scale meridional transport.

2024

A model tunnel (1:30 compared to a standard tunnel section) with a helium-air smoke mixture was used to study the vehicular blockage effect on longitudinal ventilation smoke control. The experimental results showed excellent agreement with full-scale data and confirmed that the critical velocity decreases in proportion with the blockage ratio. Nevertheless, it was found that the relative position of the fire source and the relative size of the vehicular blockage can have an opposite effect, as the vehicular blockage influenced the critical and confinement velocity. The method demonstrated the ability to provide valuable information on the effect of vehicular blockage on tunnel fire dynamics.

2024, Journal of Geophysical Research: Oceans

The evolving upper ocean response excited by the passage of hurricane Gilbert (September 14–19, 1988) was investigated using current and temperature observations acquired from the deployment of 79 airborne expendable current profilers (AXCPs) and 51 airborne expendable bathythermographs from the National Oceanic and Atmospheric Administration WP‐3D aircraft in the western Gulf of Mexico. The sea surface temperatures (SSTs), mixed layer depths, and bulk Richardson numbers were objectively analyzed to examine the spatial variability of the upper ocean response to Gilbert. Net decreases of the SSTs of 3°–4°C were observed by the profilers as well as by the airborne infrared thermometer (AIRT) along the flight tracks and advanced very high resolution radiometer (AVHRR) imagery. The AXCPs indicated a marked cooling from 29°C to about 25.5°C on September 17, 1988, which was about 1.2 inertial periods (IP) following storm passage. This pool of cooler water (3.5°) was located further downst...

2024, Thermal Science

In the aim to optimize the flow systems, a numerical simulation is carried out to investigate the entropy generation of the aiding mixed convection within two nonparallel vertical plates. The computational procedure is made by solving the laminar, steady and bi-dimensional continuity, momentum and energy equations with the finite volume method. The plates are symmetrically heated with uniform temperature. The mixed convection flow in the heated vertical channel is studied for two aiding buoyancy conditions (upstream flow between hot plates and downstream flow between cold plates). The calculations are performed for several parameters such as Richardson number or buoyancy parameter (Gr/Re 2 ), Reynolds and Bejan numbers. The results showed that the velocity, temperature and entropy generation profiles within the channel are significantly affected by the later parameters..

2024, Journal of Fluid Mechanics

Direct numerical simulations are performed to study the dynamics of an inhomogeneous stratified shear flow that models an atmospheric jet centred at the tropopause across which the density stratification is non-uniform. Small to moderate background stratifications are selected, and simulations are conducted for a range of Reynolds and Froude numbers. A spectral domain decomposition method that is particularly suitable for simulations of non-uniformly stratified shear flows is developed to simulate the desired turbulent jet, and quasi-equilibrium flow fields are obtained by long-time integration of governing equations. The structures of the mean flow and turbulence fields are calculated, which are interpreted using relevant length scales (Ozmidov, buoyancy, shear, Ellison) and Richardson number profiles. The ratios of the Ellison to buoyancy scales are much smaller than unity at the jet core and approach unity at the edges, confirming that mechanical turbulence prevails in the jet core, while nonlinear waves and stratification effects are dominating at the jet edges. The jet core is found to support sustained mechanical (active) turbulence, outside which lay a region of patchy turbulence and nonlinear gravity wave activity characterized by spatially decaying velocity fluctuations and strong temperature fluctuations. Detailed energy budgets show how energy is partitioned within the flow, including the transport of energy from the jet to its immediate vicinity by nonlinear gravity waves.

2024, Nanomedicine & Nanotechnology Open Access

The problem of MHD mixed convection is analyzed in a lid driven cavity with corrugated wavy bottom wall filled with Cu-H2 O nanofluid in presence of internal heat source. The top and right walls of the cavity are maintained with a uniform cold temperature whereas the left wall and bottom wavy wall are kept adiabatic. The top wall is moving with a constant velocity upon its lid and a rectangular heat source is placed horizontally inside the cavity. The physical problem is characterized by 2D governing partial differential equations along with proper boundary conditions and are discretized using Galerkin’s finite element formulation. The study is executed by analyzing different ranges of geometrical, physical and nondimensional parameters namely, wave number of wavy surface (0 ≤ ≤ λ 4) , the ratio of heat source height and cavity height 1 3 1 20 a l  ≤ ≤      volume fraction of nanoparticle (0 ≤ ≤ ϕ 0.09) Hartmann number (0 ≤ ≤ Ha 90) and Richardson number (0.1≤ ≤ Ri 10) . The r...

2024, Journal of Applied Meteorology

Vertical mixing of momentum and heat is investigated in turbulent stratified shear flows. It is assumed that the flow has uniform shear and stratification with homogeneous turbulence and that an equilibrium is reached between kinetic and potential energy without gravity wave oscillations. A simple model is derived to estimate vertical diffusivities for Richardson numbers in between 0 and about 1. The model is based on the budgets of kinetic and potential energy and assumes a linear relationship between dissipation, shear, and vertical velocity variance for closure. Scalar fluctuations are related to shear or buoyancy frequency depending on the Richardson number. The turbulent Prandtl number and the growth rate of kinetic energy are specified as functions of this number. Model coefficients are determined mainly from laboratory measurements. Data from large-eddy simulations are used to determine the "stationary" Richardson number with balanced shear production, dissipation, ...

2024, Journal of Thermal Engineering

This work presents a study by numerical simulation of mixed convection in a rectangular cavity with a sinusoidal temperature imposed on the right vertical wall while the other wall on the left is kept at a cold temperature. The upper and lower walls are thermally insulated, the inlet and outlet ports are respectively located on the hot wall to the bottom and on the top to the left. The enclosure represents a practical system such as an air-cooled electronic device, the heat source represents a radiator or an electronic component located in three different positions towards the left side at the bottom, the center and towards the right side at the top in such an enclosure. All calculations are made for a range of Richardson number from 0 to 10 and Reynolds numbers from 50 to 200. The influence of Richardson number, the position of the heat source and the influence of amplitude and phase deviation of the temperature imposed on the Nusselt number on the hot surface is studied. The resul...

2024, Physical science international journal

The mixed convection of heat transfer and fluid flow in a-lid drivencubical cavity filled with air is investigated numerically in this study. The computational procedure is based on the finite volume method and a full multigrid acceleration solver. The top wall of the cavity is maintained at a constant high-temperature Th, and it can move with a constant velocity U0. The bottom wall is immobile and maintained at a cold temperature Tc. While, the remaining boundary parts of the cavity are motionless and kept thermally insulated. Several numerical simulations were conducted to investigate mixed convection heat transfer in a sliding cubical cavity for a range of Reynolds numbers from 1000 to 5000 and Richardson numbers from 0.001 to 10. The influence of mixed convection parameters, Reynolds number, Richardson number, and heat transfer rate on the flow behavior was analyzed through parametric studies. The results include flow and heat transfer characteristics, iso-surfaces, and streamlines for the entire range of Richardson numbers and Reynolds numbers investigated. The study shows that as Reynolds number is increased beyond a critical value, the flow becomes unstable and bifurcates.