Meso-scale simulations of energy localization at hotspots provide closure models for multiscale frameworks of shock-to-detonation transition (SDT). Validation of such meso-scale calculations is challenging as direct comparison with experiments is constrained both by limitations of data acquisition in the experiments (e.g., of temperature fields) and modeling over-simplifications in the simulations. To address the latter problem and bring modeling closer to experiments, we advance a high-fidelity meso-scale computational framework for interface-resolved reactive calculations of shock initiation in plastic-bonded explosives (PBXs). Accurate resolution of shock and interfacial dynamics is achieved through higher-order (5th-order WENO) schemes, and sharp interface treatments are implemented for physically accurate material-material interactions. Recently obtained atomistics-consistent material models are used for HMX, with the grid resolution taken down to atomistic scale (O(nm)). The crystal geometries are obtained directly from experiments via nano-CT imaging. The impacting flyer plate, energetic crystal and binder are tracked as distinct phases, and flyer-binder impact and separatio
We derive the shock strength area rule for a Noble-Abel stiffened gas (NASG) equation of a state required in Whitham's geometrical shock dynamics approach to determine shock wave dynamics in dense gases, liquids and solids. An exact formulation requiring the solution of an ordinary differential equation is provided. Closed form solutions of various levels of approximations are also obtained as an expansion in shock wave strength. The leading order approximation recovers the geometrical acoustic limit, while higher order approximations account for the medium's compressibility. The exact shock strength area relation and the various order approximations are illustrated for shocks in liquid water. The simple closed form of the first order solution predicts the shock strength area rule up to shock pressures of approximately twice the stiffening pressure in water, i.e., approximately 1-2 GPa.
Accurate aerodynamic and aerothermodynamic predictions are crucial for numerous hypersonic applications. This paper proposes a gas-kinetic scheme (GKS) coupled with a two-temperature kinetic model, which distinguishes between the translational-rotational and vibrational modes of temperature. Compared with one-temperature model and the translational-rotational multi-temperature model, the proposed model provides a more physically accurate simulation of real gas effects when vibrational energy modes of air are excited. On the other hand, it is computationally simpler than multi-temperature model with independent translational, rotational and vibrational modes. The scheme is implemented on both structured and unstructured grids. To further improve the robustness for strong shock and rarefaction waves, the discontinuity feedback factor is employed instead of traditional limiters. Numerical verifications are conducted on one-dimensional shock structure, two-dimensional (2D) hypersonic flow over a cylinder, 2D hypersonic flow over a wedge and 2D Edney Type IV shock/shock interaction. Compared with experimental data, the reference results from direct simulation Monte Carlo (DSMC) method a
Accurate prediction of aerothermal loads in hypersonic flows is critical yet challenging due to the coupling of Shock-Wave/Boundary-Layer Interactions (SBLI) and thermal non-equilibrium. This work presents the development of a three-dimensional two-temperature Gas-Kinetic Scheme (3D 2T-GKS) on unstructured meshes. The scheme resolves translational-rotational and vibrational energy modes within a unified kinetic framework. A key innovation is the integration of a Generalized Kinetic Boundary Condition (GKBC), which physically decouples the thermal accommodation of vibrational energy from the translational-rotational mode, thereby offering a more accurate model for gas-surface interactions. Additionally, a Discontinuity Feedback Factor (DFF) is employed to capture strong shock waves with reduced numerical dissipation compared to classical limiters. The method is rigorously validated against standard experimental benchmarks, including the sharp double-cone and hollow cylinder-flare configurations. Numerical results demonstrate that the proposed solver, augmented by the GKBC, accurately captures complex wave structures, separation topologies, and surface heat flux distributions. These
Supersonic vortex plane and axisymmetric flows of non-viscous non-heatconductive gas with arbitrary thermodynamic properties in the vicinity of a steady shock wave are studied. The differential equations describing the gas flow exterior to the discontinuity surface and the dynamic compatibility conditions at this discontinuity are used. The gas flow nonuniformity in the shock vicinity is described by the space derivatives of the gasdynamic parameters at a point on the shock surface. The parameters are the gas pressure, density, velocity vector. The derivatives with respect to the directions of the streamline and normal to it, and of the shock surface and normal to it are considered. Space derivatives of all gasdynamic parameters are expressed through the flow nonisobaric factor along the streamline, the streamline curvature, and the flow vorticity and non-isoenthalpy factors. An algorithm for these factors of the gas flow downstream a shock wave determination is developed. Examples of these factors calculation for imperfect oxygen and thermodynamically perfect gas are presented. The influence coefficients of the upstream flow factors on the downstream flow factors are calculated. A
For stars, the bow shock is typically the boundary between their stellar wind and the interstellar medium. Named for the wave made by a ship as it moves through water, the bow shock wave can be created in the space when two streams of gas collide. The space is actually filled with the interstellar medium consisting of tenuous gas and dust. Stars are emitting a flow called stellar wind. Stellar wind eventually bumps into the interstellar medium, creating an interface where the physical conditions such as density and pressure change dramatically, possibly giving rise to a shock wave. Here we discuss some literature on stellar bow shocks and show observations of some of them, enhanced by image processing techniques, in particular by the recently proposed AstroFracTool software.
Self-similar shock solutions in spherically symmetric polytropic gas flows are constructed and analyzed in contexts of proto-star formation processes. Among other possible solutions, we model a similarity shock across the sonic critical curve with an inner free-fall core collapse and a simultaneous outer expansion of the extended envelope; the separation or stagnation surface between these two flow zones travels outwards in a self-similar manner at a variable speed. One readily obtains accretion shock solutions. Semi-complete self-similar solutions across the sonic critical curve either once or twice without shocks can also be constructed. Features of star formation clouds of our polytropic model include the mass density scaling in the outer flow zone $ρ\propto r^{-2/(2-γ)}$, the temperature scalings of the inner flow zone $T\propto r^{-3(γ-1)/2}$ and of the outer flow zone $T\propto r^{-2(γ-1)/(2-γ)}$, and the variable central mass accretion rate $\dot{M}=k^{3/2}t^{(3-3γ)}m_0/G$ where $γ$ is the polytropic index, $k$ is a constant, $m_0$ is the core mass, and $G$ is the gravitational constant. Spectral line profiles characteristic of the `envelope expansion with core collapse' (EE
We theoretically investigate the thermodynamic properties of a quasi-one-dimensional single-component dipolar Fermi gas at finite temperatures. A self-bound fermionic droplet can be achieved by exchange correlations with long-range dipole-dipole interactions under quasi-one-dimensional confinement, where the interaction can be tuned by tilting the dipoles along the system coordinate. Using the Hartree-Fock approximation, we show how the liquid-gas phase transition occurs in this system and elucidate the finite-temperature phase structure consisting of a gas phase, liquid phase, gas-liquid coexistence phase, and spinodal phase. We also discuss its similarity to the liquid-gas phase transition in nuclear matter through a comparison with phenomenological models. By examining the experimental conditions for realizing self-bound fermionic droplets, we find that microwave-shielded fermionic polar molecules are promising candidates. Our results will be useful for an interdisciplinary understanding of self-bound fermionic matter as well as an analog quantum simulation of nuclear systems.
Cyberattacks on pipeline operational technology systems pose growing risks to energy infrastructure. This study develops a physics-informed simulation and optimization framework for analyzing cyber-physical threats in petroleum pipeline networks. The model integrates networked hydraulic dynamics, SCADA-based state estimation, model predictive control (MPC), and a bi-level formulation for stealthy false-data injection (FDI) attacks. Pipeline flow and pressure dynamics are modeled on a directed graph using nodal pressure evolution and edge-based Weymouth-type relations, including control-aware equipment such as valves and compressors. An extended Kalman filter estimates the full network state from partial SCADA telemetry. The controller computes pressure-safe control inputs via MPC under actuator constraints and forecasted demands. Adversarial manipulation is formalized as a bi-level optimization problem where an attacker perturbs sensor data to degrade throughput while remaining undetected by bad-data detectors. This attack-control interaction is solved via Karush-Kuhn-Tucker (KKT) reformulation, which results in a tractable mixed-integer quadratic program. Test gas pipeline case st
We demonstrate that the liquid-gas transition of nuclear matter can be rigorously described with the quantum chromodynamics by combining the quark gap equation and the Faddeev equation of nucleon. Our investigation focuses on this transition at zero temperature and finite chemical potential, revealing a finite difference between the gas and liquid solution of the quark propagator. This difference emerges from the shift of the nucleon pole mass in medium, which is generated in the nucleon channel of the quark gap equation. We prove that such a difference is precisely the contour contribution from the shift of the nucleon pole. The resulting discontinuity manifests as a first-order phase transition and fundamentally determines both the nuclear binding energy and the saturation density. We then derive an analytical relation between the binding energy and the sigma term of the nucleon, yielding a binding energy of $E/A=15.9\,\textrm{MeV}$. Furthermore, by establishing the relation between the nuclear saturation density and the vector charge of nucleon in association with the binding energy, we determine the saturation density to be $n_{\textrm{B}}^{0}=0.15\,\textrm{fm}^{-3}$.
We aim to investigate the interplay between dense molecular gas, star formation, and active galactic nucleus (AGN) feedback in the luminous infrared galaxy (LIRG) NGC 1068 at sub-kiloparsec scales. We present the HCN (4-3) and HCO$^+$ (4-3) maps of NGC 1068, obtained with JCMT as part of the Mapping the dense molecular gas in the strongest star-forming galaxies (MALATANG) project, and perform spatially resolved analyses of their correlations with infrared luminosity and soft X-ray emission. Spatially resolved relations between the luminosities of infrared dust emission and dense molecular gas tracers ($L_{\rm IR}-L'_{\rm dense}$) are found to be nearly linear, without clear evidence of excess contributions from AGN activity. The spatially resolved X-ray emission ($L^{\rm gas}_{0.5-2\,\mathrm{keV}}$) displays a radially-dependent twofold correlation with the star formation rate (SFR), suggesting distinct gas-heating mechanisms between the galaxy center and the outer regions. A super-linear scaling is obtained in galactic center regions with SFR surface density ($Σ_{\rm SFR}$) $>$ 8.2 $\times$ 10$^{-6}$ $M_\odot$ yr$^{-1}$ kpc$^{-2}$: log($L^{\rm gas}_{0.5-2\,\mathrm{keV}}$/erg s$
Noble gas atoms such as Kr and Xe are byproducts of nuclear fission in nuclear plants. How to trap and confine these volatile even radioactive gases is particularly challenging. Recent studies have shown that they can be trapped in nanocages of ultrathin silica. Here, we exhibit with self-consistent phonon theory and four-phonon (4ph) scattering where the adsorption of noble gases results in an anomalous increase in lattice thermal conductivity, while the presence of Cu atoms doping leads to a reduction in lattice thermal conductivity. We trace this behavior in host-guest 2D silica to an interplay of tensile strain, rattling phonon modes, and redistribution of electrons. We also find that 4ph scatterings play indispensable roles in the lattice thermal conductivity of 2D silica. Our work illustrates the microscopic heat transfer mechanism in 2D silica nanocages with the immobilization of noble gas atoms and inspires further exploring materials with the kagome and glasslike lattice thermal conductivity.
Manipulation of intense pulse propagation in gas-filled capillaries is desirable for various high-field applications. Tuning the parameters of the driving laser pulse and the working gas is the conventional approach, and it provides limited capability of control. Here we demonstrate through numerical simulations a practical scheme to control the propagation of intense pulses. A weak ultraviolet pulse is launched into a capillary with a negative delay with respect to a main infrared pulse. The pulses begin to temporally overlap due to dispersion. As the main pulse self-compresses, the control pulse is strongly red-shifted due to cross-phase modulation. The frequency shifts of the two pulses mitigate pulse walk-off and allow an efficient coupling, substantially extending the effective interaction length. This interesting phenomenon may benefit applications such as high-order harmonic generation.
Strong mid-infrared light-matter interaction has attracted extensive attention as it opens up new frontiers in nonlinear optics. Here we observe through simulations a novel aspect of mid-infrared pulse dynamics in a high-pressure gas-filled capillary, where a pulse with a power well below the critical power for Kerr self-focusing undergoes an astonishing rise of the peak intensity following an extremely efficient spectral broadening. This intensity enhancement is attributed to the Kerr-induced focusing of the supercontinuum. Our study provides an interesting perspective for controlling the laser intensity with possible applications in nonlinear light conversion driven by mid-infrared pulses.
We design a neural network Ansatz for variationally finding the ground-state wave function of the Homogeneous Electron Gas, a fundamental model in the physics of extended systems of interacting fermions. We study the spin-polarised and paramagnetic phases with 7, 14 and 19 electrons over a broad range of densities from $r_s=1$ to $r_s=100$, obtaining similar or higher accuracy compared to a state-of-the-art iterative backflow baseline even in the challenging regime of very strong correlation. Our work extends previous applications of neural network Ansätze to molecular systems with methods for handling periodic boundary conditions, and makes two notable changes to improve performance: splitting the pairwise streams by spin alignment and generating backflow coordinates for the orbitals from the network. We illustrate the advantage of our high quality wave functions in computing the reduced single particle density matrix. This contribution establishes neural network models as flexible and high precision Ansätze for periodic electronic systems, an important step towards applications to crystalline solids.
The physical mechanisms for starburst or quenching in less massive ($M_* < 10^{10} M_{\odot}$) galaxies are unclear. The merger is one of the inescapable processes referred to as both starburst and quenching in massive galaxies. However, the effects of the merger on star formation in dwarf galaxies and their evolution results are still uncertain. We aim to explore how to trigger and quench star formation in dwarf galaxies by studying the metal-poor gas-rich dwarf mergers based on the multi-band observations at a spatial resolution of $\sim$ 460 pc. We use the archival data of ALMA (band 3, 8) and VLT/MUSE to map CO($J=$1-0), [CI]($^3$P$_1 - ^3$P$_0$), and H$α$ emission in one of the most extreme starburst merging dwarf galaxies, Haro 11. We find the molecular gas is assembled around the central two star-forming regions. The molecular/ionized gas and stellar components show complex kinematics, indicating that the gas is probably at a combined stage of collision of clouds and feedback from star formation. The peak location and distribution of [CI](1-0) strongly resemble the CO(1-0) emission, meaning that it might trace the same molecular gas as CO in such a dwarf merger starburst
We complement the MALATANG sample of dense gas in nearby galaxies with archival observations of $^{12}\rm CO$ and its isotopologues to determine scaling relations between Wide-field Infrared Survey Explorer (WISE) 12 $μ$m emission and molecular gas tracers at sub-kiloparsec scales. We find that 12 $μ$m luminosity is more tightly correlated with $^{12}\rm CO$ than it is with $^{13}\rm CO$ or dense gas tracers. Residuals between predicted and observed $^{12}\rm CO$ are only weakly correlated with molecular gas mass surface density ($Σ_{\rm mol}$) in regions where $Σ_{\rm mol}$ is very low ($\sim 10~{\rm M_{\odot}~pc^{-2}}$). Above this limit, the $^{12}\rm CO$ residuals show no correlations with physical conditions of molecular gas, while $^{13}\rm CO$ residuals depend on the gas optical depth and temperature. By analyzing differences from galaxy to galaxy, we confirm that the $^{12}\rm CO$-12 $μ$m relation is strong and statistically robust with respect to star forming galaxies and AGN hosts. These results suggest that WISE 12 $μ$m emission can be used to trace total molecular gas instead of dense molecular gas, likely because polycyclic aromatic hydrocarbons (PAHs, a major contribu
We present new high-resolution HI observations of the disk of the collisional infrared luminous (L$_{\rm IR}=2.2\times10^{11}$ L$_{\sun}$) galaxy NGC 1144, which reveal an apparent large-scale azimuthal and kinematic segregation of neutral hydrogen relative to the molecular gas distribution. Even among violently collisional galaxies, the CO/HI asymmetry in NGC 1144 is unusual, both in the inner regions, and in the outer disk. We suggest that we are observing Arp 118 at a special moment, shortly after a high-speed collision between NGC 1144 and its elliptical companion NGC 1143. HI emission with an average molecular fraction f$_{mol}$ $<$ 0.5 is observed on one side (NW) of the rotating disk of NGC 1144, while the other side (SE) is dominated by dense molecular complexes in which f$_{mol}$ is almost unity. The interface region between the warm-- and cool--cloud dominated regions, lies on a deep spiral-like dust-lane which we identify as a shock-wave responsible for the relative shift in the dominance of HI and H$_2$ gas. A strong shock being fed by diffuse HI clouds with unusually large ($>$ 400 km s$^{-1}$) rotational velocities can explain: 1) the CO/HI asymmetries, 2) a lar
Sub-femtosecond x-ray burst is powerful tool for probing and imaging electronic and concomitant atomic motion in attosecond physics. For years, x-ray source (above 2 keV) had mainly been obtained from X-ray free electron laser (XFEL) or synchrotron radiation, which are high energy consumption, high cost and huge volume. Here we propose a low-cost and small-size method to generate X-ray source. We experimentally obtained high photon energy spectrum (~ 5.2 keV) through both atom and multiple valence state ions using a near-infrared 1.45 μm driving laser interacting with krypton gas, according to our knowledge, which is the highest photon energy generated through high-order harmonic generation up to now. In our scheme, multi-keV photon energy can be achieved with a relaxed requirement on experimental conditions, and make time-resolved studies more accessible to many laboratories that are capable of producing high energy photon extending to hard x-ray region. Furthermore, our scheme minimizes the influence of X-ray fluorescence process on detection, and can also be utilized to study the quantum-optical nature of high-order harmonic generation.
We study numerically the influence of disorder and localization effects on the local spectroscopic characteristics and infrared optical properties of $\mbox{Ga}_{1-x}\mbox{Mn}_x\mbox{As}$. We treat the band structure and disorder effects at an equal level by using exact diagonalization supercell simulation method. This method accurately describes the low doping limit and gives a clear picture of the transition to higher dopings, which captures the localization effects inaccessible to other theoretical methods commonly used. Our simulations capture the rich mid-gap localized states observed in scanning tunneling microscopy studies and reproduce the observed features of the infrared optical absorption experiments. We show clear evidence of a disordered valence band model for metallic samples in which (i) there is no impurity band detached from the valence band, (ii) the disorder tends to localize and pull states near the top of the valence band into the gap region, and (iii) the Fermi energy is located deep in the delocalized region away from the mobility edge. We identify localized states deep in the gap region by visualizing the probability distribution of the quasiparticles and co