SEARCH · Search NASA
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.
Spherical bubble growth or collapse in component inviscid liquid with constant material properties through irreversible thermodynamics and kinetic theory
Explore the source record for details and available documents.
In this research, a near-wall second-order closure based on the Speziable et al.(1991) or SSG model for the pressure-strain term is proposed. Unlike the LRR model, the SSG model is quasi-nonlinear and yields better results when applied to calculate rotating homogeneous turbulent flows. An asymptotic analysis near the wall is applied to both the exact and modeled, equations so that appropriate near-wall corrections to the SSG model and the modeled dissipation-rate equation can be derived to satisfy the physical wall boundary conditions as well as the asymptotic near-wall behavior of the exact equations. Two additional model constants are introduced and they are determined by calibrating against one set of near-wall channel flow data. Once determined, their values are found to remain constant irrespective of the type of flow examined. The resultant model is used to calculate simple turbulent flows, near separating turbulent flows, complex turbulent flows and compressible turbulent flows with a freestream Mach number as high as 10. In all the flow cases investigated, the calculated results are in good agreement with data. This new near-wall model is less ad hoc, physically and mathematically more sound and eliminates the empiricism introduced by Zhang. Therefore, it is quite general, as demonstrated by the good agreement achieved with measurements covering a wide range of Reynolds numbers and Mach numbers.
The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined circular duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions at each axial location. These values are used to determine the growth rates and phase velocities of the instability waves and the near field pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, and wall impedance, and the effects of these changes on non-axisymmetric modes. The primary results indicate that the effects of the duct wall on stability characteristics diminish as the outer flow increases and as the jet azimuthal mode number increases. Also, wall reflections are reduced when using a finite impedance boundary condition at the wall; but in addition, reflections are reduced and growth rates diminished by keeping the imaginary part of the impedance negative when using the negative exponential for the harmonic dependence.
A theoretical examination of the slotted-wall flow field is conducted to determine the appropriate wall pressure drop (or boundary condition) equation. This analysis improves the understanding of the fluid physics of these types of flow fields and helps in evaluating the uncertainties and limitations existing in previous mathematical developments. It is shown that the resulting slotted-wall boundary condition contains contributions from the airfoil-induced streamline curvature and the non-linear, quadratic, slot crossflow in addition to an often neglected linear term which results from viscous shearing in the slot. Existing and newly acquired experimental data are examined in the light of this formulation and theoretical developments.
In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.
Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.
The original aim of this work was to shed some light on the physics of turbulence over rough walls using large-eddy simulations and the suboptimal-control wall boundary conditions introduced by Nicoud et al. It was hoped that, if that algorithm was used to fit the mean velocity profile of the simulations to that of a rough-walled channel, instead of to a smooth one, the wall stresses introduced by the control algorithm would give some indication of what aspects of rough walls are most responsible for the modification of the flow in real turbulence. It was similarly expected that the structure of the resulting velocity fluctuations would share some of the characteristics of rough-walled flows, thus again suggesting what is intrinsic and what is accidental in the effect of geometric wall roughness. A secondary goal was to study the effect of 'unphysical' boundary conditions on the outside flow by observing how a relatively major change of the target velocity profile, and therefore presumably of the applied wall stresses, modifies properties such as the dominant length scales of the velocity fluctuations away from the wall. As will be seen below, this secondary goal grew more important during the course of the study, which was carried out during a short summer visit of the first two authors to the CTR. It became clear that there are open questions about the way in which the control algorithm models the boundary conditions, even for smooth walls, and that these questions make the physical interpretation of the results difficult. Considerable more work in that area seems to be needed before even relatively advanced large-eddy simulations, such as these, can be used to draw conclusions about the physics of wall-bounded turbulent flows. The numerical method is the same as in Nicoud et al. The modifications introduced in the original code are briefly described in section 2, but the original paper should be consulted for a full description of the algorithm. The results are presented in section 3 and summarized in section 4. The elementary properties of turbulence over rough walls which are used in the text have been taken from recent reviews such as Raupach et al. or Jimenez.
In order to study the effects of hypodynamia on the coagulating properties of the blood, vascular wall, and myocardium, chinchilla rabbits were kept for varying periods in special cages which restricted their movements. At the end of the experiment, blood samples were taken and the animals were sacrificed. Preparations were made from the myocardium venae cavae, and layers of the aorta. Two resultant interrelated and mutually conditioned syndromes were discovered: thrombohemorrhagic in the blood and hemorrago-thrombotic in the tissues.
Category 1, Problems 1 and 2, Category 2, Problem 2, and Category 3, Problem 2 are solved computationally using the Dispersion-Relation-Preserving (DRP) scheme. All these problems are governed by the linearized Euler equations. The resolution requirements of the DRP scheme for maintaining low numerical dispersion and dissipation as well as accurate wave speeds in solving the linearized Euler equations are now well understood. As long as 8 or more mesh points per wavelength is employed in the numerical computation, high quality results are assured. For the first three categories of benchmark problems, therefore, the real challenge is to develop high quality numerical boundary conditions. For Category 1, Problems 1 and 2, it is the curved wall boundary conditions. For Category 2, Problem 2, it is the internal radiation boundary conditions inside the duct. For Category 3, Problem 2, they are the inflow and outflow boundary conditions upstream and downstream of the blade row. These are the foci of the present investigation. Special nonhomogeneous radiation boundary conditions that generate the incoming disturbances and at the same time allow the outgoing reflected or scattered acoustic disturbances to leave the computation domain without significant reflection are developed. Numerical results based on these boundary conditions are provided.
A conservative flux difference splitting is presented for the hyperbolic systems of gasdynamics. The stable robust method is suitable for wide application in a variety of schemes, explicit or implicit, iterative or direct, for marching in either time or space. The splitting is modeled on the local quasi one dimensional characteristics system for multi-dimensional flow similar to Chakravarthy's nonconservative split coefficient matrix method; but, as the result of maintaining global conservation, the method is able to capture sharp shocks correctly. The embedded characteristics formulation is cast in a primitive variable the volumetric internal energy (rather than the pressure) that is effective for treating real as well as perfect gases. Finally the relationship of the splitting to characteristics boundary conditions is discussed and the associated conservative matrix formulation for a computed blown wall boundary condition is developed as an example. The theoretical development employs and extends the notion of Roe of constructing stable upwind difference formulae by sending split simple one sided flux difference pieces to appropriate mesh sites. The developments are also believed to have the potential for aiding in the analysis of both existing and new conservative difference schemes.
The inviscid Euler equations in complicated geometries are solved using a Cartesian grid. This requires solid wall boundary conditions in the irregular grid cells near the boundary. Since these cells may be orders of magnitude smaller than the regular grid cells, stability is a primary concern. An approach to this problem is presented and its use is illustrated.
The inviscid Euler equations in complicated geometries are solved using a Cartesian grid. This requires solid wall boundary conditions in the irregular grid cells near the boundary. Since these cells may be orders of magnitude smaller than the regular grid cells, stability is a primary concern. An approach to this problem is presented and its use is illustrated.
Fluctuating wall-pressure measurements are made in shock-wave/turbulent-boundary-layer interactions generated by sharp/unswept fins at angles of attack of 16, 18, 20, 22, 24, 26, and 28 degrees at Mach 5. The experiment was conducted under approximately adiabatic wall temperature conditions. The mean and rms pressure distributions can be collapsed in conical coordinates. The wall-pressure signal near separation is intermittent for all angles of attack (16-28 deg) and is qualitatively similar to that measured in unswept flows. However, the shock frequencies are higher - about 5 kHz compared to 0.5-1 kHz. Over the range of sweepbacks examined, from 25-55 deg, the spectral content of the fluctuating pressures does not change. Thus, the increase in separation-shock frequency from 1 to 5 kHz occurs at lower interaction sweepback and is not a continuous process with increasing sweepback. Power spectra at the position of maximum rms in the intermittent region for interactions in different incoming boundary layers have the same center frequency. The maximum rms in the intermittent region correlates with interaction sweepback, not with overall inviscid pressure rise.
The concept of a correctable-interference transonic wind tunnel combining a capability for wall-interference assessment with a limited capability for wall control is introduced. Recent progress toward achieving the interference assessment capability is described. The feasibility of using experimentally measured data directly as boundary values for the assessment in lieu of more generally formulated but less accurate wall boundary conditions is demonstrated for two-dimensional subsonic flows. Also, a procedure for transonic analysis of wall interference which leads to a rational definition of the wall-induced perturbation field within the context of nonlinear transonic flow computations is outlined in principle.
Stress and stability of thin-walled conical shells for axisymmetric loads
Two new approximate boundary conditions have been applied to the large eddy simulation of channel flow with and without transpiration. These new boundary conditions give more accurate results than those previously in use, and allow significant reduction of the required CPU time over simulations in which no-slip conditions are applied. Mean velocity profiles and turbulence intensities compare well both with experimental data and with the results of resolved simulations. The influence of the approximate boundary conditions remains confined near the point of application and does not affect the turbulence statistics in the core of the flow.
Work was carried out to formulate near-wall models for the equations governing the transport of the temperature-variance and its dissipation rate. With these equations properly modeled, a foundation is laid for their extension together with the heat-flux equations to compressible flows. This extension is carried out in a manner similar to that used to extend the incompressible near-wall Reynolds-stress models to compressible flows. The methodology used to accomplish the extension of the near-wall Reynolds-stress models is examined and the actual extension of the models for the Reynolds-stress equations and the near-wall dissipation-rate equation to compressible flows is given. Then the formulation of the near-wall models for the equations governing the transport of the temperature variance and its dissipation rate is discussed. Finally, a sample calculation of a flat plate compressible turbulent boundary-layer flow with adiabatic wall boundary condition and a free-stream Mach number of 2.5 using a two-equation near-wall closure is presented. The results show that the near-wall two-equation closure formulated for compressible flows is quite valid and the calculated properties are in good agreement with measurements. Furthermore, the near-wall behavior of the turbulence statistics and structure parameters is consistent with that found in incompressible flows.