Time-averaged Field Research Articles

Periodic Rayleigh streaming vortices and Eckart flow arising from traveling-wave-based diffractive acoustic fields.

Recent studies have demonstrated that periodic time-averaged acoustic fields can be produced from traveling surface acoustic waves (SAWs) in microfluidic devices. This is caused by diffractive effects arising from a spatially limited transducer. This permits the generation of acoustic patterns evocative of those produced from standing waves, but instead with the application of a traveling wave. While acoustic pressure fields in such systems have been investigated, acoustic streaming from diffractive fields has not. In this work we examine this phenomenon and demonstrate the appearance of geometry-dependent acoustic vortices, and demonstrate that periodic, identically rotating Rayleigh streaming vortices result from the imposition of a traveling SAW. This is also characterized by a channel-spanning flow that bridges between adjacent vortices along the channel top and bottom. We find that the channel dimensions determine the types of streaming that develops; while Eckart streaming has been previously presumed to be a distinguishing feature of traveling-wave actuation, we show that Rayleigh streaming vortices also results. This has implications for microfluidic actuation, where traveling acoustic waves have applications in microscale mixing, separation, and patterning.

DEVELOPMENT OF A LIFETIME PRESSURE SENSITIVE PAINT PROCEDURE FOR HIGH-PRESSURE VANE TESTING

Abstract Optical measurements based on fast response Pressure Sensitive Paint (PSP) provide enhanced spatial resolution of the pressure field. This paper presents lifetime PSP at 20 kHz, with precise calibrations, and results from a demonstration in an annular vane cascade. The laser lifetime PSP methodology is first evaluated in a linear wind tunnel with a converging-diverging nozzle followed by a wavy surface. This test section is fully optically accessible with maximum modularity. A data reduction procedure is proposed for the PSP calibration, and optimal pixel binning is selected to reduce the uncertainty. In the annular test section, laser lifetime PSP was used to measure the time-averaged static pressure field on a section of the suction surface of a high-pressure turbine vane. Tests were performed at engine representative conditions in the Purdue Big Rig for Annular Stationary Turbine Analysis module at the Purdue Experimental Turbine Aerothermal Lab. The 2-D pressure results showed a gradual increase of pressure across the spanwise and flow directions, corroborated with local static pressure taps and computational results. The variation in PSP thickness was measured as a contribution to the uncertainty. The discrete Fourier transform of the unsteady pressure signal showed increased frequency content in wind-on compared to wind-off conditions at the mid-span and 30% span. Compared to the mid-span, the hub end wall region had an increase in frequencies and pressure amplitude. This result was anticipated given the expected presence of secondary flow structures at the near hub region.

Paleosecular Variation and the Time‐Averaged Geomagnetic Field Since 10 Ma

AbstractInvestigations into long‐term geomagnetic variations provide useful information regarding paleomagnetic field behavior. In this study, we assess the latitudinal structure of paleosecular variation (PSV) and the time‐averaged field (TAF) for the Brunhes normal and Matuyama reverse chrons, and for the 0–10 Ma period, from an updated and reviewed paleodirectional database spanning the past 10 Myr. The new database comprises 2,543 paleomagnetic sites from igneous rocks, providing improvements in the geographic and temporal distributions of high‐quality data relative to previous compilations. In addition, the new data collection differs considerably in application of strict selection criteria. Statistical analysis of the virtual geomagnetic pole (VGP) dispersion curve of Model G reveals a low latitudinal dependence of PSV for the last 10 Myr. For this period, we present a zonal TAF model based on the latitudinal distribution of inclination anomaly data. The best estimates found for axial quadrupole and octupole components were about 3% and 1% relative to axial dipole component, respectively. The new statistical models for the Brunhes and Matuyama chrons have different patterns in both PSV and TAF, in compliance with earlier studies. Our quantitative assessments indicate an apparent hemispheric PSV asymmetry, particularly in the Brunhes chron, with a stronger latitudinal signature in the southern hemisphere compared to the north. These findings suggest that equatorial PSV asymmetry, that has previously been found in modern, historical and millennial scale geomagnetic models, has persisted over the past 0.78 Ma.

Electromagnetic energy–momentum tensors in general dispersive bianisotropic media

The macroscopic electromagnetic (EM) energy–momentum tensor is one of the most important quantities characterizing the propagation and interaction of light in materials. In recent years, while exotic optical effects in various kinds of bianisotropic materials have been discovered, there still lacks a rigorous analysis of the energy and momentum of EM fields in such general cases. In this paper, using Noether’s theorem and the “Abrahamization” procedure, we obtain generalized Minkowski and Abraham EM energy–momentum tensors, applicable for both arbitrary time-dependent real EM fields and complex-valued analytic signals, in generic lossless bianisotropic media with frequency dispersion. The frequency dispersion of the materials modifies the expressions of EM energy density and Minkowski momentum, making them different from their familiar forms in nondispersive media. Our results reveal that the generalized Minkowski momenta for both real fields and analytic signals are conserved in source-free homogeneous media, while the Abraham momenta, characterizing the centroid motion of light, can change over time, which leads to the counterintuitive phenomenon that wave packets can travel along curved trajectories even in homogeneous bianisotropic media. We also show that the energy–momentum tensor for analytic signals derived from the action principle directly gives the conservation law of time-averaged fields and hence can describe the envelope evolution of waves in quasi-monochromatic approximation.

Stereoscopic particle image velocimetry in inhomogeneous refractive index fields of combustion flows.

Particle image velocimetry (PIV) measurements in reactive flows are disturbed by inhomogeneous refractive index fields, which cause measurement deviations in particle positions due to light refraction. The resulting measurement errors are known for standard PIV, but the measurement errors for stereoscopic PIV are still unknown. Therefore, for comparison, the velocity errors for standard and stereoscopic PIV are analyzed in premixed propane flames with different Reynolds numbers. For this purpose, ray-tracing simulations based on the time-averaged inhomogeneous refractive index fields of the studied non-swirled flame flows measured by the background-oriented Schlieren technique are performed to quantify the resulting position errors of the particles. In addition, the performance of the volumetric self-calibration relevant to tomographic PIV is analyzed with respect to the remaining position errors of the particles within the flames. The position errors cause significant standard PIV errors of 2% for the velocity component radial to the burner symmetry axis. Stereoscopic PIV measurements result in measurement errors of up to 3% radial to the burner axis and 13% for the velocity component perpendicular to the measurement plane. Due to the lower refractive index gradients in the axial direction, no significant velocity errors are observed for the axial velocity component. For the investigated flame configurations, the position errors and velocity errors increase with the Reynolds numbers. However, this dependence needs to be verified for other flame configurations such as swirled flame flows.

Eruptive and magma feeding system evolution of Sośnica Hill Volcano (Lower Silesia, SW Poland) revealed from Volcanological, Geophysical, and Rock Magnetic Data

The Sośnica Hill volcano is part of the Oligocene to Miocene (30.9–20.0 Ma) Strzelin volcanic field. It is located 100 km east of the Ohře Rift in the eastern part of the Fore-Sudetic Block, south of the town of Strzelin, Poland. Modern quarrying has exposed the sub-volcanic magma feeder system of the central part of the volcano and an extrusive volcanic succession that includes a 40 m thick sequence of lava flows and pyroclastic deposits that collectively form a complex scoria cone. Geophysical data (ground magnetometry and electric resistivity tomography (ERT)) reveal sharp linear anomalies that are interpreted to reflect normal faults dissecting the volcano. The ERT data map both high and low resistivity bodies, likely representing coherent clay-free dry rocks and partly argilized volcaniclastic deposits, respectively. Paleomagnetic data from 20 intrusive sites reveal two populations of reverse polarity site mean data; 11 sites are of low dispersion and yield a group mean direction that is discordant to the expected field direction, while six sites are highly scattered. Three sites did not yield interpretable results. We interpret the 11 sites as time-averaged field directions that are discordant to the expected field. The high dispersion of the remaining six sites are interpreted to indicate sub-volcanic deformation associated with the growth of the volcanic construct or multiple magma pulses over an extended period of time relative to secular variation. AMS data from 35 sites show a range of flow directions that vary across the quarry without an orderly pattern of fabric orientations. The flow pattern identified from dike paired margin data exhibits sub-vertical upward flow, sub-vertical downward, and moderately inclined northwest flow. Field observations and mapping indicate a complex development of the system in terms of styles of eruptive activity and structure of the final volcanic edifice. The activity included Strombolian and effusive phases, followed by phreatomagmatic, Hawaiian and again effusive eruptions. Such diversity of eruptive styles shows that the origin of the volcano is more complex than a simple, ‘textbook’ monogenetic scoria cone. Palaesoil on top of Strombolian deposits document a longer break in activity, after which eruptions resumed with different style; this is also not typical of monogenetic cones. The lateral variation in the volcanic succession suggests eruptions from several smaller, local vents. The complex subvolcanic magma flow patterns recorded in dikes match the variation of surface eruptive products and documents dynamically changing magma distribution paths in the near-surface and intra-cone part of the feeding system of the volcano.

Flow characteristics analysis for flow past the porous spheres: wake structures and drag force coefficients

Flow past a permeable sphere is different from that of a solid sphere due to the penetration of the fluid within porous structures, which can arise a change of flow fields. In this work, flow past porous spheres with Darcy numbers (Da) ranging from [0,10−3 ] were measured using planar Particle Image Velocimetry (PIV). The whole flow fields, including both leading edge and trailing edge, were captured at six different Reynolds numbers (Re) varying from 400 to 1400. Time-average flow fields were calculated based on instantaneous flow fields within fully-developed stages. Local minimum method was used to search for stagnation point positions. The results show positions of stagnation points are nearly proportional to the logarithm of Re. For most porous spheres, positions of stagnation points are extended to farther downstream positions than that of a solid sphere. However, at some certain Darcy numbers, ranging from 5 ∗ 10−6 to 2 ∗ 10−5, positions of stagnation points are closer to the sphere centers than that of an impermeable one.

Liquid hydrodynamics in a gas-liquid vortex reactor

A gas–liquid vortex reactor (GLVR) is considered one of the promising Process-Intensified Equipment for applications that require high mixing and mass transfer efficiency. To gain understanding of the fundamentals of GLVR operation, the gas–liquid hydrodynamics in the GLVR has been investigated by high-speed imaging, digital images analysis and cross correlation. Gas bubbles and liquid ligaments are observed in the GLVR over a wide operating window. The dynamic liquid layer thickness in the GLVR chamber varies in the range of 26–36 mm. The instantaneous and time-average velocity fields of the liquid phase are obtained with a guaranteed high signal-to-noise ratio (cross correlation peak ratio higher than 10), allowing the analysis of liquid turbulent properties including the turbulent kinetic energy and its dissipation rate. Results show that the high turbulent kinetic energy dissipation rate ranging from 23 to 114 m2/s3 contributes to a high mixing and mass transfer efficiency in the GLVR.

Visualization of the flow morphology of submerged jets in subcritical and supercritical water

Direct observation of flow behaviors is helpful for design of supercritical water reactor, so an in situ observation platform was built. Water jets under sub- and supercritical conditions (300–550 °C, 23 MPa) were observed by high-speed camera, which can be divided into four typical morphologies: mushroom momentum jets, amorphous plumes, banded plumes and cylindrical momentum jets. The time-average velocity field of each morphology was obtained through the image correlation method to quantitatively describe the flow and diffusion characteristics of jets. The influencing factors on the jet morphology, including the initial flow state of the jet, the driving force of the jet and the turbulence intensity of the ambient fluid, are analysed. Several dimensionless quantities are introduced to describe these factors quantitatively. The results show that the driving force of the jet is the decisive factor to appear to be a plume or momentum jet. The jet initial state has a large influence on the formation of the mushroom vortex. Turbulence of ambient fluid can lead to irregular deformation of the jet shape.

Towards the Large-Eddy Simulation of a full engine: Integration of a 360 azimuthal degrees fan, compressor and combustion chamber. Part II: Comparison against stand-alone simulations

Unsteady simulations of various components of a gas-turbine engine are often carried out independently and only share averaged quantities at the component interfaces. In order to study the impact and interactions between components, this work compares results from sectoral stand-alone simulations of a fan, compressor and annular combustion chamber of the DGEN-380 demonstrator engine at take-off conditions against an integrated 360 azimuthal degrees large-eddy simulation with over 2.1 billion cells of all previously listed components. Note that, at take-off conditions the compressor works at transonic conditions and generates an upstream-propagating shock that interacts with the fan modifying the shape of its wake with respect to the stand-alone simulation. Furthermore, the shock is seen as a tone in the pressure spectra at half the impeller blade passing frequency in the forward region of the engine. In the aft region, time-averaged fields are overall similar between stand-alone and integrated simulations but show a deviation in the azimuthal position of the hot-spot at the exit of the combustion chamber due to the addition of the diffuser. Pressure fluctuations generated in the compressor are captured in the combustion chamber as tones in the temperature and pressure spectra at the impeller blade-passing frequency and harmonics as well as an increase in the root-mean-square pressure.

Theory of pressure acoustics with thermoviscous boundary layers and streaming in elastic cavities.

We present an effective thermoviscous theory of acoustofluidics including pressure acoustics, thermoviscous boundary layers, and streaming for fluids embedded in elastic cavities. By including thermal fields, we thus extend the effective viscous theory by Bach and Bruus [J. Acoust. Soc. Am. 144, 766 (2018)]. The acoustic temperature field and the thermoviscous boundary layers are incorporated analytically as effective boundary conditions and time-averaged body forces on the thermoacoustic bulk fields. Because it avoids resolving the thin boundary layers, the effective model allows for numerical simulation of both thermoviscous acoustic and time-averaged fields in three-dimensional models of acoustofluidic systems. We show how the acoustic streaming depends strongly on steady and oscillating thermal fields through the temperature dependency of the material parameters, in particular the viscosity and the compressibility, affecting both the boundary conditions and spawning additional body forces in the bulk. We also show how even small steady temperature gradients ( ∼1 K/mm) induce gradients in compressibility and density that may result in very high streaming velocities ( ∼1 mm/s) for moderate acoustic energy densities ( ∼100 J/m3).

Characterization of temperature distributions in a swirled oxy-fuel coal combustor using tomographic absorption spectroscopy with fluctuation modelling

Oxy-fuel combustion promises efficient and inexpensive carbon dioxide sequestration and is therefore subject to active experimental research. However, the transfer of mature, non-intrusive diagnostic methods for temperature measurement like Coherent Anti-Stokes Raman Spectroscopy (CARS) from air-fed to oxy-fuel systems is challenging due to the deficiency in diatomic species suitable for thermometry. Although not limited to oxy-fuel atmospheres, we demonstrate the application of linear hyperspectral absorption tomography on water vapor in a swirled oxy-fuel coal combustor as a complementary diagnostic method, supplementing reference and validation data sets. Due to the burner design an axisymmetric reconstruction of the time-averaged temperature field is conducted. To compensate the temperature bias expected when evaluating time-averaged spectroscopic data we incorporate turbulent temperature fluctuations into our spectroscopic model, providing a fluctuation measure in addition to mean temperatures. The results quantitatively agree with vibrational O2-CARS measurements and qualitatively recreate spatial structures known from particle image velocimetry (PIV) flow fields.

Experimental and Numerical Investigation of Flow and Heat Transfer in Stationary Two-Pass Rectangular Duct (AR = 1:2) With Continuous and Broken V-Shaped Ribs

Abstract The combined experimental and large eddy simulations (LES) were performed in the stationary two-pass duct of aspect ratio (AR) 1:2. The experiments were conducted with three different rib arrangements, namely, 60 deg V, 60 deg V–IV, and broken 60 deg V–IV ribs, and the analysis was carried out with Reynolds numbers of 45,000, 60,000, and 75,000. The infrared thermography (IRT) technique is employed to obtain the local temperature distribution on heated smooth and ribbed surfaces. In all ribbed cases, the copper ribs are glued to the heated surface with a fixed rib height-to-hydraulic diameter (e/Dh) ratio of 0.125 and the rib pitch-to-height ratio (P/e) of 10 and 5 for continuous and broken ribs, respectively. In addition, the LES turbulence model was adopted for carrying out a simulation to understand the flow and heat transfer behavior in ducts populated with all three V-shaped ribs. The comparison of the time-averaged thermal fields generated using computations has been made with experimentally measured Nusselt numbers, friction factors, thermal performance factor (TPF), and Reynolds analogy performance parameter (RAPP) for all cases. The overall thermal performance factor was found to be quantitatively within 8.0–10.66% between experimental and numerical results. Among all the cases, the 60 deg V–IV ribbed duct provides the best TPF and RAPP than the other two ribbed ducts, whereas the smooth duct shows poor TPF.

Isothermal swirling flow characteristics and pressure drop analysis of a novel double swirl burner

Nonreacting swirling flow behavior through a novel double swirl burner has been investigated in this paper using Detached Eddy Simulation (DES) calculations and experimental measurements. The effect of mixing both the swirling flow of the inner jet (primary air or inner swirler) and the swirling flow of the annular jet (secondary air or outer swirler) on the flow structure has been analyzed. Moreover, pressure losses and flow field characteristics have been investigated at different inlet air mass flow rates and combustor diameters. Due to the drawbacks of general pressure loss to predict the pressure drop through high turbulence swirl flow as well as the complex geometry, the time-averaged viscous dissipation field function is used to predict the pressure drop through the burner. This function is based on the relationship between viscous dissipation and pressure losses. DES flow fields are compared with high-speed particle image velocimetry measurements at different inlet air mass flow rates. It is found that there is good agreement between numerical and experimental results, and both the central toroidal recirculation zone and the corner recirculation zone are well captured. In the case of a single inner swirler, the general pressure drop is very high compared to single outer and double swirl generators, and 84% of the inlet air pass through the outer swirler and just 16% of the air pass through the inner swirler in the case of a double burner (partially premixed). The time-averaged viscous dissipation field function is very important to identify the locations that lead to pressure drop inside the burner, which induces to select the appropriate design for the best performance. In addition, it is concluded that the performance of the double swirl burner is more efficient compared with the other two single swirlers in terms of pressure losses and flow structure.

Stochastic linear modes in a turbulent channel flow

Abstract

Persistent shallow magnetic inclination in the past 5 million years with implications for regional tectonics in the Philippines

Despite the growing knowledge on geomagnetic paleosecular variations, the region around the Philippine archipelago remains to be a huge data gap in this regard. This work looked into rock magnetic information from <5 Ma volcanic rocks of Luzon Island to provide data for the statistical characterization of the time-averaged field in the region.A total of 21 sites with ages from 1 ka to 4.3 Ma yielded 13 normal, 7 reversed, and 1 intermediate polarity dataset. These were aggregated with previously published work to form a combined dataset of 59 normal + reversed sites with a mean virtual geomagnetic pole at 321.1°E and 83.7°N with α95 = 3.4, and an angular standard deviation of about 14.6°. This amount of scatter is slightly larger than that predicted by Model G but relatively comparable with other near-equatorial regions such as Java and Ecuador. This work is a first examination of paleomagnetic data from the Philippines in relation to the long-term behavior of the magnetic field in this part of Southeast Asia. The long-ranging shallow inclination observed may also have implications on how paleomagnetic data are interpreted for tectonic studies in the region.

Effect of jet inflow angle on heat transportation in an impinging jet interfering with a Couette flow

The effect of the incoming angle of a jet that interferes with a Couette flow is studied using numerical simulations. When the jet flows against the direction of the Couette flow, the instantaneous vortex structure and heat transfer on the jet-impinging surface become complex and irregular. In contrast, when the jet flows in the same direction as the Couette flow, the instantaneous vortex structure has a strong and distinct periodicity about the main flow direction as a vortex street with a sine wave velocity fluctuation. This periodicity also affects the heat transfer on the jet-impinging surface. The Strouhal number calculated from these time periodicity increases with increasing jet angle. In the time-averaged flows and temperature fields, the Nusselt number exhibits a few characteristic distributions. The heat transfer concentrates around the jet impingement region at a negative jet angle whereas the area associated with significant heat transfer coefficient separates symmetrically about the flow in spanwise direction at a positive jet angle and the value of heat transfer coefficient on the center axis is low. Furthermore, a large and strong counter-rotating vortex pair is formed when the jet flows in the same direction as the Couette flow. As reported in our previous study, this pair vortex involves a main flow fluid and results in a characteristic Nusselt number distribution. Comparing the overall heat transfers on the jet-impinging surface, the lowest time-averaged heat transfer is obtained when the jet is orthogonal to the Couette main flow. Generally, the heat transfer over the surface is higher when the jet flows against the Couette flow than when the jet flows in the same direction as the Couette flow. However, since the pressure loss also increases drastically at the same time, it is found that the ratio of heat transfer over pressure loss increases with increase in jet angle. Considering the temperature and velocity vector distributions, and plots of turbulent heat flux, when the jet flows against the Couette main flow, the effect of jet impingement is dominant on heat transfer. However, when the jet flows in the same direction as the Couette main flow, most of the heat transfer is carried by the pair vortex. Although the contributing vortex type is different from similar previous studies conducted, the mechanism of generating this pair vortex is based on time-averaging the legs of the periodic vortex street.

Three-dimensional acoustic energy flow from rotating sound sources

A detailed study of the three-dimensional sound fields around a monopole point source in subsonic and supersonic rotation is investigated. We calculate the sound pressure, the acoustic velocity, as well as the instantaneous and time-averaged acoustic intensity fields with the advanced time approach. The advanced time approach is applied for sound sources at subsonic and supersonic rotation and yield comparable results as those obtained from the spherical harmonic series expansion method or the solution of the retarded time equation. Further, the direction of the time-averaged acoustic energy flow is derived from the acoustic intensity vectors. It is shown, that the direction of the energy flow path differs from the radial direction and depends on the rotational direction and rotational velocity. Those findings are useful, for example, to improve the acoustic absorption performance of sound absorbers and acoustic liners around turbomachines.

Large-Eddy Simulation of Turbulent Flow over Wavy Wall with Different Wave Steepness

Large-eddy simulation (LES) with dynamic one-equation subgrid-scale stress model was utilized to study the characteristics of turbulent flow over sea-surface generalized as two-dimensional wavy wall with different wave steepness. The statistical characteristics of turbulent field and pressure distribution were presented in detail. The simulation results showed that the separation zones induce stronger convection as the wave steepness increases. Reynolds shear stress near the wall boundary showed positive correlation with the wave steepness, while in the outer region, time-averaged turbulent fields seemed to be insensitive with the wall boundary. The equivalent velocity profiles obtained from the spatial averaging of the time-averaged velocity indicated that lower velocity occurred as the increase of wave steepness due to the increased pressure-induced form drag. The vortex structures were also visualized and showed vertically-bent characteristics which induced negative Reynolds shear stress at the stoss side as the increase of wave steepness that is equivalent to the effect of wave age.

Tomographic imaging using multi-simultaneous measurements (TIMes) for flame emission reconstructions.

The method of tomographic imaging using multi-simultaneous measurements (TIMes) for flame emission reconstructions is presented. Measurements of the peak natural CH* chemiluminescence in the flame and luminescence from different vaporised alkali metal salts that were seeded in a multi-annulus burner were used. An array of 29 CCD cameras around the Cambridge-Sandia burner was deployed, with 3 sets of cameras each measuring a different colour channel using bandpass optical filters. The three-dimensional instantaneous and time-averaged fields of the individual measured channels were reconstructed and superimposed for two new operating conditions, with differing cold flow Reynolds numbers. The contour of the reconstructed flame front followed the interface between the burnt side of the flame, where the alkali salt luminescence appears, and the cold gas region. The increased mixing between different reconstructed channels in the downstream direction that is promoted by the higher levels of turbulence in the larger Reynolds number case was clearly demonstrated. The TIMes method enabled combustion zones originating from different streams and the flame front to be distinguished and their overlap regions to be identified, in the entire volume.