A nonlinear Harmonic Navier-Stokes (HNS) framework is introduced for simulating instabilities in laminar spanwise-invariant shear layers, featuring sharp and smooth wall surface protuberances. While such cases play a critical role in the process of laminar-to-turbulent transition, classical stability theory analyses such as parabolized or local stability methods fail to provide (accurate) results, due to their underlying assumptions. The generalized incompressible Navier-Stokes (NS) equations are expanded in perturbed form, using a spanwise and temporal Fourier ansatz for flow perturbations. The resulting equations are discretized using spectral collocation in the wall-normal direction and finite-difference methods in the streamwise direction. The equations are then solved using a direct sparse-matrix solver. The nonlinear mode interaction terms are converged iteratively. The solution implementation makes use of a generalized domain transformation to account for geometrical smooth surface features, such as humps. No-slip conditions can be embedded in the interior domain to account for the presence of sharp surface features such as forward- or backward-facing steps. Common difficulties with Navier-Stokes solvers, such as the treatment of the outflow boundary and convergence of nonlinear terms, are considered in detail. The performance of the developed solver is evaluated against several cases of representative boundary layer instability growth, including linear and nonlinear growth of Tollmien-Schlichting waves in a Blasius boundary layer and stationary crossflow instabilities in a swept flat-plate boundary layer. The latter problem is also treated in the presence of a geometrical smooth hump and a sharp forward-facing step at the wall. HNS simulation results, such as perturbation amplitudes, growth rates, and shape functions, are compared to benchmark flow stability analysis methods such as Parabolized Stability Equations (PSE), Adaptive Harmonic Linearized Navier-Stokes (AHLNS), or Direct Numerical Simulations (DNS). Good agreement is observed in all cases. The HNS solver is subjected to a grid convergence study and a simple performance benchmark, namely memory usage and computational cost. The computational cost is found to be considerably lower than high-fidelity DNS at comparable grid resolutions.@en

Indirect combustion noise due to the interaction of flow inhomogeneities with a choked combustion-chamber exit is an important cause of combustion instability in solid rocket motors. Moreover, it is believed to be an issue in electrical-power generation turbines and aero-engines. If these flow inhomogeneities are essentially characterized by the fluid having a locally appreciably-different thermodynamic state, the acoustic response engendered by its interaction with the combustion-chamber exit is commonly referred to as entropy noise. In this paper, dedicated numerical-simulation results of entropy-patch choked-nozzle interactions are presented. Two types of entropy patches were considered: rectangular slugs and circular spots. Moreover, analytical-model-based analysis, of said simulation results, is presented. Based on said analysis, the authors posit the existence of three modeling regimes: the quasi-steady-modeling regime, the blended-physical-effects regime, and the inertial-modeling regime. @en

In §II. Theory, two reduced-order models are proposed, which the authors have termed: the quasi-steady model (§II.A. Quasi-steady one-dimensional model) & the inertial/hybrid model (§II.B. Quasi-one-dimensional pointmass model), respectively. N.b., in both cases time dependence isn’t explicitly modeled, i.e., technically speaking both models are quasi-steady. Ergo, in hindsight, it would have been more apt to call the model proposed in §II.A.: the matching-condition model. With that in mind, the readership is encouraged to substitute “matching-condition model/modeling regime” instead of “quasi-steady model/modeling regime,” when reading this conference paper. Moreover, the following title would have been more suitable: “Entropy-patch chokednozzle interaction: matching-condition and inertial modeling-regimes mapped”.

@en

In this paper, we build on the work of Hughes and Sangalli (2007) dealing with the explicit computation of the Fine-Scale Greens’ function. The original approach chooses a set of functionals associated with a projector to compute the Fine-Scale Greens’ function. The construction of these functionals, however, does not generalise to arbitrary projections, higher dimensions, or Spectral Element methods. We propose to generalise the construction of the required functionals by using dual functions. These dual functions can be directly derived from the chosen projector and are explicitly computable. We show how to find the dual functions for both the L² and the H₀¹ projections. We then go on to demonstrate that the Fine-Scale Greens’ functions constructed with the dual basis functions consistently reproduce the unresolved scales removed by the projector. The methodology is tested using one-dimensional Poisson and advection–diffusion problems, as well as a two-dimensional Poisson problem. We present the computed components of the Fine-Scale Greens’ function, and the Fine-Scale Greens’ function itself. These results show that the method works for arbitrary projections, in arbitrary dimensions. Moreover, the methodology can be applied to any Finite/Spectral Element or Isogeometric framework.

@en

Numerical simulations and optimisation methods, such as mesh adaptation, rely on the accurate and inexpensive use of error estimation methods. Adjoint-based error estimation is the most accurate method, and generally the most costly. A strong contributor to this cost is the need to compute a higher resolution adjoint solution. Here, it is proposed to use super-resolution neural networks to super-resolve a fine adjoint solution from a lower-cost coarse adjoint solution: a superAdjoint. The method is compared to reference error estimators on an unsteady Burgers’ equation using the method of manufactured solutions. Two forms of the superAdjoint were implemented, a twice and a four times refining super-resolution neural network. These were used to demonstrate both the computational cost reduction and the potential for the reduction of the storage footprint of the primal problem. The first, referred to as 2×CNN, was able to reconstruct the spatially enriched adjoint solution, thus providing a robust and inexpensive local output error. The second, the 4×CNN, was able to demonstrate the reconstruction ability of super-resolution neural networks for higher upscaling factors. This was leveraged in order to subsample the primal solution in space, thus reducing substantially the storage footprint of the discrete primal solution. Both superAdjoints could achieve the desired level of accuracy when compared to the reconstruction of the refined adjoint-based error estimate. Moreover, superAdjoint was shown to be able to generalise to a new QoI that was untrained for. This gives great confidence in the use of super-resolution neural networks for the reduction of both computational cost and storage requirements of adjoint-based error estimation, and goal-oriented mesh adaptation.

@en

An experimental cold-gas study of the response of a choked convergent–divergent nozzle to swirl perturbations is presented. The perturbations were obtained by means of upstream unsteady tangential injections into initially steady flows with different values of steady background swirl. The swirl perturbations induced changes in the axial mass-flow rate, due to either their ingestion or evacuation by the nozzle. This in turn caused a downstream acoustic response. For low-intensity background swirl the responses were found to be similar to those obtained without steady background swirl. Perturbations of a high-intensity background swirl led to different effects. For long injection times, the negative mass-flow rate modulation occurred in two stages. The first stage was similar to that of the background-swirl free case. The second stage occurred after a short time delay, and induced a much stronger negative acoustic response. This unexpected behavior suggests that a significant part of the tangentially injected fluid flows upstream inducing an accumulation of swirl, which is – after tangential injection is ceased – suddenly cleared out through the nozzle. A scaling rule for the amplitudes of these acoustic responses is reported. Furthermore, quasi-steady models, based on steady-state measurements are proposed. These models predict the downstream acoustic response amplitude within a factor two. Additionally, preliminary empirical evidence of the effect of swirl on the downstream acoustic response due to the interaction of entropy patches with a choked nozzle is reported. This was obtained by comparison of sound produced by abrupt radial or tangential sonic injection, upstream from the choked nozzle, of air from a reservoir at room temperature to that from a reservoir with a higher stagnation temperature. Because the mass flow through the nozzle does not increase instantaneously, the injected higher-enthalpy air accumulates upstream of the injection-port position in the main flow. This eventually induces a large downstream acoustic pulse when tangential injection is interrupted. The magnitude of the resulting sound pulse can reach that of a quasi-steady response of the nozzle to a large air patch with a uniform stagnation temperature equal to that of the upstream-injected heated air. This hypothesis is consistent with the fact that the initial indirect-sound pulse is identical to one obtained with unheated air injection. The authors posit that – given all of the insight gleaned from them in this case – acoustic measurements of indirect sound appear to be a potentially useful diagnostic tool.

@en

This work presents a transpacific airliner designed for minimal climate impact, incorporating several novel design features. These include open rotor engines, sustainable aviation fuels, natural laminar flow airfoils, and riblets. The design’s configuration and mission have been optimised simultaneously using a combination of standard preliminary techniques, experimental data, a multi-point mission analysis, and a model of average temperature response. It is demonstrated that, on an 8000 km mission, the design offers an 89.8% reduction in average temperature response relative to an Airbus A330-200, at the expense of a 7.3% increase in direct operating cost. The sensitivity of these results is investigated by comparing the performance over a range of operating conditions. In addition, several alternative designs incorporating only some of the above-mentioned features are analysed, allowing for an assessment of their individual contribution. Finally, a life-cycle average temperature response analysis is presented to place the climate impact of operation, manufacturing and end-of-life procedures in context@en

The underlying physical mechanism of the residual-based large eddy simulation (LES) based on the variational multiscale (VMS) method is clarified. Resolved large-scale energy transportation equation is mathematically derived for turbulent kinetic energy budget analysis. Firstly, statistical results of benchmark turbulent channel flow at Re_τ=180 obtained using a coarse mesh are compared with the results obtained by the classical LES with the Smagorinsky and dynamic subgrid stress (SGS) model. The present LES shows an advantage in predicting the statistical results of the incompressible turbulent flows. Secondly, the contributions of the unresolved small-scale presentation terms (Term I-IV in Eq. (10)) to the turbulent kinetic dissipation are analysed for the VMS method. The results show that the turbulent kinetic dissipation provided by the numerical diffusion in the VMS method is smaller in the inner layer, larger in the outer layer of the channel flow than those by the Smagorinsky and dynamic SGS model. The turbulent kinetic dissipation in the VMS method is mainly given by the numerical diffusion provided by one of the “cross-stress” terms (Term I, same as the stabilization term in the SUPG method) and LSIC term (Term IV). The other one of the “cross-stress” terms (Term II) gives rise to the positive turbulent kinetic energy budget, and does not dissipate the turbulent kinetic energy. The so-called “Reynolds stress” term (Term III) dissipates the turbulent energy but provides a very small numerical diffusion. Finally, on the basis of the turbulent kinetic energy dissipation analysis, a new residual-based stabilized finite element formulation is proposed by modifying the large-scale equation in the VMS method. Numerical experiments of 2D lid-driven cavity flow and 3D incompressible turbulent channel flow are tested to validate the proposed formulation. It is shown that all the stabilization terms in the proposed formulation produce additional numerical diffusions and physically increase the total turbulent kinetic dissipation. Consequently, an apparent improvement in both the first-order and second-order statistical quantities are pursued by the new stabilized finite element formulation.

@en

Data-driven parameterizations offer considerable potential for improving the fidelity of General Circulation Models. However, ensuring that these remain consistent with the governing equations while still producing stable simulations remains a challenge. In this paper, we propose a combined Variational-Multiscale (VMS) Artificial Neural Network (ANN) discretization which makes no a priori assumptions on the model form, and is only restricted in its accuracy by the precision of the ANN. Using a simplified problem, we demonstrate that good predictions of the required closure terms can be obtained with relatively compact ANN architectures. We then turn our attention to the stability of the VMS-ANN discretization in the context of a single implicit time step. It is demonstrated that the ANN parameterization introduces nonphysical solutions to the governing equations that can significantly affect or prevent convergence. We show that enriching the training data with nonphysical states from intra-time step iterations is an effective remedy. This indicates that the lack of representative ANN-induced errors in our original, exact training inputs underpin the observed instabilities. In turn, this suggests that data set enrichment might aid in resolving instabilities that develop over several time steps.

@en

The downstream acoustic response due to the interaction of main-flow oriented vorticity with a choked nozzle in a swirl-free flow was experimentally demonstrated. The response was obtained by means of impulsive radial air injection in the pipe upstream from the nozzle. The resulting downstream acoustic data are shown to obey a scaling rule that differs, from the one for swirl-nozzle interaction, which according to the literature is proportional to the square of the swirl number. In contrast, here evidence is presented that points to the scaling of main-flow oriented vorticity noise with the cross-sectional average of the square of the transversal velocity at the throat divided by the square of the critical sound speed.@en

Proper Orthogonal Decomposition (POD) plays an important role in the analysis of complex nonlinear systems governed by partial differential equations (PDEs), since it can describe the full-order system in a simplified but representative way using a handful of dominant dynamic modes. However, determining a POD from the results of complex unsteady simulations is often impractical using traditional approaches due to the need to store a large number of high-dimensional solutions. As an alternative, the incremental Singular Value Decomposition (SVD) has been developed, which can be used to avoid the storage problem by performing the POD analysis on the fly using a single-pass updating algorithm. Nevertheless, the total computing cost of incremental SVD is more than traditional approaches. In order to reduce this total cost, we incorporate POD mode truncation into the incremental procedure, leading to an enhanced algorithm for incremental SVD. Two error estimators are formulated for this enhanced incremental SVD based on an aggregated expression of the snapshot solutions, equipping the proposed algorithm with criteria for choosing the truncation number. The effectiveness of these estimators and the parallel efficiency of the enhanced algorithm are demonstrated using transient solutions from representative model problems. Numerical results show that the enhanced algorithm can significantly improve the computing efficiency for different kinds of datasets, and that the proposed algorithm is scalable in both the strong and weak sense.

@en

View Video Presentation: https://doi-org.tudelft.idm.oclc.org/10.2514/6.2021-2286.vid

Quantitative measurements of sound due to swirl-nozzle interaction are presented for the first time. In the experiment a swirl structure was generated by means of tangential injection into a steady swirl-free flow upstream from a choked convergent-divergent nozzle. Ingestion of swirl by the choked nozzle caused a mass-flow rate change, which resulted in a downstream measured acoustic response. The amplitude of this acoustic response was found to be proportional to the square of the tangential mass-flow rate used to generate swirl. This was, assuming that the upstream generated swirl intensity is proportional to the tangential injection mass-flow rate, predicted by a previously published quasi-steady model for the swirl-nozzle interaction sound source (Hirschberg, L., Hulshoff, S. J., and Bake, F., “Sound Production due to Swirl-Nozzle Interaction: Model-Based Analysis of Experiments,” AIAA Journal, Published online on Nov. 11th 2020, doi: 10.2514/1.J059669.). The tangential-injection time was varied, and found to not influence the amplitude of the acoustic response. This indicates that quasi-steady modelling remains applicable, even for smallest achievable upstream swirl structure with an axial length of ca. three upstream diameters.@en

Adaptive Mesh Refinement (AMR) is potentially an effective way to automatically generate computational meshes for high-fidelity simulations such as Large Eddy Simulation (LES). When combined with adjoint methods, which are able to localize error contributions, AMR can generate meshes that are optimal for computing a physical quantity of interest (e.g. lift or drag). In order to apply adjoint-based AMR techniques to LES, primal flow solutions are needed to solve the adjoint problem backward in time. However, the resources required to store primal flow solutions can be huge, even prohibitive, in practical problems because of the typically very fine meshes and long averaging times for computing the statistical quantities of interest. Here, a Reduced-Order Representation (ROR) based upon proper orthogonal decomposition is introduced to address this issue. We develop an Enhanced Online Algorithm (EOA) based on incremental singular value decomposition to build this ROR online, which makes adjoint-based AMR feasible for practical applications. An adjoint-based error estimation procedure is first introduced, and verified using a manufactured solution. Then a ROR-driven AMR strategy is studied using a 1D unsteady Burgers problem with a multi-frequency forcing term. This is also used to evaluate the EOA for ROR-driven AMR. Numerical results demonstrate that the enhanced online algorithm generates RORs that are sufficiently accurate for AMR, avoiding the storage of almost all of the primal solution data.

@en

Adaptive mesh refinement (AMR) is potentially an effective way to automatically generate computational meshes for high-fidelity simulations such as Large Eddy Simulation (LES). Adjoint methods, which are able to localize error contributions, can be used to optimize the mesh for computing a physical quantity of interest (e.g. lift, drag) during AMR. When adjoint-based AMR techniques are applied to LES, primal flow solutions are needed to solve the adjoint problem backward in time due to the nonlinearity of Navier-Stokes equations. However, the resources required to store primal flow solutions can be huge, even prohibitive, in practical problems because of the long averaging time for computing statistical quantities. In this paper, a Reduced-Order Model (ROM) based upon Proper Orthogonal Decomposition (POD) is introduced to circumvent this issue. First, an adjoint-based error estimation procedure is verified using a manufactured solution. Then a ROM-driven AMR strategy is studied using a LES model problem based on the 1D unsteady Burgers equation. Numerical results demonstrate that using ROMs not only lowers storage requirements, but also has no impact on the effectiveness of adjoint-based AMR.

@en

Measurements of sound due to swirl–nozzle interaction are presented. In the experiment a swirl structure was generated by means of unsteady tangential injection into a steady swirl-free flow upstream from a choked convergent–divergent nozzle. Ingestion of swirl by the choked nozzle caused a mass-flow rate change, which resulted in a downstream-measured acoustic response. The downstream acoustic pressure was found to remain negative as long as the swirl is maintained and reflections from the open downstream pipe termination do not interfere. The amplitude of this initial acoustic response was found to be proportional to the square of the tangential mass-flow rate used to generate swirl. When the tangential injection valve was closed, the mass-flow rate through the nozzle increased, resulting in an increase of the downstream acoustic pressure. This increase in signal was compared to the prediction of an empirical quasi-steady model, constructed from steady-state flow measurements. As the opening time of the valve was varied, the signal due to swirl evacuation showed an initial overshoot with respect to quasi-steady behavior, after which it gradually decayed to quasi-steady behavior for tangential injection times long compared to the convection time in the pipe upstream of the nozzle. This demonstrates that the acoustic signal can be used to obtain quantitative information concerning the time dependence of the swirl in the system. This could be useful for understanding the dynamics of flow in engines with swirl-stabilized combustion. Graphic abstract: [Figure not available: see fulltext.].

@en

Indirect noise due to the interaction of flow inhomogeneities with a choked nozzle is an important cause of combustion instability in solid rocket motors and is believed to be important in aircraft engines. A previously published experiment (Kings, N., and Bake, F., “Indirect Combustion Noise: Noise Generation by Accelerated Vorticity in a Nozzle Flow,” International Journal of Spray and Combustion Dynamics, Vol. 2, No. 3, 2010, pp. 253–266.) demonstrated that interaction of a nozzle with time-dependent axial swirl can also be a source of sound. This axial swirl was generated by intermittent tangential mass injection upstream from a choked nozzle in a so-called vortex wave generator. The present work discusses the impact of swirl–nozzle interaction in this experiment on the acoustic waves detected downstream of the nozzle. The main source of sound appears to be the reduction in mass flux through the choked nozzle, which depends quadratically on the swirl number. This effect is quantitatively predicted by a quasi-steady and quasi-cylindrical analytical model. The model, combined with empirical data for the decay of axial swirl in pipe flows, predicts the observed influence of the distance between the vortex wave generator and the nozzle. The findings presented here contradict the hypothesis found in the literature, which presumes that sound production in the aforementioned experiment is due to the acceleration of vorticity waves through the nozzle.

@en

Temporary stratospheric aerosol injection (SAI) using sulphate compounds could help to mitigate some of the adverse and irreversible impacts of global warming. Among the risks and uncertainties of SAI, the development of a delivery system presents an appreciable technical challenge. Early studies indicate that specialised aircraft appear the most feasible (McClelan et al., Aurora Flight Sciences, 2010; Smith and Wagner, Environ Res Lett 13(12), 2018). Yet, their technical design characteristics, financial cost of deployment, and emissions have yet to be studied in detail. Therefore, these topics are examined in this two- part study. This first part outlines a set of injection scenarios and proposes a detailed, feasible aircraft design. Part 2 considers the resulting financial cost and equivalent CO₂ emissions spanned by the scenarios and aircraft. Our injection scenarios comprise the direct injection of H₂SO₄ vapour over a range of possible dispersion rates and an SO₂ injection scenario for comparison. To accommodate the extreme demands of delivering large payloads to high altitudes, a coupled optimisation procedure is used to design the system. This results in an unmanned aircraft configuration featuring a large, slender, strut-braced wing and four custom turbofan engines. The aircraft is designed to carry high-temperature H₂SO₄, which is evaporated prior to injection into a single outboard engine plume. Optimised flight profiles are produced for each injection scenario, all involving an initial climb to an outgoing dispersion leg at 20 km altitude, followed by a return dispersion leg at a higher altitude of 20.5 km. All the scenarios considered are found to be technologically and logistically attainable. However, the results demonstrate that achieving high engine plume dispersion rates is of principal importance for containing the scale of SAI delivery systems based on direct H₂SO₄ injection, and to keep these competitive with systems based on SO₂ injection.

@en

To extend the application of the variational multiscale method to ocean engineering, 2D residual-based large eddy simulations of flows around a circular cylinder close to a flat plane have been considered in the present study. A generalized-α method is applied to solving variational forms of the Navier–Stokes equations. Flows around a spanned marine pipeline at a subcritical Reynolds number 1.31×10⁴ have been computed. The hydrodynamic forces are compared to the published experimental data by a PIV and numerical results by a RANS and a classical LES. Although the magnitude of the predicted hydrodynamic forces are over-predicted compared to the experimental data, they are in good agreement with those obtained by the classical LES. The effect of the gap distance between the cylinder and the seabed on the hydrodynamic forces are analyzed. It is found that the drag coefficients and the root-mean-square of the lift coefficients decrease as the gap distance decreases. The flow characteristics around the spanned pipeline are examined. It is demonstrated that the proximity of the cylinder to the seabed severely affects the wake flow of the cylinder, and that the vortex shedding is suppressed as G∕D<0.3, which is in agreement with the experimental and other numerical results. The pressure distributions along the seabed appear to coincide with the experimental data. It is significant that the present method results in more reliable shear stresses (friction velocities) along the seabed than the RANS, which is very important for the research of the sediment transportation.

@en

Indirect noise due to the interaction of flow inhomogeneities with a choked nozzle is as an important cause of combustion instability in solid rocket motors and is believed to be important in aircraft engines. A previously published experiment demonstrated that interaction of the nozzle with time-dependent axial swirl can also be a source of sound. This axial swirl was generated by intermittent circumferential mass injection upstream from a choked nozzle in a so-called Vortex Wave Generator. The present work discusses the impact of swirl-nozzle interaction in this experiment on the acoustic waves detected downstream of the nozzle. The main source of sound appears to be the reduction in mass flux through the choked nozzle, which depends quadratically on the swirl number. This effect is quantitatively predicted by a quasi-steady and quasi-cylindrical analytical model. The model, combined with empirical data for the decay of axial swirl in pipe flows, predicts the observed influence of the distance between the Vortex Wave Generator and the nozzle. The findings presented here contradict the hypothesis found in the literature, which presumes that sound production in the above-mentioned experiment is due to the acceleration of vorticity waves through the nozzle.

@en