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.
The new code, OpenSG, is established as a novel multiscale constitutive modeling approach for realistic thin-walled composite structures. Its performance in terms of accuracy and efficiency is extensively verified against VABS/SwiftComp and direct numerical simulations of a wind blade.
High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.
The research objective of this work was to improve the engineering predictions of the turbulence characteristics of flows in curved narrow channels. Such channel flows are commonly encountered in nuclear research and test reactors, with one of them being the high-flux isotope reactor (HFIR). Research reactors bear high heat fluxes, and the proper computing of turbulence is paramount for safe and reliable reactor operation. The study builds on the results of a previous direct numerical simulation of turbulence to inform a well-known Reynolds-averaged Navier–Stokes shear-stress turbulence model and improves its accuracy in simulating parallel channel flows. A new formulation of the loss term in the dissipation conservation equation is suggested. Combined with high wall distance computational grids, the new implementation provides a fast-running flow solution, suitable for engineering purposes. Model generalization for parallel channel flows, in a broader range of frictional Reynolds numbers, is suggested by introducing a new form of the model constants.
Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.
The underlying bases and developments in two techniques of detailed turbulence modeling are described where the flow is treated in the Eulerian sense, and one technique where the Lagrangian motions of vortices are followed. First, a technique is described for solving the single-point statistically averaged conservation equations. The Reynolds stresses that appear in these equations are evaluated by solving supplemental differential equations which contain terms that are modeled. A sequence of increasingly complex, but also increasingly general, modeling equations is described and computations based on these equations are compared with experimental data. The hierarchy of models described terminates with equations for the individual components of the Reynolds stress tensor. The second Eulerian technique approach to turbulence modeling is the direct numerical simulation of turbulent fields. In this approach, all three dimensional eddies between a predetermined range of sizes are computed in time within a specified volume of flow. Present day computers require a tradeoff between the size of the volume that can be considered and the degree of resolution of the turbulent eddies. Techniques of modeling the smallest eddies are described that permit enlarging the volume, or Reynolds number, that can be considered.
Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.
Voyager 2's passage through the downstream corotation region of Ganymede found disturbances in the field and particle environment. Large fluctuations in the intensities and energy spectra of ions in the 0.1-to 4.0-MeV interval were also observed with the low energy charged particle experiment. All ion species were evidently affected, up through medium (C, N, O) nuclei. Electrons intensities, while disturbed, did not behave similarly to the ions. These effects are probably associated with the wake produced by the absorption of Jovian magnetospheric corotating plasma by Ganymede. A two-dimensional theoretical model for perturbed field in the wake region is proposed and it is shown by direct numerical simulation of exact particle trajectories that such perturbations are capable of modulating the charged particle energies to about the observed amount for the observed magnetic and plasma parameters. Analysis shows that the energization process is dependent on the pitch angles of the energetic ions and the process itself is nonadiabatic (the Larmor radius is approximately the size of the disturbance region).
The time evolution of an initially coherent, sinusoidal passive-scalar disturbance is considered when the wavelength q is less than the length scale of the surrounding isotropic turbulent flow. In 64 sup 3 direct numerical simulations a Gaussian prescription for the average scalar amplitude breaks down after a timescale associated with the wavenumber of the disturbance and there is a transition to a new characteristic decay. The Gaussian prescription is given by exp(-(1/2) q-squared w(t)), where a form for w(t), the Lagrangian mean square displacement of a single fluid particle, is proposed. After the transition the decay is given by exp(-t/tau), where tau is the new characteristic timescale. If q k(sub K), then 1/tau = 1/tau(sub D) + 1/tau(sub K), where k(sub K) is the Kolmogorov wavenumber, tau(sub D) is the diffusive timescale and tau(sub K) is the Kolmogorov timescale. An experiment originally proposed by de Gennesis considered in which the evolution of a coherent laser-induced pattern is read by a diffracting laser. The theory of this experiment involves the dispersion of particle pairs, but it is shown that in a certain limit it reduces to the single Fourier-mode problem and can be described in terms of single particle diffusion. The decay of a single mode after the transition in the simulation best describes the experiment.
The applicability of active control of transition by period suction-blowing is investigated via direct numerical simulations of the Navier-Stokes equations. The time-evolution of finite-amplitude disturbances in plane channel flow is compared in detail with and without control. The analysis indicates that for relatively small three-dimensional amplitudes, a two-dimensional control effectively reduces disturbance growth rates even for linearly unstable Reynolds numbers. After the flow goes through secondary instability, three-dimensional control seems necessary to stabilize the flow. An investigation of the temperature field suggests that passive temperature contamination is operative to reflect the vorticity dynamics during transition.
A three-dimensional, time-dependent, direct numerical simulation of low-Reynolds number turbulent flow in a mildly curved channel was performed, and the results examined to determine the mechanism by which curvature affects wall-bounded turbulent shear flows. A spectral numerical method with about one-million modes was employed, and no explicit subgrid scale model was used. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides is in the Reynolds shear stress. This is accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Goertler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by enhancing the difference between the wall shear stresses.
Direct numerical simulation using the full three dimensional, time dependent Navier-Stokes equations is used to investigate V/STOL jet induced interactions. The objective of this numerical simulation is to compute accurately the details of the flow field and to achieve a better understanding of the physics of the flow, including the role of initial turbulence in the jet, the influence of forward motion on hover aerodynamics, the collision zone and fountain characteristics. Preliminary results are presented.
A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.
The applicability of active control of transition by periodic suction-blowing is investigated via direct numerical simulations of the Navier-Stokes equations. The time-evolution of finite-amplitude disturbances in plane channel flow is compared in detail with and without control. The analysis indicates that, for relatively small three dimensional amplitudes, a two dimensional control effectively reduces disturbance growth rates even for linearly unstable Reynolds numbers. After the flow goes through secondary instability, three dimensional control seems necessary to stabilize the flow. An investigation of the temperature field suggests that passive temperature contamination is operative to reflect the flow dynamics during transition.
The mechanisms by which heat release affects the fluid dynamics in a turbulent reacting mixing layer are studied by direct numerical simulation. In agreement with previous laboratory experiments, the heat release is observed to lower the rate at which the mixing layer grows and to reduce the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer are shown to produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of fluid elements. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers are shown to result from vorticity generation by baroclinic torques. The density reductions also lower the generation rates of turbulent kinetic energy and the turbulent shear stresses, resulting in less turbulent mixing of fluid elements. Calculations of the energy in the various wave numbers show that the heat release has a stabilizing effect on the growth rates of individual modes. A linear stability analysis of a simlified model problem confirms this, showing that low density fluid in the mixing region will result in a shift of the frequency of the unstable modes to lower wave numbers (longer wavelengths). The growth rates of the unstable modes decrease, contributing to the slower growth of the mixing layer.
Direct numerical simulations of sink-flow boundary layers, with acceleration parameters K between 1.5 x 10 to the -6th and 3.0 x 10 to the -6th, are presented. The three-dimensional, time-dependent Navier-Stokes equations are solved numerically, using a spectral method, with about one million degrees of freedom. The flow is assumed to be statistically steady, and self-similar. A multiple-scale approximation and periodic conditions are applied to the fluctuations. The turbulence is studied using instantaneous and statistical results. Good agreement with the experiments of Jones and Launder (1972) is observed. The two effects of the favorable pressure gradient are to extend the logarithmic layer, and to alter the energy balance of the turbulence near the edge of the boundary layer. At low Reynolds number the logarithmic layer is shortened and slightly displaced, but wall-layer streaks are present even at the lowest values of R(theta) for which turbulence can be sustained. Large quiescent patches appear in the flow. Relaminarization occurs at K = 3.0 x 10 to the -6th, corresponding to a Reynolds number R(theta) of about 330.
The structures of the vorticity fields in several homogeneous irrotational straining flows and a homogeneous turbulent shear flow were examined using a database generated by direct numerical simulation of the unsteady Navier-Stokes equations. In all cases, strong evidence was found for the presence of coherent vortical structures. The initially isotropic vorticity fields were rapidly affected by imposed mean strain and the rotational component of mean shear and developed accordingly. In the homogeneous turbulent shear-flow cases, the roll-up of mean vorticity into characteristic hairpin vortices was clearly observed, supporting the view that hairpin vortices are an important vortical structure in all turbulent shear flows; the absence of mean shear in the homogeneous irrotational straining flows precludes the presence of hairpin-like vortices.
A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of 3300, based on the mean centerline velocity and channel half-width, with about 4 million grid points. All essential turbulence scales are resolved on the computational grid and no subgrid model is used. A large number of turbulence statistics are computed and compared with the existing experimental data at comparable Reynolds numbers. Agreements as well as discrepancies are discussed in detail. Particular attention is given to the behavior of turbulence correlations near the wall. A number of statistical correlations which are complementary to the existing experimental data are reported for the first time.
A mathematical technique for analyzing local energy-transfer rates among wave-vector triads is developed and applied to the data generated by direct numerical simulations of waves with various types of initial conditions. Starting the simulation with a primary two-dimensional wave and random noise produced structures very similar to those which evolved from three-dimensional center modes. The local transfer rates determined in this case help to explain the eventual deformation of the primary two-dimensional wave which was observed in wind-tunnel experiments. Including weak streamwise vortices in the initial flow field results in large amounts of energy being transferred to the K-type modes early in the simulation. The later development of the waves (and hence the energy-transfer rates) is similar to the previous case.