Axisymmetric Jet Research Articles

The effects of jet Lorentz factor on a relativistic astrophysical jet

In this Letter, the numerical simulation of axisymmetric hydrodynamic relativistic jet propagation was performed by solving the hydrodynamic relativistic Euler equation using the computer code PLUTO [Mignone et al., Astrophys. J. Suppl. Ser. 170, 228 (2007)]. The detailed flow features involved in this relativistic jet propagation has been thoroughly discussed in this present numerical study. The effect of the jet Lorentz factor (Γj) on the shock–turbulence interaction has been studied by analyzing the divergence of the Lamb vector (L=ω×U). The strong coexistence of two layers ∇·L<0 and ∇·L>0 enhances the momentum transfer due to energy difference, causing turbulence amplification.

Spatiotemporal structure of flow field by a dielectric-barrier-discharge plasma actuator using phase-resolved particle image velocimetry

A typical dielectric-barrier-discharge plasma actuator (DBD-PA) operating in continuous mode generates a self-similar starting vortex followed by a quasi-steady wall jet. In the burst mode, it generates periodic vortices, leading to a spatially wider mean flow field. The characteristics of periodic vortices can be controlled by adjusting the duty cycle α and burst frequency fb of the actuation signal. In the present study, the flow field generated by the actuator in burst mode is studied for 50%≥α≤90%, and 10 Hz ≥fb≤90 Hz. The Reynolds number of mean flow based on the maximum induced velocity and jet half-width ranges from 171 to 524. High-speed laser-sheet visualization and time-resolved particle image velocimetry are used for flow field measurements. The flow field generated during burst mode operation shows the strong influence of the burst frequency. Periodic vortices with comparable size and convection speed of starting vortex are generated at low fb. Dipoles and small-scale vortices are formed at moderate fb similar to the behavior of a transitional wall jet. Kelvin–Helmholtz instability is observed at higher burst frequency. An ornamental vortex consisting of several small periodic vortices is formed at high values of fb during the starting period. Based on the temporal evolution of momentum, the flow field generated by the DBD-PA can be classified into (i) starting flow regime, (ii) transitional flow regime, and (iii) periodic flow regime. The self-similarity of the mean flow field generated by the DBD-PA is verified by comparing the mean velocity profile at multiple downstream locations with the self-similar solutions of laminar and turbulent wall jets. The entrainment coefficient of the mean flow field is similar to that of an axisymmetric turbulent wall jet. The bi-orthogonal analysis shows that both coherent mode energy content and the entropy are minimum (about 5% only) for the continuous mode actuation. For the burst mode actuation, both the coherent mode energy content and the entropy are higher in magnitude compared to that of continuous mode and decrease with increasing burst frequency. The vortices in the periodic flow regime are locked in with fb. The distance between subsequent vortices, λx and flow frequency, f follow a power-law relationship, λx=cf−n, where c is a constant and n can be approximated as 2/3.

A Thread-Safe Lattice Boltzmann Model for Multicomponent Turbulent Jet Simulations

In this work an optimized multicomponent lattice Boltzmann (LB) model is deployed to simulate axisymmetric turbulent jets of a fluid evolving in a quiescent, immiscible environment over a wide range of dynamic regimes. The implementation of the multicomponent LB code achieves peak performance on graphic processing units (GPUs) with a significant reduction of the memory footprint, retains the algorithmic simplicity inherent to standard LB computing, and, being based on a high-order extension of the thread-safe LB algorithm, allows us to perform stable simulations at vanishingly low viscosities. The proposed approach opens attractive prospects for high-performance computing simulations of realistic turbulent flows with interfaces on GPU-based architectures.

Numerical studies of hydrodynamic relativistic astrophysical jet: The modification of vorticity transport equation

In this Letter, we have carried out the two-dimensional numerical simulation of axisymmetric relativistic jet in cylindrical coordinates by employing higher order finite volume method in PLUTO [Mignone et al., “PLUTO: A numerical code for computational astrophysics,” Astrophys. J. Suppl. Ser. 170, 228 (2007)] solver. The modified vorticity transport equation has been proposed for relativistic flow by taking the curl of the momentum equation, which shows significant change in the baroclinic vorticity production term due to relativistic effect. Both mathematical analysis and numerical results show that the vorticity production term due to baroclinic torque is heavily influenced due to the presence of specific enthalpy gradient and square of Lorentz factor gradient in a relativistic fluid flow.

Spectral proper orthogonal decomposition analysis of a spatially developing underexpanded round jet at Re = 45 000

This study conducts a comprehensive analysis of the coherent structures within a spatially developing underexpanded axisymmetric jet at a Reynolds number of 45 000, utilizing high-fidelity implicit large eddy simulations (iLES) in conjunction with spectral proper orthogonal decomposition (SPOD). In the frequency–wavenumber space, the global SPOD analysis identifies three distinct coherent structures, corresponding to three different mechanisms, namely, Kelvin–Helmholtz (KH), Orr, and lift-up. Their salient characteristics are discussed in detail. Local SPOD analysis further explores the streamwise evolution of these coherent structures, revealing that the influence of KH mechanism is confined to the near field, while the lift-up mechanism persists and dominates the energy content beyond the potential core, with the streaks of azimuthal wavenumbers one and two being the most energetic. The reconstruction of turbulent kinetic energy (TKE) and Reynolds shear stress from SPOD modes is assessed, and the first few azimuthal modes with low wavenumbers and frequencies are found crucial for capturing the dominant features of the flow. It is found that only the m = 0 and m = 1 modes contribute to the TKE at the centerline. The Reynolds shear stress reconstruction quality is comparable to TKE, but with a negligible contribution from the m = 0 mode. The azimuthal mode m = 1 captures the slope of the actual Reynolds shear stress profile in the vicinity of centerline, while m = 2 and higher modes capture the peak location of the actual profile.

Scalar mixing and entrainment in an axisymmetric jet subjected to external turbulence

The present study aims to understand the process of turbulent entrainment into a jet, as affected by background turbulence, using scalar statistics. Planar laser-induced fluorescence was employed to capture the orthogonal cross sections of the jet at a fixed downstream station with varying background turbulence intensities and length scales. The conditional scalar profiles revealed that the thickness of the scalar turbulent/turbulent interface is greater than that of the traditional turbulent/non-turbulent interface, and the interfacial thickness is an increasing function of the background turbulence intensity. Although nibbling remains the primary entrainment mechanism in the far field, increased occurrence of concentration “holes” within the interfacial layer in the presence of ambient turbulence suggests a more significant role of large-scale engulfment in the turbulent/turbulent entrainment process (although still below 1% of the total mass flux). Enhanced contribution of the area of detached jet patches (i.e., “islands”) to that of the main jet is hypothesized to be evidence of intense detrainment events in the background turbulence. This can potentially contribute to a reduced net entrainment into the jet, which manifests as less negative values of scalar skewness within the jet core.

Multifrequency polarimetry of high-synchrotron peaked blazars probes the shape of their jets

Multifrequency polarimetry is emerging as a powerful probe of blazar jets, especially with the advent of the Imaging X-ray Polarimetry Explorer (IXPE) space observatory. We studied the polarization of high-synchrotron peaked (HSP) blazars, for which both optical and X-ray emission can be attributed to synchrotron radiation from a population of nonthermal electrons. We adopted an axisymmetric stationary force-free jet model in which the electromagnetic fields are determined by the jet shape. When the jet is nearly parabolic, the X-ray polarization degree is ΠX ∼ 15–50%, and the optical polarization degree is ΠO ∼ 5–25%. The polarization degree is strongly chromatic: ΠX/ΠO ∼ 2–9. This chromaticity is due to the softening of the electron distribution at high energies, and is much stronger than for a uniform magnetic field. The electric vector position angle (EVPA) is aligned with the projection of the jet axis on the plane of the sky. These results compare very well with multifrequency polarimetric observations of HSP blazars. When the jet is instead nearly cylindrical, the polarization degree is large and weakly chromatic (we find ΠX ∼ 70% and ΠO ∼ 60%, close to the expected values for a uniform magnetic field). The EVPA is perpendicular to the projection of the jet axis on the plane of the sky. A cylindrical geometry is therefore practically ruled out by current observations. The polarization degree and the EVPA may be less sensitive to the specific particle acceleration process (e.g., magnetic reconnection or shocks) than previously thought.

Eulerian approach for erosion induced by particle-laden impinging jets

Particle dispersion and erosion play an important role in various engineering applications, particularly in impinging jet flows. This study presents a new Eulerian method for characterizing these phenomena in axisymmetric, laminar, and turbulent jets, impinging and transiently eroding a flat wall. The Eulerian mass and momentum conservation equations are solved assuming a one-way coupling between the flow and the particles. The carrier flow boundary layer velocity profiles in the laminar case are validated with analytical solutions. The particle calculation involves two stages, incident and reflected particles, connected via a restitution model. The reflected particles’ results offer insights into secondary erosion areas but are not used for the main erosion calculations. Subsequently, erosion at the eroded surface is computed using an empirical erosion model, followed by the displacement of computational nodes using the linear-elastic, small-strain deformation equations, all of which are solved transiently. The solutions reveal the spatial particle concentrations and momentum for different Stokes numbers: smaller Stokes particles follow the carrier flow streamlines, medium Stokes particles deviate, forming a W-shaped erosion profile, while the largest Stokes particles create a U-shaped profile. This numerical model offers a computationally efficient framework for CFD-based transient erosion calculation due to its Eulerian implementation.

DMD-based spatiotemporal superresolution measurement of a supersonic jet using dual planar PIV and acoustic data

The present study applies a framework of the spatiotemporal superresolution measurement based on the total-least-squares dynamic mode decomposition, the Kalman filter and the Rauch-Tung-Striebel smoother to an axisymmetric underexpanded supersonic jet of a jet Mach number of 1.35. Dual planar particle image velocimetry was utilized, and paired velocity fields of the flow with a short time interval were obtained at a temporal resolution of 5000 Hz. High-frequency acoustic data of 200,000 Hz were simultaneously obtained. Then, the time-resolved velocity fields of the supersonic jet were reconstructed at a temporal resolution of 200,000 Hz. Also, time coefficients of dynamic modes in high temporal resolution were calculated. The correlation between time coefficients implies that the mixing promotion by screech tone causes the lift-up of the high-velocity fluid from the jet center and accelerates at the downstream side.

Flowfield and Noise Dynamics of Supersonic Rectangular Impinging Jets: Major versus Minor Axis Orientations

The current study explores the flowfield and noise characteristics of an ideally expanded supersonic (Mach 1.44) rectangular jet impinging on a flat surface. The existing literature is primarily concentrated on axisymmetric jets, known for their resonance dominance, pronounced unsteadiness, and acoustic signatures. In contrast, non-axisymmetric jets remain relatively less understood, particularly those impinging on a ground surface. By employing Schlieren imaging, high-frequency pressure measurements using high-bandwidth transducers, and particle image velocimetry (PIV), this research comprehensively examines the flow-acoustic phenomena. Schlieren imaging revealed distinct, coherent structures and strong acoustic waves, while pressure measurements at the impingement surface exhibited high-amplitude fluctuations, peaking at approximately 186 dB. Acoustic analysis identified multiple high-amplitude tones with unique directional characteristics, suggesting the potential for multiple simultaneous modes in rectangular jets. Furthermore, the PIV data elucidated differences in the jet shear layer and wall jet development attributed to the nozzle orientation. These findings contribute to a deeper understanding of non-axisymmetric jet behavior, offering insights relevant to fundamental flow physics and practical applications such as vertical takeoff and landing aircraft.

Particle Image Velocimetry And Scaling Of Mach 0.3 And 1.25 Axisymmetric Turbulent Jets

The influence of compressibility on the scaling of the Reynolds stress transport (RST) budgets is investigated for Mach 0.3 and Mach 1.25 free axisymmetric turbulent jets using particle image velocimetry (PIV). It is proposed that the RST equations can be shown in self-preserving form scaled by the local streamwise evolving spreading rate, b’. It is shown that b’ is a function of the local convective Mach number, M_c. Compressibility effects on the suppression of the Reynolds stresses and higher turbulence moments are detailed. From self-preservation principles, the budgets scaled by powers of b’. In this form, most terms in the budgets are scaled appropriately indicating self-preservation. Terms which are not properly scaled are isolated then which allows for insight on effects of compressibility on turbulence.

Move With The Flow To Track The Evolving Turbulent/Non-Turbulent Interface

This study investigates the turbulent/non-turbulent interface (TNTI) in self-similar turbulent axisymmetric jet flows, focusing on a novel approach named ’move with the flow’ where the image acquisition is space-based rather than time-based. Experiments were conducted at three different Reynolds numbers (9000, 12000, and 31000), utilizing particle image velocimetry (PIV) and laser-induced fluorescence (LIF) techniques. The core of this research was developing a traverse system specifically designed to follow the evolving flow structures at the TNTI synchronously. We succeeded in tracking large-scale events in TNTI from creation to dissolution.

Determining the local value of heat transfer coefficient (HTC) on a surface cooled by an axisymmetric water jet

This paper presents the experimental findings of cooling an EN 1.4828 steel plate with a single water jet, and the results of computing the local heat transfer coefficient (HTC) and heat flux (HF) on the plate surface cooled. The plate cool from an initial temperature of 906 °C. The plate temperature measure 2 mm under the surface cooled, using 48 sheathed thermocouples 0.5 mm in diameter. Three-dimensional inverse solutions yielded HTC and HF as a function of surface coordinates and time. The uncertainty of the measured values and the resulting uncertainty of the calculated heat transfer coefficient were determined. Temperature distribution predict by the model is compared to the data obtained from measurements. Results allow to formulate an equation of local HTC as a function of non-dimensional temperature and the non-dimensional ratio of the distance to the stagnation point to the nozzle diameter. According to the author’s knowledge, this work fills a gap in the literature by providing relationship describing local HTC as a function of surface temperature and distance from the stagnation point in the case of single water jet cooling with unsteady heat transfer.

Gamma-Ray Burst Pulses and Lateral Jet Motion

We propose that gamma-ray burst (GRB) pulses are produced when highly relativistic jets sweep across an observer’s line of sight. We hypothesize that axisymmetric jet profiles, coupled with special relativistic effects, produce the time-reversed properties of GRB pulses. Curvature resulting from rapid jet expansion is responsible for much of the observed pulse asymmetry and hard-to-soft evolution. The relative obliqueness with which the jet crosses the line of sight explains the known GRB pulse morphological types. We explore two scenarios: one in which a rigid/semirigid jet moves laterally and another in which a ballistic jet sprays material from a laterally moving nozzle. The ballistic jet model is favored based upon its consistency with standard emission mechanisms.

High-order thread-safe lattice Boltzmann model for high performance computing turbulent flow simulations

We present a highly optimized thread-safe lattice Boltzmann model in which the non-equilibrium part of the distribution function is locally reconstructed via recursivity of Hermite polynomials. Such a procedure allows the explicit incorporation of non-equilibrium moments of the distribution up to the order supported by the lattice. Thus, the proposed approach increases accuracy and stability at low viscosities without compromising performance and amenability to parallelization with respect to standard lattice Boltzmann models. The high-order thread-safe lattice Boltzmann is tested on two types of turbulent flows, namely, the turbulent channel flow at Reτ=180 and the axisymmetric turbulent jet at Re = 7000; it delivers results in excellent agreement with reference data [direct numerical simulations (DNS), theory, and experiments] and (a) achieves peak performance [∼5×1012 floating point operations (FLOP) per second and an arithmetic intensity of ∼7 FLOP/byte on a single graphic processing unit] by significantly reducing the memory footprint, (b) retains the algorithmic simplicity of standard lattice Boltzmann computing, and (c) allows to perform stable simulations at vanishingly low viscosities. Our findings open attractive prospects for high-performance simulations of realistic turbulent flows on GPU-based architectures. Such expectations are confirmed by excellent agreement among lattice Boltzmann, experimental, and DNS reference data.

An Investigation on Uncontrolled and Vortex-Generator Controlled Supersonic Jets

The present study is carried out with a motivation to investigate the axisymmetric supersonic jet both experimentally and computationally. An open jet facility was utilized to carry out the experiments, and the results were compared with computational simulations employing the K-omega SST turbulence model using ANSYS software. It is important to note that, the computational validation has been done incorporating the Rayleigh Pitot formula to match the centerline pressure for the uncontrolled jet, which has not been found in any other validation studies according to the authors’ understanding. Besides, the experimental study is extended with a focus on evaluating the impact of Vortex Generators (VGs) on Mach 1.6 supersonic jets. The aim was to enhance jet mixing, a critical factor for improving engine performance. Various nozzle geometry modifications were explored in the past, but VGs emerged as the most effective method for optimizing jet mixing efficiency. The investigation revealed a substantial decrement in the supersonic jet core length when VGs were introduced at the nozzle exit, especially under favorable pressure gradients. This reduction in the supersonic core emphasized the role of VGs in enhancing mixing efficiency. The study also confirmed that VGs significantly distort wave patterns within the supersonic core, crucial for improved jet mixing. This research signifies the importance of VGs in augmenting the mixing of Mach 1.6 jets, offering the potential for improved jet performance and reduced noise emissions in the aerospace industry.

Averaging theory for heat transfer in circular hydraulic jumps with a separation bubble

Analytical investigations of heat transfer during the vertical impingement of an unsubmerged axisymmetric liquid jet on a horizontal plate have been limited to the regions ahead of the jump. This limitation is due to the complex flow physics in the jump region arising from sudden changes in the flow field. This is addressed in here by extending the averaging theory (AT) introduced by Bohr et al. (Phys. Rev. Lett., vol. 79, issue 6, 1997, pp. 1038–1041) which was further developed by Watanabe et al. (J. Fluid Mech., vol. 480, 2003, pp. 233–265), to describe the heat transfer problem in circular hydraulic jumps including separation. The applicability of the resulting theory to determine the temperature field in the jump region is evaluated using the data available in the literature and also by means of fully resolved numerical solutions. Good agreement is observed for moderate Prandtl numbers. However, for sufficiently high Prandtl numbers, deviations become notable. The reasons for the deviations according to their relevance are (i) monotonically decreasing temperature profile inherent to the AT, whereas the fully resolved numerical solutions exhibit a local maximum in the temperature profile away from the plate; and (ii) inapplicability of the concept of dividing the flow field into a region affected and a region unaffected by heat transfer according to the thermal boundary layer thickness. This concept leads to the overestimation of the temperature close to the wall and to the existence of a threshold Prandtl number, for which the thermal boundary layer thickness does not meet the free surface anymore. Around this threshold Prandtl number, the temperature field shows a discontinuous behaviour.

Complete characterization of axisymmetric turbulent jet using background oriented schlieren and physics-informed neural network

Complete characterization of axisymmetric turbulent jet using background oriented schlieren and physics-informed neural network

Features of the dynamics of several elevated jets in a cross-flow

Using the thermoanemometric method, the effect of turbulent cross-flow on vertical low-speed elevated axisymmetric jets with a relative velocity of approximately 2.45 has been experimentally investigated. Processing and qualitative analysis of experimental data on velocity profiles for the case of two and four identical isothermal turbulent jets propagating in a uniform cross-flow are carried out. The analysis of velocity contours in different cross-sections showed that the interaction of the jets proceeds in a complex way; and in this interaction, in addition to the jets and the cross flow, the trace under the jets and tubes – the sources of the jets – are involved. It is obtained that the trajectory of a twin parallel jet shows less penetration than that of a single jet. At the same time, the trajectory of the tandem jet after the merging of the rows shows significant penetration into the transverse flow compared to the single one.

Prompt GRB polarization from non-axisymmetric jets

ABSTRACT Time-resolved linear polarization (Π) measurements of the prompt gamma-ray burst emission can reveal its dominant radiation mechanism. A widely considered mechanism is synchrotron radiation, for which linear polarization can be used to probe the jet’s magnetic-field structure, and in turn its composition. In axisymmetric jet models, the polarization angle (PA) can only change by 90°, as Π temporarily vanishes. However, some time-resolved measurements find a continuously changing PA, which requires the flow to be non-axisymmetric in at least one out of its emissivity, bulk Lorentz factor, or magnetic field. Here, we consider synchrotron emission in non-axisymmetric jets, from an ultrarelativistic thin shell, comprising multiple radially expanding mini-jets (MJs) or emissivity patches within the global jet, that yield a continuously changing PA. We explore a wide variety of possibilities with emission consisting of a single pulse or multiple overlapping pulses, presenting time-resolved and integrated polarization from different magnetic field configurations and jet angular structures. We find that emission from multiple incoherent MJs/patches reduces the net polarization due to partial cancellation in the Stokes plane. When these contain a large-scale ordered field in the plane transverse to the radial direction, Π always starts near maximal and then declines over the single pulse or shows multiple highly polarized peaks due to multiple pulses. Observing $\Pi \lesssim 40~{{\ \rm per\ cent}}$ (15 per cent) integrated over one (several) pulse(s) will instead favour a shock-produced small-scale field either ordered in the radial direction or tangled in the plane transverse to it.