SEARCH · Engineering Papers

Results for “capillary 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗

NASA NTRS2002

Fluid Mechanics of Capillary-Elastic Instabilities in Microgravity Environment

The aim of this project is to investigate the closure and reopening of lung airways due to surface tension forces, coupled with airway elasticity. Airways are liquid-lined, flexible tubes and closure of airways can occur by a Rayleigh instability of the liquid lining, or an instability of the elastic support for the airway as the surface tension of the air-liquid interface pulls the tube shut, or both. Regardless of the mechanism, the airway is closed because the liquid lining has created a plug that prevents axial gas exchange. In the microgravity environment, surface tension forces dominate lung mechanics and would lead to more prevalent, and more uniformly distributed air-way closure, thereby creating a potential for respiratory problems for astronauts. Once closed the primary option for reopening an airway is by deep inspiration. This maneuver will pull the flexible airways open and force the liquid plug to flow distally by the incoming air stream. Airway reopening depends to a large extent on this plug flow and how it may lead to plug rupture to regain the continuity of gas between the environment and the alveoli. In addition to mathematical modeling of plug flows in liquid-lined, flexible tubes, this work has involved benchtop studies of propagating liquid plugs down tube networks that mimic the human airway tree. We have extended the work to involve animal models of liquid plug propagation in rat lungs. The liquid is radio-opaque and x-ray video imaging is used to ascertain the movement and distribution of the liquid plugs so that comparisons to theory may be made. This research has other uses, such as the delivery of liquids or drugs into the lung that may be used for surfactant replacement therapy or for liquid ventilation.

Grotberg, James B.↗

NASA NTRS1996

Effect of Variable Gravity on Evaporation of Binary Fluids in a Capillary Pore Evaporator

The research project focuses on experimental investigation of the capillary-pumped evaporative heat transfer phenomenon. The objective is to examine whether the heat transfer and stability of a heated meniscus in a capillary pore can be enhanced by adding trace amounts of a non-volatile solute to a solvent and to understand the changes that occur. The experimental setup consists of a single pore evaporator connected to a reservoir which supplies liquid to the evaporator. In addition to the experiments of capillary-pumped evaporation, a parallel experimental study has been conducted to systematically investigate the effects of gravity as well as the effects of bulk composition on the heat transfer characteristics of evaporating binary thin films near the contact line region along an inclined heated surface. To investigate the buoyancy effects on evaporation along an inclined heated surface, the angle of inclination from a horizontal plane was varied fro 15 C to 90 C. An optimum concentration between 0.5% and 1% decane in pentane/decane solutions has been demonstrated at different angles of inclination. Improved heat transfer was found for the geometry with the smallest angle of inclination of 15 degrees. In addition, flow visualization has revealed that at low inclination angles effective heat transfer takes place primarily due to an extension of the thin film near the contact line. At these low inclination angles, the optimum concentration is associated with enhanced wetting characteristics and reduced thermocapillary stresses along the interface.

Girgis, Morris M.↗

NASA NTRS1999

Weakly Nonlinear Description of Parametric Instabilities in Vibrating Flows

This project focuses on the effects of weak dissipation on vibrational flows in microgravity and in particular on (a) the generation of mean flows through viscous effects and their reaction on the flows themselves, and (b) the effects of finite group velocity and dispersion on the resulting dynamics in large domains. The basic mechanism responsible for the generation of such flows is nonlinear and was identified by Schlichting [21] and Longuet-Higgins. However, only recently has it become possible to describe such flows self-consistently in terms of amplitude equations for the parametrically excited waves coupled to a mean flow equation. The derivation of these equations is nontrivial because the limit of zero viscosity is singular. This project focuses on various aspects of this singular problem (i.e., the limit C equivalent to (nu)((g)(h(exp 3)))exp -1/2 << 1,where nu is the kinematic viscosity and h is the liquid depth) in the weakly nonlinear regime. A number of distinct cases is identified depending on the values of the Bond number, the size of the nonlinear terms, distance above threshold and the length scales of interest. The theory provides a quantitative explanation of a number of experiments on the vibration modes of liquid bridges and related experiments on parametric excitation of capillary waves in containers of both small and large aspect ratio. The following is a summary of results obtained thus far.

Knobloch, E.↗

DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Additive Manufacturing and Experimental Characterization of Nickel-Titanium Shape-Memory Alloy Wick Structures and Heat Pipes for Spacecraft Thermal Control

Shape memory alloys (SMA), such as those based on nickel-titanium (NiTi), are increasingly being applied as multifunctional spacecraft components. For thermal management applications, NiTi flow tubing hinges and self-deploying loop heat pipes have been demonstrated. Emerging additive manufacturing (AM) processes are enabling more complex SMA devices than can be formed from conventional plain wire, tubing, and sheet stock materials. This paper presents our progress toward applying powder bed fusion AM to producing porous NiTi wicks and NiTi-H2O heat pipes, which could be embedded in thermally deploying radiators for spacecraft thermal management. AM near-equiatomic NiTi (55.1 wt% Ni) porous wick specimens were produced with a range of deposition parameters. Transient acetone rate-of-rise experiments were performed to estimate wick permeability (K) and average pore radius (r_pore) values and identify parameter sets with high capillary performance. Surface treatments were evaluated to achieve hydrophilic wick structures. Evaluated treatments included chemical oxide growth with H2O2, oxide and sodium titanate growth with NaOH solution, and ultrasonic cleaning with specialty detergents that can remove hydrocarbon contaminants. The most durable hydrophilic surface conditions were obtained with the NaOH treatment. High performing wick deposition parameters were used to produce a full AM NiTi heat pipe, which was treated with NaOH solution to activate the wick. This heat pipe was operated on a test stand in the inverted configuration (upper evaporator and lower condenser), and demonstrated stable nearly isothermal operation for 150 hrs.

Thermal management↗

NASA NTRSDate not supplied

Additive Manufacturing and Experimental Characterization of Nickel-Titanium Shape-Memory Alloy Wick Structures and Heat Pipes for Spacecraft Thermal Control

Thermal management↗

DOE OSTI2022

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Sockeye Heat Pipe Analysis Code Verification and Validation

Some of the most promising microreactor designs currently under development utilize heat pipe technology to transfer heat from the reactor core to the secondary side heat exchanger, due to the technology’s compactness, efficiency, passivity, and reliability. Sockeye is an engineering-scale heat pipe tool developed under the Nuclear Energy Advanced Modeling and Simulation Program to be used for the design and safety analysis of microreactors. Sockeye’s core capability lies in a 1D, two-phase, compressible flow model, used to simulate the working fluid inside the heat pipe. Sockeye is built on the Multiphysics Object-Oriented Simulation Environment framework, which allows for seamless multiphysics coupling with other Nuclear Energy Advanced Modeling and Simulation tools and can thus be used in a full-scale simulation of a microreactor assembly, which can include hundreds of heat pipes. This paper presents some initial verification and validation assessments performed for Sockeye. We demonstrate good agreement between Sockeye’s numerical results and steady-state analytic solutions for velocity and pressure drop, with differences attributable to the underlying assumptions made by the analytic solutions. We also show that Sockeye reproduces analytic predictions of key operational limits, such as the capillary limit and the sonic limit. Finally, we compare Sockeye results to experimental data for the SAFE-30 heat pipe module test.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2017

Renin-Angiotensin System (RAS) in Hematopoietic Stem/Progenitor Cells (HS/PC) Predicts Vaso-Reparative Dysfunction and Progression of Diabetic Retinopathy (DR)

Purpose: We tested the hypothesis that loss of angiotensin converting enzyme 2 (ACE2) within diabetic HS/PCs (Hematopoietic Stem/Progenitor Cells) would be detrimental to HS/PC reparative function, and alter their ability to contribute to vascular remodeling in human subjects and rodent models of DR (Diabetic Retinopathy). Methods: Subjects (n52) were recruited as controls (n13) or diabetics (n39) with either no DR, mild non-proliferative DR (NPDR), moderate NPDR, severe NPDR or proliferative DR (PDR). Fluorescein angiograms were analyzed using Vessel Generation Analysis (VESGEN) software in a cohort of subjects. CD34+ HS/PCs were isolated from peripheral blood. RAS (Renin-Angiotensin System) gene expression and migration was measured. Diabetic ACE2 knockout (KO)C57BL6-Ins2 (Akita) mice at 3, 6 and 9 months of diabetes were compared to age-matched controls. Bone marrow HS/PC populations were analyzed by flow cytometry and migration and proliferation studies performed. Results: ACE2 gene expression in human CD34+ cells from diabetics without DR was increased compared to controls (p0.0437). Mas receptor mRNA was also increased in diabetics without DR, but reduced with the onset of NPDR (p0.0002), suggesting a loss of compensation. DR was associated with CD34+ cell migratory dysfunction. By VESGEN analysis, vessel density measured by several confirming parameters in early NPDR (n3) was greater than in normal retina (n6) in both arteries and veins, which suggests active retinal remodeling. ACE2KO-Akita and Akita cohorts showed reduced retinal thickness by OCT (Optical Coherence Tomography) at 9 months of diabetes. Absence of ACE2 in 9-month Akita mice led to an accelerated increase in acellular capillaries compared to diabetic alone. Electroretinogram (ERG) in ACE2KO-Akita mice resulted in persistent deterioration of the neural retina. Reparative function studies showed that ACE2KO exacerbated diabetes-induced impairment of LK (Low Potassium) cell migration and proliferative functions as early as 3-month of diabetes (p0.0019). Conclusions: Retinopathy and adverse vascular remodeling in subjects with diabetes was associated with a loss of the protective arm of RAS in HS/PCs. Loss of ACE2 exacerbated vascular dysfunction in diabetic mice.

retina↗

NASA NTRS1987

Spacelab Life Sciences flight experiments: an integrated approach to the study of cardiovascular deconditioning and orthostatic hypotension

The microgravity environment of spaceflight produces rapid cardiovascular changes which are adaptive and appropriate in that setting, but are associated with significant deconditioning and orthostatic hypotension on return to Earth's gravity. The rapidity with which these space flight induced changes appear and disappear provides an ideal model for studying the underlying pathophysiological mechanisms of deconditioning and orthostatic hypotension, regardless of etiology. Since significant deconditioning is seen after flights of very short duration, muscle atrophy due to inactivity plays, at most, a small role. These changes in circulatory control associated with cephalad fluid shifts, rather than inactivity per se, are probably more important factors. In order to test this hypothesis in a systematic way, a multidisciplinary approach which defines and integrates inputs and responses from a wide variety of circulatory sub-systems is required. The cardiovascular experiments selected for Spacelab Life Sciences flights 1 and 2 provide such an approach. Both human and animal models will be utilized. Pre- and post-flight characterization of the payload crew includes determination of maximal exercise capacity (bicycle ergometry), orthostatic tolerance (lower body negative pressure), alpha and beta adrenergic sensitivity (isoproterenol and phenylephrine infusions), baroreflex sensitivity (ECG-gated, stepwise changes in carotid artery transmural pressure with a pneumatic neck collar), and responses to a 24 h period of 5 deg head-down tilt. Measurements of cardiac output (CO2 and C2H2 rebreathing), cardiac chamber dimensions (phased-array 2-dimensional echocardiography), direct central venous pressure, leg volume (Thornton sock), limb blood flow and venous compliance (occlusion plethysmography), blood and plasma volumes, renal plasma flow and glomerular filtration rates, and various hormonal levels including catecholamines and atrial natriuretic factor will also be obtained. The central venous catheter will be inserted immediately pre-launch and monitored with heart rate and blood pressure in-flight until cardiac output, respiratory gas exchange and quantitative 2D echocardiography measurements can be performed. In-flight hemodynamic measurements will be repeated at rest and during submaximal exercise daily and also during maximal exercise midway through the flight to document the timecourse and extent of cardiovascular changes in the payload crew. Parallel studies are planned for the animals. In addition to measurements of right atrial and aortic pressures and cardiac output, a dorsal micro-circulatory chamber will allow determinations of changes in capillary and venular architecture and function in six of the rats. The techniques and findings from many of the SLS-1 and 2 supporting studies have already yielded significant information about circulatory regulation in patients with both hypo- and hypertension. The flight experiments themselves will provide new data to test the validity of both animal and human models currently used for simulating the fluid shifts of a micro-gravity environment. The development of effective countermeasures, not only for short and long duration space travellers, but also for Earth-bound medical patients can then be physiologically based on experimental data rather than anecdote.

Non-NASA Center↗

NASA NTRS1994

Performance model of a recirculating stack nickel hydrogen cell

A theoretical model of the nickel hydrogen battery cell has been utilized to describe the chemical and physical changes during charge and overcharge in a recirculating stack nickel hydrogen cell. In particular, the movement of gas and electrolyte have been examined as a function of the amount of electrolyte put into the cell stack during cell activation, and as a function of flooding in regions of the gas screen in this cell design. Additionally, a two-dimensional variation on this model has been utilized to describe the effects of non-uniform loading in the nickel-electrode on the movement of gas and electrolyte within the recirculating stack nickel hydrogen cell. The type of nonuniform loading that has been examined here is that associated with higher than average loading near the surface of the sintered nickel electrode, a condition present to some degree in many nickel electrodes made by electrochemical impregnation methods. The effects of high surface loading were examined primarily under conditions of overcharge, since the movement of gas and electrolyte in the overcharging condition was typically where the greatest effects of non-uniform loading were found. The results indicate that significant changes in the capillary forces between cell components occur as the percentage of free volume in the stack filled by electrolyte becomes very high. These changes create large gradients in gas-filled space and oxygen concentrations near the boundary between the separator and the hydrogen electrode when the electrolyte fill is much greater than about 95 percent of the stack free volume. At lower electrolyte fill levels, these gaseous and electrolyte gradients become less extreme, and shift through the separator towards the nickel electrode. Similarly, flooding of areas in the gas screen cause higher concentrations of oxygen gas to approach the platinum/hydrogen electrode that is opposite the back side of the nickel electrode. These results illustrate the need for appropriate pore size distributions, and the maintenance of both convective electrolyte and gas flow paths through the stack, if the recirculating stack nickel hydrogen cell design is to work properly.

Zimmerman, Albert H.↗

NASA NTRS1991

Cryogenic liquid-jet breakup in two-fluid atomizers

A two-fluid atomizer was used to study the breakup of liquid-nitrogen jets in nitrogen, argon, and helium atomizing gas flows. A scattered-light scanner particle sizing instrument previously developed at NASA Lewis Research Center was further developed and used to determine characteristic drop diameters for the cryogenic sprays. In the breakup regime of aerodynamic-stripping, i.e., sonic-velocity conditions, the following correlation of the reciprocal Sauter mean diameter, D(sub 32)exp -1, with the atomizing-gas flowrate, W(g), was obtained: D(sub 32)exp -1 = k(sub c)(W(g)exp 1.33), where k(sub c) is a proportionality constant evaluated for each atomizing gas. Values of k(sub c) = 120, 220, and 1100 were obtained for argon, nitrogen, and helium gasflows, respectively. The reciprocal Sauter mean diameter and gas flowrate have the units of 1/cm and g/sec, respectively. In the regime of capillary-wave breakup, or subsonic conditions, it was found that D(sub 32)exp -1 = k(g)(W(g)exp 0.75), where k = 270, 390, and 880 for argon, nitrogen, and helium gasflows, respectively.

Ingebo, Robert D.↗

DOE OSTI2021

Represent precipitation-induced geological hazards in Earth system models using artificial intelligence

Precipitation-induced geological hazards, such as debris flow, landslides, mudflow and rockfalls (hereafter referred to as landslides), pose serious threats to public safety in many areas through the world. As residential properties and infrastructure in the US have increasingly expanded into landslide-prone areas and the drivers of landslides (e.g., wildfires and hurricanes) are predicted to intensify under climate warming, losses and fatalities from landslides are likely to increase in the future. Although our understanding of geoenvironmental factors and mechanisms contributing to landslides has greatly improved, only moderate progress has been made in predicting landslides out of well-studied watersheds. Furthermore, explicitly representing various landslide-related processes in Earth system models (ESMs), from the buckling of local bearing elements in granular materials, to frictional sliding between grains, formation of microcracks in the soil matrix, rupture of capillary bridges, or breakage of plant roots3 is still very unlikely within the next decade, even with the help of exascale computers.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2016

Demonstration of Nautilus Centripetal Capillary Condenser Technology

This paper describes the results of a proof of concept effort for development of a Nautilus Centripetal Capillary Condenser (NCCC or NC3) used for microgravity compatible water recovery from moist air with integral passive phase separation. Removal of liquid condensate from the air stream exiting a condenser is readily performed here on Earth. In order to perform this function in space however, without gravity or mechanical action, other tactics including utilization of inertial, drag and capillary forces are required. Within the NC3, liquid water forms via condensation on cold condenser surfaces as humid air passes along multiple spiral channels, each in its own plane, all together forming a stacked plate assembly. Non-mechanical inertial forces are employed to transfer condensate, as it forms, via centripetal action to the outer perimeter of each channel. A V-shaped groove, constructed on this outer edge of the spiral channel, increases local capillary forces thereby retaining the liquid. Air drag then pulls the liquid along to a collection region near the center of the device. Dry air produced by each parallel spiral channel is combined in a common orthogonal, out-of-plane conduit passing down the axial center of the stacked device. Similarly, the parallel condensate streams are combined and removed from the condenser/separator through yet another out-of-plane axial conduit. NC3 is an integration of conventional finned condenser operation, combined with static phase separation and capillary transport phenomena. A Mars' transit mission would be a logical application for this technology where gravity is absent and the use of vibrating, energy-intensive, motor-driven centrifugal separators is undesired. Here a vapor stream from either the Heat Melt Compactor or the Carbon dioxide Reduction Assembly, for example, would be dried to a dew point of 10 deg using a passive NC3 condenser/separator with the precious water condensate recycled to the water bus.

Separator↗

NASA NTRS1999

Pulsed high energy synthesis of fine metal powders

Repetitively pulsed plasma jets generated by a capillary arc discharge at high stagnation pressure (>15,000 psi) and high temperature (>10,000 K) are utilized to produce 0.1-10 .mu.m sized metal powders and decrease cost of production. The plasma jets impact and atomize melt materials to form the fine powders. The melt can originate from a conventional melt stream or from a pulsed arc between two electrodes. Gas streams used in conventional gas atomization are replaced with much higher momentum flux plasma jets. Delivering strong incident shocks aids in primary disintegration of the molten material. A series of short duration, high pressure plasma pulses fragment the molten material. The pulses introduce sharp velocity gradients in the molten material which disintegrates into fine particles. The plasma pulses have peak pressures of approximately one kilobar. The high pressures improve the efficiency of disintegration. High gas flow velocities and pressures are achieved without reduction in gas density. Repetitively pulsed plasma jets will produce powders with lower mean size and narrower size distribution than conventional atomization techniques.

Witherspoon, F. Douglas↗

DOE OSTI2025

A theoretical study of the upper bound of surface elevation variance in the Phillips initial stage during wind-wave generation

The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

58 GEOSCIENCES↗

DOE OSTI2021

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

36 MATERIALS SCIENCE↗