Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence optimization”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Use of multiple discrete wall jets for delaying boundary layer separation

The effectiveness of a spanwise array of small discrete blowing nozzles in preventing separation of a turbulent boundary layer was investigated experimentally. The spacing, axial location, and momentum flux of the nozzles were varied in a systematic way, and overall performance was measured for each combination. Extensive mean velocity profiles were measured for one selected combination. Overall diffusion achieved before separation was correlated successfully with a momentum flux excess parameter, and in terms of this parameter discrete nozzles, when advantageously placed, were found to perform somewhat better than an optimally placed two-dimensional jet slot.

Mclean, J. D.↗
NASA NTRS1995

Adaptive unstructured triangular mesh generation and flow solvers for the Navier-Stokes equations at high Reynolds number

A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.

Ashford, Gregory A.↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2024

Active Nematic Ratchet in Asymmetric Obstacle Arrays

We numerically investigate the effect of an asymmetric periodic obstacle array in a two-dimensional active nematic. We find that activity in conjunction with the asymmetry leads to a ratchet effect or unidirectional flow of the fluid along the asymmetry direction. The directional flow is still present even in the active turbulent phase when the gap between obstacles is sufficiently small. We demonstrate that the dynamics of the topological defects transition from flow-mirroring to smectic-like as the gap between obstacles is made smaller, and explain this transition in terms of the pinning of negative winding number defects between obstacles. This also leads to a non-monotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Program of Research in Laminar Flow Control in the JIAFS at NASA Langley Research Center

At high Reynolds numbers, the performance of a low drag suction LFC airplane is essentially controlled by the induced drag and the turbulent friction drag of the nonlaminarized area. The question then arises as to how the airplane cruise lift to drag ratio (L/D) sub cruise) cruise increases with increasing extent of laminar flow 0 sub lam./0 sub total (0 = airplane wetted area). In particular, the question arises as to the airplane performance in the optimum case with all laminar flow over the airplane wetted area. Design approaches of all laminar flow LFC airplanes which optimize the airplane range . R = eta sub ov times (L/D) times H times lambda n(W sub 0/W sub E) are considered.

Pfenninger, W.↗
NASA NTRS2009

Quantifying Parameter Sensitivity, Interaction and Transferability in Hydrologically Enhanced Versions of Noah-LSM over Transition Zones

We use sensitivity analysis to identify the parameters that are most responsible for shaping land surface model (LSM) simulations and to understand the complex interactions in three versions of the Noah LSM: the standard version (STD), a version enhanced with a simple groundwater module (GW), and version augmented by a dynamic phenology module (DV). We use warm season, high-frequency, near-surface states and turbulent fluxes collected over nine sites in the US Southern Great Plains. We quantify changes in the pattern of sensitive parameters, the amount and nature of the interaction between parameters, and the covariance structure of the distribution of behavioral parameter sets. Using Sobol s total and first-order sensitivity indexes, we show that very few parameters directly control the variance of the model output. Significant parameter interaction occurs so that not only the optimal parameter values differ between models, but the relationships between parameters change. GW decreases parameter interaction and appears to improve model realism, especially at wetter sites. DV increases parameter interaction and decreases identifiability, implying it is overparameterized and/or underconstrained. A case study at a wet site shows GW has two functional modes: one that mimics STD and a second in which GW improves model function by decoupling direct evaporation and baseflow. Unsupervised classification of the posterior distributions of behavioral parameter sets cannot group similar sites based solely on soil or vegetation type, helping to explain why transferability between sites and models is not straightforward. This evidence suggests a priori assignment of parameters should also consider climatic differences.

Rosero, Enrique↗
NASA NTRS1980

Application of two design methods for active flutter suppression and wind-tunnel test results

The synthesis, implementation, and wind tunnel test of two flutter suppression control laws for an aeroelastic model equipped with a trailing edge control surface are presented. One control law is based on the aerodynamic energy method, and the other is based on results of optimal control theory. Analytical methods used to design the control laws and evaluate their performance are described. At Mach 0.6, 0.8, and 0.9, increases in flutter dynamic pressure were obtained but the full 44 percent increase was not achieved. However at Mach 0.95, the 44 percent increase was achieved with both control laws. Experimental results indicate that the performance of the systems is not so effective as that predicted by analysis, and that wind tunnel turbulence plays an important role in both control law synthesis and demonstration of system performance.

Newsom, J. R.↗
NASA NTRS1984

Convergent-tapered annular seals - Analysis and testing for rotordynamic coefficients

A combined analytical-computational method has been developed to calculate the pressure field and dynamic coefficients for tapered high-pressure annular seals. Completely developed turbulent flow is assumed in both the circumferential and axial directions, according to Hirs' (1973) bulk-flow turbulent-lubrication equations. In a numerical experiment with the method, a short bearing approximation is used to derive an analytical expression for the first-order (dynamic) pressure gradient. This expression is integrated numerically to define dynamic coefficients of the seal. Numerical results of the integration are compared with previous results for straight and tapered seals. It is found that direct stiffness and leakage coefficients increase in the present seal, while the remaining coefficients decrease. An optimal taper angle is shown to exist with respect to: (1) the direct stiffness and (2) the ratio of direct stiffness to leakage. On the basis of the optimal taper angle calculations, stiffness increases on the order of 40 to 50 percent are predicted. Experimental results for three taper angles are presented, and are found to be in good agreement with the theoretical predictions.

Childs, D. W.↗
NASA NTRS2000

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗
NASA NTRSDate not supplied

Physics and Modeling of Quasi-Steady and Forced Shock Trains

The dual-mode scramjet is a propulsive concept for contemporary hypersonic vehicles promising flight cruise speeds upwards of Mach 5. Despite decades of development and a recent resurgence in interest, such engines are not routinely deployed for commercial or military purposes. Instead, they have been confined to the realm of ground-test research and limited one-off flight test missions. One of the key stumbling blocks on the critical path to deployment of a dual-mode scramjet is the optimal design and stable operation of the isolator, a typically constant area duct that separates the scramjet inlet from the combustor. During ramjet mode at Mach numbers from 3 to 5 and during the transition to scramjet mode at higher Mach numbers, the isolator contains the pseudo-shock. Necessary to meet the back-pressure in the combustor, the pseudo-shock is typically divided into a shock dominated shock train region followed by a turbulence driven mixing region. Fluctuations in combustor pressure due to vehicle maneuvering, throttling, or combustion instabilities drive the shock train to shift streamwise within the isolator. In response to sharp pressure fluctuations, the shock train is unable to stabilize within the isolator, and the scramjet unstarts by ejecting the shock train out the inlet. This work addresses three key gaps in understanding of the isolator pseudo-shock. The phenomenon is studied in three experimental facilities to investigate the quasi-steady and unsteady fluid dynamics of the shock train in fueled and un-fueled test conditions. The strong link between boundary layer morphology and shock train structure is studied in detail and used as a tool to evaluate the state-of-the-art definitions and models of the pseudo-shock. Identifying the limitations of these models, a new process driven flux-conserved pseudo-shock model is proposed. This provides a computationally efficient model grounded in the fluid dynamic mechanisms of the pseudo-shock with better predictive capability than the current state-of-the-art across all facilities tested. The unsteady behavior of the shock train is further investigated to identify and understand two distinct forced dynamic modes: mechanically induced, low-rate forcing and high-rate combustion driven unsteadiness that sheds upstream propagating shocks. The proposed process driven model is found to be effective in the the limiting case of high-rate combustion driven unsteadiness by evaluating the model in a shock stationary reference frame. Further analysis identifies the approximately rigid chain motion of the shock train structure to mechanically induced forcing. The time delay between the onset of forced dynamics is found to saturate to a constant value and maximum shock train speed along the isolator with respect to the perturbation rate of change. The result of this work provides dual-mode scramjet designers with improved tools and physical understanding to develop more efficient, robust, and capable propulsion systems for hypersonic vehicles.

JAXA↗
NASA NTRS2001

Predicting Turbulent Convective Heat Transfer in Fully Developed Duct Flows

The performance of an explicit algebraic stress model (EASM) is assessed in predicting the turbulent flow and forced heat transfer in both straight and wavy ducts, with rectangular, trapezoidal and triangular cross-sections, under fully developed conditions. A comparison of secondary flow patterns. including velocity vectors and velocity and temperature contours, are shown in order to study the effect of waviness on flow dynamics, and comparisons between the hydraulic parameters. Fanning friction factor and Nusselt number, are also presented. In all cases. isothermal conditions are imposed on the duct walls, and the turbulent heat fluxes are modeled using gradient-diffusion type models. The formulation is valid for Reynolds numbers up to 10(exp 5) and this minimizes the need for wall functions that have been used with mixed success in previous studies of complex duct flows. In addition, the present formulation imposes minimal demand on the number of grid points without any convergence or stability problems. Criteria in terms of heat transfer and friction factor needed to choose the optimal wavy duct cross-section for industrial applications among the ones considered are discussed.

Rokni, Masoud↗
NASA NTRS2024

An Extensible Perturbed Parameter Ensemble for the Community Atmosphere Model Version 6

This paper documents the methodology and preliminary results from a Perturbed Parameter Ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2-CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations are analyzed and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks. The use of various emulators is explored in the multi- dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that using PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2-CAM6 PPE is publicly available, and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g. coupling to the land surface).

Machine learning↗
DOE OSTI2024

An extensible perturbed parameter ensemble for the Community Atmosphere Model version 6

This paper documents the methodology and preliminary results from a perturbed parameter ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2–CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks are analyzed. The use of various emulators is explored in the multi-dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties, are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2–CAM6 PPE is publicly available and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g., coupling to the land surface).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

A computational study of thrust augmenting ejectors based on a viscous-inviscid approach

A viscous-inviscid interaction technique is advocated as both an efficient and accurate means of predicting the performance of two-dimensional thrust augmenting ejectors. The flow field is subdivided into a viscous region that contains the turbulent jet and an inviscid region that contains the ambient fluid drawn into the device. The inviscid region is computed with a higher-order panel method, while an integral method is used for the description of the viscous part. The strong viscous-inviscid interaction present within the ejector is simulated in an iterative process where the two regions influence each other en route to a converged solution. The model is applied to a variety of parametric and optimization studies involving ejectors having either one or two primary jets. The effects of nozzle placement, inlet and diffuser shape, free stream speed, and ejector length are investigated. The inlet shape for single jet ejectors is optimized for various free stream speeds and Reynolds numbers. Optimal nozzle tilt and location are identified for various dual-ejector configurations.

Lund, Thomas S.↗
NASA NTRS1989

Use of piezoelectric foil for flow diagnostics

A laboratory investigation was conducted to characterize two piezoelectric-film sensor configurations, a rigidly mounted sensor and a sensor mounted over an air cavity. The sensors are evaluated for sensitivity and frequency response, and methods to optimize data are presented. The cavity-mounted sensor exhibited a superior frequency response and was more sensitive to normal pressure fluctuations and less sensitive to vibrations through the structure. Both configurations were sensitive to large-scale structural vibrations. Flight-test data are shown for cavity-mounted sensors, illustrating practical aspects to consider when designing sensors for application in such harsh environments. The relation of the data to skin friction and maximum shear stress, transition detection, and turbulent viscous layers is derived through analysis of the flight data.

Carraway, Debra L.↗
NASA NTRS1995

Sailplane Glide Performance and Control Using Fixed and Articulating Winglets

An experimental study was conducted to investigate the effects of controllable articulating winglets on glide performance and yawing moments of high performance sailplanes. Testing was conducted in the Texas A&M University 7 x 10 foot Low Speed Wind Tunnel using a full-scale model of the outboard 5.6 feet of a 15 meter class high performance sailplane wing. Different wing tip configurations could be easily mounted to the wing model. A winglet was designed in which the cant and toe angles as well as a rudder on the winglet could be adjusted to a range of positions. Cant angles used in the investigation consisted of 5, 25, and 40 degrees measured from the vertical axis. Toe-out angles ranged from 0 to 22.5 degrees. A rudder on the winglet was used to study the effects of changing the camber of the winglet airfoil on wing performance and wing yawing moments. Rudder deflections consisted of-10, 0, and 10 degrees. Test results for a fixed geometry winglet and a standard wing tip are presented to show the general behavior of winglets on sailplane wings, and the effects of boundary-layer turbulators on the winglets are also presented. By tripping the laminar boundary-layer to turbulent before laminar separation occurs, the wing performance was increased at low Reynolds numbers. The effects on the lift and drag, yawing moment, pitching moment, and wing root bending moment of the model are presented. Oil flows were used on the wing model with the fixed geometry winglet and the standard wing tip to visualize flow directions and areas of boundary layer transition. A cant angle of 25 degrees and a toe-out angle of 2.5 degrees provided an optimal increase in wing performance for the cant and toe angles tested. Maximum performance was obtained when the winglet rudder remained in the neutral position of zero degrees. By varying the cant, toe, and rudder angles from their optimized positions, wing performance decreases. Although the winglet rudder proved to be more effective in increasing the yawing moment compared to varying the cant and toe angles, the amount of increased yawing moment was insignificant when compared to that produced by the vertical tail. A rudder on the winglet was determined to be ineffective for providing additional yaw control.

Colling, James David↗
NASA NTRS2003

Formation Tetrahedron Design for Phase I of the Magnetospheric Multiscale Mission

The Magnetospheric Multiscale Mission (MMS) is a NASA mission intended to make fundamental advancements in our understanding of the Earth's Magnetosphere. There are three processes that MMS will study including magnetic reconnection, charged particle acceleration, and turbulence. There are four phases in the nominal mission and this work addresses some of the outstanding issues in phase I. The nominal phase I orbit is 1.2 x 12 R(sub e) highly elliptic orbit with four spacecraft nominally forming a regular tetrahedron. In this paper we investigate the relative dynamics of the four MMS spacecraft about an assumed reference orbit. There are several tetrahedron dimensions required in Phase I of the mission and in this work we design optimal tetrahedrons for the 10 km baseline. The performance metric used in the optimization process is directly related to the science return, and is based on an extension of previous work performed by Glassmeier. The optimizer we use is a commercially available Sequential Quadratic Programming (SQP) routine. Multiple optimal solutions are found, and we characterize how the performance of the formation varies between different regions of the reference orbit.

Hughes, Steven P.↗
NASA NTRS2018

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗