Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID BOUNDARY”

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 523 records · Page 29

The effect of gravity modulation on thermosolutal convection in an infinite layer of fluid

The effect of time-periodic vertical gravity modulation on the onset of thermosolutal convection in an infinite horizontal layer with stress-free boundaries is investigated using Floquet theory for the linear stability analysis. Situations for which the fluid layer is stably stratified in either the fingering or diffusive regimes of double-diffusive convection are considered. Results are presented both with and without steady background acceleration. Modulation may stabilize an unstable base solution or destabilize a stable base solution. In addition to synchronous and subharmonic response to the modulation frequency, instability in the double diffusive system can occur via a complex conjugate mode. In the diffusive regime, where oscillatory onset occurs in the unmodulated system, regions of resonant instability occur and exhibit strong coupling with the unmodulated oscillatory frequency. The response to modulation of the fundamental instability of the unmodulated system is described both analytically and numerically; in the double-diffusive system this mode persists under subcritical conditions as a high-frequency lobe.

Saunders, B. V.↗

A continuum theory of a lubrication problem with solid particles

The governing equations for a two-dimensional lubrication problem involving the mixture of a Newtonian fluid with solid particles at an arbitrary volume fraction are developed using the theory of interacting continuua (mixture theory). The equations take the interaction between the fluid and the particles into consideration. Provision is made for the possibility of particle slippage at the boundaries. The equations are simplified assuming that the solid volume fraction varies in the sliding direction alone. Equations are solved for the velocity of the fluid phase and that of the solid phase of the mixture flow in the clearance space of an arbitrary shaped bearing. It is shown that the classical pure fluid case can be recovered as a special case of the solutions presented. Extensive numerical solutions are presented to quantify the effect of particulate solid for a number of pertinent performance parameters for both slider and journal bearings. Included in the results are discussions on the influence of particle slippage on the boundaries as well as the role of the interacting body force between the fluid and solid particles.

Dai, Fuling↗

Eigenvalue Sensitivity Computations for Linear Stability Theory

To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.

Boundary Layer Transition↗

Magnetic Field Effect on Convective Fluid Flow in a Vertical Bridgman System

We studied the magnetic field effects on convective fluid flow in a differentially heated HgCdTe Bridgman growth system with linearly stratified vertical HgTe concentration. We followed J. E. Hart's system set-up, boundary conditions, and formulation. We obtained analytical steady parallel flow solutions for fluid flow's velocity and solutal concentration. These solutions depend only on solutal Rayleigh number. Our analytical solutions had different expressions from that of Hart's, but had similar numerical simulation results for the mean fluid velocity and solutal concentration gradient as a function of solutal Rayleigh number. We studied the effects of a transverse magnetic field on the fluid flow of the Bridgman growth systems both on the ground and in the space. The simulation results and comparison with experimental results will be presented.

Wang, J. C.↗

Investigation of geomagnetic field forecasting and fluid dynamics of the core

It was established that the total absolute magnetic flux crossing the core- mantle boundary has been a constant of the core motion for the last 50 years. This provides a scalar constraint that could be added to the geometric modelling procedure. The GSFC 8 8/80 model is being evaluated. The absolute magnetic flux linking the CMB to that model was plotted as a function of time during the span covered by the data, and increasing truncation level. The inclusion of the standard error of each Gauss coefficient derived from the statistics of fit in the GSFC 9/80 model is useful. The magnitude and sense (upwelling or downe welling) of vertical fluid motion adjacent to the core-mantle boundary was calculated using the model. Standard errors were found to be sufficiently small at all but one or two of the 40 or more critical points of B sub r. They do not nearly overlap the value gamma u/gamma r = 0. It is concluded that the core is upwelling and downwelling at an observationally detectable level.

Benton, E. R.↗

Application of Micro-ramp Flow Control Devices to an Oblique Shock Interaction

Tests are planned in the 15cm x 15cm supersonic wind tunnel at NASA Glenn to demonstrate the applicability of micro-ramp flow control to the management of shock wave boundary layer interactions. These tests will be used as a database for computational fluid dynamics (CFD) validation and Design of Experiments (DoE) design information. Micro-ramps show potential for mechanically simple and fail-safe boundary layer control.

Hirt, Stefanie↗

Viscous aerodynamic design using the adjoint variable approach

The use of classical optimal control methods, in particular variational methods, to solve the airfoil optimization problem, by deriving a set of adjoint (costate) equations and boundary conditions has already been done for inviscid (potential and Euler flows) and two dimensional, steady state, incompressible flow governed by the Navier-Stokes equations. The interior and boundary terms of the volume integral have been derived (in this work) for the steady Navier-Stokes equations in three dimensions for a viscous, compressible heat conducting fluid. This can be used to derive the adjoint equations and numerical boundary conditions for general classes of problems and hence paves the way for a solution to the aerodynamic optimization problem for compressible viscous flows. The next steps to the realization of that goal are projected as below. The usual square integral pressure functional as an objective function is being replaced by a more realistic drag functional subject to a lift constraint. The feasibility of attempting the more difficult time dependent problem is being investigated. It remains to get the full system of adjoint equations and boundary conditions with the new functional. The state and adjoint equations must be discretized and coded. An appropriate optimization program must be used (steepest descents seems inadequate) and various known airfoil shapes should be recovered in test cases of the computer program.

DeRise, George↗

Delay of Turbulent Boundary Layer Detachment by Mechanical Excitation: Application to Rearward-facing Ramp

A vane oscillating about a fixed point at the inlet to a two-dimensional 20 deg rearward-facing ramp proved effective in delaying the detachment of a turbulent boundary layer. Flow-field, surface static pressure, and smoke-wire flow visualization measurements were made. Surface pressure coefficient distributions revealed that two different effects occurred with axial distance along the ramp surface. The surface pressure coefficient varied as a complex function of the vane oscillation frequency and its trailing edge displacement amplitude; that is, it varied as a function of the vane oscillation frequency throughout the entire range of frequencies covered during the test, but it varied over only a limited range of the trailing edge displacement amplitudes covered.The complexity of these findings prompted a detailed investigation, the results of which revealed a combination of phenomena that explain qualitatively how the mechanically generated, periodic, sinusoidal perturbing signal produced by the oscillating vane reacts with the fluid flow to delay the detachment of a turbulent boundary layer experiencing transitory detachment.

McKinzie, Daniel J., Jr.↗

NASTRAN hydroelastic modal studies. Volume 1: Introduction, theory, and results

The implementation of a general three-dimensional hydroelastic capability in the NASTRAN computer program is described. A method for analyzing the combined mode shapes of arbitrary fluid and structure finite element models was developed. The fluid was modeled with three-dimensional solid elements with options for tetrahedron, wedge, and hexahedron shapes. The elements were then connected to fluid grid points which defined the pressure in the fluid at the specified location. Matrices were formulated to provide efficient solutions for large-order problems. The structure matrices were processed separately and could be reduced using matrix condensation procedures or through a modal formulation using the normal modes of the empty structure as solution coordinates. The fluid/structure interface and the free surface were defined by the user with special NASTRAN boundary elements. A special purpose mesh generator program was provided to generate the actual NASTRAN data cards for the fluid, the structure, and the boundary elements for typical tank-type models.

Source record↗

Heat transfer to a full-coverage film-cooled surface with 30 degree slant-hole injection

Heat transfer behavior was studied in a turbulent boundary layer with full coverage film cooling through an array of discrete holes and with injection 30 deg to the wall surface in the downstream direction. Stanton numbers were measured for a staggered hole pattern with pitch-to-diameter ratios of 5 and 10, an injection mass flux ratio range of 0.1 to 1.3, and a range of Reynolds number Re sub x of 150,000 to 5 million. Air was used as the working fluid, and the mainstream velocity varied from 9.8 to 34.2 m/sec (32 to 112 ft/sec). The data were taken for secondary injection temperature equal to the wall temperature and also equal to the mainstream temperature. The data may be used to obtain Stanton number as a continuous function of the injectant temperature by use of linear superposition theory. The heat transfer coefficient is defined on the basis of a mainstream-to-wall temperature difference. This definition permits direct comparison of performance between film cooling and transpiration cooling. A differential prediction method was developed to predict the film cooling data base. The method utilizes a two-dimensional boundary layer program with routines to model the injection process and turbulence augmentation. The program marches in the streamwise direction, and when a row of holes is encountered, it stops and injects fluid into the boundary layer. The turbulence level is modeled by algebraically augmenting the mixing length, with the augmentation keyed to a penetration distance for the injected fluid.

Crawford, M. E.↗

Cell volume control at a surface for three-dimensional grid generation packages

An alternate method of calculating the cell size for orthogonality control in the solution of Poisson's 3D space equations is presented. The method provides the capability to enforce a better initial guess for the grid distribution required for boundary layer resolution. This grid point distribution is accomplished by enforcing grid spacing from a grid block boundary where orthogonality is required. The actual grid spacing or cell size for that boundary is determined by the two or four adjacent boundaries in the grid block definition, which are two dimensional grids. These two dimensional grids are in turn defined by the user using insight into the flow field and boundary layer characteristics. The adjoining boundaries are extended using a multifunctional blending scheme, with user control of the blending and interpolating functions to be used. This grid generation procedure results in an enhanced computational fluid dynamics calculation by allowing a quicker resolution of the configuration's boundary layer and flow field and by limiting the number of grid re-adaptations. The cell size specification calculation was applied to a variety of configurations ranging from axisymmetric to complex three-dimensional configurations. Representative grids are shown for the Space Shuttle and the Langley Lifting Body (HL-20).

Alter, Stephen J.↗

Improved Boundary Conditions for Cell-centered Difference Schemes

Cell-centered finite-volume (CCFV) schemes have certain attractive properties for the solution of the equations governing compressible fluid flow. Among others, they provide a natural vehicle for specifying flux conditions at the boundaries of the physical domain. Unfortunately, they lead to slow convergence for numerical programs utilizing them. In this report a method for investigating and improving the convergence of CCFV schemes is presented, which focuses on the effect of the numerical boundary conditions. The key to the method is the computation of the spectral radius of the iteration matrix of the entire demoralized system of equations, not just of the interior point scheme or the boundary conditions.

VanderWijngaart, Rob F.↗

Boundary layer development as a function of chamber pressure in the NASA Lewis 1030:1 area ratio rocket nozzle

Through the use of theoretical predictions of fluid properties and experimental heat transfer and thrust measurements, the zones of laminar, transitional, and turbulent boundary layer flow were defined for the NASA Lewis 1030:1 area ratio rocket nozzle. Tests were performed on the nozzle at chamber pressures from 350 to 100 psia. For these conditions, the throat diameter Reynolds numbers varied from 300,000 to 1 million. The propellants used were gaseous hydrogen and gaseous oxygen. Thrust measurements and nozzle outer wall temperature measurements were taken during the 3-sec test runs. Comparison of experimental heat transfer and thrust data with the corresponding predictions from the Two-Dimensional Kinetics (TDK) nozzle analysis program indicated laminar flow in the nozzle at a throat diameter Reynolds number of 320,000 or chamber pressure of 360 psia. Comparison of experimental and predicted heat transfer data indicated transitional flow up to and including a chamber pressure of 1000 psia. Predicted values of the axisymmetric acceleration parameter within the convergent and divergent nozzle were consistent with the above results. Based upon an extrapolation of the heat transfer data and predicted distributions of the axisymmetric acceleration parameter, transitional flow was predicted up to a throat diameter Reynolds number of 220,000 or 2600-psia chamber pressure. Above 2600-psia chamber pressure, fully developed turbulent flow was predicted.

Smith, Tamara A.↗

Boundary layer development as a function of chamber pressure in the NASA Lewis 1030:1 area ratio rocket nozzle

Through the use of theoretical predictions of fluid properties and experimental heat transfer and thrust measurements, the zones of laminar, transitional, and turbulent boundary layer flow were defined for the NASA Lewis 1039:1 area ratio rocket nozzle. Tests were performed on the nozzle at chamber pressures from 350 to 100 psia. For these conditions, the throat diameter Reynolds numbers varied from 300,000 to 1 million. The propellants used were gaseous hydrogen and gaseous oxygen. Thrust measurements and nozzle outer wall temperature measurements were taken during the 3-sec test runs. Comparison of experimental heat transfer and thrust data with the corresponding predictions from the Two-Dimensional Kinetics (TDK) nozzle analysis program indicated laminar flow in the nozzle at a throat diameter Reynolds number of 320,000 or chamber pressure of 360 psia. Comparison of experimental and predicted heat transfer data indicated transitional flow up to and including a chamber pressure of 1000 psia. Predicted values of the axisymmetric acceleration parameter within the convergent and divergent nozzle were consistent with the above results. Based upon an extrapolation of the heat transfer data and predicted distributions of the axisymmetric acceleration parameter, transitional flow was predicted up to a throat diameter Reynolds number of 220,000 or 2600-psia chamber pressure. Above 2600-psia chamber pressure, fully developed turbulent flow was predicted.

Smith, Tamara A.↗

The Effect of an Isogrid on Cryogenic Propellant Behavior and Thermal Stratification

All models for thermal stratification available in the presentation are derived using smooth, flat plate laminar and turbulent boundary layer models. This study examines the effect of isogrid (roughness elements) on the surface of internal tank walls to mimic the effects of weight-saving isogrid, which is located on the inside of many rocket propellant tanks. Computational Fluid Dynamics (CFD) is used to study the momentum and thermal boundary layer thickness for free convection flows over a wall with generic roughness elements. This presentation makes no mention of actual isogrid sizes or of any specific tank geometry. The magnitude of thermal stratification is compared for smooth and isogrid-lined walls.

Oliveira, Justin↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗