An official website of the United States government
Reduced-order models (ROMs) have achieved a lot of success in reducing the computational cost of traditional numerical methods across many disciplines. In fluid dynamics, ROMs have been successful in providing efficient and relatively accurate solutions for the numerical simulation of laminar flows. For convection-dominated (e.g., turbulent) flows, however, standard ROMs generally yield inaccurate results, usually affected by spurious oscillations. Thus, ROMs are usually equipped with numerical stabilization or closure models in order to account for the effect of the discarded modes. The literature on ROM closures and stabilizations is large and growing fast. In this paper, instead of reviewing all the ROM closures and stabilizations, we took a more modest step and focused on one particular type of ROM closure and stabilization that is inspired by large eddy simulation (LES), a classical strategy in computational fluid dynamics (CFD). These ROMs, which we call LES-ROMs, are extremely easy to implement, very efficient, and accurate. Indeed, LES-ROMs are modular and generally require minimal modifications to standard (“legacy”) ROM formulations. Furthermore, the computational overhead of these modifications is minimal. Finally, carefully tuned LES-ROMs can accurately capture the average physical quantities of interest in challenging convection-dominated flows in science and engineering applications. LES-ROMs are constructed by leveraging spatial filtering, which is the same principle used to build classical LES models. This ensures a modeling consistency between LES-ROMs and the approaches that generated the data used to train them. It also “bridges” two distinct research fields (LES and ROMs) that have been disconnected until now. This paper is a review of LES-ROMs, with a particular focus on the LES concepts and models that enable the construction of LES-inspired ROMs and the bridging of LES and reduced-order modeling. This paper starts with a description of a versatile LES strategy called evolve–filter–relax (EFR) that has been successfully used as a full-order method for both incompressible and compressible convection-dominated flows. We present evidence of this success. We then show how the EFR strategy, and spatial filtering in general, can be leveraged to construct LES-ROMs (e.g., EFR-ROM). Several applications of LES-ROMs to the numerical simulation of incompressible and compressible convection-dominated flows are presented. Finally, we draw conclusions and outline several research directions and open questions in LES-ROM development. While we do not claim this review to be comprehensive, we certainly hope it serves as a brief and friendly introduction to this exciting research area, which we believe has a lot of potential in the practical numerical simulation of convection-dominated flows in science, engineering, and medicine.
more »
« less
Data‐Driven Optimization for the Evolve‐Filter‐Relax Regularization of Convection‐Dominated Flows
ABSTRACT Numerical stabilization techniques are often employed in under‐resolved simulations of convection‐dominated flows to improve accuracy and mitigate spurious oscillations. Specifically, the evolve–filter–relax (EFR) algorithm is a framework that consists of evolving the solution, applying a filtering step to remove high‐frequency noise, and relaxing through a convex combination of filtered and original solutions. The stability and accuracy of the EFR solution strongly depend on two parameters, the filter radius and the relaxation parameter . Standard choices for these parameters are usually fixed in time, and related to the full order model setting, that is, the grid size for and the time step for . The key novelties with respect to the standard EFR approach are: (i) time‐dependent parameters and , and (ii) data‐driven adaptive optimization of the parameters in time, considering a fully‐resolved simulation as reference. In particular, we propose three different classes of optimized‐EFR (Opt‐EFR) strategies, aiming to optimize one or both parameters. The new Opt‐EFR strategies are tested in the under‐resolved simulation of a turbulent flow past a cylinder at . The Opt‐EFR proved to be more accurate than standard approaches by up to 99, while maintaining a similar computational time. In particular, the key new finding of our analysis is that such accuracy can be obtained only if the optimized objective function includes: (i) aglobalmetric (as the kinetic energy), and (ii)spatial gradients' information.
more »
« less
- Award ID(s):
- 2012253
- PAR ID:
- 10590166
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Engineering
- Volume:
- 126
- Issue:
- 9
- ISSN:
- 0029-5981
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Offline particle tracking (OPT) is a widely used tool for theanalysis of data in oceanographic research. Given the output of ahydrodynamic model, OPT can provide answers to a wide variety of researchquestions involving fluid kinematics, zooplankton transport, the dispersionof pollutants, and the fate of chemical tracers, among others. In thispaper, we introduce ROMSPath, an OPT model designed to complement theRegional Ocean Modeling System (ROMS). Based on the Lagrangian TRANSport(LTRANS) model (North et al., 2008), ROMSPath is written in Fortran90 and provides advancements in functionality and efficiency compared toLTRANS. First, ROMSPath calculates particle trajectories using the ROMSnative grid, which provides advantages in interpolation, masking, andboundary interaction while improving accuracy. Second, ROMSPath enablessimulated particles to pass between nested ROMS grids, which is anincreasingly popular scheme to simulate the ocean over multiple scales.Third, the ROMSPath vertical turbulence module enables the turbulent(diffusion) time step and advection time step to be specified separately,adding flexibility and improving computational efficiency. Lastly, ROMSPathincludes new infrastructure which enables inputting of auxiliary parameters for addedfunctionality. In particular, Stokes drift can be input and added toparticle advection. Here we describe the details of these updates andperformance improvements.more » « less
-
The crystal structures of two intermediates, 4-amino-3,5-difluorobenzonitrile, C7H4F2N2(I), and ethyl 4-amino-3,5-difluorobenzoate, C9H9F2NO2(II), along with a visible-light-responsive azobenzene derivative, diethyl 4,4′-(diazene-1,2-diyl)bis(3,5-difluorobenzoate), C18H14F4N2O4(III), obtained by four-step synthetic procedure, were studied using single-crystal X-ray diffraction. The molecules ofIandIIdemonstrate the quinoid character of phenyl rings accompanied by the distortion of bond angles related to the presence of fluorine substituents in the 3 and 5 (ortho) positions. In the crystals ofIandII, the molecules are connected by N—H...N, N—H...F and N—H...O hydrogen bonds, C—H...F short contacts, and π-stacking interactions. In crystal ofIII, only stacking interactions between the molecules are found.more » « less
-
In Part I [Appl. Opt.59,1018(2020).APOPAI0003-693510.1364/AO.381246], we used a coupled optoelectronic model to optimize a thin-film AlGaAs solar cell with a graded-bandgap photon-absorbing layer and a periodically corrugated Ag backreflector combined with localized ohmic Pd–Ge–Au backcontacts, because both strategies help to improve the performance of AlGaAs solar cells. However, the results in Part I were affected by a normalization error, which came to light when we replaced the hybridizable discontinuous Galerkin scheme for electrical computation by the faster finite-difference scheme. Therefore, we re-optimized the solar cells containing ann-AlGaAs photon-absorbing layer with either a (i) homogeneous, (ii) linearly graded, or (iii) nonlinearly graded bandgap. Another way to improve the power conversion efficiency is by using a surface antireflection texturing on the wavelength scale, so we also optimized four different types of 1D periodic surface texturing: (i) rectangular, (ii) convex hemi-elliptical, (iii) triangular, and (iv) concave hemi-elliptical. Our new results show that the optimal nonlinear bandgap grading enhances the efficiency by as much as 3.31% when then-AlGaAs layer is 400 nm thick and 1.14% when that layer is 2000 nm thick. A hundredfold concentration of sunlight can enhance the efficiency by a factor of 11.6%. Periodic texturing of the front surface on the scale of 0.5–2 free-space wavelengths provides a small relative enhancement in efficiency over the AlGaAs solar cells with a planar front surface; however, the enhancement is lower when then-AlGaAs layer is thicker.more » « less
-
Bis(triphenylsulfonium) tetrachloridomanganate(II), (C18H15S)2[MnCl4] (I), triphenylsulfonium tetrachloridoferrate(III), (C18H15S)[FeCl4] (II), and bis(triphenylsulfonium) tetrachloridocobaltate(II), (C18H15S)2[CoCl4] (III), crystallize in the monoclinic space groupsP21/n[(I) and (III)] andP21/c[(II)]. Compounds (I) and (III) each contain two crystallographically independent triphenylsulfonium (TPS+) cations in the asymmetric unit, whereas (II) has one. In all three compounds, the sulfonium centers adopt distorted trigonal–pyramidal geometries, with S—C bond lengths falling roughly in the 1.78–1.79 Å range and C—S—C angles observed at about 101 to 106°. The [MCl4]n−anions (M= Mn2+, Fe3+, Co2+;n= 2,1,2) adopt slightly distorted tetrahedral geometries, withM—Cl bond lengths in the 2.19–2.38 Å range and Cl—M—Cl angles of approximately 104–113°. Hirshfeld surface analyses shows that H...H and H...C contacts dominate the TPS+cation environments, whereas H...Cl and shortM—S interactions link each [MCl4]n−anion to the surrounding cations. In (I) and (III), inversion-centered π–π stacking further consolidates the crystal packing, while in (II) no π–π interactions are observed.more » « less
