This work extends our previous study from S. Shrestha et al. (2024) by introducing a new abstract framework for Variational Multiscale (VMS) methods at the discrete level. We introduce the concept of what we define as the optimal projector and present a discretisation approach that yields a numerical solution closely approximating the optimal projection of the infinite-dimensional continuous solution. In this approach, the infinite-dimensional unresolved scales are approximated in a finite-dimensional subspace using the numerically computed Fine-Scale Greens’ function of the underlying symmetric problem. The proposed approach involves solving the VMS problem on two separate meshes: a coarse mesh for the full PDE and a fine mesh for the symmetric part of the continuous differential operator. We consider the 1D and 2D steady advection–diffusion problems in both direct and mixed formulations as the test cases in this paper. We first present an error analysis of the proposed approach and show that the projected solution is achieved as the approximate Greens’ function converges to the exact one. Subsequently, we demonstrate the working of this method where we show that it can exponentially converge to the chosen optimal projection. We note that the implementation of the present work employs the Mimetic Spectral Element Method (MSEM), although, it may be applied to other Finite/Spectral Element or Isogeometric frameworks. Furthermore, we propose that VMS should not be viewed as a stabilisation technique; instead, the base scheme should be inherently stable, with VMS enhancing the solution quality by supplementing the base scheme.

@en

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

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

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

Efficient input data generation for reduced-order model applications to accurately predict aerodynamic performance and stability characteristics over a large part of a fighter aircraft’s flight envelope is a major challenge. In this paper, aerodynamic reduced-order models are created from two pseudorandom binary sequence (PRBS) training maneuvers. During these maneuvers, the angle of attack and pitch rate change in a periodic and deterministic manner which is characterized by white-noise-like properties. Typical PRBS signals include sudden input variations between two distinct values, such as minimum and maximum angles of attack. However, the signals used in this paper were modified to have the step changes to depend on the simulation time. In the first motion, the aircraft undergoes a signal at a constant Mach number of 0.85. In the second motion, the Mach number varies in an optimized manner from 0.1 to 0.9. The test case is a generic triple-delta wing configuration. Simulations were run using the DoD HPCMP CREATE^RM-AV/Kestrel simulation tools. A prescribed-body motion was used to vary input parameters under given freestream conditions (Mach number and angle of attack). Different reduced-order methods were applied, that comprise regression, feed-forward neural network and auto-regressive surrogate modeling techniques to predict integrated force and moment coefficients and a proper-orthogonal decomposition based neural network approach for surface pressure prediction. Once models of integrated forces and moments were created, they were used to predict static and stability derivatives at different angles of attack and Mach numbers. Models were then used to predict aerodynamic responses to arbitrary motions including pitch sinusoidal, chirp, Schroeder, and step. Model predictions were compared with actual CFD data. Overall, a good agreement was found for all models. Models to predict surface pressure data were also able to accurately predict the upper surface pressure data at different spanwise and chordwise locations at different angles of attack for both static and dynamic runs.

@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

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

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

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

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

Entropy inhomogeneities and vorticity spots induce so-called indirect combustion noise when passing through a choked nozzle; referred to as entropy noise and vorticity noise, respectively. We note that vorticity noise depends on the orientation of the vorticity; viz., oriented normal or parallel to the axial main flow. An experimental investigation of parallel component vorticity noise is presented. In the experiment a time-dependent swirling flow was induced by unsteady tangential injection in the pipe upstream of a choked convergent-divergent nozzle. As the resulting swirling flow passes through the nozzle, the axial stretching of the fluid caused an increase in rotation energy. The steady energy conservation in an isentropic flow implies a Mach number higher than unity at the throat and an associated reduction of density. Ergo, the critical mass-flow rate (for fixed reservoir pressure and temperature) decreases quadratically with increasing swirl intensity. The acoustic waves radiated downstream of the nozzle are due to the change in the mass flow through the nozzle. These are a direct measure for this mass-flow modulation. Using a semi-empirical model, this sound production mechanism is demonstrated to be quasi steady.

@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

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

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

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

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