Rheology of carbon black suspensions. V. Heat-induced gelation and its reversibility

Linear viscoelastic properties of carbon black (CB) suspensions in a mixture of a rosin-modified phenol resin-type varnish (Varnish-1)/an alkyd resin-type varnish (Varnish-2), which exhibited a sol?Cgel transition on an increase in CB concentration, were investigated from 30°C to 80°C. The viscoelastic properties were reversible from 30°C to 60°C. In contrast, at temperatures above 60°C, the storage (G??) and loss (G??) moduli were irreversible and increased significantly with increasing temperature. This increase in the moduli is due to a change of the dispersion state to agglomerated state by heating. The agglomerated state was held, when the suspensions were lowered at 30°C. However, the G?? and G?? recovered to the original values upon shearing. This heat-induced gelation should be a universal feature for suspensions of weakly attractive particles. The temperature and shearing histories of the suspensions were discussed in relation to adsorption of polymeric component in the varnish on the CB particles.     

Multiple laminar-turbulent transition cycles around a swept leading edge

Certain interesting flow features involving multiple transition/relaminarization cycles on the leading edge of a swept wing at low speeds are reported here. The wing geometry tested had a circular nose and a leading edge sweep of 60°. Tests were made at a chord Reynolds number of 1.3 × 10⁶ with model incidence α varied in the range of 3°?18° in discrete steps. Measurements made included wing chord-wise surface pressure distributions and wall shear stress fluctuations (using hot-film gages) within about 10 % of the chord in the leading edge zone. Results at α = 16° and 18° showed that several (often incomplete) transition cycles between laminar-like and turbulent-like flows occurred. These rather surprising results are attributable chiefly to the fact that the Launder acceleration parameter K (appropriately modified for swept wings) can exceed a critical range more than once along the contour of the airfoil in the leading edge region. Each such crossing results in a relaminarization followed by direct retransition to turbulence as K drops to sufficiently low values. It is further shown that the extent of each observed transition zone (of either type) is consistent with earlier data acquired in more detailed studies of direct transition and relaminarization. Swept leading edge boundary layers therefore pose strong challenges to numerical modelling.     

Direct numerical simulation of three-dimensional flow past a yawed circular cylinder of infinite length

Direct numerical simulation of flow past a stationary circular cylinder at yaw angles (α) in the range of 0–60° was conducted at Reynolds number of 1000. The three-dimensional (3-D) Navier–Stokes equations were solved using the Petrov–Galerkin finite element method. The transition of the flow from 2-D to 3-D was studied. The phenomena that were observed in flow visualization, such as the streamwise vortices, the vortex dislocation and the instability of the shear layer, were reproduced numerically. The effects of the yaw angle on wake structures, vortex shedding frequency and hydrodynamic forces of the cylinder were investigated. It was found that the Strouhal number at different yaw angles (α) follows the independence principle. The mean drag coefficient agrees well with the independence principle. It slightly increases with the increase of α and reaches a maximum value at α=60°, which is about 10% larger than that when α=0°. The root-mean-square (r.m.s.) values of the lift coefficient are noticeably dependent on α.     

Flow around four cylinders arranged in a square configuration

This paper presents an experimental study of the flow around four circular cylinders arranged in a square configuration. The Reynolds number was fixed at Re=8000, the pitch-to-diameter ratio between adjacent cylinders was varied from P/D=2 to 5 and the incidence angle was changed from α=0° (in-line square configuration) to 45° (diamond configuration) at an interval of 7.5°. The flow field was measured using digital Particle Image Velocimetry (PIV) to examine the vortex shedding characteristics of the cylinders, together with direct measurement of fluid dynamic forces (lift and drag) on each cylinder using a piezoelectric load cell. Depending on the pitch ratio, the flow could be broadly classified as shielding regime (P/D≤2), shear layer reattachment regime (2.5≤P/D≤3.5) and vortex impinging regime (P/D≥4). However, this classification is valid only in the case that the cylinder array is arranged nearly in-line with the free stream (α≈0°), because the flow is also sensitive to α. As α increases from 0° to 45°, each cylinder experiences a transition of vortex shedding pattern from a one-frequency mode to a two-frequency mode. The flow interference among the cylinders is complicated, which could be non-synchronous, quasi-periodic or synchronized with a definite phase relationship with other cylinders depending on the combined value of α and P/D. The change in vortex pattern is also reflected by some integral parameters of the flow such as force coefficients, power spectra and Strouhal numbers.     

Analysis of upstream,double-row,cylindrical holes on primary and secondary effects of endwall flow and film cooling

Flow features and film cooling performance of five configurations of double-row, cylindrical holes, upstream of an E3 vane, in a linear cascade are numerically investigated. This simulation is completed using a verified turbulence model at four blowing ratios (M = 0.5, 1.0, 1.5, 2.0). The first three configurations have two rows of cylindrical holes, each row with the same compound angle (β=-45°, 0° or 45°), while the other two have two rows with opposite compound angles (β=-45°, 45° and β=45°, -45°), which are also referred to as double-jet film cooling (DJFC) holes. The primary effects on the downstream endwall and the secondary effects on the nearby airfoil of the cooled passage are analyzed and discussed in detail. Results show that at low blowing ratios the movement of the coolant is denominated by the interaction between the jets and vortices resulting in similar film coverage on both the endwall and airfoil. The effect of vortices is reduced at high blowing ratios. It is also shown that the movement of the coolant is determined by the initial velocity direction, as well as the film cooling configuration.     

Non-symmetric bi-stable flow around the Ahmed body

The flow around the Ahmed body at varying Reynolds numbers under yawing conditions is investigated experimentally. The body geometry belongs to a regime subject to spanwise flow instability identified in symmetric flow by Cadot and co-workers (Grandemange et al., 2013b). Our experiments cover the two slant angles 25° and 35° and Reynolds numbers up to 2.784 × 10⁶. Special emphasis lies on the aerodynamics under side wind influence. For the 35° slant angle, forces and moments change significantly with the yawing angle in the range 10° ≤ |β| ≤ 15°. The lift and the pitching moment exhibit strong fluctuations due to bi-stable flow around a critical angle β of ±12.5°, where the pitching moment changes sign. Time series of the forces and moments are studied and explained by PIV measurements in the flow field near the rear of the body.     

Improved-Delayed-Detached-Eddy Simulation of cavity-induced transition in hypersonic boundary layer

Hypersonic flow transition from laminar to turbulent due to the surface irregularities, like local cavities, can greatly affect the surface heating and skin friction. In this work, the hypersonic flows over a three-dimensional rectangular cavity with length-to-width-to-depth ratio, L:W:D, of 19.9:3.57:1 at two angles of attack (AoA) were numerically studied with Improved-Delayed-Detached-Eddy Simulation (IDDES) method to highlight the mechanism of transition triggered by the cavity. The present approach was firstly applied to the transonic flow over M219 rectangular cavity. The results, including the fluctuating pressure and frequency, agreed with experiment well. In the hypersonic case at Mach number about 9.6 the cavity is seen as “open” at AoA of −10° but “closed” at AoA of −15° unconventional to the two-dimensional cavity case where the flow always exhibits closed cavity feature when the length-to-depth ratio L/D is larger than 14. For the open cavity flow, the shear layer is basically steady and the flow maintains laminar. For the closed cavity case, the external flow goes into the cavity and impinges on the bottom floor. High intensity streamwise vortices, impingement shock and exit shock are observed causing breakdown of these vortices triggering rapid flow transition.     

Testing the space-shuttle thermal-protection coating

A light-weight insulation material and its protective glassy coating will protect thespace shuttle from temperatures as high as 1250°C (2300°F). The critical performance characteristics of the brittle coating are investigated using testing techniques developed to accommodate these extreme environments and the delicate material. These include an ultimate-strain test-specimen geometry which circumvents problems created by flawed edges, as well as a tension specimen preparation and loading system with which premature failures due to excessive bending moment are avoided. Additionally, an elevated-temperature mechanical strain transducer—useable at more than 870°C (1600°F)—is described. Potential alterations to this sensor are discussed which would make it functional at up to 1600°C (3000°F).     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

Numerical investigation of the influence of incidence angle on asymmetrical turbine blade model showerhead film cooling effectiveness

The influence of various incidence angles on film cooling effectiveness of an axial turbine blade cascade with leading edge ejection from two rows of cooling holes is numerically investigated. The rows are located in the vicinity of the stagnation line. One row is located on the suction side and the other one is on the pressure side. The predicted pressure field for various blowing ratios (M = 0.7, 1.1 and 1.5) is compared with available experimental results at the design condition. Moreover, the effect of various incidence angles (?10°, ?5°, 0°, 5° and 10°) at three blowing rates is investigated by analyzing the results of both laterally averaged and area averaged values of adiabatic film cooling effectiveness. Numerical results indicate that the incidence angle can strongly affect the thermal protection of the blade at low blowing ratio but becomes less dominant at high blowing ratio. In fact, for the low blowing ratio, a small change in the incidence angle that relates to the design condition can deeply affect the thermal protection of the blade, which is evident from the laterally and area averaged film cooling effectiveness distributions.     

Statistical scale of hairpin packets in the later stage of bypass transition induced by cylinder wake

The hairpin packet's structure and its statistical scale in the later stage of bypass transition induced by a cylinder wake are investigated by time-resolved particle image velocimetry from the side and top view, respectively. Linear stochastic estimation is used to achieve the conditionally averaged velocity fields. For the side view case, the conditionally averaged structure consists of a series of swirling motions located along a line inclining at a large angle (18°) from the wall and a low-speed region occupied by the cylinder wake appearing right above them. In the (x, z)-plane at the wall-normal height y/???=?0.2, the dominant structures are shown to be the large-scale regions of low momentum elongated almost over 3?? along the streamwise. The low-speed regions are consistently bordered by small-scale counter-rotating vortice pairs organized in the streamwise with a statistical spanwise width of 0.55??. The results suggest that in the later part of the transitional zone, the upward induction of the cylinder wake enhances both the wall-normal and spanwise extent of the hairpin packets.     

Three-dimensional simulations of flow past two circular cylinders in side-by-side arrangements at right and oblique attacks

The flow past two identical circular cylinders in side-by-side arrangements at right and oblique attack angles is numerically investigated by solving the three-dimensional Navier–Stokes equations using the Petrov–Galerkin finite element method. The study is focused on the effect of flow attack angle and gap ratio between the two cylinders on the vortex shedding flow and the hydrodynamic forces of the cylinders. For an oblique flow attack angle, the Reynolds number based on the velocity component perpendicular to the cylinder span is defined as the normal Reynolds number Re_N and that based on the total velocity is defined as the total Reynolds number Re_T. Simulations are conducted for two Reynolds numbers of Re_N=500 and Re_T=500, two flow attack angles of α=0° and 45° and four gap ratios of G/D=0.5, 1, 3 and 5. The biased gap flow for G/D=0.5 and 1 and the flip-flopping bistable gap flow for G/D=1 are observed for both α=0° and 45°. For a constant normal Reynolds number of Re_N=500, the mean drag and lift coefficients at α=0° are very close to those at α=45°. The difference between the root mean square (RMS) lift coefficient at α=0° and that at α=45° is about 20% for large gap ratios of 3 and 5. From small gap ratios of 0.5 and 1, the RMS lift coefficients at α=0° and 45° are similar to each other. The present simulations show that the agreement in the force coefficients between the 0° and 45° flow attack angles for a constant normal Reynolds number is better than that for a constant total Reynolds number. This indicates that the normal Reynolds number should be used in the implementation of the independence principle (i.e., the independence of the force coefficients on the flow attack angle). The effect of Reynolds number on the bistable gap flow is investigated by simulating the flow for Re_N=100–600, α=0° and 45° and G/D=1. Flow for G/D=1 is found to be two-dimensional at Re_N=100 and weak three-dimensional at Re_N=200. While well defined biased flow can be identified for Re_N=300–600, the gap flow for Re_N=100 and 200 changes its biased direction too frequently to allow stable biased flow to develop.     

Study of the aerodynamic stability of a blunt cone with a side flap in a hypersonic flow

The stationary and time-dependent aerodynamic coefficients of a slender blunt cone with a flap located near the base section of the model are experimentally investigated. The freestream parameters (M_∞ = 6, Re _L = 0.88 × 10⁷, and γ = 1.4) ensured a turbulent regime of flow over the conical surface and the flap. At high angles of attack (α ～ 10°) laminar-turbulent transition is observable in the separation zone on the leeward side of the body. Emphasis is placed on the determination of the trimming angles of attack for different positions of the center of rotation and the static and dynamic stability coefficients (the model oscillation damping coefficient).     

Heat transfer performance comparison of steam and air in gas turbine cooling channels with different rib angles

Using steam as working fluid to replace compressed air is a promising cooling technology for internal cooling passages of blades and vanes. The local heat transfer characteristics and the thermal performance of steam flow in wide aspect ratio channels (W/H = 2) with different angled ribs on two opposite walls have been experimentally investigated in this paper. The averaged Nusselt number ratios and the friction factor ratios of steam and air in four ribbed channels were also measured under the same test conditions for comparison. The Reynolds number range is 6,000–70,000. The rib angles are 90°, 60°, 45°, and 30°, respectively. The rib height to hydraulic diameter ratio is 0.047. The pitch-to-rib height ratio is 10. The results show that the Nusselt number ratios of steam are 1.19–1.32 times greater than those of air over the range of Reynolds numbers studied. For wide aspect ratio channels using steam as the coolant, the 60° angled ribs has the best heat transfer performance and is recommended for cooling design.     

Flow characteristics and flow-induced forces of a stationary and rotating triangular cylinder with different incidence angles at low Reynolds numbers

In this paper, the problem of two-dimensional fluid flow past a stationary and rotationally oscillating equilateral triangular cylinder with a variable incident angle, Reynolds number, oscillating amplitude, and oscillating frequency is numerically investigated. The computations are carried out by using a two-step Taylor-characteristic-based Galerkin (TCBG) algorithm. For the stationary cases, simulations are conducted at various incident angles of α=0.0–60.0° and Reynolds numbers of Re=50–160. For the oscillation cases, the investigations are done at various oscillating amplitudes of θ_max=7.5–30.0° and oscillating frequencies of F_s/F_o=0.5–3.0 considering two different incidence angles (α=0.0°, 60.0°) and three different Reynolds numbers (Re=50, 100, 150). The results show that the influences of key parameters (incidence angle, Reynolds number, oscillating amplitude, and oscillating frequency) are significant on the flow pattern and hydrodynamic forces. For the stationary cases, at smaller angle of incidence (α≤30.0°), Reynolds number has a large impact on the position of the separation points. When α is between 30.0° and 60.0°, it was found that the separation points are located at the rear corners. From a topological point of view, the diagram of flow pattern is summarized, including two distinct patterns, namely, main separation and vortex merging. A deep analysis of the influence of Reynolds number and incidence angles on the mean pressure coefficient along the triangular cylinder surface is presented. Additionally, for the oscillating cases, the lock-on phenomenon is captured. The dominant flow patterns are 2S mode and P+S mode in lock-on region at α=0.0°. It is found at α=60.0°, however, that the flow pattern is predominantly 2S mode. Furthermore, except for the case of F_s/F_o=2.0, the mean drag decreases as the oscillating amplitude increases for each Reynolds number at α=0.0°. At α=60.0°, the minimum mean drag for F_s/F_o=1.5 is lower than that for stationary case, and occurs at θ_max=15.0° (Re=100) and θ_max=22.5° (Re=150), respectively. Finally, the effect of Reynolds number on a rotational oscillation cylinder is elucidated.     

Predicting crossflow induced transition with laminar kinetic energy transition model

The hypersonic laminar kinetic energy transition model is developed to be appropriate for crossflow induced boundary layer transition prediction. A crossflow timescale is constructed and incorporated in the k_T-k_L-ω transition model to reflect crossflow effect during three-dimensional boundary layer transition. The stream-wise vorticity is selected as the indicator of crossflow strength. Regarding the inviscid unstable characteristic of crossflow instability, the crossflow timescale is constructed by reference to the second mode timescale. To eliminate inappropriate development of the crossflow timescale where the effective length scale is large enough while the crossflow strength remains at a quite low level, a crossflow velocity limit function is proposed. The revised k_T-k_L-ω transition model has been applied to HIFiRE-5 and blunt cone with 1°angle of attack test cases. Results show good correspondence with the experimental data and DNS data, which demonstrates that the constructed crossflow timescale makes the revised transition model capable of reproducing crossflow induced transition behavior with a reasonable degree of accuracy.     

Wakes of elliptical cylinders at low Reynolds number

The wakes of elliptical cylinders are numerically investigated at a Reynolds number Re_D = 150. ANSYS-Fluent, based on the finite volume method, is used to simulate two-dimensional Newtonian fluid flow. The cylinder cross-sectional aspect ratio (AR) is varied from 0.25 to 1.0 (circular cylinder), and the angle of attack (α) of the cylinder is changed as α = 0° – 90°. With the changes in AR and α, three distinct wake patterns (patterns I, II, III) are observed, associated with different characteristics of fluid forces. Steady wake (pattern I) is characterised by two steady bubbles forming behind the cylinder, occurring at AR < 0.37 and α < 2.5°. Time-mean drag and fluctuating lift coefficients are small. Pattern II refers to Karman wake followed by steady wake (AR ≥ 0.37 – 0.67, depending on α) with the Karman street transitioning to two steady shear layers downstream. An inflection angle α_i is identified where the time-mean drag of the elliptical cylinder is identical to that of a circular cylinder. Pattern III is the Karman wake followed by secondary wake (AR ≤ 0.67, α > 52°), where the Karman street forming behind the cylinder is modified to a secondary vortex street with a low frequency. The Time-mean drag coefficient is maximum for this pattern.     

低Reynolds数NACA0012翼型绕流的流动特性分析

在水槽中应用PIV测速技术研究了NACA0012翼型在Reynolds数为8200时的流动特性,重点关注了翼型绕流结构中主频和扰动增长速率随迎角的变化。结果表明,分离剪切层的扰动增长符合指数规律;且随着迎角的增大,转捩过程加速,表现为扰动增长率逐渐增大,转捩的起始位置逐渐向上游移动。在所有实验迎角情况下,流场均由脱落旋涡主导,但其主导作用随着迎角的增大而削弱。     

Experimental investigation of the fluid–structure interaction in an elastic 180° curved vessel at laminar oscillating flow

Fluid–structure interaction phenomena are extremely important when laminar flows through elastic vessels such as in biomedical flow problems are considered. In general, such elastic vessels are curved which is why an elastic 180° bend at a curvature ratio \(\delta = D/D_{\rm C} = 0.\bar{2}\) defines the reference geometry in this study. It is the purpose of this study to compare the results with the steady flow through a 180° rigid pipe bend and to quantify the impact of the fluid–structure interaction on the overall flow pattern and the vessel deformation at oscillating fully developed entrance flow. The findings comprise velocity, pressure, and structure deformation measurements. The vessel dilatation amplitude was varied between 3.75 % and 7 % of the vessel diameter at Dean De and Womersley number Wo ranges of \(327\,\le\,De\,\le\,350\) and \(7\,\le\,Wo\,\le\,8.\) The flow is investigated by time-resolved stereoscopic particle-image velocimetry in five radial cross sections located in the elastic 180° bend and in the inlet pipes. The unsteady static vessel pressure is measured synchronously at these cross sections. The comparison of the steady with the unsteady flow field shows a strong change in the axial and secondary velocity distributions at periods of transition between the centrifugal forces and the unsteady inertia forces dominated regimes. These changes are characterized by asymmetric fluctuations of the centers of the counter-rotating vortex pair. The investigation of the impact of the structure deformation amplitude on these fluctuations reveals a significant attenuation at high deformation amplitudes. The spatial motion of the elastic vessel due to the forces applied by the flow exhibits amplitudes up to 15 % of the vessel diameter. Considering the fluid–structure interaction, an amplification of the volume flux amplitude by a factor of 2.1 at the vessel outlet and phase lags up to 30° occur. The static pressure distribution is characterized by a pronounced asymmetry between forward and backward flow with a 40 % higher peak magnitude at backward flow and phase lags of 35°. The results evidence that a strong distortion of the velocity distribution in the bend, which is caused by the oscillating nature of the flow, is reduced as a result of the fluid–structure interaction.     

Boundary-layer investigations on a model of the ELAC 1 configuration at high Reynolds numbers in the DNW

Results of boundary-layer investigations on the leeward side of a 1: 12 scale model of the ELAC 1 configuration of a space transportation system are presented. The configuration has the shape of a thick delta wing with a rounded leading edge. The model length is 6 m; the experiments were carried out in the 8×6 m ² low-speed German-Dutch-Windtunnel at Reynolds numbers up to Re=40·10⁶ . In a first series of experiments mean velocity profiles were determined in the turbulent boundary layer on the leeward side of the model, with a single hot-wire probe in the plane of symmetry at four positions. Comparison calculations with a numerical solution of the boundary-layer equations showed good agreement up to angles of attack α=10° . In a second series of tests the laminar-turbulent transition of the flow and its separation near the rounded leading edge were investigated at three positions with multi-sensor hot-film arrays with 40, 56, and 96 elements. These measurements demonstrated that the flow near the rounded leading edge is markedly influenced by the nose radius.