Evolution and decay of gravity wavefields in weak-rotating environments: a laboratory study

Gravity waves are prominent physical features that play a fundamental role in transport processes of stratified aquatic ecosystems. In a two-layer stratified basin, the equations of motion for the first vertical mode are equivalent to the linearised shallow water equations for a homogeneous fluid. We adopted this framework to examine the spatiotemporal structure of gravity wavefields weakly affected by the background rotation of a single-layer system of equivalent thickness \(h_{\ell }\), via laboratory experiments performed in a cylindrical basin mounted on a turntable. The wavefield was generated by the release of a diametral linear tilt of the air–water interface, \(\eta _{\ell }\), inducing a basin-scale perturbation that evolved in response to the horizontal pressure gradient and the rotation-induced acceleration. The basin-scale wave response was controlled by an initial perturbation parameter, \({\mathcal{A}}_{*} = \eta _{0}/h_{\ell }\), where \(\eta _{0}\) was the initial displacement of the air–water interface, and by the strength of the background rotation controlled by the Burger number, \({\mathcal{S}}\). We set the experiments to explore a transitional regime from moderate- to weak-rotational environments, \(0.65\le {\mathcal{S}} \le 2\), for a wide range of initial perturbations, \(0.05\le {\mathcal{A}}_{*}\le 1.0\). The evolution of \(\eta _{\ell }\) was registered over a diametral plane by recording a laser-induced optical fluorescence sheet and using a capacitive sensor located near the lateral boundary. The evolution of the gravity wavefields showed substantial variability as a function of the rotational regimes and the radial position. The results demonstrate that the strength of rotation and nonlinearities control the bulk decay rate of the basin-scale gravity waves. The ratio between the experimentally estimated damping timescale, \(T_{d}\), and the seiche period of the basin, \(T_{g}\), has a median value of \(T_{d}/T_{g}\approx 11\), a maximum value of \(T_{d}/T_{g}\approx 10^{3}\) and a minimum value of \(T_{d}/T_{g}\approx 5\). The results of this study are significant for the understanding the dynamics of gravity waves in waterbodies weakly affected by Coriolis acceleration, such as mid- to small-size lakes.     

Variations of bed elevations due to turbulence around submerged cylinder in sand beds

This paper presents the spatio-temporal variations in bed elevations and the near-bed turbulence statistics over the deformed bed generated around the submerged cylindrical piers embedded vertically on loose sediment bed at a constant flow discharge. Experiments were carried out in a laboratory flume for three blockage ratios in the range of 0.04–0.06 using three different sizes of submerged cylinders individually placed vertically at the centerline of the flume. Clear-water experimental conditions were maintained over the smooth sediment bed surface with a constant flow discharge (\(Q = 0.015\,{\rm m}^3/{\rm sec}\)), thereby giving three different cylinder Reynolds numbers \(Re_{D_c} = \frac{U_mD_c}{\nu }\) (=10200, 12750, 15300) away from the cylinder locations, where \(U_m\) is the maximum mean velocity, \(D_c\) is the cylinder diameter and \(\nu\) is the kinematic viscosity of fluid. Instantaneous sand bed elevations around the cylinders were recorded using a SeaTek 5MHz ultrasonic ranging system of net 24 transducers to estimate bed form migration, and the near-bed velocity data at transducer locations over the stable deformed bed around the pier-like structures were collected using down-looking three-dimensional (3D) Micro-acoustic Doppler velocimeter to estimate the bottom Reynolds shear stresses and the contributions of bursting events to the dominant shear stress component. The flow perturbation generated due to relatively lower flow blockage ratio favored to achieve the stable bed condition more rapidly than the others, and larger upstream scour-depth and deformed areas were noticed for greater flow blockage ratio due to larger cylinder diameter. For larger blockage ratio in the upstream of scour-hole near the bed, occurrences of probabilities of both boundary-ward interactions (Q1 and Q3) were the dominant; whereas in the downstream of the scoured region, occurrences of probabilities of second and third quadrant events (Q2 and Q4) were dominant. On the other hand, for the lower blockage ratio, quadrant (Q2) was dominant over Q4 in the downstream of scour-hole, and in the upstream of scour-hole, quadrant Q4 was the dominant.     

Effect of atmospheric boundary layer stability on the inclination angle of turbulence coherent structures

Coherent structures in the atmospheric boundary layer are fundamental to the transport of momentum and heat as well as to the production of turbulence. The present work attempts to investigate the behavior of the inclination angle of the vortex packet structures (\(\gamma\)) under different stability conditions. The data were collected from the Marine Ecosystem Research Centre (EKOMAR) site at the east coast of Peninsular Malaysia. The main measurements were conducted by placing two hotwires 3 and 12 m above ground. The two-point correlation method was used to calculate the vortex packet structure inclination angle, while the one-point correlation method was employed to calculate its length-scale. The inclination angle was found to increase under both stable and unstable conditions. As the Obukhov stability parameter (\(\zeta\)) approaches 0, the inclination angle ranged between \(\gamma = 15^\circ\) to \(\gamma = 18^\circ\) for the stable and unstable conditions, respectively, which agrees with the findings of previous research. The vertical gradient of velocity is the dominant parameter affecting the inclination angle under different stability conditions.     

Small-scale entrainment in inclined gravity currents

We investigate the effect of buoyancy on the small-scale aspects of turbulent entrainment by performing direct numerical simulation of a gravity current and a wall jet. In both flows, we detect the turbulent/nonturbulent interface separating turbulent from irrotational ambient flow regions using a range of enstrophy iso-levels spanning many orders of magnitude. Conform to expectation, the relative enstrophy isosurface velocity \(v_n\) in the viscous superlayer scales with the Kolmogorov velocity for both flow cases. We connect the integral entrainment coefficient E to the small-scale entrainment and observe excellent agreement between the two estimates throughout the viscous superlayer. The contribution of baroclinic torque to \(v_n\) is negligible, and we show that the primary reason for reduced entrainment in the gravity current as compared to the wall-jet are 1) the reduction of \(v_n\) relative to the integral velocity scale \(u_T\); and 2) the reduction in the surface area of the isosurfaces.     

On the front conditions for gravity currents in channels of general cross-section

We consider the propagation of a high-Reynolds-number gravity current in a horizontal channel with general cross-section whose width is \(f(z), 0 \le z\le H\), and the gravity acceleration g acts in \(-z\) direction. (The classical rectangular cross-section is covered by the particular case \(f(z) =\) const.) We assume a two-layer system of homogeneous fluids of constant densities \(\rho _{c}\) (current, of height \(h < H \)) and smaller \(\rho _{a}\) (ambient, filling the remaining part of the channel). We focus attention on the calculation and assessment of the nose Froude-number condition \(Fr = U/(g' h)^{1/2}\); here U is the speed of propagation of the current and \(g' = (\rho _{c}/\rho _{a}-1) g\) is the reduced gravity. We first revisit the steady-state current, and derive compact insightful expressions of Fr and energy dissipation as a function of \(\varphi \) (\(=\) area fraction occupied by the current in the cross-section). We show that the head loss \(\delta _0\) on the stagnation line is formally a degree of freedom in the determination of \(Fr(\varphi )\), and we clarify the strong connections with the head loss \(\delta \) in the ambient fluid, and with the overall rate of dissipation \(\dot{{\mathcal{D}}}\). We demonstrate that the closure \(\delta _0 = 0\) [suggested by Benjamin (J Fluid Mech 31, 209–248, 1968) for the rectangular cross-section] produces in general the smallest Fr for a given \(\varphi \); the results are valid for a significant range \([0, \varphi _{\max }]\), in which the current is dissipative, except for the point \(\varphi _{\max }\) where \(\delta = \dot{{\mathcal{D}}} = 0\). We show that imposing the closure \(\delta = \dot{{\mathcal{D}}} = 0\), which corresponds to an energy-conserving or non-dissipative current, produces in general unacceptable restrictions of the range of validity, and large values of Fr; in particular, deep currents (\(\varphi < 0.3\) say) must be excluded because they are inherently dissipative. On the other hand, the compromise closure \(\delta (\varphi ) =\delta _0(\varphi )\) produces the simple \(Fr(\varphi ) = \sqrt{2}(1 - \varphi )\) formula whose values and dissipation properties are very close, and the range of validity is identical, to these obtained with Benjamin’s closure (moreover, we show that this corresponds to circulation-conservation solutions). The results are illustrated for practical cross-section geometries (rectangle, \(\Delta \) and \(\nabla \) triangle, circle, and the general power-law \(f(z) = b z ^\alpha \) (\(b>0, \alpha \ge 0, 0< z \le H\)). Next, we investigate the connection of the steady-state results with the time-dependent current, and show that in a lock-released current the rate of dissipation of the system is equal to, or larger than, that obtained for Fr corresponding to the conditions at the nose of the current. The results and insights of this study cover a wide range of cross-section geometry and apply to both Boussinesq and non-Boussinesq systems; they reveal a remarkable robustness of Fr as a function of \(\varphi \).     

Thin-layer models for gravity currents in channels of general cross-section area,a review

We present a brief review of the recent investigations on gravity currents in horizontal channels with non-rectangular cross-section area (such as triangle, \(\bigvee \)-valley, circle/semi-circle, trapezoid) which occur in nature (e.g., rivers) and constructed environment (tunnels, reservoirs, canals). To be specific, we discuss the propagation of a gravity current (GC) in a horizontal channel along the horizontal coordinate x, with gravity g acting in the \(-z\) direction, and y the horizontal–lateral coordinate. The bottom and top of the channel are at \(z=0,H\). The “standard” problem is concerned with 2D flow in a channel with rectangular (or laterally unbounded) cross-section area (CSA). Recent investigations have successfully extended the standard knowledge to the channels of CSA given by the quite general \(-f_1(z)\le y \le f_2(z)\) for \(0 \le z \le H\). This includes the practical \(\bigvee \)-valley, triangle, circle/semi-circle and trapezoid; these geometries may be in “up” or “down” setting with respect to gravity, e.g., \(\bigtriangleup \) and \(\bigtriangledown \). The major objective of the extended theory is to predict the height of the interface \(z=h(x,t)\) and the velocity (averaged over the CSA) u(x, t), where t is time; the prediction includes the speed and position of the nose \(u_N(t), x_N(t)\). We show that the motion is governed by a set of simplified equations, called “model,” that provides versatile and insightful solutions and trends. The emphasis in on a high-Reynolds-number current whose motion is dominated by buoyancy–inertia balance; in particular a GC released from a lock, which also contains general effects such as front and internal jumps (shocks), and reflected bore. We discuss two-layer, one-layer, and box models; Boussinesq and non-Boussinesq systems; compositional and particle-driven cases; and the effect of stratification of the ambient fluid. The models are self-contained, and admit realistic initial and boundary conditions. The governing equations are amenable to analytical solutions in some special circumstances. Some salient features of the buoyancy-viscous regime, and the estimate for the length at which transition to this regime takes place, are also presented. Some experimental support to the theory, and open questions for further investigations, are also mentioned. The major conclusions are (1) The CSA geometry has significant influence on the motion of the GC; and (2) The new theory is a useful, very significant, extension of the standard two-dimensional GC problem. The standard current is just a particular case, \(f_{1,2} =\) constants, among many other covered by the new theory.     

Analysis and prediction of $$\mathrm{PM}_{10}$$ concentration levels in Tunisia using statistical learning approaches

Over the past years, the health impact of airborne particulate matter \(\mathrm{PM}_{10}\) has become a very topical subject. Thereby, a lot of research effort in the environmental sciences goes towards the modeling and the prediction of ambient \(\mathrm{PM}_{10}\) concentrations. In this paper, we are interested in the statistical classification of the daily mean \(\mathrm{PM}_{10}\) concentration in Tunisia according to the authority regulation. We consider two monitoring stations: a big industrial station and a traffic station. The main goal of this work is to determine the pertinent predictors of \(\mathrm{PM}_{10}\) concentration within a nonlinear multiclass framework. To do this, we used two popular statistical learning methods; the support vector machines (SVM) and the random forests (RF). The statistical results obtained on the real datasets, show that RF outperform SVM for the purpose of variable selection even with a reduced number of observations compared to the number of explicative variables. It was also demonstrated that the \(\mathrm{PM}_{10}\) concentration measured yesterday is the most relevant predictor of its present-day value. Moreover, we found that the more delayed values of \(\mathrm{PM}_{10}\) concentration may be crucial to get an accurate prediction.     

Wind-induced pressure at a tunnel portal

In order to properly size the mechanical ventilation system of a tunnel, it is essential to estimate the wind-driven pressure difference that might rise between its two portals. In this respect, we explore here the pressure distribution over a tunnel portal under the influence of an incident atmospheric boundary layer and, in particular, its dependency on wind direction and on tunnel geometry. Reduced scale models of generic configurations of a tunnel portal are studied in an atmospheric wind tunnel. Pressure distributions over the front section of different open cavities are measured with surface taps, which allows us to infer the influence of the tunnel aspect ratio and wind direction on a pressure coefficient \(C_{P}\), defined as a spatially and time averaged non-dimensional pressure. Experiments reveal that the magnitude of the coefficient \(C_{P}\), as a function of the wind direction, is significantly influenced by the portal height-to-width ratio and almost insensitive to its length. The experimental data set is completed by hot-wire anemometry measurements providing vertical distribution of velocity statistics. The same configurations are simulated by numerically solving the Reynolds-averaged Navier–Stokes equations, adopting the standard \(k - \varepsilon\) turbulence model. Despite some discrepancies between numerical and experimental estimates of some flow parameters (namely the turbulent kinetic energy field), the numerical estimates of the pressure coefficients \(C_{P}\) show very good agreement with experimental data. The latter is also compared to the predictions of an analytical model, based on the estimate of a spatially averaged velocity within an infinitely long street canyon. The results of the model, which takes into account varying canyon aspect ratios, are in reasonable agreement with experimental data for all cases studied. Notably, its predictions are significantly better than those provided by the simple analytical relations usually adopted as a reference in tunnel ventilation studies.     

Development of similarity relationships for energy dissipation rate and temperature structure parameter in stably stratified flows: a direct numerical simulation approach

In this study, a newly developed direct numerical simulation (DNS) solver is utilized for the simulations of numerous stably stratified open-channel flows with bulk Reynolds number (Re _b ) spanning 3400–16,900. Overall, the simulated bulk Richardson number (\(Ri_b\)) ranges from 0.08 (weakly stable) to 0.49 (very stable). Thus, both continuously turbulent and (globally) intermittently turbulent cases are represented in the DNS database. Using this comprehensive database, various flux-based and gradient-based similarity relationships for energy dissipation rate (ε) and temperature structure parameter (\(C_T^2\)) are developed. Interestingly, these relationships exhibit only minor dependency on Re _b . In order to further probe into this Re _b -effect, similarity relationships are also estimated from a large-eddy simulation (LES) run of an idealized atmospheric boundary layer (very high Re _b ) case study. Despite the fundamental differences in the estimation of ε and \(C_T^2\) from the DNS- and the LES-generated data, the resulting similarity relationships, especially the gradient-based ones, from these numerical approaches are found to be remarkably similar. More importantly, these simulated relationships are also comparable, at least qualitatively, to the traditional observational data-based ones. Since these simulated similarity relationships do not require Taylor’s hypothesis and do not suffer from mesoscale disturbances and/or measurement noise, they have the potential to complement the existing similarity relationships.     

Simulation-based optimization of in-stream structures design: rock vanes

We employ a three-dimensional coupled hydro-morphodynamic model, the Virtual Flow Simulator (VFS-Geophysics) in its Unsteady Reynolds Averaged Navier–Stokes mode closed with \(k-\omega\) model, to simulate the turbulent flow and sediment transport in large-scale sand and gravel bed waterways under prototype and live-bed conditions. The simulation results are used to carry out systematic numerical experiments to develop design guidelines for rock vane structures. The numerical model is based on the Curvilinear Immersed Boundary approach to simulate flow and sediment transport processes in arbitrarily complex rivers with embedded rock structures. Three validation test cases are conducted to examine the capability of the model in capturing turbulent flow and sediment transport in channels with mobile-bed. Transport of sediment materials is handled using the Exner equation coupled with a transport equation for suspended load. Two representative meandering rivers, with gravel and sand beds, respectively, are selected to serve as the virtual test-bed for developing design guidelines for rock vane structures. The characteristics of these rivers are selected based on available field data. Initially guided by existing design guidelines, we consider numerous arrangements of rock vane structures computationally to identify optimal structure design and placement characteristics for a given river system.     

The CPCAT as a novel tool to overcome the shortcomings of NOEC/LOEC statistics in ecotoxicology: a simulation study to evaluate the statistical power

Species reproduction is an important determinant of population dynamics. As such, this is an important parameter in environmental risk assessment. The closure principle computational approach test (CPCAT) was recently proposed as a method to derive a NOEC/LOEC for reproduction count data such as the number of juvenile Daphnia. The Poisson distribution used by CPCAT can be too restrictive as a model of the data-generating process. In practice, the generalized Poisson distribution could be more appropriate, as it allows for inequality of the population mean \(\mu\) and the population variance \(\sigma ^2\). It is of fundamental interest to explore the statistical power of CPCAT and the probability of determining a regulatory relevant effect correctly. Using a simulation, we varied between Poisson distribution (\(\mu =\sigma ^2\)) and generalized Poisson distribution allowing for over-dispersion (\(\mu <\sigma ^2\)) and under-dispersion (\(\mu >\sigma ^2\)). The results indicated that the probability of detecting the LOEC/NOEC correctly was \(\ge 0.8\) provided the effect was at least 20% above or below the mean level of the control group and mean reproduction of the control was at least 50 individuals while over-dispersion was missing. Specifically, under-dispersion increased, whereas over-dispersion reduced the statistical power of the CPCAT. Using the well-known Hampel identifier, we propose a simple and straight forward method to assess whether the data-generating process of real data could be over- or under-dispersed.     

Fitting the truncated negative binomial distribution to count data

Modeling empirical distributions of repeated counts with parametric probability distributions is a frequent problem when studying species abundance. One must choose a family of distributions which is flexible enough to take into account very diverse patterns and possess parameters with clear biological/ecological interpretations. The negative binomial distribution fulfills these criteria and was selected for modeling counts of marine fish and invertebrates. This distribution depends on a vector \(\left( K,\mathfrak {P}\right) \) of parameters, and ranges from the Poisson distribution (when \(K\rightarrow +\infty \)) to Fisher’s log-series, when \(K\rightarrow 0\). Moreover, these parameters have biological/ecological interpretations which are detailed in the literature and in this study. We compared three estimators of K, \(\mathfrak {P}\) and the parameter \(\alpha \) of Fisher’s log-series, following the work of Rao CR (Statistical ecology. Pennsylvania State University Press, University Park, 1971) on a three-parameter unstandardized variant of the negative binomial distribution. We further investigated the coherence underlying parameter values resulting from the different estimators, using both real count data collected in the Mauritanian Exclusive Economic Zone (MEEZ) during the period 1987–2010 and realistic simulations of these data. In the case of the MEEZ, we first built homogeneous lists of counts (replicates), by gathering observations of each species with respect to “typical environments” obtained by clustering the sampled stations. The best estimation of \(\left( K,\mathfrak {P}\right) \) was generally obtained by penalized minimum Hellinger distance estimation. Interestingly, the parameters of most of the correctly sampled species seem compatible with the classical birth-and-dead model of population growth with immigration by Kendall (Biometrika 35:6–15, 1948).     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Quantifying a significance of sediment particle size to hyporheic sedimentary oxygen demand with a permeable stream bed

A mechanistic model of sedimentary oxygen demand (SOD) for hyporheic flow is presented. The permeable sediment bed, e.g. sand or fine gravel, is considered with hydraulic conductivity in the range $0.1 < K < 20$  cm/s. Hyporheic pore water flow is induced by pressure fluctuations at the sediment/water interface due to near-bed turbulent coherent motions. A 2-D advection–diffusion equation is linked to the pore water flow model to simulate the effect of advection–dispersion driven by interstitial flow on oxygen transfer through the permeable sediment. Microbial oxygen uptake in the sediment is expressed as a function of the microbial growth rate, and is related to the sediment properties, i.e. the grain diameter $(d_{s})$ and porosity $(\phi )$ . The model describes the significance of sediment particle size to oxygen transfer through the sediment and microbial oxygen uptake: With increasing grain diameter $(d_{s})$ , the hydraulic conductivity $(K)$ increases so does the oxygen transfer rate, while particle surface area per volume (the available surface area for colonization by biofilms) decreases reducing the microbial oxygen uptake rate. Simulation results show that SOD increases as the hydraulic conductivity $(K)$ increases before a threshold has been reached. After that, SOD diminishes with the increment of the hydraulic conductivity $(K)$ .     

Numerical modelling of horizontal sediment-laden jets

Sediment-laden turbulent flows are commonly encountered in natural and engineered environments. It is well known that turbulence generates fluctuations to the particle motion, resulting in modulation of the particle settling velocity. A novel stochastic particle tracking model is developed to predict the particle settling out and deposition from a sediment-laden jet. Particle velocity fluctuations in the jet flow are modelled from a Lagrangian velocity autocorrelation function that incorporates the physical mechanism leading to a reduction of settling velocity. The model is first applied to study the settling velocity modulation in a homogeneous turbulence field. Consistent with basic experiments using grid-generated turbulence and computational fluid dynamics (CFD) calculations, the model predicts that the apparent settling velocity can be reduced by as much as 30 % of the stillwater settling velocity. Using analytical solution for the jet mean flow and semi-empirical RMS turbulent velocity fluctuation and dissipation rate profiles derived from CFD predictions, model predictions of the sediment deposition and cross-sectional concentration profiles of horizontal sediment-laden jets are in excellent agreement with data. Unlike CFD calculations of sediment fall out and deposition from a jet flow, the present method does not require any a priori adjustment of particle settling velocity.     

Lock-exchange gravity currents over rough bottoms

In nature, density driven currents often flow over or within a bottom roughness: a sea breeze encountering tall buildings, a shallow flow encountering aquatic vegetation, or a dense oceanic current flowing over a rough bottom. Laboratory experiments investigating the mechanisms by which bottom roughness enhances or inhibits entrainment and dilution in a lock-exchange dense gravity current have been conducted. The bottom roughness has been idealized by an array of vertical, rigid cylinders. Both spacing (sparse vs. dense configuration) and height of the roughness elements compared with the height of the current have been varied. Two-dimensional density fields have been obtained. Experimental results suggest that enhancement of the entrainment/dilution of the current can occur due to two different mechanisms. For a sparse configuration, the dense current propagates between the cylinders and the entrainment is enhanced by the vortices generated in the wake of the cylindrical obstacles. For a dense configuration, the dense current rides on top of the cylinders and the dilution is enhanced by the onset of convective instability between the dense current above the cylinders and the ambient lighter water between the cylinders. For low values of the ratio of the cylinder to lock height \(\lambda \) the dense current behavior approaches that of a current over a smooth bottom, while the largest deviations from the smooth bottom case are observed for large values of \(\lambda \).     

Nanostructured semiconducting materials for efficient hydrogen generation

Massive production of hydrogen by water decomposition triggered by a solar light active photocatalyst is a major objective in chemistry and a promising avenue to overcome the global energy crisis. The development of efficient, stable, economically viable and eco-friendly photocatalysts for hydrogen production is a challenging task. This article reviews the use of nanocomposite in three combinations: metal oxide–metal oxide semiconductor, metal–metal oxide semiconductor and metal chalcogenide–metal oxide core–shell nanostructures. These core–shell structures occur in two forms: a simple form where the photocatalyst is either in the core or the shell or in a more complex system where the core–shell structure comprises a co-catalyst deposited on a semiconducting material. We discuss the design, synthesis and development of semiconductor-based nanocomposite photocatalysts for hydrogen production. The major points are the role of catalytic active sites, the chemical nature of sacrificial agents, the effect of light sources, the variable light intensity and the energy efficiency calculation. For TiO₂-based nanocomposites, the metal oxide or metal co-catalyst loading of 1.0–3.0 wt% was optimal. TiO₂ nanotube–CuO hybrid nanocomposites produce 1,14,000 µmol h^?1 \({\text{g}}^{ - 1}_{\text{cat}}\), whereas TiO₂/Au nanocomposites display 1,60,000 µmol h^?1 \({\text{g}}^{ - 1}_{\text{cat}}\). For core–shell catalysts, a shell thickness of 2–20 nm was found for the best activity, and its performance is as follows: (a) CdS–NiO system produces around 19,949 µmol h^?1 \({\text{g}}^{ - 1}_{\text{cat}}\) and (b) CuO–Cr₂O₃ as co-catalyst immobilized on TiO₂ system produces around 82,390 µmol h^?1 \({\text{g}}^{ - 1}_{\text{cat}}\).