Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial velocity”

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 37 records · Page 2

Liquid-Phase Heat-Release Rates of the Systems Hydrazine-Nitric Acid and Unsymmetrical Dimethylhydrazine-Nitric Acid

The initial rates of heat release produced by the reactions of hydrazine and unsymmetrical dimethylhydrazine with nitric acid were determined in a bomb calorimeter under conditions of forced mixing. Fuel-oxidant weight ratio and injection velocity were varied. The rate of heat release apparently depended on the interfacial area between the propellants. Above a narrow range of injection velocities representing a critical amount of interfacial area, the rates reached a maximum and were almost constant with injection velocity. The maximum rate for hydrazine was about 70 percent greater than that for unsymmetrical dimethylhydrazine. The total heat released did not vary with mixture ratio over the range studied.

Somogyi, Dezso↗

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗

High-resolution measurements of microphysics and entrainment in marine stratocumulus clouds

High-resolution measurements from the Airborne Cloud-Turbulence Observation System (ACTOS) during the Azores Stratocumulus Measurements of Radiation, Turbulence and Aerosols (ACORES) campaign are analysed for an investigation of the vertical profiles of microphysical properties and entrainment velocity ( W e ) in marine stratocumulus clouds. Here, the vertical profiles show the transition from the cloudy layer to free troposphere with nearly linear profiles of total water mixing ratio, liquid water potential temperature and virtual potential temperature, but the thickness of entrainment interfacial layer varies significantly. Sharp transitions of cloud microphysical and optical properties within a single horizontal flight leg are found in one stratocumulus cloud system. They seem to be related to the local environmental conditions, such as the wind shear and turbulent dissipation rate. W e values estimated by three methods show consistent tendencies in general and are clearly related to the local environmental conditions, such as vertical shear of the horizontal wind and turbulence intensity. However, the magnitudes of W e values differ by up to two orders of magnitude depending on the methods, which suggests that the estimation of W e from in situ measurements is still a challenge. Analysis of the microphysical response to entrainment suggests that inhomogeneous mixing occurs dominantly. On the other hand, the analysis results for the clouds under more humid conditions indicate a higher likelihood of homogeneous mixing.

54 ENVIRONMENTAL SCIENCES↗

Interaction between the atmospheric and oceanic boundary layers

The two-layer system of an atmosphere over water bodies is reduced to a single-layer problem. Values of the interfacial quantities, such as the friction velocity, the surface velocity, the angles, alpha and beta, between the surface shear stress and the geostrophic wind velocity and the surface wind velocity, respectively, and the surface roughness, all of which depend upon external parameters, such as the geostrophic wind and stratifications, are obtained. The geostrophic drag coefficient, the geostrophic wind coefficient, and the angles alpha, and beta, of the turbulent flow at the sea-air interface are functions of a dimensionless number, mfG/kg, with S sub 1 and S sub 2 as two free stratification parameters. The surface roughness is uniquely determined from the geostrophic wind rather than from the wind profile in the boundary layer.

Yeh, G.-T.↗

The roles of fluid motion and other transport phenomena in the morphology of materials

Two crystallization problems were studied: the growth of protein crystals, in particular the influence of colloidal forces and convection, and the influence of interface resistance on the growth of dendritic crystals. The protein study involved both experimental and theoretical work; the work of dendrites was entirely theoretical. In the study of protein crystallization, experiments were carried out where crystals were grown in the presence and absence of natural convection. No evidence was found that convection retards crystal growth. The theoretical study focused on the influence of colloidal forces (electrostatic and London-van der Waals) on the interaction between a protein molecule and a flat crystal surface. It was shown that the interaction is extremely sensitive to colloidal forces and that electrostatic interactions play a strong role in deciding whether or not a molecule will find a favorable site for adsorption. In the study of dendritic growth, the role of an interfacial resistance on the selection processes was examined. Using a computational scheme, it was found that the selected velocity is strongly dependent on the magnitude of the interfacial resistance to heat transfer. This is a possible explanation for discrepancies between the theoretical and experimental results on succinonitrile.

Saville, D. A.↗

Effects of crystal-melt interfacial energy anisotropy on dendritic morphology and growth kinetics

Morphological and kinetic studies of succinonitrile, a BCC crystal with a low (0.5 percent) anisotropy and pivalic acid, and FCC crystal with relatively large (5 percent) anisotropy in solid-liquid interfacial energy, show clearly that anisotropy in the solid-liquid interfacial energy does not affect the tip radius-velocity relationship, but has a profound influence on the tip region and the rate of amplification of branching waves. Anisotropy of the solid-liquid interfacial energy may be one of the key factors by which the microstructural characteristics of cast structures reflect individual material behavior, especially crystal symmetry.

Glicksman, M. E.↗

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗

Review of critical flow rate, propagation of pressure pulse, and sonic velocity in two-phase media

For single-phase media, the critical discharge velocity, the sonic velocity, and the pressure pulse propagation velocity can be expressed in the same form by assuming isentropic, equilibria processes. In two-phase mixtures, the same concept is not valid due to the existence of interfacial transports of momentum, heat, and mass. Thus, the three velocities should be treated differently and separately for each particular condition, taking into account the various transport processes involved under that condition. Various attempts are reviewed to predict the critical discharge rate or the propagation velocities by considering slip ratio (momentum change), evaporation (mass and heat transport), flow pattern, etc. Experimental data were compared with predictions based on various theorems. The importance is stressed of the time required to achieve equilibrium as compared with the time available during the process, for example, of passing a pressure pulse.

Hsu, Y.↗

Tuning Kinetic Inductance with Doping in Superconducting Electron Gases at the KTaO 3 (111) Interface

We used coplanar waveguide resonators fabricated from the two-dimensional superconductor formed at KTaO 3 (111) interfaces to characterize the superconductor’s sheet kinetic inductance L K/□ . Upon varying the carrier density, in samples with T c ranging from 0.62 to 1.81 K, we can tune L K/□ between 1.88 and 7.42 nH/□ at the lowest temperatures, exceeding the highest values reported for granular aluminum films. The temperature dependence of L K/□ is consistent with a superconducting gap without nodes. The high L K/□ of superconducting KTaO 3 (111) interfacial electron gases combined with the high dielectric constant of KTaO 3 results in resonators with exceedingly low phase velocities (“slow light”) and small mode volumes. In addition, superconducting KTaO 3 (111) interfacial electron gases are robust in magnetic fields well above a Tesla. In conclusion, these features can enable unique device applications in superconducting electronics and quantum information science.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of particle velocity fluctuations on the inertia coupling in two-phase flow

Consistent forms for the interfacial force, the interfacial pressure, the Reynolds stresses and the particle stress have been derived for the inviscid, irrotational incompressible flow of fluid in a dilute suspension of spheres. The particles are assumed to have a velocity distribution, giving rise to an effective pressure and stress in the particle phase. The velocity fluctuations also contribute in the fluid Reynolds stress and in the (elastic) stress field inside the spheres. The relation of these constitutive equations to the force on an individual sphere is discussed.

Drew, Donald A.↗

Use of the Burton-prim-slichter Equation at High Growth Rates

The following must be assumed in order to derive the classic Burton-Prim-Slichter equation for segregation during crystal growth: isobaric, isothermal, isopotential or uncharged species, binary, planar interface, steady state, constant diffusion coefficient in fluid, no diffusion in crystal, no lateral convection within fluid film at interface with complete mixing beyond (stagnant film model), and either density, total concentration or partial molar volumes constant in the fluid phase. In addition, the effective distribution coefficient and the interfacial distribution coefficient must be defined properly. The velocity in the equation is the growth rate times a factor correcting for the difference in volumetric properties between crystal and fluid. The stagnant film thickness is found to be a function of freezing rate, with the precise functionality depending on the type of stirring occurring in the growth fluid.

Wilcox, W. R.↗

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peclet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peelet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗

Investigation of Critical Heat Flux in Reduced Gravity Using Photomicrographic Techniques

Experiments were performed to examine the effects of body force on flow boiling critical heat flux (CHF). FC-72 was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface just prior to CHF. High-speed video imaging techniques were used to identify dominant CHF mechanisms corresponding to different flow orientations and liquid velocities. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed significant sensitivity to orientation for flow velocities below 0.2 m/s, where extremely low CHF values where measured, especially with downward-facing heated wall and downflow orientations. High flow velocities dampened the effects of orientation considerably. The CHF data were used to assess the suitability of previous CHF models and correlations. It is shown the Interfacial Lift-off Model is very effective at predicting CHF for high velocities at all orientations. The flooding limit, on the other hand, is useful at estimating CHF at low velocities and for downflow orientations. A new method consisting of three dimensionless criteria is developed for determining the minimum flow velocity required to overcome body force effects on near-saturated flow boiling CHF. Vertical upflow boiling experiments were performed in pursuit of identifying the trigger mechanism for subcooled flow boiling CHF. While virtually all prior studies on flow boiling CHF concern the prediction or measurement of conditions that lead to CHF, this study was focused on events that take place during the CHF transient. High-speed video imaging and photomicrographic techniques were used to record the transient behavior of interfacial features from the last steady-state power level before CHF until the moment of power cut-off following CHF. The video records show the development of a wavy vapor layer which propagates along the heated wall, permitting cooling prior to CHF only in wetting fronts corresponding to the wave troughs. Image analysis software was developed to estimate void fraction from the individual video images. The void fraction records for subcooled flow boiling show the CHF transient is accompanied by gradual lift-off of wetting fronts culminating in some maximum vapor layer mean thickness, following which the vapor layer begins to thin down as the transition to film boiling ensues. This study proves the Interfacial Lift-off Model, which has been validated for near-saturated flow boiling CHF, is equally valid for subcooled conditions.

Mudawar, Issam↗

The coupling of interfacial instabilities and the stabilization of two-layer annular flows

In this paper the stability of annular pressure-driven parallel flows of two liquids sandwiching a free cylindrical interface is considered. For small to moderate Reynolds numbers, the interface is susceptible to capillary and interfacial wave instabilities, the latter instability caused by a jump in viscosity at the interface. It is shown that favorable velocity profiles in both liquids may stabilize capillary breakup of the interface and suppress the axisymmetric interfacial wave instability. A long-wave analysis leads to the physical mechanism responsible for stabilization of capillary breakup. This physical mechanism is a generalization of that by which capillary breakup is stabilized by interfacial shear in an annular film of a single liquid. Stabilization of intermediate wavelengths is studied with a mechanical energy analysis, which leads to a description of the energetic processes at work.

Dijkstra, Henk A.↗

Numerical and Physical Simulation of the Fluid Flow in a Beam Blank Mold Fed Through the Center of the Web

Fluid flow inside a beam blank mold fed through a three-port SEN (two lateral ports and one bottom port), positioned at the center of the mold, has been investigated. Literature survey shows that this kind of configuration is not frequently used. It had been shown that it is possible to get a symmetrical flow with this configuration. The downward inclination of the ports should influence the transient fluid flow. The greater the downward inclination, the greater is the instability of the jet and of the interface between the immiscible fluids that simulate the slag/metal behavior. The flow characterization was made by dye dispersion, PIV technique, measurement of the meniscus oscillation using ultrasonic sensors, and CFD simulations. The slag/metal behavior was simulated using water and oils with different physical properties. The results from PIV as well as from observations of the water–oil interface have been used to validate the CFD simulations. Oils with density close to water resulted in more instability at the interface with entrainment starting from fluid flow rate of 125 L/min (equivalent to casting speed of 1 m/min). Decreasing the viscosity of slag (oils) and increasing the casting velocity (water flow rate) result in reduction of the interfacial stability.

36 MATERIALS SCIENCE↗

Thermocapillary motion of deformable drops

The thermocapillary motion of initially spherical drops/bubbles driven by a constant temperature gradient in an unbounded liquid medium is simulated numerically. Effects of convection of momentum and energy, as well as shape deformations, are addressed. The method used is based on interface tracking on a base cartesian grid, and uses a smeared color or indicator function for the determination of the surface topology. Quad-tree adaptive refinement of the cartesian grid is implemented to enhance the fidelity of the surface tracking. It is shown that convection of energy results in a slowing of the drop, as the isotherms get wrapped around the front of the drop. Shape deformation resulting from inertial effects affect the migration velocity. The physical results obtained are in agreement with the existing literature. Furthermore, remarks are made on the sensitivity of the calculated solutions to the smearing of the fluid properties. Analysis and simulations show that the migration velocity depends very strongly on the smearing of the interfacial force whereas it is rather insensitive to the smearing of other properties, hence the adaptive grid.

Haj-Hariri, Hossein↗