Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 829 records · Page 46

Long, paired A'A/Pahoehoe flows of Mauna Loa: Volcanological significance and insights they provide into volcano plumbing systems

The long lava flows of Mauna Loa, Hawaii have been cited as Earth's closed analogs to the large Martian flows. It is therefore important to understand the flow mechanics and characteristics of the Mauna Loa flows and to make use of these in an attempt to gain insights into Martian eruptive processes. Two fundamentally different kinds of long lava flows can be distinguished on Hawaiian volcanoes as in Martian flows. The two kinds may have identical initial viscosities, chemical compositions, flow lengths, and flow volumes, but their flow mechanisms and thermal energy budgets are radically different. One travels a distance set by the discharge rate as envisaged by Walker and Wadge, and the other travels a distance set mainly by the eruption duration and ground slope. In the Mauna Loa lavas, yield strength becomes an important flow morphology control only in the distal part of a'a lavas. The occurrence of paired flows on Mauna Loa yields insights into the internal plumbing systems of the volcano, and it is significant that all of the volume of the a'a flow must be stored in a magma chamber before eruption, while none of the volume of the pahoehoe needs to be so stored. Differentiation between the two kinds of flows on images of Martian volcanoes is possible and hence an improved understanding of these huge structures is acquired.

Rowland, Scott K.↗

Navier-Stokes analysis and experimental data comparison of compressible flow within ducts

Many aircraft employ ducts with centerline curvature or changing cross-sectional shape to join the engine with inlet and exhaust components. S-ducts convey air to the engine compressor from the intake and often decelerate the flow to achieve an acceptable Mach number at the engine compressor by increasing the cross-sectional area downstream. Circular-to-rectangular transition ducts are used on aircraft with rectangular exhaust nozzles to connect the engine and nozzle. To achieve maximum engine performance, the ducts should minimize flow total pressure loss and total pressure distortion at the duct exit. Changes in the curvature of the duct centerline or the duct cross-sectional shape give rise to streamline curvature which causes cross stream pressure gradients. Secondary flows can be caused by deflection of the transverse vorticity component of the boundary layer. This vortex tilting results in counter-rotating vortices. Additionally, the adverse streamwise pressure gradient caused by increasing cross-sectional area can lead to flow separation. Vortex pairs have been observed in the exit planes of both duct types. These vortices are due to secondary flows induced by pressure gradients resulting from streamline curvature. Regions of low total pressure are produced when the vortices convect boundary layer fluid into the main flow. The purpose of the present study is to predict the measured flow field in a diffusing S-duct and a circular-to-rectangular transition duct with a full Navier-Stokes computer program, PARC3D, and to compare the numerical predictions with new detailed experimental measurements. The work was undertaken to extend previous studies and to provide additional CFD validation data needed to help model flows with strong secondary flow and boundary layer separation. The S-duct computation extends the study of Smith et al, and Harloff et al, which concluded that the computation might be improved by using a finer grid and more advanced turbulence models. The present study compares results for both the Baldwin-Lomas and k-epsilon turbulence models and is conducted with a refined grid. For the transition duct, two inlet conditions were considered, the first with straight flow and the second with swirling flow. The first case permits examination of the effects of the geometric transition on the flow field, while the second case includes the rotational flow effect characteristic of a gas turbine engine.

Harloff, G. J.↗

Large-volume lava flow fields on Venus: Dimensions and morphology

Of all the volcanic features identified in Magellan images, by far the most extensive and really important are lava flow fields. Neglecting the widespread lava plains themselves, practically every C1-MIDR produced so far contains several or many discrete lava flow fields. These range in size from a few hundred square kilometers in area (like those fields associated with small volcanic edifices for example), through all sizes up to several hundred thousand square kilometers in extent (such as many rift related fields). Most of these are related to small, intermediate, or large-scale volcanic edifices, coronae, arachnoids, calderas, fields of small shields, and rift zones. An initial survey of 40 well-defined flow fields with areas greater than 50,000 sq km (an arbitrary bound) has been undertaken. Following Columbia River Basalt terminology, these have been termed great flow fields. This represents a working set of flow fields, chosen to cover a variety of morphologies, sources, locations, and characteristics. The initial survey is intended to highlight representative flow fields, and does not represent a statistical set. For each flow field, the location, total area, flow length, flow widths, estimated flow thicknesses, estimated volumes, topographic slope, altitude, backscatter, emissivity, morphology, and source has been noted. The flow fields range from about 50,000 sq km to over 2,500,000 sq km in area, with most being several hundred square kilometers in extent. Flow lengths measure between 140 and 2840 km, with the majority of flows being several hundred kilometers long. A few basic morphological types have been identified.

Lancaster, M. G.↗

Control of flow separation in airfoil/wing design applications

Existing aerodynamic design methods have generally concentrated on the optimization of airfoil or wing shapes to produce a minimum drag while satisfying some basic constraints such as lift, pitching moment, or thickness. Since the minimization of drag almost always precludes the existence of separated flow, the evaluation and validation of these design methods for their robustness and accuracy when separated flow is present has not been aggressively pursued. However, two new applications for these design tools may be expected to include separated flow and the issues of aerodynamic design with this feature must be addressed. The first application of the aerodynamic design tools is the design of airfoils or wings to provide an optimal performance over a wide range of flight conditions (multipoint design). While the definition of 'optimal performance' in the multipoint setting is currently being hashed out, it is recognized that given a wide enough range of flight conditions, it will not be possible to ensure a minimum drag constraint at all conditions, and in fact some amount of separated flow (presumably small) may have to be allowed at the more demanding flight conditions. Thus a multipoint design method must be tolerant of the existence of separated flow and may include some controls upon its extent. The second application is in the design of wings with extended high speed buffet boundaries of their flight envelopes. Buffet occurs on a wing when regions of flow separation have grown to the extent that their time varying pressures induce possible destructive effects upon the wing structure or adversely effect either the aircraft controllability or the passenger comfort. A conservative approach to the expansion of the buffet flight boundary is to simply expand the flight envelope of nonseparated flow under the assumption that buffet will also thus be alleviated. However, having the ability to design a wing with separated flow and thus to control the location, extent, and severity of the separated flow regions may allow aircraft manufacturers to gain an advantage in the early design stages of an aircraft, when configuration changes are relatively inexpensive to make. Continuing the work begun last year, an airfoil design package has been modified to provide some control over the existence and extent of flow separation. This package consists of a 2-D Navier-Stokes flow solver which is coupled to the CDISC (constrained direct/iterative surface curvature) design method. The first modification is a prediction method for determining whether separation is likely based solely upon a given pressure distribution. If separation is predicted but is undesirable, the new routines will modify the pressure distribution to alleviate the problem. This new pressure distribution is then used in the design method to generate a new aerodynamic shape. Since separation may be acceptable in some cases, particularly if the separation does not extend to the trailing edge, another added logic estimates the extent of separation based upon a correlation with calculated separated flow cases. If the flow behind a shock induced separation is not predicted to reattach before the trailing edge, the logic weakens the shock strength and otherwise alters the pressure distribution in order to promote reattachment. This later addition is as yet unreliable due to secondary separation effects, but additional work is being pursued to improve the method.

Gally, Thomas A.↗

Critical Velocities in Open Capillary Flow

This paper describes the proposed research program on open capillary flow and the preliminary work performed theoretically and in drop tower experiments. The work focuses on the fundamental physical understanding of the flow through capillary bound geometries, where the circumference of the cross section of the flow path contains free surfaces. Examples for such a flow configuration are capillary vanes in surface tension tanks, flow along edges and corners and flow through liquid bridges. The geometries may be classified by their cross section areas, wetted circumferences and the radii of curvature of the free surfaces. In the streaming float zone the flow path is bound by a free surface only. The ribbon vane is a model for vane types used in surface tension tanks, where a structure in proximity to the tank wall forms a capillary gap. A groove is used in heat pipes for the transportation of the condensed working fluid to the heat source and a wedge may occur in a spaceborne experiment where fluid has to be transported by the means of surface tension. The research objectives are the determination of the maximum volume flux, the observation of the free surfaces and the liquid flow inside the flow path as well as the evaluation of the limiting capillary wave speed. The restriction of the maximum volume flux is due to convective forces (flow velocity exceeding the capillary wave speed) and/or viscous forces, i.e. the viscous head loss along the flow path must be compensated by the capillary pressure due to the curved free surface. Exceeding the maximum volume flux leads to the choking of the flow path, thus the free surface collapses and.gas ingestion occurs at the outlet. The means are ground-based experimental work with plateau tanks and in a drop tower, a sounding rocket flight, and theoretical analysis with integral balances as well as full three dimensional CFD solutions for flow with free surfaces.

Dreyer, Michael↗

The Dynamics of Rapidly Emplaced Terrestrial Lava Flows and Implications for Planetary Volcanism

The Kaupulehu 1800-1801 lava flow of Hualalai volcano and the 1823 Keaiwa flow from the Great Crack of the Kilauea southwest rift zone had certain unusual and possibly unique properties for terrestrial basaltic lava flows. Both flows apparently had very low viscosities, high effusion rates, and uncommonly rapid rates of advance. Ultramafic xenolith nodules in the 1801 flow form stacks of cobbles with lava rinds of only millimeter thicknesses. The velocity of the lava stream in the 1801 flow was extremely high, at least 10 m/s (more than 40 km/h). Observations and geological evidence suggest similarly high velocities for the 1823 flow. The unusual eruption conditions that produced these lava flows suggest a floodlike mode of emplacement unlike that of most other present-day flows. Although considerable effort has gone into understanding the viscous fluid dynamics and thermal processes that often occur in basaltic flows, the unusual conditions prevalent for the Kaupulehu and Keaiwa flows necessitate different modeling considerations. We propose an elementary flood model for this type of lava emplacement and show that it produces consistent agreement with the overall dimensions of the flow, channel sizes, and other supporting field evidence. The reconstructed dynamics of these rapidly emplaced terrestrial lava flows provide significant insights about the nature of these eruptions and their analogs in planetary volcanism.

Baloga, Stephen↗

Investigation of Flow Separation in a Transonic-fan Linear Cascade Using Visualization Methods

An extensive study into the nature of the separated flows on the suction side of modem transonic fan airfoils at high incidence is described in the paper. Suction surface.flow separation is an important flow characteristic that may significantly contribute to stall flutter in transonic fans. Flutter in axial turbomachines is a highly undesirable and dangerous self-excited mode of blade oscillations that can result in high cycle fatigue blade failure. The study basically focused on two visualization techniques: surface flow visualization using dye oils, and schlieren (and shadowgraph) flow visualization. The following key observations were made during the study. For subsonic inlet flow, the flow on the suction side of the blade is separated over a large portion of the blade, and the separated area increases with increasing inlet Mach number. For the supersonic inlet flow condition, the flow is attached from the leading edge up to the point where a bow shock from the upper neighboring blade hits the blade surface. Low cascade solidity, for the subsonic inlet flow, results in an increased area of separated flow. For supersonic flow conditions, a low solidity results in an improvement in flow over the suction surface. Finally, computational results modeling the transonic cascade flowfield illustrate our ability to simulate these flows numerically.

Lepicovsky, Jan↗

Flow Visualization Techniques in Wind Tunnel Tests of a Full-Scale F/A-18 Aircraft

The proposed paper presents flow visualization performed during experiments conducted on a full-scale F/A-18 aircraft in the 80- by 120-Foot Wind-Tunnel at NASA Ames Research Center. The purpose of the flow-visualization experiments was to document the forebody and leading edge extension (LEX) vortex interaction along with the wing flow patterns at high angles of attack and low speed high Reynolds number conditions. This investigation used surface pressures in addition to both surface and off-surface flow visualization techniques to examine the flow field on the forebody, canopy, LEXS, and wings. The various techniques used to visualize the flow field were fluorescent tufts, flow cones treated with reflective material, smoke in combination with a laser light sheet, and a video imaging system for three-dimension vortex tracking. The flow visualization experiments were conducted over an angle of attack range from 20 deg to 45 deg and over a sideslip range from -10 deg to 10 deg. The various visualization techniques as well as the pressure distributions were used to understand the flow field structure. The results show regions of attached and separated flow on the forebody, canopy, and wings as well as the vortical flow over the leading-edge extensions. This paper will also present flow visualization comparisons with the F-18 HARV flight vehicle and small-scale oil flows on the F-18.

Lanser, Wendy R.↗

Evaluation of Sulfur Flow Emplacement on Io from Galileo Data and Numerical Modeling

Galileo images of bright lava flows surrounding Emakong Patera have been analyzed and numerical modeling has been performed to assess whether these flows could have resulted from the emplacement of sulfur lavas on Io. Images from the solid-state imaging.(SSI) camera show that these bright, white to yellow Emakong flows are up to 370 km long and contain dark, sinuous features that are interpreted to be lava conduits, approx. 300-500 m wide and > 100 km long. Near-Infrared Mapping Spectrometer (NIMS) thermal emission data yield a color temperature estimate of 344 K +/- 60 K (less than or equal to 131 C) within the Emakong caldera. We suggest that these bright flows likely resulted from either sulfur lavas or silicate lavas that have undergone extensive cooling, pyroclastic mantling, and/or alteration with bright sulfurous materials. The Emakong bright flows have estimated volumes of approx. 250-350 cu km, similar to some of the smaller Columbia River Basalt flows. If the Emakong flows did result from effusive sulfur eruptions, then they are orders of magnitude greater in volume than any terrestrial sulfur flows. Our numerical modeling results show that sulfur lavas on Io could have been emplaced as turbulent flows, which were capable of traveling tens to hundreds of kilometers, consistent with the predictions of Sagan [ 19793 and Fink et al. [ 19831. Our modeled flow distances are also consistent with the measured lengths of the Emakong channels and bright flows. Modeled thermal erosion rates are approx. 1-4 m/d for flows erupted at approx. 140-180 C, which are consistent with the melting rates of Kieffer et al. [2000]. The Emakong channels could be thermal erosional in nature; however, the morphologic signatures of thermal erosion channels cannot be discerned from available images. There are planned Galileo flybys of Io in 2001 which provide excellent opportunities to obtain high-resolution morphologic and color data of Emakong Patera. Such observations could, along with further modeling, provide additional information to better constrain whether sulfur lavas produced the Emakong flows.

Williams, David A.↗

Laminar Flow About a Rotating Body of Revolution in an Axial Airstream

We have set ourselves the problem of calculating the laminar flow on a body of revolution in an axial flow which simultaneously rotates about its axis. The problem mentioned above, the flow about a rotating disk in a flow, which we solved some time ago, represents the first step in the calculation of the flow on the rotating body of revolution in a flow insofar as, in the case of a round nose, a small region about the front stagnation point of the body of revolution may be replaced by its tangential plane. In our problem regarding the rotating body of revolution in a flow, for laminar flow, one of the limiting cases is known: that of the body which is in an axial approach flow but does not rotate. The other limiting case, namely the flow in the neighborhood of a body which rotates but is not subjected to a flow is known only for the rotating circular cylinder, aside from the rotating disk. In the case of the cylinder one deals with a distribution of the circumferential velocity according to the law v = omega R(exp 2)/r where R signifies the cylinder radius, r the distance from the center, and omega the angular velocity of the rotation. The velocity distribution as it is produced here by the friction effect is therefore the same as in the neighborhood of a potential vortex. When we treat, in what follows, the general case of the rotating body of revolution in a flow according to the calculation methods of Prandtl's boundary-layer theory, we must keep in mind that this solution cannot contain the limiting case of the body of revolution which only rotates but is not subjected to a flow. However, this is no essential limitation since this case is not of particular importance for practical purposes.

Schlichting, H.↗

Investigation of Unsteady Flow Behavior in Transonic Compressor Rotors with LES and PIV Measurements

In the present study, unsteady flow behavior in a modern transonic axial compressor rotor is studied in detail with large eddy simulation (LES) and particle image velocimetry (PIV). The main purpose of the study is to advance the current understanding of the flow field near the blade tip in an axial transonic compressor rotor near the stall and peak-efficiency conditions. Flow interaction between the tip leakage vortex and the passage shock is inherently unsteady in a transonic compressor. Casing-mounted unsteady pressure transducers have been widely applied to investigate steady and unsteady flow behavior near the casing. Although many aspects of flow have been revealed, flow structures below the casing cannot be studied with casing-mounted pressure transducers. In the present study, unsteady velocity fields are measured with a PIV system and the measured unsteady flow fields are compared with LES simulations. The currently applied PIV measurements indicate that the flow near the tip region is not steady even at the design condition. This self-induced unsteadiness increases significantly as the compressor rotor operates near the stall condition. Measured data from PIV show that the tip clearance vortex oscillates substantially near stall. The calculated unsteady characteristics of the flow from LES agree well with the PIV measurements. Calculated unsteady flow fields show that the formation of the tip clearance vortex is intermittent and the concept of vortex breakdown from steady flow analysis does not seem to apply in the current flow field. Fluid with low momentum near the pressure side of the blade close to the leading edge periodically spills over into the adjacent blade passage. The present study indicates that stall inception is heavily dependent on unsteady behavior of the flow field near the leading edge of the blade tip section for the present transonic compressor rotor.

Hah, Chunill↗

Initial Observations of Lunar Impact Melts and Ejecta Flows with the Mini-RF Radar

The Mini-RF radar on the Lunar Reconnaissance Orbiter's spacecraft has revealed a great variety of crater ejecta flow and impact melt deposits, some of which were not observed in prior radar imaging. The craters Tycho and Glushko have long melt flows that exhibit variations in radar backscatter and circular polarization ratio along the flow. Comparison with optical imaging reveals that these changes are caused by features commonly seen in terrestrial lava flows, such as rafted plates, pressure ridges, and ponding. Small (less than 20 km) sized craters also show a large variety of features, including melt flows and ponds. Two craters have flow features that may be ejecta flows caused by entrained debris flowing across the surface rather than by melted rock. The circular polarization ratios (CPRs) of the impact melt flows are typically very high; even ponded areas have CPR values between 0.7-1.0. This high CPR suggests that deposits that appear smooth in optical imagery may be rough at centimeter- and decimeter- scales. In some places, ponds and flows are visible with no easily discernable source crater. These melt deposits may have come from oblique impacts that are capable of ejecting melted material farther downrange. They may also be associated with older, nearby craters that no longer have a radar-bright proximal ejecta blanket. The observed morphology of the lunar crater flows has implications for similar features observed on Venus. In particular, changes in backscatter along many of the ejecta flows are probably caused by features typical of lava flows.

Carter, Lynn M.↗

A regional comparison of sub-daily flow variability in regulated and unregulated rivers in the United States

Regulating rivers for hydropower or other purposes can dramatically alter river flow patterns, including creating substantial changes in flow over short, minutes-to-hours-long timespans known as sub-daily flow variability (SDFV). The impacts of flexible hydropower production on flow and aquatic organisms are increasingly documented in research. However, the degree to which flow alteration relates to different hydropower operational modes in distinct geographical regions and seasons is not well understood. This study offers a methodology for regional- and species-appropriate evaluations of potential impacts of flow on fish based on sub-daily flow characteristics of hydropower operational modes. We analyzed 15-min discharge data between 2018 and 2021 from 69 USGS stream gages to compare SDFV in hydropeaking, run-of-river, and unregulated systems in the US Southeast and Pacific Northwest. Regulated systems exhibited significant SDFV downstream from hydropower facilities relative to unregulated systems, but specific impacts differed between regions. Regulated systems in the Southeast were characterized by high flow coefficients of variation and ratios (hydropeaking only) and extended durations of daily upramping flow phases. Regulated systems in the Pacific Northwest were characterized by many short flow phases per day and large portions of the day spent upramping. Pacific Northwest unregulated systems displayed the strongest seasonal flow patterns while Southeastern hydropeaking systems displayed the greatest SDFV. Given that SDFV impacts multiple dimensions of fish ecology, region-specific sub-daily flow signatures have important implications for understanding and mitigating potential community-, species-, and age-specific effects on fish in different parts of the country.

Fish↗

Noncontact Flow Rate Using Laser Ultrasonics

Several types of advanced nuclear reactors are cooled with high-temperature liquid metal or molten salt flows. There is a critical need to measure flow velocity in flow channels for test purposes, and eventually in operational reactors. In conventional ultrasonic flow sensors, ultrasonic waves traveling in both the upstream and downstream directions are generated and detected by transducers that must contact the flow channels. A shift in the frequency or transit time between the two ultrasonic waves is measured to determine flow velocity. We describe here an initial effort to apply that sensing concept when the contact transducers are replaced by laser-based generation and detection instrumentation. This noncontact sensing avoids many practical problems associated with contact transducers when implemented on flow channels at high temperature. Laser-based flow monitoring can also be applied to hot-process piping in the geothermal energy, chemical-processing and petroleum-refining industries. Our effort has included theoretical simulation of noncontact laser-based flow monitoring, indicating capability of measuring flow velocities relevant to reactor cooling. It also included a room-temperature experimental demonstration using water as the flow liquid, and indicating capability of measuring flow velocity at a responsivity roughly consistent with simulation predictions. Plans were made for an experimental demonstration at Oregon State University using liquid metal at 110 degrees C.

02 PETROLEUM↗

Graphical and analytical methods for the determination of a flow of a compressible fluid around an obstacle

Chaplygin introduced the hodograph method in the theory of compressible fluid flows and developed a method for constructing stream functions of such flows. This method, which has been extensively used in investigation of compressible fluid flows, is limited in certain respects. The expression for the stream function obtained in this manner can represent only certain types of flow patterns. In general, flow patterns obtained in this way cannot represent the whole flow around an obstacle, but only a part of such a flow, and therefore several expressions are needed in order to obtain the whole flow. On the other hand, in many instances it is important to have a single expression representing the whole flow. Recently Von Karman and Tsien constructed more general types of stream functions, but only by replacing the true pressure density relation by the linear pressure-specific volume relation so that their method is essential limited to flows the maximum Mach number of which is not too large. In a companion report the author derived a new formula for stream functions based on the true pressure density relation. It is not subject to the limitations of the Chaplygin method. In the present report this formula is employed to construct two-dimensional subsonic compressible fluid flows around a body similar in shape to a given symmetric obstacle. The methods described in the report are illustrated by numerical examples.

FLOW, COMPRESSIBLE-MEASUREMENT↗

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗

Efficient Thermal Energy Storage with Radial Flow in Packed Beds

Thermal energy storage (TES) is a way to store excess heat in order to generate power later. These types of systems are valuable in the solar industry, where power production can still occur when sunlight is not available. This not only increases system flexibility, it drives down the cost of electricity. One option to store thermal energy is with a packed bed where the storage media resides in a cylindrical container. Flow from one end of the cylinder to the other deposits or recovers heat (axial flow). While a promising technology in terms of energy storage, it exhibits a high pressure drop that lowers overall system efficiency. This project seeks to flow the heat transfer fluid through the storage media in the radial direction. This system offers the potential to retain reasonable thermal performance while substantially lowering the pressure drop. The overall goal was to assess the performance of radial flow experimentally and numerically. Three main designs were considered in this work. The first system utilized a central pipe and four receiving pipes near the wall. The second system utilized a central pipe and an annulus with holes near the wall. The last approach considered segmenting the bed so only select zones would receive flow in the radial direction. In a high-aspect ratio system, the hole patterns in the piping are important for ensuring even flow into the bed, and a major design effort was testing holes patterns to promote even flow. The four pipe system was shown to be not feasible, as preferential flow paths occur that do not allow even and full heating of the packed bed. The annular system can lead to even flow, but the results show this is better accomplished through larger piping instead of considering variations in the hole pattern along the axial length of the pipe. Segmenting the bed can lead to similar exergetic performance when compared to axial or radial flow. However, the pressure drop in the segmented design is likely too high for practical implementation. A baseline commercial design showed exergetic efficiencies for axial, radial, and radial segments of 81.6%, 82.8%, and 80.2%, respectively. Pressure drop for the axial and radial results were 2.36 psi and 2.45 psi, respectively, with segments being nearly an order of magnitude higher. At these large scales, the aspect ratio of the system is important. From a baseline of 0.64, an aspect ratio of 0.32 for radial flow showed an exergetic efficiency of 86% but a pressure drop of 5.39 psi. All of these results provide new insights into packed bed thermal energy storage with radial flow. Competing effects must be considered when designing a radial system, and results show a radial design can show strong thermal results at the expense of system efficiencies from pressure drop. The COMSOL models used in this analysis are available to the public and other researchers. Public benefits include a potential increase in the efficiency of packed bed thermal energy storage. Higher efficiency in storage promotes green energy technologies while reducing costs.

14 SOLAR ENERGY↗

Shear rate dependency on flowing granular biomass material

The commercialization of bioenergy has been significantly limited by various material handling issues due to the poor flowability of granular biomass materials. A good understanding of flow physics and robust constitutive models to predict flow behavior across multiple regimes are essential to address these issues. In this study, we investigated the multi-regime flow behavior of loblolly pine chips, a widely used bioenergy feedstock, through comprehensive inclined plane flow experiments and simulations. A quasi-static hypoplastic model and a cross-regime Drucker-Prager-µ(I) model were calibrated and validated against the physical experiments to investigate the quasi-static and dense flow behavior. The results show that for granular biomass, 1) plane flow (iso-thickness along the plane) exists within a smaller range of inclination while heap flow (varying thickness along the plane) exists within a broader range of inclination, as compared with conventional granular materials (e.g., glass beads); 2) the scaling law of granular biomass flowing on an inclined plane (Froude Number versus dimensionless thickness) forms a bi-linear trend with the turning point governed by the quasi-static and dense flow regimes; 3) the multi-regime DP-µ(I) model can capture the flow behavior in both regimes well at the cost of extra calibration. In conclusion, these findings advance the scientific understanding of the multi-regime flow behavior of granular biomass materials and shed light on formulating novel constitutive models to assist granular biomass handling in the bioenergy industry.

09 BIOMASS FUELS↗