A hybrid, design-order sliding mesh algorithm, which uses a control volume finite element method (CVFEM), in conjunction with a discontinuous Galerkin (DG) approach at non-conformal interfaces, is ...outlined in the context of a low-Mach fluid dynamics equation set. This novel hybrid DG approach is also demonstrated to be compatible with a classic edge-based vertex centered (EBVC) scheme. For the CVFEM, element polynomial, P, promotion is used to extend the low-order P=1 CVFEM method to higher-order, i.e., P=2. An equal-order low-Mach pressure-stabilized methodology, with emphasis on the non-conformal interface boundary condition, is presented. A fully implicit matrix solver approach that accounts for the full stencil connectivity across the non-conformal interface is employed. A complete suite of formal verification studies using the method of manufactured solutions (MMS) is performed to verify the order of accuracy of the underlying methodology. The chosen suite of analytical verification cases range from a simple steady diffusion system to a traveling viscous vortex across mixed-order non-conformal interfaces. Results from all verification studies demonstrate either second- or third-order spatial accuracy and, for transient solutions, second-order temporal accuracy. Significant accuracy gains in manufactured solution error norms are noted even with modest promotion of the underlying polynomial order. The paper also demonstrates the CVFEM/DG methodology on two production-like simulation cases that include an inner block subjected to solid rotation, i.e., each of the simulations include a sliding mesh, non-conformal interface. The first production case presented is a turbulent flow past a high-rate-of-rotation cube (Re, 4000; RPM, 3600) on like and mixed-order polynomial interfaces. The final simulation case is a full-scale Vestas V27 225 kW wind turbine (tower and nacelle omitted) in which a hybrid topology, low-order mesh is used. Both production simulations provide confidence in the underlying capability and demonstrate the viability of this hybrid method for deployment towards high-fidelity wind energy validation and analysis.
We examine the application of neural network-based methods to improve the accuracy of large eddy simulations of incompressible turbulent flows. The networks are trained to learn a mapping between ...flow features and the subgrid scales, and applied locally and instantaneously—in the same way as traditional physics-based subgrid closures. Models that use only the local resolved strain rate are poorly correlated with the actual subgrid forces obtained from filtering direct numerical simulation data. We see that highly accurate models in a priori testing are inaccurate in forward calculations, owing to the preponderance of numerical errors in implicitly filtered large eddy simulations. A network that accounts for the discretization errors is trained and found to be unstable in a posteriori testing. We identify a number of challenges that the approach faces, including a distribution shift that affects networks that fail to account for numerical errors.
•Framework presented to estimate model-form uncertainty in LES closures - Independent of initial model form - Uncertainty in terms of magnitude, shape and orientation (eigenspace-based) - Physically ...plausible bounds derived for each degree of freedom - Computationally efficient and suitable to general solvers•Computational dataset of turbulent jet validated against experimental data - First- and second-order statistics present complete agreement with experiment - Base LES of turbulent axisymmetric jet using WALE (eddy-viscosity) model•Important differences between reference and modeled SGS stresses identified - Magnitude of SGS tends to be underpredicted - Extremely low anisotropy correlation between reference and modeled tensors•Eigensensitivity analysis performed to quantify main modeling uncertainties - Perturbations in tensor anisotropy (especially) produce larger relative impacts - Uncertainty estimates are able to envelope the reference solution•Future work focused on more complex flows and extending the framework
The study of complex turbulent flows by means of large-eddy simulation approaches has become increasingly popular in many scientific and engineering applications. The underlying filtering operation of the approach enables to significantly reduce the spatial and temporal resolution requirements by means of representing only large-scale motions. However, the small-scale stresses and their effects on the resolved flow field are not negligible, and therefore require additional modeling. As a consequence, the assumptions made in the closure formulations become potential sources of model-form uncertainty that can impact the quantities of interest. The objective of this work, thus, is to perform a model-form sensitivity analysis in large-eddy simulations of an axisymmetric turbulent jet following an eigenspace-based strategy recently proposed. The approach relies on introducing perturbations to the decomposed subgrid-scale stress tensor within a range of physically plausible values. These correspond to discrepancy in magnitude (trace), anisotropy (eigenvalues) and orientation (eigenvectors) of the normalized, small-scale stresses with respect to a given tensor state, such that propagation of their effects can be assessed. The generality of the framework with respect to the six degrees of freedom of the small-scale stress tensor makes it also suitable for its application within data-driven techniques for improved subgrid-scale modeling.
Motivated by the sizable increase of available computing resources, large-eddy simulation of complex turbulent flow is becoming increasingly popular. The underlying filtering operation of this ...approach enables to represent only large-scale motions. However, the small-scale fluctuations and their effects on the resolved flow field require additional modeling. As a consequence, the assumptions made in the closure formulations become potential sources of incertitude that can impact the quantities of interest. The objective of this work is to introduce a framework for the systematic estimation of structural uncertainty in large-eddy simulation closures. In particular, the methodology proposed is independent of the initial model form, computationally efficient, and suitable to general flow solvers. The approach is based on introducing controlled perturbations to the turbulent stress tensor in terms of magnitude, shape and orientation, such that propagation of their effects can be assessed. The framework is rigorously described, and physically plausible bounds for the perturbations are proposed. As a means to test its performance, a comprehensive set of numerical experiments are reported for which physical interpretation of the deviations in the quantities of interest are discussed.
•Design-order, variable density low-Mach methods have been established on atypical unstructured mesh topologies.•Credible, high-quality low-Mach large-eddy simulation results are achieved using ...generalized unstructured topologies.•Modest polynomial promotion provides substantial reduction in numerical error.•Large-eddy simulation studies can benefit from the usage of hybrid and tetrahedral meshing approaches.
An implicit, low-dissipation, low-Mach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for large-eddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between low-order hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a third-order, unstructured hexahedral topology. Using smooth analytical and manufactured low-Mach solutions, design-order convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energy-stable methodologies previously deployed within a finite-difference context. A wide range of simulations demonstrate that low-order hexahedral- and wedge-based element topologies behave nearly identically in both computed numerical errors and overall simulation timings. Moreover, low-order tetrahedral and pyramid element topologies also display nearly the same numerical characteristics. Although the superiority of the hexahedral-based topology is clearly demonstrated for trivial laminar, principally-aligned flows, e.g., a 1x2x10 channel flow with specified pressure drop, this advantage is reduced for non-aligned, turbulent flows including the Taylor–Green Vortex, turbulent plane channel flow (Reτ395), and buoyant flow past a heated cylinder. With the order of accuracy demonstrated for both homogenous and hybrid meshes, it is shown that solution verification for the selected complex flows can be established for all topology types. Although the number of elements in a mesh of like spacing comprised of tetrahedral, wedge, or pyramid elements increases as compared to the hexahedral counterpart, for wall-resolved large-eddy simulation, the increased assembly and residual evaluation computational time for non-hexahedral is offset by more efficient linear solver times. Finally, most simulation results indicate that modest polynomial promotion provides a significant increase in solution accuracy.
A high-fidelity large-eddy simulation and unsteady flamelet combustion model construct is deployed to numerically investigate the effects of crosswind magnitude and pool fire shape on large-scale ...pool fire attributes. These include general flame dynamics, flame shape and radiative flux magnitude in and around the fire. Three pool fire shapes at a nominal length scale of 10 m are subjected to four crosswind magnitudes between 0 and 20 m s
$^{-1}$
. The pool shapes studied are circular, square and rectangular. The study includes the sensitivity of parameters to mesh and time step refinement. Results demonstrate that the rectangular shape, under crosswind, has low-levels of vertical velocity induction, resulting in a plume that is closer to the ground. In the quiescent regime, under-resolved meshes provide a higher radiative heat flux prediction compared with the most refined mesh. However, as crosswind increases, low mesh resolutions underpredicted radiative flux. This is due to the coarse mesh resolution not capturing small-scale vortical features that increased mixing and combustion efficiency. A transition of peak radiative flux with respect to crosswind occurs from the leeward- to windward-side of the pool, while sharp pool features result in larger radiative heat fluxes concentrated in regions of high scalar dissipation rate.
The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used ...to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.
A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission ...probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen transport from a synthetic cough emanating from a kneeling humanoid. For an outdoor configuration that mimics the recent open-space social distance strategy of San Francisco, maximum primary droplet deposition distances are shown to approach 8.1 m in a moderate wind configuration with the aerosol plume transported in excess of 15 m. In quiescent conditions, the aerosol plume extends to approximately 4 m before the emanating pulsed jet becomes neutrally buoyant. A dose-response model, which is based on previous SARS coronavirus (SARS-CoV) data, is exercised on the high-fidelity aerosol transport database to establish relative risk at eighteen virtual receptor probe locations.