Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Utilizing Direct Numerical Simulations of Transition and Turbulence in Design Optimization

Design optimization methods that use the Reynolds-averaged Navier-Stokes equations with the associated turbulence and transition models, or other model-based forms of the governing equations, may result in aerodynamic designs with actual performance levels that are noticeably different from the expected values because of the complexity of modeling turbulence/transition accurately in certain flows. Flow phenomena such as wake-blade interaction and trailing edge vortex shedding in turbines and compressors (examples of such flows) may require a computational approach that is free of transition/turbulence models, such as direct numerical simulations (DNS), for the underlying physics to be computed accurately. Here we explore the possibility of utilizing DNS data in designing a turbine blade section. The ultimate objective is to substantially reduce differences between predicted performance metrics and those obtained in reality. The redesign of a typical low-pressure turbine blade section with the goal of reducing total pressure loss in the row is provided as an example. The basic ideas presented here are of course just as applicable elsewhere in aerodynamic shape optimization as long as the computational costs are not excessive.

Direct numerical simulations of transition and tur↗
NASA NTRS1980

Numerical prediction of wall temperatures for near-critical para-hydrogen in turbulent upflow inside vertical tubes

Presented herein is a method of including density fluctuations in the equations of turbulent transport. Results of a numerical analysis indicate that the method may be used to predict heat transfer for the case of near-critical para-hydrogen in turbulent upflow inside vertical tubes. Wall temperatures, heat transfer coefficients, and velocities obtained by coupling the equations of turbulent momentum and heat transfer with a perturbed equation of state show good agreement with experiment for inlet reduced pressures of 1.28-5.83.

Bellmore, C. P.↗
DOE OSTI2024

One-to-one aeroservoelastic validation of operational loads and performance of a 2.8 MW wind turbine model in OpenFAST

Abstract. This article presents a validation study of the popular aeroservoelastic code suite OpenFAST leveraging weeks of measurements obtained during normal operation of a 2.8 MW land-based wind turbine. Measured wind conditions were used to generate one-to-one turbulent flow fields (i.e., comparing simulation to measurement in 10 min increments, or bins) through unconstrained and constrained assimilation methods using the kinematic turbulence generators TurbSim and PyConTurb. A total of 253 bins of 10 min of normal turbine operation were selected for analysis, and a statistical comparison in terms of performance and loads is presented. We show that successful validation of the model was not strongly dependent on the type of inflow assimilation method used for mean quantities of interest, which had median modeling errors per wind-speed interval generally within 5 %–10 % of the measurement. The type of inflow assimilation method did have a larger effect on the fatigue predictions for blade-root flapwise and tower-base fore–aft quantities, which surprisingly saw larger errors from the assumed higher-fidelity assimilation methods. Avenues for further work are discussed and include possible improvements to the aerodynamic, structural, and controller modeling that may offer insight on the origin of the up to ∼ 40 % median overprediction of fatigue for these quantities.

17 WIND ENERGY↗
NASA NTRS1970

Calculation of the planar supercritical flow over a NASA supercritical profile

An unsteady finite difference procedure was used to calculate the steady inviscid flow over an 11% thickness ratio NASA supercritical profile of LWP 505 at M sub infinity = 0.80 and alpha = 0. An attempt is made to include the viscous effects using a modified form of Head's entrainment method to calculate the turbulent boundary layer. The attempt to predict the viscous effects using the compressible form of Head's integral method with a modified auxiliary equation for the form factors was unsatisfactory. Though a reasonably separation bubble was established on the lower surface, a grossly exaggerated displacement effect resulted downstream of the shock on the upper surface. There clearly is substantial further effort required to evolve a satisfactory boundary layer procedure, which must then be coupled in a still unproven manner with the inviscid procedure.

Yoshihara, H.↗
NASA NTRS1989

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.↗
NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testingon the CobraMRV Mars Entry Vehicle Concept

The CFD solver OVERFLOW was used to simulate the CobraMRV undergoing supersonic retropropulsion (SRP) in the Langley Unitary Plan Wind Tunnel as part of a pre-test study. Aerodynamics of the CobraMRV are summarized and the sensitivity of vehicle loads to CFD parameters at a subset of operating conditions are investigated. Specific numerical methods, boundary conditions, and turbulence modeling options have been selected after considering the complex phenomena in SRP flows. Proper shock-capturing methods and dynamic grid adaption are necessary to correctly capture vehicle loads and dynamics. Present data indicates that at particular conditions, upstream asymmetries originating from the tunnel inflow plane contribute to measurable asymmetries on the CobraMRV heatshield.

OVERFLOW↗
NASA NTRS1979

Venera-11 and Venera 12: Preliminary estimates for the wind speed and turbulence in the atmosphere of Venus

The methods and results of measurements for wind speed and atmospheric turbulence in the clouds of Venus are described, and compared with earlier results. The distribution of wind speed obtained from the data of Venera 12 is in good conformity with the data of the preceding Venera and Pioneer probes, indicating the existence of a constant and powerful zonal movement of the troposphere.

Kerzhanovich, V. V.↗
NASA NTRS2009

Turbulent CO2 Flux Measurements by Lidar: Length Scales, Results and Comparison with In-Situ Sensors

The vertical CO2 flux in the atmospheric boundary layer (ABL) is investigated with a Doppler differential absorption lidar (DIAL). The instrument was operated next to the WLEF instrumented tall tower in Park Falls, Wisconsin during three days and nights in June 2007. Profiles of turbulent CO2 mixing ratio and vertical velocity fluctuations are measured by in-situ sensors and Doppler DIAL. Time and space scales of turbulence are precisely defined in the ABL. The eddy-covariance method is applied to calculate turbulent CO2 flux both by lidar and in-situ sensors. We show preliminary mean lidar CO2 flux measurements in the ABL with a time and space resolution of 6 h and 1500 m respectively. The flux instrumental errors decrease linearly with the standard deviation of the CO2 data, as expected. Although turbulent fluctuations of CO2 are negligible with respect to the mean (0.1 %), we show that the eddy-covariance method can provide 2-h, 150-m range resolved CO2 flux estimates as long as the CO2 mixing ratio instrumental error is no greater than 10 ppm and the vertical velocity error is lower than the natural fluctuations over a time resolution of 10 s.

Gilbert, Fabien↗
NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗
NASA NTRS1974

A new technique for inducing a turbulent boundary layer in a supersonic flow

The feasibility of using slot injection to establish a turbulent boundary layer corresponding to a known Reynolds number was investigated. The basic concept consists of injecting secondary air through a slot. The air is injected at the same local static pressure at a selected stagnation temperature to simulate a given up-stream heat conduction. The mass of the injected air is controlled to match a desired Reynolds number. Secondary air was injected through a supersonic nozzle over a flat plate model. Tests were conducted with different injection air conditions. Heat transfer rates were measured behind the slot and results were compared with theoretical estimates. Static pressure measurements were made over the surface in the streamwise direction. Velocity and Mach number profiles were determined downstream of the slot. Laminar and turbulent profiles were calculated to compare with the results of measurements. Experimental results indicate that this is an effective method for establishing a turbulent boundary layer having a prescribed value of R sub e Theta.

Miyazawa, M.↗
NASA NTRS1985

Large-eddy substitution via vortex cancellation for wall turbulence control

A system of co-rotating longitudinal vortices was used to introduce streamline (as opposed to wall) curvature into a turbulent wall flow. Two methods of vortex cancellation, unwinding and self-annihilation, were tested as a means of removing the vortices once they had processed most of the incoming turbulent boundary layer. Vortex unwinding, which uses vorticity of the opposite sign, was shown to be a viable method for cancelling the co-rotating vortices. Vortex self-annihilation, caused by interference effects resulting from a close initial spanwise vortex spacing, eliminated the vortices within 60 delta downstream. In each case, reductions in boundary layer entrainment were found once the vortices were cancelled.

Mcginley, C. B.↗
NASA NTRS1981

Scale effects on turbulent boundary layer development and flow separation around V/STOL inlets at high incidence

Numerical methods for calculating laminar and turbulent boundary layers development around V/STOL engine inlets at high incidence angles, along with the procedures for predicting flow separation, are presented. Results of scale-effects, which are obtained by a numerical scaling procedure on the boundary layer characteristics and incidence angles at onset of separation are discussed. The interesting 'cross-over' phenomena, where the full-scale models actually exhibit earlier separation than the scaled-models, is illustrated for a typical V/STOL inlet at a certain operating condition. Some of the numerical results are compared with the existing wind-tunnel test data for a 1/6 scale inlet model to demonstrate the validity of the numerical approach.

Chou, D. C.↗
NASA NTRS1988

A study of turbulent flow computations in an angled duct with a step

This paper presents a method for effective computation of turbulent flows in an angled duct which has a step at the angled section by using a grid generation technique. An effect of flow curvature was analyzed on the separating and reattaching flows. The computations are based on the Reynolds averaged transport equations which are transformed for use on a generalized curvilinear coordinate system. A special care was taken on the treatment of pressure corrections. Wall function boundary conditions for k-epsilon model were developed.

Maruszewski, J. P.↗
NASA NTRS1992

The present state and future directions of PDF methods

The objectives of the workshop are presented in viewgraph format, as is this entire article. The objectives are to discuss the present status and the future direction of various levels of engineering turbulence modeling related to Computational Fluid Dynamics (CFD) computations for propulsion; to assure that combustion is an essential part of propulsion; and to discuss Probability Density Function (PDF) methods for turbulent combustion. Essential to the integration of turbulent combustion models is the development of turbulent model, chemical kinetics, and numerical method. Some turbulent combustion models typically used in industry are the k-epsilon turbulent model, the equilibrium/mixing limited combustion, and the finite volume codes.

Pope, S. B.↗
NASA NTRS1976

Turbulent boundary layers over nonstationary plane boundaries

Methods of predicting integral parameters and skin-friction coefficients of turbulent boundary layers developing over moving-ground-planes are evaluated using test information from three different wind tunnel facilities at the NASA Langley Research Center. These data include test information from the VSTOL tunnel which is presented for the first time. The three methods evaluated were: (1) relative integral parameter method, (2) relative power law method, and (3) modified law of the wall method. Methods (1) and (2) can be used to predict moving-ground-plane shape factors with an expected accuracy of + or - 10%. They may also be used to predict moving-ground-plane displacement and momentum thicknesses with lower expected accuracy. This decrease in accuracy can be traced to the failure of approximations upon which these methods are based to prove universal when compared with VSTOL tunnel test results.

Roper, A. T.↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
NASA NTRS1992

Ship viscous flow: A report on the 1990 SSPA-IIHR Workshop

To assess the state of the art in ship viscous flow computation a Workshop was organized in 1990 by three organizations: SSPA Maritime Consulting AB, Chalmers University of Technology, and the Iowa Institute of Hydraulic Research. Two test cases were specified by the organizers and sent out to all interested research groups, which were asked to submit results in a prescribed format. In September 1990 a meeting was held at Chalmers University of Technology. All results had then been collected and presented in a common format, and the theories based on responses to a questionnaire sent out earlier. During the meeting, each research group was first given the opportunity to briefly introduce their method and results. Thereafter, a considerable time was spent on general discussions on the performance of the different methods considering the differences in the underlying theories. Specific items that were addressed were grid generation, governing equations, boundary conditions, turbulence modelling, and numerical methods. Practical aspects of the results, for instance from the point of view of propeller design, were also discussed. The Workshop Proceedings contain a description of the participating methods, and the results of both test cases. In the present paper, a summary of the Workshop and its results is presented.

Patel, Virendra C.↗