To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
This study is devoted to the analysis of capillary oscillations of a gas bubble in a liquid with an insoluble surfactant adsorbed on the surface. The influence of the Gibbs elasticity, the viscosities of the liquid and gas, as well as the shear and dilatational surface viscosities, on the damping of free oscillations is examined. Dependences of the frequency shift and the damping rate on the parameters of the problem are determined. In the limit of small viscosities and neglecting the surfactant surface diffusion, a simplified dispersion relation is obtained, which includes finite parameters of surface viscosities and Gibbs elasticity. From this relation, conditions are identified under which the damping of capillary oscillations can occur with a small frequency. Numerical solutions of the full dispersion relation demonstrate that a non-oscillatory regime is impossible for the considered configuration. An additional mode associated with Gibbs elasticity is discovered, characterized as a rule by low natural frequency and damping rate. Approximate relations for the complex natural frequency of bubble oscillations in a low-viscosity liquid in the presence of a surfactant are derived, including an estimate of the contribution of the gas inside the bubble to viscous dissipation. An original Lagrangian–Eulerian method is proposed and used to perform direct numerical simulations based on the full nonlinear Navier–Stokes equations and natural boundary conditions at the interface, accounting for shear and dilatational viscosities. The numerical data on the damping process confirm the results of the linear theory.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
As IBMs have gained popularity, their use has expanded to multiphysics problems in which the Navier-Stokes equations are only one among many other possibilities. In this chapter, a list of advanced applications is described in which IBMs are used to solve heat transfer, phase change and chemical reaction problems. These examples are intended as suggestions to extend the application of immersed boundary methods to complex physics problems.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
The various forcing strategies to be implemented in the governing equations are described in this chapter. Two big categories are first introduced, namely continuous forcing and discrete forcing methods. The various techniques are then detailed and the steps needed to implement them into an existing flow solver are described.
As any immersed boundary method has to be coupled with a solution algorithm for the governing equations, pseudo-compressibility and fractional-step methods are described in detail and some issues related to their combination with IBMs illustrated.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
With this chapter, the technical part of immersed boundary methods is initiated. Here it is explained how to define in the most convenient way a complex geometry object and how, after having immersed it in a computational grid, it is possible to determine the position (tagging) of the Eulerian nodes with respect to the boundary of the body.
Several computational geometry theorems are used to design an efficient computational algorithm which makes possible the tagging step within limited CPU time even when the computational grid contains tens of millions of nodes and the immersed object is described by hundreds of thousands of elements. This efficiency is key in problems involving moving bodies, deformable objects or fluid-structure interaction problems.
A comprehensive review on using different transmission lines for feeding ME dipole antennas and arrays is presented, including the SIW, ridge gap waveguide, packaged microstrip line, and substrate-integrated coaxial line feeds. In addition, the developments of low profile of ME dipole arrays, filtering ME dipoles, and all-metal ME dipole arrays for high-power applications are summarized. Some other recent applications are briefly reported. Hopefully, our readers can appreciate the attractiveness of the ME dipoles for future wireless applications at millimeter-wave and terahertz frequencies.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
This chapter is devoted to the application of IBMs to problems with moving boundaries. Specific adaptations of the algorithms are needed in order to cope with the Eulerian nodes at the interface that change position from inside to outside the body within one time step.
In turn, the boundary reconstruction of the solution is also affected and the necessary changes to the method are described.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
In this chapter it is explained how to compute the hydrodynamic loads produced by pressure and viscous stresses over an immersed surface. Several procedures are illustrated that entail different computational costs and degree of precision. The choice depends on whether only the resultant of the forces is needed or if the local values of the loads are needed. Finally, a simple validation of the discussed methods for a body with prescribed kinematics is shown.
Roberto Verzicco, Università degli Studi di Roma ‘Tor Vergata’, Gran Sasso Science Institute, L’Aquila, and University of Twente, Enschede,Marco D. de Tullio, Politecnico di Bari,Francesco Viola, Gran Sasso Science Institute, L’Aquila
Turbulent Rayleigh–Bénard convection in an extended layer of square cross-section with moderate aspect ratio $L/H=8.6$ ($L$ is the length of the cell, $H$ is its height) is studied numerically for Rayleigh numbers in the range ${\textit{Ra}}= 10^6{-}10^8$. We focus on the influence of different types of boundary conditions, including asymmetrical ones, on the characteristics of Rayleigh–Bénard convection with and without an immersed freely floating body. Convection without a floating body is characterised by the formation of stable thermal superstructures with preferred location. The crucial role of the symmetry of the boundary conditions is revealed. In the case of thermal boundary conditions of different types at the upper and lower boundaries, the flow pattern in Rayleigh–Bénard convection has a regular shape. The immersed body makes random wanderings and actively mixes the fluid, preventing the formation of superstructures. The mean flow structure with an immersed body is similar for all combinations of boundary conditions except for the case of a fixed heat flux at both boundaries. The floating disk does not change the tendency of turbulent convection to form a circulation of the maximal available scale under symmetric Neumann-type conditions. The type of boundary conditions has a weak influence on the Nusselt and Reynolds number values, significantly changing the ratio of the mean and fluctuating components of the heat flux. As the Rayleigh number increases, the motions of the body become more intensive and intermittent. The increase of $Ra$ also changes the structure of the mean flow without the body but the additional mixing provided by the floating body preserves the flow structure.
We conduct three-dimensional numerical simulations on centrifugal convection (CC) in a closed annular container, incorporating gravity and no-slip top and bottom boundaries, to systematically investigate rotation-induced secondary flow. The Stewartson layer, identified by an elongated circulation in mean vertical velocity plots, emerges near the inner and outer cylinders only beyond a critical gravitational forcing. Quantitative analyses confirm that the layer thickness scales as $\delta _{\,\!\textit{st}}\sim {\textit{Ek}}^{1/3}$ due to rotational effects, consistent with results from rotating Rayleigh–Bénard convection, where $Ek$ represents the Ekman number. The internal circulation strength, however, is determined by both gravitational buoyancy and rotational effects. We propose that gravitational buoyancy drives the internal flow, which balances against viscous forces to establish a terminal velocity. Through theoretical analysis, the vertical velocity amplitude follows $W_{\,\!\textit{st}}\sim {\textit{Ek}}^{5/3}\,Ro^{-1}\,{\textit{Ra}}_g\,Pr^{-1}$, showing good agreement with simulation results across a wide parameter range. Here, $Ro^{-1}$ represents the inverse Rossby number, ${\textit{Ra}}_g$ is the gravitational Rayleigh number, and ${\textit{Pr}}$ is the Prandtl number. The theoretical predictions match simulations well, demonstrating that the Stewartson layer is gravity-induced and rotationally constrained through geostrophic balance in the CC system. These findings yield fundamental insights into turbulent flow structures and heat transfer mechanisms in the CC system, offering both theoretical advances and practical engineering applications.
A ϕ2.5 m Gregorian antenna has been designed, analyzed, developed, and tested for a contoured beam for India and its Islands. The downlink band is 11.7–12.2 GHz. Uplink band is 17.3–17.8 GHz. The frequency band ratio is 1:1.52. The feed system consists of a radial corrugated horn and a turnstile orthomode transducer. Since there is a common feed system for Ku-Tx and Ku-Rx bands, it causes a large variation in the phase center for Tx and Rx bands. This large disparity in Tx and Rx phase center makes shaping challenging to achieve higher edge of coverage (EoC) gain, over widely separated transmit and receive bands. The optimization is carried out to achieve enhanced EoC gain and compliance of cross-polar isolation (XPI) through surface shaping of main and sub-reflectors. We have also optimized the feed coordinate to achieve the goal. The scattered near field at feed aperture is also minimized to get good XPI. Minimum radius of curvature of the surfaces has been controlled, which is required for the manufacturing of sub- and main-reflectors with minimum fabrication error. A new technique has been devised for the accommodation of a Gregorian antenna on a spacecraft with suitable radio frequency (RF) clearance. The impact of photogrammetry targets on RF performance is also brought out in the article.
In this paper, we showcase how flow obstruction by a deformable object can lead to symmetry breaking in curved domains subject to angular acceleration. Our analysis is motivated by the deflection of the cupula, a soft tissue located in the inner ear that is used to perceive rotational motion as part of the vestibular system. The cupula is understood to block the rotation-induced flow in a toroidal region with the flow-induced deformation of the cupula used by the brain to infer motion. By asymptotically solving the governing equations for this flow, we characterise regimes for which the sensory system is sensitive to either angular velocity or angular acceleration. Moreover, we show the fluid flow is not symmetric in the latter case. Finally, we extend our analysis of symmetry breaking to understand the formation of vortical flow in cavernous regions within channels. We discuss the implications of our results for the sensing of rotation by mammals.
Vertical thermal convection exhibits weak turbulence and spatio-temporally chaotic behaviour. For this configuration, we report seven new equilibria and 26 new periodic orbits. These orbits, together with four previously studied in Zheng et al. (J. Fluid Mech., 2024b, vol. 1000, p. A29) bring the number of periodic-orbit branches computed so far to 30, all solutions to the fully nonlinear three-dimensional Navier–Stokes equations. These new and unstable invariant solutions capture intricate spatio-temporal flow patterns including straight, oblique, wavy, skewed and distorted convection rolls, as well as bursts and defects. These interesting and important fluid mechanical processes in a small flow unit are shown to also appear locally and instantaneously in a chaotic simulation in a large domain. Most of the solution branches show rich spatial and/or spatio-temporal symmetries. The bifurcation-theoretic organisation of these solutions is discussed; the bifurcation scenarios include Hopf, pitchfork, saddle-node, period-doubling, period-halving, global homoclinic and heteroclinic bifurcations, as well as isolas. Furthermore, these orbits are shown to be able to reconstruct statistically the core part of the attractor, so that these results may contribute to a quantitative description of transitional fluid turbulence using periodic orbit theory.
BaLa2Cu1−xBaxTi2O9 (x = 0.00, 0.15, and 0.30) ceramics were synthesized in polycrystalline form via the conventional solid-state reaction techniques in air. The crystal structure of the title compositions was characterized by room-temperature X-ray powder diffraction and analyzed using the Rietveld refinement method. All the compositions crystallize in the tetragonal symmetry of space group I4/mcm (No. 140) with cell volumes: 249.43(1) Å3 for x = 0.00, 249.42(1) Å3 for x = 0.15, and 250.05(1) Å3 for x = 0.30. The tilt system of the MO6 octahedra (M = Cu(Ba2)/Ti) corresponds to the notation a0a0c−. The MO6 octahedra share the corners via oxygen atoms in 3D. Along the c-axis, the octahedra are connected by O(1) atoms of (0, 0, 1/4) positions; while in the ab-plane, they are linked by O(2) atoms of (x, x + 1/2, 0) positions. The bond angle of M–O2–M is 168.6(7)° for x = 0.00, 168.6(6)° for x = 0.15, and 166.8(6)° for x = 0.30, whereas the bond angle of M–O1–M is constrained to be 180° by space group I4/mcm.
Fluid mechanics, solid state diffusion and heat conduction are deeply interconnected through the mathematics and physical principles that define them. This concise and authoritative book reveals these connections, providing a detailed picture of their important applications in astrophysics, plasmas, energy systems, aeronautics, chemical engineering and materials science. This sophisticated and focused text offers an alternative to more expansive volumes on heat, mass and momentum transfer and is ideal for students and researchers working on fluid dynamics, mass transfer or phase transformations and industrial scientists seeking a rigorous understanding of chemical or materials processes. Accessible yet in depth, this modern treatment distills the essential theory and application of these closely related topics, includes numerous real world applications and can be used for teaching a range of related courses in physics, engineering and materials science departments.
In this paper, we study experimentally the dispersion of colloids in a two-dimensional, time-independent, Rayleigh–Bénard flow in the presence of salt gradients. Due to the additional scalar, the colloids do not follow exactly the Eulerian flow field, but have a (small) extra velocity $\boldsymbol{v}_{{dp}} = D_{{dp}}\, \boldsymbol{\nabla }\log C_s$, where $D_{{dp}}$ is the phoretic constant, and $C_s$ is the salt concentration. Such a configuration is motivated by the theoretical work by Volk et al. (2022, J.FluidMech., vol. 948, A42), which predicted enhanced transport or blockage in a stationary cellular flow depending on the value of a blockage coefficient. By means of high dynamical range light-induced fluorescence, we study the evolution of the colloids concentration field at large Péclet number. We find good agreement with the theoretical work, although a number of hypotheses are not satisfied, as the experiment is non-homogeneous in space, and intrinsically transient. In particular, we observe enhanced transport when salt and colloids are injected at both ends of the Rayleigh–Bénard chamber, and blockage when colloids and salt are injected together and phoretic effects are strong enough.
We present experiments of settling and dissolving sugar grains continuously sieved above a water tank with varying grain size and mass flux. Through drag and dissolution, grains force a downward flow whose dynamics are analysed in a laser sheet through particle image velocimetry and the use of home-made fluorescent sugar to track the negatively buoyant sugary water. We reveal different regimes, mostly controlled by the grain size, from a particle-constrained laminar flow at large grain size, to a turbulent plume with an effectively fluid-like behaviour when grains are small. The transitions between regimes are predicted from dimensionless numbers quantifying fluid–particle coupling, collective effects between grains and the possible onset of a Rayleigh–Taylor instability at the source. When a quasi-steady regime is reached, all grains dissolve above a finite depth, below which the flow is exclusively driven by dissolved sugar. We derive simple idealised models based on the source properties that predict the depth of this dissolution layer as well as the characteristic flow velocity.