Search NASA⌕ Search

Engineering topics

Tryggvason, Gretar

Publications and source records attributed to Tryggvason, Gretar.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Shape effects on the local dynamics of suspensions of spheroidal particles

The effect of shape on the dynamics of suspensions of non-spherical heavy particles is examined by fully resolved numerical simulations of oblate and prolate spheroids, as well as spheres, for a density ratio of ten, volume fractions ranging from 0.5% to 5%, and Reynolds numbers between 20 and 30. The dynamics is determined both by the interactions of the particles with the fluid as well as by collisions, with the number and importance of collisions increasing with volume fractions. A single isolated oblate or prolate spheroid falling under gravity generally falls broadside on, for the governing parameters examined here, and at low-volume fractions, the majority of particles in a suspension fall that way. At higher-volume fractions, the orientation is more random. The slip velocity decreases as the volume fraction increases for all shapes, as expected, but the effect of the shape is much less than seen for a single particle. Furthermore, this seems to be due to two effects. For all volume fractions, the most deformed particles cluster more than spheres and less deformed particles, which increases their slip velocity. As the concentration increases, the increased particle interactions also causes more particles to fall short side-on, which reduces the frontal area and the resulting drag, increasing the slip velocity. This second effect is, of course, absent for spherical particles.

42 ENGINEERING↗
DOE OSTI2023

The dynamics of suspensions of prolate spheroidal particles—Effects of volume fraction

Here, the effect of volume fraction on the dynamics of a gravity driven suspension of prolate spheroidal solid particles in a fully periodic domain is examined by fully resolved numerical simulations, for relatively modest Reynolds numbers (around 20, depending on the volume fraction). Three systems are examined, for volume fractions of 2.2% (20 particles), 5.5% (50 particles), and 9.9% (90 particles). The results show a transition from flow at low volume fractions dominated by the hydrodynamic interactions between the particles and the fluid, although modified by collisions, to flow dominated by collisions at higher volume fractions. The distribution of the particles with respect to each other, as measured by the probability distribution of nearest distances, is nearly random at low volume fractions but at high volume fractions, there is more clustering than for a random distribution. At lower volume fractions most of the particles fall broadside on whereas at the highest volume fraction, their orientation is essentially random.

42 ENGINEERING↗
DOE OSTI2021

Effect of topology changes on the breakup of a periodic liquid jet

Here, the breakup of a periodic jet is examined computationally, using a front-tracking/finite-volume method, where the interface is represented by connected marker points moving with the fluid, while the governing equations are solved on a fixed grid. Tracking the interface allows control of whether topology changes take place or not. The Reynolds and Capillary numbers are kept relatively low ($Re = 150$ and $Ca = 2$) so most of the flow is well resolved. The effect of topology changes is examined by following the jet until it has mostly disintegrated, for different “coalescence criterion,” based on the thickness of thin films and threads. The evolution of both two-dimensional and fully three-dimensional flows is examined. It is found that although there is a significant difference between the evolution when no breakup takes place and when it does, once breakup takes place the evolution is relatively insensitive to exactly how it is triggered for a range of coalescence criterion, and any differences are mostly confined to the smallest scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Verification of bubble tracking method and DNS examinations of single- and two-phase turbulent channel flows

Direct numerical simulation (DNS) has been regarded as a reliable data source for the development and validation of turbulence models along with experiments. The realization of DNS usually involves a very fine mesh that should be able to resolve all relevant turbulence scales down to Kolmogorov scale [1]. As the most computationally expensive approach compared to other CFD techniques, DNS applications used to be limited to flow studies at very low Reynolds numbers. Thanks to the tremendous growth of computing power over the past decades, the simulation capability of DNS has now started overlapping with some of the most challenging engineering problems. One of those examples in nuclear engineering is the turbulent coolant flow inside reactor cores. Coupled with interface tracking methods (ITM), the simulation capability of DNS can be extended to more complicated two-phase flow regimes. Departure from nucleate boiling (DNB) is the limiting critical heat flux phenomena for the majority of accidents that are postulated to occur in pressurized water reactors (PWR) [2]. As one of the major modeling and simulation (M&S) challenges pursued by CASL, the prediction capability is being developed for the onset of DNB utilizing multiphase-CFD (M-CFD) approach. DNS (coupled with ITM) can be employed to provide closure law information for the multiphase flow modeling at CFD scale. In the presented work, research groups at NCSU and UND will focus on applying different ITM to different geometries. Higher void fraction flow analysis at reactor prototypical conditions will be performed, and novel analysis methods will be developed, implemented and verified for the challenging flow conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Using statistical learning to close two-fluid multiphase flow equations for a simple bubbly system

Direct numerical simulations of bubbly multiphase flows are utilized to find closure terms for a simple model of the average flow, using Neural Networks (NNs). The flow considered consists of several nearly spherical bubbles rising in a periodic domain where the initial vertical velocity and the average bubble density are homogeneous in two directions but non-uniform in one of the horizontal directions. After an initial transient motion the average void fraction and vertical velocity become approximately uniform. The NN is trained on a dataset from one simulation and then used to simulate the evolution of other initial conditions. As a whole, the resulting model predicts the evolution of the various initial conditions reasonably well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Investigations of Bubbly Flows using Direct Numerical Simulations

Direct numerical simulations are used to examine the behavior of relatively high void fraction bubbly flows. The full Navier-Stokes equations are solved for both the bubbles and the ambient liquid and the flow near the moving and deforming phase boundaries is resolved fully. This yields detailed data for the evolution of the system that can be used to extract information about the basic bubble-bubble and bubble-flow interactions as well as about the averaged properties of the system. Such information are critical for the generation of averaged engineering models of industrial size systems as well as for novel applications of bubbly flows where experiments are difficult. The project has yielded detailed information about many aspects of bubbly flows. The results for flows in vertical channels are particularly significant and have yielded an essentially complete understanding of many aspects of the flows. The important insight is that the flow has a particularly simple structure, consisting of a homogeneous core in hydrostatic balance and a wall layer whose composition depends of whether the liquid is flowing upward or downward. These results hold for both laminar and turbulent flows and have allowed us to derive an analytical model to predict the flow. We have also shown that bubble deformability results in a completely different flow structure, but also giving rise to relatively simple predictive models. Preliminary studies of transient migration of a non-homogeneous bubble distribution and coalescence induced flow regime transitions in vertical channels were also done.

42 ENGINEERING↗
NASA NTRS2004

Spray Cooling Processes for Space Applications

The present paper reports ongoing work to develop numerical and modeling tools used to design efficient and effective spray cooling processes and to determine characteristic non-dimensional parametric dependence for practical fluids and conditions. In particular, we present data that will delineate conditions towards control of the impingement dynamics of droplets upon a heated substrate germane to practical situations.

Kizito, John P.↗
NASA NTRS2004

Droplet-Surface Impingement Dynamics for Intelligent Spray Design

Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1) By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface; 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling; and 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions. Determining the splash/non-splash boundary is important for many practical applications. Coating and cooling processes would each benefit from near-term empirical relations and subsequent models. Such demonstrations can guide theoretical development by providing definitive testing of its predictive capabilities. Thus, empirical relations describing the boundary between splash and non-splash are given for drops impinging upon a dry solid surface and upon a thin fluid film covering a similar surface. Analytical simplification of the power laws describing the boundary between the splash and non-splash regions yields insight into the engineering parameters governing the splash and non-splash outcomes of the fluid droplets. The power law correlation is shown separating the splashing versus non-splashing regions as developed for droplets impinging upon a dry solid surface. Splashing upon a dry surface is reasonably described by Ca greater than 0.85, reflecting the competing roles of surface tension and viscosity. The power law correlation is shown separating the splashing versus non-splashing regions as developed for droplets impinging upon a thin fluid film covering the solid surface. Splashing upon a thin fluid film, as described by v (pd/s) greater than 63, is governed by fluid density and surface tension, but is rather independent of viscosity. Finally, the data presented here suggests that a more direct dependence upon the surface tension and viscosity, given a better understanding of their interplay, would allow accurate description of the droplet-surface impacts for more complicated situations involving non-Newtonian fluids, specifically those exhibiting viscoelastic behavior.

Wal, Randy L. Vander↗
NASA NTRS2004

Droplet-surface Impingement Dynamics for Intelligent Spray Design

Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1. By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling. 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions.

VanderWal, Randy L.↗
NASA NTRS2002

Splashing Droplets

Current data on droplet breakup is scarce for the sizes and velocities typical of practical applications such as in spray combustion processes and coating processes. While much more representative of practical applications, the small spatial scales and rapid time-scales prevent detailed measurement of the internal fluid dynamics and liquid property gradients produced by impinging upon surfaces. Realized through the extended spatial and temporal scales afforded by a microgravity environment, an improved understanding of drop breakup dynamics is sought to understand and ultimately control the impingement dynamics of droplets upon surfaces in practical situations. The primary objective of this research will be to mark the onset of different 'splashing modes' and to determine their temperature, pressure and angle dependence for impinging droplets representative of practical fluids. In addition, we are modeling the evolution of droplets that do not initially splash but rather undergo a 'fingering' evolution observed on the spreading fluid front and the transformation of these fingers into splashed products. An example of our experimental data is presented below. These images are of Isopar V impacting a mirror-polished surface. They were acquired using a high-speed camera at 1000 frames per second. They show the spreading of a single droplet after impact and ensuing finger instabilities. Normal gravity experimental data such as this will guide low gravity measurements in the 2.2 second drop tower and KC-135 aircraft as available. Presently we are in the process of comparing the experimental data of droplet shape evolution to numerical models, which can also capture the internal fluid dynamics and liquid property gradients such as produced by impingement upon a heated surface. To-date isothermal numerical data has been modeled using direct numerical simulations of representative splashing droplets. The data obtained so far indicates that the present model describes well the droplet wall interactions to a point in time just before splash. Additional information is included in the original extended abstract.

VanderWal, Randall L.↗
NASA NTRS2001

Concentration Field of Reactants and Products Species in a Reacting Vortex Ring

The proposed paper will present experimental and numerical results on the concentration fields of both reactants and products species in a reacting vortex ring that is generated from the interaction between a diffusion flame and a laminar vortex ring. Flame-vortex interactions are canonical configurations used to study the underlying processes occurring in complicated turbulent reacting flows. This type of configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. The current configuration has been studied experimentally by Chen and Dahm and Chen et al. under microgravity conditions, and by Park and Shin, and You et al. under normal gravity conditions. This configuration is similar to that used in the analyses of Karagozian and Manda of their 2-D vortex pair in which both fuel and entrained oxidizer are present. The vortex ring used in this study is generated by issuing methane into an air environment through the exit of an axisymmetric nozzle. The experiments were conducted under microgravity conditions in order to remove the undesirable effects of buoyancy that can affect both the flame structure and ring dynamics resulting in possibly asymmetric and nonrepeatable interactions. The experimental technique of diode laser wavelength modulation spectroscopy (WMS) is used to measure concentration fields of reactants, CH4 and O2, products, H2O, CO2, OH, and temperature fields which can be inferred from either line pairs of O2 or OH lines. This technique has been investigated previously by Silver and Bomse et al. This is the first time that the technique has been applied to reacting vortex rings under microgravity conditions. The effect of ring circulation and fuel volume on the species concentration fields will be investigated. The experimental results will be compared to the current numerical results, and used to validate the numerical studies. In addition, the existence of burned cores during the interactions will be determined, and the increase in reactant consumption with increased ring circulation will be examined. Numerical studies were also conducted by solving the Navier-Stokes and mixture fraction equations with the assumptions of unity Lewis and Schmidt numbers. Equilibrium chemistry and flamelet libraries were used to obtain the temperature and species mass fraction fields. The numerical results will serve as guidelines in conducting the experimental studies. Ring circulation and fuel volume effects on the interactions and species concentration fields will be investigated and compared to experimental results.

Chen, Shin-Juh↗
NASA NTRS1999

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗
NASA NTRS1997

Fundamentals of Mold Free Casting: Experimental and Computational Studies

Researchers are developing the technology of 'Ballistic Particle Manufacturing' (BPM) in which individual drops are precisely layered onto a substrate, and the drops are deposited so as to prevent splatting. These individual drops will ultimately be combined to form a net-shape, three-dimensional object. Our understanding of controlled drop deposition as applied to BPM is far from complete. Process parameters include the size and temperature of the liquid metal drop, its impact velocity and trajectory, and the condition and temperature of the substrate. Quantitative knowledge of the fluid mechanics and heat transfer of drop deposition and solidification are necessary to fully optimize the manufacturing process and to control the material microstructure of the final part. The object of this study is to examine the dynamics of liquid metal drops as they impinge upon a solid surface and solidify under conditions consistent with BPM (i.e. conditions which produce non-splatting drops). A program of both numerical simulations and experiments will be conducted. Questions this study will address include the following: How do the deformation and solidification of the drop depend on the properties of the fluid drop and the solid substrate? How does the presence of previously deposited drops affect the impingement and solidification process? How does the impingement of the new drop affect already deposited material? How does the cooling rate and solidification of the drops influence the material microstructure?

Tryggvason, Gretar↗