Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

At least 145 records · Page 8

Tetragonal Lysozyme Nucleation and Crystal Growth: The Role of the Solution Phase

Lysozyme, and most particularly the tetragonal form of the protein, has become the default standard protein for use in macromolecule crystal nucleation and growth studies. There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests this proteins crystal nucleation and growth is by addition of associated species that are preformed by standard reversible concentration-driven self association processes in the bulk solution. The evidence includes high resolution AFM studies of the surface packing and of growth unit size at incorporation, fluorescence resonance energy transfer measurements of intermolecular distances in dilute solution, dialysis kinetics, and modeling of the growth rate data. We have developed a selfassociation model for the proteins crystal nucleation and growth. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. Further, it indicates that nucleation and crystal growth are not distinct mechanistically, but identical, with the primary difference being the probability that the particle will continue to grow or dissolve. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small (M.W. = 14,400), a structured octamer in the 4(sub 3) helix configuration (the proposed average sized growth unit) would have a M.W. = 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm commonly attach to the crystal. These measurements were made at approximately saturation conditions, and they reflect the sizes of species that both added or desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. Flow-driven effects resulting in misorientation upon addition to and incorporation into the crystal need only be a small fraction of a percentage to significantly affect the resulting crystal. One Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species will be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc L.↗

An Alternative Hypothesis for How Microgravity Improves Macromolecular Crystal Quality

There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests that for many proteins crystal nucleation and growth is by addition of associated species that are preformed by reversible concentration-driven self association processes in the bulk solution. We have developed a self-association model for the crystal nucleation and growth of the protein chicken egg lysozyme. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small an octamer in the 43 helix configuration (the proposed average sized growth unit) would have a M.W. approx. 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm (a 24-mer) commonly attach to the crystal. AFM results from Weichmann et al. (Ultramicroscopy 86, 159-166, 2001) suggest that associated species of up to 40-mers in size add to the (101) faces. These measurements reflect the sizes of units that both added and desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. On Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species is more likely to be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc↗

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interfacial Area Transport through Vertical-Downward Elbows in Air-water Two-phase Flow

In the study of two-phase flow, models are often developed by the data acquired in a test section with fixed orientation and with no flow restrictions. In view of practical nuclear reactor systems, however, two-phase flow through flow restrictions is important. Flow restrictions and orientations can alter both global and local two-phase flow parameters such as: cause additional pressure loss; alter the flow regime transition; change bubble distributions; affect the bubble velocity and the bubble interaction mechanisms. However, there is insufficient database and predictive models that can account for such effects. In view of above, the objectives of the current study are: (1) to establish database for local and global two-phase flow parameters downstream of a horizontal-to-vertical glass elbow with a non-dimensional radius of curvature, R{sub c}/D, of three, (2) to develop predictive models in the elbow region, and (3) to apply the developed model to the interfacial area transport equation (IATE) and evaluate it in elbow regions. (authors)

42 ENGINEERING↗

Stability of oscillatory two phase Couette flow

We investigate the stability of two phase Couette flow of different liquids bounded between plane parallel plates. One of the plates has a time dependent velocity in its own plane, which is composed of a constant steady part and a time harmonic component. In the absence of time harmonic modulations, the flow can be unstable to an interfacial instability if the viscosities are different and the more viscous fluid occupies the thinner of the two layers. Using Floquet theory, we show analytically in the limit of long waves, that time periodic modulations in the basic flow can have a significant influence on flow stability. In particular, flows which are otherwise unstable for extensive ranges of viscosity ratios, can be stabilized completely by the inclusion of background modulations, a finding that can have useful consequences in many practical applications.

Coward, Adrian V.↗

Structural Gating Enhances Long-Distance Light-Driven Interfacial Electron Transfer

Structural gating provides a molecular means to transfer electrons preferentially in one desired vectorial direction, a behavior needed for applications in artificial photosynthesis. At the interfaces utilized herein, visible-light absorption by a transition metal complex opens a “structural gate” by planarization of otherwise rotating phenyl rings in p-phenylene ethynylene (PE) bridge units. Planarization provides a conjugated pathway for electron flow toward a conductive oxide surface. Interfacial electron transfer to the oxide restores rotation and closes the gate to the unwanted recombination reaction. This structural gating results in nearly quantitative long-distance (>20 Å) interfacial electron transfer that occurs ~1000 times faster than transfer in the opposite direction. A comparative kinetic study of these complexes with those that contain ionic bridge units, without gating function, as a function of the applied potential and hence –ΔG° provided a physical basis for the structural gating. A small distance-dependent reorganization energy with weak electronic coupling underlies the success of this gate that enables efficient long-distance electron transfer and slow recombination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement of the Surface Dilatational Viscosity of an Insoluble Surfactant Monolayer at the Air/Water Interface Using a Pendant Drop Apparatus

When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Lorenzo, Jose↗

Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with Jacobian-free Newton–Krylov method

This study presents a numerical investigation on using the Jacobian-free Newton–Krylov (JFNK) method to solve the two-phase flow four-equation drift flux model with realistic constitutive correlations (‘closure models’). The drift flux model is based on Isshi and his collaborators’ work. Additional constitutive correlations for vertical channel flow, such as two-phase flow pressure drop, flow regime map, wall boiling and interfacial heat transfer models, were taken from the RELAP5-3D Code Manual and included to complete the model. The staggered grid finite volume method and fully implicit backward Euler method was used for the spatial discretization and time integration schemes, respectively. The Jacobian-free Newton–Krylov method shows no difficulty in solving the two-phase flow drift flux model with a discrete flow regime map. In addition to the Jacobian-free approach, the preconditioning matrix is obtained by using the default finite differencing method provided in the PETSc package, and consequently the labor-intensive implementation of complex analytical Jacobian matrix is avoided. Extensive and successful numerical verification and validation have been performed to prove the correct implementation of the models and methods. Code-to-code comparison with RELAP5-3D has further demonstrated the successful implementation of the drift flux model.

42 ENGINEERING↗

Current–voltage characteristics of organic heterostructure devices with insulating spacer layers

The dark current density in donor/acceptor organic planar heterostructure devices at a given forward voltage bias can either increase or decrease when an insulating spacer layer is added between the donor and acceptor layers. The dominant current flow process in these systems involves the formation and subsequent recombination of interfacial exciplex states. If the exciplex recombination rate limits current flow, an insulating interface layer decreases the dark current. However, if the exciplex formation rate limits the current, an insulating interface layer may increase the dark current. As a result, we present a device model to describe this behavior, and we discuss relevant experimental data.

25 ENERGY STORAGE↗

Thermocapillary convection in two immiscible liquid layers with free surface

Thermocapillary convection is studied in two immiscible liquid layers with one free surface, one liquid/liquid interface, and differential heating applied parallel to the interfaces. An analytical solution is introduced for infinite horizontal layers. The defining parameter for the flow pattern is lambda, the ratio of the temperature coefficient of the interfacial tension to that of the surface tension. Four different flow patterns exist under zero gravity conditions. 'Halt' conditions which halt the fluid motion in the lower encapsulated liquid layer have been found. A numerical experiment is carried out to study effects of vertical end walls on the double layer convection in a 2D cavity. The halt condition obtained from the analytical study is found to be valid in the limit of small Reynolds numbers. The flow in the encapsulated liquid layer can be suppressed substantially.

Doi, Takao↗

Simulation and analysis of slip flow of water at hydrophobic silica surfaces of nanometer slit pores

Water flow in nanopores is controlled fundamentally by pore wettability and interfacial water features. In the research herein reported, Molecular dynamics (MD) was used to simulate water flow in nanometer slit pores of silica surfaces with variable wetting characteristics. Compared to Couette flow, the simulated slit-pores of the talc (001) and octadecyltrichlorosilane (OTS) monolayer surfaces exhibited water flow enhancement and a positive slip length at the interface, whereas water flow velocity at the quartz (001) and amorphous silica surfaces was zero, adhering to the no slip condition. More intensive interactions between the interfacial water and hydrophilic silica surfaces were confirmed by MD simulation and sum-frequency vibrational spectroscopy (SFVS) results. A water exclusion zone of approximately 3 Å was found at the hydrophobic silica surfaces by quantitative MD interfacial water analysis and validated by SFVS. This exclusion zone dimension was similar in magnitude to the slip length (4–5 Å) as determined from simulated water flow in hydrophobic silica slit pores. The simulated slip lengths at the hydrophobic surfaces (0.2–1 nm) were significantly less than the experimental slip lengths (>40 nm) reported in the literature. In the case of the experimental results, this difference appears to be due to nanobubbles expected to be present at hydrophobic surfaces.

42 ENGINEERING↗

Stochastic level-set variational implicit-solvent approach to solute-solvent interfacial fluctuations

Recent years have seen the initial success of a variational implicit-solvent model (VISM), implemented with a robust level-set method, in capturing efficiently different hydration states and providing quantitatively good estimation of solvation free energies of biomolecules. The level-set minimization of the VISM solvation free-energy functional of all possible solute-solvent interfaces or dielectric boundaries predicts an equilibrium biomolecular conformation that is often close to an initial guess. In this work, we develop a theory in the form of Langevin geometrical flow to incorporate solute-solvent interfacial fluctuations into the VISM. Such fluctuations are crucial to biomolecular conformational changes and binding process. We also develop a stochastic level-set method to numerically implement such a theory. We describe the interfacial fluctuation through the “normal velocity” that is the solute-solvent interfacial force, derive the corresponding stochastic level-set equation in the sense of Stratonovich so that the surface representation is independent of the choice of implicit function, and develop numerical techniques for solving such an equation and processing the numerical data. We apply our computational method to study the dewetting transition in the system of two hydrophobic plates and a hydrophobic cavity of a synthetic host molecule cucurbit[7]uril. Numerical simulations demonstrate that our approach can describe an underlying system jumping out of a local minimum of the free-energy functional and can capture dewetting transitions of hydrophobic systems. In the case of two hydrophobic plates, we find that the wavelength of interfacial fluctuations has a strong influence to the dewetting transition. In addition, we find that the estimated energy barrier of the dewetting transition scales quadratically with the inter-plate distance, agreeing well with existing studies of molecular dynamics simulations. Our work is a first step toward the inclusion of fluctuations into the VISM and understanding the impact of interfacial fluctuations on biomolecular solvation with an implicit-solvent approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent motions and time scales that control heat and mass transfer at wind-swept water surfaces

Forecast of the heat and chemical budgets of lakes, rivers, and oceans requires improved predictive understanding of air-water interfacial transfer coefficients. Here in this paper we present laboratory observations of the coherent motions that occupy the air-water interface at wind speeds (U 10 ) 1.1–8.9 m/s. Spatiotemporal near-surface velocity data and interfacial renewal data are made available by a novel flow tracer method. The relative activity, velocity scales, and time scales of the various coherent interfacial motions are measured, namely for Langmuir circulations, streamwise streaks, nonbreaking wind waves, parasitic capillary waves, nonturbulent breaking wind waves, and turbulence-generating breaking wind waves. Breaking waves exhibit a sudden jump in streamwise interfacial velocity wherein the velocity jumps up to exceed the wave celerity and destroys nearby parasitic capillary waves. Four distinct hydrodynamic regimes are found to exist between U 10 =0 and 8.9 m/s, each with a unique population balance of the various coherent motions. The velocity scales, time scales, and population balance of the different coherent motions are input to a first-principles gas transfer model to explain the waterside transfer coefficient (k w ) as well as experimental patterns of temperature and gas concentration. The model mixes concepts from surface renewal and divergence theories and requires surface divergence strength (β), the Lagrangian residence time inside the upwelling zone (t Lu ), and the total lifetime of new interface before it is downwelled (t LT ). The model’s output agrees with time-averaged measurements k w , patterns of temperature in infrared photographs, and spatial patterns of gas concentration and k w from direct numerical simulations. Finally, several nondimensional parameters, e.g. βt Lu and τ s t LT where ss is the interfacial shear rate, determine the effectiveness of a particular type of coherent motion for affecting k w .

16 TIDAL AND WAVE POWER↗

Impact of KOH on the interfacial precipitation rates of C-S-H during the early hydration of C3S

In this work we investigated the influence of alkalis on the hydration of two triclinic samples of C{sub 3}S. One pure C{sub 3}S and one Al-doped C{sub 3}S (0.5 wt.‐% Al{sub 2}O{sub 3}). We added differently concentrated KOH solutions representing alkalis in the pore solution, which lead to an increase of the interfacial precipitation rates of C-S-H during the hydration of C{sub 3}S compared to the hydration without KOH. Isothermal calorimetry was conducted for the determination of the rates of reaction during the hydration process and ICP-OES measurements for the determination of the pore solution composition at certain times of hydration. The increasing effect of KOH solution on the calculated interfacial precipitation rates from pore solution data can explain the enhancement of the rates of reaction recorded by heat flow calorimetry during the early hydration. Therefore, the interfacial precipitation rates of C-S-H can accurately describe the acceleration effect of alkalis on the early hydration of C{sub 3}S.

36 MATERIALS SCIENCE↗

Investigation of Body Force Effects on Flow Boiling Critical Heat Flux

The bubble coalescence and interfacial instabilities that are important to modeling critical heat flux (CHF) in reduced-gravity systems can be sensitive to even minute body forces. Understanding these complex phenomena is vital to the design and safe implementation of two-phase thermal management loops proposed for space and planetary-based thermal systems. While reduced gravity conditions cannot be accurately simulated in 1g ground-based experiments, such experiments can help isolate the effects of the various forces (body force, surface tension force and inertia) which influence flow boiling CHF. In this project, the effects of the component of body force perpendicular to a heated wall were examined by conducting 1g flow boiling experiments at different orientations. FC-72 liquid was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface at conditions approaching CHF. High-speed video imaging was employed to capture dominant CHF mechanisms. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed great sensitivity to orientation for flow velocities below 0.2 m/s, where very small CHF values where measured, especially with downflow and downward-facing heated wall orientations. High flow velocities dampened the effects of orientation considerably. Figure I shows representative images for the different CHF regimes. The Wavy Vapor Layer regime was dominant for all high velocities and most orientations, while all other regimes were encountered at low velocities, in the downflow and/or downward-facing heated wall orientations. The Interfacial Lift-off model was modified to predict the effects of orientation on CHF for the dominant Wavy Vapor Layer regime. The photographic study captured a fairly continuous wavy vapor layer travelling along the heated wall while permitting liquid contact only in wetting fronts, located in the troughs of the interfacial waves. CHF commenced when wetting fronts near the outlet were lifted off the wall. The Interfacial Lift-off model is shown to be an effective tool for predicting the effects of body force on CHF at high velocities.

Zhang, Hui↗

Liquid-phase turbulence measurements in air-water two-phase flows using particle image velocimetry

Liquid-phase turbulence plays a vital role in determining various gas-liquid two-phase flow parameters, such as void fraction distribution, bubble morphology, bubble-bubble interactions, and interfacial area concentration. In this study, a two-phase flow database including both the gas- and liquid-phases measurements was developed focusing on three bubbly flow conditions in an air-water two-phase flow loop with a vertical one-inch diameter circular pipe test section. A particle image velocimetry (PIV) system integrating an optical phase separation method, i.e., the planar laser-induced fluorescence (PLIF) technique using fluorescent particles and optical filtration, was applied to measure the liquid-phase turbulence information, including the time-averaged velocity, Reynolds stress, and turbulent kinetic energy for the liquid phase. The PIV measurements were taken at three ports along the test section at 14.5, 51.5, and 88.5 pipe inner diameters downstream of a bubble injector. In addition, a double-sensor conductivity probe was used to measure radial distributions of the local time-averaged void fraction and gas velocity. The measured liquid-phase turbulence was used to benchmark Sato's turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. Lastly, in the two bubbly flows tested that have low void fractions being less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction.

42 ENGINEERING↗