Modeling of Cavitation Bubble Cloud with Discrete Lagrangian Tracking

In this paper, a Lagrangian-Eulerian (LE) two-way coupling model is developed to numerically study the cavitation bubble cloud. In this model, the gas-liquid mixture is treated directly as a continuous and compressible fluid and the governing equations are solved by methods in Eulerian descriptions. An isobaric closure exhibiting better consistency properties is applied to evaluate the pressure of gas-liquid mixture. The dispersed gas/vapor bubbles are tracked in a Lagrangian fashion, and their compression and expansion are described by a modified Rayleigh-Plesset equation, which considers the close-by flow properties other than these of the infinity for each bubble. The performance of the present method is validated by a number of benchmark tests. Then, this model is applied to study how the bubble cloud affects the shape and propagation of a pressure wave when the pressure pulse travels through. In the end, a three-dimensional simulation of a vapor cloud’s Rayleigh collapse is carried out, and the induced extreme pressure is discussed in detail. The total bubble number’s influence on the extreme collapse pressure and the size distribution of bubbles during the collapse are also analyzed.

Lagrangian eddy tracking reveals the Eratosthenes anticyclonic attractor in the eastern Levantine Basin

Statistics of anticyclonic eddy activity and eddy trajectories in the Levantine Basin over the 2000–2018 period are analyzed using the DYNED-Atlas database, which links automated mesoscale eddy detection by the Angular Momentum Eddy Detection and Tracking Algorithm (AMEDA) algorithm to in situ oceanographic observations. This easternmost region of the Mediterranean Sea, delimited by the Levantine coast and Cyprus, has a complex eddying activity, which has not yet been fully characterized. In this paper, we use Lagrangian tracking to investigate the eddy fluxes and interactions between different subregions in this area. The anticyclonic structure above the Eratosthenes Seamount is identified as hosting an anticyclone attractor, constituted by a succession of long-lived anticyclones. It has a larger radius and is more persistent (staying in the same position for up to 4 years with successive merging events) than other eddies in this region. Quantification of anticyclone flux shows that anticyclones that drift towards the Eratosthenes Seamount are mainly formed along the Israeli coast or in a neighboring area west of the seamount. The southeastern Levantine area is isolated, with no anticyclone transfers to or from the western part of the basin, defining the effective attraction basin for the Eratosthenes anticyclone attractor. Co-localized in situ profiles inside eddies provide quantitative information on their subsurface physical anomaly signature, whose intensity can vary greatly with respect to the dynamical surface signature intensity. Despite interannual variability, the so-called Eratosthenes anticyclone attractor stores a larger amount of heat and salt than neighboring anticyclones, in a deeper subsurface anomaly that usually extends down to 500 m. This suggests that this attractor could concentrate heat and salt from this subbasin, which will impact the properties of intermediate water masses created there.

‪Impulse wave generation: a comparison of landslides of block and granular masses by coupled Lagrangian tracking using VOF over a set mesh

This paper presents the numerical results of impulsive waves generated by landslides of solid block, granular materials and heavy block sinking. An impulse product parameter P is developed and a wide range of effective parameters are studied. The volume-of-fluid (VOF) and overset mesh methods have been used to study landslide-generated tsunamis. Also, a Lagrangian tracking approach coupled with the VOF to simulate the granular movement. The effect of the water reservoir depth, the landslide height, the landslide density and the geometrical parameters on the wave height (elevation) has been investigated using the open-source OpenFOAM software. The results have been presented for dimensionless distances and the normalized geometry of the landslide in the range 5–7, 1–2, respectively. These numbers have been normalized the height of the landslide (a). According to the results of simulations, the tsunami formation process is divided into three stages, which were analyzed in detail by considering the interactions between the solid and the water reservoir. The Scott Russell wave has the highest impulse product parameter among the impulse wave mechanisms which is 58.6% of the total impulse production. In addition, the duration of the wave propagation has been computed based on the wave height.

Deposition Characteristics of Firebrands on and Around Rectangular Cubic Structures

The focus of the present work is on the deposition of firebrands in a flow over a rectangular cubic block representative of a structure in wildland-urban interface (WUI). The study was carried out by physics based modeling where the wind flow turbulence was dealt with by large eddy simulation (LES) and firebrands were treated by Lagrangian tracking. The Lagrangian equations coupled with the flow solver, accounted for both translational and rotational motions as well as thermochemical degradation of firebrands, assumed to be cylindrical. The dimensions of the structure were varied from 3 to 9 m in the simulations for a parametric study. The simulations were carried out by tracking many firebrands randomly released with a uniform distribution from a horizontal plane 35 m above the ground into the computational domain. The coordinates of the deposited firebrands were used to calculate their normalized number density (number of landed firebrands per unit surface area) to quantify their deposition pattern. On the leewardside of the block, an area, referred to as the safe zone, was identified right behind the structure where firebrands never deposit. The size of the safe zone in the direction perpendicular to the wind was nearly identical to the width of the structure. The length of the safe zone in the wind direction was proportional to the height of the structure. The leeward face of the blocks was never hit by a firebrand. The windward face was hit by many more firebrands than the lateral faces but much less than the top face. The distribution of the number density of the deposited firebrands on the top face was found to be correlated with the flow separation and reattachment on this face.

Future changes in European windstorm severities and impacts

<p>Extratropical cyclones that impact Europe are significant natural hazards that can cause severe damages and lead to large socio-economic losses. There are large uncertainties associated with changes in storm number, and storm intensity, in future climate projections. Here we use a Lagrangian tracking algorithm applied to reanalysis data and to historical and two future scenario simulations of a number of CMIP6 models to investigate future changes in characteristics of windstorms over Europe. As well as storm frequencies and peak wind speeds, we also quantify changes in two versions of a storm severity index (SSI), one of which is population weighted. These metrics are calculated using the footprints of cyclones as they pass over the European continent.</p><p>The models show differing abilities to represent the historical SSI compared to ERA5. Future changes in SSI are somewhat uncertain, but tend to show an increase in severities over central and northwest Europe, and a decrease over lower latitudes and the Mediterranean, with responses tending to be larger for the more extreme climate change scenario. These changes are associated with the changes in the extreme winds over land.</p><p>By considering the parameters of population density and wind intensity threshold we could explore the relevance of future socio-economic and adaptive changes on the windstorm impacts. For the population-weighted SSI, smaller increases are found in the future cases where population densities do not change and/or adaptation to increases in extreme wind speed climatologies occur.</p>

Suppression of Baroclinic Eddies by Strong Jets

<p>The midlatitude storm tracks are one of the most prominent features of extratropical climate. Despite the theoretical expectation, based on baroclinic instability theory that baroclinic eddy strength correlates with jet intensity, there is a decrease in storm-track activity during midwinter over the Pacific compared to the shoulder seasons. Recent studies suggest this phenomenon is a result of the general circulation effect on the storm-track through interaction with the jet-stream. To isolate the effect of jet strength, we conduct a series of GCM experiments with a systematically varied jet intensity. The simulations are analyzed using Lagrangian tracking to understand the response from a single eddy perspective. The results of the Lagrangian analysis show that while the response of upper-level eddies is dominated by a reduction in the amount of tracked features, the lower-level eddies' response is also affected by a reduction in their lifetime. Analyzing the effect of the jet strength on the pairing between the upper- and lower-level eddies, we show how the jet intensification break the baroclinic wave structure and limits its growth. Furthermore, we show that these results can be settled with linear baroclinic instability models if the eddies' spatial scale is considered. The intensification of the jet and increase in the deformation radius shift the preferred scale for growth from the synoptic-scale toward the planetary-scale, consistent with the reduction in storm activity. This mechanism potentially explains the midwinter suppression of storm activity over the Pacific and the difference from the response over the Atlantic.</p>

Lagrangian eddy tracking reveals the Eratosthenes anticyclonicattractor in the eastern Levantine basin

<p>Statistics of anticyclone activity and trajectories in the southeastern Mediterranean sea over the period 2000-2018<br>is created using the DYNED atlas, which links the automated mesoscale eddy detection by the AMEDA algorithm with in<br>situ oceanographic observations. This easternmost region of the Mediterranean sea, delimited by the Levantine coast and<br>Cyprus, has a complex eddying activity, which has not yet been fully characterized. Using Lagrangian tracking<br>to investigate the eddy fluxes and interactions between different subregions in this area, we find that the southeastern Levantine<br>area is isolated, with no anticyclone exchanges with the western part of the basin. Moreover the anticyclonic structure above<br>the Eratosthenes seamount is identified as being an anticyclone attractor, differentiated from other anticyclones and staying<br>around this preferred position up to four years with successive mergings. Colocalized in situ profiles inside eddies provide<br>quantitative information on their subsurface structure and show that similar surface signatures correspond to very different<br>physical properties. Despite interannual variability, the so-called "Eratosthenes attractor" stores a larger amount of heat and<br>salt than neighboring anticyclones, in a deeper subsurface anomaly that usually extend down to 500 m. This suggests that this<br>attractor could concentrate heat and salt from this sub-basin, which will impact the properties of intermediate water masses<br>created there.</p>

Lagrangian eddy tracking reveals the Eratosthenes anticyclonic attractor in the eastern Levantine basin

Statistics of anticyclone activity and trajectories in the southeastern Mediterranean sea over the period 2000–2018 is created using the DYNED atlas, which links the automated mesoscale eddy detection by the AMEDA algorithm with in situ oceanographic observations. This easternmost region of the Mediterranean sea, delimited by the Levantine coast and Cyprus, has a complex eddying activity, which has not yet been fully characterized. In this paper we use Lagrangian tracking to investigate the eddy fluxes and interactions between different subregions in this area. We find that the southeastern Levantine area is isolated, with no anticyclone exchanges with the western part of the basin. Moreover the anticyclonic structure above the Eratosthenes seamount is identified as being an anticyclone attractor, differentiated from other anticyclones and staying around this preferred position up to four years with successive mergings. Colocalized in situ profiles inside eddies provide quantitative information on their subsurface structure and show that similar surface signatures correspond to very different physical properties. Despite interannual variability, the so-called Eratosthenes attractor stores a larger amount of heat and salt than neighbouring anticyclones, in a deeper subsurface anomaly that usually extend down to 500 m. This suggests that this attractor could concentrate heat and salt from this sub-basin, which will impact the properties of intermediate water masses created there.

Computational Protocol for Spray Flow Simulations Including Primary Atomization

Primary atomization is the key element in spray flow simulations. We have, in our previous work, used and validated the integral form of the conservation equations, leading to the “quadratic formula” for determination of the drop size during spray atomization in various geometry. A computational protocol has been developed where this formulation is adapted to existing computational frameworks for continuous and dispersed (droplet) liquid phase, for simulations of pressure-atomized sprays with and without swirl. In principle, this protocol can be applied to any spray geometry, with appropriate modifications in the atomization criterion. The preatomization continuous liquid motion (e.g., liquid column or sheet) is computed using volume-of-fluid (VOF) or similar methods, then the velocity data from this computation is input to the quadratic formula for determination of the local drop size. This initial drop size, along with the local liquid velocities from VOF, is then used in a Lagrangian tracking algorithm for the postatomization dispersed droplet calculations. This protocol can be implemented on coarse-grid, time-averaged simulations of spray flows, and produces convincing results when compared with experimental data for pressure-atomized sprays with and without swirl. This approach is general, and can be adapted in any spray geometries for complete and efficient computations of spray flows.