Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solutal buoyancy”

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.

229 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Generation of Bubbly Suspensions in Low Gravity

Generation of a uniform monodisperse bubbly suspension in low gravity is a rather difficult task because bubbles do not detach as easily as on Earth. Under microgravity, the buoyancy force is not present to detach the bubbles as they are formed from the nozzles. One way to detach the bubbles is to establish a detaching force that helps their detachment from the orifice. The drag force, established by flowing a liquid in a cross or co-flow configuration with respect to the nozzle direction, provides this additional force and helps detach the bubbles as they are being formed. This paper is concerned with studying the generation of a bubbly suspension in low gravity in support of a flight definition experiment titled "Behavior of Rapidly Sheared Bubbly Suspension." Generation of a bubbly suspension, composed of 2 and 3 mm diameter bubbles with a standard deviation <10% of the bubble diameter, was identified as one of the most important engineering/science issues associated with the flight definition experiment. This paper summarizes the low gravity experiments that were conducted to explore various ways of making the suspension. Two approaches were investigated. The first was to generate the suspension via a chemical reaction between the continuous and dispersed phases using effervescent material, whereas the second considered the direct injection of air into the continuous phase. The results showed that the reaction method did not produce the desired bubble size distribution compared to the direct injection of bubbles. However, direct injection of air into the continuous phase (aqueous salt solution) resulted in uniform bubble-diameter distribution with acceptable bubble-diameter standard deviation.

Nahra, Henry K.↗
NASA NTRS1999

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗
NASA NTRS2003

Simulating the Effect of Space Vehicle Environments on Directional Solidification of a Binary Alloy

Space microgravity missions are designed to provide a microgravity environment for scientific experiments, but these missions cannot provide a perfect environment, due to vibrations caused by crew activity, on-board experiments, support systems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps. Therefore, it is necessary to predict the impact of these vibrations on space experiments, prior to performing them. Simulations were conducted to study the effect of the vibrations on the directional solidification of a dendritic alloy. Finite element ca!cu!attie?ls were dme with a simd2titcr based on a continuum model of dendritic solidification, using the Fractional Step Method (FSM). The FSM splits the solution of the momentum equation into two steps: the viscous intermediate step, which does not enforce continuity; and the inviscid projection step, which calculates the pressure and enforces continuity. The FSM provides significant computational benefits for predicting flows in a directionally solidified alloy, compared to other methods presently employed, because of the efficiency gains in the uncoupled solution of velocity and pressure. finite differences, arises when the interdendritic liquid reaches the eutectic temperature and concentration. When a node reaches eutectic temperature, it is assumed that the solidification of the eutectic liquid continues at constant temperature until all the eutectic is solidified. With this approach, solidification is not achieved continuously across an element; rather, the element is not considered solidified until the eutectic isotherm overtakes the top nodes. For microgravity simulations, where the convection is driven by shrinkage, it introduces large variations in the fluid velocity. When the eutectic isotherm reaches a node, all the eutectic must be solidified in a short period, causing an abrupt increase in velocity. To overcome this difficulty, we employed a scheme to numerically predict a more accurate value for the rate of change of fraction of liquid as the liquid in an element solidifies. The new method enables us to contrast results of simulations in which the alloy is subjected to no gravity or a steady-state acceleration versus simulations when the alloy is subjected to vibration disturbances; therefore, the effect of vibration disturbances can be assessed more accurately. To assess the impact of these vibration-perturbations, transient accelerometer data from a space shuttle mission are used as inputs for the simulation model. These on-orbit acceleration data were obtained from the Microgravity Science Division at Glenn Research Center (GRC- MSD) and are applied to the buoyancy term of the momentum equation in a simulation of a Pb-5.8 wt. % Sb alloy that solidifies in a thermal gradient of 4000 K/m and a translation velocity of 3 p d s . Figure 2 shows the vertical velocity of a node that begins in the all-liquid region and subsequently solidifies; the vibrations are applied at 5000 seconds in this simulation. An important difficulty, common to all solidification models based on finite elements or 2 The magnitudes of the velocity oscillations that are vibration-induced are very small and acceptable. The biggest concern is whether the concentration of the liquid near the dendrite tips is distorted because of the vibration-induced perturbations. Results for this case show no concentration oscillations present in the all-liquid region.

Westra, D. G.↗
NASA NTRS2001

Remobilizing the Interface of Thermocapillary Driven Bubbles Retarded By the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow that propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillary to direct the movement of bubbles in space is the fact that surfactant impurities, which are unavoidably present in the continuous phase, can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. (This can also be the case if kinetic desorption of surfactant at the back end of the bubble is much slower than convection.) For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like and bubbles translate as solid particles). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature-induced tension gradient, and can decrease to near zero the bubble's thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have be retarded by the adsorption onto the bubble surface of a surfactant impurity. Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble surface much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition, the rapid kinetic exchange and high bulk concentration maintain a saturated surface with uniform surface concentrations. This prevents retarding surface tension gradients and keeps the thermocapillary velocity high. In our reports over the first 2 years, we presented numerical simulations of the bubble motion and surfactant transport which verified theoretically the concept of remobilization, and the development of an apparatus to track and measure the velocity of rising bubbles in a glycerol/water surfactant solution. This year, we detail experimental observations of remobilization. Two polyethylene oxide surfactants were studied, C12E6 (CH3(CH2)11(OCH2)6OH) and C10E8 (CH3(CH2)4(OCH2CH2)8OH). Measurements of the kinetic exchange for these surfactants show that the one with the longer hydrophobe chain C12E6 has a lower rate of kinetic exchange. In addition, this surfactant is much less soluble in the glycerol/water mixture because of the shorter ethoxylate chain. As a result, we found that C12E6 had only a very limited ability to remobilize rising bubbles because of the limited kinetic exchange and reduced solubility. However, C10E8, with its higher solubility and more rapid exchange was found to dramatically remobilize rising bubbles. We also compared our theoretical calculations to the experimental measurements of velocity for both the non-remobilizing and remobilizing surfactants and found excellent agreement. We further observed that for C10E8 at high concentrations, which exceeded the critical micelle concentrations, additional remobilization was measured. In this case the rapid exchange of monomer between micelle and surfactant provides an additional mechanism for maintaining a uniform surface concentrations.

Palaparthi, Ravi↗
NASA NTRS1992

The formation of ore mineral deposits on the Moon: A feasibility study

Most of the ore deposits on Earth are the direct result of formation by hydrothermal solutions. Analogous mineral concentrations do not occur on the Moon, however, because of the absence of water. Stratified ore deposits form in layered instrusives on Earth due to fractional crystallization of magma and crystal settling of high-density minerals, particularly chromium in the mineral chromite. We have evaluated the possibility of such mineral deposition on the Moon, based upon considerations of 'particle settling velocities' in lunar vs. terrestrial magmas. A first approximation of Stoke's Law would seem to indicate that the lower lunar gravity (1/6 terrestrial) would result in slower crystal settling on the Moon. However, the viscosity of the silicate melt is the most important factor affecting the settling velocity. The viscosities of typical lunar basaltic melts are 10-100 times less than their terrestrial analogs. These lower viscosities result from two factors: (1) lunar basaltic melts are typically higher in FeO and lower in Al2O3, Na2O, and K2O than terrestrial melts; and (2) lunar igneous melts and phase equilibria tend to be 100-150 C higher than terrestrial, largely because of the general paucity of water and other volatile phases on the Moon. Therefore, particle settling velocities on the Moon are 5-10 times greater than those on Earth. It is highly probable that stratiform ore deposits similar to those on Earth exist on the Moon. The most likely ore minerals involved are chromite, ilmenite, and native FeNi metal. In addition, the greater settling velocities of periodotite in lunar magmas indicate that the buoyancy effects of the melt are less than on Earth. Consequently, the possibility is considerably less than on Earth of deep-seated volcanism transporting upper mantle/lower crustal xenoliths to the surface of the Moon, such as occurs in kimberlites on Earth.

Taylor, Lawrence A.↗
NASA NTRS2003

Influence of Solutocapillary Convection on Macrovoid Defect Formation in Polymeric Membranes

The focus of this research project involved the dry-cast process for polymeric membrane formation, whereby evaporation of solvent from an initially homogeneous polymer/solvent/ nonsolvent solution results in phase separation and the formation of polymer-rich and polymer-lean phases. Under certain conditions the polymer-lean phase gives rise to very large and usually undesirable, tear-drop-shaped pores (size approx. 10 - 50 microns) termed macrovoids (MVs). Although in many cases the presence of MV pores has deleterious effects on membrane performance, there are a number of innovative applications where the presence of such pores is highly desirable. Although researchers have proposed a variety of mechanisms for MV formation over the past three decades, two main hypotheses are currently favored: one asserts that MV growth can be attributed solely to diffusion (the diffusive growth hypothesis), whereas the other states that solutocapillary convection (the SC hypothesis) at the MV interface contributes to growth. The overall goal of this research was to obtain a more comprehensive understanding of the fundamental mechanism of MV growth. This research incorporates a coupled modeling and experimental approach to test a solutocapillary convection hypothesis for the growth of macrovoid pores in polymeric membranes. Specifically, we utilized a modification of the first principles model developed by two of the PIs (ARG and WBK) for dry-cast CA membranes. For the experimental component, two separate and mutually complementary approaches were used to study MV growth. In the first, membranes cast in a zero-g environment aboard the NASA KC-135 aircraft were compared with those cast on the ground to assess the effect of the buoyancy force on membrane morphology and MV size and shape. In the second approach, videomicroscopy flow visualization (VMFV) was utilized to observe MV formation and growth in real time and to assess the effect of surface tension on the MV growth dynamics. As a result of these fundamental studies, our research group advanced a new hypothesis for MV pore development in polymeric membranes.

Greenberg, Alan R.↗
NASA NTRS2014

Submerged AUV Charging Station

Autonomous Underwater Vehicles (AUVs) are becoming increasingly important for military surveillance and mine detection. Most AUVs are battery powered and have limited lifetimes of a few days to a few weeks. This greatly limits the distance that AUVs can travel underwater. Using a series of submerged AUV charging stations, AUVs could travel a limited distance to the next charging station, recharge its batteries, and continue to the next charging station, thus traveling great distances in a relatively short time, similar to the Old West “Pony Express.” One solution is to use temperature differences at various depths in the ocean to produce electricity, which is then stored in a submerged battery. It is preferred to have the upper buoy submerged a reasonable distance below the surface, so as not to be seen from above and not to be inadvertently destroyed by storms or ocean going vessels. In a previous invention, a phase change material (PCM) is melted (expanded) at warm temperatures, for example, 15 °C, and frozen (contracted) at cooler temperatures, for example, 8 °C. Tubes containing the PCM, which could be paraffin such as pentadecane, would be inserted into a container filled with hydraulic oil. When the PCM is melted (expanded), it pushes the oil out into a container that is pressurized to about 3,000 psi (approx equals 20.7 MPa). When a valve is opened, the high-pressure oil passes through a hydraulic motor, which turns a generator and charges a battery. The low-pressure oil is finally reabsorbed into the PCM canister when the PCM tubes are frozen (contracted). Some of the electricity produced could be used to control an external bladder or a motor to the tether line, such that depth cycling is continued for a very long period of time. Alternatively, after the electricity is generated by the hydraulic motor, the exiting low-pressure oil from the hydraulic motor could be vented directly to an external bladder on the AUV, such that filling of the bladder causes the AUV to rise, and emptying of the bladder allows the AUV to descend. This type of direct buoyancy control is much more energy efficient than using electrical pumps in that the inefficiencies of converting thermal energy to electrical energy to mechanical energy is avoided. AUV charging stations have been developed that use electricity produced by waves on floating buoys and that use electricity from solar photovoltaics on floating buoys. This is the first device that has absolutely no floating or visible parts, and is thus impervious to storms, inadvertent ocean vessel collisions, or enemy sabotage.

Jones, Jack A.↗
NASA NTRS1995

Resolution of the buoyancy in the 8-foot high temperature tunnel combustor

Currently, the 8-Foot High Temperature Tunnel (8-Ft. HTT) combustor produces a good profile at only one point (2000 psia and 3650 R with oxygen enrichment). Air is enriched with oxygen (liquid) so that the combustor product gas will contain the volumetric amount of oxygen normally found in air. The oxygen enriched air has a large fraction that is not reacted and flows through the outer periphery of the fuel injector. This ring of cold air in addition to the relatively cold walls of the combustor set up buoyancy forces that produce a segregation of relatively cool gases at the bottom of the combustor exit. The basic problem is to produce a test gas that has uniform properties at all combustor conditions. The combustor temperature may be as high as 3700 R or as low as 2000 R. Combustor pressures can be as high as 3500 psia (no oxygen enrichment) and as low as 600 psia. The segregation is most severe with oxygen enriched air, since its temperature is lower and its density is high. The combustor is lined with nickel 201 and can be operated at about 1600 R maximum. A global mixing process is desired that produces an acceptable profile of temperature, species, and velocity at the exit of the combustor. The ultimate goal is a temperature profile with about 100 R variance and about 2 percent variance in oxygen. The exit total temperature must not be lowered significantly by the mixing apparatus or mechanisms employed. If immersed bodies are used, they must also be kept very hot. All combustor wall modifications must be able to survive the heat and structural conditions of the varied operating conditions. Our approach to resolving this issue is being conducted in three stages: (1) Consider mixing exclusively, (2) Resolve the heat transfer concerns resulting from the chosen mixing strategy, and (3) Solve the material and structural problems resulting from stages (1) and (2). Since the 8-Ft. HTT is unavailable for experimentation, the study is conducted exclusively with computational fluid dynamic (CPD) codes (Fluent/Uns and Rampant 3.1) using unstructured grid through body fitted coordinates. Both CFD codes are general purpose Navier-Stokes solution packages that can solve integral conservation equations for conservation of mass, momentum and energy. The governing equations are discretized using a control-volume finite-element method on unstructured triangular 2-D grids. In the interest of time, a 3-D tetrahedral grid was used to check the 2-D results on one mixing strategy and the 2-D results were confirmed. Preliminary results indicate that excellent mixing can be achieved with a body placed in the center of the flowing hot fluid with a minor modification to the combustor wall similar to a model positioned in the test section of a wind tunnel. The concept here, is to create longitudinal vortices strong enough to bring both fluids into intimate contact with each other near the wall where the cold fluid resides. However, there may be a trade-off in the choice of mixing strategy, heat transfer and structural requirements. Currently, we are examining the more promising geometries for heat transfer concerns and developing strategies for the material of construction of the center body.

Loney, Norman W.↗
NASA NTRS1995

Solid surface combustion experiment flame spread in a quiescent, microgravity environment implications of spread rate and flame structure

A unique environment in which flame spreading, a phenomenon of fundamental, scientific interest, has importance to fire safety is that of spacecraft in which the gravitational acceleration is low compared with that of the Earth, i.e., microgravity. Experiments aboard eight Space Shuttle missions between October 1990 and February 1995 were conducted using the Solid Surface Combustion Experiment (SSCE) payload apparatus in an effort to determine the mechanisms of gas-phase flame spread over solid fuel surfaces in the absence of any buoyancy induced or externally imposed oxidizer flow. The overall SSCE effort began in December of 1984. The SSCE apparatus consists of a sealed container, approximately 0.039 cu m, that is filled with a specified O2/N2 mixture at a prescribed pressure. Five of the experiments used a thin cellulosic fuel, ashless filter paper, 3 cm wide x 10 cm long, 0.00825 cm half-thickness, ignited in five different ambient conditions. Three of the experiments, the most recent, used thick polymethylmethacrylate (PMMA) samples 0.635 cm wide x 2 cm long, 0.32 cm half-thickness. Three experiments, STS 41, 40 and 43, were designed to evaluate the effect of ambient pressure on flame spread over the thin cellulosic fuel while flights STS 50 and 47 were at the same pressure as two of the earlier flights but at a lower oxygen concentration in order to evaluate the effect of ambient oxygen level on the flame spread process at microgravity. For the PMMA flights, two experiments, STS 54 and 63, were at the same pressure but different oxygen concentrations while STS 64 was at the same oxygen concentration as STS 63 but at a higher pressure. Two orthogonal views of the experiments were recorded on 16 mm cine-cameras operating at 24 frames/s. In addition to filmed images of the side view of the flames and surface view of the burning samples, solid- and gas-phase temperatures were recorded using thermocouples. The experiment is battery powered and follows an automated sequence upon activation by the Shuttle Crew. In this study we separate the SSCE data into two groups according to the fuel type: (1) thin cellulose; and (2) thick PMMA. The experimental spread rates are compared with prediction from a number of models in an effort to uncover the important physics that characterize microgravity flame spread. Both steady and unsteady solutions are employed to explore the flame evolution, especially for thick fuels. Finally, the flame structure in downward spread is compared with the microgravity flame structure and modeling results to delineate the difference between the two configurations and the influence of normal gravity.

Bundy, Matthew↗
DOE OSTI2021

Powering the Blue Economy: A Survey of Station-Keeping Methods for Mooringless Platforms

The term “ocean platform” is used to reference everything from stationary, typically moored, buoys to mobile water vehicles, whether they operate on the ocean’s surface or underwater. For certain applications for which relatively stationary station-keeping conditions are desired, the use of mooring systems is not always a viable alternative either for economic, environmental, regulatory, or otherwise practical reasons, or a combination thereof, (e.g., short deployments, sensitive ecosystems, very deep project sites). Maintaining a platform at a single waypoint or reference location without being moored would require additional control systems and a power source to counteract the drift forces that would naturally displace it. Mobile platforms, which are usually untethered except for remotely operated vehicles, typically require energy input to power their station-keeping capabilities so that they hold or control their location in the ocean. Currently, most of these platforms use combustion engines or batteries for this purpose, which, depending on the specific systems, may be costly, pollute the environment, or create limitations on the length of the deployment. However, powering this kind of platforms with surrounding renewable resources (waves, currents, winds, or sun) has been identified as a promising solution to expand their application. The intent of this report is to investigate station-keeping methods for various ocean platforms that are not moored or otherwise anchored to the ocean floor, or another platform or vessel, paying particular interest to technologies that use marine renewable resources to power their operation, because that is of particular interest to the U.S. Department of Energy’s Powering the Blue Economy (PBE) initiative. As a first step, 72 articles and technical reports related to mooringless station-keeping methods were collected for review. The preliminary literature review provided a broad overview of common themes across the literature from which a descriptive methodology for analyzing various platforms was developed. That is, station-keeping methods were categorized based on their predominant energy source and consumption (renewable, nonrenewable, or hybrid if the platform uses renewable and nonrenewable resources equally), and their localization strategy (drift reduction, “path-planning or “waypoint-holding”). In addition, platform types were segregated into the following groups: buoys, surface drifters, and unoccupied surface vehicles (USVs); offshore renewable energy systems; and unoccupied underwater vehicles (UUVs). The main types of station-keeping methods encountered in this report achieve their intended localization strategy by means of drift mitigation, steering, and/or propulsion. Drift mitigation is commonly accomplished via drogues and sea anchors. Stand-along steering subsystems use control surfaces (e.g., ship rudder, wing sail, etc.) that react to ocean currents, waves, or winds to provide varying-degrees of course adjustments. Combined steering and propulsion subsystems include differential thrusters, directional thrusters separate from a primary thruster that cause the platform to pitch up/down or yaw clockwise/counterclockwise, or vectored thrusters that direct the propulsion in a range of directions relative to the platform’s local coordinate system. Propulsion is often achieved by running a motor and applying active control strategies but can also involve buoyancy shifts and using sails to generate lifting forces that propel a platform in a desired direction. Future research is primarily expected to take place in the form of a technoeconomic analysis that would aim to determine the technological viability, cost, and added value of mooringless station-keeping use cases identified through this research, including docking for UUV recharging or for georeferencing drifter buoys, deep-sea floating wind farms, U.S. Navy sonar arrays, and a Pacific Ocean wave buoy network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1999

Experimental Studies of the Interaction Between a Parallel Shear Flow and a Directionally-Solidifying Front

It has long been recognized that flow in the melt can have a profound influence on the dynamics of a solidifying interface and hence the quality of the solid material. In particular, flow affects the heat and mass transfer, and causes spatial and temporal variations in the flow and melt composition. This results in a crystal with nonuniform physical properties. Flow can be generated by buoyancy, expansion or contraction upon phase change, and thermo-soluto capillary effects. In general, these flows can not be avoided and can have an adverse effect on the stability of the crystal structures. This motivates crystal growth experiments in a microgravity environment, where buoyancy-driven convection is significantly suppressed. However, transient accelerations (g-jitter) caused by the acceleration of the spacecraft can affect the melt, while convection generated from the effects other than buoyancy remain important. Rather than bemoan the presence of convection as a source of interfacial instability, Hurle in the 1960s suggested that flow in the melt, either forced or natural convection, might be used to stabilize the interface. Delves considered the imposition of both a parabolic velocity profile and a Blasius boundary layer flow over the interface. He concluded that fast stirring could stabilize the interface to perturbations whose wave vector is in the direction of the fluid velocity. Forth and Wheeler considered the effect of the asymptotic suction boundary layer profile. They showed that the effect of the shear flow was to generate travelling waves parallel to the flow with a speed proportional to the Reynolds number. There have been few quantitative, experimental works reporting on the coupling effect of fluid flow and morphological instabilities. Huang studied plane Couette flow over cells and dendrites. It was found that this flow could greatly enhance the planar stability and even induce the cell-planar transition. A rotating impeller was buried inside the sample cell, driven by an outside rotating magnet, in order to generate the flow. However, it appears that this was not a well-controlled flow and may also have been unsteady. In the present experimental study, we want to study how a forced parallel shear flow in a Hele-Shaw cell interacts with the directionally solidifying crystal interface. The comparison of experimental data show that the parallel shear flow in a Hele-Shaw cell has a strong stabilizing effect on the planar interface by damping the existing initial perturbations. The flow also shows a stabilizing effect on the cellular interface by slightly reducing the exponential growth rate of cells. The left-right symmetry of cells is broken by the flow with cells tilting toward the incoming flow direction. The tilting angle increases with the velocity ratio. The experimental results are explained through the parallel flow effect on lateral solute transport. The phenomenon of cells tilting against the flow is consistent with the numerical result of Dantzig and Chao.

Zhang, Meng↗
DOE OSTI2022

Development of Numerical Model of Metal Foam with PCM for the Estimation of Effective Thermal Conductivity

Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result

25 ENERGY STORAGE↗
NASA NTRS1999

Control of Meridional Flow in Circular Cylinders by a Travelling Axial Magnetic Field

Convective flow in a Bridgman or float zone configuration significantly affects the interface shape and segregation phenomena. While the primary causative factor for this flow is buoyancy induced convection in an enclosed Bridgman melt, the presence of a free surface gives rise to surface tension driven flows in the floating zone processing of melts. It is of interest to curtail these flows in order to realize near quiescent growth conditions that have shown to result in crystals with good longitudinal and radial homogeneity and thereby of better overall quality. While buoyancy effects can be reduced by careful processing in a low gravity (space) environment, the reduction of Marangoni flows due to surface tension variations is not that straight forward. Attempts have been made with some limited success with the use of external fields to affect the melt thermo-fluid behavior. The use of a static magnetic field that reduces convective contamination through the effects of a non-intrusively induced, dissipative Lorentz force in an electrically conducting melt is one such approach. Experiments have shown that axial fields of the order of 5 Tesla can significantly eliminate convection and yield close to diffusion limited crystal growth conditions. The generation and use of such high magnetic fields require substantial hardware and incur significant costs for its operation. Lately, the use of rotating magnetic fields has been tested in semiconductor crystal growth. The method is fairly well known and commonly used in metal processing but its adaptation to crystal growth of semiconductors is fairly recent. The elegance of the technique rests in its low power requirement (typically 10-20 milli-Tesla at 50-400 Hz) and its efficacy in curtailing deleterious temperature fluctuations in the melt. A rotating magnetic field imposes a rotational force and thereby induces a circulation within the melt that tends to dominate other sporadic convective effects. Thus a known low level of convective flow is introduced into the system. A new novel variation of the Lorentz force mechanism is proposed and investigated in this study. Since one of the desired process conditions in melt crystal growth is the minimization of convective effects, this investigation examines the use of an external field of magnetic origin to counteract existing convective flow within the melt. This is accomplished by utilizing a running or traveling axial magnetic wave in the system. The concept is similar to the use of vibrational means in order to induce streaming flows that oppose buoyant or surface tension driven convection in the system. The rotation direction as well as the magnitude (strength) of this circulation can be easily controlled by external inputs thus affording a direct means of controlling the developing shape of the crystallizing front (interface). The theoretical model of this technique is fully developed and presented in this paper. Results from the solution of the developed governing equations and boundary conditions are also presented. An experimental demonstration of the concept is presented through the suppression of natural convective flow in a mercury column. Implications to crystal growth systems will be fully explored in the final manuscript.

Mazuruk, K.↗