Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary layer 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 811 records · Page 45

Improved ceramic slip casting technique

A primary concern in modern fluid dynamics research is the experimental verification of computational aerothermodynamic codes. This research requires high precision and detail in the test model employed. Ceramic materials are used for these models because of their low heat conductivity and their survivability at high temperatures. To fabricate such models, slip casting techniques were developed to provide net-form, precision casting capability for high-purity ceramic materials in aqueous solutions. In previous slip casting techniques, block, or flask molds made of plaster-of-paris were used to draw liquid from the slip material. Upon setting, parts were removed from the flask mold and cured in a kiln at high temperatures. Casting detail was usually limited with this technique -- detailed parts were frequently damaged upon separation from the flask mold, as the molded parts are extremely delicate in the uncured state, and the flask mold is inflexible. Ceramic surfaces were also marred by 'parting lines' caused by mold separation. This adversely affected the aerodynamic surface quality of the model as well. (Parting lines are invariably necessary on or near the leading edges of wings, nosetips, and fins for mold separation. These areas are also critical for flow boundary layer control.) Parting agents used in the casting process also affected surface quality. These agents eventually soaked into the mold, the model, or flaked off when releasing the case model. Different materials were tried, such as oils, paraffin, and even an algae. The algae released best, but some of it remained on the model and imparted an uneven texture and discoloration on the model surface when cured. According to the present invention, a wax pattern for a shell mold is provided, and an aqueous mixture of a calcium sulfate-bonded investment material is applied as a coating to the wax pattern. The coated wax pattern is then dried, followed by curing to vaporize the wax pattern and leave a shell mold of the calcium sulfate-bonded investment material. The shell mold is cooled to room temperature, and a ceramic slip is poured therein. After a ceramic shell of desired thickness has set up in the shell mold, excess ceramic slip is poured out. While still wet, the shell mold is peeled from the ceramic shell to expose any delicate or detailed parts, after which the ceramic shell is cured to provide a complete, detailed, precision ceramic article without parting lines.

Buck, Gregory M.↗

Efficiency and reliability enhancements in propulsion flowfield modeling

The implementation of traditional CFD algorithms in practical propulsion related flowfields often leads to dramatic reductions in efficiency and/or robustness. The present research is directed at understanding the reasons for this deterioration and finding methods to circumvent it. Work to date has focussed on low Mach number regions, viscous dominated regions, and high grid aspect ratios. Time derivative preconditioning, improved definition of the local time stepping, and appropriate application of boundary conditions are employed to decrease the required time to obtain a solution, while maintaining accuracy. A number of cases having features typical of rocket engine flowfields are computed to demonstrate the improvement over conventional methods. These cases include laminar and turbulent high Reynolds number flat plate boundary layers, flow over a backward-facing step, a diffusion flame, and wall heat-flux calculations in a turbulent converging-diverging nozzle. Results from these cases show convergence that is virtually independent of the local Mach number and the grid aspect ratio, which translates to a convergence speed-up of up to several orders of magnitude over conventional algorithms. Current emphasis is in extending these results to three-dimensional flows with highly stretched grids.

Buelow, Philip E. O.↗

Relationship between transition and modes of instability in supersonic boundary layers

The relationship between the predicted transition location and the first and second modes of instability in two-dimensional supersonic boundary-layer flow on a flat plate is examined. Linear stability theory and the N-factor criterion are used to predict transition location. The effect of heat transfer is also studied; the results demonstrate that the transition reversal phenomenon can be explained by the opposite effect of cooling on the first and second modes of instability. Compressibility of destabilizing at free-stream Mach numbers of 2 to 3.5. The predicted transition location is due to the oblique first modes of instability, up to free-stream Mach numbers between 6 and 6.5. At higher Mach numbers, the predicted transition location is due to a combination of two-dimensional first and second modes of instability.

Masad, Jamal A.↗

Performance of Renormalization Group Algebraic Turbulence Model on Boundary Layer Transition Simulation

The RNG-based algebraic turbulence model, with a new method of solving the cubic equation and applying new length scales, is introduced. An analysis is made of the RNG length scale which was previously reported and the resulting eddy viscosity is compared with those from other algebraic turbulence models. Subsequently, a new length scale is introduced which actually uses the two previous RNG length scales in a systematic way to improve the model performance. The performance of the present RNG model is demonstrated by simulating the boundary layer flow over a flat plate and the flow over an airfoil.

Ahn, Kyung H.↗

Boundary-layer receptivity due to distributed surface imperfections of a deterministic or random nature

Acoustic receptivity of a Blasius boundary layer in the presence of distributed, 2D surface irregularities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface irregularities, only a narrow band of Fourier components can lead to an efficient conversion of the acoustic input at any given frequency to an unstable eigenmode of the boundary-layer flow. Random irregularities which are spatially homogeneous in nature, and also possess a continuous spectrum in the streamwise direction, lead to instability amplitudes that are intermediate to those caused by the periodic and isolated irregularities, respectively, being, in fact, of the same order as the geometrical mean of the amplitudes in the latter two cases. A physical explanation for these asymptotic scalings is given, in addition to providing an analytical expression for the expected value of the instability amplitude for an ensemble of statistically irregular surfaces with random phase distributions.

Choudhari, Meelan↗

Scalability study of parallel spatial direct numerical simulation code on IBM SP1 parallel supercomputer

The implementation and the performance of a parallel spatial direct numerical simulation (PSDNS) code are reported for the IBM SP1 supercomputer. The spatially evolving disturbances that are associated with laminar-to-turbulent in three-dimensional boundary-layer flows are computed with the PS-DNS code. By remapping the distributed data structure during the course of the calculation, optimized serial library routines can be utilized that substantially increase the computational performance. Although the remapping incurs a high communication penalty, the parallel efficiency of the code remains above 40% for all performed calculations. By using appropriate compile options and optimized library routines, the serial code achieves 52-56 Mflops on a single node of the SP1 (45% of theoretical peak performance). The actual performance of the PSDNS code on the SP1 is evaluated with a 'real world' simulation that consists of 1.7 million grid points. One time step of this simulation is calculated on eight nodes of the SP1 in the same time as required by a Cray Y/MP for the same simulation. The scalability information provides estimated computational costs that match the actual costs relative to changes in the number of grid points.

Hanebutte, Ulf R.↗

The Earth's magnetosphere is 165 R(sub E) long: Self-consistent currents, convection, magnetospheric structure, and processes for northward interplanetary magnetic field

The subject of this paper is a self-consistent, magnetohydrodynamic numerical realization for the Earth's magnetosphere which is in a quasi-steady dynamic equilibrium for a due northward interplanetary magnetic field (IMF). Although a few hours of steady northward IMF are required for this asymptotic state to be set up, it should still be of considerable theoretical interest because it constitutes a 'ground state' for the solar wind-magnetosphere interaction. Moreover, particular features of this ground state magnetosphere should be observable even under less extreme solar wind conditions. Certain characteristics of this magnetosphere, namely, NBZ Birkeland currents, four-cell ionospheric convection, a relatively weak cross-polar potential, and a prominent flow boundary layer, are widely expected. Other characteristics, such as no open tail lobes, no Earth-connected magnetic flux beyond 155 R(sub E) downstream, magnetic merging in a closed topology at the cusps, and a 'tadpole' shaped magnetospheric boundary, might not be expected. In this paper, we will present the evidence for this unusual but interesting magnetospheric equilibrium. We will also discuss our present understanding of this singular state.

Fedder, J. A.↗

CFD Study of Turbo-Ramjet Interactions in Hypersonic Airbreathing Propulsion System

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.

Chang, Ing↗

Water Vapor Effects on Silica-Forming Ceramics

Silica-forming ceramics such as SiC and Si3N4 are proposed for applications in combustion environments. These environments contain water vapor as a product of combustion. Oxidation of silica-formers is more rapid in water vapor than in oxygen. Parabolic oxidation rates increase with the water vapor partial pressure with a power law exponent value close to one. Molecular water vapor is therefore the mobile species in silica. Rapid oxidation rates and large amounts of gases generated during the oxidation reaction in high water vapor pressures may result in bubble formation in the silica and nonprotective scale formation. It is also shown that silica reacts with water vapor to form Si(OH)4(g). Silica volatility has been modeled using a laminar flow boundary layer controlled reaction equation. Silica volatility depends on the partial pressure of water vapor, the total pressure, and the gas velocity. Simultaneous oxidation and volatilization reactions have been modeled with paralinear kinetics.

Opila, E. J.↗

On the Connection Between One-and Two-Equation Models of Turbulence

A formalism will be presented that allows the transformation of two-equation eddy viscosity turbulence models into one-equation models. The transformation is based on an assumption that is widely accepted over a large range of boundary layer flows and that has been shown to actually improve predictions when incorporated into two-equation models of turbulence. Based on that assumption, a new one-equation turbulence model will be derived. The new model will be tested in great detail against a previously introduced one-equation model and against its parent two-equation model.

Menter, F. R.↗

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.↗

Numerical Simulation of Roughness-Induced Transient Growth in a Laminar Boundary Layer

Numerical simulations are used to examine the roughness-induced transient growth in a laminar boundary-layer flow. Based on the spectral element method, these simulations model the stationary disturbance field associated with a nonsmooth roughness geometry, such as the spanwise periodic array of circular disks used by White and co-workers during a series of wind tunnel experiments at Case Western Reserve University. Besides capturing the major trends from the recent measurements by White and Ergin, the simulations provide additional information concerning the relative accuracy of the experimental findings derived from two separate wall-finding procedures. The paper also explores the dependence of transient growth on geometric characteristics of the roughness distribution, including the height and planform shape of the roughness element and the ratio of roughness due to spacing between an adjacent pair of elements. Results are used for a preliminary assessment of the differences between recently reported theoretical results of Tumin and Reshotko and the measurements by White and Ergin.

Fischer, Paul↗

Scalability of Parallel Spatial Direct Numerical Simulations on Intel Hypercube and IBM SP1 and SP2

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube and IBM SP1 and SP2 parallel computers is documented. Spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows are computed with the PSDNS code. The feasibility of using the PSDNS to perform transition studies on these computers is examined. The results indicate that PSDNS approach can effectively be parallelized on a distributed-memory parallel machine by remapping the distributed data structure during the course of the calculation. Scalability information is provided to estimate computational costs to match the actual costs relative to changes in the number of grid points. By increasing the number of processors, slower than linear speedups are achieved with optimized (machine-dependent library) routines. This slower than linear speedup results because the computational cost is dominated by FFT routine, which yields less than ideal speedups. By using appropriate compile options and optimized library routines on the SP1, the serial code achieves 52-56 M ops on a single node of the SP1 (45 percent of theoretical peak performance). The actual performance of the PSDNS code on the SP1 is evaluated with a "real world" simulation that consists of 1.7 million grid points. One time step of this simulation is calculated on eight nodes of the SP1 in the same time as required by a Cray Y/MP supercomputer. For the same simulation, 32-nodes of the SP1 and SP2 are required to reach the performance of a Cray C-90. A 32 node SP1 (SP2) configuration is 2.9 (4.6) times faster than a Cray Y/MP for this simulation, while the hypercube is roughly 2 times slower than the Y/MP for this application. KEY WORDS: Spatial direct numerical simulations; incompressible viscous flows; spectral methods; finite differences; parallel computing.

Joslin, Ronald D.↗

Direct Simulation of Evolution and Control of Three-Dimensional Instabilities in Attachment-Line Boundary Layers

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic- source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in at-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.

Joslin, Ronald D.↗

Optimal disturbances in boundary layers subject to streamwise pressure gradient

An analysis of the optimal 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-Scan 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.

Ashpis, David E.↗

Turbulent Aeroheating Testing of Mars Science Laboratory Entry Vehicle in Perfect-Gas Nitrogen

An experimental investigation of turbulent aeroheating on the Mars Science Laboratory entry vehicle heat shield has been conducted in the Arnold Engineering Development Center Hypervelocity Wind Tunnel No. 9. Testing was performed on a 6-in. (0.1524 m) diameter MSL model in pure N2 gas in the tunnel s Mach 8 and Mach 10 nozzles at free stream Reynolds numbers of 4.1x10(exp 6)/ft to 49x10(exp 6)/ft (1.3x10(exp 7)/m to 16x10(exp 7)/m) and 1.2x10(exp 6)/ft to 19x10(exp 6)/ft (0.39x10(exp 7)/m to 62x10(exp 7)/m), respectively. These conditions were sufficient to span the regime of boundary-layer flow from completely laminar to fully-developed turbulent flow over the entire forebody. A supporting aeroheating test was also conducted in the Langley Research Center 20-Inch Mach 6 Air Tunnel at free stream Reynolds number of 1x10(exp 6)/ft to 7x10(exp 6)/ft (0.36x10(exp 7)/m to 2.2x10(exp 7)/m) in order to help corroborate the Tunnel 9 results. A complementary computational fluid dynamics study was conducted in parallel to the wind tunnel testing. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins on predictions for aeroheating environments during entry into the Martian atmosphere. Data from both wind tunnel tests and comparisons with the predictions are presented herein. It was concluded from these comparisons that for perfect-gas conditions, the computational tools could predict fully-laminar or fully-turbulent heating conditions to within 10% of the experimental data

Hollis, Brian R.↗

Heterogeneous fuel for hybrid rocket

Heterogeneous fuel compositions suitable for use in hybrid rocket engines and solid-fuel ramjet engines, The compositions include mixtures of a continuous phase, which forms a solid matrix, and a dispersed phase permanently distributed therein. The dispersed phase or the matrix vaporizes (or melts) and disperses into the gas flow much more rapidly than the other, creating depressions, voids and bumps within and on the surface of the remaining bulk material that continuously roughen its surface, This effect substantially enhances heat transfer from the combusting gas flow to the fuel surface, producing a correspondingly high burning rate, The dispersed phase may include solid particles, entrained liquid droplets, or gas-phase voids having dimensions roughly similar to the displacement scale height of the gas-flow boundary layer generated during combustion.

Stickler, David B.↗

Turbulent Aeroheating Testing of Mars Science Laboratory Entry Vehicle

An experimental investigation of turbulent aeroheating on the Mars Science Laboratory entry vehicle heat shield has been conducted in the Arnold Engineering Development Center Hypervelocity Wind Tunnel No. 9. Testing was performed on a 6-in. (0.1524 m) diameter MSL model in pure N2 gas in the tunnel's Mach 8 and Mach 10 nozzles at free stream Reynolds numbers of 4.1 x 10(exp 6)/ft to 49 x 10(exp 6)/ft (1.3 x 10(exp 7)/m to 19 x 10(exp 6/ft) and 1.2 x 10(exp 6)/ft to 19 x 10(exp 6)/ft (0.39 x 10(exp 7)/m to 62 x 10(exp 7)/m), respectively. These conditions were sufficient to span the regime of boundary-layer flow from completely laminar to fully-developed turbulent flow over the entire forebody. A supporting aeroheating test was also conducted in the Langley Research Center 20-Inch Mach 6 Air Tunnel at free stream Reynolds number of 1 x 10(exp 6)/ft to 7 x 10(exp 6)/ft (0.36 x 10(exp 7)/m to 2.2 x 10(exp 7)/m) in order to help corroborate the Tunnel 9 results. A complementary computational fluid dynamics study was conducted in parallel to the wind tunnel testing. Laminar and turbulent predictions were generated for the wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins on predictions for aeroheating environments during entry into the Martian atmosphere. Data from both wind tunnel tests and comparisons with the predictions are presented herein. It was concluded from these comparisons that for perfect-gas conditions, the computational tools could predict fully-laminar or fully-turbulent heating conditions to within 12% or better of the experimental data.

Hollis, Brian R.↗