This paper brings to fruition recently found \cite{migdal2023exact} exact reduction of decaying turbulence in the Navier-Stokes equation in $3 + 1$ dimensions to a Number Theory problem of finding the statistical limit of the Euler ensemble. We reformulate the Euler ensemble as a Markov chain and show the equivalence of this formulation to the \textbf{quantum} statistical theory of $N$ fermions on a ring, with an external field related to the random fractions of $\pi$. We find the solution of this system in the turbulent limit $N\to \infty, \nu \to 0$ in terms of a complex trajectory (instanton) providing a saddle point to the path integral over the density of these fermions. This results in an analytic formula for the observable correlation function of vorticity in wavevector space. This is a full solution of decaying turbulence from the first principle without approximations or fitted parameters. The energy spectrum decays as $k^{-\frac{7}{2}}$. We compute resulting integrals in \Mathematica{} and present effective index $n(t) = -t\partial_t \log E $ for the energy decay as a function of time Fig.\ref{fig::NPlot}. The asymptotic value of the effective index in energy decay $n(\infty) = \frac{5}{4}$, but the universal function $n(t)$ is neither constant nor linear.

This manuscript presents a novel approach to estimating the potential within Schrödinger equations, with a particular focus on applications to the Navier-Stokes problem in fluid dynamics. By establishing new theoretical estimates, we delve into the intricate dynamics of fluid flow, aiming to unveil previously obscured aspects of the Navier-Stokes equations. Through a rigorous mathematical framework, we explore how these new potential estimates can provide a fresh perspective on fluid mechanics, contributing to the ongoing quest to solve some of its most persistent challenges.

This study delves into the Schrödinger equation's new invariants, shedding light on their crucial implications for theoretical physics and beyond. Specifically, it explores the behavior of scattering amplitudes and bound states concerning the choice of coordinate systems. The research unveils that while the energy of bound states remains invariant under coordinate transformations, the phase of scattering amplitudes undergoes variations, underscoring its pivotal role in theories reliant on phase normalization. Moreover, these findings have played a pivotal role in constructing estimates for three-dimensional Navier-Stokes equations, enhancing our ability to model complex fluid dynamics with greater precision and reliability. Additionally, in seismic exploration, tomography, and ultrasound imaging, these properties allow researchers to optimize phase selection, facilitating the most effective interpretation of measurement results. This strategic phase manipulation ensures clearer insights into subsurface structures, tissue composition, and fluid dynamics, with the flexibility to seamlessly transition back to the original coordinate system post-interpretation. This interdisciplinary synergy underscores the profound impact of fundamental theoretical research on advancing practical applications across diverse scientific domains.

This work presents a novel approach for the study of the movement of droplets on inclined surfaces under the influence of gravity and chemical heterogeneities. The developed numerical methodology uses data-driven modeling to extend the applicability limits of an analytically derived reduced-order model for the contact line velocity. More specifically, while the reduced-order model is able to capture the effects of the chemical heterogeneities to a satisfactory degree, it does not account for gravity. To alleviate this shortcoming, datasets generated from direct numerical simulations are used to train a data-driven model for the contact line velocity, which is based on the Fourier neural operator and corrects the reduced-order model predictions to match the reference solutions. This hybrid surrogate model, which comprises of both analytical and data-driven components, is then integrated in time to simulate the droplet movement, offering a speedup of five orders of magnitude compared to direct numerical simulations. The performance of this hybrid model is quantified and assessed in different wetting scenarios, by considering various inclination angles and values for the Bond number, demonstrating the accuracy of the predictions as long as the adopted parameters lie within the ranges considered in the training dataset.

We delve into the intricate dynamics governing anisotropic fluids within the cosmological paradigm, elucidating the nuanced influence of conformal anomalies upon the stress-energy tensor. Through the presentation of a sophisticated theoretical framework, we seamlessly intertwine unconventional fluidic properties with profound cosmological considerations, thereby augmenting and refining our comprehension of the primordial universe.

In addition to their importance in statistical thermodynamics, probabilistic entropy measurements are crucial for understanding and analyzing complex systems, with diverse applications in time series and one-dimensional profiles. However, extending these methods to two- and three-dimensional data still requires further development. In this study, we present a new method to classify spatiotemporal processes based on entropy measurements. To test and validate the method, we selected four classes of similar processes related to the evolution of random patterns: dynamic colored noises ((i)~white and (ii)~red); (iii)~weak turbulence from reaction-diffusion; (iv)~hydrodynamic fully developed turbulence, and (v)~plasma turbulence from MHD. Considering seven possible ways to measure entropy from a matrix, we present the method as a parameter space composed of the two best separating measures of the five selected classes. The results highlight better combined performance of Shannon Permutation Entropy ($S^p_H$) and a new approach based on Tsallis Spectral Permutation Entropy ($S^s_q$). Notably, our observations reveal the segregation of reaction terms in this $S^p_H \times S^s_q$ space, a result that identifies specific sectors for each class of dynamic process, and can be used to train machine learning models for automatic classification of complex spatiotemporal patterns.

We derive an ab initio equilibrium statistical mechanics of the gas-liquid-solid contact angle on planar periodic, monodisperse, textured surfaces subject to electrowetting. To that end, we extend an earlier theory that predicts how advance or recession of the contact line amount to distinct first-order phase transitions of the filling state in the ensemble of nearby surface cavities. Upon calculating the individual capacitance of a cavity subject to the influence of its near neighbors, we show how hysteresis, which is manifested by different advancing and receding contact angles, is affected by electrowetting. The analysis reveals nine distinct regimes characterizing contact angle behavior, three of which arise only when a voltage is applied to the conductive liquid drop. As the square voltage is progressively increased, the theory elucidates how the drop occasionally undergoes regime transitions triggering jumps in the contact angle, possibly changing its hysteresis, or saturating it at a value weakly dependent on further voltage growth. To illustrate these phenomena and validate the theory, we confront its predictions with four data sets. A benefit of the theory is to forsake trial-and-error when designing textured surfaces with specific contact angle behavior.

This work is an addition to the previously developed two-dimensional model of the shock-plasma interaction extending it into the third dimension. The model allows tracing the evolution of the state of the hypersonic flow and the shock front refracted at a thermal discontinuity. The advantages of using the spherical coordinate system for this type of the problems include more transparency in interpreting the solution and a shortened calculation procedure, because all the changes to the front become reduced to one distortion component only. Although the vorticity generation triggered at the interface is the consequence of the refraction and tied to the steep changes of the front, it is shown here that it is not because of an instant parameter jump at the interface due to refraction itself.

This review is a first attempt at bringing together various concepts from research on wall- and magnetically-bounded turbulent flows. Brief reviews of both fields are provided: The main similarities identified are coherent (turbulent) structures, flow generation and transport barriers. Examples are provided and discussed.

The paper is devoted to the study of a high-pressure short-arc xenon discharge, which is widely used as a light source. Cathodes of such lamps are often doped with thorium for better cathode operation. As a result, thorium can evaporate into a discharge. The study is based on a model of a short-arc xenon discharge plasma, which includes presence of thorium atoms in the plasma. The data obtained show that thorium atoms strongly affect the characteristics of the plasma, reducing its temperature, changing the ionic composition near the cathode, energy balance and generation of optical radiation. Especially interesting is the case when the discharge electrodes have a different shape of the working surface: the anode has a convex or concave shape; the cathode has a different cone angle. The obtained data make it possible to propose a new way to increase the luminous efficiency of xenon lamps due to the appropriate choice of the electrodes working surface shape, which turned out to be an unexpectedly strong factor of affecting plasma characteristics. This, in turn, allows us to develop xenon lamps with the desired characteristics providing optical radiation in the required range of the spectrum and better use of energy.

We study electrohydrodynamic (EHD) linear (in)stability of microfluidic channel flows, i.e. the stability of interface between two-shearing viscous (perfect) dielectrics exposed to an electric field in large aspect ratio microchannels. We then apply our results to particular microfluidic systems known as electroosmotic (EO) pumps. Our results are detailed analytical expressions for the growth rate of two-dimensional EHD modes in Couette-Poiseuille flows in the limit of small Reynolds number ($R$); the expansions to both zeroth- and first-order-$R$ are considered. The growth rates are complicated functions of viscosity-, height-, density- and dielectric-constant ratio, as well as of wavenumbers and voltages, and to our knowledge have not been presented before in literature. To make the results more useful, e.g., for voltage-control EO pump operations, we also derive equations for the impending voltages of the neutral stability curves that divide stable from unstable regions in voltage-wavenumber stability diagrams. The voltage equations and the stability diagrams are given for all wavenumbers. We finally outline the flow regimes in which our first-order-$R$ voltage corrections could potentially be experimentally measured. Our work gives the insight into the coupling mechanism between electric field and shear flow in parallel-planes channel flows, correcting an erroneous attempt from literature. Our general analysis enables us also to refine the renowned instability mechanism due to viscosity stratification, i.e., the case of pure shear instability, when the first-order-$R$ voltage correction equals zero.

This is the second paper in a cycle investigating the exact solution of loop equations in decaying turbulence. We perform numerical simulations of the Euler ensemble, suggested in the previous work, as a solution to the loop equations. We designed novel algorithms for simulation, which take a small amount of computer RAM so that only the CPU time grows linearly with the size of the system and its statistics. This algorithm allows us to simulate systems with billions of discontinuity points on our fractal curve, dual to the decaying turbulence. For the vorticity correlation function, we obtain \textbf{quantum fractal laws} with regime changes between universal discrete values of indexes. The traditional description of turbulence with fractal or multifractal scaling laws does not apply here. The measured conditional probabilities of fluctuating variables are smooth functions of the logarithm of scale, with statistical errors being negligible at our large number of random samples $T = 167,772,160$ with $N = 200,000,000$ points on a fractal curve. The quantum jumps arise from the analytical distribution of scaling variable $X = \frac{\cot^2(\pi p/q)}{q^2}$ where $\frac{p}{q}$ is a random fraction. This distribution is related to the totient summatory function and is a discontinuous function of $X$. In particular, the energy spectrum is a devil's staircase with tilted uneven steps with slopes between $-1$ and $-2.03$. The logarithmic derivative of energy decay $n(t)$ as a function of time is jumping down the stairs of universal levels between $1.50$ and $1.53$. The quantitative verification of this quantization would require more precise experimental data.

This study explores how the properties of heater materials impact nucleate pool boiling using a comprehensive simulation model. The model integrates well-known heat transfer mechanisms, allowing us to assess the effects of two distinct heater materials. Findings indicate that materials with higher thermal conductivity, such as copper, notably improve refrigeration efficiency in nucle- ate boiling. The study provides insights into the relationship between bubble growth, microlayer recovery beneath a bubble, temperature fluctuations, and heater properties. Comparisons between copper and silicon oxide underscore variations in bubble frequency, attributed to differences in bubble growth time, microlayer recovery time, and material-dependent behavior. The influence of neighboring bubbling sites, especially pronounced in silicon oxide due to its low thermal conductivity and diffusivity, was observed. Temperature variations in such materials became highly visible due to their significant inertia in recovery. Simulation results align well with empirical correlations, confirming the model’s versatility in capturing the intricate phenomena of nucleate pool boiling.

In this paper, we explore the use of superluminal ionization fronts to accelerate and amplify electromagnetic radiation. These fronts are defined as optical boundaries between two regions of a gas, the neutral region and the plasma region, characterized by two different values of the refractive index. For that reason, the front velocity is not necessarily related with the motion of material particles, such as neutral atoms, ions and electrons, which can stay at rest. The fronts can therefore become superluminal without violating causality. In recent years, different experimental configurations, such as the flying-focus, showed that it is possible to create superluminal fronts in the laboratory. These fronts can easily be described theoretically in a special reference frame, called the time-frame, which is used here. In this frame, superluminal fronts reduce to time-refraction, a process that is symmetrical to the well-known optical refraction. We show that propagation through such fronts can lead to considerable frequency shift and energy amplification of probe laser beams. This could eventually be used to develop new sources of turnable radiation.

Runaway electrons (RAEs) generation in high-pressure gases is an important physical phenomenon that influences significantly discharge shapes and properties of initiated plasma. The diffuse discharges formed due to RAEs in air and other gases at atmospheric pressure find wide application. In the present review, theoretical and experimental results that explain the reason for the RAEs occurrence at high pressures are analyzed and recommendations are given on the implementation of conditions under which the runaway electron beam (RAEB) with the highest current can be obtained at atmospheric pressure. Experimental results were obtained using subnanosecond, nanosecond, and submicrosecond generators, including those specially developed for runaway electron generation. The RAEB were recorded by oscilloscopes and collectors with picosecond time resolution. To describe theoretically the phenomenon of continuous electron acceleration, the method of physical kinetics has been used based on the Boltzmann kinetic equation that takes into account the minimum but sufficient number of elementary processes, including shock gas ionization and elastic electron scattering. Results of modeling allowed the main factors to be established that control the RAE appearance, the most important of which is electron scattering on neutral atoms and/or molecules. Theoretical modeling has allowed the influence of various parameters (including the voltage, pressure, gas type, and geometrical characteristics of the discharge gap) to be taken into account. The results of research presented here allow the RAEs accelerators with desirable parameters to be developed and the possibility of obtaining diffuse discharges to be accessed under various conditions. The review consists of the Introduction, 5 Sections, the Conclusion, and the References.

In this paper, for the Reynolds-averaged Navier-Stokes equations, a self-closed turbulence model without any adjustable parameter is formulated and a simple model is proposed. The validation rule for self-closed turbulence model is rigorously derived from the Reynolds-averaged Navier-Stokes equation. The rule is not effected by turbulence modelling on the Reynolds stresses. As an application, two dimensional turbulent boundary layers, and plate turbulent boundary layers have been formulated.

The phenomenon of droplet impact on moving surfaces is widely observed in fields such as transportation, rotating machinery, and inkjet printing. Droplets exhibit non-axisymmetric behavior due to the motion of solid surfaces which significantly determines core parameters such as contact time, maximum spreading radius, and bounding velocity, thereby affecting the efficiency of related applications. In this study, we focus on kinetics and morphology of the non-axisymmetric bouncing behaviors for droplets impacting on a moving superhydrophobic surface (SHPS) within the normal (Wen) and tangential (Wet) Weber numbers. Considering the influences of moving surface on the contact area and contact time, the previous scaling formula for the horizontal velocity of droplets has been improved. Based on the velocity superposition hypothesis, we establish a theoretical model for the ratio of the maximum spreading radius at both ends depending on Wen and Wet. This research provides both experimental and theoretical evidence for understanding and controlling the non-axisymmetric behavior of droplets impacting on moving surfaces.

The influence of a background vortex flow on clustering of floating tracer is addressed. The vortex flow considered is induced by an ellipsoidal vortex evolving in a deformation. The system exhibits various vortex motion regimes: (i) steady state, (ii) oscillation and (iii) rotation of the ellipsoidal vortex core. The latter two induce unsteady velocity field for the tracer thus leading to irregular (chaotic) tracer motion. Superimposing a stochastic divergent velocity field onto the deterministic vortex flow allows us to observe significantly different tracer evolution. We use the method of characteristics to integrate the floating tracer density evolution equation and the Euler-Ito scheme for obtaining the floating tracer trajectories with the random velocity field. The cluster area and cluster mass from the statistical topography were used as the quantitative diagnostics of floating-tracer’s clustering. For the case of the steady ellipsoidal vortex embedded into the deformation flow with a random velocity field component we found that clustering characteristics are weakened by the steady vortex. For the cases of an unsteady ellipsoidal vortex, we observed clustering in the floating-tracer density field if the contribution of the divergent component is higher or equal to that of the rotational (nondivergent) component. Even when the initial floating tracer patch was set on the boundary of oscillating ellipsoidal vortex, we observed formation of clusters. In the case of a rotating ellipsoidal vortex, we also observed pronounced clustering. Thus, we argue that unsteady ellipsoidal vortex regimes (oscillation and rotation), which induce chaotic motion of the nearby passive tracer’s trajectories, are still conducive to clustering of floating tracer observed in the density field despite the intense deformation introduced by strain and shear.

The possibility of a shock wave recovery at a discrete closed interface with a heated gas has been investigated. A two-dimensional model applied to conditions of optical discharges featuring spherical, elliptical, and drop-like configurations demonstrated that non-symmetry in the shock refraction contributes to the specific mechanism of recovery other than simply its compensation. Even though the full restoration of the hypersonic flow state doesn’t occur in a strict sense of it, clear reverse changes toward the initial shape of the shock front eventually take place, thus creating an appearance of a full recovery seen in experiments. From analysis for different interfaces symmetries, the factors determining the recovery dynamics are identified. The results are directly applicable to the problem of energy deposition into a hypersonic flow, however it can be useful anywhere else where the flow modifications following the interaction are important. The dimensionless form of the equations allows applications on any scale other than that demonstrated for the optical discharges.

Numerical simulations of standard situations in the transition region from gas kinetics to fluid dynamics at small Mach numbers indicate a clear dependence of the simulation results on the underlying kinetic model (here: nonlinear and linearized Boltzmann collision operator vs. BGK relaxation model). We develop an improved mathematical framework (trace theory) to explain these differences. In particular we reveal certain deficiencies for the classical BKG system as well as for the standard Navier Stokes approach.

Practically important properties of real liquids, including seawater, are their nonlinear properties, which include a nonlinear acoustic parameter, as well as cavitation strength – a rupture of the continuity of a liquid at high intensities in an acoustic wave. The manifestation of nonlinear effects is greatly facilitated by the presence of various nuclei in the liquid – gas bubbles, foreign particles, and other inclusions of various origins. For practical applications, it is important to study the properties of real liquids with inclusions in them. Due to the difficulties of measuring cavitation strength by acoustic methods, other methods are being sought. It is proposed to use an optoacoustic method associated with the use of laser radiation that causes optical breakdown - optical cavitation, that accompanied by a strong sound generation effect. The connection between optical breakdown thresholds and acoustic cavitation thresholds has been established, which can later be used to identify cavitation thresholds at high static pressures (at great depths in the sea), where the use of acoustic emitters is extremely difficult. The paper studies the bubble size distribution function in seawater, combined with studies of the acoustic nonlinearity parameter and cavitation strength at various depths. The interrelation of these characteristics for seawater is shown.

The studies of laminar unsteady boundary layer flows is crucial for understanding laminar-turbulence transition and origins of turbulence. However, the task of finding its solutions poses a significant challenge. In this paper, we propose a novel approach by introducing a similar transformation to convert the 2D unsteady laminar boundary layer equations into a single partial differential equation. By applying this transformation, we are able to obtain the "exact" similarity solutions for the 2D unsteady laminar boundary layer equations, specifically for the case of flat plate boundary flow. Notably, this is the first time that such an exact solution has been obtained. Application of this exact solution is being used to solve Stokes' first problem in 2D boundary layer. Perspectives on the transition fron Rayleigh solution to Blasius solution is provided.

The discovery of solitonlike coherent structure (SCS) in boundary layer flows is crucial for understanding turbulence origins and in particular for laminar-turbulence transition. However, the task of finding solutions for solitonlike coherent structure from the Navier-Stokes equations poses a significant challenge. In this paper, based on the author's previous work [B.H. Sun, Exact similarity solutions of unsteady laminar boundary layer flows, Preprints 2023, 2023092117], we are able to study convergent flow boundary layers, whose solution encompasses both shock wave and solitary wave solutions, and their superposition gives rise to solitary-like waves, namely solitonlike coherent structure. It is found that the solitonlike coherent structure can only be obtained by combining the Navier-Stokes equations and mass conservation, since the combined equations will have the third order derivatives respect to coordinate $x$, in other words, without the mass conservation condition, the Navier-Stokes equations does not contain solitonlike coherent structure by itself. Finally proved solitonlike coherent structure is a universal motion form of viscous fluid.

We presents results obtained from the numerical simulation of the gas-phase chemical kinetics in atmospheric pressure air non-equilibrium plasmas. In particular we have addressed the effect of pulsed operation mode of a planar dielectric barrier discharge. As it was conjectured, the large difference in the time scales involved in the fast dissociation of molecules in plasmas and their subsequent reactions to produce stable chemical species, makes the presence of a continuously repeated plasma production stage unnecessary and a waste of electrical power and efficiency. Results on NOx remediation, ozone production, water vapor and ammonia dissociation are discussed. A few comparisons with experimental findings in a dielectric barrier discharge reactor already used for applications is briefly addressed too. Our results clearly indicate a pattern for the optimization of the discharge using a carefully designed repetition rate and duty-cycle.

The presented work studies the processes of synthesis of nitrogen-containing structures of ZnO using atmospheric pressure microwave nitrogen plasma and investigates their photocatalytic activity in the processes of degradation of 2,4-dinitrophenol and the antibiotic ciprofloxacin when irradiated with sunlight. The work proposes an effective method for formation of photosensitive ZnO powders. Due to the features of plasma treatment in the open atmosphere of zinc metal microparticles, ZnO structures are formed with sizes from hundreds of nanometers to several micrometers with various micromorphologies. High photoactivity was demonstrated (rate constants 0.036 min-1 and 0.051 min-1) of synthesized ZnO structures during photo-degradation of 2,4-dinitrophenol and ciprofloxacin, respectively, when exposed to solar radiation. Photo-active structures of ZnO synthesized using microwave nitrogen plasma can find application in processes of mineralization of toxic organic compounds.

In previous papers we have shown how Schrödinger’s equation which includes an electromagnetic field interaction can be deduced from a fluid dynamical Lagrangian of a charged potential flow that interacts with an electromagnetic field. The quantum behaviour was derived from Fisher information terms which were added to the classical Lagrangian. It was thus shown that a quantum mechanical system is driven by information and not only electromagnetic fields. This program was applied also to Pauli’s equations by removing the restriction of potential flow and using the Clebsch formalism. Although the analysis was quite successful there were still terms that did not admit interpretation, some of them can be easily traced to the relativistic Dirac theory. Here we repeat the analysis for a relativistic flow, pointing to a new approach for deriving relativistic quantum mechanics.

We study the generalized χ and η cross helicities for non-ideal non-barotropic magnetohydrodynamics (MHD). χ and η, the additional label translation symmetry group, have been used to generalize the cross helicity in ideal flows. Both new helicities are additional topological invariants of ideal MHD. To study there behaviour in non-ideal MHD, we derive the time derivative of both helicities using the non-ideal MHD equations taking viscosity, finite resistivity and heat conduction into account. Physical variables are divided into ideal and non-ideal quantities separately during the mathematical analysis for simplification. The analytical results indicate that χ and η cross helicities are not strict constant of motion in the non-ideal MHD and show the rate of dissipation which is compared to the dissipation of other topological constants of motion.

The paper reports an analytical study of two problems of Rayleigh-Bénard convection(RBC): one in a Newtonian liquid occupying an axi-symmetric cylindrical enclosure(CE) and second in a two-dimensional rectangular enclosure(RE). Linear stability analysis is employed to obtain information on the onset of convection in terms of a critical Rayleigh number, Ra, as a function of the aspect ratio, Ro. The nature and the flow structures of the convective cells, at the onset of convection, are predicted theoretically by plotting the streamlines for different values of Ro. Not so shallow(Ro<1 or Ro≪1), and very shallow enclosures(Ro≫1) are considered. Concepts of cell-width and the wave number are respectively used in these two types of enclosures. An explicit analytical expression for the number of cells as a function of Ro is obtained theoretically. Exact prediction of the critical value of Ra, viz., Rac, as a function of Ro is also made. These are new findings in the case of both CE and RE. The fundamental differences between the problems of CE and RE are also highlighted.

RANS simulations have been broadly used to investigate turbulence in the oceans and atmosphere. Within these environments there are a multitude of tracers undergoing reactions (eg. phytoplankton growth, chemical reactions). The distribution of these reactive tracers is strongly influenced by turbulent mixing. With a 50 member ensemble of two dimensional Rayleigh-Taylor induced turbulent mixing, we show that the dynamics of a reactive tracer growing according to Fisher’s equation are poorly captured by the mean. A fluctuation-dependent sink introduced by Reynolds averaging Fisher’s equation transfers concentration from the mean to the fluctuations. We compare the dynamics of the reactive tracer with those of a passive tracer. The reaction causes the concentration of reactive tracer to increase thereby increasing Fickian diffusion and allowing the reactive tracer to diffuse into turbulent structures that the passive tracer cannot reach. A positive feedback between turbulent mixing and fluctuation growth is identified. We show that eddy viscosity and diffusivity parameterizations fail to capture the bulk trends of the system and identify a need for negative eddy diffusivities. One must therefore be cautious when interpreting RANS results for reactive tracers.

Localized turbulence is common in geophysical flows, where the roles of rotation and stratification are paramount. In this study, we investigate the evolution of a stratified turbulent cloud under rotation. Recognizing that a turbulent cloud is composed of vortices of varying scales and shapes, we start our investigation with a single eddy. Compared to an eddy under pure rotation, the stratified eddy shows the physical manifestation of a known potential vorticity mode, appearing as a static stable vortex. In addition, the expected shift from inertial waves to inertial-gravity waves is observed. In our numerical simulations of the turbulent cloud, carried out at a constant Rossby number over a range of Froude numbers, stratification causes columnar structures to deviate from vertical alignment. This deviation increases with increasing stratification, slowing the expansion rate of the cloud. The observed characteristics of these columnar structures are consistent with the predictions of linear theory, particularly in their tilt angles and vertical growth rates, suggesting a significant influence of inertial-gravity waves. Using Lagrangian particle tracking, we have identified regions where wave activity dominates over turbulence. In scenarios of milder stratification, these inertial-gravity waves are responsible for a significant energy transfer away from the turbulent cloud, a phenomenon that attenuates with increasing stratification.

We present a mixing length based algebraic turbulence model calibrated to pipe flow; the main purpose of the model is to capture the increasing turbulence production-to-dissipation ratio observed in connection with the high Reynolds number transition region. The model includes the mixing length description by Gersten and Herwig which takes the observed variation of the von Kármán number with Reynolds number into account. Pipe wall roughness effects are included in the model. Results are presented for area-averaged (integral) quantities which can be used both as a self-contained model and as initial inlet boundary conditions for computational fluid dynamics simulations.

Steady-state one-dimensional flows of five-component air behind a normal shock wave are considered with one-temperature model. A mathematical model is formulated to describe the relaxation of a five-component air mixture with a one-temperature non-equilibrium approximation. A numerical study of non-equilibrium flows of a reacting five-component air mixture behind shock waves at different heights and velocities of free flow is performed. The contribution of different types of reactions to the overall relaxation of the mixture is discussed, and the distributions of macro-parameters of the flow behind the shock wave front are calculated. The lengths of the relaxation zones behind the shock wave front are compared at different initial conditions.

Maxwell equations governing electromagnetic effects are being shown to be equivalent to the compressible inviscid Navier-Stokes equations applicable in fluid dynamics and representing conservation of mass and linear momentum. The latter applies subject to a generalized Beltrami condition to be satisfied by the magnetic field. This equivalence indicates that the compressible inviscid Navier-Stokes equations are Lorentz invariant as they derive directly from the Lorentz invariant Maxwell equations subject to the same Beltrami condition. In addition, the derivation results provide support for the claim that electromagnetic potentials have physical significance as demonstrated by Aharonov-Bohm effect, and are not only a convenient mathematical formulation.

The decay of the low-mode internal tide due to the superharmonic energy cascade is investigated in a realistically forced global Hybrid Coordinate Ocean Model (HYCOM) simulation with 1/25 degrees (4 km) horizontal grid spacing. Time-mean and depth-integrated supertidal kinetic energy is found to be largest near low-latitude internal tide generation sites, such as the Bay of Bengal, Amazon Shelf, and Mascarene Ridge. The supertidal kinetic energy can make up to 50% of the total internal tide kinetic energy several hundred kilometers from the generation sites. As opposed to the tidal flux divergence, the supertidal flux divergence does not correlate with the barotropic to baroclinic energy conversion. Instead, the time-mean and depth-integrated supertidal flux divergence correlates with the nonlinear kinetic energy transfers from (sub)tidal to supertidal frequency bands as estimated with a novel coarse-graining approach. The regular spaced banding patterns of the surface-intensified nonlinear energy transfers are attributed to semidiurnal mode 1 and mode 2 internal waves that interfere constructively at the surface. This causes patches where both surface tidal KE and nonlinear energy transfers are elevated. The simulated internal tide off the Amazon shelf steepens significantly near these patches, generating solitary-like waves in good agreement with Synthetic Aperture Radar (SAR) imagery. Globally, we find that regions of high supertidal energy flux also show a high correlation with observed instances of solitary nonlinear internal waves.

Registration of the flow pattern fine structure and the matter distribution of a free falling liquid drop in a target fluid at rest in the impact mode of coalescence, when the kinetic energy of the drop exceeds its available potential surface energy (APSE), was carried out by photo and video recording. The main attention was paid to the study of the flow structure at the initial stage of the cavity formation. To carry out color registration, the flow pattern was illuminated by several matrix LED and fiber optic sources of constant light. The planning of experiments and interpretation of the results were based on the properties of the complete solutions of the fundamental equations of fluid mechanics system, including the transfer and conversion of energy processes. Complete solutions of the system of equations describe large-scale flow components that are waves or vortices as well as thin jets (ligaments, filaments, fibers, trickles). In experiments, the jets are accelerated by the converted available potential surface energy (APSE) when the free surfaces of merging fluids were eliminated. The experiments were performed with the coalescence of water, solutions of alizarin ink, potassium permanganate, copper sulphate or iron sulphate drops in deep water. In all cases, at the initial contact, the drop begins to lose its continuity and breaks up into a thin veil and jets, the velocity of which exceeds the drop contact velocity. Small droplets, the size of which grows with time, are thrown into the air from spikes at the jet tops. On the surface of the liquid, the fine jets leave colored traces that form linear and reticular structures. Part of the jets penetrating through the bottom and wall of the cavity forms an intermediate covering layer. The forming layer jets are separated by interfaces of the target fluid. The processes of molecular diffusion equalize the density differences and form an intermediate layer with sharp boundaries in the target fluid. All noted structural features of the flow are also visualized when a fresh water drop isothermally spreads in the same tap water. The fast-changing and finely structured diffuse boundary of merging fluids, which at the initial stage has a complex and irregular shape, is gradually smoothed out by molecular diffusion processes. The similar flow patterns were observed in all performed experiments, however the geometric features of the flow depend on the individual thermodynamic and kinetic parameters of the contacting fluids.

Conversion of electromagnetic energy into magnetohydrodynamic energy occurs when the electric conductivity changes from negligible to finite values. This process is relevant during the epoch of reheating of the early Universe at the end of inflation and before the emergence of the radiation-dominated era. We find that conversion into kinetic and thermal energies is primarily the result of electric energy dissipation and that the magnetic energy plays only a secondary role in this process. This means that, since electric energy dominates over magnetic energy during inflation and reheating, significant amounts of electric energy can be converted into magnetohydrodynamic energy when conductivity emerges early enough, before the relevant length scales become stable.

The nonlocality (superdiffusion) of turbulence is expressed in the empiric Richardson t3 scaling law for the mean square of the mutual separation of a pair of particles in a fluid or gaseous medium. The development of the theory of nonlocality of various processes in physics and other sciences based on the concept of Lévy flights resulted in the idea of Shlesinger and colleagues about the possibility of describing the nonlocality of turbulence using a linear integro-differential equation with a slowly falling kernel. The close approach developed by us made it possible to establish the closeness of the superdiffusion parameter of plasma density fluctuations moving across a strong magnetic field in a tokamak to the Richardson law. In this paper, we show the possibility of a universal description of the characteristics of nonlocality of transfer in a stochastic medium (including turbulence of gases and fluids) using the Biberman-Holstein approach to the transfer of excitation of a medium by photons, generalized to take into account the finiteness of the velocity of excitation carriers. This approach enables us to propose a scaling that generalizes Richardson's t3 scaling law to the combined regime of Lévy flights and Lévy walks.

An experimental platform for laser-driven ion (sub-MeV) acceleration and potential applications was recently commissioned at the HiLASE laser facility. The auxiliary beam of the Bivoj laser system operating at GW peak power (~10 J in 5-10 ns) and 1-10 Hz repetition rate enabled a sta-ble production of high-current ion beams of multiple species (Al, Ti, Fe, Si, Cu, Sn). The pro-duced laser-plasma ion sources were fully characterized against the laser intensity on target (1013-1015 W/cm2) by varying the laser energy, focal spot size, and pulse duration. This al-lowed to provide accurate scaling laws of the maximum ion energy for the different target ma-terials investigated. Such experimental scaling laws are presented for the first time in the inves-tigated laser intensity range and for ns-class laser pulses, and allow to provide a qualitative in-terpretation of the laser-plasma interaction underpinning physics, thus to tune the main features of the accelerated ion beams (energy, temperature, and current). Such a detailed study was facil-itated by the large amount of data acquired at high repetition rate (1-10 Hz) provided by the Bivoj laser system. The versatility and tuneability of such high-repetition-rate laser-plasma ion sources are of po-tential interest for multidisciplinary user applications.

Nonlinear thermal transport of non-Newtonian polymer flows is an increasingly important area in materials engineering. Motivated by new developments in this area which entail more re-fined and mathematical frameworks, the present analysis investigates the boundary layer ap-proximation and heat transfer persuaded by a symmetrical cylindrical surface positioned hori-zontally. To simulate thermal relaxation impacts the bioconvection cross nanofluid flow Buongiorno model is deployed. The impact of the magnetic field applied on the nanofluid using the heat generation and melting phenomena are inspected. The nonlinear effect of thermosolutal buoyant forces is incorporated into the proposed model. The novel mathematical equations comprise thermophoresis and Brownian diffusion effects. Via robust transformation techniques, the primitive resulting partial equations for momentum, energy, concentration, and motile living microorganism are rendered into nonlinear ordinary equations with convective boundary postulates. An explicit and efficient numerical solver procedure in Mathematica 11.0 programming platform is developed to engage the nonlinear equations. The consequence of multiple governing parameters on dimensionless fluid profiles is inspected through plotted visuals and tables. Finally, outcomes declared on the surface drag force, heat, and mass transfer coefficients are exhibited through 3D visuals via different influential parameters.

Wind energy is a rapidly expanding renewable energy technique. Wind farm developers need to understand the interaction between wind farms and the atmospheric flow over complex terrain. Large-eddy simulations provide valuable data for gaining further insight into the impact of rough topography on wind-farm performance. In this research, we investigate the influence of spatial heterogeneity on wind turbine performance. We conducted numerical simulations of a 12×5 wind-turbine array on various rough topographies. First, we evaluated our LES method through mesh convergence analysis, using mean vertical profiles, vertical friction velocity, and resolved and subgrid-scale kinetic energy. Next, we analyze the effects of surface roughness and dispersive stresses on the performance of fully developed large wind farms. Our results demonstrate that the ground roughness element’s flow resistance boosts large wind-farm power production by almost 68% in fully aerodynamic rough surface compared to flat terrain. Dispersive stress analysis revealed that the primary degree of spatial heterogeneity in the wind farm is in the streamwise direction, which is the “wake-occupied” region, and the relative contribution of dispersive shear stress is almost 45% to the overall drag. We also observed that the power performance of the wind farm in complex terrain outperforms the drag. Our study has implications for improving the design of wind turbines and wind farms in complex terrain to increase their efficiency and energy output.

We consider open non-twist Hamiltonian systems represented by an area-preserving two-dimensional map describing incompressible planar flows in the reference frame of a propagating wave, and possessing exits through which map orbits can escape. The corresponding escape basins have a fractal nature that can be revealed by the so-called basin entropy, a novel concept developed to quantify final-state uncertainty in dynamical systems. Since the map considered violates locally the twist condition, there is a shearless barrier that prevents global chaotic transport. In this paper we show that it is possible to determine the shearless barrier breakup by considering the variation of the escape basin entropy with a tunable parameter.

The effect of coal hydrophilic particles in water-glycerol drops on the maximum diameter of spreading along a hydrophobic solid surface is experimentally studied by analyzing the velocity of internal flows by Particle Image Velocimetry. The grinding fineness of coal particles was 45-80 μm and 120-140 μm, their concentration was 0.06 wt.% and 1 wt.%. The impact of particle-laden drops on a solid surface occurred at Weber numbers (We) from 30 to 120. It is revealed the interrelated influence of We and the concentration of coal particles on changes in the maximum absolute velocity of internal flows in a drop within the kinetic and spreading phases of the drop-wall impact. A physical model is formulated for internal convective flows in the longitudinal section of a drop parallel to the plane of the solid wall. The kinetic energy of translational motion of coal particles in a spreading drop compensates for the energy expended by the drop on sliding friction along the wall. At We=120, the inertia-driven spreading of the particle-laden drop is mainly determined by the dynamics of the deformable Taylor rim. An increase in We contributes to more noticeable differences in the convection velocities in spreading drops. When the drop spreading diameter rises at the maximum velocity of internal flows, a growth of the maximum spreading diameter occurs. The presence of coal particles causes a general tendency to reduce the drop spreading.

Mathematical modeling is the most important tool for constructing the theory of different kinds of two-phase flows. The review is devoted to the analysis of the advent of the approach of mathematical modeling of two-phase flows, where solid particles act mainly as the dispersed phase. The main problems and features of the study of gas-solid two-phase flows are excluded. The main characteristics of gas flows with solid particles are given and the classification of two-phase flows is developed on their basis. The Lagrangian and Euler approaches to modeling the motion of a dispersed phase (particles) are described. Much attention is paid to the consideration of numerical simulation methods that describe descriptions of turbulent gas flow at different hierarchical levels (RANS, LES, DNS), different levels of description of interphase interactions (one-way coupling (OWC), two-way coupling (TWC) and four-way coupling (FWC)), as well as at different levels of interface resolution (partial-point (PP) and particle-resolved (PR)). Excluded are examples of some studies carried out on the basis of the identified approaches, as well as their exclusion for mathematical modeling of various classes of gas-solid two-phase flows.

This article examines the wavelet transforms and the machine learning methods in accelerating large eddy simulations of turbulent flows. An overarching challenge of large eddy simulations is to accurately represent the cascade of kinetic energy across the cut-off length scale. Taylor hypothesis suggests that the energy cascade occurs through the process of vortex stretching. However, Richardson hypothesis assumes the self-amplification of the strain field and a hierarchical break down of large eddies lead the energy cascade. The large eddy simulations typically employs the self-amplification of the strain in formulating subgrid models. However, several studies also proposed subgrid models based on vortex stretching. The wavelet-based coherent vortex simulation of turbulence directly accounts for vortex stretching in overall forward scatter of energy, while allowing local backscatter. The wavelet-based large eddy simulation adapts the grid to capture the creation of small-scale eddies, while adopting subgrid models based on the self-amplification of strain. The advancement of artificial intelligence in turbulence modelling is currently evolving around accelerating the numerical simulations of turbulent flow. However, there is a clear connection between the application of wavelet transfroms and neural networks for directly solving the Navier-Stokes equation, indicating some potential benefits of wavelet methods over the neural networks.

We report a preliminary exploration on the low-pressure Ar-O2 plasma generated by low-frequency alternating current (AC) power supply in the low-temperature Plasma Experimental Power Supply (CTP-2000S) for the first time. The electron density 〖(n〗_e), the electron temperature 〖(T〗_e) and the electron energy distribution function (EEDF) are measured by the single Langmuir probe in this experiment. Experimental results show that the ne increases with gas pressure rising, while the ne decreases with the oxygen concentrations rising. And, a plasma density peak is also observed with the AC power supply increase. A comparison between this work and the previous results by RF source suggests that the low-frequency AC power supply may be very fruitful exploration to generate the low-pressure plasma.

Magnetic helicity plays a tremendously important role when it is different from zero on average. Most notably, it leads to the phenomenon of an inverse cascade. Here, we consider decaying magnetohydrodynamic turbulence as well as some less common examples of magnetic evolution under the Hall effect and ambipolar diffusion, as well as cases in which the magnetic field evolution is constrained by the presence of an asymmetry in the number density of chiral fermions, whose spin is systematically either aligned or anti-aligned with its momentum. In all those cases, there is a new conserved quantity: the Hosking integral. We present quantitative scaling results for the magnetic integral scale as well as the magnetic energy density and its spectrum. We also compare with cases were also a magnetic version of the Saffman integral is initially finite.

Features of nanosecond discharge development in a non-uniform electric field are studied experimentally. High spatial resolution imaging showed that thin luminous tracks of great length with a cross-section of a few microns are observed against the background of discharge glow in air and argon. It has been established that the detected tracks are adjacent to brightly luminous white spots on the electrodes or to the vicinity of these spots, and are associated with the flight of small particles. It is shown that the traces have various shapes and change from pulse to pulse. The particle tracks may look like curvy or straight lines. In some photos, they can change the direction of movement to the opposite. The particle trace was found to abruptly cut off and a bright flash is visible at the trail break point. The color of the tracks differs from that of the spark leaders, while the bands of the second positive nitrogen system dominate in the gap spectra during the existence of a diffuse discharge. Areas of blue light are evident near the electrodes, as well. Development of glow in the gap during its breakdown is revealed using an ICCD camera. Physical reasons for the observed phenomena are discussed.

The common dandelion uses a bundle of drag-enhancing bristles (the pappus) that enables seed dispersal over formidable distances; however, the scaling laws of aerodynamic drag underpinning pappus-mediated flight remains unresolved. In this paper, we study the aerodynamic shape of dandelion, derive the scaling law of resistance, determine the Vogel exponent. In particular, we find that the total drag coefficient is proportional to the -2/3 power of the dandelion pappus Reynolds number, and obtain the terminal velocity of the dandelion seed under gravitation field.

Shear thickening fluids (STF) are a class of fluids whose viscosity significantly rises under external loads. The research on these fluids has been advancing in recent years regarding prospective practical applications, including developing impact-absorbing composites. Following the green agendas governments and legislators advocate, this study investigates the combination of STF and other sustainable materials. Cork is a naturally occurring cellular material with a negative carbon footprint and superior energy absorption properties. With varying concentrations of STF material, cork agglomerates were formed. Notably, the analysed blends are homogeneous, contrasting with previous literary works. A series of experiments were conducted in quasi-static and dynamic conditions to determine the various mechanical responses. STF appears to influence the mechanical behaviour of cork agglomerates by providing softer deceleration rates and dispersing more energy through disaggregation mechanisms.

This article provides insight into flow measurement techniques in water injection wells in oil production fields, with a particular focus on the initial phases of operation. Consequently, the method created by Ramey in 1962, originally intended for estimating the temperature of the injection fluid, has been adapted to calculate the flow rate. In this technique, the calculation This procedure is based on the correlation between the thermal flux formed in the well. The discrepancy between the temperatures of the injected liquid and the geothermal temperature of the reservoir is the main source of the systematic errors in Ramey’s technique. To a lesser extent, but still significant, failure to observe the injection time in a fluid variation also results in an error problem that needs the failure to adhere to the scheduled injection time for a fluid alteration also yields a notable error dilemma that needs to be fixed. The reduction of listed systematic errors is the product of the main part of this article.

Accurate measurement of low-level thoron gas and high-accuracy calibration of thoron measurement devices are essential for assessing and preventing thoron radiological risks. This study aimed to develop a thoron gas measurement technique using an airflow-through scintillation cell for both low-level measurement and high-accuracy calibration. To achieve this, a compartment model was developed to estimate the influence of progeny deposition and accumulation on the wall of the scintillation cell to prevent overestimation of thoron. A self-developed scintillation cell was utilized to implement and validate this technique. The lower detection limit and measurement uncertainty were then evaluated to assess the feasibility of the technique for low-level measurement and high-accuracy calibration. The results showed that the compartment model effectively addressed the influence of the progeny deposition. The measurement technique achieved a lower detection limit below 100 Bq m-3 even with the coexistence of 100 Bq m-3 of radon and attained a measurement uncertainty (k = 2) below 10% when the concentration of thoron exceeded 1,000 Bq m-3. In summary, this study developed a reliable and practical thoron gas measurement technique using an airflow-through scintillation cell with consideration for progeny deposition, and is expected to contribute to the assessment and prevention of thoron radiological risk.

Projected future sea level rise and marine storminess is a serious threat for beaches as they induce beach flooding and erosion. Among other factors, wave runup play an important role in beach evolution and must be robustly assessed. However, little attention has been paid to the uncertainties associated to its characterization and how do they compare to other sources of uncertainty. As a first step towards this goal, we have quantified the impact of several sources of error in the estimation of wave runup on sandy beaches. Understanding what factors are more influential in the quality of the results will help to determine the main sources of uncertainty in beach flooding projections. To reach that goal, a calibrated state-of-the-art numerical modelling system has been setup for two beaches in the Mallorca islands (NW Mediterranean). The system has been forced with the best available information of nearshore incoming waves and has been validated against observations to define the benchmark quality. To determine the key factors affecting the quality of the system’s results, different systems configurations have been tested with different degrees of complexity.Our results show that using the most sophisticated modelling system with the best information on boundary conditions, bottom bathymetry, and submerged vegetation leads to a swash RMSE comparable to the standard deviation of the observed swash. We have also found that the choice of lateral boundary conditions (i.e., source of information for the incoming waves) can double the RMSE and induce large biases. Our results also show that using a simple empirical approach usually underestimates the wave runup. However, in locations with vegetated seabed there is a compensation error, and the empirical approach can lead to acceptable results if forced by nearshore waves. In addition, we have compared the error estimates obtained with uncertainties associated to projected sea level rise. Our findings suggest that the uncertainty associated with wave runup modelling should be considered in the assessment of the total uncertainty of future beach flooding, since the analysis performed indicates that this uncertainty can account for between 16% to 60% of the uncertainties linked to mean sea level projections.

Anatomical airways parameters as length, diameter and angles have a strong effect on the flow dynamics. Aiming to explore the effect of variations of the bifurcation angle (BA) and carina rounding radius (CRR) of lower human airways on the respiratory processes, numerical simulations of airflow during inhalation and exhalation were performed using a synthetic bifurcation models. Geometries for the airways models were parameterized based on a set of different BA’s and several CRR’s. A range of Reynolds numbers (Re) relevant to the human breathing process were selected to analyzed the airflow behaviour. The numerical results show a significant influence of BA and the CRR on the development of the airflow within the airways, and therefore affecting some relevant features of the flow, namely the deformation of velocity profiles, alterations of pressure drop, the secondary flow patterns, and finally enhance or attenuation of wall shear stresses (WSS) appearing during the regular respiratory process. A particular relationship between velocity profiles, vorticity and the secondary flow patterns is also discussed.

The evolution in the post shock nonequilibrium relaxation in a hypersonic plasma flow was investigated during a shock’s reflection off a thermal discontinuity. Within a transitional period, the relaxation zone parameters past both, the reflected and transmitted waves, evolve differently compared to that in the incident wave. In a numerical example for the non-dissociating N₂ gas heated to 5000 K/10,000 K across the interface and M = 3.5, the relaxation time for the transmitted wave is up to 50% shorter and the relaxation depth for both waves is significantly reduced, thus resulting in a weakened wave structure. The results of extension into larger values of heating strength and the shock Mach numbers are discussed. The findings can be useful in the areas of research involving strong shocks interacting with optical discharges or other heated media on the scale where the shock structure becomes important.

The turbulent Poiseuille flow between two parallel plates is one of the simplest possible physical situations, and it has been studied intensively. In this paper, we propose a modified Prandtl-van Driest mixing length that satisfies both boundary conditions and wall damping effects. With our new formulations, we numerically solve the problem and, moreover, propose an approximate analytical solution of mean velocity. As applications of our solution, an approximate analytical friction coefficient of turbulent Poiseuille flow is proposed.

A formula for the turbulent Prandtl number is derived in terms of local variables available from two-equations turbulence models. The derivation is a direct consequence of the expected square additivity of the molecular and turbulent parameters defining the effective viscosity and the effective conductivity. The formula does not degenerate and leads to a Kays like formulation if approximated. Its utility and its limitation are discussed.

We present large eddy simulations of a mid-latitude open ocean front using a modified state of the art computational fluid dynamics code. We investigate the energy and information fluxes at the submesoscale/small-scale range in the absence of any atmospheric forcing. We find submesoscale conditions (Ro∼1, Ri∼1) near surface, within baroclinic structures related with partially imbalanced frontogenetic activity. Near surface, the simulations show a significant scale coupling on scales larger than ∼103(m). This is manifested as a strong direct energy cascade and intense mutual communication between scales, where the latter was evaluated using an estimator based on Mutual Information Theory. At scales smaller than ∼103(m) the results show near-zero energy flux, however, at this scale range, the estimator of mutual communication still shows values corresponding with a significant level of communication between them. This fact motivates to investigate the nature of the self-organized turbulent motion at this scale range with weak energetic coupling but where communication between scales is still significant and to inquire into the existence of synchronization or functional relationships between scales with emphasis on eventual underlying non-local processes.

Heavy Ion Beam Probing (HIBP) diagnostic is a powerful tool for electric field studies in hot dense plasma of modern day toroidal magnetic confinement devices. On the TUMAN-3M to-kamak, the HIBP have been used in regimes with improved plasma confinement to clear up the role of radial electric field in the transition to good confinement regimes. Recently, a moderni-zation of the TUMAN-3M HIBP diagnostics was performed aiming to reconfigure it for a work with a reversed plasma current direction and improvement of overall stability of the diagnostic. The results of first measurements of plasma potential in co-NBI scenario are reported and dis-cussed.

We investigated experimentally, analytically and numerically the formation process of double emulsion formations under dripping regime in a tri-axial co-flow capillary device. The results show that mismatches of core and shell droplets under a given flow condition can be captured both experimentally and numerically. We propose a semi-analytical model using the match ratio between the pinch-off length of the shell droplet and the product of the core growth rate and its pinch-off time. The mismatch issue can be avoided if the match ratio is lower than unity. We considered a model with the wall effect to predict the size of the matched double emulsion. The model shows slight deviations with experimental data if the Reynolds number of continuous phase is lower than 0.06, but asymptotically approaches to good agreement if the Reynolds number increases from 0.06 to 0.14. The numerical simulation generally agrees with the experiments under various flow conditions.

When existing at emulsion interface and continuous phase, dispersibility of water-soluble flexible polymer chains has obvious effect on rheology and dielectric properties of whole emulsion. CNF Pickering emulsion is a good system to research these properties with respect to their microscopic phase structure, dielectric and rheological properties by using cellulose nanofibers (CNF) as water-soluble Pickering emulsifier, liquid Paraffin as oil phase and DDAB as a cationic auxiliary surfactant. The CNF- and DDAB-contents were systematically varied while the water to paraffin oil ratio remained constant to discern the influence of the Pickering emulsifiers. Polarized optical microscopic images revealed that the droplets have tendency to become smaller size for higher CNF content, but bigger size for higher DDAB content, which was proved by the fluorescent analysis for CNF dispersibility with varying DDAB content. The dielectric damping exhibits a minimum, whose value decreases with increasing DDAB- and CNF-content. Increasing the DDAB-content promotes solubilization of CNFs in the aqueous phase and will increase the overall viscosity and yield points. Similarly, a higher CNF-content leads to a higher viscosity and yield point, but at high DDAB-contents the viscosity function exhibits an S-shape at intermediate CNF-content. To evaluate the results further they were compared with CNF-dispersions (without oil-phase), which showed a surfactant effect slightly on maximum stress but strongly on yield stress τy, indicating DDAB can promote the formation of CNF network rather than the viscosity of whole system.

We study in this work a dilute granular gas immersed in a thermal bath made of smaller particles with nonnegligible masses as compared with the granular ones. The kinetic theory for this system is developed and described by an Enskog--Fokker--Planck equation for the one-particle velocity distribution function. Granular particles are assumed to have inelastic and hard interactions, losing energy in collisions as accounted by a constant coefficient of normal restitution. The interaction with the thermal bath is based on a nonlinear drag force plus a white-noise stochastic force. To get explicit results of the temperature aging and steady states, Maxwellian and first Sonine approximations are developed. The latter takes into account the coupling of the excess kurtosis with the temperature. Theoretical predictions are compared with direct simulation Monte Carlo and event-driven molecular dynamics simulations. While good results for the granular temperature are obtained from the Maxwellian approximation, a much better agreement, especially as inelasticity and drag nonlinearity increase, is found when using the first Sonine approximation. The latter approximation is, additionally, crucial to account for memory effects like Mpemba and Kovacs-like ones.

The unsteady flow around the monopile under the influence of wave and vibration is studied by numerical method. Based on VOF model and SST Turbulence model, combined with Intermittent Transition model and dynamic mesh method, wave generation and pile vibration in numerical flume are simulated. The flow field around vibrating monopile in a calm water flume and a wave flume is simulated and the simulation results are compared. The calculation results show that the monopile vibration has a great disturbance to the unsteady flow around the pile, and the disturbance degree intensifies with the increase of the vibration amplitude of the monopile: Under the action of inertia, the vibration of the monopile leads to the reattachment of the shedding vortices; Under the influence of the reattachment vortices, the asymmetry of the vorticity on the side wall of the monopile intensifies; The intensification of the disturbance causes the reduction of the phase lead angle of the bed shear stress on the side of the single pile, and an additional high bed shear stress appears on the side of the monopile. In addition, the vortex structure formed alternately around the monopile during the vibration process leads to the periodic change of the high bed shear stress.

An overview is presented of several diverse branches of work in the area of effectively 2D fluid equilibria which have in common that they are constrained by an infinite number of conservation laws. Broad concepts, and the enormous variety of physical phenomena that can be explored, are highlighted. These span, roughly in order of increasing complexity, Euler flow, nonlinear Rossby waves, 3D axisymmetric flow, shallow water dynamics, and 2D magnetohydrodynamics. The classical field theories describing these systems bear some resemblance to perhaps more familiar fluctuating membrane and continuous spin models, but the fluid physics drives these models into unconventional regimes exhibiting large scale jet and eddy structures. From a dynamical point of view these structures are the end result of various conserved variable forward and inverse cascades. The resulting balance between large scale structure and small scale fluctuations is controlled by the competition between energy and entropy in the system free energy, in turn highly tunable through setting the values of the conserved integrals. Although the statistical mechanical description of such systems is fully self-consistent, with remarkable mathematical structure and diversity of solutions, great care must be taken because the underlying assumptions, especially ergodicity, can be violated or at minimum lead to exceedingly long equilibration times. Generalization of the theory to include weak driving and dissipation (e.g., non-equilibrium statistical mechanics and associated linear response formalism) could provide additional insights, but has yet to be properly explored.

Mesoscale imperfections, such as pores and voids, can strongly modify the properties and the mechanical response of materials under extreme conditions. Tracking the material response and microstructure evolution during void collapse is crucial for understanding its performance. In particular, imperfections in ablator materials, such as voids, can limit the efficiency of the fusion reaction and ultimately hinder ignition. To characterize how voids influence the response of materials during dynamic loading and seed hydrodynamic instabilities, we have developed a tailored fabrication procedure for designer targets with voids at specific locations. Our procedure uses SU-8 as a proxy for ablators materials and hollow silica microspheres as proxy for voids and pores. By using photolithography to design the targets’ geometry, we demonstrate precise and highly reproducible placement of a single void within the sample, which is key for a detailed understanding of its behavior under shock compression. This fabrication technique will benefit high-repetition rate experiments at x-ray and laser facilities. Insight from shock compression experiments will provide benchmarks for the next generation of microphysics modelling.

A century-old scientific conundrum is solved in this paper. The Prandtl mixing length modelled plane boundary turbulent flow is described by: $\frac{du^+}{d{y^+} }+\kappa^2{y^+} ^2\left( \frac{du^+}{d{y^+} }\right)^2=1$, together with boundary condition $ {y^+} =0:\, u^+=0$. Only an approximate solution to this nonlinear ordinary differential equation (ODE) has been sought so far, however, the exact solution to this ODE has not been obtained. By introducing a transformation, $2\kappa y^+=\sinh \xi$, I successfully find the exact solution of the ODE as follows: $u^+=\frac{1}{\kappa}\ln(2\kappa {y^+} +\sqrt{1+4\kappa^2{y^+} ^2})-\frac{2{y^+} }{1+\sqrt{1+4\kappa^2{y^+} ^2}}$.

In this paper, we present the extension of the pressure-based solver designed for the simulation of compressible and/or incompressible two-phase flows of viscous fluids. The core of the numerical scheme is based on the hybrid Kurganov— Noele — Petrova/PIMPLE algorithm. The governing equations are discretized in the conservative form and solved for velocity and pressure, with the density evaluated by an equation of state. The acoustic-conservative interface discretization technique helps to prevent the unphysical instabilities on the interface. The solver was validated on various cases in wide range of Mach number, both for single-phase and two-phase flows. The numerical algorithm was implemented on the basis of the well-known open-source Computational Fluid Dynamics library OpenFOAM in the solver called interTwoPhaseCentralFoam. The source code and the pack of test cases are available on GitHub: https://github.com/unicfdlab/hybridCentralSolvers

Using an original method for measuring weak radiothermal electromagnetic signals, experimental data were obtained on the analysis of dispersions of noise attractors of the radiobrightness temperature of aqueous solutions of NaCl, Al (NO₃)₃ at dilutions C1 (1:100) – C12 (12:100). The values of dispersions of noise attractors significantly differ from those for the control experiment with water. On the basis of the changed values of the noise attractors of the radiobrightness temperature, the values of the absorption coefficients of the solutions are calculated. The difference between the values of the solutions absorption coefficients and those for water is shown with a successive increase in their values with an increase in the degree of dilution. The difference between the values of the absorption coefficients of solutions of high dilutions and those for water is interpreted as a mechanism of cooperative interaction of ions with water molecules located outside the first hydration shell. A network of water molecules perturbed by ions contributes to the dielectric permittivity was obtained due to new properties at supramolecular interactions on the basis of dispersion interaction.

The exact similarity solutions (also called as special exact solutions) of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, no similarity solutions for two dimensional turbulent boundary layers have been reported in the literature. With the help of dimensional analysis and invariance principle, Prandtl mixing length $\ell=\kappa y$ for one dimensional turbulent boundary layer is extended to $\ell(x,y)=\kappa y (1-\frac{y}{\delta})\sqrt{\frac{\nu}{U\delta}} $ for the two dimensional turbulent boundary layers, furthermore, with a similarity transformation, we successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation $f'''+ ff''+\beta(1-f'^2)+\kappa^2[\eta^2(1-\eta)^2f''|f''|]'=0$. As an application, similarity solutions of the two dimensional turbulent boundary layer on a flat plate at zero incidence have been studied in detail. To solve the ordinary differential equation numerically, a complete Maple code is provided.

The exact similarity solutions of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, for two dimensional turbulent boundary layers, no similarity solutions (special exact solutions) have ever been found. In the light of Blasius' pioneer works, we extend Blasius similarity transformation to the two dimensional turbulent boundary layers, and successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation. By author's Maple code, we numerically solve the ordinary differential equation and produce some useful quantities.

The aim of the studies was to test the hydraulic pressure losses along the length, local pressure losses, absolute and relative roughness values of the pipeline walls made of modern polymer material. The article presents the results of the experimental studies of hydraulic resistance of pipelines from modern polymer materials, the method of theoretical determination of resistance coefficients in the VALTEC ISO 21003 metal polymer pipeline, the classification and method of determining the values of local hydraulic resistance coefficients of T-Junctions are given, the dependence for determining the Darcy λ hydraulic friction coefficient associated with the value of the equivalent hydraulic roughness coefficient of the inner surface of pipes and the Reynolds number. The experimental data of the hydraulic coefficient of friction Darcy λ was compared with the known theoretical relationships of A.D. Altschul, Colbrooke for smooth pipes, Colbrooke and White, an alternative to the Colbrooke-White, Offengenden equation for low-shear pipes and pipes from polymer materials with Reynolds numbers Re>10^4. Dependence of hydraulic coefficient of friction on relative equivalent roughness of pipeline walls and Reynolds number, specific energy loss per one linear meter on total flow rate in system of hydraulically short pipeline is obtained. The empirical dependence of values of coefficients of local resistance ζ=f(Q2/Q3) from a ratio of the divided expenses of Q2/Q3 of the equal-pass T-Junction α=90 the pipeline to a branch up, on pass in the direct T-Junction and division of streams at various internal pipeline diameters of VALTEC.

We report on the study of ultrafast laser-induced plasma expansion dynamics in a gas microjet. To this purpose, we focused femtosecond laser pulses on a nitrogen jet produced through a homemade De Laval micronozzle. The laser excitation leads to plasma excitation with a characteristic spectral line emission at 391 nm. By following the emitted signal with a detection system based on an Intensified Charge-Coupled Device (ICCD) we captured the two-dimensional spatial evolution of the photo-excited nitrogen ions with a temporal resolution on the nanosecond time scale. We fabricated the micronozzle on fused silica substrate by femtosecond laser micromachining. This technique enables high accuracy and three-dimensional capabilities, thus providing an ideal platform for developing glass-based microfluidic structures for application to plasma physics and ultrafast spectroscopy.

This study revisits the Reynolds-averaged Navier--Stokes equations (RANS) and finds that the existing literature is erroneous regarding the primary unknowns and the number of independent unknowns in the RANS. The literature claims that the Reynolds stress tensor has six independent unknowns, but in fact the six unknowns can be reduced to three that are functions of the three velocity fluctuation components, because the Reynolds stress tensor is simply an integration of a second-order dyadic tensor of flow velocity fluctuations rather than a general symmetric tensor. This difficult situation is resolved by returning to the time of Reynolds in 1895 and revisiting Reynolds' averaging formulation of turbulence. The study of turbulence modeling could focus on the velocity fluctuations instead of on the Reynolds stress. An advantage of modeling the velocity fluctuations is, from both physical and experimental perspectives, that the velocity fluctuation components are observable whereas the Reynolds stress tensor is not.

At INFN-LNS, and in collaboration with the ATOMKI laboratories, an innovative multi-diagnostic system with advanced analytical methods has been designed and implemented. This is based on several detectors and techniques (Optical Emission Spectroscopy, RF systems, Interfero-polarimetry, X-ray detectors) and here we focus on high resolution spatially-resolved X-ray spectroscopy, performed by means of a X-ray pin-hole camera setup operating in the 0.5−20 keV energy domain. The diagnostic system was installed at a 14 GHz Electron Cyclotron Resonance (ECR) ion source (ATOMKI, Debrecen), enabling high precision X-ray spectrally-resolved imaging of ECR plasmas heated by hundreds of Watts. The achieved spatial and energy resolution were 0.5 mm and 300 eV at 8 keV, respectively. We here present the innovative analysis algorithm that we properly developed for obtaining Single Photon-Counted (SPhC) images providing the local plasma emitted spectrum in a High-Dynamic-Range (HDR) mode, by distinguishing fluorescence lines of the materials of the plasma chamber (Ti, Ta) from plasma (Ar). This method allows a quantitative characterization of warm electrons population in the plasma (and its 2D distribution) which are the most important for ionization, and also to estimate local plasma density and spectral temperatures. The developed post-processing analysis is also able to remove the readout noise, that is often observable at very low exposure times (msec). The setup is now under update including fast shutters and trigger systems in order to allow simultaneously space and time-resolved plasma spectroscopy during transients, stable and turbulent regimes.

In large eddy simulation (LES) of turbulent flows, the most critical dynamical processes to be considered by dynamic subgrid models to account for an average cascade of kinetic energy from the largest to the smallest scales of the flow is not fully clear. Furthermore, evidence of vortex stretching being the primary mechanism of the cascade is not out of the question. In this article, we study some essential statistical characteristics of vortex stretching and its role in dynamic approaches of modeling subgrid-scale turbulence. We have compared the interaction of subgrid stresses with the filtered quantities among four models using invariants of the velocity gradient tensor. This technique is a single unified approach to studying a wide range of length scales in the turbulent flow. In addition, it also provides a rational basis for the statistical characteristics a subgrid model must serve in physical space to ensure an appropriate cascade of kinetic energy. Results indicate that the stretching mechanism extracts energy from the large-scale straining motion and passes it onto small-scale stretched vortices.

This paper considers the general ocean circulation within the thermodynamical closure of our climate theory, which aims to deduce the generic climate state from first principles. The preceding papers of the theory have reduced planetary fluids to warm/cold masses and determined their bulk thermal properties, which provide prior constraints for the derivation of the upper-bound circulation when the potential vorticity is homogenized in moving masses. In a companion paper on the atmosphere, this upper bound is seen to reproduce the prevailing wind, forsaking therefore previous discordant explanations of the easterly trade and the polar jet stream. In this paper on the ocean, we again show that this upper bound may replicate broad features of the observed circulation, including a western-intensified subtropical gyre and a counter-rotating tropical gyre feeding the equatorial undercurrent. Together, we posit that PV homogenization may provide a unifying dynamical principle of the large-scale planetary circulation, which may be interpreted as the maximum macroscopic motion extractable by microscopic stirring --- within the confine of the thermal differentiation.

The dynamics of a low-energy charged particle in an axis-symmetric magnetic field is known to be a regular superposition of periodic–although possibly incommensurate–motions. The projection of the particle orbit along the two non-ignorable coordinates (x,y) may be expressed in terms of each other: y=y(x), yet–to our knowledge–such a functional relation has never been directly produced in literature, but only by way of a detour: first, equations of motion are solved, yielding x=x(t),y=y(t), and then one of the two relations is inverted, x(t)→t(x). In this paper we present a closed-form functional relation which allows to express coordinates of the particle’ orbit without the need to pass through the hourly law of motion.

The current paper is devoted to the introduction of simpler Eulerian variational principles from which all the relevant equations of non-barotropic stationary magnetohydrodynamics can be derived for magnetic fields which lie on surfaces. A variational principle is given in terms of three independent variables for stationary non-barotropic magnetohydrodynamic flows. This is a smaller number of variables than the eight variables which appear in the standard equations of non-barotropic magnetohydrodynamics which are the magnetic field, the velocity field, the specific entropy and the density. We further investigate the case in which the flow along magnetic lines is not ideal.

The objective of this work is to identify the mechanism of dropwise condensation on a smooth solid surface. We investigate the stable dropwise condensation that occurs at a droplet growth rate of 1 µm/s in diameter on a gold-coated glass surface. Additionally, we present our observations on unstable dropwise condensation, i.e., degradation of dropwise condensation of steam on a gold surface. The Surface Plasmon Resonance Imaging (SPRi) approach used in this study has lateral resolutions of 4-10 µm, thickness resolutions of 0.1-1 nm, and temporal resolutions of 200-10,000 frames per second (FPS). SPRi is used to evaluate the existence and structure of thin films and occurrence of initial nuclei during condensation. Visualization of the onset of stable dropwise condensation suggests droplets form at heterogeneous nucleation sites and that no film greater than a monolayer exists on the surface before the formation of droplets. Observation of the unstable dropwise condensation of steam shows the existence of water films that are several nanometers thick between droplets. This work shows that neither the nucleation theory or film rupture theory can individually explain the physics of dropwise condensation. Therefore, there is a need for a more comprehensive theory that can explain the mechanism of dropwise condensation.

Natural selection is commonly seen not just as an explanation for adaptive evolution, but as the inevitable consequence of “heritable variation in fitness among individuals”. Although it remains embedded in biological concepts, such a formalisation makes it tempting to explore whether this precondition may be met not only in life as we know it, but also in other physical systems. This would imply that these systems are subject to natural selection and may perhaps be investigated in a biological framework, where properties are typically examined in light of their putative functions. Here we relate the major questions that were debated during a three-day workshop[1] devoted to discussing whether natural selection may take place in non-living physical systems. We start this report with a brief overview of research fields dealing with “life-like” or “proto-biotic” systems, where mimicking evolution by natural selection in test tubes stands as a major objective. We contend the challenge may be as much conceptual as technical. Taking the problem from a physical angle, we then discuss the framework of dissipative structures. Although life is viewed in this context as a particular case within a larger ensemble of physical phenomena, this approach does not provide general principles from which natural selection could be derived. Turning back to evolutionary biology, we ask to what extent the most general formulations of the necessary conditions or signatures of natural selection may be applicable beyond biology. In our view, such a cross-disciplinary jump is in large part impeded by reliance on individuality as a central yet implicit and loosely defined concept. Overall, these discussions thus lead us to conjecture that understanding, in physico-chemical terms, how individuality emerges and how it can be recognised, will be essential in the search for instances of evolution by natural selection outside of living systems. [1] Natural Selection Beyond Life? Observing the physico-chemical world through Darwinian glasses; 12-15 November 2019 - Institut d'Etudes Scientifiques, Cargèse, France

We present the results of direct numerical simulations of power spectral densities for kinetic energy, convective entropy and heat flux for unsteady Rayleigh-Bénard magnetoconvection in the frequency space. For larger values of frequency, the power spectral densities for all the global quantities vary with frequency f as f⁻². The scaling exponent is independent of Rayleigh number, Chandrasekhar’s number and thermal Prandtl number.

We consider the developed turbulence of capillary waves on shallow water. Analytic theory shows that an isotropic cascade spectrum is unstable with the respect to small angular perturbations, in particular, to spontaneous breakdown of the reflection symmetry and generation of nonzero momentum. By computer modeling we show that indeed a random pumping, generating on average zero momentum, produces turbulence with a nonzero total momentum. A strongly anisotropic large-scale pumping produces turbulence whose degree of anisotropy decreases along a cascade. It tends to saturation in the inertial interval and then further decreases in the dissipation interval. Surprisingly, neither the direction of the total momentum nor the direction of the compensated spectrum anisotropy is locked by our square box preferred directions (side or diagonal) but fluctuate.

The general coupling between particle transport and ionization-recombination processes in hot plasma is considered on the key concept of equilibrium charge state (CS) transport. A theoretical interpretation of particle and CS transport is gained in terms of a two-dimensional (2D) Markovian stochastic (random) processes, a discrete 2D Fokker-Plank-Kolmogorov equation (in charge and space variables) and generalized 2D coronal equilibrium between atomic processes and particle transport. The basic tool for analysis of CS equilibrium and transport is the equilibrium cell (EC) (two states on charge and two on space), which presents (i) a unit phase volume, (ii) the characteristic scale of local equilibrium, (iii) a comprehensive solution for the simplest nonlinear relations between transport and atomic processes. The approach opens up new perspectives on transport studies: (i) the direct modelling of equilibrium and transport of impurity using the atomic data base, (ii) recovery of the complete recombination rate profile based on knowledge of density profiles and ionization rate profiles, (iii) the local transport analysis, based on the reduction of the equilibrium set to the single EC (in particular, central or edge), (iv) analysis of the reduced transport coefficients (diffusion and convection) on the density profile measurements.

The nature of fluid mechanics makes experimentation an important part of a course taught in the subject. Presented here is the application of a novel, large scale multidisciplinary model of practical education in a fluids engineering laboratory. Advantages of this approach include efficiencies through economy of scale leading to better pedagogy for students. The scale justifies dedicated academic resources to focus on developing laboratory classes and giving specific attention to designing activities that meet learning outcomes. Four examples of applying this approach to fluids mechanics experiments are discussed, illustrating tactics that have been developed and honed through many repeated instances of delivery. “The measurement lab” uses a flow measurement context to teach identifying and managing general experimental uncertainty. New students, unfamiliar with fluid mechanics are guided through a process to gain understanding that can be applied to all future experimental activities. The “pressure loss in pipes” lab discusses the advantage and process for sharing equipment and teaching resources between multiple cohorts. The provision for students is adapted for context, such as the degree program or year of study. The “weirs big and small” lab provides a methodology to teach the power of dimensional analysis to mechanical engineers using a field of fluid mechanics that is outside their usual theoretical studies. Finally, the “spillway design” lab discusses mechanisms to deliver student independent, open ended experiments at scale, without excessive staff resource requirement.

Hasselmann and coauthors proposed the discrete interaction approximation (DIA) as the best tool replacing the nonlinear evolution term in a numerical wind-wave model. Much later, Polnikov and Farina radically improved the original DIA by means of location all the interacting four wave vectors, used in the DIA configuration, exactly at the nodes of the numerical frequency-angular grid. This provides nearly two-times enhancing the speed of numerical calculation for the nonlinear evolution term in a wind-wave model. For this reason, the proposed version of the DIA was called as the fast DIA (FDIA). In this paper we demonstrate all details of the FDIA concept for several frequency-angular numerical grids of high resolution, with the aim of active implementation the FDIA in modern versions of world-wide used wind-wave models.

This work investigates the utility of the National Oceanic and Atmospheric Administration’s National Water Model (NWM) for water management operations by assessing the total inflow into a select number of reservoirs across the Central and Western U.S. Total inflow is generally an unmeasured quantity, though critically important for anticipating both floods and shortages in supply over a short-term (hourly) to sub-seasonal (monthly) time horizon. The NWM offers such information at over 5,000 reservoirs across the U.S., however, its skill at representing inflow processes is largely unknown. The goal of this work is to understand the drivers for both well performing and poor performing NWM inflows such that managers can get a sense of the capability of NWM to capture natural hydrologic processes and in some cases, the effects of upstream management. We analyzed the inflows for a subset of Bureau of Reclamation (BoR) reservoirs within the NWM over the long-term simulations (retrospectively, seven years) and for short, medium and long-range operational forecast cycles over a one-year period. We utilize ancillary reservoir characteristics (e.g. physical and operational) to explain variation in inflow performance across the selected reservoirs. In general, we find that NWM inflows in snow-driven basins outperform those in rain-driven, and that assimilated basin area, upstream management, and calibrated basin area all influence the NWM’s ability to reproduce daily reservoir inflows. The final outcome of this work proposes a framework for how the NWM reservoir inflows can be useful for reservoir management, linking reservoir purposes with the forecast cycles and retrospective simulations.

A solution for the Prandtl-Blasius equation is essential to all kinds of boundary layer problems. This paper revisits this classic problem and presents a general Maple code as its numerical solution. The solutions were obtained from the Maple code, using the Runge-Kutta method. The study also considers convergence radius expanding and an approximate analytic solution is proposed by curve fitting. Similarly, the study resolves some boundary layer related problems and provide relevant Maple codes for these.

Metal foams are widely studied as possible tools for the enhancement of heat transfer from hot bodies. The basic idea is that a metal foam tends to increase significantly the heat exchange area between the hot solid body and the external cooling fluid. For this reason, this class of porous materials are considered as good candidates for an alternative to finned surfaces, with different pros and cons. Among the pros, we mention the generally wider area of contact between solid and fluid. Among the cons is the difficulty of producing different specimens with the same inner structure, with the consequence that their performance may be significantly variable. This paper will offer a survey of the literature with a focus on the main heat transfer characteristics of the metal foams. Then, a numerical simulation of the heat transfer at the pore scale level for an artificial foam with a spatially periodic structure will be discussed. Finally, these numerical results will be employed to assess the macroscopic modelling of the flow and heat transfer in a metal foam.

The simultaneous evaluation of heat/mass transference procures momentous significance in the polymer industry and engineering activities. Extensive utilization of such aspects includes propylene flares, heat exchangers, energy transfiguration in chilling towers, and amputation of post fortuitous heat in nuclear reactors. Owing to such prospective demands, we interpreted the Jeffrey liquid reactive flow under non-linear convection. We scrutinized the transference of heat/mass under generalized Fourier-Fick relations. Thermal conductivity depends on temperature while mass diffusivity is dependent on concentration. Besides, heat source along with first-order chemical reaction aspects are accounted for. Relevant transformations are exerted to achieve non-linear differential systems that are solved through the homotopy scheme. Influences of divers factors are exhibited via graphical benchmark.

For the solution existence condition of the Navier-Stokes equation, we propose a conjecture as follows: "\emph{The Navier-Stokes equation has a solution if and only if the determinant of flow velocity gradient is not zero, namely $\det (\bm \nabla \bm v)\neq 0$.}"

The material derivative is important in continuum physics. This Letter shows that the expression $\frac{d }{dt}=\frac{\partial }{\partial t}+(\bm v\cdot \bm \nabla)$, used in most literature and textbooks, is incorrect. The correct expression $ \frac{d (:)}{dt}=\frac{\partial }{\partial t}(:)+\bm v\cdot [\bm \nabla (:)]$ is formulated. The solution existence condition of Navier-Stokes equation has been proposed from its form-solution, the conclusion is that "\emph{The Navier-Stokes equation has a solution if and only if the determinant of flow velocity gradient is not zero, namely $\det (\bm \nabla \bm v)\neq 0$.}"

Previously, a conductivity invariance phenomena (CIP) has been discovered – at a certain lift-off, the inductance change of the sensor due to a test sample is immune to conductivity variations, i.e. the inductance – lift-off curve passes through a common point at a certain lift-off, termed as conductivity invariance lift-off. However, this conductivity invariance lift-off is fixed for a particular sensor setup, which is not convenient for various sample conditions. In this paper, we propose using two parameters in the coil design – the horizontal and vertical distances between the transmitter and the receiver to control the conductivity invariance lift-off. The relationship between these two parameters and the conductivity invariance lift-off is investigated by simulation and experiments and it has been found that there is an approximate linear relationship between these two parameters and the conductivity invariance lift-off. This is useful for applications where the measurements have restrictions on lift-off, e.g. uneven coating thickness which limits the range of the lift-off of probe during the measurements. Therefore, based on this relationship, it can be easier to adjust the configuration of the probe for a better inspection of the test samples.

This paper presents an investigation into the influence of drag-reducing polymer (DRP) on the flow of water inside and around U-bends. The DRPs tested were partially hydrolysed polyacrylamide (HPAM) and two different molecular weight polyethylene oxides (PEO). The Reynolds number and in-situ concentration ranges of water and DRP were and respectively. The test section was made up of clear polyvinyl chloride (CPVC) pipes that contained two U-bends with radii of curvature (R) of 100 and 200 mm as well as straight sections of 19-mm ID. Higher pressure gradient was recorded in the under-developed flow section after the bend. Drag reduction (DR) increased with in-situ concentration of the DRP and the flowrate up to some threshold values. The concentration and flowrate thresholds were about 20 ppm and 1.07 m/s respectively. DR differed among the different sections considered with the highest values recorded for developed flows upstream of the bend and the least values recorded in the bend. DR at station 1 (immediately after the bend) was generally lower than that at station 2 (further downstream). P-K Plots showed that, the DR envelop for developed flow in straight pipes did not apply to the bend and the under-developed flow section. DR increased with polymer molecular weight and ionic strength in straight sections. However, molecular weight (and polymer shearing) had limited effect on DR in the bends. DR characteristics at the bend suggested that there was a different or additional mechanism. The effect of bend curvature ratio on DR was predominant after the U-bend. The larger the ratio of pipe diameter to bend diameter, the lower the DR in the redeveloping sections after the bend. For R = 100 mm, the highest DR recorded for fully developed flows before bend, in the bend, at stations 1 and 2 were 57, 31, 36 and 33% respectively while that for R = 200 mm were 64, 28, 42 and 49% respectively.

Pressure drop and flow pattern of oil-water flows were investigated in a 19 mm ID clear polyvinyl chloride pipe consisting of U-bend with radius of curvature of 100 mm. The range for oil and water superficial velocities tested were and respectively. Measurements were carried out under different flow conditions in a test section that consisted of four different parts: upstream of the bend, at the bend and at two redeveloping flow locations after the bend. The result indicated that the bend had limited influence on downstream flow patterns. However, the shear forces imposed by the bend caused some shift flow pattern transition and bubble characteristics in the redeveloping flow section after the bend relative to develop flow before the bend. Generally, pressure gradient at all the test sections increased with both oil fraction and water superficial velocity and there was a sharp change of pressure gradient profile during phase inversion. The transition point where phase inversion occurred was always within the range of . Pressure losses differed at the various test sections and the difference was strongly linked to the superficial velocity of the phases and the flow pattern. At high mixture velocity, pressure losses at the redeveloping section after the bend were higher than that at the bend and that for fully developed flows. At low mixture velocity, pressure losses at the bend are higher than in the straight sections. Pressure drop generally decreased with level of flow development downstream of the bend.

In 1895, Korteweg and de Vries (KdV), derived their celebrated equation describing the motion of waves of long wavelength in shallow water. In doing so they made a number of quite reasonable assumptions, incompressibility of the water and irrotational fluid. The resulting equation, the celebrated KdV equation, has been shown to be a very reasonable description of real water waves. However there are other phenomena which have an impact on the shape of the wave, that of vorticity and viscosity. This paper examines how a constant vorticity affects the shape of waves in electrohydrodynamics. For constant vorticity, the vertical component of the velocity obeys a Laplace equation and also has the usual lower boundary condition. In making the vertical component of the velocity take central stage, the Burns condition can be thus bypassed.

Based on recent computational and experimental studies, hesperidin, a bioactive flavonoid abundant in citrus peel, stands out for its high binding affinity to the main cellular receptors of SARS-CoV-2, outperforming drugs already recommended for clinical trials. Thus, it is very promising for prophylaxis and treatment of COVID-19, along with other coexistent flavonoids such as naringin, which could help restraining the pro-inflammatory overreaction of the immune system. Controlled hydrodynamic cavitation processes showed the highest speed, effectiveness and efficiency in the integral and green aqueous extraction of flavonoids, essential oils and pectin from citrus peel waste. After freeze-drying, the extracted pectin showed high quality and excellent antioxidant and antibacterial activities, attributed to flavonoids and essential oils adsorbed and concentrated on its surface. This study reviews the recent evidence about hesperidin as a promising molecule, and proposes a feasible and affordable process based on hydrodynamic cavitation for the integral aqueous extraction of citrus peel waste resulting in hesperidin-rich products, either aqueous extracts or pectin tablets. The uptake of this process on a relevant scale is urged, in order to achieve large-scale production and distribution of hesperidin-rich products. Meanwhile, experimental and clinical studies could determine the effective doses either for therapeutic and preventive purposes.

The need to understand the process by which particles, including solar wind and coronal ions as well as pickup ions, are accelerated to high energies (ultimately to become anomalous cosmic rays) motivate a multi-fluid shock wave model which includes kinetic effects (e.g. ion acceleration) in an electromagnetically self-consistent framework. Particle reflection at the cross-shock potential leads to ion acceleration in the motional electric field and thus anisotropic heating and pressure in the shock layer, with important consequences for the multi-fluid dynamics. This motivates development of a multi-fluid model of solar wind electrons and ions treated as fluid, coupled self-consistently with a small population of kinetically treated ions (e.g. pickup ions.) Consideration of both the time dependent and steady state regimes, indicate that such a multi-fluid approach is necessary for resolving the, Debye scale, particle reflecting cross-shock potential and subsequent dynamics. To study charge separation effects in narrow, supersonic wave layers we consider a reduction of the system to the steady state for cold ions and hot electrons and find two types of solitary waves inherent to the reduced two-fluid system in this limiting case.

Beam dumps are indispensable components for particle accelerator facilities to absorb or dispose beam kinetic energy in a safe way. However, the design of beam dumps based on conventional technology, i.e. the energy deposition via beam-dense matter interaction, makes the beam dump facility complicated and large in size, partly due to nowadays’ high beam intensities and energies achieved. In addition, these high-power beams generate radioactive hazards, which need specific methods to deal with. On the other hand, the EuPRAXIA project can advance the laser-plasma accelerator significantly by achieving 1-5 GeV high quality electron beam in a compact layout. Nevertheless, the beam dump based on conventional technique will still produce radiation hazards and make the overall footprint less compact. Here, we propose to implement a plasma beam dump to absorb the kinetic energy from the EuPRAXIA beam. In doing so, the overall compactness of the EuPRAXIA layout will not be impacted, and the radioactivity generated by the facility can be mitigated. In this paper, results from particle-in-cell (PIC) simulations are presented for plasma beam dumps based on EuPRAXIA beam parameters.

Medical institutions, where several patients are treated and medical workers engaged, are always exposed to secondary viral and bacterial infections. It is critical to prevent infection transmission by indirect as well as direct contact through air or splash. The infections of most diseases can be transmitted through the air. HEPA filters installed in air conditioning equipment are used to prevent infection transmission through air in medical institutions, but air circulation takes a long time in a large space. Virus and bacteria smaller than 0.3 μm cannot be removed by the HEPA filter; hence, those microbes remain alive throughout the air ventilation. A plasma sterilizer has the capability to provide environmental friendly sterilization by employing reactive oxide species and reactive nitrogen species at a low cost. We developed an excellent plasma sterilizer by using a non-thermal atmospheric-pressure biocompatible plasma (NBP). Ozone concentration in plasma sources has been derived by Kuhn et al. [1]. The diffusion coefficients inside (D₀) and outside (D₁) the plasma sterilizer have been calculated to be 0.0641 m² s^-1 and 0.717 m² s^-1, respectively. To sustain high O₃ concentrations over 121 ppm inside the plasma source and low O₃ concentrations below 0.05 ppm outside the sterilizer, it is necessary to keep O₃ concentrations at the exit of plasma sterilizer below 0.28 ppm. so that diffusion coefficient D₁ has been designed to be as large as 11 times of D₀.

With improved measurement and modelling technology, variability has emerged as an essential feature in non-equilibrium processes. While traditionally, mean values and variance have been heavily used, they are not appropriate in describing extreme events where a significant deviation from mean values often occurs. Furthermore, stationary Probability Density Functions (PDFs) miss crucial information about the dynamics associated with variability. It is thus critical to go beyond a traditional approach and deal with time-dependent PDFs. Here, we consider atmospheric data from the Whole Atmosphere Community Climate Model (WACCM) and calculate time-dependent PDFs and the information length from these PDFs, which is the total number of statistically different states that a system passes through in time. Time-dependent PDFs are shown to be non-Gaussian in general, and the information length calculated from these PDFs shed us a new perspective of understanding variabilities, correlation among different variables and regions. Specifically, we calculate time-dependent PDFs and information length and show that the information length tends to increase with the altitude albeit in a complex form. This tendency is more robust for flows/shears than temperature. Also, much similarity among flows and shears in the information length is found in comparison with the temperature. This means a stronger correlation among flows/shears because of a strong coupling through gravity waves in this particular WACCM model. We also find the increase of the information length with the latitude and interesting hemispheric asymmetry for flows/shears/temperature, a stronger anti-correlation (correlation) between flows/shears and temperature at a higher (low) latitude. These results also suggest the importance of high latitude/altitude in the information budge in the Earth's atmosphere, the spatial gradient of the information as a useful proxy for the transport of physical quantities.

In this paper we develop from first principles a unique law pertaining to the flow of fluids through closed conduits. This law, which we call “Quinn’s Law”, may be described as follows: When fluids are forced to flow through closed conduits under the driving force of a pressure gradient, there is a linear relationship between the normalized dimensionless pressure gradient, P_Q, and the normalized dimensionless fluid current, C_Q. The relationship is expressed mathematically as: P_Q = k₁ +k₂C_Q. This linear relationship remains the same whether the conduit is filled with or devoid of solid obstacles. The law differentiates, however, between a packed and an empty conduit by virtue of the tortuosity of the fluid path, which is seamlessly accommodated within the normalization framework of the law itself. When movement of the fluid is very close to being at rest, i.e., very slow, this relationship has the unique minimum constant value of k₁, and as the fluid acceleration increases, it varies with a slope of k₂ as a function of normalized fluid current. Quinn’s Law is validated herein by applying it to the data from published classical studies of measured permeability in both packed and empty conduits, as well as to the data generated by home grown experiments performed in the author’s own laboratory.