2024-08-15 | | Total: 6

The paper aims at assessing a hypothesis that resolution required to evaluate fuel consumption and heat release rates by directly (i.e., without a subgrid model of unresolved influence of small-scale turbulent eddies on the local flame) processing filtered fields of density, temperature, and species mass fractions should be significantly finer than resolution required to directly compute flame surface density by processing the same filtered fields. For this purpose, box filters of various widths are applied to three-dimensional Direct Numerical Simulation data obtained earlier by Dave et al. (Combust. Flame 196 (2018) 386-399) from a statistically one-dimensional and planar, moderately lean H2/air complex-chemistry flame propagating in a box under conditions of sufficiently intense small-scale turbulence (Karlovitz number is larger than unity and a ratio of laminar flame thickness to Kolmogorov length scale is about 20). Results confirm this hypothesis and show that the mean flame surface density and area can be predicted with acceptable accuracy by processing filtered combustion progress variable fields computed using a sufficiently wide filter. Such an approach does not require a model of the influence of subgrid turbulent eddies on flame surface density provided that filter width and laminar flame thickness are of the same order of magnitude. Good performance of this approach is attributed to inability of small-scale (when compared to the thickness) turbulent eddies to substantially change the local flame structure, which, nevertheless, is significantly perturbed by larger turbulent eddies that strain the local flame.

Barchans are eolian dunes of crescent shape found on Earth, Mars and other celestial bodies. Among the different types of barchan-barchan interaction, there is one, known as chasing, in which the dunes remain close but without touching each other. In this paper, we investigate the origins of this barchan-barchan dune repulsion by carrying out grain-scale numerical computations in which a pair of granular heaps is deformed by the fluid flow into barchan dunes that interact with each other. In our simulations, data such as position, velocity and resultant force are computed for each individual particle at each time step, allowing us to measure details of both the fluid and grains that explain the repulsion. We show the trajectories of grains, time-average resultant forces, and mass balances for each dune, and that the downstream barchan shrinks faster than the upstream one, keeping, thus, a relatively high velocity although in the wake of the upstream barchan. In its turn, this fast shrinkage is caused by the flow disturbance, which induces higher erosion on the downstream barchan and its circumvention by grains leaving the upstream dune. Our results help explaining the mechanisms behind the distribution of barchans in dune fields found on Earth and Mars.

In this work, we investigate the mixing of active scalars in two dimensions by the stirring action of stochastically generated shock waves. We use direct numerical simulations (DNS) of the interaction of shock waves with two non-reacting species to analyse the mixing dynamics for different Atwood numbers (At). Unlike passive scalars, the presence of density gradients in active scalars makes the species diffusion nonlinear, introducing a concentration gradient-driven term and a density gradient-driven nonlinear dissipation term in the concentration evolution equation. We show that the direction of the concentration gradient causes the interface across which molecular diffusion occurs to expand outward or inward, even without any stirring action. Shock waves enhance the mixing process by increasing the perimeter of the interface and by sustaining concentration gradients. Negative Atwood number mixtures sustain concentration gradients for longer time than positive Atwood number mixtures due to the so-called nonlinear dissipation terms. We estimate the time till when the action of stirring is dominant over molecular mixing. We also highlight the role of baroclinicity in increasing the interface perimeter in the stirring dominant regime. We compare the stirring effect of shock waves on mixing of passive scalars with active scalars and show that the vorticity generated by baroclinicity is responsible for the folding and stretching of the interface in the case of active scalars. We conclude by showing that lighter mixtures with denser inhomogeneities (At < 0) take longer time to homogenise than the denser mixtures with lighter inhomogeneities (At > 0).

We present a novel hybrid incompressible flow/material point method solver for simulating the combustion of flammable solids. Our approach utilizes a sparse grid representation of solid materials in the material point method portion of the solver and a hybrid Eulerian/FLIP solver for the incompressible portion. We utilize these components to simulate the motion of heated air and particulate matter as they interact with flammable solids, causing combustion-related damage. We include a novel particle sampling strategy to increase Eulerian flow accuracy near regions of high temperature. We also support control of the flame front propagation speed and the rate of solid combustion in an artistically directable manner. Solid combustion is modeled with temperature-dependent elastoplastic constitutive modeling. We demonstrate the efficacy of our method on various real-world three-dimensional problems, including a burning match, incense sticks, and a wood log in a fireplace.

A Stokes wave is a traveling free-surface periodic water wave that is constant in the direction transverse to the direction of propagation. In 1981 McLean discovered via numerical methods that Stokes waves are unstable with respect to transverse perturbations. In \cite{CreNguStr} for the case of infinite depth we proved rigorously that the spectrum of the water wave system linearized at small Stokes waves, with respect to transverse perturbations, contains unstable eigenvalues lying approximately on an ellipse. In this paper we consider the case of finite depth and prove that the same spectral instability result holds for all but finitely many values of the depth. The computations and some aspects of the theory are considerably more complicated in the finite depth case.

Hemodynamic quantities are valuable biomedical risk factors for cardiovascular pathology such as atherosclerosis. Non-invasive, in-vivo measurement of these quantities can only be performed using a select number of modalities that are not widely available, such as 4D flow magnetic resonance imaging (MRI). In this work, we create a surrogate model for hemodynamic flow field estimation, powered by machine learning. We train graph neural networks that include priors about the underlying symmetries and physics, limiting the amount of data required for training. This allows us to train the model using moderately-sized, in-vivo 4D flow MRI datasets, instead of large in-silico datasets obtained by computational fluid dynamics (CFD), as is the current standard. We create an efficient, equivariant neural network by combining the popular PointNet++ architecture with group-steerable layers. To incorporate the physics-informed priors, we derive an efficient discretisation scheme for the involved differential operators. We perform extensive experiments in carotid arteries and show that our model can accurately estimate low-noise hemodynamic flow fields in the carotid artery. Moreover, we show how the learned relation between geometry and hemodynamic quantities transfers to 3D vascular models obtained using a different imaging modality than the training data. This shows that physics-informed graph neural networks can be trained using 4D flow MRI data to estimate blood flow in unseen carotid artery geometries.