No Arabic abstract
In the present study, detailed numerical simulations are performed to investigate the primary breakup of a gasoline surrogate jet under non-evaporative Spray G operating conditions. The Spray G injector and operating conditions, developed by the Engine Combustion Network (ECN), represent the early phase of spray-guided gasoline injection. To focus the computational resources on resolving the primary breakup, simplifications have been made on the injector geometry. The effect of the internal flow on the primary breakup is modeled by specifying a nonzero injection angle at the inlet. The nonzero injection angle results in an increase of the jet penetration speed and also a deflection of the liquid jet. A parametric study on the injection angle is performed, and the numerical results are compared to the experimental data to identify the injection angle that best represents the Spray G conditions. The nonzero injection angle introduces an azimuthally non-uniform velocity in the liquid jet, which in turn influences the instability development on the jet surfaces and also the deformation and breakup of the jet head. The asymmetric primary breakup dynamics eventually lead to an azimuthal variation of droplet size distributions. The number of droplets varies significantly with the azimuthal angle, but interestingly, the probability density functions (PDF) of droplet size for different azimuthal angles collapse to a self-similar profile. Analysis has also been conducted to estimate the percentage and statistics of the tiny droplets that are under resolved in the present simulation. The PDF of the azimuthal angle is also presented, which is also shown to exhibit a self-similar form that varies little over time. Finally, a model is developed to predict the droplet number as a function of droplet diameter, azimuthal angle where a droplet is located, and time.
This article presents an original methodology for the prediction of steady turbulent aerodynamic fields. Due to the important computational cost of high-fidelity aerodynamic simulations, a surrogate model is employed to cope with the significant variations of several inflow conditions. Specifically, the Local Decomposition Method presented in this paper has been derived to capture nonlinear behaviors resulting from the presence of continuous and discontinuous signals. A combination of unsupervised and supervised learning algorithms is coupled with a physical criterion. It decomposes automatically the input parameter space, from a limited number of high-fidelity simulations, into subspaces. These latter correspond to different flow regimes. A measure of entropy identifies the subspace with the expected strongest non-linear behavior allowing to perform an active resampling on this low-dimensional structure. Local reduced-order models are built on each subspace using Proper Orthogonal Decomposition coupled with a multivariate interpolation tool. The methodology is assessed on the turbulent two-dimensional flow around the RAE2822 transonic airfoil. It exhibits a significant improvement in term of prediction accuracy for the Local Decomposition Method compared with the classical method of surrogate modeling for cases with different flow regimes.
Despite their fundamental and applied importance, a general model to predict the natural breakup length of steady capillary jets has not been proposed yet. In this work, we derive a scaling law with two universal constants to calculate that length as a function of the liquid properties and operating conditions. These constants are determined by fitting the scaling law to a large set of experimental and numerical measurements, including previously published data. Both the experimental and numerical jet lengths conform remarkably well to the proposed scaling law. This law is explained in terms of the growth of perturbations excited by the jet breakup itself.
Turbulent flows under transcritical conditions are present in regenerative cooling systems of rocker engines and extraction processes in chemical engineering. The turbulent flows and the corresponding heat transfer phenomena in these complex processes are still not well understood experimentally and numerically. The objective of this work is to investigate the turbulent flows under transcritical conditions using DNS of turbulent channel flows. A fully compressible solver is used in conjunction with a Peng-Robinson real-fluid equation of state to describe the transcritical flows. A channel flow with two isothermal walls is simulated with one heated and one cooled boundary layers. The grid resolution adopted in this study is slightly finer than that required for DNS of incompressible channel flows. The simulations are conducted using both fully (FC) and quasi-conservative (QC) schemes to assess their performance for transcritical wall-bounded flows. The instantaneous flows and the statistics are analyzed and compared with the canonical theories. It is found that results from both FC and QC schemes qualitatively agree well with noticeable difference near the top heated wall, where spurious oscillations in velocity can be observed. Using the DNS data, we then examine the usefulness of Townsend attached eddy hypothesis in the context of flows at transcritical conditions. It is shown that the streamwise energy spectrum exhibits the inverse wavenumber scaling and that the streamwise velocity structure function follows a logarithmic scaling, thus providing support to the attached eddy model at transcritical conditions.
Cardiovascular diseases, specifically cerebral aneurysms, represent a major cause of morbidity and mortality, having a significant impact on the cost and overall status of health care. In the present work, we employ a haemorheological blood model originally proposed by Owens to investigate the haemodynamics of blood flow through an aneurytic channel. This constitutive equation for whole human blood is derived using ideas drawn from temporary polymer network theory to model the aggregation and disaggregation of erythrocytes in normal human blood at different shear rates. To better understand the effect of rheological models on the haemodynamics of blood flow in cerebral aneurysms we compare our numerical results with those obtained with other rheological models such as the Carreau-Yasuda (C-Y) model. The results show that the velocity profiles for the Newtonian and the Owens models are approximately similar but differ from those of the C-Y model. In order to stabilize our numerical simulations, we propose two new stabilization techniques, the so-called N-Owens and I-Owens methods. Employing the N-Owens stabilization method enables us to capture the effect of erythrocyte aggregation in blood flow through a cerebral aneurysm at higher Weissenberg (We) and Reynolds (Re) numbers than would otherwise be possible.
Sand traps are used to measure Aeolian flux. Since they modify the surrounding wind velocity field their gauging represents an important challenge. We use numerical simulations under the assumption of homogeneous turbulence based on FLUENT to systematically study the flow field and trapping efficiency of one of the most common devices based on a hollow cylinder with two slits. In particular, we investigate the dependence on the wind speed, the Stokes number, the permeability of the membrane on the slit and the saltation height.