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 721 records · Page 40

Comparison of Boundary Layer Similarity Transformations for High Mach Number Flows

Reynolds averaged Navier-Stokes (RANS) simulations of the turbulent boundary layer flow on a flat plate are performed across a Mach number range of 0.3 up to 16 using the computational fluid dynamics (CFD) software, VULCAN-CFD. The simulation results are used to evaluate the ability of the Van Driest transformation and the transformation of Trettel and Larsson to collapse boundary layer velocity profiles under various flow conditions across the range of Mach numbers. The flat platewall boundary conditions and leading edge geometry are varied from adiabatic to isothermal and sharp to blunted to reveal the physics of how these effects impact the performance of the chosen flow transformations. Results indicate that the transformation of Trettel and Larsson produces a better collapse of velocity profile data than the Van Driest transformation regardless of Mach number, wall boundary condition, or leading edge geometry. In addition, the velocity gradients obtained from the transformation of Trettel and Larsson match the normalized untransformed velocity gradients much more closely than those obtained using the Van Driest transformation.

DiGregorio, Nicholas J.↗

The 3D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a three-dimensional (3D) Navier-Stokes flow solver. A two-equation turbulence model, and a porous bleed model based on unchoked bleed hole discharge coefficients were used. Comparisons were made with experimental data, inviscid theory, and two-dimensional Navier-Stokes analyses. The main objective was to gain insight into the inlet fluid dynamics. Examination of the computational results along with the experimental data suggest that the cowl shock-sidewall boundary layer interaction near the leading edge caused a substantial separation in the wind tunnel inlet model. As a result, the inlet performance may have been compromised by increased spillage and higher bleed mass flow requirements. The internal flow contained substantial waves that were not in the original inviscid design. 3D effects were fairly minor for this inlet at on-design conditions. Navier-Stokes analysis appears to be an useful tool for gaining insight into the inlet fluid dynamics. It provides a higher fidelity simulation of the flowfield than the original inviscid design, by taking into account boundary layers, porous bleed, and their interactions with shock waves.

Mizukami, M.↗

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗

CFD Simulations of Molten Salts Within Plasma-Facing Components with Twisted Tape Inserts Under One-Sided–Heated Conditions

For liquid immersion blankets (LIBs), the heat removal within different plasma-facing components (PFCs) requires substantial heat removal capabilities that could benefit from the usage of heat transfer enhancements. Of consideration are twisted tape insert heat transfer enhancements that are ribbon-like metal pieces inserted within cooling channels that increase fluid mixing, resulting in higher heat transfer and pressure drop. Unlike traditional applications, PFCs such as divertor coolant tubes experience a significantly higher heat flux on one side versus the other or entirely one-sided heating. These conditions are not well studied for molten salts where the higher Prandtl numbers characteristic of the fluids are significantly higher than water, gases, and liquid metals. As an effort to study this knowledge gap, large eddy simulations were done with Prandtl numbers representative of FLiBe at proposed operating conditions. In this study, the effects of applying a surface heat flux to one side of a tube versus the entire tube with a twisted tape insert were investigated for the heat transfer coefficient and effects on the temperature profile. The findings of the study suggest that for Prandtl number fluids above unity, the Nusselt number was largely unaffected by the change in boundary condition. Both boundary conditions showed good agreement with the Manglik and Bergles Nusselt number correlation for tight-fitting twisted tapes. Finally, the outcomes of this study will be useful for future design studies of LIBs to provide fundamental guidance implementing twisted tape inserts and their performance in one-sided–heated conditions.

molten salts↗

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Periaux, J.↗

Vortex simulation of an inviscid shear layer

The accuracy of the vortex-blob method was tested by simulating a free-shear-layer instability, Kirchhoff's elliptical vortex, and a circular vortex. The main numerical parameters in the vortex-blob method are the density of the vortices, and the distribution of vorticity within each vortex core. The growth rate of a periodic unstable mode of the shear layer was calculated numerically and compared with the exact result. The error is only a few percent for about 10 rows of vortex blobs. The error is reduced by decreasing the spacing between vortices and, correspondingly, the core size. In the simulation of the motion of the elliptical vortex, the rotation of the boundary, without change of shape, and the circular particle paths of the vortical fluid were well simulated. For the circular vortex, optimum sets of parameters were obtained by comparing them with the exact velocity. The results are consistent with convergence theories of the vortex-blob method. In particular, second-order convergence is observed with a Gaussian core from velocity calculation.

Nakamura, Y.↗

Application of the MHD energy principle to magnetostatic atmospheres

The MHD energy principle is applied to the stability of a magnetized atmosphere which is bounded below by much denser fluid, as is the solar corona. The two fluids are treated as ideal; the approximation is consistent with the energy principle, and the dynamical conditions that must hold at a fluid-fluid interface are used to show that if vertical displacements of the lower boundary are premitted, then the lower atmosphere must be perturbed as well. However, displacements which do not perturb the coronal boundary can be properly treated as isolated perturbations of the corona alone.

Zweibel, E. G.↗

Application of the MHD energy principle to magnetostatic atmospheres

The MHD energy principle is applied to the stability of a magnetized atmosphere which is bounded below by much denser fluid, as is the solar corona. The two fluids are treated as ideal; the approximation is consistent with the energy principle, and the dynamical conditions that must hold at a fluid-fluid interface are used to show that if vertical displacements of the lower boundary are permitted, then the lower atmosphere must be perturbed as well. However, displacements which do not perturb the coronal boundary can be properly treated as isolated perturbations of the corona alone.

Zweibel, E. G.↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The influence of bulk turbulence and boundary layer thickness on the secondary flow development in a square, 90 degree turning duct was investigated. A three-dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity at six cross-planes in the duct. The results from this investigation, with entrance boundary layer thickness of 20 percent, were compared with the thin boundary layer results documented in NASA CR-174811. The axial velocity profiles, cross-flow velocities, and turbulence intensities were compared and evaluated with regard to the influence of bulk turbulence intensity and boundary layer thickness, and the influence was significant. The results of this investigation expand the 90 degree curved duct experimental data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, R. A.↗

Noble metals in mid-ocean ridge volcanism: A significant fractionation of gold with respect to platinum group metals

Hydrothermal precipitates, black smoker particulate, and massive sulphide dredge samples from the Explorer Ridge on the Juan de Fuca Plate and the TAG hydrothermal area on the Mid-Atlantic Ridge were analyzed for selected noble metals including Au, Ir and Pd by radiochemical neutron activation analysis. The preliminary results indicate that gold contents may reach the ppm range although values in the neighborhood of 100 to 200 ppb are more typical. The platinum group elements (PGE) represented by Ir and Pd are typically less than 0.02 ppb and less than 2 ppb respectively. These abundances represent a significant enrichment of gold relative to the PGE in comparison with average noble metal abundances in mid-ocean ridge basalts (MORB). A partial explanation of this distinctive fractionation can be found in the concepts of sulfur-saturation of basic magma in mid-ocean ridge (MOR) settings, and the origin of MOR hydrothermal fluids. Experimental and petrological data suggest that MORBs are sulfur-saturated at the time of magma generation and that an immiscible sulfide component remains in the mantle residue. Hence, MORBs are noble metal-poor, particularly with respect to PGE. Consequently, black smoker fluids can be expected to reflect the low Ir and Pd contents of the rock column. The average Au content of MORB is 1.3 ppb, and so the rock column is not significantly enriched in Au. The generation of fluids which precipitate solids with 200 ppb Au is apparently dependent on highly efficient fluid chemistry to mobilize Au from the rock column, high Au solubility in seawater hydrothermal fluids and efficient precipitation mechanisms to coprecipitate Au on Fe, Zn and Cu sulfides. Significant differences in these parameters appear to be the ultimate cause of the strong Au-PGE fractionation in the MOR setting. It does not appear from the current data base that MOR hydrothermal fluids are significant contributors to the Ir enrichment seen in Cretaceous-Tertiary boundary sediments.

Crocket, James H.↗

Considerations on the moving contact-line singularity, with application to frictional drag on a slender drop

It has previously been shown that the no-slip boundary conditions leads to a singularity at a moving contact line and that this presumes some form of slip. Present considerations on the energetics of slip due to shear stress lead to a yield stress boundary condition. A model for the distortion of the liquid state near solid boundaries gives a physical basis for this boundary condition. The yield stress condition is illustrated by an analysis of a slender drop rolling down an incline. That analysis provides a formula for the frictional drag resisting the drop movement. With the present boundary condition, the length of the slip region becomes a property of the fluid flow.

Durbin, P. A.↗

Simulating transitional flow and heat transfer over the flat plate and circular cylinder using a K-epsilon turbulence model

A boundary layer computer code, called TEXSTAN, has been used to simulate external, transitional, fluid flow and heat transfer over the flat plate and circular cylinder. Turbulence modeling was accomplished by using a low-Reynolds number K-epsilon turbulence model. The study was made for free-stream turbulence intensities less than 15 percent. The comparison of the flat plate heat transfer predictions to experimental data showed good agreement, for the location of transition as well as in the heat transfer coefficient. The typical discrepancy in predicting the heat transfer coefficient was less than 5 percent of the measured value.

Sullivan, T. J.↗

Unsteady panel method for flows with multiple bodies moving along various paths

A potential flow based three-dimensional panel method was modified to treat time dependent conditions in which several submerged bodies can move within the fluid along different trajectories. This modification was accomplished by formulating the momentary solution in an inertial frame-of-reference, attached to the undisturbed stationary fluid. Consequently, the numerical interpretation of the multiple-body, solid-surface boundary condition and the viscous wake rollup was considerably simplified. The unsteady capability of this code was validated by comparing computed and experimental results for a finite wing undergoing pitch oscillations. In order to demonstrate the multicomponent capability, computations were made for two wings following closely intersecting paths (e.g., to avoid mid air collisions) and for a flow field with relative rotation (e.g., helicopter-rotor/fuselage interaction). Results were compared to experimental data when such data was available.

Richason, Thomas F.↗

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗

Unsteady panel method for flows with multiple bodies moving along various paths

A potential flow based three-dimensional panel method was modified to treat time-dependent conditions in which several submerged bodies can move within the fluid along different trajectories. This modification was accomplished by formulating the momentary solution in an inertial frame of reference, attached to the undisturbed stationary fluid. Consequently, the numerical interpretation of the multiple-body, solid-surface boundary condition and the viscous wake rollup was considerably simplified. The usteady capability of this code was calibrated and validated by comparing computed results with closed-form analytical results available for an airfoil, which was impulsively set into a constant speed forward motion. To demonstrate the multicomponent capability, computations were made for two wings following closely intersecting paths (i.e., simulations aimed at avoiding mid-air collisions) and for a flowfield with relative rotation (i.e., the case of a helicopter rotor rotating relative to the fuselage). Computed results for the cases were compared to experimental data, when such data was available.

Richardson, Thomas F.↗

An Experimental Study of the Near Field Region of a Free Jet with Passive Mixing Tabs

An experimental study was performed to determine the flow characteristics of a tabbed free jet. Results were acquired in the near field (nominally 2 tab widths upstream to 2 tab widths downstream of the exit plane) of a tabbed jet. Upstream pressure results showed static pressure distributions in both the x-and y-directions along the top surface of the tunnel. Hot-wire measurements showed rapid expansion of the core fluid into the ambient region. Two counter rotating regions of streamwise vorticity were shown on each side of the primary tab. An enhancement of the tabbed jet concept was proposed and tested. Specifically, two tabs, half the scale of the primary tab, were added to the primary tab to provide attachment surfaces for the normally occurring ejection of fluid. The secondary tabs caused a slight increase in the streamwise vorticity created from the upstream static pressure gradient while significantly increasing the re-oriented boundary layer vorticity. The combined pumping effect of the two counter rotating regions of vorticity caused a significant increase in the transport of the jet core fluid into the surrounding region.

Bohl, D. G.↗

Heat Transfer in a Liquid Metal Flowing Turbulently through a Channel with a Step Function Boundary Temperature

An analytical heat transfer solution is derived and evaluated for the general case of a turbulently flowing liquid metal which suddenly encounters a step-function boundary temperature in a channel system. Local Nusselt moduli, dimensionless mixed-mean fluid temperatures, and arithmetic-mean Nusselt moduli are given as functions of Reynolds and Prandtl moduli and a dimensionless axial-distance modulus. These solutions are compared with known solutions of more specific systems as well as with a set of experimental liquid-metal heat transfer data for a thermal entrance region.

Poppendiek, H. F.↗