Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stokes parameters”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Navier-Stokes Simulation of the Canard-Wing-Body Longitudinal Dynamic Stability Characteristics

Many modern aircraft are canard-configured for aircraft control and improved aerodynamic performance. Canards can often enhance aircraft cruise performance, maneuverability and agility. For close-coupled canard configurations, the aerodynamic interaction between the canard and wing significantly changes the flow characteristics of the wing. In unsteady flow, such changes in the flow structure and performance of wings can be quite pronounced. Accurate modeling of the unsteady aerodynamics is essential for potential CFD design and analysis of such configurations. A time-accurate numerical simulation is performed to study the unsteady aerodynamic interaction between a canard and wing with emphasis on the effects of the canard on the configuration's dynamic response characteristics. The thin-layer Reynolds-averaged Navier-Stokes Equations with various turbulence models are used in this study. Computations are made on a generic, analytically-defined, close-coupled canard-wing-body configuration which has been the subject of numerous previously published experimental studies during the 1970's to mid-80's. More recently, a series of steady-flow simulations has been performed and published by the author. In the current study, the configuration is given prescribed ramp and oscillatory motions in order to predict characteristics such as the damping-in-pitch and oscillatory longitudinal stability parameters. The current computations are made at high-subsonic and transonic Mach numbers, moderate angles-of- attack from -4 to 20 degrees, and at various pitch rates and reduced frequencies. Comparisons of pressures and integrated force quantities (e.g. lift, drag, pitching moment and selected dynamic parameters) are made with other published computational results and available experimental data. Results showing the unsteady effects of the canard on surface pressures, integrated forces, canard-wing vortex interaction and vortex breakdown will be presented.

Tu, Eugene L.↗
NASA NTRS1996

A Survey of Reynolds Number and Wing Geometry Effects on Lift Characteristics in the Low Speed Stall Region

This paper presents a survey of the effects of Reynolds number on the low- speed lift characteristics of wings encountering separated flows at their leading and side edges, with emphasis on the region near the stall. The influence of leading-edge profile and Reynolds number on the stall characteristics of two- dimensional airfoils are reviewed first to provide a basis for evaluating three- dimensional effects associated with various wing planforms. This is followed by examples of the effects of Reynolds number and geometry on the lift characteristics near the stall for a series of three-dimensional wings typical of those suitable for high-speed aircraft and missiles. Included are examples of the effects of wing geometry on the onset and spanwise progression of turbulent reseparation near the leading edge and illustrations of the degree to which simplified theoretical approaches can be useful in defining the influence of the various geometric parameters. Also illustrated is the manner in which the Reynolds number and wing geometry parameters influence whether the turbulent reseparation near the leading edge results in a sudden loss of lift, as in the two-dimensional case, or the formation of a leading-edge vortex with Rs increase in lift followed by a gentle stall as in the highly swept wing case. Particular emphasis is placed on the strong influence of 'induced camber' on the development of turbulent reseparation. R is believed that the examples selected for this report may be useful in evaluating viscous flow solutions by the new computational methods based on the Navier-Stokes equations as well as defining fruitful research areas for the high-Reynolds-number wind tunnels.

Polhamus, Edward C.↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗
NASA NTRS2005

CARS Temperature and Species Concentration Measurements in a Supersonic Combustor with Normal Injection

The dual-pump coherent anti-Stokes Raman spectroscopy (CARS) method was used to measure temperature and the absolute mole fractions of N2, O2 and H2 in a supersonic combustor. Experiments were conducted in the NASA Langley Direct-Connect Supersonic Combustion Test Facility. CARS measurements were performed at the facility nozzle exit and at three planes downstream of fuel injection. Processing the CARS measurements produced maps of the mean temperature, as well as quantitative N2 and O2 and qualitative H2 mean mole fraction fields at each plane. The CARS measurements were also used to compute correlations between fluctuations of the different simultaneously measured parameters. Comparisons were made between this 90 degree angle fuel injection case and a 30 degree fuel injection case previously presented at the 2004 Reno AIAA Meeting.

Tedder, S. A.↗
NASA NTRS1992

Automatic differentiation of advanced CFD codes for multidisciplinary design

Automated multidisciplinary design of aircraft and other flight vehicles requires the optimization of complex performance objectives with respect to a number of design parameters and constraints. The effect of these independent design variables on the system performance criteria can be quantified in terms of sensitivity derivatives which must be calculated and propagated by the individual discipline simulation codes. Typical advanced CFD analysis codes do not provide such derivatives as part of a flow solution; these derivatives are very expensive to obtain by divided (finite) differences from perturbed solutions. It is shown that sensitivity derivatives can be obtained accurately and efficiently using the ADIFOR source translator for automatic differentiation. In particular, it is demonstrated that the 3-D, thin-layer Navier-Stokes, multigrid flow solver called TLNS3D is amenable to automatic differentiation in the forward mode even with its implicit iterative solution algorithm and complex turbulence modeling. It is significant that by using computational differentiation, consistent discrete nongeometric sensitivity derivatives have been obtained from an aerodynamic 3-D CFD code in a relatively short time, e.g., O(man-week) not O(man-year).

Bischof, C.↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1993

A computational investigation of fuel mixing in a hypersonic scramjet

A parabolized, Navier-Stokes code, SHIP3D, is used to numerically investigate the mixing between air injection and hydrogen injection from a swept ramp injector configuration into either a mainstream low-enthalpy flow or a hypervelocity test flow. The mixing comparisons between air and hydrogen injection reveal the importance of matching injectant-to-mainstream mass flow ratios. In flows with the same injectant-to-mainstream dynamic pressure ratio, the mixing definition was altered for the air injection cases. Comparisons of the computed results indicate that the air injection cases overestimate the mixing performance associated with hydrogen injection simulation. A lifting length parameter, to account for the time a fluid particle transverses through the mixing region, is defined and used to establish a connection of injectant mixing in hypervelocity flows, based on nonreactive, low-enthalpy flows.

Fathauer, Brett W.↗
NASA NTRS1976

Atmospheric flow over two-dimensional bluff surface obstructions

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

Bitte, J.↗
NASA NTRS1999

Development of a Multiobjective Optimization Procedure for Sonic Boom Minimization

A design optimization procedure for improved sonic boom and aerodynamic performance of high speed aircraft is presented. The multiobjective optimization procedure simultaneously minimizes the sonic boom at a given distance from the aircraft and the drag-to-lift ratio (C(sub D)/C(sub L)) Of the aircraft. Upper and lower bounds are also imposed on the lift coefficient. The Kreisselmeier - Steinhauser function is used for the multiobjective optimization formulation. A discrete semi-analytical aerodynamic sensitivity analysis procedure coupled with an analytical grid sensitivity analysis technique is used for evaluating design sensitivities. The use of the semi-analytical sensitivity analysis techniques results in significant computational savings. The flow equations are solved using a three-dimensional parabolized Navier-Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. A nonlinear programming technique and an approximate analysis procedure are used for the optimization. The optimization procedure developed is applied to the design of two high speed configurations, namely, a doubly swept wing-body configuration and a delta wing-body configuration. For the two sweep case only, minimization of the first peak in the pressure signature is performed first by optimizing only the nose radius and length of the aircraft. Minimization of the second peak in the pressure signature is performed next by optimizing only the wing geometric parameters. Significant improvements are obtained in the sonic boom characteristics and the aerodynamic performance of the wing-body configurations.

Narayan, J. R.↗
DOE OSTI2023

Nanoparticle Diffusion in Miscible Polymer Nanocomposite Melts

X-ray photon correlation spectroscopy measurements were used to quantify the dynamics of bare and bimodal grafted silica nanoparticles mixed with PEO melts of different molecular weights. In dilute polymer nanocomposite (PNC) samples, we find diffusive NP behavior as described by the Stokes–Einstein relationship so long as the adsorbed PEO polymer layer is taken into account in determining both the effective NP size and its role on composite viscosity. The size of this bound layer was found to be approximately 2R g , where R g is the chain radius of gyration. We also expanded our system to investigate how the dynamics of grafted NPs differ from bare NPs with an adsorbed layer. We showed that the dynamics again can be determined by an effective NP radius at a scale smaller than the effective interparticle spacing; however, at larger length scales, the morphology and grafting parameters play a major role in the system dynamics. These results allow us to quantify NP ordering driven by polymer crystallization. It has previously been speculated that behavior is controlled by the relative ratio of time scale of crystal growth and the diffusive time scale of the NPs, a Peclet number. When the former time scale is longer, then the NPs are expected to be segregated into the interlamellar amorphous zones, while the NPs are too slow to be reorganized in the opposite case. We show here that this conjecture is quantitatively correct and the demarcation in behavior occurs for Pe = 1. Thus, we provide a way to estimate a critical spherulite growth rate for any semicrystalline PNC, at which a given NP can be ordered. Together, the results of this study permit us to tunably design PNCs through directed dispersion of NPs in a polymer matrix.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1987

Numerical study of ribbon-induced transition in Blasius flow

The early three-dimensional stages of transition in the Blasius boundary layer are studied by numerical solution of the Navier-Stokes equations. A finite-amplitude two-dimensional wave and low-amplitude three-dimensional random disturbances are introduced. Rapid amplification of the three-dimensional components is observed and leads to transition. For intermediate amplitudes of the two-dimensional wave the breakdown is of subharmonic type, and the dominant spanwise wavenumber increases with the amplitude. For high amplitudes the energy of the fundamental mode is comparable to the energy of the subharmonic mode, but never dominates it; the breakdown is of mixed type. Visualizations, energy histories, and spectra are presented. The sensitivity of the results to various physical and numerical parameters is studied. The agreement with experimental and theoretical results is discussed.

Spalart, Philippe R.↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Numerical simulation of boundary layers. Part 2: Ribbon-induced transition in Blasius flow

The early three-dimensional stages of transition in Blasius boundary layers are studied by numerical solution of the Navier-Stokes equations. A finite-amplitude two-dimensional wave and random low-amplitude three-dimensional disturbances are introduced. Rapid amplification of the three-dimensional components is observed and leads to transition. For intermediate amplitudes of the two-dimensional wave the breakdown is of subharmonic type, and the dominant spanwise wave number increases with the amplitude. For high amplitudes the energy of the fundamental mode is comparable to the energy of the subharmonic mode, but never dominates it; the breakdown is of mixed type. Visualizations, energy histories, and spectra are presented. The sensitivity of the results to various physical and numerical parameters is studied. Agreement with experimental and theoretical results is discussed.

Spalart, P.↗
DOE OSTI2020

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗
NASA NTRS1989

Estimations of ABL fluxes and other turbulence parameters from Doppler lidar data

Techniques for extraction boundary layer parameters from measurements of a short-pulse CO2 Doppler lidar are described. The measurements are those collected during the First International Satellites Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE). By continuously operating the lidar for about an hour, stable statistics of the radial velocities can be extracted. Assuming that the turbulence is horizontally homogeneous, the mean wind, its standard deviations, and the momentum fluxes were estimated. Spectral analysis of the radial velocities is also performed from which, by examining the amplitude of the power spectrum at the inertial range, the kinetic energy dissipation was deduced. Finally, using the statistical form of the Navier-Stokes equations, the surface heat flux is derived as the residual balance between the vertical gradient of the third moment of the vertical velocity and the kinetic energy dissipation. Combining many measurements would normally reduce the error provided that, it is unbiased and uncorrelated. The nature of some of the algorithms however, is such that, biased and correlated errors may be generated even though the raw measurements are not. Data processing procedures were developed that eliminate bias and minimize error correlation. Once bias and error correlations are accounted for, the large sample size is shown to reduce the errors substantially. The principal features of the derived turbulence statistics for two case studied are presented.

Gal-Chen, Tzvi↗
NASA NTRS1999

Design and Predictions for a High-Altitude (Low-Reynolds-Number) Aerodynamic Flight Experiment

A sailplane being developed at NASA Dryden Flight Research Center will support a high-altitude flight experiment. The experiment will measure the performance parameters of an airfoil at high altitudes (70,000 to 100,000 ft), low Reynolds numbers (200,000 to 700,000), and high subsonic Mach numbers (0.5 and 0.65). The airfoil section lift and drag are determined from pitot and static pressure measurements. The locations of the separation bubble, Tollmien-Schlichting boundary layer instability frequencies, and vortex shedding are measured from a hot-film strip. The details of the planned flight experiment are presented. Several predictions of the airfoil performance are also presented. Mark Drela from the Massachusetts Institute of Technology designed the APEX-16 airfoil, using the MSES code. Two-dimensional Navier-Stokes analyses were performed by Mahidhar Tatineni and Xiaolin Zhong from the University of California, Los Angeles, and by the authors at NASA Dryden.

Greer, Donald↗
NASA NTRS2005

Navier-Stokes Simulation of Airconditioning Facility of a Large Modem Computer Room

NASA recently assembled one of the world's fastest operational supercomputers to meet the agency's new high performance computing needs. This large-scale system, named Columbia, consists of 20 interconnected SGI Altix 512-processor systems, for a total of 10,240 Intel Itanium-2 processors. High-fidelity CFD simulations were performed for the NASA Advanced Supercomputing (NAS) computer room at Ames Research Center. The purpose of the simulations was to assess the adequacy of the existing air handling and conditioning system and make recommendations for changes in the design of the system if needed. The simulations were performed with NASA's OVERFLOW-2 CFD code which utilizes overset structured grids. A new set of boundary conditions were developed and added to the flow solver for modeling the roomls air-conditioning and proper cooling of the equipment. Boundary condition parameters for the flow solver are based on cooler CFM (flow rate) ratings and some reasonable assumptions of flow and heat transfer data for the floor and central processing units (CPU) . The geometry modeling from blue prints and grid generation were handled by the NASA Ames software package Chimera Grid Tools (CGT). This geometric model was developed as a CGT-scripted template, which can be easily modified to accommodate any changes in shape and size of the room, locations and dimensions of the CPU racks, disk racks, coolers, power distribution units, and mass-storage system. The compute nodes are grouped in pairs of racks with an aisle in the middle. High-speed connection cables connect the racks with overhead cable trays. The cool air from the cooling units is pumped into the computer room from a sub-floor through perforated floor tiles. The CPU cooling fans draw cool air from the floor tiles, which run along the outside length of each rack, and eject warm air into the center isle between the racks. This warm air is eventually drawn into the cooling units located near the walls of the room. One major concern is that the hot air ejected to the middle isle might recirculate back into the cool rack side and cause thermal short-cycling. The simulations analyzed and addressed the following important elements of the computer room: 1) High-temperature build-up in certain regions of the room; 2) Areas of low air circulation in the room; 3) Potential short-cycling of the computer rack cooling system; 4) Effectiveness of the perforated cooling floor tiles; 5) Effect of changes in various aspects of the cooling units. Detailed flow visualization is performed to show temperature distribution, air-flow streamlines and velocities in the computer room.

Source record↗
DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗