Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient flow”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Analysis of separated boundary-layer flows

A method for predicting the strong interaction between the viscous and inviscid flows which occurs in flow separation is reviewed. An inverse boundary-layer procedure approximately accounts for normal pressure gradients that may be important in strongly interacting flows. Transformed boundary-layer equations are written in which the pressure gradient is set equal to the inviscid pressure gradient. As the boundary-layer edge is approached and the viscous shear and heat conduction terms vanish, the viscous flow solution is required to asymptotically approach the inviscid solution over the generalized displacement body. Attention is then focused on viscous-inviscid interacting flows with a first-order viscous formulation and constant pressure across the boundary layer. Results obtained with this procedure are presented for: (1) transitional separation bubbles near an airfoil leading edge, (2) subsonic boattail separated turbulent flow, and (3) transonic turbulent shock wave boundary-layer interaction on an axisymmetric bump configuration.

Carter, J. E.↗
NASA NTRS1998

The Role of Marangoni Convection for the FZ-Growth of Silicon

Fluctuations of the electrical resistivity due to inhomogeneous dopant distribution are still a serious problem for the industrial processing yield of doped silicon crystals. In the case of silicon floating-zone growth, the main sources of these inhomogeneities are time- dependent flows in the liquid phase during the growth process. Excluding radio frequency (RF) induced convection, buoyancy and thermocapillary (Marangoni) convection are the two natural reasons for fluid flow. Both originate from temperature/concentration gradients in the melt, buoyancy convection through thermal/concentrational volume expansion, and thermocapillary convection through the temperature/concentration dependence of the surface tension. To improve the properties of grown crystals, knowledge of the strength, the characteristic, and the relation of these two flow mechanisms is essential. By the use of microgravity, the effect and the strength of buoyancy (gravity dependent) and thermocapillary (gravity independent) convection can be separated and clarified. Applying magnetic fields, both convective modes can be influenced: fluid flow can either be damped (static magnetic fields) or overlaid by a regular flow regime (rotating magnetic fields). Two complementary approaches have been pursued: Silicon full zones (experiments on the German sounding rockets TEXUS 7, 12, 22, 29, and 36) with the maximum temperature at half of the zone height and silicon half zones (experiments on the Japanese sounding rockets TR-IA4 and 6) with the maximum temperature at the top of the melt. With the full zone arrangement, the intensity and the frequency of the dopant striations could be determined and the critical Marangoni number could be identified. The half zone configuration is suited to classify the flow pattern and to measure the amplitude and the frequency of temperature fluctuations in the melt by inserting thermocouples or temperature sensors into the melt. All experiments have been carried out in monoellipsoid mirror furnaces. Typical zone geometries are approx. 8 to 14 mm in diameter and height. The crystals grown under microgravity are compared to crystals grown in static axial magnetic fields (B<5 tesla) and in transversal rotating magnetic fields (B<7.5 mT / f=50 Hz). Experimental results are completed by 3D numerical simulations: the obtained temperature and concentration distribution in the melt confirm the damping effect of rotating magnetic fields and explain the change in the radial segregation under static magnetic fields.

Dold, P.↗
NASA NTRS2003

Optimal Disturbances in Boundary Layers Subject to Streamwise Pressure Gradient

An analysis of the non-modal growth of perturbations in a boundary layer in the presence of a streamwise pressure gradient is presented. The analysis is based on PSE equations for an incompressible fluid. Examples with Falkner- Skan profiles indicate that a favorable pressure gradient decreases the non-modal growth while an unfavorable pressure gradient leads to an increase of the amplification. It is suggested that the transient growth mechanism be utilized to choose optimal parameters of tripping elements on a low-pressure turbine (LPT) airfoil. As an example, a boundary-layer flow with a streamwise pressure gradient corresponding to the pressure distribution over a LPT airfoil is considered. It is shown that there is an optimal spacing of the tripping elements and that the transient growth effect depends on the starting point. The amplification is found to be small at the LPT s very low Reynolds numbers, but there is a possibility to enhance the transient energy growth by means of wall cooling.

Ashpis, David E.↗
NASA NTRS2004

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗
NASA NTRS1985

A class of airfoils having finite trailing edge pressure gradients

Some new developments relevant to the design of single-element airfoils using potential flow methods are presented. In particular, the ramifications of the unbounded trailing edge pressure gradients generally present in the potential flow solution for the flow over an airfoil are examined, and the conditions necessary to obtain a class of airfoils having finite trailing edge pressure gradients developed. The incorporation of these conditions into the inverse method of Eppler for the design of low-speed airfoils is discussed, and designs generated using the modified scheme are presented for consideration. A detailed viscous analysis of one of these airfoils demonstrates a significant reduction in the strong inviscid-viscid interactions generally present near the trailing edge. These reductions offer the possibility of improved airfoil performance, as well as the possibility of improved accuracy in the methods of airfoil design and analysis.

Ormsbee, A. I.↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Solutal convection during directional solidification

During directional solidification of a binary alloy at constant velocity, buoyancy-driven fluid flow may occur due to the solute gradients generated by the solidification process. Numerical calculations of the solute and fluid flow fields in the melt have been carried out using finite differences in a two-dimensional, time-dependent model that assumes a planar crystal-melt interface and allows time-dependent gravitational accelerations. The container walls are rigid and perfectly insulating for solute. For constant vertical gravitational accelerations, as the solutal Rayleigh number is varied, multiple steady states and time-dependent states may occur. The bifurcation from the quiescent state may be subcritical or transcritical, depending on the aspect ratio of the container. Calculations have also been performed for a gravitational acceleration that is assumed to be uniform in magnitude with its direction rotating uniformly. Numerical results have been obtained for a Schmidt number of 10 and a gravitational acceleration of 0.0001 G. The maximum variation in the solute concentration at the crystal-melt interface is calculated for various values of the rotation rate of the gravitational acceleration.

Mcfadden, G. B.↗
NASA NTRS1985

Juncture flow measurements using laser velocimetry

The turbulent flow in a simulated wing-fuselage juncture has been investigated using Laser Velocimetry. The juncture was simulated by a splitter plate with an unswept wing model mounted perpendicular to it. The use of sharp wing leading edge results in the elimination of separation ahead of leading edge, and hence the separation vortex. Secondary flow induced by Reynolds stress gradients is the main source of vortical motion in the juncture formed by the plate and the flat surface of the wing, whereas in the juncture formed by the plate and the circular arc surface of the wing, streamwise vorticity is generated because of the presence of lateral curvature in the flow field. The results of this study indicate that fillets in the juncture can improve the flow characteristics near the wing trailing edge.

Kubendran, L. R.↗
NASA NTRS1994

Schlieren Imaging Of An Aircraft In Flight

Technique for making schlieren images of airplanes and missiles in supersonic flight devised to help understand physics of compressible aerodynamic flows about complicated aircraft shapes. Technique also used to study far-field sonic booms. Data obtained from schlieren images useful in optimizing designs of prototype aircraft. Technique incorporates elements of focusing schlieren photography, astronomical photography, and streak photography. Using sun or moon as source of light, apparatus forms image revealing gradients of density in air flow.

Weinstein, Leonard M.↗
DOE OSTI2024

Energy Transport and Conversion Above a Bright Discrete Auroral Arc

Magnetically connected observations of particle distributions and luminosity from the Reimei spacecraft are used to examine energy transport and conversion occurring above a discrete auroral arc. By combining imaging and in situ measurements it is shown how transverse electromagnetic and kinetic energy fluxes measured along the spacecraft trajectory converge across geomagnetic field-lines into the acceleration region. It is shown how cross-field energy transport is facilitated by the formation of vortices along the length of the arc. From an integration over the vertical extent of the acceleration region it is shown that the transverse and field-aligned flow of energy into the arc locally supports the dissipation needed to power the electron acceleration observed. Estimates of gradients in electromagnetic and kinetic energy flows show how field-aligned electron energization in the acceleration region is supported by the divergence of Poynting flux along and across the background magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1970

Laminar wave train structure of collisionless magnetic slow shocks

The laminar wave train structure of collisionless magnetic slow shocks is investigated using two fluid hydromagnetics with ion cyclotron radius dispersion. For shock strengths less than the maximally strong switch-off shock, in the shock leading edge dispersive steepening forms a magnetic field gradient, while in the downstream flow dispersive propagation forms a trailing wave train; dispersion scale lengths are the ion inertial length if beta is smaller than 1 and the ion cyclotron radius if beta is greater than 1. In the switch-off slow shock leading edge, dispersion only produced rotations of the magnetic field direction; the gradient of the magnetic field magnitude, and hence the shock steepening length, is determined solely by resistive diffusion. The switch-off shock structure consists of a long trailing of magnetic rotations which are gradually damped by resistivity.

Coroniti, F. V.↗
NASA NTRS1988

Leakage predictions for Rayleigh-step, helium-purge seals

Rayleigh-step, helium purge, annular shaft seals, studied for use in liquid oxygen turbopumps, generate a hydrodynamic force that enables the seal to follow shaft perturbations. Hence, smaller clearances can be used to reduce seal leakage. FLOWCAL, a computer code developed by Mechanical Technology Incorporated, predicts gas flow rate through an annular seal with an axial pressure gradient. Analysis of a 50-mm Rayleigh-step, helium-purge, annular seal showed the flow rate increased axial pressure gradient, downstream pressure, and eccentricity ratio. Increased inlet temperature reduced leakage. Predictions made at maximum and minimum clearances (due to centrifugal and thermal growths, machining tolerances and + or - 2 percent uncertainty in the clearance measurement) placed wide boundaries on expected flow rates. The widest boundaries were set by thermal growth conditions. Predicted flow rates for a 50-mm Rayleigh-step, helium-purge, annular seal underestimated measured flow rates by three to seven times. However, the analysis did accurately predict flow rates for choked gas flow through annular seals when compared to flow rates measured in two other independent studies.

Proctor, Margaret P.↗
DOE OSTI2025

Vertical Characterization of Flow and Contaminant Concentration in New Hanford Extraction Wells

Effective groundwater remediation should focus on high-concentration, low-permeability zones within the aquifer, as these areas often present significant challenges to contaminant removal. Characterizing the vertical distribution of contamination in an aquifer is therefore critical for optimizing and monitoring pump-and-treat (P&T) systems to meet target cleanup timelines. Flow patterns driven by hydraulic gradients and variations in permeability are described to provide insight into contaminant mass removal from long-screened wells. The objective of this multi-phase work is to provide actionable recommendations for improving remedy performance monitoring and P&T operations, which ongoing evaluations will address. This report summarizes the preliminary results of vertical flow and contaminant distribution monitoring in new extraction wells at the Hanford Site.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Observations of Laminar Surface Shear Stress Patterns on the Sidewall of a Rectangular, Mach 3 Quiet Tunnel Using Liquid Crystal Coatings

A shear-stress--sensitive liquid crystal coating (LCC) was used to visualize the surface shear stress distribution on the flat sidewall of the MSU quiet supersonic tunnel as a function of flow stagnation pressure. Under conditions of quiet operation, the LCC color-change response indicated the existence of a nonuniform surface shear stress distribution. This shear pattern was characterized by an elongated, down st ream-point ing triangular region of relatively low shear, with its apex on and its axis coincident with, the sidewall centerline. This low-shear zone was bounded symmetrically by two regions of relatively higher shear; these high-shear zones originated within the concave-curvature portion of the nozzle, in the corners between the flat sidewall and the contoured upper and lower nozzle surfaces. A 3-D Navier-Stokes code was used to compute the pressure and surface shear distributions on the sidewall. Flow-expansion-induced transverse pressure gradients on the nozzle sidewall generated symmetric inflows from the corners towards the sidewall centerline; these inflows caused a thickening of the sidewall boundary layer along the centerline, resulting in lower shear stresses consistent with the liquid crystal results. Peripherally nonuniform laminar boundary layer development, and the associated stability of such complex 3-D flows, must be considered in quiet-tunnel applications using rectangular nozzles. A color video will be shown.

Reda, D.↗
NASA NTRS2017

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method↗
NASA NTRS1994

The Eady problem for a basic state with zero PV gradient but beta does not equal 0

The classic Eady problem is modified to include beta does not equal 0, but with the basic distributions of temperature and zonal flow adjusted to preserve zero meridional gradients of basic-state potential vorticity in the fluid interior. Much of the mathematical simplicity of the classic problem is retained; however, the results differ in important ways. Specifically, the instability now has a long-wave cutoff in addition to the traditional short-wave cutoff. The former is associated with the fact that the phase speeds of the edge waves begin to differ so much as wavenumber is reduced that the two edge waves can no longer interact in order to form unstable modes. For the unstable modes, this manifests itself in that the steering level for unstable modes is always below the middle of the fluid and approaches the lower boundary near the long-wave cutoff. Relatedly, the amplitude of the unstable geopotential perturbations is larger at the upper boundary than at the lower boundary. Finally, below the long-wave cutoff, one of the neutral waves has a phase speed that becomes increasingly easterly as wavenumber decreases. This allows a resonant response to planetary-scale stationary forcing.

Lindzen, Richard S.↗
NASA NTRS1980

Analysis and calculation of macrosegregation in a casting ingot. MPS solidification model. Volume 1: Formulation and analysis

The physical and numerical formulation of a model for the horizontal solidification of a binary alloy is described. It can be applied in an ingot. The major purpose of the model is to calculate macrosegregation in a casting ingot which results from flow of interdendritic liquid during solidification. The flow, driven by solidification contractions and by gravity acting on density gradients in the interdendritic liquid, was modeled as flow through a porous medium. The symbols used are defined. The physical formulation of the problem leading to a set of equations which can be used to obtain: (1) the pressure field; (2) the velocity field: (3) mass flow and (4) solute flow in the solid plus liquid zone during solidification is presented. With these established, the model calculates macrosegregation after solidification is complete. The numerical techniques used to obtain solution on a computational grid are presented. Results, evaluation of the results, and recommendations for future development of the model are given. The macrosegregation and flow field predictions for tin-lead, aluminum-copper, and tin-bismuth alloys are included as well as comparisons of some of the predictions with published predictions or with empirical data.

Maples, A. L.↗