Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow forming”

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 451 records · Page 25

Formation of perched lava ponds on basaltic volcanoes: Interaction between cooling rate and flow geometry allows estimation of lava effusion rates

Perched lava ponds are infrequent but distinctive topographic features formed during some basaltic eruptions. Two such ponds, each approximately 150 m in diameter, formed during the 1968 eruption at Napau Crater and the 1974 eruption of Mauna Ulu, both on Kilauea Volcano, Hawaii. Each one formed where a channelized, high volume flux lava flow encountered a sharp reduction of slope: the flow spread out radially and stalled, forming a well-defined terminal levee enclosing a nearly circular lava pond. We describe a model of how cooling limits the motion of lava spreading radially into a pond and compare this with the case of a channelized flow. The difference in geometry has a major effect, such that the size of a pond is a good indicator of the volume flux of the lava forming it. Lateral spreading on distal shallow slopes is a major factor limiting the lengths of lava flows.

Wilson, L.↗

The response of a laminar boundary layer in supersonic flow to small amplitude progressive waves

The effect of a small amplitude progressive wave on the laminar boundary layer on a semi-infinite flat plate, due to a uniform supersonic freestream flow, is considered. The perturbation to the flow divides into two streamwise zones. In the first, relatively close to the leading edge of the plate, on a transverse scale comparable to the boundary layer thickness, the perturbation flow is described by a form of the unsteady linearized compressible boundary layer equations. In the freestream, this component of flow is governed by the wave equation, the solution of which provides the outer velocity conditions for the boundary layer. This system is solved numerically, and also the asymptotic structure in the far downstream limit is studied. This reveals a breakdown and a subsequent second streamwise zone, where the flow disturbance is predominantly inviscid. The two zones are shown to match in a proper asymptotic sense.

Duck, Peter W.↗

Analysis of the fluid flow and heat transfer in a thin liquid film in the presence and absence of gravity

The hydrodynamic and thermal behavior of a thin liquid film flowing over a solid horizontal surface is analyzed for both plane and radially spreading flows. The situations where the gravitational force is completely absent and where it is significant are analyzed separately and their practical relevance to a micro-gravity environment is discussed. In the presence of gravity, in addition to Reynolds number, the Froude number of the film is found to be an important parameter that determines the supercritical and subcritical flow regimes and any associated hydraulic jump. A closed-form solution is possible under some flow situations, whereas others require numerical integration of ordinary differential equations. The approximate analytical results are found to compare well with the available two-dimensional numerical solutions.

Rahman, M. M.↗

On steady field-aligned double-adiabatic flow

The paper summarizes the four general double-adiabatic invariants governing ideal MHD flow in a narrow tube and renders them in a form suitable for the calculation of flow in tubes of nonconstant cross-sectional area. Results for converging-diverging flux tubes are developed and compared to single-adiabatic results which are identical to those for ordinary gasdynamic flow. It is found that in certain parameter regimes the governing equation for the vector potential is elliptic and can be reduced to Laplace's equation by simple stretching of one coordinate. In other regimes the equation is hyperbolic, indicating the presence of wave patterns. The relevance of these results to flow over magnetopause undulations, including flux transfer events, is discussed.

Sonnerup, B. U. O.↗

Energy Equation Approximation in Fluid Mechanics

There is some confusion in the literature of fluid mechanics in regard to the correct form of the energy equation for the study of the flow of nearly incompressible fluids. Several forms of the energy equation and their use are therefore discussed in this note.

Goldstein, Arthur W.↗

Parametric study of grid size, time step and turbulence modeling on Navier-Stokes computations over airfoils

An upwind-biased implicit approximate factorization algorithm is applied to several steady and unsteady turbulent flows. The thin layer form of the compressible Navier-Stokes equation is used. Both the flux vector splitting and flux difference splitting methods are used to determine fluxes, and the results are compared. Flux difference splitting predicts results more accurately than flux vector splitting on a given mesh size, but, in its present implementation, is more severely limited by the maximum CFL number for unsteady time accurate flows. Physical aspects of the computations are also examined. An equilibrium turbulent boundary layer model computes generally better steady and unsteady results than a nonequilibrium model when there is little to no boundary layer separation. Conversely, when a significant region of separation exists, the nonequilibrium model performs in better agreement with experiment.

Rumsey, Christopher L.↗

Shaped Recess Flow Control

One or more embodiments of techniques or systems for shaped recess flow control are provided herein. A shaped recess or cavity can be formed on a surface associated with fluid flow. The shaped recess can be configured to create or induce fluid effects, temperature effects, or shedding effects that interact with a free stream or other structures. The shaped recess can be formed at an angle to a free stream flow and may be substantially "V" shaped. The shaped recess can be coupled with a cooling channel, for example. The shaped recess can be upstream or downstream from a cooling channel and aligned in a variety of manners. Due to the fluid effects, shedding effects, and temperature effects created by a shaped recess, lift-off or separation of cooling jets of cooling channels can be mitigated, thereby enhancing film cooling effectiveness.

Shyam, Vikram↗

Steady vibrations of wing of circular plan form

This paper treats the problem of determining the lift, moment, and induced drag of a thin wing of circular plan form in uniform incompressible flow on the basis of linearized theory. As contrasted to a similar paper by Kinner, in which the acceleration potential method was used, the present paper utilizes the concept of the velocity potential. Calculations of the lift and moment are presented for several deformed shapes. It is shown that considerable deviations exist between the strip theory analysis and the more exact theory. The lift, moment, and induced drag are also determined for a harmonically oscillatory circular plan form wing. As contrasted to a similar paper by Schade, in which the acceleration potential method was used, the present paper utilizes the concept of the velocity potential. Expressions for lift, moment, and induced drag are given and finally specialized to the case of a slowly oscillating circular wing.

FLUTTER - WINGS↗

A theory for turbulent curved wall jets

An attempt is made to formulate the simplest possible model for the flow of a wall jet emanating from a two dimensional source into quiescent surroundings in the presence of a wall of arbitrary shape. The method used self similar profiles for the mean velocity together with a simple eddy viscosity model. The streamwise and radial momentum equations are integrated across the wall jet flow to give an expression for the momentum balance including the effect of the pressure gradient induced by the rate of change of surface curvature. The streamwise momentum equation is also evaluated at the point of maximum velocity to provide a second equation and thereby permit a solution for the two unknown quantities, the jet half width and the maximum velocity. This approach provides approximate closed form solutions for the flow of the wall jet over surfaces of various shapes, and in particular permits a direct comparison with the available experimental results for plane, cylindrical or logarithmic spiral surfaces.

Roberts, Leonard↗

Active Flow Control for Trailing Edge Flap Separation

Active flow control (AFC) in the form of sweeping jet (SWJ) excitation and discrete steady jet excitation is used to control the flow separation on an NACA 0015 semispan wing with a deflected, simple-hinged, trailing edge flap. This geometry has been the focus of several recent publications that investigated methods to improve the efficiency of sweeping jet actuators. In the current study, the interaction of the AFC excitation with the separated flowfields present at several flap deflection angles was examined. Previous studies with this model have been limited to a maximum flap deflection angle of 40°. The flap deflection range was extended to 60! because systems studies have indicated that a high-lift system with simple-hinged flaps may require larger flap deflections than the Fowler flaps found on most high-lift systems. The results obtained at flap deflection angles of 20°, 40°, and 60° are presented and compared. Force and moment data, Particle Image Velocimetry (PIV) data, and steady and unsteady surface pressure data are used to describe the flowfield with and without AFC. With a flap deflection of 60°, increasing the SWJ actuator momentum at the flap shoulder increased lift due to an increase in circulation but did not completely eliminate the recirculation region above the flap surface. AFC using the discrete steady jet actuators of this study increased lift as well but required more mass flow than the SWJ actuators and had a detrimental effect on lift at the highest mass flow level tested. PIV results showed that the angle between the excitation and the flap surface was not optimal for attaching the separated shear layer.

Melton, LaTunia Pack↗

Concepts for radically increasing the numerical convergence rate of the Euler equations

Integral equation and finite difference methods have been developed for solving transonic flow problems using linearized forms of the transonic small disturbance and Euler equations. A key element is the use of a strained coordinate system in which the shock remains fixed. Additional criteria are developed to determine the free parameters in the coordinate straining; these free parameters are functions of the shock location. An integral equation analysis showed that the shock is located by ensuring that no expansion shocks exist in the solution. The expansion shock appears as oscillations in the solution near the sonic line, and the correct shock location is determined by removing these oscillations. A second objective was to study the ability of the Euler equation to model separated flow.

Nixon, David↗

Intrepretation of surface features and surface processes on Mars

Eolian erosion and deposition on earth was studied in order to interpret the eolian land forms of Mars. Emphasis of the wind tunnel studies was on the flow field around models of eolian forms. Areas of the wind tunnel studies include: simulation of the atmospheric boundary layer; velocity profile measurements around different models in the desert boundary layer, and estimating shear stress distributions on the model surfaces; flow visualization techniques; streamline mapping using tuft photographs; and roughness contrasts.

Gad-El-hak, M.↗

Numerical simulation of flows in curved diffusers with cross-sectional transitioning using a three-dimensional viscous analysis

A three dimensional analysis for fully viscous, subsonic, compressible flow is evaluated. An approximate form of the Navier Stokes equations is solved by an implicit spatial marching technique. Calculations were made for flow in a circular S duct and in the F 16 inlet duct. The computed total pressure contours and secondary flow velocity vectors are presented. Qualitative comparisons with experiment are shown for both ducts. The analysis is used to show how the cross section transitioning in the F 16 inlet suppresses the development of a secondary flow vortex.

Towne, C. E.↗

Generation of Martian chaos and channels by debris flows

A debris flow mechanism is proposed to account for the formation of chaos and the large channels debouching into Crysae Planitia from the adjacent southern uplands of Mars. Based on considerations of the juxtaposition of individual channel environments, the morphological assemblages within each environment and flow dynamics, it is suggested that the debris flows were triggered by the large-scale failure of subsurface sediments, possibly initiated by a seismic event. During the initial, slow-moving phase of the flow, the debris would have formed gently sinuous channels with multiple side-wall slumps, grooves and ridges, and elongate erosional remnants. The flow would have gained mobility as the debris moved downslope, producing travel distances greatly in excess of those characteristic of terrestrial examples, and eroded, streamlined remnants at the distal reaches of the channel. Finally, due to internal and boundary friction, the flow would have been slowed down once it entered the Chryse plains, resulting in a thin debris blanket with no depositional relief.

Nummedal, D.↗

Systems Modeling to Implement Integrated System Health Management Capability

ISHM capability includes: detection of anomalies, diagnosis of causes of anomalies, prediction of future anomalies, and user interfaces that enable integrated awareness (past, present, and future) by users. This is achieved by focused management of data, information and knowledge (DIaK) that will likely be distributed across networks. Management of DIaK implies storage, sharing (timely availability), maintaining, evolving, and processing. Processing of DIaK encapsulates strategies, methodologies, algorithms, etc. focused on achieving high ISHM Functional Capability Level (FCL). High FCL means a high degree of success in detecting anomalies, diagnosing causes, predicting future anomalies, and enabling health integrated awareness by the user. A model that enables ISHM capability, and hence, DIaK management, is denominated the ISHM Model of the System (IMS). We describe aspects of the IMS that focus on processing of DIaK. Strategies, methodologies, and algorithms require proper context. We describe an approach to define and use contexts, implementation in an object-oriented software environment (G2), and validation using actual test data from a methane thruster test program at NASA SSC. Context is linked to existence of relationships among elements of a system. For example, the context to use a strategy to detect leak is to identify closed subsystems (e.g. bounded by closed valves and by tanks) that include pressure sensors, and check if the pressure is changing. We call these subsystems Pressurizable Subsystems. If pressure changes are detected, then all members of the closed subsystem become suspect of leakage. In this case, the context is defined by identifying a subsystem that is suitable for applying a strategy. Contexts are defined in many ways. Often, a context is defined by relationships of function (e.g. liquid flow, maintaining pressure, etc.), form (e.g. part of the same component, connected to other components, etc.), or space (e.g. physically close, touching the same common element, etc.). The context might be defined dynamically (if conditions for the context appear and disappear dynamically) or statically. Although this approach is akin to case-based reasoning, we are implementing it using a software environment that embodies tools to define and manage relationships (of any nature) among objects in a very intuitive manner. Context for higher level inferences (that use detected anomalies or events), primarily for diagnosis and prognosis, are related to causal relationships. This is useful to develop root-cause analysis trees showing an event linked to its possible causes and effects. The innovation pertaining to RCA trees encompasses use of previously defined subsystems as well as individual elements in the tree. This approach allows more powerful implementations of RCA capability in object-oriented environments. For example, if a pressurizable subsystem is leaking, its root-cause representation within an RCA tree will show that the cause is that all elements of that subsystem are suspect of leak. Such a tree would apply to all instances of leak-events detected and all elements in all pressurizable subsystems in the system. Example subsystems in our environment to build IMS include: Pressurizable Subsystem, Fluid-Fill Subsystem, Flow-Thru-Valve Subsystem, and Fluid Supply Subsystem. The software environment for IMS is designed to potentially allow definition of any relationship suitable to create a context to achieve ISHM capability.

Figueroa, Jorge F.↗

Thermodynamic constraints on the textural evolution of eucrite, EET 90020

Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.

Eucrites↗

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analyses were performed to identify the mechanisms by whichtwo-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for threeboundary-layer transition trip configurations:a two-dimensional smooth hump,a two-dimensional rectangular-shapedrod and a three-dimensional zig-zagtape. The smooth hump was used by Park et al.1in the linear and nonlinear parabolized stability (PSE) computations. Reynolds numbersbased on the tripheights are 710, 683,and 300 for the two-and three-dimensional trips, respectively. Simulations withthe smooth hump showed that the flow becomes stable upstream and strongly unstable downstream of the hump. Simulations with the two-dimensional rectangular trip showed that the flow separates at the top of the trip and forms a long recirculation zonedownstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff2 experiment. Three-dimensional simulations with secondary disturbances of fundamental, subharmonic, and oblique types produced turbulent flows shortlydownstream. The three-dimensional zig-zag trip with a small roughness Reynolds number lead to early breakdown to a turbulent flow.

Stability↗

Experimental and numerical investigation of supersonic turbulent flow in an annular duct

Experimental and numerical results are presented for developing supersonic turbulent flow in an annular duct formed by a circular centerbody and outer shroud. The experimental results are based on data taken in a new flow facility that was designed to generate a shock-free, supersonic annular flow. Numerical computations were performed using the Baldwin-Lomax turbulence model for comparison with experimentally measured profiles. The results demonstrate that computed and measured profiles are in excellent agreement, so that studies can now be conducted of shock wave/boundary layer interaction phenomena within the duct, such as those induced by changes in downstream duct geometry or by the placement of struts between the duct walls.

Williams, K. E.↗