Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hele-shaw flows”

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.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Flamelet Formation In Hele-Shaw Flow

A Hele-Shaw flow apparatus constructed at Michigan State University (MSU) produces conditions that reduce influences of buoyancy-driven flows. In addition, in the MSU Hele-Shaw apparatus it is possible to adjust the heat losses from the fuel sample (0.001 in. thick cellulose) and the flow speed of the approaching oxidizer flow (air) so that the "flamelet regime of flame spread" is entered. In this regime various features of the flame-to-smolder (and vice versa) transition can be studied. For the relatively wide (approx. 17.5 cm) and long (approx. 20 cm) samples used, approximately ten flamelets existed at all times. The flamelet behavior was studied mechanistically and statistically. A heat transfer analysis of the dominant heat transfer mechanisms was conducted. Results indicate that radiation and conduction processes are important, and that a simple 1-D model using the Broido-Shafizadeh model for cellulose decomposition chemistry can describe aspects of the flamelet spread process. Introduction

Wichman, I. S.↗
NASA NTRS1995

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗
NASA NTRS1999

Premixed-Gas Flame Propagation in Hele-Shaw Cells

It is well known that buoyancy and thermal expansion affect the propagation ra and shapes of premixed gas flames. The understanding of such effects is complicated by the large density ratio between the reactants and products, which induces a baroclinic production of vorticity due to misalignment of density and pressure gradients at the front, which in turn leads to a complicated multi-dimensional flame/flow interaction. The Hele-Shaw cell, i.e., the region between closely-spaced flat parallel plates, is probably the simplest system in which multi-dimensional convection is presents consequently, the behavior of fluids in this system has been studied extensively (Homsy, 1987). Probably the most important characteristic of Hele-Shaw flows is that when the Reynolds number based on gap width is sufficiently small, the Navier-Stokes equations averaged over the gap reduce to a linear relation, namely a Laplace equation for pressure (Darcy's law). In this work, flame propagation in Hele-Shaw cells is studied to obtain a better understanding of buoyancy and thermal expansion effects on premixed flames. This work is also relevant to the study of unburned hydrocarbon emissions produced by internal combustion engines since these emissions are largely a result of the partial burning or complete flame quenching in the narrow, annular gap called the "crevice volume" between the piston and cylinder walls (Heywood, 1988). A better understanding of how flames propagate in these volumes through experiments using Hele-Shaw cells could lead to identification of means to reduce these emissions.

Sharif, J.↗
NASA NTRS2005

The Dynamics of Miscible Interfaces: Simulations

This research project focuses on the dynamics of interfacial regions between miscible fluids. While much attention has focused on immiscible interfaces in the past, miscible interfaces have been explored to a much lesser degree, so that there are many open questions regarding their dynamics at this time. Among the more pressing issues is the role that nonconventional stresses can play in such interfacial regions. Such stresses are typically not accounted for in efforts to model the dynamics of miscible flows. Our research aims to clarify under which circumstances these stresses do have to be taken into account, and what quantitative approaches are most suitable in this regard. In order to address these issues, we have focused on conducting linear stability analyses and nonlinear simulations for capillary tube and Hele-Shaw flows, and to compare the results with corresponding experiments performed in the labs of our co-investigators Prof. Maxworthy at USC, and Dr. Balasubramaniam at NASA. Over the duration of the project we have, among other things, focused on the effects of variable diffusion coefficients in such flows, and specifically on their influence in the growth of instabilities. Furthermore, our three-dimensional spectral element simulations have made good progress, so that we have come to a point where we can conduct more detailed comparisons with experimental observations. We are currently focusing our efforts on reproducing the tip-splitting instability observed by Maxworthy. Finally, we have discovered a new core-annular flow instability in the Stokes flow regime during the last year. This represents a significant finding, as this instability does not have an immiscible counterpart.

Meiburg, Eckart↗
DOE OSTI2026

Quantum Solver Using Singular Value Decomposition for Computational Fluid Dynamics

Numerical solutions for fluid flow problems are challenging and have been focus of Computational Fluid Dynamics (CFD) research for past several decades. The advent of quantum computing promises exponential speedup in comparison to existing classical methods and alleviate computational constraints posed by CFD problems. Although solutions for most problems of interest in fluid dynamics using quantum computing are distant, recent advances in algorithms, software and hardware provide a path towards realizing this goal. Quantum linear solver algorithms (QLSA) such as Harrow–Hassidim–Lloyd (HHL) and Variational Quantum Linear Solver (VQLS) have been successfully implemented to solve for canonical problems such as Hele-Shaw flow. However, these algorithms still suffer to scale and address problems with ill-conditioned Jacobians. In the current paper, we alleviate these restrictions with a new quantum solver based on Singular Value Decomposition (SVD) and simulate flow past a 2D cylinder. The fidelity of the SVD based quantum solver in predicting the flow past 2D cylinder is computed along with an assessment of errors. Classical and quantum solutions for the flow are compared for different resolutions. Finally, we discuss variation in the solutions based on number of shots used.

Gottiparthi, Kalyan [ORNL] (ORCID:0000000213540255↗
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↗
NASA NTRS2002

The Dynamics of Miscible Interfaces: Simulations

The goal of this experimental/computational investigation (joint with Prof Maxworthy at USC) has been to study the dynamics of miscible interfaces, both from a scientific and a practical point of view, and to prepare a related experiment to be flown on the International Space Station. In order to address these effects, we have focused experimental and computational investigations on miscible displacements in cylindrical capillary tubes, as well as in Hele-Shaw cells. Regarding the flow in a capillary tube, the question was addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations and experiments for miscible flows in capillary tubes. An estimate of the vorticity and streamfunction fields induced by the Kortewegs stresses was derived, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, were seen to be very small, and they have a negligible influence on the tip velocity. As a result, it can be concluded that they are not responsible for the discrepancies between experiments and conventional Stokes simulations. A further part of our investigation focussed on the development of high-accuracy three-dimensional spectral element simulation techniques for miscible flows in capillary tubes, including the effects of variable density and viscosity. Towards this end, the conservation equations are treated in cylindrical coordinates.

Meiburg, Eckart↗
NASA NTRS1999

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗
DOE OSTI2025

Invariant Forms of Dissolution Fingers

Dissolution of fractured and porous media introduces a positive feedback between fluid flow and reactant transport, leading to the emergence of pronounced, fingerlike channels. We investigate the formation of these structures using a microfluidic Hele-Shaw cell with a soluble bottom. Our experiments show that the shape of dissolution fingers is invariant and reveals itself over time as the fingers extend into the system. Here, by combining reactive-transport theory and conformal mapping techniques we derive these invariant forms. We relate these results to natural dissolution fingers in karst landscapes, and illustrate how to determine the groundwater flow rate responsible for their formation based on the finger shape.

Flows in porous media↗
NASA NTRS1992

Gallery of fluid motion

Winning photographs from the ninth Annual Fluid Mechanics Photo Contest held by the Division of Fluid Dynamics of the American Physical Society in November 1991 are displayed. Pictures are presented of vortex flows, a periodic axisymmetric vortex breakdown in a cylinder with a rotating end wall, impacting water drops, a vortex structure behind a disk started from rest, dynamic stall, blue waves, isotherms in turbulent Rayleigh-Benard convection, a square jet, 'vapor lines' emanating from water droplets, viscosity effects on the directional solidification of NH4Cl solution in a Hele-Shaw cell, and stationary cross-flow vortices.

Reed, Helen L.↗
DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗
NASA NTRS1996

Marangoni Effects on the Bubble Dynamics in a Pressure Driven Flow

The motion of air bubbles and water drops in a Hele-Shaw cell filled with a silicone oil has been studied experimentally and theoretically. By adding a predetermined amount of a surfactant to the water drops we attempted to investigate the surfactant influence systematically. While the motion of air bubbles was in reasonable agreement with the predictions of Taylor and Saffman, water drops behaved quite differently in that the translational velocities were smaller by an order of magnitude and their shapes were very unusual as observed previously by Kopf-Sill and Homsy. Assuming that the surrounding fluid wets the solid wall and the bubble (or the drop) surface is rigid due to the surfactant influence, we have estimated the translational velocity of an elliptic bubble. The calculated velocities were in good agreement with the observations indicating that the surfactant influence could retard the bubble motion significantly. The present study also indicates that the unusual bubble shapes are also due to the surfactant influence.

Park, Chang-Won↗
NASA NTRS2001

Design and Fabrication of a Hele-Shaw Apparatus for Observing Instabilities of Diffusion Flames

Examinations of flame fronts spreading over solid fuels in an opposed flow of oxidizer have shown that the flame front fragments into smaller (cellular) flames. These 'flamelets' will oscillate, recombine, or extinguish, indicating that they are in the near extinction limit regime (i.e., to one side of the quenching branch of the flammability map). Onset of unstable cellular flamelet formation for flame spread over thin fuels occurs when a heat-sink substrate is placed a small distance from the underside of the fuel. This heat-sink substrate (or backing) displaces the quenching branch of the flammability map in a direction that causes the instabilities to occur at higher air velocities. Similar near-limit behavior has been observed in other works using different fuels, thus suggesting that these dynamic mechanisms are fuel-independent and therefore fundamental attributes of flames in this near-limit flame spread regime. The objective of this project is to determine the contributions of the hydrodynamic and thermodiffusive mechanisms to the observed formation of flame instabilities. From this, a model of diffusion flame instabilities shall be generated. Previously, experiments were conducted in NASA drop towers, thereby limiting observation time to O(1-5 sec). The NASA tests exhibited flamelet survival for the entire drop time, suggesting that flamelets (i.e., small cellular flames) might exist, if permitted, for longer time periods. By necessity, experiments were limited to thermally thin cellulose fuels (approximately 0.001 in thick): instabilities could form by virtue of faster spread rates over thin fuels. Unstable behavior was unlikely in the short drop time for thicker fuels. In the International Space Station (ISS), microgravity time is unlimited, so both thin and thick fuels can be tested.

Wichman, I. S.↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗
DOE OSTI2022

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Surface Tension and Fingering of Miscible Interfaces

Experiments on miscible, buoyantly unstable reaction-diffusion fronts and non-reacting displacement fronts in Hele-Shaw cells show a fingering-type instability whose wavelengths (lambda*) are consistent with an interfacial tension (sigma) at the front caused by the change in chemical composition, even though the solutions are miscible in all proportions. In conjunction with the Saffman-Taylor model, the relation sigma = K/tau, where tau is the interface thickness and K approximately equal 4 +/- 2 x 10(exp -6) dyne, enables prediction of our measured values of lambda* as well as results from prior experiments on miscible interfaces. These results indicate that even for miscible fluids, surface tension is generally a more significant factor than diffusion in interfacial stability and flow characteristics.

Abib, Mohammed↗
NASA NTRS2013

Multiphase Transport in Porous Media: Gas-Liquid Separation Using Capillary Pressure Gradients International Space Station (ISS) Flight Experiment Development

Advances in the understanding of multiphase flow characteristics under variable gravity conditions will ultimately lead to improved and as of yet unknown process designs for advanced space missions. Such novel processes will be of paramount importance to the success of future manned space exploration as we venture into our solar system and beyond. In addition, because of the ubiquitous nature and vital importance of biological and environmental processes involving airwater mixtures, knowledge gained about fundamental interactions and the governing properties of these mixtures will clearly benefit the quality of life here on our home planet. The techniques addressed in the current research involving multiphase transport in porous media and gas-liquid phase separation using capillary pressure gradients are also a logical candidate for a future International Space Station (ISS) flight experiment. Importantly, the novel and potentially very accurate Lattice-Boltzmann (LB) modeling of multiphase transport in porous media developed in this work offers significantly improved predictions of real world fluid physics phenomena, thereby promoting advanced process designs for both space and terrestrial applications.This 3-year research effort has culminated in the design and testing of a zero-g demonstration prototype. Both the hydrophilic (glass) and hydrophobic (Teflon) media Capillary Pressure Gradient (CPG) cartridges prepared during the second years work were evaluated. Results obtained from ground testing at 1-g were compared to those obtained at reduced gravities spanning Martian (13-g), Lunar (16-g) and zero-g. These comparisons clearly demonstrate the relative strength of the CPG phenomena and the efficacy of its application to meet NASAs unique gas-liquid separation (GLS) requirements in non-terrestrial environments.LB modeling software, developed concurrently with the zero-g test effort, was shown to accurately reproduce observed CPG driven gas-liquid separation phenomena. The design and fabrication of a micropost plate-lamina Hele-Shaw (HS) cell was performed which served as a computationally attainable geometric structure facilitating direct comparison between physical phenomena observed in our laboratory and the LB software predictions.

Wheeler, Richard R., Jr.↗