The geometry was tested in the Sandia Trisonic Wind Tunnel (TWT). The geometry is an insert into the floor/ceiling of the wind tunnel. It consists of a cavity (rectangular cut-out) with some “complex” features. The features and dimensions are presented here. Wall pressure data has already been presented at International open conferences. This is being used for validation purposes by international researchers. The purpose of this document is to make the geometry available for release to such groups.
We demonstrate a Bayesian method that can be used to calibrate computationally expensive 3D RANS models with complex response surfaces. Such calibrations, conditioned on experimental data, can yield turbulence model parameters as probability density functions (PDF), concisely capturing the uncertainty in the estimation. Methods such as Markov chain Monte Carlo construct the PDF by sampling, and consequently a quick-running surrogate is used instead of the RANS simulator. The surrogate can be very difficult to design if the model’s response i.e., the dependence of the calibration variable (the observable) on the parameters being estimated is complex. We show how the training data used to construct the surrogate models can also be employed to isolate a promising and physically realistic part of the parameter space, within which the response is well-behaved and easily modeled. We design a classifier, based on treed linear models, to model the “well-behaved region”. This classifier serves as a prior in a Bayesian calibration study aimed at estimating 3 k-ε parameters C = (Cμ, Cε2, Cε1) from experimental data of a transonic jet-in-crossflow interaction. The robustness of the calibration is investigated by checking its predictions of variables not included in the calibration data. We also check the limit of applicability of the calibration by testing at an off-calibration point.
The flow over aircraft bays exhibits many characteristics of cavity flows, namely resonant pressures that can create high structural loading. An extensive dataset of pressure measurements within both simple and complex cavities was previously obtained and analyzed using power-spectral densities, coherence levels, and cross correlations between sensor pairs within the cavity. More in-depth analysis of the flow structure is studied here using modal decomposition techniques. Both Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) were applied to the experimental and computational results within a simple rectangular cavity. POD was able to show that the cavity modes are coherent across the cavity width. Only higher modes that were associated with more turbulent fluctuations exhibited spanwise variations. These were concentrated at the aft end of the cavity. DMD was able to isolate structures associated with single frequencies in the flow. At the Rossiter frequencies, coherent structures across the front cavity width were found, while more complex shapes were observed at the cavity rear, consistent with the POD analysis. Additional DMD modes in between the dominant Rossiter frequencies also appeared. These additional modes were associated with a low-frequency modulation of the cavity tones.
This work examines simulation requirements for ensuring accurate predictions of compressible cavity flows. Lessons learned from this study will be used in the future to study the effects of complex geometric features, representative of those found on real weapons bays, on compressible flow past open cavities. A hybrid RANS/LES simulation method is applied to a rectangular cavity with length-to-depth ratio of 7, in order to first validate the model for this class of flows. Detailed studies of mesh resolution, absorbing boundary condition formulation, and boundary zone extent are included and guidelines are developed for ensuring accurate prediction of cavity pressure fluctuations.
This work examines simulation requirements for ensuring accurate predictions of compressible cavity flows. Lessons learned from this study will be used in the future to study the effects of complex geometric features, representative of those found on real weapons bays, on compressible flow past open cavities. A hybrid RANS/LES simulation method is applied to a rectangular cavity with length-to-depth ratio of 7, in order to first validate the model for this class of flows. Detailed studies of mesh resolution, absorbing boundary condition formulation, and boundary zone extent are included and guidelines are developed for ensuring accurate prediction of cavity pressure fluctuations.
An approach for building energy-stable Galerkin reduced order models (ROMs) for linear hyperbolic or incompletely parabolic systems of partial differential equations (PDEs) using continuous projection is developed. This method is an extension of earlier work by the authors specific to the equations of linearized compressible inviscid flow. The key idea is to apply to the PDEs a transformation induced by the Lyapunov function for the system, and to build the ROM in the transformed variables. For linear problems, the desired transformation is induced by a special inner product, termed the "symmetry inner product", which is derived herein for several systems of physical interest. Connections are established between the proposed approach and other stability-preserving model reduction methods, giving the paper a review flavor. More specifically, it is shown that a discrete counterpart of this inner product is a weighted L2 inner product obtained by solving a Lyapunov equation, first proposed by Rowley et al. and termed herein the "Lyapunov inner product". Comparisons between the symmetry inner product and the Lyapunov inner product are made, and the performance of ROMs constructed using these inner products is evaluated on several benchmark test cases.
We demonstrate a Bayesian method that can be used to calibrate computationally expensive 3D RANS (Reynolds Av- eraged Navier Stokes) models with complex response surfaces. Such calibrations, conditioned on experimental data, can yield turbulence model parameters as probability density functions (PDF), concisely capturing the uncertainty in the parameter estimates. Methods such as Markov chain Monte Carlo (MCMC) estimate the PDF by sampling, with each sample requiring a run of the RANS model. Consequently a quick-running surrogate is used instead to the RANS simulator. The surrogate can be very difficult to design if the model's response i.e., the dependence of the calibration variable (the observable) on the parameter being estimated is complex. We show how the training data used to construct the surrogate can be employed to isolate a promising and physically realistic part of the parameter space, within which the response is well-behaved and easily modeled. We design a classifier, based on treed linear models, to model the "well-behaved region". This classifier serves as a prior in a Bayesian calibration study aimed at estimating 3 k - ε parameters ( C μ, C ε2 , C ε1 ) from experimental data of a transonic jet-in-crossflow interaction. The robustness of the calibration is investigated by checking its predictions of variables not included in the cal- ibration data. We also check the limit of applicability of the calibration by testing at off-calibration flow regimes. We find that calibration yield turbulence model parameters which predict the flowfield far better than when the nomi- nal values of the parameters are used. Substantial improvements are still obtained when we use the calibrated RANS model to predict jet-in-crossflow at Mach numbers and jet strengths quite different from those used to generate the ex- perimental (calibration) data. Thus the primary reason for poor predictive skill of RANS, when using nominal values of the turbulence model parameters, was parametric uncertainty, which was rectified by calibration. Post-calibration, the dominant contribution to model inaccuraries are due to the structural errors in RANS.
The datasets being released consist of cavity configurations for which measurements were made in the Sandia Trisonic Wind Tunnel (TWT) facility. The cavities were mounted on the walls (ceiling/floor) of the wind tunnel, with the approach flow boundary layer thickness dictated by the run-length from the settling chamber of the tunnel. No measurements of the boundary layer for the different cases were made explicitly. However, prior measurements of the boundary layer have been made and simulations of the tunnel from the settling chamber on have shown that this method yields the correct boundary layer thickness at the leading edge of the cavity. The measurements focused on the cavity flow field itself and the cavity wall pressures. For each of the cases, the stagnation conditions are prescribed in order to obtain the correct inflow conditions upstream of the cavity. The wind tunnel contours have been approved for public release and will be made available also.
This report describes work performed from June 2012 through May 2014 as a part of a Sandia Early Career Laboratory Directed Research and Development (LDRD) project led by the first author. The objective of the project is to investigate methods for building stable and efficient proper orthogonal decomposition (POD)/Galerkin reduced order models (ROMs): models derived from a sequence of high-fidelity simulations but having a much lower computational cost. Since they are, by construction, small and fast, ROMs can enable real-time simulations of complex systems for onthe- spot analysis, control and decision-making in the presence of uncertainty. Of particular interest to Sandia is the use of ROMs for the quantification of the compressible captive-carry environment, simulated for the design and qualification of nuclear weapons systems. It is an unfortunate reality that many ROM techniques are computationally intractable or lack an a priori stability guarantee for compressible flows. For this reason, this LDRD project focuses on the development of techniques for building provably stable projection-based ROMs. Model reduction approaches based on continuous as well as discrete projection are considered. In the first part of this report, an approach for building energy-stable Galerkin ROMs for linear hyperbolic or incompletely parabolic systems of partial differential equations (PDEs) using continuous projection is developed. The key idea is to apply a transformation induced by the Lyapunov function for the system, and to build the ROM in the transformed variables. It is shown that, for many PDE systems including the linearized compressible Euler and linearized compressible Navier-Stokes equations, the desired transformation is induced by a special inner product, termed the “symmetry inner product”. Attention is then turned to nonlinear conservation laws. A new transformation and corresponding energy-based inner product for the full nonlinear compressible Navier-Stokes equations is derived, and it is demonstrated that if a Galerkin ROM is constructed in this inner product, the ROM system energy will be bounded in a way that is consistent with the behavior of the exact solution to these PDEs, i.e., the ROM will be energy-stable. The viability of the linear as well as nonlinear continuous projection model reduction approaches developed as a part of this project is evaluated on several test cases, including the cavity configuration of interest in the targeted application area. In the second part of this report, some POD/Galerkin approaches for building stable ROMs using discrete projection are explored. It is shown that, for generic linear time-invariant (LTI) systems, a discrete counterpart of the continuous symmetry inner product is a weighted L2 inner product obtained by solving a Lyapunov equation. This inner product was first proposed by Rowley et al., and is termed herein the “Lyapunov inner product“. Comparisons between the symmetry inner product and the Lyapunov inner product are made, and the performance of ROMs constructed using these inner products is evaluated on several benchmark test cases. Also in the second part of this report, a new ROM stabilization approach, termed “ROM stabilization via optimization-based eigenvalue reassignment“, is developed for generic LTI systems. At the heart of this method is a constrained nonlinear least-squares optimization problem that is formulated and solved numerically to ensure accuracy of the stabilized ROM. Numerical studies reveal that the optimization problem is computationally inexpensive to solve, and that the new stabilization approach delivers ROMs that are stable as well as accurate. Summaries of “lessons learned“ and perspectives for future work motivated by this LDRD project are provided at the end of each of the two main chapters.
We propose a Bayesian method to calibrate parameters of a k-∈ RANS model to improve its predictive skill in jet-in-crossflow simulations. The method is based on the hypotheses that (1) informative parameters can be estimated from experiments of flow configurations that display the same, strongly vortical features of jet-in-crossflow interactions and (2) one can construct surrogates of RANS models for judiciously chosen outputs which serve as calibration observables. We estimate three k - ∈ parameters, (Cμ,C∈2,C∈1), from Reynolds stress measurements obtained from an incompressible flow-over-a-square-cylinder experiment. The k - ∈ parameters are estimated as a joint probability density function. Jet-in-crossflow simulations performed with (Cμ,C∈2,C∈) samples drawn from this distribution are seen to provide far better predictions than those obtained with nominal parameter values. We also find a (Cμ,C∈2,C∈1) combination which provides less than 15% error in a number of performance metrics. In contrast, the errors obtained with nominal parameter values may exceed 60%.
Simulations of a rectangular cavity containing a model captive store are performed using a Hybrid Reynolds-averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) model. The fluid flow simulations are coupled to a structural dynamics finite element model using a one-way pressure transfer procedure. Simulation results for pressure fluctuation spectra and store acceleration are compared to measurements made on the same configuration in a tri-sonic wind tunnel at Mach numbers of 0.60, 0.80, and 1.47. The simulation results are used to calculate unsteady integrated forces and moments acting on the store. Spectra of the forces and moments reveal that a complex relationship exists between the unsteady integrated forces and the measured resonant cavity modes as indicated in the cavity wall pressure measurements. Predictions of the store accelerations from the coupled model show some success in predicting both forced and natural modal responses of the store within the cavity environment, while also highlighting some challenges in obtaining statistically converged results for this class of problems.
A newly-developed computational fluid-structure interaction framework for simulation of stores in captive carriage environments is validated. The computational method involves one-way coupling, with pressure loads calculated by a hybrid RANS-LES CFD model transferred to a structural dynamics solver. Validation is performed at several levels. First, the ability of the CFD model to accurately predict the flow-field and resulting aerodynamic loads in an empty cavity is assessed against wind tunnel data. In parallel, the structural dynamics model for a simulated store is calibrated and then validated against a shaker table experiment. Finally, predictions of aerodynamic loads and store vibrations from the coupled simulation model are compared to new wind tunnel experimental data for a model captive carriage configuration.
We develop a novel calibration approach to address the problem of predictive ke RANS simulations of jet-incrossflow. Our approach is based on the hypothesis that predictive ke parameters can be obtained by estimating them from a strongly vortical flow, specifically, flow over a square cylinder. In this study, we estimate three ke parameters, C%CE%BC, Ce2 and Ce1 by fitting 2D RANS simulations to experimental data. We use polynomial surrogates of 2D RANS for this purpose. We conduct an ensemble of 2D RANS runs using samples of (C%CE%BC;Ce2;Ce1) and regress Reynolds stresses to the samples using a simple polynomial. We then use this surrogate of the 2D RANS model to infer a joint distribution for the ke parameters by solving a Bayesian inverse problem, conditioned on the experimental data. The calibrated (C%CE%BC;Ce2;Ce1) distribution is used to seed an ensemble of 3D jet-in-crossflow simulations. We compare the ensemble's predictions of the flowfield, at two planes, to PIV measurements and estimate the predictive skill of the calibrated 3D RANS model. We also compare it against 3D RANS predictions using the nominal (uncalibrated) values of (C%CE%BC;Ce2;Ce1), and find that calibration delivers a significant improvement to the predictive skill of the 3D RANS model. We repeat the calibration using surrogate models based on kriging and find that the calibration, based on these more accurate models, is not much better that those obtained with simple polynomial surrogates. We discuss the reasons for this rather surprising outcome.