The article presents an original spectroscopic method for determining the initial stage of the clustering process in supersonic jets. The technique has been tested on a supersonic argon jet excited by electrons at a distance of 30 millimetres from the nozzle outlet, where the clustering process is practically complete, and the influence of secondary processes on the intensity of the observed emissions is significantly suppressed. It is proposed to use the continuum emitted by neutral excited clusters with a wavelength of 127 nm as an indicator of cluster formation in a supersonic argon jet under various flow conditions. Analysis of the temperature dependence of the intensity of this continuum recorded at multiple argon pressures at the nozzle inlet allowed us to establish that the parameters of supersonic jet flow corresponding to the onset of crystallisation are related by the empirical expression. The calculated value of the constant for argon in the studied range of pressures and temperatures was 0.011.
This paper concerns the structural stability of smooth cylindrically symmetric supersonic Euler-Poisson flows in nozzles. Both three-dimensional and axisymmetric perturbations are considered. On one hand, we establish the existence and uniqueness of three-dimensional smooth supersonic solutions to the potential flow model of the steady Euler-Poisson system. On the other hand, the existence and uniqueness of smooth supersonic flows with nonzero vorticity to the steady axisymmetric Euler-Poisson system are proved. The problem is reduced to solve a nonlinear boundary value problem for a hyperbolic-elliptic mixed system. One of the key ingredients in the analysis of three-dimensional supersonic irrotational flows is the well-posedness theory for a linear second order hyperbolic-elliptic coupled system, which is achieved by using the multiplier method and the reflection technique to derive the energy estimates. For smooth axisymmetric supersonic flows with nonzero vorticity, the deformation-curl-Poisson decomposition is utilized to reformulate the steady axisymmetric Euler-Poisson system as a deformation-curl-Poisson system together with several transport equations, so that one can desi
Turbulence plays a critical role in the atmosphere, oceans, engineering, and astrophysics. The dissipation (heating) induced by turbulent flows is particularly important for the thermodynamics and chemistry of interstellar clouds, yet its structure and statistics remain poorly understood. Using high-resolution turbulence simulations with controlled explicit viscosity, we study the kinetic energy dissipation rate, $\varepsilon_{\mathrm{kin}}$, across subsonic and supersonic regimes. We find that dissipation lags large-scale kinetic energy injection events by $1.64\pm0.21$ and $0.48\pm0.07$ turbulent turnover times in subsonic and supersonic turbulence, respectively. Correlations show $\varepsilon_{\mathrm{kin}}\propto\vert abla\times\mathbf{v}\vert^2$ (vorticity squared) in the subsonic regime, where density fluctuations are negligible, while in the supersonic regime dissipation is primarily correlated with density, $\varepsilon_{\mathrm{kin}}\proptoρ^{3/2}$. A spectral analysis demonstrates that achieving numerical convergence of $\varepsilon_{\mathrm{kin}}$ across all scales is challenging, especially in the subsonic case, even at $2048^3$ resolution. Nonetheless, subsonic dissipa
This paper concerns the structural stability of smooth cylindrically symmetric supersonic spiral flows with large angular velocity for the steady Euler-Poisson system in a concentric cylinder. We establish the existence and uniqueness of some smooth supersonic Euler-Poisson flows with nonzero angular velocity and vorticity including both cylindrical spiral flows and axisymmetric spiral flows. The deformation-curl-Poisson decomposition for the steady Euler-Poisson system is utilized to deal with the hyperbolic-elliptic mixed structure in the supersonic region. For smooth cylindrical supersonic spiral flows, the key point lies on the well-posedness of a boundary value problem for a linear second order hyperbolic-elliptic coupled system, which is achieved by finding an appropriate multiplier to obtain the important basic energy estimates. The nonlinear structural stability is established by designing a two-layer iteration and combining the estimates for the hyperbolic-elliptic system and the transport equations. For smooth axisymmetric supersonic spiral flows, we use the special structure of the steady Euler-Poisson system to derive a priori estimates of the linearized second order el
A new method for determining the condensate fraction and cluster density in absolute units has been proposed and tested for a supersonic argon jet, which can also be applied to supersonic jets of other gases. The method is based on measuring the absolute intensities of the Ar II resonance lines (93.2 and 92.0 nm) when the supersonic jet is excited by an electron beam with constant current density. Knowing the absolute intensities of the lines and their emission cross sections, we determined the density of the noncondensed atomic component in the supersonic argon jet and the evolution of the condensate fraction over the whole temperature range investigated.
The increasing computational demand from growing data rates and complex machine learning (ML) algorithms in large-scale scientific experiments has driven the adoption of the Services for Optimized Network Inference on Coprocessors (SONIC) approach. SONIC accelerates ML inference by offloading it to local or remote coprocessors to optimize resource utilization. Leveraging its portability to different types of coprocessors, SONIC enhances data processing and model deployment efficiency for cutting-edge research in high energy physics (HEP) and multi-messenger astrophysics (MMA). We developed the SuperSONIC project, a scalable server infrastructure for SONIC, enabling the deployment of computationally intensive tasks to Kubernetes clusters equipped with graphics processing units (GPUs). Using NVIDIA Triton Inference Server, SuperSONIC decouples client workflows from server infrastructure, standardizing communication, optimizing throughput, load balancing, and monitoring. SuperSONIC has been successfully deployed for the CMS and ATLAS experiments at the CERN Large Hadron Collider (LHC), the IceCube Neutrino Observatory (IceCube), and the Laser Interferometer Gravitational-Wave Observat
This paper concerns supersonic flows with nonzero vorticity governed by the steady Euler-Poisson system, under the coupled effects of the electric potential and the geometry of a convergent nozzle. By the coordinate rotation, the existence of radially symmetric supersonic flows is proved. We then establish the structural stability of these background supersonic flows under multi-dimensional perturbations of the boundary conditions. One of the crucial ingredients of the analysis is the reformulation of the steady Euler-Poisson system into a deformation-curl-Poisson system and several transport equations via the deformation-curl-Poisson decomposition. Another one is to obtain the well-posedness of the boundary value problem for the associated linearized hyperbolic-elliptic coupled system, which is achieved through a delicate choice of multiplier to derive a priori estimates. The result indicates that the electric field force in compressible flows can counteract the geometric effects of the convergent nozzle and thereby stabilize key physical features of the flow.
Supersonic combustion plays a vital role in various applications, including scramjets, dual-mode ramjets, and pulse detonation engines. However, the flame characteristics can vary significantly, depending on the application. To model supersonic combustion, large eddy simulation coupled with a partially stirred reactor (PaSR) is commonly used. This method assumes that reactions occur at turbulent fine flame structure in a computational cell, making it less effective for flames governed by not only turbulent mixing but also homogenous autoignition, which is common under supersonic flow. To address this limitation, this study proposes a novel dynamic combustion model to enable more versatile modeling of supersonic flames. The model utilizes a two-delta probability density function, which represents the sub-grid composition variation and models the filtered sub-grid reaction rate. The corresponding weights are dynamically modeled based on the level of cell composition inhomogeneity. The new dynamic model is tested in a supersonic combustion case study with a strut-cavity flame holder. Results demonstrate that the new dynamic model can properly recover the limits of supersonic flames th
This paper concerns the structural stability of supersonic flows with a contact discontinuity in a finitely long curved nozzle for the two-dimensional steady compressible rotating Euler system. Concerning the effect of Coriolis force, we first establish the existence of supersonic shear flows with a contact discontinuity in the flat nozzle. Then we consider the stability of these background supersonic shear flows with a contact discontinuity when the incoming supersonic flow and the upper and lower nozzle walls are suitably perturbed. The problem can be formulated as an initial boundary value problem with a contact discontinuity as a free boundary. To deal with the free boundary value problem, the Lagrangian transformation is introduced to straighten and fix the contact discontinuity. The rotating Euler system is reduced to a first order hyperbolic system for the Riemann invariants. We design an iteration scheme and derive some estimates for the solution to the hyperbolic system. Finally, by using the inverse Lagrangian transformation, we prove the original free boundary problem admits two layers of smooth supersonic flows separated by a smooth contact discontinuity.
Software optimization refines programs for resource efficiency while preserving functionality. Traditionally, it is a process done by developers and compilers. This paper introduces a third option, automated optimization at the source code level. We present Supersonic, a neural approach targeting minor source code modifications for optimization. Using a seq2seq model, Supersonic is trained on C/C++ program pairs ($x_{t}$, $x_{t+1}$), where $x_{t+1}$ is an optimized version of $x_{t}$, and outputs a diff. Supersonic's performance is benchmarked against OpenAI's GPT-3.5-Turbo and GPT-4 on competitive programming tasks. The experiments show that Supersonic not only outperforms both models on the code optimization task but also minimizes the extent of the change with a model more than 600x smaller than GPT-3.5-Turbo and 3700x smaller than GPT-4.
Understanding the transition mechanism of boundary layer flows is of great significance in physics and engineering, especially due to the current development of supersonic and hypersonic aircraft. In this paper, we construct multiple unstable acoustic modes so-called Mack modes, which play a crucial role during the early stage of transition in the supersonic boundary layer. To this end, we develop an inner-outer gluing iteration to solve a hyperbolic-elliptic mixed type and singular system.
A novel supersonic jet oscillating method is investigated both experimentally and numerically. A rectangular primary supersonic jet is issued into a confined chamber with sudden enlargement. Secondary control jets are issued from the top, and bottom backwards-facing step regions formed due to sudden enlargement. The secondary jet also expands in the confined chamber shrouding the primary jet from the top and bottom sides. The primary jet is oscillated in the transverse direction by blowing the secondary jets in the streamwise direction in a pulsating manner with a phase shift. The out-of-phase secondary jet blowing causes the primary jet to periodically adhere to the upper and lower part of the confined chamber, causing flapping of the primary jet and acting as a supersonic fluidic oscillator. The supersonic jet oscillation characteristics are experimentally investigated using shadowgraph type flow visualization technique and steady and unsteady pressure measurements. Quantitative analysis of the shadowgraph images using the construction of $y-t$ and $y-f$ plots reveals the presence of periodic jet oscillation with a discrete dominant frequency similar to the secondary jet excitati
The chromosphere above active regions (ARs) on the Sun hosts magnetized supersonic downflows. Studies of these supersonic downflows help to decipher the magnetic fine structure and dynamics of the chromosphere. We perform a statistical analysis of the magnetized supersonic downflows in a number of ARs and survey their characteristics. We analyze spectro-polarimetric scans of parts of 13 ARs obtained in the infrared He I 10830 Å triplet formed in the upper chromosphere recorded with the GREGOR Infrared Spectrograph (GRIS) mounted at the GREGOR solar telescope. We retrieve the line-of-sight velocities and the magnetic field vector using the HeLIx+ inversion code that assumes Milne-Eddington atmospheres. We find magnetized supersonic downflows in all the ARs, with larger area coverage by such flows in ARs observed during their emerging phase. The fact that supersonic downflows were detected in all scans, though they cover only a small fraction, 0.2--6.4%, of the observed field-of-view, suggests that they are a common phenomenon in the upper chromospheres of ARs. The supersonic downflows are found to be associated with many AR features such as pores, sunspot umbrae, sunspot penumbrae,
The study performs large-eddy simulations of supersonic free jet flows using the Discontinuous Galerkin Spectral Element Method (DGSEM). The main objective of the present work is to assess the resolution requirements for adequate simulation of such flows with the DGSEM approach. The study looked at the influence of the mesh and the spatial discretization accuracy on the simulation results. The present analysis involves four simulations, incorporating three different numerical meshes and two different orders of spatial discretization accuracy. The numerical meshes are generated with distinct mesh topologies and refinement levels. Detailed descriptions of the grid generation and refinement procedures are presented. The study compares flow property profiles and power spectral densities of velocity components with experimental data. The results show a consistent improvement in the computed data as the simulation resolution increases. This investigation revealed a trade-off between mesh and polynomial refinement, striking a balance between computational cost and the accuracy of large-eddy simulation results for turbulent flow analyses.
We consider parallel splitting of a strip composed of two different chains. As a waveguide, the dissimilar-chain structure radically differs from the well-studied identical-chain system. It is characterized by three speeds of the long waves, two for the separate chains, and one for the intact strip where the chains are connected. Accordingly, there exist three ranges, the subsonic for both chains, the intersonic and the supersonic. The speed in the latter range is supersonic for the intact strip and at the same time, it is subsonic for the separate higher-speed chain. This fact allows the splitting wave to propagate in the strip supersonically. We derive steady-state analytical solutions and find that the splitting can propagate steadily only in two of these speed ranges, the subsonic and the supersonic, whereas the intersonic regime is forbidden. In the case of considerable difference in the chain stiffness, the lowest dynamic threshold corresponds to the supersonic regime. The peculiarity of the supersonic mode is that the supersonic energy delivery channel, being initially absent, is opening with the moving splitting point. Based on the discrete and related continuous models we
The resolvent formulation of McKeon & Sharma (2010) is applied to supersonic turbulent boundary layers to study the validity of Morkovin's hypothesis, which postulates that high-speed turbulence structures in zero pressure-gradient turbulent boundary layers remain largely the same as its incompressible counterpart. Supersonic zero-pressure-gradient turbulent boundary layers with adiabatic wall boundary conditions at Mach numbers ranging from 2 to 4 are considered. Resolvent analysis highlights two distinct regions of the supersonic turbulent boundary layer in the wave parameter space: the relatively supersonic region and the relatively subsonic region. In the relatively supersonic region, where the flow is supersonic relative to the freestream, resolvent modes display structures consistent with Mach wave radiation that are absent in the incompressible regime. In the relatively subsonic region, we show that the low-rank approximation of the resolvent operator is an effective approximation of the full system and that the response modes predicted by the model exhibit universal and geometrically self-similar behaviour via a transformation given by the semi-local scaling. Moreover,
Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposition from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions. Our study include jet--clump interaction, transient jets, and magnetised jets. We find that collimated supersonic jets do not excite supersonic motions far from the vicinity of the jet. Supersonic fluctuations are damped quickly and do not spread into the parent cloud. Instead subsonic, non-compressional modes occupy most of the excited volume. This is a generic feature which can not be fully circumvented by overdense jets or magnetic fields. Nevertheless, jets are able to leave strong imprints in their cloud structure and can disrupt dense clumps. Our results question the ability of collimated jets to sustain supersonic turbulence in molecular clouds.
We consider a supersonic flow past an airfoil in the context of the steady, isentropic and irrotational compressible Euler equations. We show that for appropriate data, either a shock forms, in which case certain derivatives of the solution blow up, or a sonic line forms, which translates to the fact that the hyperbolic character of the equations degenerates.
The presence of giant shells or loops in giant HII regions are clear witness of the mechanical energy input from massive stars. Here we evaluate the impact that winds may have on the structure of giant nebulae and on their supersonic velocity dispersion. We follow the suggestion from Chu \& Kennicutt (1994) to see if a combination of a large number of unresolved wind-driven shells caused by massive stars could produce the integrated broad Gaussian profiles typical of giant HII regions. The results, accounting for a wide range of energies, densities and velocity or age of the expanding shells, show that supersonic Gaussian profiles may arise only from a collection of unresolved wind-driven shells if the shells present a peculiar velocity distribution which implies a strongly peaked age distribution leading to an awkward star formation history. On the other hand, a uniform distribution of ages originates profiles with a flat-topped core defined by the terminal shell velocity and a steep decay as $v^{-6}$ up to the largest detectable shell speed. Thus, supersonic profiles can arise only if the final speed of the unresolved shells is supersonic. This implies an equally supersonic r
We develop a dynamical model of the propagating nonlinear localized excitations, supersonic kinks, in the cation layer in a silicate mica crystal. We start from purely electrostatic Coulomb interaction and add the Ziegler-Biersack-Littmark short-range repulsive potential and the periodic potential produced by other atoms of the lattice. This approach allows the construction of supersonic kinks which can propagate in the lattice within a large range of energies and velocities. The interparticle distances in the lattice kinks with high energy are physically reasonable values. The introduction of the periodic lattice potential results in the important feature that the kinks propagate with a single velocity and a single energy which are independent on the excitation conditions. The found kinks are ultra-discrete and can be described with the "magic wave number" $q\simeq 2π/3a$, which was previously revealed in the nonlinear sinusoidal waves and supersonic kinks in the Fermi-Pasta-Ulam lattice. The extreme discreteness of the supersonic kinks, with basically two particles moving at the same time, allows the interpretation of their double-kink structure. The energy of the supersonic kink