2024-12-03 | | Total: 17
We present analytical and simulation study of twisted terahertz (THz) radiation generation via propagation of a circularly polarized Laguerre Gaussian (LG) laser pulse in homogeneous plasma embedded in an axial magnetic field. Analytical formulation is based on perturbation technique and quasistatic approximation. Longitudinal and transverse wakefields generated via laser plasma interactions are evaluated using Lorentz force and Maxwells equations in the mildly nonlinear regime. It is observed that two linearly polarized twisted terahertz (THz) radiation beams are generated in mutually perpendicular planes. Superposition of the two beams result in a single linearly polarized twisted THz radiation beam with modified amplitude and polarization direction. Three dimensional (3D) particle in cell (PIC) simulations are performed for this configuration using FBPIC code. Graphical comparison of amplitude of the resultant THz beam obtained via analytical and simulation studies is presented.
Turbulence grounds aircraft and combating it in flight requires energy, yet volant wildlife fly effortlessly even on windy days. The nature of the interactions between soaring birds and transient turbulent gusts is not clear, especially when compared with our understanding of flight in larger and steadier airflows during thermal or dynamic soaring. We show that soaring golden eagles (Aquila chrysaetos) experienced short upward accelerations indicative of preferential engagement with strong and intermittent turbulent updrafts. The vertical accelerations reflect changes in lift that were as large as 25 standard deviations from the mean, or more than three times the acceleration of gravity, and so large as not to be consistent with gust mitigation or avoidance. These extreme events occurred in short bursts that mimic movement with turbulent vortices. The burst statistics and their symmetries approach those of turbulence toward longer timescales. On the shortest timescales, the bursts break the symmetry of small-scale turbulence in favor of upward accelerations that are more intermittent than turbulence. We introduce a simple nonlinear model that predicts the scale at which symmetry breaks and the stronger intermittency on the smaller scales. These findings suggest a ratcheting mechanism on turbulent gusts and constitute the first quantitative evidence in favor of turbulent gust harvesting by wildlife. An implications is that turbulence is so strong and pervasive as to make unsteady and nonlinear aerodynamics an intrinsic and beneficial aspect of both flapping and soaring flight in the atmospheric boundary layer - one that we need to incorporate in our understanding of the energetics of flight.
Machine learning based surrogate models offer researchers powerful tools for accelerating simulation-based workflows. However, as standard datasets in this space often cover small classes of physical behavior, it can be difficult to evaluate the efficacy of new approaches. To address this gap, we introduce the Well: a large-scale collection of datasets containing numerical simulations of a wide variety of spatiotemporal physical systems. The Well draws from domain experts and numerical software developers to provide 15TB of data across 16 datasets covering diverse domains such as biological systems, fluid dynamics, acoustic scattering, as well as magneto-hydrodynamic simulations of extra-galactic fluids or supernova explosions. These datasets can be used individually or as part of a broader benchmark suite. To facilitate usage of the Well, we provide a unified PyTorch interface for training and evaluating models. We demonstrate the function of this library by introducing example baselines that highlight the new challenges posed by the complex dynamics of the Well. The code and data is available at https://github.com/PolymathicAI/the_well.
In this paper, we introduce the proper latent decomposition (PLD) as a generalization of the proper orthogonal decomposition (POD) on manifolds. PLD is a nonlinear reduced-order modeling technique for compressing high-dimensional data into nonlinear coordinates. First, we compute a reduced set of intrinsic coordinates (latent space) to accurately describe a flow with fewer degrees of freedom than the numerical discretization. The latent space, which is geometrically a manifold, is inferred by an autoencoder. Second, we leverage tools from differential geometry to develop numerical methods for operating directly on the latent space; namely, a metric-constrained Eikonal solver for distance computations. With this proposed numerical framework, we propose an algorithm to perform PLD on the manifold. Third, we demonstrate results for a laminar flow case and the turbulent Kolmogorov flow. For the laminar flow case, we are able to identify a semi-analytical expression for the solution of Navier-Stokes; in the Kolmogorov flow case, we are able to identify a dominant mode that exhibits physical structures, which are compared with POD. This work opens opportunities for analyzing autoencoders and latent spaces, nonlinear reduced-order modeling and scientific insights into the structure of high-dimensional data.
We are concerned with the compressible barotropic Navier--Stokes equations for a $\gamma$-law gas with density-dependent bulk viscosity coefficient $\lambda=\lambda(\rho)=\rho^\beta$ on the two-dimensional periodic domain $\mathbb{T}^2$. The global existence of weak solutions with initial density bounded away from zero and infinity for $\beta>3$, $\gamma>1$ has been established by Va\uıgant--Kazhikhov [Sib. Math. J. 36 (1995), 1283--1316]. When $\gamma=\beta>3$, the large-time behaviour of the weak solutions and, in particular, the absence of formation of vacuum and concentration of density as $t \to \infty$, has been proved by Perepelitsa [SIAM J. Math. Anal. 39 (2007/08), 1344--1365]. Huang--Li [J. Math. Pures Appl. 106 (2016), 123--154] extended these results by establishing the global existence of weak solutions and large-time behaviour under the assumptions $\beta >3/2$, $1< \gamma<4\beta-3$, and that the initial density stays away from infinity (but may contain vacuum). Improving upon the works listed above, we prove that in the regime of parameters as in Huang--Li, namely that $\beta >3/2$ and $1< \gamma<4\beta-3$, if the density has no vacuum or concentration at $t=0$, then it stays uniformly away from zero and infinity at all later time $t \in ]0,\infty]$. Moreover, under the mere assumption that $\beta>1$ and $\gamma>1$, we establish the global existence of weak solutions, thus pushing the global existence theory of the barotropic Navier--Stokes equations on $\mathbb{T}^2$ to the most general setting to date. One of the key ingredients of our proof is a novel application -- motivated by the recent work due to Danchin--Mucha [Comm. Pure Appl. Math. 76 (2023), 3437--3492] -- of Desjardins' logarithmic interpolation inequality.
Manipulating particles, such as cells and tissues, in a flowing liquid environment is crucial for life science research. Traditional contactless tweezers, although widely used for single-cell manipulation, face several challenges. These include potential damage to the target, restriction to static environments, complex excitation setups, and interference outside the target area. To address these issues, we propose an ``invisible hydrodynamic tweezer'' utilizing near-zero index hydrodynamic metamaterials. This metamaterial-based device creates an equipotential resistance zone, effectively immobilizing particles in flowing fluids without disturbing the external flow field and without causing damage to the targets. Unlike traditional active control methods, our tweezer passively captures and releases particles by adjusting the flow channel, eliminating the need for continuous and stable excitation devices, thereby significantly simplifying the setup complexity. Furthermore, these tweezers can be modularly designed in different sizes to flexibly accommodate various application needs. Simulations and experimental validations demonstrated the non-interfering, stable trapping, and precise movement capabilities of these tweezers. This proposed technique holds significant potential for applications in biomedicine, microfluidics, and environmental monitoring.
Shock tubes are instrumental in studying high-temperature kinetics and simulating high-speed flows. They swiftly elevate the thermodynamic conditions of test gases, making them ideal for examining rapid chemical reactions and generating high-enthalpy flows for aerodynamic research. However, non-ideal effects, stemming from factors like diaphragm opening processes and viscous effects, can significantly influence thermodynamic conditions behind the shock wave. This study investigates the impact of various diaphragm opening patterns on the shock parameters near the driven section end-wall. Experiments were conducted using helium and argon as driver and driven gases, respectively, at pressures ranging from 1.32 to 2.09 bar and temperatures from 1073 to 2126 K behind the reflected shock. High-speed imaging captured different diaphragm rupture profiles, classified into four distinct types based on their dynamics. Results indicate that the initial stages of diaphragm opening, including the rate and profile of opening, play crucial roles in resulting incident shock Mach number and test time. A sigmoid function was employed to fit the diaphragm opening profiles, allowing for accurate categorization and analysis. New correlations were developed to predict the incident shock attenuation rate and post-shock pressure rise, incorporating parameters such as diaphragm opening time, rupture profile constants, and normalized experimental Mach number. The results emphasize the importance of considering diaphragm rupture dynamics in shock tube experiments to achieve accurate predictions of shock parameters.
The stochastic dynamics of small elastic objects in fluid are central to many important and emerging technologies. It is now possible to measure and use the higher modes of motion of elastic structures when driven by Brownian motion alone. Although theoretical descriptions exist for idealized conditions, computing the stochastic multimodal dynamics for the complex conditions of experiment is very challenging. We show that this is possible using deterministic finite element calculations with the fluctuation dissipation theorem by exploring the multimodal stochastic dynamics of a doubly-clamped nanobeam. We use a very general, and flexible, finite-element computational approach to quantify the stochastic dynamics of multiple modes simultaneously using only a single deterministic simulation. We include the experimentally relevant features of an intrinsic tension in the beam and the influence of a nearby rigid boundary on the dynamics through viscous fluid interactions. We quantify the stochastic dynamics of the first eleven flexural modes of the beam when immersed in air or water. We compare the numerical results with theory, where possible, and find excellent agreement. We quantify the limitations of the computational approach and describe its range of applicability. These results pave the way for computational studies of the stochastic dynamics of complex 3D elastic structures in a viscous fluid where theoretical descriptions are not available.
Liquid composites moulding is an important manufacturing technology for fibre reinforced composites, due to its cost-effectiveness. Challenges lie in the optimisation of the process due to the lack of understanding of key characteristic of textile fabrics - permeability. The problem of computing the permeability coefficient can be modelled as the well-known Stokes-Brinkman equation, which introduces a heterogeneous parameter $\beta$ distinguishing macropore regions and fibre-bundle regions. In the present work, we train a Fourier neural operator to learn the nonlinear map from the heterogeneous coefficient $\beta$ to the velocity field $u$, and recover the corresponding macroscopic permeability $K$. This is a challenging inverse problem since both the input and output fields span several order of magnitudes, we introduce different regularization techniques for the loss function and perform a quantitative comparison between them.
Simulating complex hydraulic conditions, particularly two-phase flows over spillway chutes, can be achieved with high accuracy using three-dimensional numerical models. This study investigates the potential for vacuum generation and cavitation phenomena on the Aghchai Dam service spillway through numerical simulations conducted in Flow-3D. The analysis focuses on two specific flow rates, 4400 and 1065 cubic meters per second, as determined by experimental data. The Volume of Fluid (VOF) method is employed to accurately calculate the free surface flow. Simulation results at a discharge rate of 4400 cubic meters per second indicate a high likelihood of cavitation at critical locations, including the ogee curve and the angle transition in the chute channel. These areas require specific mitigation measures to prevent cavitation-induced damage. In contrast, at the lower flow rate of 1065 cubic meters per second, the risk of cavitation is minimal due to reduced flow velocity and the absence of flow separation from the bed. The numerical findings align closely with empirical observations, demonstrating the reliability of the simulation approach in predicting cavitation behavior.
We propose a one-dimensional mixing model to investigate the impact of Ablative Rayleigh-Taylor Instability (ARTI) on compression, addressing the limitations of high-dimensional simulations. In this model, the scale of the mixed region is predicted using an ablative buoyancy drag model and the mass-weighted average field distributions in the direction of acceleration are obtained by interpolating between the initial and final isothermal molecule mixing states based on the mixing extent, defined by the percentage of maximum specific turbulent kinetic energy and a free multiplier f, which varies from chunk mixing (0) to atomic mixing (1). After validation through two-dimensional simulations, the model is integrated into a Lagrangian framework and applied to spherical implosion scenarios. The results indicate that the compression of the inflight shell is diminished, leading to a reduced time interval between the convergence of the main shock and stagnation, potentially providing a measurable signal in single-shot experiments.
A Mesoscale Convective System (MCS) is a collection of thunderstorms that function as a system, representing a widely discussed phenomenon in both the natural sciences and visual effects industries, and embodying the untamed forces of nature.In this paper, we present the first interactive, physically inspired mesoscale thunderstorms simulation model that integrates Grabowski-style cloud microphysics with atmospheric electrification processes. Our model simulates thunderclouds development and lightning flashes within a unified meteorological framework, providing a realistic and interactive approach for graphical applications. By incorporating key physical principles, it effectively links cloud formation, electrification, and lightning generation. The simulation also encompasses various thunderstorm types and their corresponding lightning activities.
This work presents a comprehensive framework for the efficient implementation of finite-volume-based reacting flow solvers, specifically tailored for high speed propulsion applications. Using the exascale computing project (ECP) based AMReX framework, a compressible flow solver for handling high-speed reacting flows is developed. This work is complementary to the existing PeleC solver, emphasizing specific applications that include confined shock-containing flows, stationary and moving shocks and detonations. The framework begins with a detailed exposition of the numerical methods employed, emphasizing their application to complex geometries and their effectiveness in ensuring accurate and stable numerical simulations. Subsequently, an in-depth analysis evaluates the solver's performance across canonical and practical geometries, with particular focus on computational cost and efficiency. The solver's scalability and robustness are demonstrated through practical test cases, including flow path simulations of scramjet engines and detailed analysis of various detonation phenomena.
Linear stability of a locally parallel annular swirling jet laden with particles in a swirl flow combustor is considered. At low Stokes numbers, the eigenspectra of the particle-laden jet with uniform particle concentration shows three unstable modes namely centre, sinuous and varicose modes. As the Stokes number is increased to unity, the growth rates of the centre and shear layer modes reduces compared to that of the unladen swirling jet. The magnitude of the velocity eigenmodes peaks in the vortex core and decays radially outward. The variation in particle concentration occurs mostly in the vortex core and almost none in the shear layer. The strength of flow reversal at the jet centreline is given by the backflow parameter. An increase in the backflow parameter increases the growth rate of the centre mode. Non-uniformity in the base-state particle concentration is introduced using a Gaussian function varying in the radial direction and a reduction in the growth rate of the centre mode is seen compared to the uniform particle concentration profile. When the location of the peak of the base-state particle concentration profile is inside the vortex core, the centre modes are stable. Linearized vorticity budget analysis reveals that this is accompanied by a decrease in the net generation of perturbation vorticity in the axial direction and increased radial and azimuthal perturbation vorticity.
When numerically solving partial differential equations, for a given problem and operating condition, adaptive mesh refinement (AMR) has proven its efficiency to automatically build a discretization achieving a prescribed accuracy at low cost. However, with continuously varying operating conditions, such as those encountered in uncertainty quantification, adapting a mesh for each evaluated condition becomes complex and computationally expensive. To enable more effective error and cost control, this work introduces a novel approach to mesh adaptation. The method consists in building a unique adapted mesh that aims at minimizing the average error for a continuous set operating conditions. In the proposed implementation, this unique mesh is built iteratively, informed by an estimate of the local average error over a reduced set of sample conditions. The effectiveness and performance of the method are demonstrated on a one-dimensional Burgers equation and a two-dimensional Euler scramjet shocked flow configurations.
Rheology of bubble suspensions is critical for the prediction and control of bubbly flows in a wide range of industrial processes. It is well-known that the bubble suspension exhibits a shear-thinning behavior due to the bubble shape deformation under pure shear, but how the shear rheology response to dilatation remains unexplored. Here, we report a compression-thinning behavior that the bubble suspension exhibits a decreasing shear viscosity upon compressing. This peculiar rheological behavior is microscopically due to that a shrinking bubble surface effectively weakens the flow resistance of the surrounding liquid. We theoretically propose a constitutive equation for dilute bubble suspensions considering both shear and dilatation effects, and demonstrate that the contribution of dilatation effect on the shear viscosity can be significant at a changing pressure.
In this study, we designed a spiral-shaped milli-reactor with a T-junction microchannel to generate digital droplets for studying and observing the digital freezing process of droplets. During the study of the recalescence and solidification processes of digital droplets dynamically moving in microchannels, we found that although the digital generation of droplets in our channel aligns well with the literature, achieving the digitalization of the droplet freezing process is very challenging. Even the initial phase of freezing (the recalescence process) exhibits significant randomness. A key feature of the randomness in the freezing process is the nucleation position of droplets within the channel, which significantly impacts the digital characteristics and hinders digital freezing. During the investigation of freezing randomness, we identified five distinct nucleation profiles, which largely determine the evolution of the freezing front and the duration of the recalescence phase. However, upon studying the motion velocity of the freezing front, we found that these velocities are temperature-dependent. This aligns with the results of our phase-field simulations and experimental findings, indicating that the release of latent heat during the recalescence process is stable. Additionally, the randomness in freezing may also stem from the deformation of droplets during the solidification process. In this study, we identified two distinct solidification modes during the freezing phase: one initiating from the droplet's head or tail and the other starting from the middle, with the latter causing significant droplet deformation. Through statistical analysis, we further explored the influence of flow rate variation on the digital clustering of droplet freezing and discovered flow rate parameters that optimize freezing digitalization.