Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Anisotropic flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Flow Physics at a Large Rotor Tip Gap in a Multi-Stage Axial Compressor

The flow physics in a large rotor tip gap in a 1.5-stage axial compressor is investigated in the current study. The flow structure in the rotor tip region is complex with several dominant vortical structures of opposite rotation, resulting in inhomogeneous and highly anisotropic turbulence. Earlier measurements show that eddy viscosity is negative over large parts of the tip region and eddy viscosity varies among stress/strain components. The present study aims to understand how the complex nature of rotor tip leakage flow affects compressor performance when the tip gap size is greater than 4–5% of the rotor span, which is typical of advanced small core engines. Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) techniques are applied to study flow physics in a large rotor tip gap (5.5% of rotor span) in a 1.5-stage axial compressor. Calculated flow fields from the two different approaches are compared with available measurements and examined in detail. LES calculates the pressure rise in the present compressor fairly well, while URANS with a standard two-equation turbulence closure underpredicts the pressure rise by 15-20% of the measured values. The current study shows that URANS with the current turbulence closure produces much higher all-positive eddy viscosity in the tip-gap region compared to measurements and LES. The distribution of eddy viscosity in the URANS simulation is also wrong. Consequently, the flow in the tip region is highly damped with significantly larger blockage generation, which results in the tip leakage vortex (TLV) staying closer to the blade suction side compared to the measurement. When the TLV stays closer to the blade, both flow turning and the pressure rise across the compressor are reduced compared to the measurements. It appears that this effect is amplified by a large rotor tip gap.

tip gap↗

Optimizing microfluidic flow cell geometry for in situ resonant soft X-ray characterization of molecular nanostructures

Liquid-phase resonant soft X-ray scattering (LP-RSoXS) is an emerging label-free technique to probe chemically resolved nanostructures of molecular or hybrid materials in liquid environments. Still, quantitative analysis is hindered by the pressure-induced deformation of thin silicon nitride (SiN) membranes used as windows in microfluidic flow cells, which attenuates the signal in nonlinear ways, making experimental optimization difficult. Here, in this work, we directly characterize this deformation under experimental conditions for a variety of cell configurations. We use this to develop a predictive model that combines transmission effects of SiN bowing, incident X-ray beam profiles, and material-dependent resonant scattering cross sections to simulate the effective scattering intensity at the detector across the carbon K-edge. Maps of the total signal across the flow cell window reveal that increasing the window width and polymer concentration shifts the anisotropic intensity distributions from the center toward the edges of the window. It was determined that an optimal SiN thickness of 50 nm, with a window aperture of 104 μm, maximizes the total signal for typical solute concentrations and energies across the carbon K-edge. Our results overturn the assumption that corner regions dominate the scattering signal, offering explicit design guidelines for maximizing LP-RSoXS signals and significantly advancing the quantitative application of this technique to the characterization of molecular and hybrid nanostructured materials in liquids.

Grabner, Devin [Washington State Univ., Pullman, W↗

Effect of anisotropic turbulence on aerodynamic noise

A model based on Lighthill's theory for predicting aerodynamic noise from a turbulent shear flow is developed. This model is a generalization of the one developed by Ribner. It does not require that the turbulent correlations factor into space and time-dependent parts. It replaces his assumption of isotropic turbulence by the more realistic one of axisymmetric turbulence. In the course of the analysis, a hierarchy of equations is developed wherein each succeeding equation involves more assumptions than the preceding equation but requires less experimental information for its use. The implications of the model for jet noise are discussed. It is shown that for the particular turbulence data considered anisotropy causes the high-frequency self-noise to be beamed downstream.

Goldstein, M.↗

Multiple-scale turbulence model in confined swirling jet predictions

A recently developed multiple-scale turbulence model which attempts to circumvent the deficiencies of earlier models by taking nonequilibrium spectral energy transfer into account is presented. The model's validity is tested by predicting the confined swirling coaxial jet flow in a sudden expansion. It is noted that, in order to account for anisotropic turbulence, a full Reynolds stress model is required.

Chen, C. P.↗

Tetrahedral and Hexahedral Mesh Adaptation for CFD Problems

This paper presents two unstructured mesh adaptation schemes for problems in computational fluid dynamics. The procedures allow localized grid refinement and coarsening to efficiently capture aerodynamic flow features of interest. The first procedure is for purely tetrahedral grids; unfortunately, repeated anisotropic adaptation may significantly deteriorate the quality of the mesh. Hexahedral elements, on the other hand, can be subdivided anisotropically without mesh quality problems. Furthermore, hexahedral meshes yield more accurate solutions than their tetrahedral counterparts for the same number of edges. Both the tetrahedral and hexahedral mesh adaptation procedures use edge-based data structures that facilitate efficient subdivision by allowing individual edges to be marked for refinement or coarsening. However, for hexahedral adaptation, pyramids, prisms, and tetrahedra are used as buffer elements between refined and unrefined regions to eliminate hanging vertices. Computational results indicate that the hexahedral adaptation procedure is a viable alternative to adaptive tetrahedral schemes.

Biswas, Rupak↗

Sketch-to-Solution: A Case Study in RCS Aerodynamic Interaction

Thanks to recent advances in the fields of anisotropic grid adaptation, error estimation, and geometry modeling, a sketch-to-solution work flow is now possible for viscous computational fluid dynamic (CFD) simulations. With this workflow, a CFD application engineer provides geometry, boundary conditions, and flow parameters; and the sketch-to-solution process yields a CFD simulation through automatic, error-based, grid adaptation. To explore the benefits of this nascent capability, a conventional manual grid generation work flow is compared to this new automatic grid generation work flow for a given engineering question: What are the aerodynamic interactions caused by the reaction control system (RCS) on an entry vehicle? This case study indicates that while the automatic grid generation sketch-to- solution process is not yet mature, it is preferred over a manual grid generation work flow because it greatly reduces manual labor, eliminates many opportunities for human error, and provides grid sensitivity information.

Bill Kleb↗

Three-dimensional turbulent-mixing-length modeling for discrete-hole coolant injection into a crossflow

Three dimensional mixing length models of a flow field immediately downstream of coolant injection through a discrete circular hole at a 30 deg angle into a crossflow were derived from the measurements of turbulence intensity. To verify their effectiveness, the models were used to estimate the anisotropic turbulent effects in a simplified theoretical and numerical analysis to compute the velocity and temperature fields. With small coolant injection mass flow rate and constant surface temperature, numerical results of the local crossflow streamwise velocity component and surface heat transfer rate are consistent with the velocity measurement and the surface film cooling effectiveness distributions reported in previous studies.

Wang, C. R.↗

Computation of supersonic and low subsonic cascade flows using an explicit Navier-Stokes technique and the kappa-epsilon turbulence model

A fully explicit two-dimensional flow solver, based on a four-stage Runge-Kutta scheme, was developed and utilized to predict two-dimensional viscous flow through turbomachinery cascades for which experimental data is available. The formulation is applied to the density averaged Navier-Stokes equations. Several features of the technique improve the ability of the code to predict high Reynolds number flows on highly stretched grids. These include a low Reynolds number compressible form of the k-epsilon turbulence model, anisotropic scaling of artificial dissipation terms and locally varying timestep evaluation based on hyperbolic and parabolic stability considerations. Comparisons between computation and experiment are presented for both a supersonic and a low-subsonic compressor cascade. These results indicate that the code is capable of predicting steady two-dimensional viscous cascade flows over a wide range of Mach numbers in reasonable computational times.

Kunz, R. F.↗

Theory and observation of energetic ions in the lobes of the geomagnetic tail

The Owen et al. (1990) model, which attempts to explain the nature of the pitch angle distributions of energetic ions within the lobes of the distant geomagnetic tail is briefly reviewed. Energetic ion data from the ISEE-3 spacecraft obtained during early 1983, when the spacecraft made several traversals of the distant geomagnetic tail are then presented. The data demonstrate that during quiet periods in which the spacecraft is continuously located in the tail lobes, the pitch angle distribution is observed to be highly anisotropic, being peaked closely perpendicular to the magnetic field direction, but with a small net flow in the antisunward direction, in agreement with the model results. Further predictions of the model, concerning the variation of the lobe energetic ion distributions with position in the lobes, are compared to observations made as the spacecraft performed a traversal of the lobe. Finally, since the model indicates that a more isotropic distribution should exist in the tail lobe during solar particle events, data is presented from such a period for further comparison. In all the above cases, good agreement is demonstrated between the data and the expections of the model.

Owen, C. J.↗

Failure of Anisotropic Unstructured Mesh Adaption Based on Multidimensional Residual Minimization

An automated anisotropic unstructured mesh adaptation strategy is proposed, implemented, and assessed for the discretization of viscous flows. The adaption criteria is based upon the minimization of the residual fluctuations of a multidimensional upwind viscous flow solver. For scalar advection, this adaption strategy has been shown to use fewer grid points than gradient based adaption, naturally aligning mesh edges with discontinuities and characteristic lines. The adaption utilizes a compact stencil and is local in scope, with four fundamental operations: point insertion, point deletion, edge swapping, and nodal displacement. Evaluation of the solution-adaptive strategy is performed for a two-dimensional blunt body laminar wind tunnel case at Mach 10. The results demonstrate that the strategy suffers from a lack of robustness, particularly with regard to alignment of the bow shock in the vicinity of the stagnation streamline. In general, constraining the adaption to such a degree as to maintain robustness results in negligible improvement to the solution. Because the present method fails to consistently or significantly improve the flow solution, it is rejected in favor of simple uniform mesh refinement.

Wood, William A.↗

The case for anisotropic mass loss from T Tauri stars

Observational evidence and theoretical arguments are assembled in favor of directed, rather than isotropic, mass flows from T Tauri stars. These arguments include the structure of T Tauri nebulosities, radio continuum maps of the stars, the motions of Herbig-Haro objects away from their exciting T Tauri stars, bipolar CO flows, and particularly the spatial coincidence of an optical jet and radio emission associated with an embedded, probably T Tauri star in the dark cloud L1551.

Cohen, M.↗

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗

Anisotropic ion heating and parallel O(+) acceleration in regions of rapid E x B convection

A numerical solution to the 20-moment set of transport equations has been found in order to study subauroral ionospheric outflows during periods of enhanced perpendicular ion drifts. The numerical model solves the time-dependent O(+) density, momentum, and both the parallel and perpendicular energy and heat flow equations in the 200-6000 km altitude range. Assuming perpendicular drifts of 3 km/s relative to the neutral atmosphere, we have found that anisotropic heating of O(+) (a result of ion-neutral collisions) leads to a temperature anisotropy, with perpendicular temperatures exceeding 8000 K and parallel temperatures greater than 5000 K (near 200 km altitude). Above approximately 2000 km, transport processes dominate the effects of collisions and wavelike oscillations in O(+) velocity, temperature and heat flux were noted.

Korosmezey, A.↗

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.↗

Noise produced by turbulent flow into a rotor: Theory manual for atmospheric turbulence prediction and mean flow and turbulence contraction prediction

Prediction of helicopter main rotor noise due to ingestion of atmospheric turbulence was analyzed. The analysis combines several different models that describe the fluid mechanics of the turbulence and the ingestion process. Two models, atmospheric turbulence, and mean flow and turbulence contraction were covered. The third model, covered in a separate report, describes the rotor acoustic mode. The method incorporates the atmospheric turbulence model and a rapid distortion turbulence contraction description to determine the statistics of the anisotropic turbulence at the rotor plane. The analytical basis for a module was provided which was incorporated in NASA's ROTONET helicopter noise prediction program. The mean flow and turbulence statistics associated with the atmospheric boundary layer were modeled including effects of atmospheric stability length, wind speed, and altitude. The turbulence distortion process is modeled as a deformation of vortex filaments (which represent the turbulence field) by a mean flow field due to the rotor inflow.

Simonich, J. C.↗

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗

Self-Diffusion of Drops in a Dilute Sheared Emulsion

Self-diffusion coefficients that describe cross-flow migration of non-Brownian drops in a dilute sheared emulsion were obtained by trajectory calculations. A boundary integral formulation was used to describe pairwise interactions between deformable drops; interactions between undeformed drops were described with mobility functions for spherical drops. The results indicate that drops have large anisotropic self-diffusivities which depend strongly on the drop viscosity and modestly on the shear-rate. Pairwise interactions between drops in shear-flow do not appreciably promote drop breakup.

Loewenberg, Michael↗