Tandem Cylinder Noise Predictions

In an effort to better understand landing-gear noise sources, we have been examining a simplified configuration that still maintains some of the salient features of landing-gear flow fields. In particular, tandem cylinders have been studied because they model a variety of component level interactions. The present effort is directed at the case of two identical cylinders spatially separated in the streamwise direction by 3.7 diameters. Experimental measurements from the Basic Aerodynamic Research Tunnel (BART) and Quiet Flow Facility (QFF) at NASA Langley Research Center (LaRC) have provided steady surface pressures, detailed off-surface measurements of the flow field using Particle Image Velocimetry (PIV), hot-wire measurements in the wake of the rear cylinder, unsteady surface pressure data, and the radiated noise. The experiments were conducted at a Reynolds number of 166 105 based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent shedding process and simulate the effects of a high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The current calculations further explore the influence of the grid resolution and spanwise extent on the flow and associated radiated noise. Extensive comparisons with the experimental data are used to assess the ability of the computations to simulate the details of the flow. The results show that the pressure fluctuations on the upstream cylinder, caused by vortex shedding, are smaller than those generated on the downstream cylinder by wake interaction. Consequently, the downstream cylinder dominates the noise radiation, producing an overall directivity pattern that is similar to that of an isolated cylinder. Only calculations based on the full length of the model span were able to capture the complete decay in the spanwise correlation, thereby producing reasonable noise radiation levels

Aeroacoustic Simulations of Tandem Cylinders with Subcritical Spacing

Tandem cylinders are being studied because they model a variety of component level interactions of landing gear. The present effort is directed at the case of two identical cylinders with their centroids separated in the streamwise direction by 1.435 diameters. Experiments in the Basic Aerodynamic Research Tunnel and Quiet Flow Facility at NASA Langley Research Center have provided an extensive experimental database of the nearfield flow and radiated noise. The measurements were conducted at a Mach number of 0.1285 and Reynolds number of 1.66x10(exp 5) based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent flow separation and, hence, to simulate a major aspect of high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The experiments exhibited an asymmetry in the surface pressure that was persistent despite attempts to eliminate it through small changes in the configuration. To model the asymmetry, the simulations were run with the cylinder configuration at a nonzero but small angle of attack. The computed results and experiments are in general agreement that vortex shedding for the spacing studied herein is weak relative to that observed at supercritical spacings. Although the shedding was subdued in the simulations, it was still more prominent than in the experiments. Overall, the simulation comparisons with measured near-field data and the radiated acoustics are reasonable, especially if one is concerned with capturing the trends relative to larger cylinder spacings. However, the flow details of the 1.435 diameter spacing have not been captured in full even though very fine grid computations have been performed. Some of the discrepancy may be associated with the simulation s inexact representation of the experimental configuration, but numerical and flow modeling errors are also likely contributors to the observed differences

Boundary-Layer Stability Analysis of the Mean Flows Obtained Using Unstructured Grids

Boundary-layer stability analyses of mean flows extracted from unstructured-grid Navier- Stokes solutions have been performed. A procedure has been developed to extract mean flow profiles from the FUN3D unstructured-grid solutions. Extensive code-to-code validations have been performed by comparing the extracted mean ows as well as the corresponding stability characteristics to the predictions based on structured-grid solutions. Comparisons are made on a range of problems from a simple at plate to a full aircraft configuration-a modified Gulfstream-III with a natural laminar flow glove. The future aim of the project is to extend the adjoint-based design capability in FUN3D to include natural laminar flow and laminar flow control by integrating it with boundary-layer stability analysis codes, such as LASTRAC

Impact of Detergent on the Biophysical Properties and Immune Response of the IpaDB Fusion Protein, a Candidate Subunit Vaccine against Shigella spp.

Shigella spp. are causative agents of bacillary dysentery, a human illness with high global morbidity levels, particularly among elderly and infant populations. Shigella infects via the fecal-oral route, and its virulence is dependent upon a type III secretion system (T3SS). Two components of the exposed needle tip complex of the Shigella T3SS, invasion plasmid antigen D (IpaD) and IpaB, have been identified as broadly protective antigens in the mouse lethal pneumonia model. A recombinant fusion protein (DB fusion) was created by joining the coding sequences of IpaD and IpaB. The DB fusion is coexpressed with IpaB\u27s cognate chaperone, IpgC, for proper recombinant expression. The chaperone can then be removed by using the mild detergents octyl oligooxyethelene (OPOE) or N,N-dimethyldodecylamine N-oxide (LDAO). The DB fusion in OPOE or LDAO was used for biophysical characterization and subsequent construction of an empirical phase diagram (EPD). The EPD showed that the DB fusion in OPOE is most stable at neutral pH below 55°C. In contrast, the DB fusion in LDAO exhibited remarkable thermal plasticity, since this detergent prevents the loss of secondary and tertiary structures after thermal unfolding at 90°C, as well as preventing thermally induced aggregation. Moreover, the DB fusion in LDAO induced higher interleukin-17 secretion and provided a higher protective efficacy in a mouse challenge model than did the DB fusion in OPOE. These data indicate that LDAO might introduce plasticity to the protein, promoting thermal resilience and enhanced protective efficacy, which may be important in its use as a subunit vaccine

Spatial linear global instability analysis of the HIFiRE-5 elliptic cone model flow

The linear instability of the three-dimensional boundary-layer over the HIFiRE-5 flight test geometry, i.e. a rounded-tip 2:1 elliptic cone, at Mach 7, has been analyzed through spatial BiGlobal analysis, in a effort to understand transition and accurately predict local heat loads on next-generation ight vehicles. The results at an intermediate axial section of the cone, Re x = 8x10 5, show three different families of spatially amplied linear global modes, the attachment-line and cross- ow modes known from earlier analyses, and a new global mode, peaking in the vicinity of the minor axis of the cone, termed \center-line mode". We discover that a sequence of symmetric and anti-symmetric centerline modes exist and, for the basic ow at hand, are maximally amplied around F* = 130kHz. The wavenumbers and spatial distribution of amplitude functions of the centerline modes are documente

The Effect of Acoustic Forcing on Instabilities and Breakdown in Swept-Wing Flow over a Backward-Facing Step

Instability interaction and breakdown were experimentally investigated in the flow over a swept backward-facing step. Acoustic forcing was used to excite the Tollmien-Schlichting (TS) instability and to acquire phase-locked results. The phase-averaged results illustrate the complex nature of the interaction between the TS and stationary cross flow instabilities. The weak stationary cross flow disturbance causes a distortion of the TS wavefront. The breakdown process is characterized by large positive and negative spikes in velocity. The positive spikes occur near the same time and location as the positive part of the TS wave. Higher-order spectral analysis was used to further investigate the nonlinear interactions between the TS instability and the traveling cross flow disturbances. The results reveal that a likely cause for the generation of the spikes corresponds to nonlinear interactions between the TS, traveling cross flow, and stationary cross flow disturbances. The spikes begin at low amplitudes of the unsteady and steady disturbances (2-4% U (sub e) (i.e. boundary layer edge velocity)) but can achieve very large amplitudes (20-30 percent U (sub e) (i.e. boundary layer edge velocity)) that initiate an early, though highly intermittent, breakdown to turbulence

Maximizing Happiness in Graphs of Bounded Clique-Width

Clique-width is one of the most important parameters that describes structural complexity of a graph. Probably, only treewidth is more studied graph width parameter. In this paper we study how clique-width influences the complexity of the Maximum Happy Vertices (MHV) and Maximum Happy Edges (MHE) problems. We answer a question of Choudhari and Reddy '18 about parameterization by the distance to threshold graphs by showing that MHE is NP-complete on threshold graphs. Hence, it is not even in XP when parameterized by clique-width, since threshold graphs have clique-width at most two. As a complement for this result we provide a $n^{\mathcal{O}(\ell \cdot \operatorname{cw})}$ algorithm for MHE, where $\ell$ is the number of colors and $\operatorname{cw}$ is the clique-width of the input graph. We also construct an FPT algorithm for MHV with running time $\mathcal{O}^*((\ell+1)^{\mathcal{O}(\operatorname{cw})})$, where $\ell$ is the number of colors in the input. Additionally, we show $\mathcal{O}(\ell n^2)$ algorithm for MHV on interval graphs.Comment: Accepted to LATIN 202

Closed-loop control of boundary layer streaks induced by free-stream turbulence

The central aim of the paper is to carry out a theoretical and numerical study of active wall transpiration control of streaks generated within an incompressible boundary layer by free-stream turbulence. The disturbance flow model is based on the linearized unsteady boundary-region (LUBR) equations, studied by Leib, Wundrow, and Goldstein [J. Fluid Mech. 380, 169 (1999)], which are the rigorous asymptotic limit of the Navier-Stokes equations for low-frequency and long-streamwise wavelength. The mathematical formulation of the problem directly incorporates the random forcing into the equations in a consistent way. Due to linearity, this forcing is factored out and appears as a multiplicative factor. It is shown that the cost function (integral of kinetic energy in the domain) is properly defined as the expectation of a random quadratic function only after integration in wave number space. This operation naturally introduces the free-stream turbulence spectral tensor into the cost function. The controller gains for each wave number are independent of the spectral tensor and, in that sense, universal. Asymptotic matching of the LUBR equations with the free-stream conditions results in an additional forcing term in the state-space system whose presence necessitates the reformulation of the control problem and the rederivation of its solution. It is proved that the solution can be obtained analytically using an extension of the sweep method used in control theory to obtain the standard Riccati equation. The control signal consists of two components, a feedback part and a feed-forward part (that depends explicitly on the forcing term). Explicit recursive equations that provide these two components are derived. It is shown that the feed-forward part makes a negligible contribution to the control signal. We also derive an explicit expression that a priori (i.e., before solving the control problem) leads to the minimum of the objective cost function (i.e., the fundamental performance limit), based only on the system matrices and the initial and free-stream boundary conditions. The adjoint equations admit a self-similar solution for large spanwise wave numbers with a scaling which is different from that of the LUBR equations. The controlled flow field also has a self-similar solution if the weighting matrices of the objective function are chosen appropriately. The code developed to implement this algorithm is efficient and has modest memory requirements. Computations show the significant reduction of energy for each wave number. The control of the full spectrum streaks, for conditions corresponding to a realistic experimental case, shows that the root-mean-square of the streamwise velocity is strongly suppressed in the whole domain and for all the frequency ranges examined

Numerical Evaluation of P-Multigrid Method for the Solution of Discontinuous Galerkin Discretizations of Diffusive Equations

This paper describes numerical experiments with P-multigrid to corroborate analysis, validate the present implementation, and to examine issues that arise in the implementations of the various combinations of relaxation schemes, discretizations and P-multigrid methods. The two approaches to implement P-multigrid presented here are equivalent for most high-order discretization methods such as spectral element, SUPG, and discontinuous Galerkin applied to advection; however it is discovered that the approach that mimics the common geometric multigrid implementation is less robust, and frequently unstable when applied to discontinuous Galerkin discretizations of di usion. Gauss-Seidel relaxation converges 40% faster than block Jacobi, as predicted by analysis; however, the implementation of Gauss-Seidel is considerably more expensive that one would expect because gradients in most neighboring elements must be updated. A compromise quasi Gauss-Seidel relaxation method that evaluates the gradient in each element twice per iteration converges at rates similar to those predicted for true Gauss-Seidel