No Arabic abstract
Miniature heaters are immersed in flows of quantum fluid and the efficiency of heat transfer is monitored versus velocity, superfluid fraction and time. The fluid is $^4$He helium with a superfluid fraction varied from 71% down to 0% and an imposed velocity up to 3 m/s, while the characteristic sizes of heaters range from 1.3 $mu$m up to few hundreds of microns. At low heat fluxes, no velocity dependence is observed. In contrast, some velocity dependence emerges at larger heat flux, as reported previously, and three non-trivial properties of heat transfer are identified. First, at the largest superfluid fraction (71%), a new heat transfer regime appears at non-null velocities and it is typically 10% less conductive than at zero velocity. Second, the velocity dependence of the mean heat transfer is compatible with the square-root dependence observed in classical fluids. Surprisingly, the prefactor to this dependence is maximum for an intermediate superfluid fraction or temperature (around 2 K). Third, the heat transfer time series exhibit highly conductive short-lived events. These textit{cooling glitches} have a velocity-dependent characteristic time, which manifest itself as a broad and energetic peak in the spectrum of heat transfer time series, in the kHz range. After showing that the velocity dependence can be attributed to the breaking of superfluidity within a thin shell surrounding heaters, an analytical model of forced heat transfer in a quantum flow is developed to account for the properties reported above. We argue that large scale flow patterns must form around the heater, having a size proportional to the heat flux (here two decades larger than the heater diameter) and resulting in a turbulent wake. The observed spectral peaking of heat transfer is quantitatively consistent with the formation of a Von Karman vortex street in the wake of a bluff body.
In this paper, the problem of compressible flow over a thin airfoil located near the ground is studied. A singular integral equation, also known as Possio equation, that relates the pressure jump along the airfoil to its downwash is derived. The derivation of the equation utilizes Laplace transform, Fourier transform, method of images, and theory of Mikhlin multipliers. The existence and uniqueness of solution to the Possio equation is verified for the steady state case and an approximate solution is obtained. The aerodynamic loads are then calculated based on the approximate solution. Moreover, the divergence speed of a continuum wing structure located near the ground is obtained based on the derived expressions for the aerodynamic loads.
The movement of subaqueous sediment in laminar shearing flow is numerically investigated by the coupled lattice Boltzmann and discrete element methods. First, the numerical method is validated by comparing the phase diagram proposed by Ouriemi {it et al.} ({it J. Fluid Mech}., vol. 636, 2009, pp. 321-336). Second, a detailed study on sediment movement is performed for sediment with varying solid volume fractions, and a nonlinear relationship between the normalised thickness of the mobile layer and the normalised fluid flow rate is observed for a densely-packed sediment. Third, an independent investigation on the effective viscosity and friction coefficient of the sediment under different fluid flow rates is conducted in a shear cell; and substitution of these two critical parameters into a theoretical expression proposed by Aussillous {it et al.} ({it J. Fluid Mech}., vol. 736, 2013, pp. 594-615) provides consistent predictions of bedload thickness with the simulation results of sediment movement. Therefore, we conclude that the non-Newtonian behaviour of densely-packed sediment leads to the nonlinear relationship between the normalised thickness of the mobile layer and the normalised fluid flow rate.
In this work, the static stability of plates with fixed trailing edges in axial airflow is studied using the framework of Possio integral equation. First, we introduce a new derivation of a Possio integral equation that relates the pressure jump along thin plates to their downwash based on the linearization of the governing equations of an ideal compressible fluid. The steady state solution to the Possio equation is used to account for the aerodynamic forces in the steady state plate governing equation resulting in a singular differential-integral equation which is transformed to an integral equation. Next, we verify the solvability of the integral equation based on the Fredholm alternative for compact operators in Banach spaces and the contraction mapping theorem. Then, we derive explicit formulas for the characteristic equations of free-clamped and free-pinned plates. The minimum solutions to the characteristic equations are the divergence speeds which indicate when static instabilities start to occur. We show analytically that free-pinned plates are statically unstable. After that, we move to derive analytically flow speed intervals that correspond to static stability regions for free-clamped plates. We also resort to numerical computations to obtain an explicit formula for the divergence speed of free-clamped plates. Finally, we apply the obtained results on piezoelectric plates and we show that free-clamped piezoelectric plates are statically more stable than conventional free-clamped plates due to the piezoelectric coupling.
A sensitive porosity adjoint method (SPAM) for optimizing the topology of fluid machines has been proposed. A sensitivity function with respect to the porosity has been developed. In the first step of the optimization process, porous media are introduced into the flow regime according to the sensitivity function. Then the optimized porous media are transformed to solid walls. The turbulent flow in porous media is accounted for by a modified eddy-viscosity based turbulence model. Its influence on the adjoint equations is nevertheless neglected, which refers to the so called frozen turbulence assumption. A test case of application in terms of the turbulent rough wall channel flow shows that a considerable reduction of the objective function can be obtained by this method. The transformation from porous media to solid walls may have important effect on the optimization results.
We present numerical simulations of laminar and turbulent channel flow of an elastoviscoplastic fluid. The non-Newtonian flow is simulated by solving the full incompressible Navier-Stokes equations coupled with the evolution equation for the elastoviscoplastic stress tensor. The laminar simulations are carried out for a wide range of Reynolds numbers, Bingham numbers and ratios of the fluid and total viscosity, while the turbulent flow simulations are performed at a fixed bulk Reynolds number equal to 2800 and weak elasticity. We show that in the laminar flow regime the friction factor increases monotonically with the Bingham number (yield stress) and decreases with the viscosity ratio, while in the turbulent regime the the friction factor is almost independent of the viscosity ratio and decreases with the Bingham number, until the flow eventually returns to a fully laminar condition for large enough yield stresses. Three main regimes are found in the turbulent case, depending on the Bingham number: for low values, the friction Reynolds number and the turbulent flow statistics only slightly differ from those of a Newtonian fluid; for intermediate values of the Bingham number, the fluctuations increase and the inertial equilibrium range is lost. Finally, for higher values the flow completely laminarises. These different behaviors are associated with a progressive increases of the volume where the fluid is not yielded, growing from the centerline towards the walls as the Bingham number increases. The unyielded region interacts with the near-wall structures, forming preferentially above the high speed streaks. In particular, the near-wall streaks and the associated quasi-streamwise vortices are strongly enhanced in an highly elastoviscoplastic fluid and the flow becomes more correlated in the streamwise direction.