Scalar Forcing Methodology for Direct Numerical Simulations of Turbulent Stratified Mixture Combustion

Scalar forcing in the context of turbulent stratified flame simulations aims to maintain the fuel-air inhomogeneity in the unburned gas. With scalar forcing, stratified flame simulations have the potential to reach a statistically stationary state with a prescribed mixture fraction distribution and root-mean-square value in the unburned gas, irrespective of the turbulence intensity. The applicability of scalar forcing for Direct Numerical Simulations of stratified mixture combustion is assessed by considering a recently developed scalar forcing scheme, known as the reaction analogy method, applied to both passive scalar mixing and the imperfectly mixed unburned reactants of statistically planar stratified flames under low Mach number conditions. The newly developed method enables statistically symmetric scalar distributions between bell-shaped and bimodal to be maintained without any significant departure from the specified bounds of the scalar. Moreover, the performance of the newly proposed scalar forcing methodology has been assessed for a range of different velocity forcing schemes (Lundgren forcing and modified bandwidth forcing) and also without any velocity forcing. It has been found that the scalar forcing scheme has no adverse impact on flame-turbulence interaction and it only maintains the prescribed root-mean-square value of the scalar fluctuation, and its distribution. The scalar integral length scale evolution is shown to be unaffected by the scalar forcing scheme studied in this paper. Thus, the scalar forcing scheme has a high potential to provide a valuable computational tool to enable analysis of the effects of unburned mixture stratification on turbulent flame dynamics.

CFD Analysis of a Large Marine Engine Scavenging Process

The scavenging process is an important part of the two-stroke engine operation. Its efficiency affects the global engine performance such as power, fuel consumption, and pollutant emissions. Slow speed marine diesel engines are uniflow scavenged, which implies inlet scavenging ports on the bottom of the liner and an exhaust valve on the top of the cylinder. A CFD model of such an engine process was developed with the OpenFOAM software tools. A 12-degree sector of the mesh was used corresponding to one of the 30 scavenging ports. A mesh sensitivity test was performed, and the cylinder pressure was compared to experimental data for the analyzed part of the process. The scavenging performances were analyzed for real operation parameters. The influence of the scavenge air pressure and inlet ports geometric orientation was analyzed. The scavenging process is analyzed by means of a passive scalar representing fresh air in the cylinder. Isosurfaces that show the concentration of fresh air were presented. The variation of oxygen and carbon dioxide with time and the axial and angular momentum in the cylinder were calculated. Finally, the scavenging performance for the various operation parameters was evaluated by means of scavenging efficiency, charging efficiency, trapping efficiency, and delivery ratio. It was found that the scavenging efficiency decreases with the engine load due to the shorter time for the process. The scavenging efficiency increases with the pressure difference between the exhaust and scavenging port, and the scavenging efficiency decreases with the increase in the angle of the scavenging ports. It was concluded that smaller angles than the industry standard of 20° could be beneficial to the scavenging efficiency. In the investigation, the charging efficiency ranged from 0.91 to over 0.99, the trapping efficiency ranged from 0.54 to 0.83, the charging efficiency ranged from 0.78 to 0.92, and the delivery ratio ranged from 1.21 to 2.03.

Numerical Flow Modeling of Artificial Ocean Upwelling

Artificial Upwelling (AU) of nutrient-rich Deep Ocean Water (DOW) to the ocean's sunlit surface layer has recently been put forward as a means of increasing marine CO2 sequestration and fish production. AU and its possible benefits have been studied in the context of climate change mitigation as well as food security for a growing human population. However, extensive research still needs to be done into the feasibility, effectiveness and potential risks, and side effects associated with AU to be able to better predict its potential. Fluid dynamic modeling of the AU process and the corresponding inorganic nutrient transport can provide necessary information for a better quantification of the environmental impacts of specific AU devices and represents a valuable tool for their optimization. Yet, appropriate capture of all flow phenomena relevant to the AU process remains a challenging task that only few models are able to accomplish. In this paper, simulation results obtained with a newly developed numerical solution method are presented. The method is based on the open-source modeling environment OpenFOAM. It solves the unsteady Reynolds-Averaged Navier-Stokes (RANS) equations with additional transport equations for energy, salinity, and inorganic nutrients. The method aims to be widely applicable to oceanic flow problems including temperature- and salinity-induced density stratification and passive scalar transport. The studies presented in this paper concentrate on the direct effects of the AU process on nutrient spread and concentration in the ocean's mixed surface layer. Expected flow phenomena are found to be captured well by the new method. While it is a known problem that cold DOW that is upwelled to the surface tends to sink down again due to its high density, the simulations presented in this paper show that the upwelled DOW settles at the lower boundary of the oceans mixed surface layer, thus keeping a considerable portion of the upwelled nutrients available for primary production. Comparative studies of several design variants, with the aim of maximizing the amount of nutrients that is retained inside the mixed surface layer, are also presented and analyzed.

Simulation of the Scalar Transport above and within the Amazon Forest Canopy

The parallelized large-eddy simulation model (PALM) was used to understand better the turbulent exchanges of a passive scalar above and within a forested region located in the central Amazon. Weak (2 ms−1) and strong (6 ms−1) wind conditions were simulated. A passive scalar source was introduced to the forest floor for both simulations. The simulations reproduced the main characteristics of the turbulent flow and of the passive scalar transport between the forest and the atmosphere. Noteworthily, strong and weak wind conditions presented different turbulence structures that drove different patterns of scalar exchange both within and above the forest. These results show how passive scalar concentration is influenced by the wind speed at the canopy top. Additionally, higher wind speeds are related to stronger sweep and ejection regimes, generating more intense plumes that are able to reduce the passive scalar concentration inside the forest canopy. This work was the first that used PALM to investigate scalar transport between the Amazon rainforest and the atmosphere.

Investigation of internal intermittency by way of higher-order spectral moments

The analysis of turbulence by way of higher-order spectral moments is uncommon, despite the relatively frequent use of such statistical analyses in other fields of physics and engineering. In this work, higher-order spectral moments are used to investigate the internal intermittency of the turbulent velocity and passive-scalar (temperature) fields. This study first introduces the theory behind higher-order spectral moments as they pertain to the field of turbulence. Then, a short-time Fourier-transform-based method is developed to estimate these higher-order spectral moments and provide a relative, scale-by-scale measure of intermittency. Experimental data are subsequently analysed and consist of measurements of homogeneous, isotropic, high-Reynolds-number, passive and active grid turbulence over the Reynolds-number range $35\leq R_{\lambda } \leq ~731$ . Emphasis is placed on third- and fourth-order spectral moments using the definitions formalised by Antoni (Mech. Syst. Signal Pr., vol. 20 (2), 2006, pp. 282–307), as such statistics are sensitive to transients and provide insight into deviations from Gaussian behaviour in grid turbulence. The higher-order spectral moments are also used to investigate the Reynolds (Péclet) number dependence of the internal intermittency of velocity and passive-scalar fields. The results demonstrate that the evolution of higher-order spectral moments with Reynolds number is strongly dependent on wavenumber. Finally, the relative levels of internal intermittency of the velocity and passive-scalar fields are compared and a higher level of internal intermittency in the inertial subrange of the scalar field is consistently observed, whereas a similar level of internal intermittency is observed for the velocity and passive-scalar fields for the high-Reynolds-number cases as the Kolmogorov length scale is approached.

On mixing enhancement by secondary baroclinic vorticity in a shock–bubble interaction

To investigate the intrinsic mechanism for mixing enhancement by variable-density (VD) behaviour, a canonical VD mixing extracted from a supersonic streamwise vortex protocol, a shock–bubble interaction (SBI), is numerically studied and compared with a counterpart of passive-scalar (PS) mixing. It is meaningful to observe that the maximum concentration decays much faster in a VD SBI than in a PS SBI regardless of the shock Mach number ( $Ma=1.22 - 4$ ). The quasi-Lamb–Oseen-type velocity distribution in the PS SBI is found by analysing the azimuthal velocity that stretches the bubble. Meanwhile, for the VD SBI, an additional stretching enhanced by the secondary baroclinic vorticity (SBV) production contributes to the faster-mixing decay. The underlying mechanism of the SBV-enhanced stretching is further revealed through the density and velocity difference between the light shocked bubble and the heavy ambient air. By combining the SBV-accelerated stretching model and the initial shock compression, a novel mixing time estimation for VD SBI is theoretically proposed by solving the advection–diffusion equation under a deformation field of an axisymmetric vortex with the additional SBV-induced azimuthal velocity. Based on the mixing time model, a mixing enhancement number, defined by the ratio of VD and PS mixing time further, reveals the contribution from the VD effect, which implies a better control of the density distribution for mixing enhancement in a supersonic streamwise vortex.