Search NASA⌕ Search

SEARCH · Search NASA

Results for “evolving boundary layer”

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 163 records · Page 9

Coherent structures in a boundary layer and shear layer of a turbulent backward-facing step flow

A wind tunnel experiment has been carried out at the NASA Ames Research Center to analyze the evolution of coherent structures from a boundary layer to a shear layer in a turbulent, backward-facing, step flow. A miniature X-wire/cold-wire probe has been used in conjunction with two arrays of cold wires, one aligned in the plane of main shear and the other in the spanwise direction of the flow, to detect and characterize delta-scale organized structures in the outer regions of the flow and to provide detailed information concerning these structures. Kinematic features of the events associated with the large scale structures were analyzed and topological pictures of the evolving flow, as well as the contributions to the Reynolds shear stress components are presented.

Jovic, Srba↗

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗

Reconstruction framework advancements to support streaming for the ePIC detector at the EIC

The ePIC collaboration adopted the JANA2 framework to manage its reconstruction algorithms. This framework has since evolved substantially in response to ePIC’s needs. There have been three main design drivers: integrating cleanly with the Podio-based data models and other layers of the key4hep stack, enabling external configuration of existing components, and supporting timeframe splitting for streaming readout. The result is a unified component model featuring a new declarative interface for specifying inputs, outputs, parameters, services, and resources. This interface enables the user to instantiate, configure, and wire components via an external file. One critical new addition to the component model is a hierarchical decomposition of data boundaries into levels such as Run, Timeframe, PhysicsEvent, and Subevent. Two new component abstractions, Folder and Unfolder, are introduced in order to traverse this hierarchy, e.g. by splitting or merging. The pre-existing components can now operate at different event levels, and JANA2 will automatically construct the corresponding parallel processing topology. This means that a user may write an algorithm once, and configure it at runtime to operate on timeframes or on physics events. Overall, these changes mean that the user requires less knowledge about the framework internals, obtains greater flexibility with configuration, and gains the ability to reuse the existing abstractions in new streaming contexts.

Brei, Nathan [Thomas Jefferson National Accelerato↗

Evolution of giant gaseous protoplanets embedded in the primitive solar nebula

In contrast to earlier evolutionary calculations, which assumed that protoplanets are isolated, the present spherically symmetric approximation for a protoplanet of one Jovian mass during the early phase of quasi-static contraction invokes a time-dependent surface boundary condition that simulates physical conditions in an evolving primitive solar nebula. Assuming in a first set of calculations that the protoplanet is surrounded by a thermal bath whose temperature varies with time, and whose pressure is small and constant, results show evaporation and complete dispersal of the object. A second set of calculations varies both temperature and pressure at the surface with time, according to solar nebula models, and yields an acceleration or retardation of evolution, relative to that of an isolated protoplanet, depending on the relative entropy of the nebula and the protoplanet's outer layers. Processes by which terrestrial planets can form in the cores of giant gaseous protoplanets are discussed.

Cameron, A. G. W.↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

Simulations of curved turbulent boundary layers

The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Orlandi, Paolo↗

The development of concentration gradients in a suspension of chemotactic bacteria

When a suspension of bacterial cells of the species Bacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e., they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development when N and C depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

NASA Discipline Cell Biology↗

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

An evolving MHD vortex street model for quasi-periodic solar wind fluctuations

Magnetohydrodynamic (MHD) simulation are used to provide a dynamical basis for the 'vortex street' model of the quasi-periodic meridional flow observed by Voyager 2 in the outer heliosphere. Various observations suggest the existence near the current sheet at solar minimum, of a vorticity distribution of two opposite shear layers with an antisymmetric staggered velocity pattern due to structured high-speed wind surrounding low-speed equatorial flow. It is shown that this flow pattern leads to the formation of a highly stable vortex street through the nonlinear interaction of the two shear layers. Spatial profiles of various simulated parameters (velocity, density, meridional flow angle and the location of magnetic sector boundaries) and their relative locations in the quasi-steady vortex street are generally in good agreement with the observations.

Siregar, Edouard↗

Toward a Turbulence Constitutive Relation for Rotating Flows

In rapidly rotating turbulent flows the largest scales of the motion are in approximate geostrophic balance. Single-point turbulence closures, in general, cannot attain a geostrophic balance. This article addresses and resolves the possibility of constitutive relation procedures for single-point second order closures for a specific class of rotating or stratified flows. Physical situations in which the geostrophic balance is attained are described. Closely related issues of frame-indifference, horizontal nondivergence, Taylor-Proudman theorem and two-dimensionality are, in the context of both the instantaneous and averaged equations, discussed. It is shown, in the absence of vortex stretching along the axis of rotation, that turbulence is frame-indifferent. A derivation and discussion of a geostrophic constraint which the prognostic equations for second-order statistics must satisfy for turbulence approaching a frame-indifferent limit is given. These flow situations, which include rotating and nonrotating stratified flows, are slowly evolving flows in which the constitutive relation procedures are useful. A nonlinear non-constant coefficient representation for the rapid-pressure strain covariance appearing in the Reynolds stress and heat flux equations consistent with the geostrophic balance is described. The rapid-pressure strain model coefficients are not constants determined by numerical optimization but are functions of the state of the turbulence as parameterized by the Reynolds stresses and the turbulent heat fluxes. The functions are valid for all states of the turbulence attaining their limiting values only when a limit state is achieved. These issues are relevant to strongly vortical flows as well as flows such as the planetary boundary layers, in which there is a transition from a three-dimensional shear driven turbulence to a geostrophic or horizontal turbulence.

Ristorcelli, J. R.↗

Three-dimensional compressible convection at low Prandtl numbers

Numerical simulations are used to study fully compressible thermal convection at large Rayleigh numbers. Results are presented from a sequence of three-dimensional simulations that reveal a transition from gradually-evolving laminar convection to nearly turbulent convection as the Prandtl number is reduced from a value of unity to one-tenth. The convective flows form irregular cellular patterns near the upper surface, possesing a network of fast downflow at cell peripheries and gentler upflow at cell centers. At greater depths the curving sheets of downflow collapse into plumes which may twist and possess substantial vertical vorticity. For the lowest Prandtl number, the convection near the bottom of the layer appears to be turbulent, yet the rapidly varying small-scale flow structure there is accompanied by more ordered sites of wavering upflow, with the latter able to penetrate all the way to the upper boundary. Thus a significant component of the flow is able to extend over multiple density scale heights, in contrast to what is argued in formulating mixing-length models for stellar convection. Results are also shown from two-dimensional simulations carried out with very high spatial resolution, which reveal that supersonic convection with fluttering shock systems can be realized.

Toomre, Juri↗

Development of buried wire gages for measurement of wall shear stress in Blastane experiments

Buried Wire Gages operated from a Constant Temperature Anemometer System are among the special types of instrumentation to be used in the Boundary Layer Apparatus for Subsonic and Transonic flow Affected by Noise Environment (BLASTANE). These Gages are of a new type and need to be adapted for specific applications. Methods were developed to fabricate Gage inserts and mount those in the BLASTANE Instrumentation Plugs. A large number of Gages were prepared and operated from a Constant Temperature Anemometer System to derive some of the calibration constants for application to fluid-flow wall shear-stress measurements. The final stage of the calibration was defined, but could not be accomplished because of non-availability of a suitable flow simulating apparatus. This report provides a description of the Buried Wire Gage technique, an explanation of the method evolved for making proper Gages and the calibration constants, namely Temperature Coefficient of Resistance and Conduction Loss Factor.

Murthy, S. V.↗

Development and calibration of buried wire gages for wall shear stress measurements in fluid flow

Special methods were developed to arrange 'Buried Wire Gage' inserts flush to the contoured flow surfaces of instrument plugs of a boundary-layer flow apparatus. The fabrication process was aimed at producing proper bonding of the sensor wire to the substrate surface, without causing excessive surface waviness. A large number of gages were built and first calibrated for the resistance-temperature characteristics. The gages were then installed in a flow calibration apparatus and operated from a constant temperature anemometer system for a series of flow settings to derive the calibration constants of each of the gages. The flow settings included a range of subsonic freestream Mach numbers in order to help establish the gage calibration characteristics for compressible flow fields. This paper provides a description of the buried wire gage technique, an explanation of the method evolved for making proper gages, the procedure for calibrating the gages and the results of measurements performed for determining the calibration constants.

Murthy, Sreedhara V.↗

Farming in space: environmental and biophysical concerns

The colonization of space will depend on our ability to routinely provide for the metabolic needs (oxygen, water, and food) of a crew with minimal re-supply from Earth. On Earth, these functions are facilitated by the cultivation of plant crops, thus it is important to develop plant-based food production systems to sustain the presence of mankind in space. Farming practices on earth have evolved for thousands of years to meet both the demands of an ever-increasing population and the availability of scarce resources, and now these practices must adapt to accommodate the effects of global warming. Similar challenges are expected when earth-based agricultural practices are adapted for space-based agriculture. A key variable in space is gravity; planets (e.g. Mars, 1/3 g) and moons (e.g. Earth's moon, 1/6 g) differ from spacecraft orbiting the Earth (e.g. Space stations) or orbital transfer vehicles that are subject to microgravity. The movement of heat, water vapor, CO2 and O2 between plant surfaces and their environment is also affected by gravity. In microgravity, these processes may also be affected by reduced mass transport and thicker boundary layers around plant organs caused by the absence of buoyancy dependent convective transport. Future space farmers will have to adapt their practices to accommodate microgravity, high and low extremes in ambient temperatures, reduced atmospheric pressures, atmospheres containing high volatile organic carbon contents, and elevated to super-elevated CO2 concentrations. Farming in space must also be carried out within power-, volume-, and mass-limited life support systems and must share resources with manned crews. Improved lighting and sensor technologies will have to be developed and tested for use in space. These developments should also help make crop production in terrestrial controlled environments (plant growth chambers and greenhouses) more efficient and, therefore, make these alternative agricultural systems more economically feasible food production systems. c2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Review↗

Genesis of Twin Tropical Cyclones as Revealed by a Global Mesoscale Model: The Role of Mixed Rossby Gravity Waves

In this study, it is proposed that twin tropical cyclones (TCs), Kesiny and 01A, in May 2002 formed in association with the scale interactions of three gyres that appeared as a convectively coupled mixed Rossby gravity (ccMRG) wave during an active phase of the Madden-Julian Oscillation (MJO). This is shown by analyzing observational data, including NCEP reanalysis data and METEOSAT 7 IR satellite imagery, and performing numerical simulations using a global mesoscale model. A 10-day control run is initialized at 0000 UTC 1 May 2002 with grid-scale condensation but no sub-grid cumulus parameterizations. The ccMRG wave was identified as encompassing two developing and one non-developing gyres, the first two of which intensified and evolved into the twin TCs. The control run is able to reproduce the evolution of the ccMRG wave and thus the formation of the twin TCs about two and five days in advance as well as their subsequent intensity evolution and movement within an 8-10 day period. Five additional 10-day sensitivity experiments with different model configurations are conducted to help understand the interaction of the three gyres, leading to the formation of the TCs. These experiments suggest the improved lead time in the control run may be attributed to the realistic simulation of the ccMRG wave with the following processes: (1) wave deepening (intensification) associated with a reduction in wavelength and/or the intensification of individual gyres, (2) poleward movement of gyres that may be associated with boundary layer processes, (3) realistic simulation of moist processes at regional scales in association with each of the gyres, and (4) the vertical phasing of low- and mid-level cyclonic circulations associated with a specific gyre.

forecast↗

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗