Incompressible Flows Free Surfaces
Energy Science and Technology Software Center (ESTSC)
1992-02-01
NASA-VOF3D is a three-dimensional, transient, free surface, incompressible fluid dynamics program. It is specifically designed to calculate confined flows in a low gravity environment in which surface physics must be accurately treated. It allows multiple free surfaces with surface tension and wall adhesion and includes a partial cell treatment that allows curved boundaries and internal obstacles. Variable mesh spacing is permitted in all three coordinate directions. Boundary conditions available are rigid free-slip wall, rigid no-slipmore » wall, continuative, periodic, and specified pressure outflow boundary.« less
Retention/Diffusivity Studies in Free-Surface Flowing Liquid Lithium
R.A. Stubbers; G.H. Miley; M. Nieto; W. Olczak; D.N. Ruzic; A. Hassanein
2004-12-14
FLIRE was designed to measure the hydrogen and helium retention and diffusivity in a flowing stream of liquid lithium, and it has accomplished these goals. Retention coefficients for helium in the flowing liquid stream were 0.1-2% for flow speeds of 44 cm/s and implantation energies between 500 and 2000 eV. The energy dependence of retention is linear for the energy range considered, as expected, and the dependence of retention on flow velocity fits the expected square-root of flow speed dependence. Estimates of the helium diffusion coefficient in the flowing lithium stream were {approx} 4 x 10{sup -7} cm{sup 2}/s, and are independent of implantation energy. This value is much lower than expected, which could be due to several factors, such as mixing, bubble formation or surface film formation. In the case of hydrogen, long term retention and release mechanisms are of greatest importance, since this relates to tritium inventory in flowing lithium PFCs for fusion applications. The amount of hydride formation was measured for flowing lithium exposed to neutral deuterium gas. Thermal desorption spectroscopy (TDS) measurements indicate that the hydride concentration was between 0.1 and 0.2% over a wide range of pressures (6.5 x 10{sup -5} to 1 Torr). This result implies that the deuterium absorption rate is limited by the surface dissociation rate, since deuterium (hydrogen/tritium) is absorbed in its atomic form, not its molecular form.
Computing Incompressible Flows With Free Surfaces
NASA Technical Reports Server (NTRS)
Kothe, D.
1994-01-01
RIPPLE computer program models transient, two-dimensional flows of incompressible fluids with surface tension on free surfaces of general shape. Surface tension modeled as volume force derived from continuum-surface-force model, giving RIPPLE both robustness and accuracy in modeling surface-tension effects at free surface. Also models wall adhesion effects. Written in FORTRAN 77.
Natural Elements Method for Free Surface Flow
NASA Astrophysics Data System (ADS)
Darbani, M.; Ouahsine, A.; Villon, P.
2009-09-01
The Natural Element Method (NEM) is used to simulate a 2D shallow water flow in presence of free surface and a varying bathymetry. This meshless method used a fully Lagrangian formulation and natural neighbors, which remain a very striking problem related the boundary conditions. The method was succefully used to simulate dam-break flows by solving the fully nonlinear Shallow Water Equations (SWE) and by using an implicit scheme under a transient flow and the Coriolis effect.
Finite element solutions of free surface flows
NASA Technical Reports Server (NTRS)
Zarda, P. R.; Marcus, M. S.
1977-01-01
A procedure is presented for using NASTRAN to determine the flow field about arbitrarily shaped bodies in the presence of a free surface. The fundamental unknown of the problem is the velocity potential which must satisfy Laplace's equation in the fluid region. Boundary conditions on the free surface may involve second order derivatives in space and time. In cases involving infinite domains either a tractable radiation condition is applied at a truncated boundary or a series expansion is used and matched to the local finite elements. Solutions are presented for harmonic, transient, and steady state problems and compared to either exact solutions or other numerical solutions.
Free-surface and coating flows
NASA Astrophysics Data System (ADS)
Ashmore, Jacqueline
In thus thesis we present results on four problems which fall into the class of free-surface flows. Two of the problems relate to coating flows, in which a thin film of fluid is deposited on a moving substrate. The third problem concerns the motion of a sphere down an inclined plane that is coated with a thin film of viscous fluid. Finally, we study the Rayleigh capillary instability of a cylindrical fluid thread surrounded by a second viscous fluid in a rotating frame. One useful tool for analysis of coating flows is the method of matched asymptotics, which enables us to obtain a solution to the nonlinear third-order differential equation that describes the free-surface shape in two of the problems we study. The first coating flow problem we consider is the interface shape of fluid inside a horizontal cylinder rotating about its axis with a small fraction of its volume filled with viscous Newtonian fluid. By accounting for surface tension effects we find a new axially uniform steady solution valid at low rotation rates. Predictions for the scalings based on matched asymptotic methods are confirmed numerically. Secondly, theoretical predictions of the film thickness that coats a substrate when it is withdrawn from a bath of viscoelastic fluid are presented, for a number of substrate geometries. In the third problem studied, we use lubrication analysis to characterize the steady translational velocity and rotation rate of a sphere that moves down an inclined plane, which is coated with a thin film of fluid, under the action of gravity. In this case, an understanding of the fluid flow must be used to determine the forces that act on the sphere, which set the translational speed and rotation rate of the sphere. Finally, we study the effect of the Coriolis force on the Rayleigh capillary instability of a cylindrical column of viscous liquid surrounded by a second viscous liquid in a rotating frame. Our study comprises numerical results of a temporal linear stability
PIV measurements of flow around an arbitrarily moving free surface
NASA Astrophysics Data System (ADS)
Park, Jinsoo; Im, Sunghyuk; Sung, Hyung Jin; Park, Jun Sang
2015-03-01
We present an image preprocessing method for particle image velocimetry (PIV) measurements of flow around an arbitrarily moving free surface. When performing PIV measurements of free surface flows, the interrogation windows neighboring the free surface are vulnerable to a lack, or even an absence, of seeding particles, which induces less reliable measurements of the velocity field. In addition, direct measurements of the free surface velocity using PIV have been challenging due to the intermittent appearance of the arbitrarily moving free surface. To address the aforementioned limitations, the PIV images with a curvilinear free surface can be treated to be suitable for a structured interrogation window arrangement in a Cartesian grid. The proposed image preprocessing method is comprised of a free surface detection method and an image transform process. The free surface position was identified using a free surface detection method based on multiple textons. The detected free surface points were used to transform PIV images of a curvilinear free surface into images with a straightened free surface using a cubic Hermite spline interpolation scheme. After the image preprocessing, PIV algorithms can be applied to the treated PIV images. The fluid-only region velocities were measured using standard PIV method with window deformation, and the free surface velocities were resolved using PIV/interface gradiometry method. The velocity field in the original PIV images was constructed by inverse transforming that in the transformed images. The accuracy of the proposed method was quantitatively evaluated with two sets of synthetic PIV images, and its applicability was examined by applying the present method to free surface flow images, specifically sloshing flow images.
RIPPLE - A new model for incompressible flows with free surfaces
NASA Technical Reports Server (NTRS)
Kothe, D. B.; Mjolsness, R. C.
1991-01-01
A new free surface flow model, RIPPLE, is summarized. RIPPLE obtains finite difference solutions for incompressible flow problems having strong surface tension forces at free surfaces of arbitrarily complex topology. The key innovation is the continuum surface force model which represents surface tension as a (strongly) localized volume force. Other features include a higher-order momentum advection model, a volume-of-fluid free surface treatment, and an efficient two-step projection solution method. RIPPLE's unique capabilities are illustrated with two example problems: low-gravity jet-induced tank flow, and the collision and coalescence of two cylindrical rods.
Free-Surface Induced Axial Flows in Oscillating Cylinder Wakes
NASA Astrophysics Data System (ADS)
Voorhees, A. V.; Wei, T.
1999-11-01
The problem to be discussed is that of the effect of a free surface on flow along the cores of Karman vortices shed from an inverted pendulum type cylinder. The motivation for this research was to determine the nature and extent of free surface effects in flow induced vibrations of structures aligned perpendicular to and penetrating a free surface. The 2.54 cm diameter cylinder studied had a low mass ratio and was attached at the bottom end by a leaf spring; the free, upper end protruded through the free surface of a large water tunnel. It has been observed that the vortex induced oscillations of the cylinder resulted in very strong axial flows upward toward the free surface, even for small amplitude motions. Approaching the free surface, however, there was an equal likelyhood that the axial flows would be directed down away from the surface. The connection between the direction of axial flow and the cylinder motion will be described using single and two view flow visualization techniques.
Application of the PTT model to axisymmetric free surface flows
NASA Astrophysics Data System (ADS)
Merejolli, R.; Paulo, G. S.; Tomé, M. F.
2013-10-01
This work is concerned with numerical simulation of axisymmetric viscoelastic free surface flows using the Phan-Thien-Tanner (PTT) constitutive equation. A finite difference technique for solving the governing equations for unsteady incompressible flows written in Cylindrical coordinates on a staggered grid is described. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are applied. The numerical method is verified by comparing numerical predictions of fully developed flow in a pipe with the corresponding analytic solutions. To demonstrate that the numerical method can simulate axisymmetric free surface flows governed by the PTT model, numerical results of the flow evolution of a drop impacting on a rigid dry plate are presented. In these simulations, the rheological effects of the parameters ɛ and ξ are investigated.
Energy-decomposition analysis for viscous free-surface flows.
Colagrossi, Andrea; Bouscasse, Benjamin; Marrone, Salvatore
2015-11-01
This work is dedicated to the energy decomposition analysis of a viscous free-surface flow. In the presence of a free surface, the viscous dissipation for a Newtonian liquid can be decomposed into two terms: an enstrophy component and a free-surface deformation component. The latter requires the evaluation of volume and surface integrals in the meshless framework. The analysis is based on the weakly compressible smoothed particle hydrodynamics formalism. The behavior of the energy terms is studied in standing wave problems by changing the viscosity and the wave amplitude. Finally, an analysis of a complex shallow water breaking wave case is provided. It is shown that in presence of intense breaking phenomena the two energy components are always comparable, whereas generally the free surface component is dominant on the viscous dissipation of gravity waves. PMID:26651775
Local mesh refinement for incompressible fluid flow with free surfaces
Terasaka, H.; Kajiwara, H.; Ogura, K.
1995-09-01
A new local mesh refinement (LMR) technique has been developed and applied to incompressible fluid flows with free surface boundaries. The LMR method embeds patches of fine grid in arbitrary regions of interest. Hence, more accurate solutions can be obtained with a lower number of computational cells. This method is very suitable for the simulation of free surface movements because free surface flow problems generally require a finer computational grid to obtain adequate results. By using this technique, one can place finer grids only near the surfaces, and therefore greatly reduce the total number of cells and computational costs. This paper introduces LMR3D, a three-dimensional incompressible flow analysis code. Numerical examples calculated with the code demonstrate well the advantages of the LMR method.
Numerical study of free surface flow around large obstacles
NASA Astrophysics Data System (ADS)
Zhang, Yanming
In this thesis a numerical model was developed to study three-dimensional turbulent flows around large obstacles in an open channel. With this numerical model, a series of numerical tests was carried out, and the properties of turbulent flows around a single obstacle or a cluster of obstacles were investigated. The origin of this study was to study the flow properties around fish habitat structures. Actually, the numerical model can be applied to the study of general turbulent flows under free surfaces. In the numerical model the three-dimensional Reynolds-averaged Navier-Stokes equations in conjunction with k-epsilon turbulence model were solved in a free surface fitted coordinate system. First, different forms of governing equations for turbulent flow were investigated, and a concise form of fully transformed governing equations in a general curvilinear coordinate system was derived. In the numerical solution the FAVOR (Fractional Area/Volume Obstacle Representation) technique was extended into the free surface fitted coordinate system. With this feature the problem of complex turbulent flow with a free surface and general shaped obstacles could be solved efficiently. To locate the free surface, a method based on integrating the momentum equation in the vertical direction was developed. After study and tests of several popular difference schemes, a QUICK scheme with UMIST limiter was adopted in this numerical model. Several test cases were presented to demonstrate the present numerical model. The first test case was to simulate a submerged hydraulic jump. The calculated velocity, free surface profile and turbulence properties of the flow showed a close match with the experimental data. The second test was a submerged hydraulic jump with a baffle sill. The comparison between numerical and experimental data indicated that the current numerical model could catch the general flow structures of the submerged hydraulic jumps. The last two test cases were flows around a
The surface layer for free-surface turbulent flows
NASA Astrophysics Data System (ADS)
Shen, Lian; Zhang, Xiang; Yue, Dick K. P.; Triantafyllou, George S.
1999-05-01
Direct numerical simulation (DNS) is used to examine low Froude number free-surface turbulence (FST) over a two-dimensional mean shear flow. The Navier Stokes equations are solved using a finite-difference scheme with a grid resolution of 1283. Twenty separate simulations are conducted to calculate the statistics of the flow. Based on the velocity deficit and the vertical extent of the shear of the mean flow, the Reynolds number is 1000 and the Froude number is 0.7. We identify conceptually and numerically the surface layer, which is a thin region adjacent to the free surface characterized by fast variations of the horizontal vorticity components. This surface layer is caused by the dynamic zero-stress boundary conditions at the free surface and lies inside a thicker blockage (or ‘source’) layer, which is due to the kinematic boundary condition at the free surface. The importance of the outer blockage layer is manifested mainly in the redistribution of the turbulence intensity, i.e. in the increase of the horizontal velocity fluctuations at the expense of the vertical velocity fluctuation. A prominent feature of FST is vortex connections to the free surface which occur inside the surface layer. It is found that as hairpin-shaped vortex structures approach the free surface, their ‘head’ part is dissipated quickly in the surface layer, while the two ‘legs’ connect almost perpendicularly to the free surface. Analysis of the evolution of surface-normal vorticity based on vortex surface-inclination angle shows that both dissipation and stretching decrease dramatically after connection. As a result, vortex structures connected to the free surface are persistent and decay slowly relative to non-connected vorticities. The effects of surface and blockage layers on the turbulence statistics of length scales, Reynolds-stress balance, and enstrophy dynamics are examined, which elucidate clearly the different turbulence mechanisms operating in the respective near
Free surface and flow problem for a viscous liquid
Zaytsev, M. L. Akkerman, V. B.
2011-10-15
An exact closed system of equations is proposed for describing the shape of the free surface of a viscous steady-state liquid in the 2D case in terms of the surface itself. A method that lowers the dimensionality in the Navier-Stokes equation is demonstrated, and its application in problems of steady-state flow past solids is considered.
Wetting and free surface flow modeling for potting and encapsulation.
Brooks, Carlton, F.; Brooks, Michael J. (Los Alamos National Laboratory, Los Alamos, NM); Graham, Alan Lyman; Noble, David F. ); Notz, Patrick K.; Hopkins, Matthew Morgan; Castaneda, Jaime N.; Mahoney, Leo James; Baer, Thomas A.; Berchtold, Kathryn (Los Alamos National Laboratory, Los Alamos, NM); Adolf, Douglas Brian; Wilkes, Edward Dean; Rao, Rekha Ranjana; Givler, Richard C.; Sun, Amy Cha-Tien; Cote, Raymond O.; Mondy, Lisa Ann; Grillet, Anne Mary; Kraynik, Andrew Michael
2007-06-01
As part of an effort to reduce costs and improve quality control in encapsulation and potting processes the Technology Initiative Project ''Defect Free Manufacturing and Assembly'' has completed a computational modeling study of flows representative of those seen in these processes. Flow solutions are obtained using a coupled, finite-element-based, numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. In addition, two commercially available codes, ProCAST and MOLDFLOW, are also used on geometries representing encapsulation processes at the Kansas City Plant. Visual observations of the flow in several geometries are recorded in the laboratory and compared to the models. Wetting properties for the materials in these experiments are measured using a unique flowthrough goniometer.
Unbounded wall flow with free surface waves and horizontal shear
NASA Astrophysics Data System (ADS)
Lapham, Gary; McHugh, John
2015-11-01
Free surface waves in the presence of a non-uniform shear flow are treated. The shear flow of interest varies with both the transverse and vertical coordinates, U (y , z) . Initial results treat a mean flow varying only with the transverse, U (y) . The domain is bounded on one side by a flat rigid vertical wall and is unbounded on the other side. The mean flows considered here are nonzero near the vertical wall and approach zero far from the wall, e.g. U =e-γy . The flowfield is treated as inviscid but rotational. Linear solutions are obtained using a nonuniform coordinate transformation that converts the free surface boundary condition into a modified Bessel equation. Velocity components are expanded in modified Bessel functions of the first kind of purely imaginary order. The dispersion relation for steady waves are found with wavespeeds outside the range of U, matching previous results for a flow bounded on both sides. Corresponding eigenvectors show a sequence of wave profiles of increasing complexity near the wall. The wave amplitude approaches zero far from the wall.
Flow past a cylinder near a free surface
NASA Astrophysics Data System (ADS)
Delaney, Keegan; Vanella, Marcos; Balaras, Elias; Riaz, Amir
2013-11-01
Flow past a cylinder close to a free surface gives rise flow phenomena that are very different from ones in cases where the cylinder is fully submerged. In this study we will report resolved large-eddy simulations, where we examine the effects of various parameters that have been shown to effect the flow phenomena in previous experimental studies. In all computations a Navier-Stokes solver for multiphase incompressible flows with immersed boundaries and Adaptive Mesh Refinement is utilized. It employs level-set techniques to sharply define the interface between different phases. A fractional step method is used to solve the momentum and continuity equations, which results in a variable coefficient Poisson pressure equation. Proper jump conditions are applied to the Poisson pressure equation to accurately capture the jump in pressure that results from surface tension between different phases. Scalability and efficiency were placed at a premium during development of the solver, which has been tested to core counts on the order of 10,000. We will present details on the interactions between the free surface and vortices shed from the cylinder and their impact in the structure of the immediate wake and air entrainment. Supported by ONR N000141110588 monitored by Dr. Ki-Han Kim.
Free-surface flow over a semi-circular obstruction
NASA Astrophysics Data System (ADS)
Lowery, Kristen; Liapis, Stergios
1999-05-01
The fully non-linear free-surface flow over a semi-circular bottom obstruction was studied numerically in two dimensions using a mixed Eulerian-Lagrangian formulation. The problem was solved in the time domain that allows the prediction of a number of transient phenomena, such as the generation of upstream advancing solitary waves, as well as the simulation of wave breaking. A parametric study was performed for a range of values of the depth-based Froude number up to 2.5 and non-dimensional obstacle heights, up to 0.9. When wave breaking does not occur, three distinct flow regimes were identified: subcritical, transcritical and supercritical. When breaking occurs it may be of any type: spilling, plunging or surging. In addition, for values of the Froude number close to 1, the upstream solitary waves break. A systematic study was undertaken to define the boundaries of each type of breaking and non-breaking pattern and to determine the drag and lift coefficients, free-surface profile characteristics and transient behavior. Copyright
JOSEPHINE: A parallel SPH code for free-surface flows
NASA Astrophysics Data System (ADS)
Cherfils, J. M.; Pinon, G.; Rivoalen, E.
2012-07-01
JOSEPHINE is a parallel Smoothed Particle Hydrodynamics program, designed to solve unsteady free-surface flows. The adopted numerical scheme is efficient and has been validated on a first case, where a liquid drop is stretched over the time. Boundary conditions can also be modelled, as it is demonstrated in a second case: the collapse of a water column. Results show good agreement with both reference numerical solutions and experiments. The use of parallelism allows significant reduction of the computational time, even more with large number of particles. JOSEPHINE has been written so that any untrained developers can handle it easily and implement new features. Catalogue identifier: AELV_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AELV_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 5139 No. of bytes in distributed program, including test data, etc.: 22 833 Distribution format: tar.gz Programming language: Fortran 90 and OpenMPI Computer: All shared or distributed memory parallel processors, tested on a Xeon W3520, 2.67 GHz. Operating system: Any system with a Fortran 90 compiler and MPI, tested on Debian Linux. Has the code been vectorised or parallelised?: The code has been parallelised but has not been explicitly vectorised. RAM: Dependent upon the number of particles. Classification: 4.12 Nature of problem:JOSEPHINE is designed to solve unsteady incompressible flows with a free-surface and large deformations. Solution method:JOSEPHINE is an implementation of Smoothed Particle Hydrodynamics. SPH is a Lagrangian mesh free particle method, thus, no explicit tracking procedure is required to catch the free surface. Incompressibility is satisfied using a weakly compressible model. Boundary conditions at walls are enforced by means of the ghost particles
Finite element simulation of temperature dependent free surface flows
NASA Technical Reports Server (NTRS)
Engelman, M. S.; Sani, R. L.
1985-01-01
The method of Engelman and Sani (1984) for a finite-element simulation of incompressible surface flows with a free and/or moving fluid interface, such as encountered in crystal growth and coating and polymer technology, is extended to temperature-dependent flows, including the effect of temperature-dependent surface tension. The basic algorithm of Saito and Scriven (1981) and Ruschak (1980) has been generalized and implemented in a robust and versatile finite-element code that can be employed with relative ease for the simulation of free-surface problems in complex geometries. As a result, the costly dependence on the Newton-Raphson algorithm has been eliminated by replacing it with a quasi-Newton iterative method, which nearly retains the superior convergence properties of the Newton-Raphson method.
A kinematic conservation law in free surface flow
NASA Astrophysics Data System (ADS)
Gavrilyuk, Sergey; Kalisch, Henrik; Khorsand, Zahra
2015-06-01
The Green-Naghdi system is used to model highly nonlinear weakly dispersive waves propagating at the surface of a shallow layer of a perfect fluid. The system has three associated conservation laws which describe the conservation of mass, momentum, and energy due to the surface wave motion. In addition, the system features a fourth conservation law which is the main focus of this note. It will be shown how this fourth conservation law can be interpreted in terms of a concrete kinematic quantity connected to the evolution of the tangent velocity at the free surface. The equation for the tangent velocity is first derived for the full Euler equations in both two and three dimensional flows, and in both cases, it gives rise to an approximate balance law in the Green-Naghdi theory which turns out to be identical to the fourth conservation law for this system. It is also shown that the conservation equation for the tangent velocity at the free surface appears as an endpoint case of a more general conservation equation for tangent velocities along material surfaces in the body of the fluid.
Aging and free surface flow of a thixotropic fluid
NASA Astrophysics Data System (ADS)
Huynh, H. T.; Roussel, N.; Coussot, P.
2005-03-01
Free surface flows of thixotropic fluids such as paints, self-compacting concrete, or natural mudflows are of noticeable practical interest. Here we study the basic characteristics of the uniform flow of a layer of thixotropic fluid under gravity. A theoretical approach relying on a simple thixotropy constitutive equation shows that after some time at rest over a small slope angle the fluid layer should start to flow rather abruptly beyond a new, larger, critical slope angle. The theory also predicts that the critical time at which the layer velocity should significantly increase is proportional to the duration of the preliminary rest and tends to infinity when the new slope approaches the critical slope. Experiments carried out with different suspensions show that the qualitative trends of the flows are in very good agreement with the theoretical predictions, except that the critical time for flow start appears to be proportional to a power 0.6 of the time of rest whereas the theory predicts a linear dependence. We show that this indicates a restructuration process at rest differing from the restructuration process under flow.
Energy flow and energy dissipation in a free surface.
NASA Astrophysics Data System (ADS)
Goldburg, Walter; Cressman, John
2005-11-01
Turbulent flows on a free surface are strongly compressible [1] and do not conserve energy in the absence of viscosity as bulk fluids do. Despite violation of assumptions essential to Kolmogorov's theory of 1941 (K41) [2, 3], surface flows show strong agreement with Kolmogorov scaling, though intermittency is larger there. Steady state turbulence is generated in a tank of water, and the spatially averaged energy flux is measured from the four-fifth's law at each instant of time. Likewise, the energy dissipation rate as measured from velocity gradients is also a random variable in this experiment. The energy flux - dissipation rate cross-correlation is measured to be correlated in incompressible bulk flows, but strongly anti-correlated on the surface. We argue that the reason for this discrepancy between surface and bulk flows is due to compressible effects present on the surface. [1] J. R. Cressman, J. Davoudi, W. I. Goldburg, and J. Schumacher, New Journal of Physics, 6, 53, 2004. [2] U. Frisch. Turbulence: The legacy of A. N. Kolmogorov, Cambridge University Press, Cambridge, 1995. [3] A. N. Kolmogorov, Doklady Akad. Nauk SSSR, 32, 16, 1941.
Rotating polygon instability of a swirling free surface flow.
Tophøj, L; Mougel, J; Bohr, T; Fabre, D
2013-05-10
We explain the rotating polygon instability on a swirling fluid surface [G. H. Vatistas, J. Fluid Mech. 217, 241 (1990) and Jansson et al., Phys. Rev. Lett. 96, 174502 (2006)] in terms of resonant interactions between gravity waves on the outer part of the surface and centrifugal waves on the inner part. Our model is based on potential flow theory, linearized around a potential vortex flow with a free surface for which we show that unstable resonant states appear. Limiting our attention to the lowest order mode of each type of wave and their interaction, we obtain an analytically soluble model, which, together with estimates of the circulation based on angular momentum balance, reproduces the main features of the experimental phase diagram. The generality of our arguments implies that the instability should not be limited to flows with a rotating bottom (implying singular behavior near the corners), and indeed we show that we can obtain the polygons transiently by violently stirring liquid nitrogen in a hot container. PMID:23705710
Velocity measurements on highly turbulent free surface flow using ADV
NASA Astrophysics Data System (ADS)
Cea, L.; Puertas, J.; Pena, L.
2007-03-01
The 3D instantaneous velocity recorded with an acoustic Doppler velocimeter (ADV) in a highly turbulent free surface flow is analysed using several filters in order to eliminate the corrupted data from the sample. The filters used include the minimum/maximum threshold, the acceleration threshold, and the phase-space threshold. Following some ideas of the phase-space filter, a new method based on the 3D velocity cross-correlation is proposed and tested. A way of computing the constants of the acceleration threshold method is proposed, so no parameters need to be fixed by the user, which makes the filtering process simpler, more objective and more efficient. All the samples analysed are highly turbulent. Nevertheless, the turbulence intensity and the air entrainment vary widely in the flow under study, which produces data records of different quality depending on the measurement point. The performance of the filtering methods when applied to samples of different quality, and the effects of the filtering process in the mean velocity, turbulent kinetic energy and frequency spectra are discussed.
On the flow structure of cloud cavitating flow around an axisymmetric body near the free surface
NASA Astrophysics Data System (ADS)
Wang, Yiwei; Wu, Xiaocui; Huang, Chenguang; Yu, XianXian
2015-12-01
The influence of the free surface on the cavitating flow is an important issue involved in the design of high speed surface vehicles. In the present paper, unsteady cavitating turbulent flow around an axisymmetric body near the free surface was investigated by both launching experiment and LES simulation. The vortex motion induced by cavity shedding under the effect of the free surface is emphatically analyzed by comparing with the submerged condition. The vortex shedding process around the projectile is not synchronized, while the asymmetric characteristic in collapse process is more remarkable, with the generation of multiple vortex ring structures.
Measurement of the interaction between the flow and the free surface of a liquid
Okamoto, Koji; Schmidl, W.D.; Philip, O.G.
1995-09-01
The interaction between the flow and free surface was evaluated measuring the velocity distribution and surface movement simultaneously. The test section was a rectangular tank having a free surface. A rectangular nozzle was set near the free surface, causing the wavy free surface condition. The flow under the free surface was visualized by a laser light sheet and small tracer particles. With image processing techniques, the movement of the free surface and the movement of the particles were simultaneously measured from the recorded images, resulting in the velocity distributions and surface locations. Then, the interactions between the flow and free surface were evaluated using the form of turbulent energy and surface-related turbulent values. By increasing the turbulent energy near the free surface, the fluctuations of the free surface height and the inclination of the free surface were increased. The higher fluctuation of horizontal velocity was related to the higher surface position and negative inclination. The image processing technique is found to be very useful to evaluate the interaction between free surface and flow.
Free surface flow simulation with application to bluff body flow control
NASA Astrophysics Data System (ADS)
Kocabiyik, S.; Bozkaya, C.
2015-03-01
To better understand the interaction of a free surface wave motion with moving bluff bodies, a two-dimensional numerical study of the forced streamwise oscillation of a circular cylinder beneath a free surface is conducted based on a two-fluid model. Computations are carried out at a Reynolds number of R = 200, a fixed displacement amplitude, A = 0.13 and the forcing frequency-to-natural shedding frequency ratios, f/ f 0 = 1.5,2.5,3.5. Finite volume discretization of the special integral form of two-dimensional continuity and unsteady Navier-Stokes equations (when a solid body is present) are performed on a fixed Cartesian grid. Improved volume-of-fluid method is used to discretize the free surface. The laminar asymmetric flow regimes in the near wake region and the fluid forces are analyzed at a fixed Froude number of Fr = 0.4 and for submergence depths at h = 0.25,0.5,0.75. A comparison of the present results with the case in the absence of a free surface is also included to illustrate the effects of inclusion of a free surface. The code validation in special cases shows good comparisons with previous numerical and experimental results. Flow regime analyses include free surface physics-based analysis, and results confirm findings of a recent work of Brøns et al. [25].
Method of driving liquid flow at or near the free surface using magnetic microparticles
Snezhko, Oleksiy; Aronson, Igor; Kwok, Wai-Kwong; Belkin, Maxim V.
2011-10-11
The present invention provides a method of driving liquid flow at or near a free surface using self-assembled structures composed of magnetic particles subjected to an external AC magnetic field. A plurality of magnetic particles are supported at or near a free surface of liquid by surface tension or buoyancy force. An AC magnetic field traverses the free surface and dipole-dipole interaction between particles produces in self-assembled snake structures which oscillate at the frequency of the traverse AC magnetic field. The snake structures independently move across the free surface and may merge with other snake structures or break up and coalesce into additional snake structures experiencing independent movement across the liquid surface. During this process, the snake structures produce asymmetric flow vortices across substantially the entirety of the free surface, effectuating liquid flow across the free surface.
Selected flows with free surfaces: Streams and drops
NASA Astrophysics Data System (ADS)
Kowalewski, Tomasz A.
1995-03-01
The basic purpose of the research described in this article was to develop a non contact method for diagnosing the physical parameters of the free surface of a liquid using drop oscillation analysis. In particular, the purpose is to measure the temperature of an evaporating surface. The realization of this goal has led to the development of new experimental techniques which make it possible to record fast processes using video and digital imaging equipment. Experimental studies of the process of the formation of drops as a result of the controlled breakup of a stream revealed the existence of an additional phase in the process based on the formation of microstreams and microsatellites with micrometer-like dimensions. A comparison of measurement results with Eggers' asymptotic model (23) confirmed the model's basic assumption of the local nature of the final phase in the disintegration of the stream, which at the same time points to the existence of a number of discrepancies which provide evidence of the limitations of this approximation. The next part of the article presents the results of observations of the instability of streams of liquid caused by its evaporation. In an attempt to analyze the mechanisms which initiate the turbulence of the evaporating surface, the author focused on surface tension gradients as an essential factor in the destabilization of small-diameter streams. The author also described the occurrence of a number of new phenomena in the destabilization of a stream, including the separation of surface fragments, their stabilization by the flow of vapor, and a quasistable change in the trajectory of the stream. The author also developed an experimental method which makes it possible to detect and produce a precise description of the deformation of drops. Measurements of the oscillations of small drops in the air led to the development of a complete non-linear model of the oscillations of a viscous drop and made it possible to verify simplified
Computation of effective free surfaces in two phase flows
NASA Astrophysics Data System (ADS)
Yapalparvi, R.; Protas, B.
2012-08-01
In this investigation we revisit the concept of "effective free surfaces" arising in the solution of the time-averaged fluid dynamics equations in the presence of free boundaries. This work is motivated by applications of the optimization and optimal control theory to problems involving free surfaces, where the time-dependent formulations lead to many technical difficulties which are, however, alleviated when steady governing equations are used instead. By introducing a number of precisely stated assumptions we develop and validate an approach in which the interface between the different phases, understood in the time-averaged sense, is sharp. In the proposed formulation, the terms representing the fluctuations of the free boundaries and of the hydrodynamic quantities appear as boundary conditions on the effective surface and require suitable closure models. As a simple model problem we consider impingement of free-falling droplets onto a fluid in a pool with a free surface, and a simple algebraic closure model is proposed for this system. The resulting averaged equations are of the free-boundary type and an efficient computational approach based on shape optimization formulation is developed for their solution. The computed effective surfaces exhibit consistent dependence on the problem parameters and compare favorably with the results obtained when the data from the actual time-dependent problem are used in lieu of the closure model.
An experimental study of free surface deformation in oscillatory thermocapillary flow
NASA Technical Reports Server (NTRS)
Lin, Jianlian; Kamotani, Yasuhiro; Ostrach, Simon
1995-01-01
Free surface deformation of a test fluid induced by steady and oscillatory thermocapillary flow in a small cylindrical container (diameter = 4.8 mm) is studied experimentally. The fluid (2 Cs oil) is heated by a cylindrical wire (diameter = 0.48 mm) placed along the centerline of the cylinder. The relative displacement and oscillation frequency of free surface are measured in a radial cross-section for various applied temperature differences. It is found that the amplitude of the free surface height oscillations is large near the heater and the cold wall and relatively small in the middle part. The frequency of the free surface oscillations is same as the frequency of the temperature oscillations. Based on these data, the free surface motion during oscillations is delineated.
Improving the Dupuit-Forchheimer Approximation for Free Surface Flow in an Unconfined Aquifer
NASA Astrophysics Data System (ADS)
Knight, J. H.
2003-12-01
The classical Dupuit-Forchheimer (DF) approximation for groundwater free surface flow in an unconfined aquifer assumes that the vertical component of the seepage velocity is zero. This assumption is expected to be least accurate when there is non-zero accretion at the free surface. The DF approximation leads to a nonlinear diffusion equation satisfied by the height of the free surface. The general principles of integral methods used by Yves Parlange are to assume some simple approximate shape for some unknown function, and then to choose the parameters of this function to satisfy some known integral relation of the flow system. The DF approximation is improved by assuming that the vertical velocity component is zero at the impermeable horizontal base, and increases linearly to its unknown value at the free surface. The well known Guirinsky potential which depends only on the free surface height corresponds to the DF assumptions. Youngs used an integral relation to define a new potential which depends on the free surface height and also on the vertical velocity component, and which for steady flow satisfies a Poisson equation in the horizontal coordinates. We use the assumption of linear variation of vertical velocity to calculate an approximation to the Youngs potential. In some simple flow systems such as the classical dam problem this leads to a simple differential equation for the free surface height, which can be solved numerically. ln some cases simple explicit approximations can be found for quantities of interest, such as the maximum free surface height between drainage ditches.
Free-surface flow around an appended hull
NASA Astrophysics Data System (ADS)
Lungu, A.; Pacuraru, F.
2010-08-01
The prediction of the total drag experienced by an advancing ship is a complicated problem which requires a thorough understanding of the hydrodynamic forces acting on the hull, the physical processes from which these forces arise as well as their mutual interaction. A general numerical method to predict the hydrodynamic performance of a twin-propeller combatant ship is presented in the paper, which describes the solution of a RANS solver coupled with a body force method as an attempt in investigating the flow features around the ship hull equipped with rotating propellers and rudders. A special focus is made on the propeller non-symmetrical inflow field, aimed at obtaining the necessary data for the propulsive performances evaluation as well as for the propeller final design. The reported work allows not only the performance evaluation for the overall performances of a hull, but also leads to the development, implementation and validation of new concepts in modeling the turbulent vortical flows, with direct connection to the ship propulsion problem.
Nonlinear dynamics and breakup of free-surface flows
Eggers, J.
1997-07-01
Surface-tension-driven flows and, in particular, their tendency to decay spontaneously into drops have long fascinated naturalists, the earliest systematic experiments dating back to the beginning of the 19th century. Linear stability theory governs the onset of breakup and was developed by Rayleigh, Plateau, and Maxwell. However, only recently has attention turned to the nonlinear behavior in the vicinity of the singular point where a drop separates. The increased attention is due to a number of recent and increasingly refined experiments, as well as to a host of technological applications, ranging from printing to mixing and fiber spinning. The description of drop separation becomes possible because jet motion turns out to be effectively governed by one-dimensional equations, which still contain most of the richness of the original dynamics. In addition, an attraction for physicists lies in the fact that the separation singularity is governed by universal scaling laws, which constitute an asymptotic solution of the Navier-Stokes equation before and after breakup. The Navier-Stokes equation is thus continued uniquely through the singularity. At high viscosities, a series of noise-driven instabilities has been observed, which are a nested superposition of singularities of the same universal form. At low viscosities, there is rich scaling behavior in addition to aesthetically pleasing breakup patterns driven by capillary waves. The author reviews the theoretical development of this field alongside recent experimental work, and outlines unsolved problems. {copyright} {ital 1997} {ital The American Physical Society}
Effect of free surface on near-wake flow of elliptic cylinders with different aspect ratios
NASA Astrophysics Data System (ADS)
Lee, Sang Joon; Daichin, -
2003-04-01
The flow fields behind elliptic cylinders with different aspect ratios adjacent to a free surface were investigated experimentally in a circulating water channel. The elliptic cylinders tested in this study have same cross section area. For each elliptic cylinder, the experiments were carried out under different conditions by varying the submergence depth of the cylinder beneath the free surface. The flow fields were measured using a single-frame double-exposure PIV system. For each experimental condition, 350 instantaneous velocity fields were captured and ensemble-averaged to obtain the mean flow field information and spatial distribution of turbulent statistics. The near-wakes can be basically classified into three typical patterns, which are formation of Coanda effect, generation of substantial jet-like flow, and attachment of jet flow to the free surface. The general flow structures behind the elliptic cylinder are similar to previous results for a circular submerged near to a free surface. However, the wake width and the angle of downward deflection of the shear layer developed from the lower surface of the elliptic cylinder are smaller than those for the circular cylinder. These trends are enhance with increasing of the cylinder aspect ratios.
Free-surface potential flow of an ideal fluid due to a singular sink
NASA Astrophysics Data System (ADS)
Mestnikova, A. A.; Starovoitov, V. N.
2016-06-01
A two-dimensional problem of a potential free-surface flow of an ideal incompressible fluid caused by a singular sink is considered. The sink is placed at the horizontal bottom of the channel. By employing a conformal map, the problem is equivalently rewritten in the unit circle. After that, it is investigated by the Levi — Civita technique with the extraction of the singular part of the flow that corresponds to the sink. We derive a Nekrasov type equation that describes exactly the form of the free boundary. This equation is studied at first numerically and then by an exact mathematical technique. It is shown that for the Froude number greater than some particular value, there exists a unique solution of the problem such that the free surface decreases monotonically when moving from the infinity to the sink. At the point over the sink, the free surface has a cusp.
The capturing of free surfaces in incompressible multi-fluid flows
NASA Astrophysics Data System (ADS)
Pan, Dartzi; Chang, Chih-Hao
2000-05-01
By treating it as a contact discontinuity in the density field, a free surface between two immiscible fluids can be automatically captured by the enforcement of conservation laws. A surface-capturing method of this kind requires no special tracking or fitting treatment for the free surface, thereby offering the advantage of algorithm simplicity over the surface-tracking or the surface-fitting method. A surface-capturing method based on a new multi-fluid incompressible Navier-Stokes formulation is developed. It is applied to a variety of free-surface flows, including the Rayleigh-Taylor instability problem, the ship waves around a Wigley hull and a model bubble-rising problem to demonstrate the validity and versatility of the present method. Copyright
Phase-detection measurements in free-surface turbulent shear flows
NASA Astrophysics Data System (ADS)
Chanson, Hubert
2016-04-01
High-velocity self-aerated flows are described as ‘white waters’ because of the entrained air bubbles. The air entrainment induces a drastic change in the multiphase flow structure of the water column and this leads to significant bubble-turbulence interactions, turbulence modulation and associated mixing processes impacting on the bulk flow properties. In these high-velocity free-surface turbulent flows, the phase-detection needle probe is a most reliable instrumentation. The signal processing of a phase-detection probe is re-visited herein. It is shown that the processing may be performed on the raw probe signal as well as the thresholded data. The latter yields the time-averaged void fraction, the bubble count rate, the particle chord time distributions and the particle clustering properties within the particulate flow regions. The raw probe signal analysis gives further the auto-correlation time scale and the power spectrum density function. Finally dimensional considerations are developed with a focus on the physical modelling of free-surface flows in hydraulic structures. It is argued that the notion of scale effects must be defined in terms of some specific set of air-water flow properties within well-defined testing conditions, while a number of free-surface flow characteristics are more prone to scale effects than others, even in large-size physical facilities.
Laboratory Study of MHD Effects on Stability of Free-surface Liquid Metal Flow
NASA Astrophysics Data System (ADS)
Burin, M. J.; Ji, H.; McMurtry, K.; Peterson, L.; Giannakis, D.; Rosner, R.; Fischer, P.
2006-10-01
The dynamics of free-surface MHD shear flows is potentially important to both astrophysics (e.g. in the mixing of dense plasma accreted upon neutron star surfaces) and fusion reactors (e.g. in liquid metal ‘first walls’). To date however few relevant experiments exist. In order to study the fundamental physics of such flows, a small-scale laboratory experiment is being built using a liquid gallium alloy flowing in an open- channel geometry. The flow dimensions are nominally 10cm wide, 1cm deep, and 70cm long under an imposed magnetic field up to 7kG, leading to maximum Hartman number of 2000 and maximum Reynolds number of 4x10^5. Two basic physics issues will ultimately be addressed: (1) How do MHD effects modify the stability of the free surface? For example, is the flow more stable (through the suppression of cross-field motions), or less stable (through the introduction of new boundary layers)? We also investigate whether internal shear layers and imposed electric currents can control the surface stability. (2) How do MHD effects modify free-surface convection driven by a vertical and/or horizontal temperature gradient? We discuss aspects of both of these issues, along with detailed descriptions of the experimental device. Pertinent theoretical stability analyses and initial hydrodynamic results are presented in companion posters. This work is supported by DoE under contract #DE-AC02-76-CH03073.
A particle filter to reconstruct a free-surface flow from a depth camera
NASA Astrophysics Data System (ADS)
Combés, Benoit; Heitz, Dominique; Guibert, Anthony; Mémin, Etienne
2015-10-01
We investigate the combined use of a kinect depth sensor and of a stochastic data assimilation (DA) method to recover free-surface flows. More specifically, we use a weighted ensemble Kalman filter method to reconstruct the complete state of free-surface flows from a sequence of depth images only. This particle filter accounts for model and observations errors. This DA scheme is enhanced with the use of two observations instead of one classically. We evaluate the developed approach on two numerical test cases: a collapse of a water column as a toy-example and a flow in an suddenly expanding flume as a more realistic flow. The robustness of the method to depth data errors and also to initial and inflow conditions is considered. We illustrate the interest of using two observations instead of one observation into the correction step, especially for unknown inflow boundary conditions. Then, the performance of the Kinect sensor in capturing the temporal sequences of depth observations is investigated. Finally, the efficiency of the algorithm is qualified for a wave in a real rectangular flat bottomed tank. It is shown that for basic initial conditions, the particle filter rapidly and remarkably reconstructs the velocity and height of the free surface flow based on noisy measurements of the elevation alone.
SPH simulation of free surface flow over a sharp-crested weir
NASA Astrophysics Data System (ADS)
Ferrari, Angela
2010-03-01
In this paper the numerical simulation of a free surface flow over a sharp-crested weir is presented. Since in this case the usual shallow water assumptions are not satisfied, we propose to solve the problem using the full weakly compressible Navier-Stokes equations with the Tait equation of state for water. The numerical method used consists of the new meshless Smooth Particle Hydrodynamics (SPH) formulation proposed by Ferrari et al. (2009) [8], that accurately tracks the free surface profile and provides monotone pressure fields. Thus, the unsteady evolution of the complex moving material interface (free surface) can been properly solved. The simulations involving about half a million of fluid particles have been run in parallel on two of the most powerful High Performance Computing (HPC) facilities in Europe. The validation of the results has been carried out analysing the pressure field and comparing the free surface profiles obtained with the SPH scheme with experimental measurements available in literature [18]. A very good quantitative agreement has been obtained.
Coupling of Smoothed Particle Hydrodynamics with Finite Volume method for free-surface flows
NASA Astrophysics Data System (ADS)
Marrone, S.; Di Mascio, A.; Le Touzé, D.
2016-04-01
A new algorithm for the solution of free surface flows with large front deformation and fragmentation is presented. The algorithm is obtained by coupling a classical Finite Volume (FV) approach, that discretizes the Navier-Stokes equations on a block structured Eulerian grid, with an approach based on the Smoothed Particle Hydrodynamics (SPH) method, implemented in a Lagrangian framework. The coupling procedure is formulated in such a way that each solver is applied in the region where its intrinsic characteristics can be exploited in the most efficient and accurate way: the FV solver is used to resolve the bulk flow and the wall regions, whereas the SPH solver is implemented in the free surface region to capture details of the front evolution. The reported results clearly prove that the combined use of the two solvers is convenient from the point of view of both accuracy and computing time.
Algorithm for the simulation of transient viscoelastic flows with free surfaces
Keunings, R.
1986-01-01
We propose a numerical procedure for solving a class of transient viscoelastic flows with free surfaces. It is based on a Galerkin/Finite Element technique on deforming elements combined with a predictor-corrector scheme. The method is applied to the analysis of jet breakup caused by capillary forces. Non-linear effects known to experimentalists are predicted and a detailed comparison with asymptotic results is carried out.
Algorithm for the simulation of transient viscoelastic flows with free surfaces
Keunings, R.
1984-10-01
We propose a numerical procedure for solving a class of transient viscoelastic flows with free surfaces. It is based on a Galerkin/Finite Element technique on deforming elements combined with a predictor-corrector scheme. The method is applied to the analysis of jet breakup caused by capillary forces. Non-linear effects known to experimentalists are predicted and a detailed comparison with asymptotic results is carried out.
An improved parallel SPH approach to solve 3D transient generalized Newtonian free surface flows
NASA Astrophysics Data System (ADS)
Ren, Jinlian; Jiang, Tao; Lu, Weigang; Li, Gang
2016-08-01
In this paper, a corrected parallel smoothed particle hydrodynamics (C-SPH) method is proposed to simulate the 3D generalized Newtonian free surface flows with low Reynolds number, especially the 3D viscous jets buckling problems are investigated. The proposed C-SPH method is achieved by coupling an improved SPH method based on the incompressible condition with the traditional SPH (TSPH), that is, the improved SPH with diffusive term and first-order Kernel gradient correction scheme is used in the interior of the fluid domain, and the TSPH is used near the free surface. Thus the C-SPH method possesses the advantages of two methods. Meanwhile, an effective and convenient boundary treatment is presented to deal with 3D multiple-boundary problem, and the MPI parallelization technique with a dynamic cells neighbor particle searching method is considered to improve the computational efficiency. The validity and the merits of the C-SPH are first verified by solving several benchmarks and compared with other results. Then the viscous jet folding/coiling based on the Cross model is simulated by the C-SPH method and compared with other experimental or numerical results. Specially, the influences of macroscopic parameters on the flow are discussed. All the numerical results agree well with available data, and show that the C-SPH method has higher accuracy and better stability for solving 3D moving free surface flows over other particle methods.
Simulation of gravity driven free-surface flow in fractured geological media
Kordilla, Jannes; Geyer, Tobias; Tartakovsky, Alexandre M.
2011-08-25
Abstract—Simulation of unsaturated free-surface flow in fractured geological media represents a challenge due to the highly heterogeneous flow field induced by extensive faults, joints and fissures. Free-surface flow in unsaturated media leads to highly intermittent flow regimes and flow velocities well above those assumed for the bulk volume. However, common modeling approaches relying on volume-averaged effective equations fail to capture this flow feature. In this work we present micro-scale flow simulations using a three-dimensional multiphase SPH code. Pairwise fluid-fluid and solid-fluid interaction forces are used to simulate a wide range of wetting conditions encountered on rock surfaces. It is shown that static contact angles for sessiles droplets are independent of the model discretization, i.e. the total amount of particles. Thus, computation times can be reduced without sacrificing qualitative or quantitative information. Furthermore we show that our model is in accordance with general scaling laws for droplet flow.
Free surface flow through rock-fill dams analyzed by FEM with level set approach
NASA Astrophysics Data System (ADS)
Sharif, N. H.; Wiberg, N.-E.; Levenstam, M.
A stabilized-finite element formulation is coupled with a level set technique for computations of incompressible non-linear flow with interfaces between two immiscible fluids. An interface capturing formulation (ICF) for non-linear, free surface, seepage flow in rock-fill dams is proposed. The formulation is derived for two- and three-dimensional flow within a fixed mesh domain. The resulting formulation is general and applicable for various steady and transient two-phase flow problems. FE-refinement is processed for the entire fixed mesh domains. A general solver is also reviewed for large and non-symmetric non-positive definite linear system of equations with the GMRES-update technique based on a Newton-iterative method. The computational procedure has been implemented in MATLAB. A comparison is performed between the 2-D computed test problem for coarse and refined meshes together with some proposed analytical solutions for nonlinear seepage flow with free surface in rock-fill dams. An expansion of the 2-D program code to a 3-D one for a rectangular rock-fill dam is also developed and simulated in MATLAB. The performance of the computations in 3-D is very promising and its opening the future for possible industrial applications using the same simple technique. Computations for a simple 3-D seepage flow problem with free surface in rock-fill dam are included in present paper. A general mesh generator and solver for large scale and complex 3-D flow problems in a real embankment dam is also under construction in C++.
Three-dimensional low Reynolds number flows with a free surface
NASA Technical Reports Server (NTRS)
Degani, D.; Gutfinger, C.
1977-01-01
The two-dimensional leveling problem (Degani, Gutfinger, 1976) is extended to three dimensions in the case where the flow Re number is very low and attention is paid to the free surface boundary condition with surface tension effects included. The no-slip boundary condition on the wall is observed. The numerical solution falls back on the Marker and Cell (MAC) method (Harlow and Welch, 1965) with the computation region divided into a finite number of stationary rectangular cells (or boxes in the 3-D case) and fluid flow traverses the cells (or boxes).
Magnetohydrodynamics and heat transfer in a free-surface, flowing liquid metal experiment
NASA Astrophysics Data System (ADS)
Rhoads, John Rayburn
This thesis examines magnetohydrodynamic (MHD) effects and heat transfer in free-surface, liquid metal flow. The interaction of the flow with the applied magnetic field introduces an anisotropy to the fluid motion, resulting in starkly different behavior depending on the orientation and strength of the applied magnetic field. The transition from a hydrodynamic state to an MHD state was experimentally found to be characterized by the interaction parameter, N = sigma L B20/rho v0 , crossing above 0.5. The experiment consisted of a channel 1 m in length, up to 16 cm in width and 1 to 2.5 cm in depth situated within an electromagnet capable of producing a uniform magnetic field up to 2.7 kG. The fluid velocity was measured with an array of 25 potential probes that resolved the spanwise velocity profile, while waves on the surface of the metal were monitored via a diagnostic that tracked the motion of a laser beam reflected off of the free surface. Lastly, the temperature response of the fluid was recorded with an array of 32 thermocouples embedded in the bottom of the channel and a mid-wavelength infrared camera used to image the free surface. Three distinct sets of experiments were performed to investigate the effects of the applied magnetic field on the dynamics within the flow, studying: 1) changes to vortices in the wake of a cylinder, 2) behavior of surface fluctuations due to the turbulent flow, and 3) vertical heat transfer resulting from heat deposition on the free surface. The first set of experiments examined the wake of a cylinder inserted into the flow with its axis parallel to the magnetic field. Measurements with the potential probes indicated that vortices in the wake became laminarized as the field was increased due to the alignment of the injected vorticity with the applied field. The second set of experiments showed that surface fluctuations were suppressed at high magnetic fields, but the precise nature of the damping depended on the orientation of
Nanometer-scale free surface flow of molten polyethylene from a heated atomic force microscope tip
NASA Astrophysics Data System (ADS)
Ewoldt, Randy; Felts, Jonathan; Somnath, Suhas; King, William
2012-11-01
We experimentally investigate nanometer-scale free surface flow of molten polyethylene from a heated atomic force microscope (AFM) cantilever, a nanofabrication process known as thermal dip-pen nanolithography (tDPN). Fluid is deposited from the AFM tip onto non-porous substrates whether the tip is moving or fixed. We find that polymer flow depends on surface capillary forces and not on shear between tip and substrate. The polymer mass flow rate is sensitive to the temperature-dependent polymer viscosity. Additionally, the flow rate increases when a temperature gradient exists between the tip and substrate. We hypothesize that the polymer flow is governed by thermal Marangoni forces and non-equilibrium wetting dynamics caused by a solidification front within the feature.
Motion of a cylinder adjacent to a free-surface: flow patterns and loading
NASA Astrophysics Data System (ADS)
Zhu, Q.; Lin, J.-C.; Unal, M. F.; Rockwell, D.
The flow structure and loading due to combined translatory and sinusoidal motion of a cylinder adjacent to a free-surface are characterized using a cinema technique of high-image-density particle image velocimetry and simultaneous force measurements. The instantaneous patterns of vorticity and streamline topology are interpreted as a function of degree of submergence beneath the free-surface. The relative magnitudes of the peak vorticity and the circulation of vortices formed from the upper and lower surfaces of the cylinder, as well as vortex formation from the free-surface, are remarkably affected by the nominal submergence. The corresponding streamline topology, interpreted in terms of foci, saddle points, and multiple separation and reattachment points also exhibit substantial changes with submergence. All of these features affect the instantaneous loading of the cylinder. Calculation of instantaneous moments of vorticity and the incremental changes in these moments during the cylinder motion allow identification of those vortices that contribute most substantially to the instantaneous lift and drag. Furthermore, the calculated moments are in general accord with the time integrals of the measured lift and drag acting on the cylinder for sufficiently large submergence.
NASA Astrophysics Data System (ADS)
Halfon, S.; Arenshtam, A.; Kijel, D.; Paul, M.; Weissman, L.; Aviv, O.; Berkovits, D.; Dudovitch, O.; Eisen, Y.; Eliyahu, I.; Feinberg, G.; Haquin, G.; Hazenshprung, N.; Kreisel, A.; Mardor, I.; Shimel, G.; Shor, A.; Silverman, I.; Tessler, M.; Yungrais, Z.
2014-05-01
The free-surface Liquid-Lithium Target, recently developed at Soreq Applied Research Accelerator Facility (SARAF), was successfully used with a 1.9 MeV, 1.2 mA (2.3 kW) continuous-wave proton beam. Neutrons (˜2 × 1010 n/s having a peak energy of ˜27 keV) from the 7Li(p,n)7Be reaction were detected with a fission-chamber detector and by gold activation targets positioned in the forward direction. The setup is being used for nuclear astrophysics experiments to study neutron-induced reactions at stellar energies and to demonstrate the feasibility of accelerator-based boron neutron capture therapy.
Halfon, S.; Feinberg, G.; Arenshtam, A.; Kijel, D.; Weissman, L.; Aviv, O.; Berkovits, D.; Dudovitch, O.; Eisen, Y.; Eliyahu, I.; Haquin, G.; Hazenshprung, N.; Kreisel, A.; Mardor, I.; Shimel, G.; Shor, A.; Silverman, I.; Yungrais, Z.; Paul, M. Tessler, M.
2014-05-15
The free-surface Liquid-Lithium Target, recently developed at Soreq Applied Research Accelerator Facility (SARAF), was successfully used with a 1.9 MeV, 1.2 mA (2.3 kW) continuous-wave proton beam. Neutrons (∼2 × 10{sup 10} n/s having a peak energy of ∼27 keV) from the {sup 7}Li(p,n){sup 7}Be reaction were detected with a fission-chamber detector and by gold activation targets positioned in the forward direction. The setup is being used for nuclear astrophysics experiments to study neutron-induced reactions at stellar energies and to demonstrate the feasibility of accelerator-based boron neutron capture therapy.
NASA Astrophysics Data System (ADS)
Xu, Xiaoyang; Deng, Xiao-Long
2016-04-01
In this paper, an improved weakly compressible smoothed particle hydrodynamics (SPH) method is proposed to simulate transient free surface flows of viscous and viscoelastic fluids. The improved SPH algorithm includes the implementation of (i) the mixed symmetric correction of kernel gradient to improve the accuracy and stability of traditional SPH method and (ii) the Rusanov flux in the continuity equation for improving the computation of pressure distributions in the dynamics of liquids. To assess the effectiveness of the improved SPH algorithm, a number of numerical examples including the stretching of an initially circular water drop, dam breaking flow against a vertical wall, the impact of viscous and viscoelastic fluid drop with a rigid wall, and the extrudate swell of viscoelastic fluid have been presented and compared with available numerical and experimental data in literature. The convergent behavior of the improved SPH algorithm has also been studied by using different number of particles. All numerical results demonstrate that the improved SPH algorithm proposed here is capable of modeling free surface flows of viscous and viscoelastic fluids accurately and stably, and even more important, also computing an accurate and little oscillatory pressure field.
Efficient Computation of N-S Equation with Free Surface Flow Around an ACV on ShirazUCFD Grid
NASA Astrophysics Data System (ADS)
Sheikhalishahi, Seyyed Mehdi; Alizadehrad, Davood; Dastghaibyfard, Gholamhossein; Alishahi, Mohammad Mehdi; Nikseresht, Amir Hossein
This paper presents the application of a parallel high accuracy simulation code for Incompressible Navier-Stokes solution with free surface flow around an ACV (Air Cushion Vehicle) on ShirazUCFD Grid environment. The parallel finite volume code is developed for incompressible Navier-Stokes solver on general curvilinear coordinates system for modeling free surface flows. A single set of dimensionless equations is derived to handle both liquid and air phases in viscous incompressible free surface flow in general curvilinear coordinates. The volume of fluid (VOF) method with lagrangian propagation in computational domain for modeling the free surface flow is implemented. The parallelization approach uses a domain decomposition method for the subdivision of the numerical grid, the SPMD program model and MPICH-G2 as the message passing environment is used to obtain a portable application.
Time-dependent liquid metal flows with free convection and free surfaces
McClelland, M.A.
1990-11-01
A finite element analysis is given for time-dependent liquid metal flows with free convection and free surfaces. Consideration is given to a two-dimensional shallow trough with vertical walls maintained at different temperatures. The spatial formulation incorporates mixed Lagrangian approximations to the velocity, pressure, temperature, and interface position. The time integration method is performed using the Trapezoid Rule with step-size control. The Galerkin method is employed to reduce the problem to a set of nonlinear algebraic equations which are solved with the Newton-Raphson method. Calculations are performed for conditions relevant to the electron beam vaporization of refractory metals. The Prandtl number is 0.015, and Grashof numbers are in the transition region between laminar and turbulent flow. The results reveal the effects of flow intensity, surface-tension gradients, and mesh and time-step refinement.
NASA Technical Reports Server (NTRS)
Orlin, W James; Lindner, Norman J; Butterly, Jack G
1947-01-01
The theory of the hydraulic analogy -- that is, the analogy between water flow with a free surface and two-dimensional compressible gas flow -- and the limitations and conditions of the analogy are discussed. A test was run using the hydraulic analogy as applied to the flow about circular cylinders of various diameters at subsonic velocities extending into the supercritical range. The apparatus and techniques used in this application are described and criticized. Reasonably satisfactory agreement of pressure distributions and flow fields existed between water and air flow about corresponding bodies. This agreement indicated the possibility of extending experimental compressibility research by new methods.
Incompressible SPH Model for Simulating Violent Free-Surface Fluid Flows
NASA Astrophysics Data System (ADS)
Staroszczyk, Ryszard
2014-06-01
In this paper the problem of transient gravitational wave propagation in a viscous incompressible fluid is considered, with a focus on flows with fast-moving free surfaces. The governing equations of the problem are solved by the smoothed particle hydrodynamics method (SPH). In order to impose the incompressibility constraint on the fluid motion, the so-called projection method is applied in which the discrete SPH equations are integrated in time by using a fractional-step technique. Numerical performance of the proposed model has been assessed by comparing its results with experimental data and with results obtained by a standard (weakly compressible) version of the SPH approach. For this purpose, a plane dam-break flow problem is simulated, in order to investigate the formation and propagation of a wave generated by a sudden collapse of a water column initially contained in a rectangular tank, as well as the impact of such a wave on a rigid vertical wall. The results of simulations show the evolution of the free surface of water, the variation of velocity and pressure fields in the fluid, and the time history of pressures exerted by an impacting wave on a wall.
Instabilities in free-surface Hartmann flow at low magnetic Prandtl numbers
NASA Astrophysics Data System (ADS)
Giannakis, Dimitrios
2009-06-01
Free-surface Hartmann flow is the parallel flow of a viscous, electrically conducting, capillary fluid on a planar surface, subject to gravity and a flow- normal magnetic field. This type of flow arises in a variety of industrial and astrophysical contexts, including liquid-metal walls in fusion devices, heavy- ion accelerator targets, and surface layers of white dwarfs and neutron stars. Typically, the Reynolds number, Re >10 4 , is high, and the background magnetic field is strong ( Ha >100, where the Hartmann number, Ha , measures the square root of the ratio of electromagnetic to viscous forces). On the other hand, the magnetic Prandtl number, Pm (the ratio of viscous to magnetic diffusivity), of laboratory fluids is small (e.g., Pm <10 -4 for liquid metals), as is the case in a number of astrophysical models. When the background magnetic field is zero, free-surface Hartmann flow exhibits the so-called soft and hard instability modes; the former being a surface wave destabilized by viscous stresses acting on the free surface, whereas the latter is a shear mode destabilized by positive Reynolds stress associated with an internal critical layer. We study in detail the influence of the external magnetic field on these two instabilities, working in the regime Pm <10^-4. We also consider flows in the inductionless limit, Pr [arrow right]0, where magnetic field perturbations diffuse infinitely fast, and the sole MHD effect is a Lorentz force arising from currents induced by the perturbed fluid motion within the background magnetic field. We have developed a spectral Galerkin method to solve the coupled Orr- Sommerfeld and induction equations, which, in conjunction with suitable stress conditions at the free surface and continuity conditions for the magnetic field, govern the linear stability of free-surface Hartmann flow. Our scheme's discrete bases for the velocity and magnetic fields consist of linear combinations of Legendre polynomials, chosen according to the
NASA Technical Reports Server (NTRS)
Orlin, W James; Lindner, Norman J; Bitterly, Jack G
1947-01-01
The theory of hydraulic analogy, that is, the analogy between water flow with a free surface and two-dimensional compressible gas flow and the limitations and conditions of the analogy are discussed. A test run was made using the hydraulic analogy as applied to the flow about circular cylinders at various diameters at subsonic velocities extending to the super critical range. The apparatus and techniques used in this application are described and criticized. Reasonably satisfactory agreement of pressure distributions and flow fields existed between water and airflow about corresponding bodies. This agreement indicated the possibility of extending experimental compressibility research by new methods.
Strongly coupled partitioned approach for fluid structure interaction in free surface flows
NASA Astrophysics Data System (ADS)
Facci, Andrea Luigi; Ubertini, Stefano
2016-06-01
In this paper we describe and validate a methodology for the numerical simulation of the fluid structure interaction in free surface flows. Specifically, this study concentrates on the vertical impact of a rigid body on the water surface, (i.e. on the hull slamming problem). The fluid flow is modeled through the volume of fluid methodology, and the structure dynamics is described by the Newton's second law. An iterative algorithm guarantees the tight coupling between the fluid and solid solvers, allowing the simulations of lightweight (i.e. buoyant) structures. The methodology is validated comparing numerical results to experimental data on the free fall of different rigid wedges. The correspondence between numerical results and independent experimental findings from literature evidences the reliability and the accuracy of the proposed approach.
Investigation of ALEGRA shock hydrocode algorithms using an exact free surface jet flow solution.
Hanks, Bradley Wright.; Robinson, Allen Conrad
2014-01-01
Computational testing of the arbitrary Lagrangian-Eulerian shock physics code, ALEGRA, is presented using an exact solution that is very similar to a shaped charge jet flow. The solution is a steady, isentropic, subsonic free surface flow with significant compression and release and is provided as a steady state initial condition. There should be no shocks and no entropy production throughout the problem. The purpose of this test problem is to present a detailed and challenging computation in order to provide evidence for algorithmic strengths and weaknesses in ALEGRA which should be examined further. The results of this work are intended to be used to guide future algorithmic improvements in the spirit of test-driven development processes.
Axial Flows in Kármán Vortices Due to Amplitude Modulation and a Free Surface
NASA Astrophysics Data System (ADS)
Voorhees, A.; Benaroya, H.; Wei, T.
2000-11-01
Fluid-structure interaction experiments were conducted in a large free-surface water tunnel facility using a low-mass ratio circular cylinder. The 2.54-cm cylinder was attached at the lower end to the tunnel floor by a leaf spring and the upper end protruded through the free surface. In this manner, the cylinder was free to oscillate as an inverted pendulum in response to the Kármán vortex shedding phenomenon. Far from the free-surface, strong axial flows directed toward the free-surface were observed along the cores of Kármán vortices shed from the oscillating cylinder; such flows were not observed for matched Reynolds number flows where the cylinder was held stationary. However, near the free-surface fluid motion along vortex cores traveled both up and down, i.e. to and from the free-surface. It has been found that this phenomenon is directly linked to an observed beating of the cylinder when the oscillation frequency approaches the cylinder's natural frequency. In this study, the Reynolds number based on cylinder diameter and free stream velocity was 3800. The problem to be discussed is that of the interaction between the amplitude modulated cylinder motion, the ensuing vortex-street, and the free-surface. These interactions will be described using DPIV measurements taken over a range of levels beneath the free-surface. A phenomenological description of the axial vortex induction process will be developed using Kármán vortex strength as a function of distance from the free-surface.
Theory for source-responsive and free-surface film modeling of unsaturated flow
Nimmo, J.R.
2010-01-01
A new model explicitly incorporates the possibility of rapid response, across significant distance, to substantial water input. It is useful for unsaturated flow processes that are not inherently diffusive, or that do not progress through a series of equilibrium states. The term source-responsive is used to mean that flow responds sensitively to changing conditions at the source of water input (e.g., rainfall, irrigation, or ponded infiltration). The domain of preferential flow can be conceptualized as laminar flow in free-surface films along the walls of pores. These films may be considered to have uniform thickness, as suggested by field evidence that preferential flow moves at an approximately uniform rate when generated by a continuous and ample water supply. An effective facial area per unit volume quantitatively characterizes the medium with respect to source-responsive flow. A flow-intensity factor dependent on conditions within the medium represents the amount of source-responsive flow at a given time and position. Laminar flow theory provides relations for the velocity and thickness of flowing source-responsive films. Combination with the Darcy-Buckingham law and the continuity equation leads to expressions for both fluxes and dynamic water contents. Where preferential flow is sometimes or always significant, the interactive combination of source-responsive and diffuse flow has the potential to improve prediction of unsaturated-zone fluxes in response to hydraulic inputs and the evolving distribution of soil moisture. Examples for which this approach is efficient and physically plausible include (i) rainstorm-generated rapid fluctuations of a deep water table and (ii) space- and time-dependent soil water content response to infiltration in a macroporous soil. ?? Soil Science Society of America.
Surface Ripples Generated in a Couette Flow with a Free Surface
NASA Astrophysics Data System (ADS)
Masnadi, N.; Washuta, N.; Duncan, J. H.
2014-11-01
Free surface ripples created by subsurface turbulence in the gap between a vertical surface-piercing moving wall and a parallel fixed wall are studied experimentally. The moving wall is created with the aide of a meter-wide stainless steel belt that travels horizontally in a loop around two rollers with vertically oriented axes, which are separated by 7.5 meters. One of the two 7.5-m-long belt sections between the rollers is in contact with the water in a large open-surface water tank and forms the moving wall. The fixed wall is an acrylic plate located 4 cm from the belt surface. The water surface ripples are measured in a plane normal to the belt using a cinematic LIF technique. Measurements are done at a location about 100 gap widths downstream of the leading edge of the fixed plate in order to have a fully developed flow condition. It is found that the overall RMS surface fluctuations increase linearly with belt speed. The frequency-domain spectra of the surface height fluctuation and its temporal derivative are computed at locations across the gap width and are used to explore the physics of the free surface motions. The support of the Office of Naval Research is gratefully acknowledged.
Oscillating line source in a shear flow with a free surface: critical layer-like contributions
NASA Astrophysics Data System (ADS)
Ellingsen, Simen Å.; Tyvand, Peder A.
2016-07-01
The linearized water-wave radiation problem for an oscillating submerged line source in an inviscid shear flow with a free surface is investigated analytically at finite, constant depth in the presence of a shear flow varying linearly with depth. The surface velocity is taken to be zero relative to the oscillating source, so that Doppler effects are absent. The radiated wave out from the source is calculated based on Euler's equation of motion with the appropriate boundary and radiation conditions, and differs substantially from the solution obtained by assuming potential flow. To wit, an additional wave is found in the downstream direction in addition to the previously known dispersive wave solutions; this wave is non-dispersive and we show how it is the surface manifestation of a critical layer-like flow generated by the combination of shear and mass flux at the source, passively advected with the flow. As seen from a system moving at the fluid velocity at the source's depth, streamlines form closed curves in a manner similar to Kelvin's cat's eye vortices. A resonant frequency exists at which the critical wave resonates with the downstream propagating wave, resulting in a downstream wave pattern diverging linearly in amplitude away from the source.
Choi, Hyunhee; Takezoe, Hideo
2016-01-14
We demonstrate circular flow formation at a surface in homeotropically oriented nematic liquid crystals with a free surface using focused laser beam irradiation. Under a weak laser power, a pit together with an associated circular bulge is formed: the Marangoni effect. Here a diverging molecular flow from the pit (thermocapillary flow) also induces director tilt in the radial direction. Upon increasing the laser power, the pit becomes deeper, and eventually evolves into a circular flow associated with a deeper pit and a subsidiary circular bulge or valley structure. This phenomenon is induced by escaping from excess deformation energy due to a bend deformation of the director. Actually, we confirmed that the circular flow is never formed in the isotropic phase. The handedness of the vortex cannot be controlled by circular polarisation, but is controllable by doping with chiral molecules. This rotational motion (a nematic micro-rotor) is a unique phenomenon only exhibited by anisotropic liquids, and is expected to be applied for novel devices. PMID:26482229
NASA Astrophysics Data System (ADS)
Wieland, M.; Gray, J. M. N. T.; Hutter, K.
1999-08-01
A series of laboratory experiments and numerical simulations have been performed to investigate the rapid fluid-like flow of a finite mass of granular material down a chute with partial lateral confinement. The chute consists of a section inclined at 40° to the horizontal, which is connected to a plane run-out zone by a smooth transition. The flow is confined on the inclined section by a shallow parabolic cross-slope profile. Photogrammetric techniques have been used to determine the position of the evolving boundary during the flow, and the free-surface height of the stationary granular deposit in the run-out zone. The results of three experiments with different granular materials are presented and shown to be in very good agreement with numerical simulations based on the Savage Hutter theory for granular avalanches. The basal topography over which the avalanche flows has a strong channelizing effect on the inclined section of the chute. As the avalanche reaches the run-out zone, where the lateral confinement ceases, the head spreads out to give the avalanche a characteristic ‘tadpole’ shape. Sharp gradients in the avalanche thickness and velocity began to develop at the interface between the nose and tail of the avalanche as it came to rest, indicating that a shock wave develops close to the end of the experiments.
NASA Astrophysics Data System (ADS)
Daichin, K. V.; Lee, Sang Joon
The flow fields behind elliptic cylinders adjacent to a free surface were investigated experimentally in a circulating water channel. A range of cylinder aspect ratios (AR=2, 3, 4) were considered, while the cross-sectional area of the elliptical cylinder was kept constant. The main objective of this study was to investigate the effect of cylinder aspect ratio and a free surface on the flow structure in the near-wake behind elliptic cylinders. For each elliptic cylinder, the flow structure was analyzed for various values of the submergence depth of the cylinder beneath the free surface. The flow fields were measured using a single-frame double-exposure PIV (Particle Image Velocimetry) system. For each experimental condition, 350 instantaneous velocity fields were obtained and ensemble-averaged to obtain the mean velocity field and spatial distribution of the mean vorticity statistics. The results show that near-wake can be classified into three typical flow patterns: formation of a Coanda flow, generation of substantial jet-like flow, and attachment of this jet flow to the free surface. The general flow structure observed behind the elliptic cylinders resembles the structure previously reported for a circular cylinder submerged near a free surface. However, the wake width and the angle of downward deflection of the shear layer developed from the lower surface of the elliptic cylinder differ from those observed for a circular cylinder. These trends are enhanced as cylinder aspect ratio is increased. In addition, the free surface distortion is also discussed in the paper.
NASA Astrophysics Data System (ADS)
Banerjee, Arindam; Kolekar, Nitin
2015-11-01
The current experimental investigation aims at understanding the effect of free surface proximity and associated blockage on near-wake flow-field and performance of a three bladed horizontal axis marine hydrokinetic turbine. Experiments were conducted on a 0.14m radius, three bladed constant chord turbine in a 0.61m ×0.61m test section water channel. The turbine was subjected to various rotational speeds, flow speeds and depths of immersion. Experimental data was acquired through a submerged in-line thrust-torque sensor that was corrected to an unblocked dataset with a blockage correction using measured thrust data. A detailed comparison is presented between blocked and unblocked datasets to identify influence of Reynolds number and free surface proximity on blockage effects. The percent change in Cp was found to be dependent on flow velocity, rotational speed and free surface to blade tip clearance. Further, flow visualization using a stereoscopic particle image velocimetry was carried out in the near-wake region of turbine to understand the mechanism responsible for variation of Cp with rotational speed and free surface proximity. Results revealed presence of slower wake at higher rotational velocities and increased asymmetry in the wake at high free surface proximity.
Free surface flow impact on a vertical wall: a numerical assessment
NASA Astrophysics Data System (ADS)
Pugliese Carratelli, Eugenio; Viccione, Giacomo; Bovolin, Vittorio
2016-03-01
The sudden impact of a free surface flow upon a solid wall is a common occurrence in many situations in nature and technology. The design of marine structures is probably the most obvious example, but also river and dam hydraulics as well as the necessity of understanding flood and debris flow-induced damage have led to theoretical and experimental work on the mechanism of fluid slamming loads. This is therefore a very old and rich research field, which has not yet reached full maturity, so that semi-empirical methods in design practice are still the rule in many sectors. Up-to-date CFD technology with both Eulerian and Lagrangian approaches is employed to investigate highly non-stationary fluid impact on a solid wall. The development of the pressure wave produced by the impact is examined as it propagates and interacts with the fluid boundaries, as well as the subsequent build-up of high-pressure gradients of high fluid velocities. The geometry and the velocity field of the problem considered are very simple, but the results seem to provide new insight, in particular, into the connection between phenomena with different timescales.
Granular-front formation in free-surface flow of concentrated suspensions.
Leonardi, Alessandro; Cabrera, Miguel; Wittel, Falk K; Kaitna, Roland; Mendoza, Miller; Wu, Wei; Herrmann, Hans J
2015-11-01
A granular front emerges whenever the free-surface flow of a concentrated suspension spontaneously alters its internal structure, exhibiting a higher concentration of particles close to its front. This is a common and yet unexplained phenomenon, which is usually believed to be the result of fluid convection in combination with particle size segregation. However, suspensions composed of uniformly sized particles also develop a granular front. Within a large rotating drum, a stationary recirculating avalanche is generated. The flowing material is a mixture of a viscoplastic fluid obtained from a kaolin-water dispersion with spherical ceramic particles denser than the fluid. The goal is to mimic the composition of many common granular-fluid materials, such as fresh concrete or debris flow. In these materials, granular and fluid phases have the natural tendency to separate due to particle settling. However, through the shearing caused by the rotation of the drum, a reorganization of the phases is induced, leading to the formation of a granular front. By tuning the particle concentration and the drum velocity, it is possible to control this phenomenon. The setting is reproduced in a numerical environment, where the fluid is solved by a lattice-Boltzmann method, and the particles are explicitly represented using the discrete element method. The simulations confirm the findings of the experiments, and provide insight into the internal mechanisms. Comparing the time scale of particle settling with the one of particle recirculation, a nondimensional number is defined, and is found to be effective in predicting the formation of a granular front. PMID:26651686
Granular-front formation in free-surface flow of concentrated suspensions
NASA Astrophysics Data System (ADS)
Leonardi, Alessandro; Cabrera, Miguel; Wittel, Falk K.; Kaitna, Roland; Mendoza, Miller; Wu, Wei; Herrmann, Hans J.
2015-11-01
A granular front emerges whenever the free-surface flow of a concentrated suspension spontaneously alters its internal structure, exhibiting a higher concentration of particles close to its front. This is a common and yet unexplained phenomenon, which is usually believed to be the result of fluid convection in combination with particle size segregation. However, suspensions composed of uniformly sized particles also develop a granular front. Within a large rotating drum, a stationary recirculating avalanche is generated. The flowing material is a mixture of a viscoplastic fluid obtained from a kaolin-water dispersion with spherical ceramic particles denser than the fluid. The goal is to mimic the composition of many common granular-fluid materials, such as fresh concrete or debris flow. In these materials, granular and fluid phases have the natural tendency to separate due to particle settling. However, through the shearing caused by the rotation of the drum, a reorganization of the phases is induced, leading to the formation of a granular front. By tuning the particle concentration and the drum velocity, it is possible to control this phenomenon. The setting is reproduced in a numerical environment, where the fluid is solved by a lattice-Boltzmann method, and the particles are explicitly represented using the discrete element method. The simulations confirm the findings of the experiments, and provide insight into the internal mechanisms. Comparing the time scale of particle settling with the one of particle recirculation, a nondimensional number is defined, and is found to be effective in predicting the formation of a granular front.
NASA Astrophysics Data System (ADS)
Freydier, Perrine; Chambon, Guillaume; Naaim, Mohamed
2016-04-01
Rheological studies concerning natural muddy debris flows have shown that these materials can be modelled as non-Newtonian viscoplastic fluids. These complex flows are generally represented using models based on a depth-integrated approach (Shallow Water) that take into account closure terms depending on the shape of the velocity profile. But to date, there is poor knowledge about the shape of velocity profiles and the position of the interface between sheared and unsheared regions (plug) in these flows, especially in the vicinity of the front. In this research, the internal dynamics of a free-surface viscoplastic flow down an inclined channel is investigated and compared to the predictions of a Shallow Water model based on the lubrication approximation. Experiments are conducted in an inclined channel whose bottom is constituted by an upward-moving conveyor belt with controlled velocity, which allows generating and observing gravity-driven stationary surges in the laboratory frame. Carbopol microgel has been used as a homogeneous and transparent viscoplastic fluid. High-resolution measurements of velocity field is performed through optical velocimetry techniques both in the uniform zone and within the front zone where flow thickness is variable and where recirculation takes place. Specific analyses have been developed to determine the position of the plug within the surge. Flow height is accessible through image processing and ultrasonic sensors. Sufficiently far from the front, experimental results are shown to be in good agreement with theoretical predictions regarding the velocity profiles and the flow height evolution. In the vicinity of the front, however, analysis of measured velocity profiles shows an evolution of the plug different from that predicted by lubrication approximation. Accordingly, the free surface shape also deviates from the predictions of the classical Shallow Water model. These results highlight the necessity to take into account higher
Detection of liquid-metal, free-surface flow using the DLP measurement technique
NASA Astrophysics Data System (ADS)
Hillenbrand, M. P.; Stieglitz, R.; Neitzel, G. P.
2012-01-01
Novel accelerator applications favor free-surface liquid-metal flows, in which the liquid acts both as a target producing secondary particles but also to remove efficiently the heat deposited. A crucial aspect for the operation is the continuous monitoring of both shape and position of the liquid's surface. This demands, in a nuclear environment, a non-intrusive measurement technique with high temporal and spatial resolution. In this context, the double-layer projection (DLP) technique based on geometric optics has been developed, allowing one to detect either point-wise or area-wise the shape and position of the nearly totally reflecting liquid-metal surface. The DLP technique employs a laser beam projected through a coplanar glass plate to the surface from which it is reflected to the plate again. Beam locations captured by means of a camera permit the position and shape of the surface to be reconstructed. The parameters affecting the resolution and performance of the DLP technique are discussed. Additionally, validation studies using static and moving objects of pre-defined shape are conducted, exhibiting spatial and temporal resolutions of 300 μm and 100 Hz, respectively. Finally, the DLP system is applied to perform measurements of a circular hydraulic jump (CHJ) in a liquid metal. The DLP system has proved the capability to measure the jump both qualitatively and quantitatively. Additionally, the experiments identified, at high Reynolds numbers, the existence of a two-step jump. The analysis of spectral data of the DLP surface measurements shows clearly that, at the outer radius, gravity waves occur. Also, contributions from the pump oscillations were found, demonstrating the high performance of the DLP system.
A new numerical framework to simulate viscoelastic free-surface flows with the finite-volume method
NASA Astrophysics Data System (ADS)
Comminal, R.; Spangenberg, J.; Hattel, J. H.
2015-04-01
A new method for the simulation of 2D viscoelastic flow is presented. Numerical stability is obtained by the logarithmic-conformation change of variable, and a fully-implicit pure-streamfunction flow formulation, without use of any artificial diffusion. As opposed to other simulation results, our calculations predict a hydrodynamic instability in the 4:1 contraction geometry at a Weissenberg number of order 4. This new result is in qualitative agreement with the prediction of a non-linear subcritical elastic instability in Poiseuille flow. Our viscoelastic flow solver is coupled with a volume-of-fluid solver in order to predict free- surfaces in extrusion.
On the solution of transient free-surface flow problems in porous media by a fixed-domain method
NASA Astrophysics Data System (ADS)
Lin, Chin-Shing; Bruch, John C.; Sloss, James M.; Comincioli, Valeriano
1984-03-01
The numerical model presented here is the solution to the exact initial-boundary-value problems arising in recharge fluid flow through porous media having a free surface. Since the problem is nonlinear and includes a moving boundary, a numerical solution is obtained by using a fixed-domain approach. The Baiocchi transformation and method are used to develop a boundary-value problem which is then solved by an iterative method of successive over-relaxation type. Transient free-surface seepage through a two-dimensional dam with accretion is presented as an example problem. The effects of hydraulic conductivity, specific storativity and accretion on the seepage flow are studied. Two modified cases are satisfactorily compared with published results.
Wang, Xiao-dan; Zhai, Zhen-hua; Zhao, Shuang; Li, Rong-qi; Ma, Wen-lin; Li, Yan-hong
2009-01-01
To achieve the effects of artificial wetland on the bacterial diversity, the culturable bacteria and total cell counts of three wetland cells, including sewage pond (SP), free surface wetland (SF) and subsurface flow wetland (SSF), were investigated using the traditional culture-dependent approach and flow cytometry method, based on the detecting the water quality. The bacterial diversity and dominant groups were also compared by PCR-DGGE profiles and 16S rDNA library technique based on its V3 region. Results show that SF and SSF cells can remove the nutrients effectively, the highest removal ratio of COD, total nitrogen, and total phosphorus reach to 42.33%, 52.92% and 41.4%, respectively; The total microbes are increased continuously with the treatment by SF and SSF, and the culturable bacteria clones are decreased after treatment by SF, and increased after further train by SSF. The Shannon-Weaver index is increased to 3.2850 from 3.0819 while the water flowing through SF, but decreased to 3.0181 after flowing through SSF; The dominant groups in SP include Actinobacteria, Cyanobacteria and alpha-Proteobacteria, reach to 38%, 18% and 18%, respectively; but the most dominant bacteria is changed to beta-Proteobacteria with the ratio of 32% and 44%, after treatment by SF and SSF, respectively. Cytophagal Flexibacter/Bacteroides (CFB) phylum is also increased to 24% finally. Therefore, while the Cuihu Wetland removing the nutrients,the bacterial counts, diversity and dominant groups are also changed,some beneficial bacteria in beta-Proteobacteria and CFB phylum increased, and part of those deleterious bacteria in Actinobacteria and Cyanobacteria decreased. PMID:19353894
An explicit Lagrangian finite element method for free-surface weakly compressible flows
NASA Astrophysics Data System (ADS)
Cremonesi, Massimiliano; Meduri, Simone; Perego, Umberto; Frangi, Attilio
2016-07-01
In the present work, an explicit finite element approach to the solution of the Lagrangian formulation of the Navier-Stokes equations for weakly compressible fluids or fluid-like materials is investigated. The introduction of a small amount of compressibility is shown to allow for the formulation of a fast and robust explicit solver based on a particle finite element method. Newtonian and Non-Newtonian Bingham laws are considered. A barotropic equation of state completes the model relating pressure and density fields. The approach has been validated through comparison with experimental tests and numerical simulations of free surface fluid problems involving water and water-soil mixtures.
NASA Technical Reports Server (NTRS)
Lai, Chun-Liang; Greenberg, Paul S.; Chai, An-Ti
1988-01-01
To study thermocapillary flows in a two-dimensional thin liquid layer with heat fluxes imposed on the free surface experimentally, a long tray configuration was employed to simulate the infinite layer. The surface temperature distribution due to thermocapillary convection for different flow regimes was measured and compared with theorectical predictions. A short tray configuration was also employed to study the end wall effects (insulating or conducting). The results show that, for a strong convection flow with an insulating wall as the boundary, the surface temperature distribution became quite uniform. Consequently, the thermocapillary driving force was greatly reduced. On the other hand, a strong fluid motion always existed adjacent to the conducting wall because of the large surface temperature gradient near the wall.
NASA Technical Reports Server (NTRS)
Lai, Chun-Liang; Greenberg, Paul S.; Chai, An-Ti
1988-01-01
To study thermocapillary flows in a two-dimensional thin liquid layer with heat fluxes imposed on the free surface experimentally, a long tray configuration was employed to simulate the infinite layer. The surface temperature distribution due to thermocapillary convective for different flow regimes was measured and compared with theoretical predictions. A short tray configuration was also employed to study the end wall effects (insulating or conducting). The results show that for a strong convection flow with an insulating wall as the boundary the surface temperature distribution became quite uniform. Consequently, the thermocapillary driving force was greatly reduced. On the other hand, a strong fluid motion always existed adjacent to the conducting wall because of the large surface temperature gradient near the wall.
Baer, T.A.; Cairncross, R.A.; Rao, R.R.; Sackinger, P.A.; Schunk, P.R.
1999-01-29
To date, few researchers have solved three-dimensional free-surface problems with dynamic wetting lines. This paper extends the free-surface finite element method described in a companion paper [Cairncross, R.A., P.R. Schunk, T.A. Baer, P.A. Sackinger, R.R. Rao, "A finite element method for free surface flows of incompressible fluid in three dimensions, Part I: Boundary-Fitted mesh motion.", to be published (1998)] to handle dynamic wetting. A generalization of the technique used in two dimensional modeling to circumvent double-valued velocities at the wetting line, the so-called kinematic paradox, is presented for a wetting line in three dimensions. This approach requires the fluid velocity normal to the contact line to be zero, the fluid velocity tangent to the contact line to be equal to the tangential component of web velocity, and the fluid velocity into the web to be zero. In addition, slip is allowed in a narrow strip along the substrate surface near the dynamic contact line. For realistic wetting-line motion, a contact angle which varies with wetting speed is required because contact lines in three dimensions typically advance or recede a different rates depending upon location and/or have both advancing and receding portions. The theory is applied to capillary rise of static fluid in a corner, the initial motion of a Newtonian droplet down an inclined plane, and extrusion of a Newtonian fluid from a nozzle onto a moving substrate. The extrusion results are compared to experimental visualization. Subject Categories
NASA Astrophysics Data System (ADS)
Freydier, Perrine; Chambon, Guillaume; Naaim, Mohamed
2015-04-01
Debris flows constitute one of the most important natural hazards throughout the mountainous regions of the world, causing significant damages and economic losses. These mass are composed of particles of all sizes from clay to boulders suspended in a viscous fluid. An important goal resides in developing models that are able to accurately predict the hydraulic properties of debris flows. First, these flows are generally represented using models based on a momentum integral approach that consists in assuming a shallow flow and in depth averaging the local conservation equations. These models take into account closure terms depending on the shape of the velocity profile inside the flow. Second, the specific migration mechanisms of the suspended particles, which have a strong influence on the propagation of the surges, also depend on the internal dynamics within the flow. However, to date, few studies concerning the internal dynamics in particular in the vicinity of the front, of such flows have been carried out. The aim of this study is to document the internal dynamics in free-surface viscoplastic flows down an inclined channel. The rheological studies concerning natural muddy debris flows, rich in fine particles, have shown that these materials can be modeled, at least as a first approximation as non-Newtonian viscoplastic fluids. Experiments are conducted in an inclined channel whose bottom is constituted by an upward-moving conveyor belt with controlled velocity. Carbopol microgel has been used as a homogeneous transparent viscoplastic fluid. This experimental setup allows generating and monitoring stationary gravity-driven surges in the laboratory frame. We use PIV technique (Particle Image Velocimetry) to obtain velocity fields both in the uniform zone and within the front zone where flow thickness is variable and where recirculation takes place. Experimental velocity profiles and determination of plug position will be presented and compared to theoretical
Continuation in a parameter - Experience with viscous and free surface flows
NASA Technical Reports Server (NTRS)
Kheshgi, H. S.; Basaran, O. A.; Benner, R. E.; Kistler, S. F.; Scriven, L. E.
1983-01-01
The results of modifications in continuation methods applied to obtain solutions to the Navier-Stokes systems of equations for incompressible, two-dimensional, steady flows are reported. It is shown that parameter continuation permits prediction of accurate, initial estimates for iterative processing of nonlinear finite difference and finite element equations of motions. The new parameter steps are derived from values of the preceding parameter steps. The accuracy of the estimates is ensured with appropriate choices of the step size. The continuation predictor/iterative corrector is demonstrated to trace the branches of parameter space along which steady flow states are found, and techniques are available for tracing multiply branching paths. The techniques are applied to solving the Navier-Stokes equations for flow through a rotating square channel, the formation of a falling liquid curtain, and gyrostatic equilibria of rotating cylindrical drops.
Analysis of free-surface flows through energy considerations: Single-phase versus two-phase modeling
NASA Astrophysics Data System (ADS)
Marrone, Salvatore; Colagrossi, Andrea; Di Mascio, Andrea; Le Touzé, David
2016-05-01
The study of energetic free-surface flows is challenging because of the large range of interface scales involved due to multiple fragmentations and reconnections of the air-water interface with the formation of drops and bubbles. Because of their complexity the investigation of such phenomena through numerical simulation largely increased during recent years. Actually, in the last decades different numerical models have been developed to study these flows, especially in the context of particle methods. In the latter a single-phase approximation is usually adopted to reduce the computational costs and the model complexity. While it is well known that the role of air largely affects the local flow evolution, it is still not clear whether this single-phase approximation is able to predict global flow features like the evolution of the global mechanical energy dissipation. The present work is dedicated to this topic through the study of a selected problem simulated with both single-phase and two-phase models. It is shown that, interestingly, even though flow evolutions are different, energy evolutions can be similar when including or not the presence of air. This is remarkable since, in the problem considered, with the two-phase model about half of the energy is lost in the air phase while in the one-phase model the energy is mainly dissipated by cavity collapses.
Al-Baldawi, Israa Abdul Wahab; Abdullah, Siti Rozaimah Sheikh; Suja, Fatihah; Anuar, Nurina; Mushrifah, Idris
2013-11-30
Two types of flow system, free surface flow (FSF) and sub-surface flow (SSF), were examined to select a better way to remove total petroleum hydrocarbons (TPH) using diesel as a hydrocarbon model in a phytotoxicity test to Scirpus grossus. The removal efficiencies of TPH for the two flow systems were compared. Several wastewater parameters, including temperature (T, °C), dissolved oxygen (DO, mgL(-1)), oxidation-reduction potential (ORP, mV), and pH were recorded during the experimental runs. In addition, overall plant lengths, wet weights, and dry weights were also monitored. The phytotoxicity test using the bulrush plant S. grossus was run for 72 days with different diesel concentrations (1%, 2%, and 3%) (Vdiesel/Vwater). A comparison between the two flow systems showed that the SSF system was more efficient than the FSF system in removing TPH from the synthetic wastewater, with average removal efficiencies of 91.5% and 80.2%, respectively. The SSF system was able to tolerate higher diesel concentrations than was the FSF system. PMID:24113536
Marrone, Salvatore; Colagrossi, Andrea; Di Mascio, Andrea; Le Touzé, David
2016-05-01
The study of energetic free-surface flows is challenging because of the large range of interface scales involved due to multiple fragmentations and reconnections of the air-water interface with the formation of drops and bubbles. Because of their complexity the investigation of such phenomena through numerical simulation largely increased during recent years. Actually, in the last decades different numerical models have been developed to study these flows, especially in the context of particle methods. In the latter a single-phase approximation is usually adopted to reduce the computational costs and the model complexity. While it is well known that the role of air largely affects the local flow evolution, it is still not clear whether this single-phase approximation is able to predict global flow features like the evolution of the global mechanical energy dissipation. The present work is dedicated to this topic through the study of a selected problem simulated with both single-phase and two-phase models. It is shown that, interestingly, even though flow evolutions are different, energy evolutions can be similar when including or not the presence of air. This is remarkable since, in the problem considered, with the two-phase model about half of the energy is lost in the air phase while in the one-phase model the energy is mainly dissipated by cavity collapses. PMID:27300984
Instability on the Free Surface of Superfluid He-II Induced by a Steady Heat Flow in Bulk
NASA Astrophysics Data System (ADS)
Remizov, I. A.; Levchenko, A. A.; Mezhov-Deglin, L. P.
2016-06-01
We report observations of the onset of irregular motion on a free surface of superfluid He-II induced by a quasi-stationary heat flow in a rectangular container. The container open from the top is mounted inside an optical cell partly filled with superfluid He-II. Three holes in the container walls provide free circulation of the normal and superfluid components inside and outside the container. The results of measurements are discussed in terms of the Korshunov theory (Eurphys Lett 16:673, 1991; JETP Lett 75:423, 2002) of the Kelvin-Helmholtz instability on an initially flat He-II surface induced by a relative motion of superfluid and normal components of the liquid along the surface when the counterflow velocity exceeds the threshold value. The experimental data are qualitatively consistent with the theoretical predictions (Korshunov in JETP Lett 75:423, 2002) taking into account the finite viscosity of He-II.
Distribution of resistive body-force in curved free-surface flow
NASA Technical Reports Server (NTRS)
Sivakumaran, N. S.; Dressler, R. F.
1986-01-01
The customary procedure for including resistive effects in turbulent hydraulic and stratified atmospheric flows is to integrate the empirically-known boundary shears over the entire wetted boundary of a thin fluid slab. A resistive body-force is then assumed to exist everywhere in each slab to replace the boundary shearing force. For the classical Saint-Venant (1871) model, this body-force can be shown to have a constant distribution in the vertical direction, and therefore can be evaluated for use in the momentum differential equation. In the newer Dressler theory (1978), however, for unsteady flow over curved beds, it is proved here that a constant body-force distribution is not possible. Its variable distribution is determined as well as its magnitude for use in the curved-flow equations. This variable distribution acts to produce an equal resultant in every thin layer of fluid parallel to the bed in an angular wedge over the curved channel bed. The new curved-flow equations are therefore extended to include resistive effects.
Plane Stokes flow driven by capillarity on the free surfaces of a doubly-connected region
Hopper, R.W.
1991-01-07
The free creeping viscous incompressible plant flow a finite region, bounded by a simple smooth closed curve and driven solely by surface tension, was analyzed previously. The shape evolution was described in terms of a time-dependent mapping function z = {Omega}({zeta},t) of the unit circle, conformal on {vert bar}{zeta}{vert bar} {le} 1. An equation giving the time evolution of the mapping, typically in parametric form, was derived. The theory is here extended to doubly-connected regions mapped from an annulus. It is found that the conjectured parametric map must satisfy three conditions. The elementary problem of the collapse of a concentric circular annulus is easily recovered, but no non-trivial flows have been discovered. The theory is used to show that confocal elliptic ring does not collapse through a sequence of confocal elliptic rings. 5 refs., 1 fig.
Analysis of pipe flow with free surface. Part II. Theoretical analysis and experiment
NASA Astrophysics Data System (ADS)
Tanaka, Amane; Takaki, Ryuji
1994-05-01
Flow field near the front of an incompressible viscous fluid pushed into a circular pipe is analyzed theoretically and observed experimentally. In the theory, an approximated stream function for a steady state near the axis of the pipe is obtained by use of the Stokes equation. In the experiment, the shape of the surface was observed by a video camera. The theoretical velocity profile and the surface shape near the axis coincide with those from computation (Part I) and experiment.
Interfacial stress balances in structured continua and free surface flows in ferrofluids
NASA Astrophysics Data System (ADS)
Chaves, Arlex; Rinaldi, Carlos
2014-04-01
Interfacial linear and internal angular momentum balances are obtained for a structured continuum and for the special case of a ferrofluid, a suspension of magnetic nanoparticles in a Newtonian fluid. The interfacial balance equations account for the effects of surface tension and surface tension gradient, magnetic surface excess forces, antisymmetric stresses, and couple stresses in driving interfacial flows in ferrofluids. Application of the interfacial balance equations is illustrated by obtaining analytical expressions for the translational and spin velocity profiles in a thin film of ferrofluid on an infinite flat plate when a rotating magnetic field is applied with axis of rotation parallel to the ferrofluid/air interface. The cases of zero and non-zero spin viscosity are considered for small applied magnetic field amplitude. Expressions for the maximum translational velocity, slope of the translational velocity profile at the ferrofluid/air interface, and volumetric flow rate are obtained and their use to test the relevance of spin viscosity and couple stresses in the flow situation under consideration is discussed.
Simulation of the 3D viscoelastic free surface flow by a parallel corrected particle scheme
NASA Astrophysics Data System (ADS)
Jin-Lian, Ren; Tao, Jiang
2016-02-01
In this work, the behavior of the three-dimensional (3D) jet coiling based on the viscoelastic Oldroyd-B model is investigated by a corrected particle scheme, which is named the smoothed particle hydrodynamics with corrected symmetric kernel gradient and shifting particle technique (SPH_CS_SP) method. The accuracy and stability of SPH_CS_SP method is first tested by solving Poiseuille flow and Taylor-Green flow. Then the capacity for the SPH_CS_SP method to solve the viscoelastic fluid is verified by the polymer flow through a periodic array of cylinders. Moreover, the convergence of the SPH_CS_SP method is also investigated. Finally, the proposed method is further applied to the 3D viscoelastic jet coiling problem, and the influences of macroscopic parameters on the jet coiling are discussed. The numerical results show that the SPH_CS_SP method has higher accuracy and better stability than the traditional SPH method and other corrected SPH method, and can improve the tensile instability. Project supported by the Natural Science Foundation of Jiangsu Province, China (Grant Nos. BK20130436 and BK20150436) and the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province, China (Grant No. 15KJB110025).
Interfacial stress balances in structured continua and free surface flows in ferrofluids
Chaves, Arlex; Rinaldi, Carlos
2014-04-15
Interfacial linear and internal angular momentum balances are obtained for a structured continuum and for the special case of a ferrofluid, a suspension of magnetic nanoparticles in a Newtonian fluid. The interfacial balance equations account for the effects of surface tension and surface tension gradient, magnetic surface excess forces, antisymmetric stresses, and couple stresses in driving interfacial flows in ferrofluids. Application of the interfacial balance equations is illustrated by obtaining analytical expressions for the translational and spin velocity profiles in a thin film of ferrofluid on an infinite flat plate when a rotating magnetic field is applied with axis of rotation parallel to the ferrofluid/air interface. The cases of zero and non-zero spin viscosity are considered for small applied magnetic field amplitude. Expressions for the maximum translational velocity, slope of the translational velocity profile at the ferrofluid/air interface, and volumetric flow rate are obtained and their use to test the relevance of spin viscosity and couple stresses in the flow situation under consideration is discussed.
SPH-DCDEM model for arbitrary geometries in free surface solid-fluid flows
NASA Astrophysics Data System (ADS)
Canelas, Ricardo B.; Crespo, Alejandro J. C.; Domínguez, Jose M.; Ferreira, Rui M. L.; Gómez-Gesteira, Moncho
2016-05-01
A unified discretization of rigid solids and fluids is introduced, allowing for resolved simulations of fluid-solid phases within a meshless framework. The numerical solution, attained by Smoothed Particle Hydrodynamics (SPH) and a variation of Discrete Element Method (DEM), the Distributed Contact Discrete Element Method (DCDEM) discretization, is achieved by directly considering solid-solid and solid-fluid interactions. The novelty of the work is centred on the generalization of the coupling of the DEM and SPH methodologies for resolved simulations, allowing for state-of-the-art contact mechanics theories to be used in arbitrary geometries, while fluid to solid and vice versa momentum transfers are accurately described. The methods are introduced, analysed and discussed. Initial validations on the DCDEM and the fluid coupling are presented, drawing from test cases in the literature. An experimental campaign serves as a validation point for complex, large scale solid-fluid flows, where a set of blocks in several configurations is subjected to a dam-break wave. Blocks are tracked and positions are then compared between experimental data and the numerical solutions. A Particle Image Velocimetry (PIV) technique allows for the quantification of the flow field and direct comparison with numerical data. The results show that the model is accurate and is capable of treating highly complex interactions, such as transport of debris or hydrodynamic actions on structures, if relevant scales are reproduced.
A new three-dimensional terrain-following tidal model of free-surface flows
NASA Astrophysics Data System (ADS)
Lu, Fuqiang; Zhang, Zhuo; Song, Zhiyao; Yue, Songshan; Wen, Yongning
2015-12-01
A three-dimensional hydrodynamic model is presented which combines a terrain-following vertical coordinate with a horizontally orthogonal curvilinear coordinate system to fit the complex bottom topography and coastlines near estuaries, continental shelves, and harbors. To solve the governing equations more efficiently, we improve the alternating direction implicit method, which is extensively used in the numerical modeling of horizontal two-dimensional shallow water equations, and extend it to a three-dimensional tidal model with relatively little computational effort. Through several test cases and realistic applications, as presented in the paper, it can be demonstrated that the model is capable of simulating the periodic to-and-fro currents, wind-driven flow, Ekman spirals, and tidal currents in the near-shore region.
Asymptotic and numerical analysis of free surface flows: Lump solitons and wave turbulence
NASA Astrophysics Data System (ADS)
Berger, Kurt Michael
Three-dimensional solitary waves or lump solitons are known to be solutions to the Kadomtsev-Petviashvili I (KP-I) equation, which models small-amplitude shallow-water waves when the Bond number is greater than ⅓. Recently, Pego and Quintero presented a proof of the existence of such waves for the Benney-Luke equation with surface tension. Here we establish an explicit connection between the lump solitons of these two equations and numerically compute the Benney-Luke lump solitons and their speed-amplitude relation. Furthermore, we numerically collide two Benney-Luke lump solitons to illustrate their soliton wave character. Finally, we study the flow over an obstacle near the linear shallow-water speed and show that three-dimensional lump solitons are periodically generated. In the second part of this dissertation, our goal is to study, numerically, the statistics of a large number of interacting finite-depth gravity surface waves. The weak- or wave-turbulence problem consists of finding statistical states with constant flux of energy in wavenumber space. These states are obtained by forcing and dissipating the conservative water wave problem at disparate scales and predicting the spectrum, often as a Kolmogorov-like power law, at intermediate scales. Majda, McLaughlin, and Tabak started the numerical investigation of the predictions of weak turbulence theory using a nonlinear dispersive NLS model equation. Here we investigate wave turbulence in a manner similar to Majda et al, but for an equation modeling gravity water waves, which is the original context in which the problem was posed. We perform long time computations on the one-dimensional, finite-depth Benney-Luke equation and compute various statistical quantities of interest. To validate this model, we first show, analytically and numerically, that Benney-Luke equations correctly predict the main deterministic aspects of resonant gravity wave interactions: resonant quartets, Benjamin-Feir type wave
Plane stokes flow driven by capillarity on a free surface: Overview and status
Hopper, R.W.
1991-10-25
This article summarizes some recent work on the creeping viscous incompressible plane flow in a region, bounded by a simple smooth closed curve and driven solely by surface tension. Such problems are self-contained in that the applied tractions are intrinsic in the geometry. The objective is to determine exactly the time evolution of the shape of the region. The problems are fundamentally nonlinear due to the large changes in shape, and it is emphasized that no mathematical approximations are made. For finite regions, the shape in the complex z-plane is described in terms of a time-dependent conformal mapping function {Omega}({zeta}, t) on the fixed region {vert bar} {zeta} {vert bar} {le} 1 of the complex {zeta}-plane. In practice, it has been necessary to conjecture an explicit parametric form {Omega}({zeta};{alpha}{sub 1}(t), {alpha}{sub 2}(t), {hor ellipsis}) whose validity must be verified using the shape- evolution equation. Suitable parameterizations are not always obvious. When the conjectured form holds and the equations can be solved, the evolution of the shape with time is obtained in simple, exact and closed form except for the value of the time, which requires a mere quadrature. Internal velocity and stress fields can be obtained.
Plane stokes flow driven by capillarity on a free surface: Overview and status. Revision 2
Hopper, R.W.
1991-10-25
This article summarizes some recent work on the creeping viscous incompressible plane flow in a region, bounded by a simple smooth closed curve and driven solely by surface tension. Such problems are self-contained in that the applied tractions are intrinsic in the geometry. The objective is to determine exactly the time evolution of the shape of the region. The problems are fundamentally nonlinear due to the large changes in shape, and it is emphasized that no mathematical approximations are made. For finite regions, the shape in the complex z-plane is described in terms of a time-dependent conformal mapping function {Omega}({zeta}, t) on the fixed region {vert_bar} {zeta} {vert_bar} {le} 1 of the complex {zeta}-plane. In practice, it has been necessary to conjecture an explicit parametric form {Omega}[{zeta};{alpha}{sub 1}(t), {alpha}{sub 2}(t), {hor_ellipsis}] whose validity must be verified using the shape- evolution equation. Suitable parameterizations are not always obvious. When the conjectured form holds and the equations can be solved, the evolution of the shape with time is obtained in simple, exact and closed form except for the value of the time, which requires a mere quadrature. Internal velocity and stress fields can be obtained.
NASA Astrophysics Data System (ADS)
Zhong, Qiang; Chen, Qigang; Wang, Hao; Li, Danxun; Wang, Xingkui
2016-05-01
Long streamwise-elongated high- and low-speed streaks are repeatedly observed near the free surface in open channel flows in natural rivers and lab experiments. Super-streamwise vortex model has been proposed to explain this widespread phenomenon for quite some time. However, statistical evidence of the existence of the super-streamwise vortices as one type of coherent structures is still insufficient. Correlation and proper orthogonal decomposition (POD) analysis based on PIV experimental data in the streamwise-spanwise plane near the free surface in a smooth open channel flow are employed to investigate this topic. Correlation analysis revealed that the streaky structures appear frequently near the free surface and their occurrence probability at any spanwise position is equal. The spanwise velocity fluctuation usually flows from low-speed streaks toward high-speed streaks. The average spanwise width and spacing between neighboring low (or high) speed streaks are approximately h and 2h respectively. POD analysis reveals that there are streaks with different spanwise width in the instantaneous flow fields. Typical streamwise rotational movement can be sketched out directly based on the results from statistical analyses. Point-by-point analysis indicates that this pattern is consistent everywhere in the measurement window and is without any inhomogeneity in the spanwise direction, which reveals the essential difference between coherent structures and secondary flow cells. The pattern found by statistical analysis is consistent with the notion that the super-streamwise vortices exist universally as one type of coherent structure in open channel flows.
NASA Technical Reports Server (NTRS)
Preiswerk, Ernst
1940-01-01
The application is treated in sufficient detail to facilitate as much as possible its application by the engineer who is less familiar with the subject. The present work was undertaken with two objects in view. In the first place, it is considered as a contribution to the water analogy of gas flows, and secondly, a large portion is devoted to the general theory of the two-dimensional supersonic flows.
NASA Astrophysics Data System (ADS)
Furuichi, Mikito; May, Dave A.
2015-07-01
We present implicit time integration schemes suitable for modeling free surface Stokes flow dynamics with marker in cell (MIC) based spatial discretization. Our target is for example thermal convection surrounded by deformable surface boundaries to simulate the long term planetary formation process. The numerical system becomes stiff when the dynamical balancing time scale for the increasing/decreasing load by surface deformation is very short compared with the time scale associated with thermal convection. Any explicit time integration scheme will require very small time steps; otherwise, serious numerical oscillation (spurious solutions) will occur. The implicit time integration scheme possesses a wider stability region than the explicit method; therefore, it is suitable for stiff problems. To investigate an efficient solution method for the stiff Stokes flow system, we apply first (backward Euler (BE)) and second order (trapezoidal method (TR) and trapezoidal rule-backward difference formula (TR-BDF2)) accurate implicit methods for the MIC solution scheme. The introduction of implicit time integration schemes results in nonlinear systems of equations. We utilize a Jacobian free Newton Krylov (JFNK) based Newton framework to solve the resulting nonlinear equations. In this work we also investigate two efficient implicit solution strategies to reduce the computational cost when solving stiff nonlinear systems. The two methods differ in how the advective term in the material transport evolution equation is treated. We refer to the method that employs Lagrangian update as "fully implicit" (Imp), whilst the method that employs Eulerian update is referred to as "semi-implicit" (SImp). Using a finite difference (FD) method, we have performed a series of numerical experiments which clarify the accuracy of solutions and trade-off between the computational cost associated with the nonlinear solver and time step size. In comparison with the general explicit Euler method
NASA Astrophysics Data System (ADS)
Tomé, M. F.; Bertoco, J.; Oishi, C. M.; Araujo, M. S. B.; Cruz, D.; Pinho, F. T.; Vynnycky, M.
2016-04-01
This work is concerned with the numerical solution of the K-BKZ integral constitutive equation for two-dimensional time-dependent free surface flows. The numerical method proposed herein is a finite difference technique for simulating flows possessing moving surfaces that can interact with solid walls. The main characteristics of the methodology employed are: the momentum and mass conservation equations are solved by an implicit method; the pressure boundary condition on the free surface is implicitly coupled with the Poisson equation for obtaining the pressure field from mass conservation; a novel scheme for defining the past times t‧ is employed; the Finger tensor is calculated by the deformation fields method and is advanced in time by a second-order Runge-Kutta method. This new technique is verified by solving shear and uniaxial elongational flows. Furthermore, an analytic solution for fully developed channel flow is obtained that is employed in the verification and assessment of convergence with mesh refinement of the numerical solution. For free surface flows, the assessment of convergence with mesh refinement relies on a jet impinging on a rigid surface and a comparison of the simulation of a extrudate swell problem studied by Mitsoulis (2010) [44] was performed. Finally, the new code is used to investigate in detail the jet buckling phenomenon of K-BKZ fluids.
NASA Technical Reports Server (NTRS)
Preiswerk, Ernst
1940-01-01
In this paper an introduction to shock polar diagrams is given which then leads into an examination of water depths in hydraulic jumps. Energy loss during these jumps is considered along with an extended look at elementary solutions of flow. An experimental test set-up is described and the results presented.
NASA Astrophysics Data System (ADS)
Adjali, Saadia; Belkadi, Mustapha; Aounallah, Mohammed; Imine, Omar
2015-05-01
Accurate simulation of turbulent free surface flows around surface ships has a central role in the optimal design of such naval vessels. The flow problem to be simulated is rich in complexity and poses many modeling challenges because of the existence of breaking waves around the ship hull, and because of the interaction of the two-phase flow with the turbulent boundary layer. In this paper, our goal is to estimate the lift and drag coefficients for NACA 0012 of hydrofoil advancing in calm water under steady conditions with free surface and emerged NACA 0015. The commercial CFD software FLUENT version 14 is used for the computations in the present study. The calculated grid is established using the code computer GAMBIT 2.3.26.The shear stress k-ωSST model is used for turbulence modeling and the volume of fluid technique is employed to simulate the free-surface motion. In this computation, the second order upwind scheme is used for discretizing the convection terms in the momentum transport equations, the Modified HRIC scheme for VOF discretisation. The results obtained compare well with the experimental data.
Lithium-Based High Energy Density Flow Batteries
NASA Technical Reports Server (NTRS)
Bugga, Ratnakumar V. (Inventor); West, William C. (Inventor); Kindler, Andrew (Inventor); Smart, Marshall C. (Inventor)
2014-01-01
Systems and methods in accordance with embodiments of the invention implement a lithium-based high energy density flow battery. In one embodiment, a lithium-based high energy density flow battery includes a first anodic conductive solution that includes a lithium polyaromatic hydrocarbon complex dissolved in a solvent, a second cathodic conductive solution that includes a cathodic complex dissolved in a solvent, a solid lithium ion conductor disposed so as to separate the first solution from the second solution, such that the first conductive solution, the second conductive solution, and the solid lithium ionic conductor define a circuit, where when the circuit is closed, lithium from the lithium polyaromatic hydrocarbon complex in the first conductive solution dissociates from the lithium polyaromatic hydrocarbon complex, migrates through the solid lithium ionic conductor, and associates with the cathodic complex of the second conductive solution, and a current is generated.
Reserve, flowing electrolyte, high rate lithium battery
NASA Astrophysics Data System (ADS)
Puskar, M.; Harris, P.
Flowing electrolyte Li/SOCl2 tests in single cell and multicell bipolar fixtures have been conducted, and measurements are presented for electrolyte flow rates, inlet and outlet temperatures, fixture temperatures at several points, and the pressure drop across the fixture. Reserve lithium batteries with flowing thionyl-chloride electrolytes are found to be capable of very high energy densities with usable voltages and capacities at current densities as high as 500 mA/sq cm. At this current density, a battery stack 10 inches in diameter is shown to produce over 60 kW of power while maintaining a safe operating temperature.
NASA Astrophysics Data System (ADS)
Talmage, Gita; Walker, John S.; Brown, Samuel H.; Sondergaard, Neal A.; Burt, Patricia E.
1990-11-01
Fully developed, viscous liquid-metal velocity profiles and induced magnetic field contours were studied for Hartmann numbers of M=2 and 10 and for different load currents for a particular rectangular channel configuration (two-dimensional Couette flow). The rectangular channel was assumed to have a homogeneous external (axial) magnetic field parallel to the moving, perfectly conducting top wall and the stationary, perfectly conducting bottom wall. The two stationary side walls were also perfect conductors. The small gap between the moving wall and each side wall was an insulating, free surface. The method of weighted residuals was used to obtain truncated series solutions for the variables of interest. The heavy load currents across the channel were obtained by simulating an external potential to the conducting moving wall. The load currents in each case were opposed by the induced electric field. Since there is no pressure gradient, the flow along the channel is driven by the viscous effects of the moving wall and the Lorentz body force and is retarded by the stationary walls. In the case where no load current is applied across the channel, the current circulates in the channel. The circulation is driven by the generator that is due to the axial variation of velocity in an axial magnetic field. The numerical results presented show that the radial gap and the free surface region represent electrical resistances in parallel between the perfectly conducting stationary wall and the perfectly conducting moving wall. The numerical results also show that the resistance of the radial gap increases as M2 while that of the free surface increases by M or M1/2. Thus, as M increases, the division of current shifts to the free surface region and the current density in the radial gap decreases as M-1. The theoretical magnetohydrodynamic model presented here was developed to provide numerical parameters to help in the design of liquid-metal current collectors. Numerical results
Lithium-Polysulfide Flow Battery Demonstration
Zheng, Wesley
2014-07-16
In this video, Stanford graduate student Wesley Zheng demonstrates the new low-cost, long-lived flow battery he helped create. The researchers created this miniature system using simple glassware. Adding a lithium polysulfide solution to the flask immediately produces electricity that lights an LED. A utility version of the new battery would be scaled up to store many megawatt-hours of energy.
Lithium-Polysulfide Flow Battery Demonstration
Zheng, Wesley
2014-06-30
In this video, Stanford graduate student Wesley Zheng demonstrates the new low-cost, long-lived flow battery he helped create. The researchers created this miniature system using simple glassware. Adding a lithium polysulfide solution to the flask immediately produces electricity that lights an LED. A utility version of the new battery would be scaled up to store many megawatt-hours of energy.
On the effect of standard PFEM remeshing on volume conservation in free-surface fluid flow problems
NASA Astrophysics Data System (ADS)
Franci, Alessandro; Cremonesi, Massimiliano
2016-07-01
The aim of this work is to analyze the remeshing procedure used in the particle finite element method (PFEM) and to investigate how this operation may affect the numerical results. The PFEM remeshing algorithm combines the Delaunay triangulation and the Alpha Shape method to guarantee a good quality of the Lagrangian mesh also in large deformation processes. However, this strategy may lead to local variations of the topology that may cause an artificial change of the global volume. The issue of volume conservation is here studied in detail. An accurate description of all the situations that may induce a volume variation during the PFEM regeneration of the mesh is provided. Moreover, the crucial role of the parameter α used in the Alpha Shape method is highlighted and a range of values of α for which the differences between the numerical results are negligible, is found. Furthermore, it is shown that the variation of volume induced by the remeshing reduces by refining the mesh. This check of convergence is of paramount importance for the reliability of the PFEM. The study is carried out for 2D free-surface fluid dynamics problems, however the conclusions can be extended to 3D and to all those problems characterized by significant variations of internal and external boundaries.
NASA Astrophysics Data System (ADS)
Talmage, Gita; Walker, John S.; Brown, Samuel H.; Sondergaard, Neal A.; Burt, Patricia E.
1989-11-01
Fully developed, viscous liquid-metal velocity profiles and induced magnetic field contours were studied for Hartmann number of M=2 and 10 and for different load currents for a particular rectangular channel configuration (two-dimensional Couette flow). The rectangular channel was assumed to have a homogeneous external (axial) magnetic field parallel to the moving, perfectly conducting top wall and the stationary, perfectly conducting bottom wall. The two stationary side walls were also perfect conductors. The small gap between the moving wall and each side wall was an insulating, free surface. The method of weighted residuals was used to obtain truncated series solutions for the variables of interest. The heavy load currents across the channel were obtained by simulating an external potential to the conducting moving wall. The load currents in each case were opposed by the induced electric field. Since there is no pressure gradient, the flow along the channel is driven by the viscous effects of the moving wall and the Lorentz body force and is retarded by the stationary walls. The circulation is driven by the generator that is due to the axial variation of velocity in an axial magnetic field. The numerical data presented show that the radial gap and the free surface region represent electrical resistances in parallel between the perfectly conducting stationary wall and the perfectly conducting moving wall.
NASA Astrophysics Data System (ADS)
Wan, Gang; Jin, Yong; Li, Haiyuan; Li, Baoming
2016-03-01
Active boundary layer flow control and boundary layer manipulation in the channel flow that was based on low temperature plasma were studied by means of a lattice Boltzmann method. Two plasma actuators were placed in a row to obtain the influence rule of their separation distance on the velocity profile at three locations and maximum velocity in the flow field. Two plasma actuators were placed symmetrically inside a channel to examine the effect of channel height and voltage on the velocity profile and flow rate. It was found that the channel height controls the distribution of flow velocity, which affected the flow rate and its direction. Increasing plasma voltage had a negative effect on the flow rate due to the generation of a larger and stronger flow vortex.
Free-surface wave-induced separation
Zhang, Z.J.; Stern, F.
1996-09-01
Free-surface wave-induced separation is studied for a surface-piercing NACA 0024 foil over a range of Froude numbers (0, .2, .37, .55) through computational fluid dynamics of the unsteady Reynolds-averaged Navier-Stokes and the continuity equations with the Baldwin-Lomax turbulence model, exact nonlinear kinematic and approximate dynamic free-surface boundary conditions, and a body/free-surface conforming grid. The flow conditions and uncertainty analysis are discussed. A topological rule for a surface-piercing body is derived and verified. Steady-flow results are presented and analyzed with regard to the wave and viscous flow and the nature of the separation.
NASA Astrophysics Data System (ADS)
Yamamoto, T.; Takagi, Y.; Okano, Y.; Dost, S.
2016-03-01
NASA astronaut Pettit has conducted thermocapillary flow experiments in water films suspended in a solid ring onboard the International Space Station (ISS) in 2003 and 2011. In one of these experiments, an oscillatory thermocapillary flow was observed. The developed flow broke its symmetry along the centerline of the film. To the best of our knowledge, there are no studies on such oscillatory thermocapillary flows in thin films, and the flow-mechanism giving rise to such oscillatory flows is also not well understood. In order to shed light on the subject, we have carried out a numerical simulation study. The simulation results have shown that the water film geometry (film surface shape; being concave) is an important parameter and give rise to three oscillatory flow structures in the film, namely, a hydrothermal wave developing near the heated section, a symmetric oscillatory flow due to temperature variations, and a symmetry breaking flow due to the hydrodynamic instability along the free boundary layer (mixing layer) and the development of the hydrothermal waves. Simulation results show that the symmetry-breaking phenomenon observed in the thin film experiment on the ISS can be explained by the hydrodynamic instability and the development of hydrothermal waves.
Lithium mass flow control for high power Lorentz Force Accelerators
NASA Astrophysics Data System (ADS)
Kodys, Andrea D.; Emsellem, Gregory; Cassady, Leonard D.; Polk, James E.; Choueiri, Edgar Y.
2001-02-01
A lithium feeding system has been developed to measure and control propellant flow for 30-200 kW Lithium Lorentz Force Accelerators (LiLFAs). The new, mechanically actuated, liquid lithium feed system has been designed and tested as a central component of a campaign to obtain basic data and establish scaling laws and performance relations for these thrusters. Calibration data are presented which demonstrate reliable and controllable feed of liquid lithium to the vaporizer hollow cathode of the thruster at flow rates between 10 and 120 mg/s. The ability to thermally track the liquid lithium through the system by the use of external temperature measurements is demonstrated. In addition, recent developments are presented in the establishment and successful testing of a lithium handling facility and safety procedures allowing for the in-house loading of the feed system and the neutralization, cleaning and disposal of up to 300 g of lithium. .
NASA Technical Reports Server (NTRS)
Faghri, Amir; Swanson, Theodore D.
1988-01-01
The results of a numerical computation and theoretical analysis are presented for the flow of a thin liquid film in the presence and absence of a gravitational body force. Five different flow systems were used. Also presented are the governing equations and boundary conditions for the situation of a thin liquid emanating from a pressure vessel; traveling along a horizontal plate with a constant initial height and uniform initial velocity; and traveling radially along a horizontal disk with a constant initial height and uniform initial velocity.
A low diffusive Lagrange-remap scheme for the simulation of violent ai-water free-surface flows
NASA Astrophysics Data System (ADS)
Bernard-Champmartin, Aude; De Vuyst, Florian
2014-10-01
In 2002, Després and Lagoutière [17] proposed a low-diffusive advection scheme for pure transport equation problems, which is particularly accurate for step-shaped solutions, and thus suited for interface tracking procedure by a color function. This has been extended by Kokh and Lagoutière [28] in the context of compressible multifluid flows using a five-equation model. In this paper, we explore a simplified variant approach for gas-liquid three-equation models. The Eulerian numerical scheme has two ingredients: a robust remapped Lagrange solver for the solution of the volume-averaged equations, and a low diffusive compressive scheme for the advection of the gas mass fraction. Numerical experiments show the performance of the computational approach on various flow reference problems: dam break, sloshing of a tank filled with water, wate-water impact and finally a case of Rayleigh-Taylor instability. One of the advantages of the present interface capturing solver is its natural implementation on parallel processors or computers. wave formation and wave breaking; wall wave impacts, local pressure peaks and pressure loadings; formation of air pockets; ejection, fragmentation of liquid droplets; Archimedes buoyancy effect with rising of bubbles and fall of droplets; effects of gas compressibility inducing a gas-to-liquid response by a pressure wave, etc. In this paper, we consider immiscible gas-liquid two-phase flow problems. The strong ratio of mass density between gas and liquid (typically 1:1000) is known to be a source of numerical stiffness and numerical instability. Therefore robust computational approaches supporting high density ratio have to be considered. Among the family of conservative Finite Volume methods (FVM), the Lagrange-remapped solvers (see e.g. [42,45,6,4,25,2]) provide both robustness and stability with achievement of mathematical properties of positiveness and entropy compatibility.Lagrange-remap numerical schemes (also referred to as Eule
A low diffusive Lagrange-remap scheme for the simulation of violent ai-water free-surface flows
NASA Astrophysics Data System (ADS)
Bernard-Champmartin, Aude; De Vuyst, Florian
2014-10-01
In 2002, Després and Lagoutière [17] proposed a low-diffusive advection scheme for pure transport equation problems, which is particularly accurate for step-shaped solutions, and thus suited for interface tracking procedure by a color function. This has been extended by Kokh and Lagoutière [28] in the context of compressible multifluid flows using a five-equation model. In this paper, we explore a simplified variant approach for gas-liquid three-equation models. The Eulerian numerical scheme has two ingredients: a robust remapped Lagrange solver for the solution of the volume-averaged equations, and a low diffusive compressive scheme for the advection of the gas mass fraction. Numerical experiments show the performance of the computational approach on various flow reference problems: dam break, sloshing of a tank filled with water, wate-water impact and finally a case of Rayleigh-Taylor instability. One of the advantages of the present interface capturing solver is its natural implementation on parallel processors or computers. wave formation and wave breaking; wall wave impacts, local pressure peaks and pressure loadings; formation of air pockets; ejection, fragmentation of liquid droplets; Archimedes buoyancy effect with rising of bubbles and fall of droplets; effects of gas compressibility inducing a gas-to-liquid response by a pressure wave, etc. In this paper, we consider immiscible gas-liquid two-phase flow problems. The strong ratio of mass density between gas and liquid (typically 1:1000) is known to be a source of numerical stiffness and numerical instability. Therefore robust computational approaches supporting high density ratio have to be considered. Among the family of conservative Finite Volume methods (FVM), the Lagrange-remapped solvers (see e.g. [42,45,6,4,25,2]) provide both robustness and stability with achievement of mathematical properties of positiveness and entropy compatibility.Lagrange-remap numerical schemes (also referred to as Eule
Free surface calculations in mantle convection
NASA Astrophysics Data System (ADS)
Rose, I.; Buffett, B. A.; Heister, T.
2015-12-01
Geodynamic simulations increasingly rely on simulations with a true free surface to investigate questions of dynamic topography, tectonic deformation, gravity perturbations, and global mantle convection. However, implementations of free surface boundary conditions have proven challenging from a standpoint of accuracy, robustness, and stability. In particular, free surfaces tend to suffer from sloshing instabilities, also known as the "drunken sailor" instability, which severely limit time step sizes. Several schemes have been proposed in the literature to deal with these instabilities. Here we analyze the problem of creeping viscous flow with a free surface and discuss the origin of these instabilities. We demonstrate their cause and how existing stabilization schemes work to damp them out. Our analysis is based on formulating a generalized eigenvalue problem for the relaxation spectra of the linearized free surface problem. We also propose a new scheme for removing instabilites from free surface calculations. It does not require modifications to the system matrix, nor additional variables, but is instead an explicit scheme based on nonstandard finite differences. It relies on a single stabilization parameter which may be identified with the smallest relaxation timescale of the free surface. We analyze the stability and accuracy of the nonstandard finite difference scheme, and describe its implementation in the open source mantle convection software Aspect. We also provide comparisons between the nonstandard finite difference scheme and the quasi-implicit scheme proposed by Kaus, Muhlhaus, and May (2010).
Fluid Dynamics with Free Surfaces
Energy Science and Technology Software Center (ESTSC)
1992-02-01
RIPPLE is a two-dimensional, transient, free surface incompressible fluid dynamics program. It allows multiple free surfaces with surface tension and wall adhesion forces and has a partial cell treatment which allows curved boundaries and interior obstacles.
NASA Astrophysics Data System (ADS)
Furuichi, M.; Nakagawa, T.; May, D.
2013-12-01
Stabilizing a numerical oscillation in free surface treatment is chagrining topic for a geodynamics simulation [e.g. Kaus et al. 2010, Duretz et al., 2011]. It is especially important for the Stokes flow simulation under the self-gravitating field based on 'Spherical Cartesian' method [Gerya et al., 2007], which is useful for simulating a long time scale dynamics of sinking metal rich materials to construct planetary core. The conventional explicit time stepping algorithm, which solves Stokes flow equation for a given material distribution at a previous time step, however has a difficulty for simulating dynamics such as a thermal convection, after the construction of layered structure in the planetary interior because of numerical oscillation. One effective approach for such numerically problematic behavior is an implicit treatment of advection term. In this study, three types of implicit strategy are discussed. First is the full implicit treatment with iterative non-linear solver which uses transported density by maker-in-cell method as nonlinear update. The maker-in-cell method is commonly used as low diffusive advection method, but is computationally expensive with makers to mesh interpolation. Second approach uses semi-Lagrangian method for nonlinear update instead of the maker-in-cell method to reduce computational cost. Third approach is to solve the Stokes flow equation combined with the linearized advection term in central-difference discretization to avoid the nonlinear update by the transport. In the second and third algorithms, physical value at the next time step is still transported by low diffusive maker-in-cell method. These three types of implicit method are examined by numerical experiment.
Han, L.H. Hu, X.Y. Adams, N.A.
2015-01-01
In this paper we present a scale separation approach for multi-scale modeling of free-surface and two-phase flows with complex interface evolution. By performing a stimulus-response operation on the level-set function representing the interface, separation of resolvable and non-resolvable interface scales is achieved efficiently. Uniform positive and negative shifts of the level-set function are used to determine non-resolvable interface structures. Non-resolved interface structures are separated from the resolved ones and can be treated by a mixing model or a Lagrangian-particle model in order to preserve mass. Resolved interface structures are treated by the conservative sharp-interface model. Since the proposed scale separation approach does not rely on topological information, unlike in previous work, it can be implemented in a straightforward fashion into a given level set based interface model. A number of two- and three-dimensional numerical tests demonstrate that the proposed method is able to cope with complex interface variations accurately and significantly increases robustness against underresolved interface structures.
NASA Astrophysics Data System (ADS)
Rosatti, Giorgio; Zugliani, Daniel
2015-03-01
In a two-phase free-surface flow, the transition from a mobile-bed condition to a fixed-bed one (and vice versa) occurs at a sharp interface across which the relevant system of partial differential equations changes abruptly. This leads to the possibility of conceiving a new type of Riemann Problem (RP), which we have called Composite Riemann Problem (CRP), where not only the initial constant values of the variables but also the system of equations change from left to right of a discontinuity. In this paper, we present a strategy for solving a CRP by reducing it to a standard RP of a single, composite system of equations. This can be obtained by combining the two original systems by means of a suitable weighting function, namely the erodibility variable, and the introduction of an appropriate differential equation for this quantity. In this way, the CRP problem can be analyzed theoretically with standard methods, and the features of the solutions can be clearly identified. In particular, a stationary contact wave is able to correctly describe the sharp transition between mobile- and fixed-bed conditions. A finite volume scheme based on the Multiple Averages Generalized Roe approach (Rosatti and Begnudelli (2013) [22]) was used to numerically solve the fixed-mobile CRP. Several test cases demonstrate the effectiveness, exact well balanceness and high accuracy of the scheme when applied to problems that fall within the physical range of applicability of the relevant mathematical model.
A free surface sharpening strategy using optimization method
NASA Astrophysics Data System (ADS)
Song, Hongchao; Ji, Lucheng; Tu, Shuangzhang
2015-11-01
VOF method which consists in transporting a discontinuous marker variable is widely used to capture the free surface in computational fluid dynamics. There is numerical dissipation in simulations involving the transport of the marker. Numerical dissipation makes the free surface lose its physical nature. A free surface sharpening strategy based on optimization method is presented in the paper. The strategy can keep the location of the free surface and local mass conservation at both time, and can also keep free surface in a constant width. It is independent on the types of solvers and meshes. Two famous cases were chosen for verifying the free surface sharpening strategy performance. Results show that the strategy has a very good performance on keeping local mass conservation. The efficiency of prediction of the free surface is improved by applying the strategy. Accurate modeling of flow details such as drops can also be captured by this method.
Free surface mixing with heat transfer
Raad, P.E.; Fadda, D.
1997-05-23
This study investigates the thermal mixing phenomenon as external rocking excitations are applied to a container partially filled with a thermally nonhomogeneous fluid. The aim of this work is to determine the effects of the rocking frequency on thermal mixing in free surface flows. Numerical results show that optimum mixing is achieved while rocking at the wave natural frequency of the tank where the free surface displacements are greatest. At this optimal frequency the strong mixing is observed to be highly concentrated in the area near the free surface. The interest is motivated by the desire to understand the fundamental physics as well as by the importance of current engineering applications in fields such as solar energy (collectors), environmental engineering (ponds, estuaries) aerospace (fuel tanks), electronics (thermal management), manufacturing (sand casting), and air conditioning (residential enclosures).
Simulation of Nematic Free Surfaces
NASA Astrophysics Data System (ADS)
de Miguel, Enrique; Martín del Río, Elvira
Molecular dynamics and Monte Carlo methods are applied to study the liquid free surfaces in model liquid crystals. The simulation results suggest that the attractive interactions promote parallel alignment of the molecules at the nematic free surface in the Gay-Berne model, in agreement with theoretical predictions. A change in the orientation from planar to homeotropic is observed and explained in terms of a competing effect between attractive and repulsive interactions. Finally, the simulation results give clear evidence that the hard-core repulsions favor homeotropic orientation at the nematic free surface, in agreement with most theories.
NASA Astrophysics Data System (ADS)
Malgarinos, Ilias; Nikolopoulos, Nikolaos; Gavaises, Manolis
2015-11-01
This study presents the implementation of an interface sharpening scheme on the basis of the Volume of Fluid (VOF) method, as well as its application in a number of theoretical and real cases usually modelled in literature. More specifically, the solution of an additional sharpening equation along with the standard VOF model equations is proposed, offering the advantage of "restraining" interface numerical diffusion, while also keeping a quite smooth induced velocity field around the interface. This sharpening equation is solved right after volume fraction advection; however a novel method for its coupling with the momentum equation has been applied in order to save computational time. The advantages of the proposed sharpening scheme lie on the facts that a) it is mass conservative thus its application does not have a negative impact on one of the most important benefits of VOF method and b) it can be used in coarser grids as now the suppression of the numerical diffusion is grid independent. The coupling of the solved equation with an adaptive local grid refinement technique is used for further decrease of computational time, while keeping high levels of accuracy at the area of maximum interest (interface). The numerical algorithm is initially tested against two theoretical benchmark cases for interface tracking methodologies followed by its validation for the case of a free-falling water droplet accelerated by gravity, as well as the normal liquid droplet impingement onto a flat substrate. Results indicate that the coupling of the interface sharpening equation with the HRIC discretization scheme used for volume fraction flux term, not only decreases the interface numerical diffusion, but also allows the induced velocity field to be less perturbed owed to spurious velocities across the liquid-gas interface. With the use of the proposed algorithmic flow path, coarser grids can replace finer ones at the slight expense of accuracy.
Effect of Energetic Plasma Flux on Flowing Liquid Lithium Surfaces
NASA Astrophysics Data System (ADS)
Kalathiparambil, Kishor; Jung, Soonwook; Christenson, Michael; Fiflis, Peter; Xu, Wenyu; Szott, Mathew; Ruzic, David
2014-10-01
An operational liquid lithium system with steady state flow driven by thermo-electric magneto-hydrodynamic force and capable of constantly refreshing the plasma exposed surface have been demonstrated at U of I. To evaluate the system performance in reactor relevant conditions, specifically to understand the effect of disruptive plasma events on the performance of the liquid metal PFCs, the setup was integrated to a pulsed plasma generator. A coaxial plasma generator drives the plasma towards a theta pinch which preferentially heats the ions, simulating ELM like flux, and the plasma is further guided towards the target chamber which houses the flowing lithium system. The effect of the incident flux is examined using diagnostic tools including triple Langmuir probe, calorimeter, rogowski coils, Ion energy analyzers, and fast frame spectral image acquisition with specific optical filters. The plasma have been well characterized and a density of ~1021 m-3, with electron temperature ~10 - 20 eV is measured, and final plasma velocities of 34 - 74 kms-1 have been observed. Calorimetric measurements using planar molybdenum targets indicate a maximum plasma energy (with 6 kV plasma gun and 20 kV theta pinch) of 0.08 MJm-2 with plasma divergence effects resulting in marginal reduction of 40 +/- 23 J in plasma energy. Further results from the other diagnostic tools, using the flowing lithium targets and the planar targets coated with lithium will be presented. DOE DE-SC0008587.
Free surface dynamics of nematic liquid crystal
NASA Astrophysics Data System (ADS)
Cummings, Linda; Kondic, Lou; Lam, Michael; Lin, Te-Sheng
2014-11-01
Spreading thin films of nematic liquid crystal (NLC) are known to behave very differently to those of isotropic fluids. The polar interactions of the rod-like molecules with each other, and the interactions with the underlying substrate, can lead to intricate patterns and instabilities that are not yet fully understood. The physics of a system even as simple as a film of NLC spreading slowly over a surface (inclined or horizontal) are remarkably complex: the outcome depends strongly on the details of the NLC's behavior at both the substrate and the free surface (so-called ``anchoring'' effects). We will present a dynamic flow model that takes careful account of such nematic-substrate and nematic-free surface interactions. We will present model simulations for several different flow scenarios that indicate the variety of behavior that can emerge. Spreading over a horizontal substrate may exhibit a range of unstable behavior. Flow down an incline also exhibits intriguing instabilities: in addition to the usual transverse fingering, instabilities can be manifested behind the flowing front in a manner reminiscent of Newtonian flow down an inverted substrate. NSF DMS-1211713.
Two-Dimensional Optical Measurement of Waves on Liquid Lithium Jet Simulating IFMIF Target Flow
Kazuhiro Itoh; Hiroyuki Koterazawa; Taro Itoh; Yutaka Kukita; Hiroo Kondo; Nobuo Yamaoka; Hiroshi Horiike; Mizuho Ida; Hideo Nakamura; Hiroo Nakamura; Takeo Muroga
2006-07-01
Waves on a liquid-lithium jet flow, simulating a proposed high-energy beam target design, have been measured using an optical technique based on specular reflection of a single laser beam on the jet surface. The stream-wise and spanwise fluctuations of the local free-surface slope were least-square fitted with a sinusoidal curve to makeup the signals lost due to the constriction in the optical arrangement. The waveform was estimated with an assumption that wave phase speed can be calculated using the dispersion relation for linear capillary-gravity waves. The direction of propagation on the jet surface was also evaluated so that the wave amplitudes, calculated by integral of slope angle signal, agree consistently in stream-wise and spanwise direction. These measurements and analyses show that the waves at the measurement location for a jet velocity of 1.2 m/s can best be represented by oblique waves with an inclination of 1.23 rad, a wavelength of 3.8 mm and a wave amplitude of about 0.05 mm. (authors)
Jaskula, B.W.
2011-01-01
In 2010, lithium consumption in the United States was estimated to have been about 1 kt (1,100 st) of contained lithium, a 23-percent decrease from 2009. The United States was estimated to be the fourth largest consumer of lithium. It remained the leading importer of lithium carbonate and the leading producer of value-added lithium materials. Only one company, Chemetall Foote Corp. (a subsidiary of Chemetall GmbH of Germany), produced lithium compounds from domestic resources. In 2010, world lithium consumption was estimated to have been about 21 kt (22,000 st) of lithium contained in minerals and compounds, a 12-percent increase from 2009.
Jaskula, B.W.
2010-01-01
In 2009, lithium consumption in the United States was estimated to have been about 1.2 kt (1,300 st) of contained lithium, a 40-percent decrease from 2008. The United States was estimated to be the fourth largest consumer of lithium, and remained the leading importer of lithium carbonate and the leading producer of value-added lithium materials. Only one company, Chemetall Foote Corp. (a subsidiary of Chemetall GmbH of Germany), produced lithium compounds from domestic resources. In 2009, world lithium consumption was estimated to have been about 18.7 kt (20,600 st) of lithium contained in minerals and compounds.
NASA Astrophysics Data System (ADS)
Caviedes-Voullième, Daniel; Juez, Carmelo; Murillo, Javier; García-Navarro, Pilar
2014-12-01
Avalanches, debris flows and other types of gravity-driven granular flows are a common hazard in mountainous regions. These regions often have human settlements in the lower parts of valleys, with human structures dangerously exposed to the destructive effects of these geophysical flows. Therefore a scientific effort has been made to understand, model and simulate geophysical granular flows. In order for computer models and simulations to be of predictive value they need to be validated under controlled, yet nature-like conditions. This work presents an experimental study of granular flow over a simplified mountain slope and valley topography. The experimental facility has a rough bed with very high slope at the upstream end and adverse slope on the downstream end, following a parabolic profile. Obstacles are present in the lower regions. Transient measurements of the moving granular surfaces were taken with a consumer-grade RGB-D sensor, providing transient 2D elevation fields around the obstacles. Three experimental configurations were tested, with semispheres of different diameters and a square dike obstacle. The experimental results are very consistent and repeatable. The quantitative, transient and two-dimensional data for all three experiments constitute excellent benchmarking tests for computational models, such as the one presented in a companion paper.
NASA Astrophysics Data System (ADS)
Juez, C.; Caviedes-Voullième, D.; Murillo, J.; García-Navarro, P.
2014-12-01
Dense granular flows are present in geophysics and in several industrial processes, which has lead to an increasing interest for the knowledge and understanding of the physics which govern their propagation. For this reason, a wide range of laboratory experiments on gravity-driven flows have been carried out during the last two decades. The present work is focused on geomorphological processes and, following previous work, a series of laboratory studies which constitute a further step in mimicking natural phenomena are described and simulated. Three situations are considered with some common properties: a two-dimensional configuration, variable slope of the topography and the presence of obstacles. The setup and measurement technique employed during the development of these experiments are deeply explained in the companion work. The first experiment is based on a single obstacle, the second one is performed against multiple obstacles and the third one studies the influence of a dike on which overtopping occurs. Due to the impact of the flow against the obstacles, fast moving shocks appear, and a variety of secondary waves emerge. In order to delve into the physics of these types of phenomena, a shock-capturing numerical scheme is used to simulate the cases. The suitability of the mathematical models employed in this work has been previously validated. Comparisons between computed and experimental data are presented for the three cases. The computed results show that the numerical tool is able to predict faithfully the overall behavior of this type of complex dense granular flow.
Krynkin, A; Horoshenkov, K V; Nichols, A; Tait, S J
2014-11-01
In this paper, the directivity of the airborne sound field scattered by a dynamically rough free flow surface in a flume is used to determine the mean roughness height for six hydraulic conditions in which the uniform depth of the turbulent flow. The nonlinear curve fitting method is used to minimize the error between the predicted directivity and directivity data. The data fitting algorithm is based on the averaged solution for the scattered sound pressure as a function of angle which is derived through the Kirchhoff integral and its approximations. This solution takes into account the directivity of the acoustic source. For the adopted source and receiver geometry and acoustic frequency it is shown that the contribution from the stationary phase point (single specular point on the rough surface) yields similar results to those which can be obtained through the full Kirchhoff's integral. The accuracy in the inverted mean roughness height is comparable to that achieved with an array of conductive wave probes. This method enables non-invasive estimation of the flow Reynolds number and uniform flow depth. PMID:25430137
Jaskula, B.W.
2012-01-01
In 2011, world lithium consumption was estimated to have been about 25 kt (25,000 st) of lithium contained in minerals and compounds, a 10-percent increase from 2010. U.S. consumption was estimated to have been about 2 kt (2,200 st) of contained lithium, a 100-percent increase from 2010. The United States was estimated to be the fourth-ranked consumer of lithium and remained the leading importer of lithium carbonate and the leading producer of value-added lithium materials. One company, Chemetall Foote Corp. (a subsidiary of Chemetall GmbH of Germany), produced lithium compounds from domestic brine resources near Silver Peak, NV.
Ober, J.A.
2006-01-01
In 2005, lithium consumption in the United States was at 2.5 kt of contained lithium, nearly 32% more than the estimate for 2004. World consumption was 14.1 kt of lithium contained in minerals and compounds in 2003. Exports from the US increased slightly compared with 2004. Due to strong demand for lithium compounds in 2005, both lithium carbonate plants in Chile were operating at or near capacity.
A study of free-surface electrohydrodynamics
Kaiser, K.L.
1989-01-01
Free-surface electrohydrodynamics (EHD) is a complex phenomena that is difficult to describe, explain, and predict. An analytical and numerical study is performed to assist in the understanding of free-surface EHD. Rayleigh-type critical charge relationships are derived for spherical droplets of perfect and imperfect insulative liquids. Also, mathematical models and a simulator are developed to analyze free-surface EHD in two dimensions. Several configurations are investigated with the simulator. Analytical results indicate that if charge is uniformly dispersed throughout a spherical droplet of perfectly insulating fluid, breakup is possible. If the dielectric constant of the liquid is approximately equal to that of the surrounding air, the charge limit is approximately 10% less than that of a conducting fluid. If ionizable impurities are in the fluid but no transfer of charge occurs between the electrode and the fluid, breakup can still occur. If the liquid's dielectric constant is large, the total charge on the electrode is not much different from the Rayleigh number. A unique donor-cell method was developed to approximate the charge flow equation, properly conserving charge in each cell. A pseudo-surface charge method was also developed. The results were shown to be good for low curvature, simple geometries. Simulation results indicate that the effects of inertia can be substantial in the rise of dielectrophoretic fluid between two electrified plates. The initial stages of the radial breakup of a charged conducting fluid cylinder consisted of wave-like oscillations along the surface. The shape of an insulating jet from a nozzle was seen to be highly dependent on the electric field configuration.
Code of Federal Regulations, 2014 CFR
2014-10-01
... 46 Shipping 1 2014-10-01 2014-10-01 false Free surface. 28.540 Section 28.540 Shipping COAST GUARD... Stability § 28.540 Free surface. (a) When doing the stability calculations required by this subpart, the... calculating the following— (1) For each type of consumable liquid, the maximum free surface effect of a...
Code of Federal Regulations, 2011 CFR
2011-10-01
... 46 Shipping 1 2011-10-01 2011-10-01 false Free surface. 28.540 Section 28.540 Shipping COAST GUARD... Stability § 28.540 Free surface. (a) When doing the stability calculations required by this subpart, the... calculating the following— (1) For each type of consumable liquid, the maximum free surface effect of a...
Code of Federal Regulations, 2010 CFR
2010-10-01
... 46 Shipping 1 2010-10-01 2010-10-01 false Free surface. 28.540 Section 28.540 Shipping COAST GUARD... Stability § 28.540 Free surface. (a) When doing the stability calculations required by this subpart, the... calculating the following— (1) For each type of consumable liquid, the maximum free surface effect of a...
Code of Federal Regulations, 2013 CFR
2013-10-01
... 46 Shipping 1 2013-10-01 2013-10-01 false Free surface. 28.540 Section 28.540 Shipping COAST GUARD... Stability § 28.540 Free surface. (a) When doing the stability calculations required by this subpart, the... calculating the following— (1) For each type of consumable liquid, the maximum free surface effect of a...
Code of Federal Regulations, 2012 CFR
2012-10-01
... 46 Shipping 1 2012-10-01 2012-10-01 false Free surface. 28.540 Section 28.540 Shipping COAST GUARD... Stability § 28.540 Free surface. (a) When doing the stability calculations required by this subpart, the... calculating the following— (1) For each type of consumable liquid, the maximum free surface effect of a...
Transient shear flow of model lithium lubricating greases
NASA Astrophysics Data System (ADS)
Delgado, M. A.; Franco, J. M.; Valencia, C.; Kuhn, E.; Gallegos, C.
2009-03-01
This paper deals with the analysis of the transient shear flow behavior of lithium lubricating greases differing in soap concentration and base oil viscosity. The shear-induced evolution of grease microstructure has been studied by means of stress-growth experiments. With this aim, different lubricating grease formulations were manufactured by modifying the concentration of lithium 12-hydroxystearate and the viscosity of the base oil, according to a RSM statistical design. Moreover, atomic force microscopy (AFM) observations were carried out. The transient stress response can be successfully described by the generalized Leider-Bird model based on two exponential terms. Different rheological parameters, related to both the elastic response and the structural breakdown of greases, have been analysed. In this sense, it has been found that the elastic properties of lithium lubricating greases were highly influenced by soap concentration and oil viscosity. The stress overshoot, τ max , depends linearly on both variables in the whole shear rate range studied, although the effect of base oil viscosity on this parameter is opposite at low and high shear rates. Special attention has been given to the first part of the stress-growth curve. In this sense, it can be deduced that the “yielding” energy density not only depends on grease composition, but also on shear rate. Moreover, an interesting asymptotic tendency has been found for both the “yielding” energy density and the stress overshoot by increasing shear rate. The asymptotic values of these parameters have been correlated to the friction coefficient obtained in a ball-disc tribometer.
Wei, Ze-Jun; Xie, Jian-Ping; Huang, Yu-Ming
2012-11-01
Many previous studies demonstrated that the performance of the subsurface constructed wetlands (SSCW) for wastewater treatment was superior to that of the free flow surface constructed wetlands (FFSCW). However, our results indicated that the performance of FFSCW derived from the evolution of SSCW due to clogging for COD, TOC, total nitrogen (TN), and total phosphor (TP) removal was higher than those of SSCW with the same substrate and plant. The laboratory culture experiments were adopted to evaluate the effect of the constructed wetland evolution on the organic matter mineralization, nitrification/denitrification as well as removal of nitrogen and phosphor. It was shown that, after evolution of SSCW into FFSCW, the mineralization rate for organic matter (as TOC) was 1.82 mg x h(-1), and it was 1.49 mg x h(-1) for SSCW. The removal efficiency for NO3(-) was 96.8%, and it was 58.1% for SSCW. The abiotic denitrification removal efficiency was 40%, and it was 28.2% for SSCW. In addition, the maximum equilibrium adsorption capacity of the substrate after evolution for phosphor (as P) was 160 mg x kg(-1), and it was 140 mg x kg(-1) for SSCW substrate. The organic coverage of the substrate was found to be beneficial to phosphor removal. The nitrification ability decreased after evolution. These results suggest the important effect of constructed wetland evolution on its performance. PMID:23323410
Ober, J.
1998-01-01
The lithium industry can be divided into two sectors: ore concentrate producers and chemical producers. Ore concentrate producers mine lithium minerals. They beneficiate the ores to produce material for use in ceramics and glass manufacturing.
A numerical procedure for three-dimensional transient free surface seepage
NASA Astrophysics Data System (ADS)
Desai, C. S.; Lightner, J. G.; Somasundaram, S.
A three-dimensional procedure based on the finite element method is proposed for transient free surface seepage. It involves solution of the governing equations by using a time integration scheme. The procedure is applied for solution of confined, and transient free surface flow; the latter includes verification with respect to test results from a laboratory model. It is also applied to free surface flow through a dam with a crack.
Particle accumulation on periodic orbits by repeated free surface collisions
NASA Astrophysics Data System (ADS)
Hofmann, Ernst; Kuhlmann, Hendrik C.
2011-07-01
The motion of small particles suspended in cylindrical thermocapillary liquid bridges is investigated numerically in order to explain the experimentally observed particle accumulation structures (PAS) in steady two- and time-dependent three-dimensional flows. Particles moving in this flow are modeled as perfect tracers in the bulk, which can undergo collisions with the free surface. By way of free-surface collisions the particles are transferred among different streamlines which represents the particle trajectories in the bulk. The inter-streamline transfer-process near the free surface together with the passive transport through the bulk is used to construct an iterative map that can describe the accumulation process as an attraction to a stable fixed point which represents PAS. The flow topology of the underlying azimuthally traveling hydrothermal wave turns out to be of key importance for the existence of PAS. In a frame of reference exactly rotating with the hydrothermal wave the three-dimensional flow is steady and exhibits co-existing regular and chaotic streamlines. We find that particles are attracted to accumulation structures if a closed regular streamline exists in the rotating frame of reference which closely approaches the free surface locally. Depending on the closed streamline and the particle radius PAS can arise as a specific trajectory which winds about the closed regular streamline or as the surface of a particular stream tube containing the closed streamline.
Numerical evaluation of gas core length in free surface vortices
NASA Astrophysics Data System (ADS)
Cristofano, L.; Nobili, M.; Caruso, G.
2014-11-01
The formation and evolution of free surface vortices represent an important topic in many hydraulic intakes, since strong whirlpools introduce swirl flow at the intake, and could cause entrainment of floating matters and gas. In particular, gas entrainment phenomena are an important safety issue for Sodium cooled Fast Reactors, because the introduction of gas bubbles within the core causes dangerous reactivity fluctuation. In this paper, a numerical evaluation of the gas core length in free surface vortices is presented, according to two different approaches. In the first one, a prediction method, developed by the Japanese researcher Sakai and his team, has been applied. This method is based on the Burgers vortex model, and it is able to estimate the gas core length of a free surface vortex starting from two parameters calculated with single-phase CFD simulations. The two parameters are the circulation and the downward velocity gradient. The other approach consists in performing a two-phase CFD simulation of a free surface vortex, in order to numerically reproduce the gas- liquid interface deformation. Mapped convergent mesh is used to reduce numerical error and a VOF (Volume Of Fluid) method was selected to track the gas-liquid interface. Two different turbulence models have been tested and analyzed. Experimental measurements of free surface vortices gas core length have been executed, using optical methods, and numerical results have been compared with experimental measurements. The computational domain and the boundary conditions of the CFD simulations were set consistently with the experimental test conditions.
Emission of dimers from a free surface of heated water
NASA Astrophysics Data System (ADS)
Bochkarev, A. A.; Polyakova, V. I.
2014-09-01
The emission rate of water dimers from a free surface and a wetted solid surface in various cases was calculated by a simplified Monte Carlo method with the use of the binding energy of water molecules. The binding energy of water molecules obtained numerically assuming equilibrium between the free surface of water and vapor in the temperature range of 298-438 K corresponds to the coordination number for liquid water equal to 4.956 and is close to the reference value. The calculation results show that as the water temperature increases, the free surface of water and the wetted solid surface become sources of free water dimers. At a temperature of 438 K, the proportion of dimers in the total flow of water molecules on its surface reaches 1%. It is found that in the film boiling mode, the emission rate of dimers decreases with decreasing saturation vapor. Two mechanisms of the emission are described.
Free surfaces in open capillary channels—Parallel plates
NASA Astrophysics Data System (ADS)
Bronowicki, P.; Canfield, P.; Grah, A.; Dreyer, M.
2015-01-01
This paper is concerned with forced flow through partially open capillary channels under microgravity conditions. The investigated channel consists of two parallel plates and is bounded by free liquid surfaces along the open sides. The curvature of the channel's gas-liquid interface, which is exposed to the ambient pressure, adjusts to the pressure difference across the interface in accordance with the Young-Laplace equation. Flow within the channel becomes unstable when the free surface collapses and gas ingestion into the flow path occurs—a process that is also referred to as the "choking" phenomenon. During stable flow, the behavior of the free surface is influenced by flow conditions, geometric properties of the channel, and the pre-defined system pressure. In this work, a previously published stability theory is verified for a wide range of model parameters. A detailed study is provided for stable flow in capillary channels, including static and dynamic solutions. The results of the Capillary Channel Flow (CCF) experiment are evaluated and are found to agree well with numerical predictions. A clear limit is determined between stable and unstable flows. It is shown that the model can predict the shape of the free surface under various flow conditions. A numerical tool is employed to exploit the mathematical model, and the general behavior of free surfaces in said capillary channels is studied. Studies are conducted in both viscous and convective flow regimes and in the transition area between the two. The validity of the model is confirmed for a wide range of geometrical configurations and parameters.
... bipolar disorder (manic-depressive disorder; a disease that causes episodes of depression, episodes of mania, and other abnormal moods). Lithium ... Lithium is also sometimes used to treat depression, schizophrenia (a mental ... emotions), disorders of impulse control (inability to resist the urge ...
Bombardment of thin lithium films with energetic plasma flows
NASA Astrophysics Data System (ADS)
Gray, Travis Kelly
The Divertor Erosion and Vapor Shielding Experiment (DEVEX) has been constructed in the Center for Plasma-Material Interactions at the University of Illinois at Urbana-Champaign. It consists of a conical theta-pinch connected to a 60 kV, 36 muF capacitor bank which is switched with a rise time of 3.5 mus. This results in a peak current of 300 kA for a 30 kV charge on the capacitor bank. The resulting plasma is created and compressed under the theta-pinch coil and then expelled axially towards a target chamber due to the conical taper of the theta-coil. The plasma that reaches the target chamber is dense, 1021 m-3 and cool, 10--20 eV. For the purposes of this study, a thin stainless steel target, sputter coated with a lithium magnetron is the target/material of interest. Both computational [A. Hassanein, Fus. Eng. Des. 60: 527546 (2002)] and experimental [M.L. Apicella, et al., J. Nuc. Mater. 386--388:821823 (2009)] studies have shown that lithium, under fusion relevant plasma bombardment, maintains a much lower surface temperature than other plasma facing materials such as tungsten or carbon. This is believed clue to the strong evaporation and/or sputtering of lithium under these conditions. Subsequently a vapor cloud is formed in front of the plasma-facing surface. The lithium vapor interacts with the incident plasma stream absorbing a fraction of the incident plasma energy via the lithium-plasma interactions. Here, we present experimental verification of substantially reduced target surface temperature with the use of thin lithium films on surface of the target as compared to a bare target. Furthermore, optical measurements are made to determine the density and temperature of the lithium vapor cloud as it expands away from the target surface. A collisional-radiative model for both neutrals and singly ionized lithium is used to model the lithium vapor and is found to correlate well with the optical measurements. The vapor cloud electron temperature is found to
Glassy Dynamics Altered by a Free Surface
NASA Astrophysics Data System (ADS)
Tsui, Ophelia
Studies of polymer dynamics in thin films showed that a highly mobile region exists at the free surface of most if not all polymers. In this talk, I shall review some of these observations, with highlights given to the recent findings that chain flexibility and connectivity may on occasions be necessary for the free surface to exercise its influence. Afterward, I shall ponder on how the influence of the free surface may penetrate as far as several polymer radii of gyration into the inner region, as found both in experiments and simulations. Near the glass transition temperature, our MD simulations showed that the dynamics consist mainly of string-like particle hopping motions, as found by others. Importantly, as the temperature decreases, the hopping motions become increasingly repetitive and back-and-forth, contributing no structural relaxations. We propose that structural relaxations are then brought about by pair-interactions between strings. Near the free surface, however, similar repetitive hopping motions are only observed sufficiently far removed from the free surface. We propose that the free surface induces a penetrating surface mobile region by breaking the memory in the particle dynamics. A possible mechanism based on string interactions will be discussed. We are grateful to the support of NSF through Project DMR-1310536 and Hong Kong GRF Grant 15301014.
A finite difference model for free surface gravity drainage
Couri, F.R.; Ramey, H.J. Jr.
1993-09-01
The unconfined gravity flow of liquid with a free surface into a well is a classical well test problem which has not been well understood by either hydrologists or petroleum engineers. Paradigms have led many authors to treat an incompressible flow as compressible flow to justify the delayed yield behavior of a time-drawdown test. A finite-difference model has been developed to simulate the free surface gravity flow of an unconfined single phase, infinitely large reservoir into a well. The model was verified with experimental results in sandbox models in the literature and with classical methods applied to observation wells in the Groundwater literature. The simulator response was also compared with analytical Theis (1935) and Ramey et al. (1989) approaches for wellbore pressure at late producing times. The seepage face in the sandface and the delayed yield behavior were reproduced by the model considering a small liquid compressibility and incompressible porous medium. The potential buildup (recovery) simulated by the model evidenced a different- phenomenon from the drawdown, contrary to statements found in the Groundwater literature. Graphs of buildup potential vs time, buildup seepage face length vs time, and free surface head and sand bottom head radial profiles evidenced that the liquid refills the desaturating cone as a flat moving surface. The late time pseudo radial behavior was only approached after exaggerated long times.
Heat transfer to a thin liquid film with a free surface
NASA Technical Reports Server (NTRS)
Rahman, M. M.; Faghri, A.; Hankey, W. L.; Swanson, T. D.
1989-01-01
The numerically computed flow field and heat transfer coefficient are presented for the free surface flow of a thin liquid film in the presence or absence of a gravitational body force. The results of an approximate analysis using the Pohlhausen integral method is also discussed. The flow systems studied here include both plane and radial film flows in the presence or absence of a gravitational body force. The heating conditions include isothermal and uniformly heated surfaces. The transport conditions considered at the free surface are an adiabatic condition and an evaporative free surface maintained at its saturation temperature. The height of the free surface, flow field and heat transfer coefficient, were found to be strongly affected by the gravitational body force. They were also found to depend on the Reynolds number and Froude number of the incoming fluid. In the presence of gravity, a hydraulic jump was found to occur under some flow conditions.
Sulphur-impregnated flow cathode to enable high-energy-density lithium flow batteries.
Chen, Hongning; Zou, Qingli; Liang, Zhuojian; Liu, Hao; Li, Quan; Lu, Yi-Chun
2015-01-01
Redox flow batteries are promising technologies for large-scale electricity storage, but have been suffering from low energy density and low volumetric capacity. Here we report a flow cathode that exploits highly concentrated sulphur-impregnated carbon composite, to achieve a catholyte volumetric capacity 294 Ah l(-1) with long cycle life (>100 cycles), high columbic efficiency (>90%, 100 cycles) and high energy efficiency (>80%, 100 cycles). The demonstrated catholyte volumetric capacity is five times higher than the all-vanadium flow batteries (60 Ah l(-1)) and 3-6 times higher than the demonstrated lithium-polysulphide approaches (50-117 Ah l(-1)). Pseudo-in situ impedance and microscopy characterizations reveal superior electrochemical and morphological reversibility of the sulphur redox reactions. Our approach of exploiting sulphur-impregnated carbon composite in the flow cathode creates effective interfaces between the insulating sulphur and conductive carbon-percolating network and offers a promising direction to develop high-energy-density flow batteries. PMID:25565112
Sulphur-impregnated flow cathode to enable high-energy-density lithium flow batteries
NASA Astrophysics Data System (ADS)
Chen, Hongning; Zou, Qingli; Liang, Zhuojian; Liu, Hao; Li, Quan; Lu, Yi-Chun
2015-01-01
Redox flow batteries are promising technologies for large-scale electricity storage, but have been suffering from low energy density and low volumetric capacity. Here we report a flow cathode that exploits highly concentrated sulphur-impregnated carbon composite, to achieve a catholyte volumetric capacity 294 Ah l-1 with long cycle life (>100 cycles), high columbic efficiency (>90%, 100 cycles) and high energy efficiency (>80%, 100 cycles). The demonstrated catholyte volumetric capacity is five times higher than the all-vanadium flow batteries (60 Ah l-1) and 3-6 times higher than the demonstrated lithium-polysulphide approaches (50-117 Ah l-1). Pseudo-in situ impedance and microscopy characterizations reveal superior electrochemical and morphological reversibility of the sulphur redox reactions. Our approach of exploiting sulphur-impregnated carbon composite in the flow cathode creates effective interfaces between the insulating sulphur and conductive carbon-percolating network and offers a promising direction to develop high-energy-density flow batteries.
Flow Boiling Heat Transfer to Lithium Bromide Aqueous Solution in Subcooled Region
NASA Astrophysics Data System (ADS)
Furukawa, Masahiro; Kaji, Masao; Nishizumi, Takeharu; Ozaki, Shinji; Sekoguchi, Kotohiko
To improve the thermal performance of high temperature generator of absorption chiller/heater, heat transfer characteristics of flow boiling of lithium bromide aqueous solution in the subcooled region were experimentally investigated. Experiments were made for water and lithium bromide aqueous solution flowing in a rectangular channel (5 mm × 20 mm cross section) with one side wall heated. Boiling onset quality of lithium bromide aqueous solution is greater than that of water. The heat transfer coefficient of lithium bromide aqueous solution is about a half of that of water under the same experimental conditions of inlet velocity and heat flux. The experimental data of heat transfer coefficient for water are compared with the empirical correlation of Thom et al.11) and a fairly good agreement is obtained. The predictive calculations by the method of Sekoguchi et al.12) are compared with the data for water and lithium bromide aqueous solution. Agreement between them is good for water, while the results for lithium bromide aqueous solution are not satisfactory.
A numerical study of three-dimensional viscous interactions of vortices with a free surface
NASA Astrophysics Data System (ADS)
Dommermuth, D.; Yue, D.
Semi-implicit and explicit numerical methods are developed for the direct simulation of the three-dimensional Navier-Stokes equations with a free surface. The efficiency of a novel multigrid flow solver permits the simulation of three-dimensional flows with free surfaces at low Reynolds and Froude numbers. The numerical schemes are used to study vortex rings and tubes interacting with walls and free surfaces. In the case of vortex rings interacting with a no-slip wall observations of the formation of secondary and tertiary vortex rings agree qualitatively with experimental measurements. When a free surface is present, the results are sensitive to the Froude number. For sufficiently low Froude numbers, the free surface behaves like a free-slip wall, which agrees qualitatively with experimental observations of vortex rings interacting with clean free surfaces. At intermediate Froude numbers, the normal incidence of a vortex ring with a clean free surface results in the formation of secondary vortex rings. Numerical studies of vortex tubes interacting with free surfaces show two possible mechanisms for the reconnection of vorticity with a free surface including primary and secondary vorticity reconnections. One type of primary vorticity reconnection should result in a cusp pattern on the free surface and secondary vorticity reconnections should manifest themselves as paired dimples on the free surface. The essential stages of the reconnection of secondary vorticity with the free surface are as follows: generation of helical vortex sheets by the primary vortex tube, stripping of the helical vortex sheets due to self-induced straining flows, attachment of the helical vortex sheets to the separated free-surface boundary layer, wrapping of U-shaped vortices around the primary vortex tube, feeding of boundary-layer vorticity into the U-vortices, and reconnection of U-vortices with the free surface. We provide evidence which suggests that the striations that may be observed
... mania (frenzied, abnormally excited mood) in people with bipolar disorder (manic-depressive disorder; a disease that causes episodes of depression, episodes of mania, and other abnormal moods). Lithium is in a ... antimanic agents. It works by decreasing abnormal activity in the brain.
Effect of free surface shape on combined thermocapillary and natural convection
NASA Technical Reports Server (NTRS)
Kamotani, Yasuhiro; Platt, Jonathan
1992-01-01
Combined thermocapillary and natural convection in an open square cavity with differentially-heated side walls is studied numerically as well as experimentally. The test fluid is silicone oil with Prandtl number of 105. The shape of fluid-free surface is made either flat or curved to study its effect on the flow. A finite difference scheme to deal with a curved free surface is developed. The experimental results shown agree with the numerical results. With the curved-free surface, the flow and local heat transfer rate are reduced in the corner regions, and a sharp peak in heat transfer rate at the top edge of the cold wall disappears.
Can a dike "feel" a free surface?
NASA Astrophysics Data System (ADS)
Pansino, S.; Taisne, B.
2015-12-01
This study uses analogue experiments to examine the various factors that could influence the dynamic or arrest of a dike prior to an eruption. Dike propagation has been well-studied from the perspective of analogue experiments and it has been shown how various competing forces interact (buoyancy, viscosity, elastic properties, fracture toughness,…). It has been observed that dikes experience acceleration as they approach a free surface, which theoretically is due to increasing stress intensity as distance to the free surface decreases. Other unrelated studies analysed the constant-volume propagation of low-surface-energy liquids (e.g. hydrophilic liquids) in a semi-infinite medium and found the propagation rate to decelerate, as liquid gradually transferred from the head region to the tail and buoyancy abated. In contrast, high-surface-energy fluids transfer negligible quantities of their volume into the tail region, essentially maintain their volume, and experience no such deceleration. This study focuses on the aforementioned factors influencing dike propagation, acceleration due to a free surface and deceleration due to head volume loss, to determine what their combined effect is on magma eruptability. As previous studies used hydrophobic liquids to investigate free-surface effects, and as magma is expected to have a low surface energy with the crust (they are both silicate), hydrophilic liquids (e.g. glucose syrup) are used in this study. Each experiment is arranged to have a free surface condition occupying half of the apparatus, and a no-slip condition occupying the other half, which allows for side-by-side comparison. The results are analysed to determine the minimum volume needed for a dike to erupt; failure to erupt implies that the dike would instead come to arrest in the crust.
A liquid metal flume for free surface magnetohydrodynamic experiments
NASA Astrophysics Data System (ADS)
Nornberg, M. D.; Ji, H.; Peterson, J. L.; Rhoads, J. R.
2008-09-01
We present an experiment designed to study magnetohydrodynamic effects in free surface channel flow. The wide aspect ratio channel (the width to height ratio is about 15) is completely enclosed in an inert atmosphere to prevent oxidization of the liquid metal. A custom-designed pump reduces entrainment of oxygen, which was found to be a problem with standard centrifugal and gear pumps. Laser Doppler velocimetry experiments characterize velocity profiles of the flow. Various flow constraints mitigate secondary circulation and end effects on the flow. Measurements of the wave propagation characteristics in the liquid metal demonstrate the surfactant effect of surface oxides and the damping of fluctuations by a cross-channel magnetic field.
A Liquid Metal Flume for Free Surface Magnetohydrodynamic Experiments
Nornberg, M.D.; Ji, H.; Peterson, J.L.; Rhoads, J.R.
2008-08-27
We present an experiment designed to study magnetohydrodynamic effects in free-surface channel flow. The wide aspect ratio channel (the width to height ratio is about 15) is completely enclosed in an inert atmosphere to prevent oxidization of the liquid metal. A custom-designed pump reduces entrainment of oxygen, which was found to be a problem with standard centrifugal and gear pumps. Laser Doppler Velocimetry experiments characterize velocity profiles of the flow. Various flow constraints mitigate secondary circulation and end effects on the flow. Measurements of the wave propagation characteristics in the liquid metal demonstrate the surfactant effect of surface oxides and the damping of fluctuations by a cross-channel magnetic field.
NASA Astrophysics Data System (ADS)
Longo, Sandro
2011-01-01
This paper is a companion paper to a study devoted to the analysis of experimental instantaneous fluid levels and three-component fluid velocity measurements in a stationary flow field generated by a Crump weir in a laboratory flume, using an ultrasonic distance sensor and a three-probe arrangement of an ultrasonic Doppler velocity profiler (UVP) (Longo in Exp Fluids, doi: 10.1007/s00348-010-0881-5, 2010). Whereas Longo (Exp Fluids, doi: 10.1007/s00348-010-0881-5, 2010) deals with a general overview of the problem, the description of the experiments and the analysis of the free surface statistics and relevant scales, the present manuscript is devoted to a detailed analysis of the turbulence and the correlation with the free surface. The data are elaborated by obtaining the macroturbulence Reynolds tensor, using conditional averages based on free surface-fluctuation statistics. We also compute the two-point correlations of fluid velocity, the micro-scale and the integral scale, the correlation between free surface and the turbulence beneath. A free surface-boundary layer was detected having a thickness proportional to the root mean square of the free surface-height series and with a velocity scale that correlates well with the free surface-elevation time gradient. Most of the relevant state variables, such as the mean velocity and Reynolds stress components, collapse to a single curve if scaled appropriately. There are many indicators that a specific regime occurs that has an optimal tuning between the free surface and turbulence. In this regime, the length scales are considered as an indicator.
A flowing liquid lithium limiter for the Experimental Advanced Superconducting Tokamak
Ren, J.; Zuo, G. Z.; Hu, J. S.; Sun, Z.; Yang, Q. X.; Li, J. G.; Xie, H.; Chen, Z. X.; Zakharov, L. E.
2015-02-15
A program involving the extensive and systematic use of lithium (Li) as a “first,” or plasma-facing, surface in Tokamak fusion research devices located at Institute of Plasma Physics, Chinese Academy of Sciences, was started in 2009. Many remarkable results have been obtained by the application of Li coatings in Experimental Advanced Superconducting Tokamak (EAST) and liquid Li limiters in the HT-7 Tokamak—both located at the institute. In furtherance of the lithium program, a flowing liquid lithium (FLiLi) limiter system has been designed and manufactured for EAST. The design of the FLiLi limiter is based on the concept of a thin flowing film which was previously tested in HT-7. Exploiting the capabilities of the existing material and plasma evaluation system on EAST, the limiter will be pre-wetted with Li and mechanically translated to the edge of EAST during plasma discharges. The limiter will employ a novel electro-magnetic pump which is designed to drive liquid Li flow from a collector at the bottom of limiter into a distributor at its top, and thus supply a continuously flowing liquid Li film to the wetted plasma-facing surface. This paper focuses on the major design elements of the FLiLi limiter. In addition, a simulation of incoming heat flux has shown that the distribution of heat flux on the limiter surface is acceptable for a future test of power extraction on EAST.
A flowing liquid lithium limiter for the Experimental Advanced Superconducting Tokamak.
Ren, J; Zuo, G Z; Hu, J S; Sun, Z; Yang, Q X; Li, J G; Zakharov, L E; Xie, H; Chen, Z X
2015-02-01
A program involving the extensive and systematic use of lithium (Li) as a "first," or plasma-facing, surface in Tokamak fusion research devices located at Institute of Plasma Physics, Chinese Academy of Sciences, was started in 2009. Many remarkable results have been obtained by the application of Li coatings in Experimental Advanced Superconducting Tokamak (EAST) and liquid Li limiters in the HT-7 Tokamak-both located at the institute. In furtherance of the lithium program, a flowing liquid lithium (FLiLi) limiter system has been designed and manufactured for EAST. The design of the FLiLi limiter is based on the concept of a thin flowing film which was previously tested in HT-7. Exploiting the capabilities of the existing material and plasma evaluation system on EAST, the limiter will be pre-wetted with Li and mechanically translated to the edge of EAST during plasma discharges. The limiter will employ a novel electro-magnetic pump which is designed to drive liquid Li flow from a collector at the bottom of limiter into a distributor at its top, and thus supply a continuously flowing liquid Li film to the wetted plasma-facing surface. This paper focuses on the major design elements of the FLiLi limiter. In addition, a simulation of incoming heat flux has shown that the distribution of heat flux on the limiter surface is acceptable for a future test of power extraction on EAST. PMID:25725839
A flowing liquid lithium limiter for the Experimental Advanced Superconducting Tokamak
NASA Astrophysics Data System (ADS)
Ren, J.; Zuo, G. Z.; Hu, J. S.; Sun, Z.; Yang, Q. X.; Li, J. G.; Zakharov, L. E.; Xie, H.; Chen, Z. X.
2015-02-01
A program involving the extensive and systematic use of lithium (Li) as a "first," or plasma-facing, surface in Tokamak fusion research devices located at Institute of Plasma Physics, Chinese Academy of Sciences, was started in 2009. Many remarkable results have been obtained by the application of Li coatings in Experimental Advanced Superconducting Tokamak (EAST) and liquid Li limiters in the HT-7 Tokamak—both located at the institute. In furtherance of the lithium program, a flowing liquid lithium (FLiLi) limiter system has been designed and manufactured for EAST. The design of the FLiLi limiter is based on the concept of a thin flowing film which was previously tested in HT-7. Exploiting the capabilities of the existing material and plasma evaluation system on EAST, the limiter will be pre-wetted with Li and mechanically translated to the edge of EAST during plasma discharges. The limiter will employ a novel electro-magnetic pump which is designed to drive liquid Li flow from a collector at the bottom of limiter into a distributor at its top, and thus supply a continuously flowing liquid Li film to the wetted plasma-facing surface. This paper focuses on the major design elements of the FLiLi limiter. In addition, a simulation of incoming heat flux has shown that the distribution of heat flux on the limiter surface is acceptable for a future test of power extraction on EAST.
Free surface deformation and heat transfer by thermocapillary convection
NASA Astrophysics Data System (ADS)
Fuhrmann, Eckart; Dreyer, Michael; Basting, Steffen; Bänsch, Eberhard
2016-04-01
Knowing the location of the free liquid/gas surface and the heat transfer from the wall towards the fluid is of paramount importance in the design and the optimization of cryogenic upper stage tanks for launchers with ballistic phases, where residual accelerations are smaller by up to four orders of magnitude compared to the gravity acceleration on earth. This changes the driving forces drastically: free surfaces become capillary dominated and natural or free convection is replaced by thermocapillary convection if a non-condensable gas is present. In this paper we report on a sounding rocket experiment that provided data of a liquid free surface with a nonisothermal boundary condition, i.e. a preheated test cell was filled with a cold but storable liquid in low gravity. The corresponding thermocapillary convection (driven by the temperature dependence of the surface tension) created a velocity field directed away from the hot wall towards the colder liquid and then in turn back at the bottom towards the wall. A deformation of the free surface resulting in an apparent contact angle rather different from the microscopic one could be observed. The thermocapillary flow convected the heat from the wall to the liquid and increased the heat transfer compared to pure conduction significantly. The paper presents results of the apparent contact angle as a function of the dimensionless numbers (Weber-Marangoni and Reynolds-Marangoni number) as well as heat transfer data in the form of a Nusselt number. Experimental results are complemented by corresponding numerical simulations with the commercial software Flow3D and the inhouse code Navier.
Thermocapillary convection in two immiscible liquid layers with free surface
NASA Technical Reports Server (NTRS)
Doi, Takao; Koster, Jean N.
1993-01-01
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.
Flow Boiling Heat Transfer to Lithium Bromide Aqueous Solution in Subcooled Region
NASA Astrophysics Data System (ADS)
Kaji, Masao; Furukawa, Masahiro; Nishizumi, Takeharu; Ozaki, Shinji; Sekoguchi, Kotohiko
A theoretical prediction model of the boiling heat transfer coefficient in the subcooled region for water and lithium bromide aqueous solution flowing in a rectangular channel is proposed. In the present heat transfer model, a heat flux is assumed to consist of both the forced convective and the boiling effect components. The forced convective component is evaluated from the empirical correlation of convective heat transfer coefficient for single-phase flow considering the effect of increase of liquid velocity due to net vapor generation. Empirical correlations for determining the heat flux due to the boiling effect and the quality at the onset point of net vapor generation are obtained from the data presented in the first report1). Agreement between the present theoretical prediction and the experimental data is satisfactorily good both for water and lithium bromide aqueous solution.
NASA Astrophysics Data System (ADS)
Hu, J. S.; Zuo, G. Z.; Ren, J.; Yang, Q. X.; Chen, Z. X.; Xu, H.; Zakharov, L. E.; Maingi, R.; Gentile, C.; Meng, X. C.; Sun, Z.; Xu, W.; Chen, Y.; Fan, D.; Yan, N.; Duan, Y. M.; Yang, Z. D.; Zhao, H. L.; Song, Y. T.; Zhang, X. D.; Wan, B. N.; Li, J. G.; EAST Team
2016-04-01
As an alternative choice of solid plasma facing components (PFCs), flowing liquid lithium can serve as a limiter or divertor PFC and offers a self-healing surface with acceptable heat removal and good impurity control. Such a system could improve plasma performance, and therefore be attractive for future fusion devices. Recently, a continuously flowing liquid lithium (FLiLi) limiter has been successfully designed and tested in the EAST superconducting tokamak. A circulating lithium layer with a thickness of <0.1 mm and a flow rate ~2 cm3 s-1 was achieved. A novel in-vessel electro-magnetic pump, working with the toroidal magnetic field of the EAST device, was reliable to control the lithium flow speed. The flowing liquid limiter was found to be fully compatible with various plasma scenarios, including high confinement mode plasmas heated by lower hybrid waves or by neutral beam injection. It was also found that the controllable lithium emission from the limiter was beneficial for the reduction of recycling and impurities, for the reduction of divertor heat flux, and in certain cases, for the improvement of plasma stored energy, which bodes well application for the use of flowing liquid lithium PFCs in future fusion devices.
A carbon-free lithium-ion solid dispersion redox couple with low viscosity for redox flow batteries
NASA Astrophysics Data System (ADS)
Qi, Zhaoxiang; Koenig, Gary M.
2016-08-01
A new type of non-aqueous redox couple without carbon additives for flow batteries is proposed and the target anolyte chemistry is demonstrated. The so-called "Solid Dispersion Redox Couple" incorporates solid electroactive materials dispersed in organic lithium-ion battery electrolyte as its flowing suspension. In this work, a unique and systematic characterization approach has been used to study the flow battery redox couple in half cell demonstrations relative to a lithium electrode. An electrolyte laden with Li4Ti5O12 (LTO) has been characterized in multiple specially designed lithium half cell configurations. The flow battery redox couple described in this report has relatively low viscosity, especially in comparison to other flow batteries with solid active materials. The lack of carbon additive allows characterization of the electrochemical properties of the electroactive material in flow without the complication of conductive additives and unambiguous observation of the electrorheological coupling in these dispersed particle systems.
Thermal analysis and two-directional air flow thermal management for lithium-ion battery pack
NASA Astrophysics Data System (ADS)
Yu, Kuahai; Yang, Xi; Cheng, Yongzhou; Li, Changhao
2014-12-01
Thermal management is a routine but crucial strategy to ensure thermal stability and long-term durability of the lithium-ion batteries. An air-flow-integrated thermal management system is designed in the present study to dissipate heat generation and uniformize the distribution of temperature in the lithium-ion batteries. The system contains of two types of air ducts with independent intake channels and fans. One is to cool the batteries through the regular channel, and the other minimizes the heat accumulations in the middle pack of batteries through jet cooling. A three-dimensional anisotropic heat transfer model is developed to describe the thermal behavior of the lithium-ion batteries with the integration of heat generation theory, and validated through both simulations and experiments. Moreover, the simulations and experiments show that the maximum temperature can be decreased to 33.1 °C through the new thermal management system in comparison with 42.3 °C through the traditional ones, and temperature uniformity of the lithium-ion battery packs is enhanced, significantly.
High-energy density nonaqueous all redox flow lithium battery enabled with a polymeric membrane.
Jia, Chuankun; Pan, Feng; Zhu, Yun Guang; Huang, Qizhao; Lu, Li; Wang, Qing
2015-11-01
Redox flow batteries (RFBs) are considered one of the most promising large-scale energy storage technologies. However, conventional RFBs suffer from low energy density due to the low solubility of the active materials in electrolyte. On the basis of the redox targeting reactions of battery materials, the redox flow lithium battery (RFLB) demonstrated in this report presents a disruptive approach to drastically enhancing the energy density of flow batteries. With LiFePO4 and TiO2 as the cathodic and anodic Li storage materials, respectively, the tank energy density of RFLB could reach ~500 watt-hours per liter (50% porosity), which is 10 times higher than that of a vanadium redox flow battery. The cell exhibits good electrochemical performance under a prolonged cycling test. Our prototype RFLB full cell paves the way toward the development of a new generation of flow batteries for large-scale energy storage. PMID:26702440
High–energy density nonaqueous all redox flow lithium battery enabled with a polymeric membrane
Jia, Chuankun; Pan, Feng; Zhu, Yun Guang; Huang, Qizhao; Lu, Li; Wang, Qing
2015-01-01
Redox flow batteries (RFBs) are considered one of the most promising large-scale energy storage technologies. However, conventional RFBs suffer from low energy density due to the low solubility of the active materials in electrolyte. On the basis of the redox targeting reactions of battery materials, the redox flow lithium battery (RFLB) demonstrated in this report presents a disruptive approach to drastically enhancing the energy density of flow batteries. With LiFePO4 and TiO2 as the cathodic and anodic Li storage materials, respectively, the tank energy density of RFLB could reach ~500 watt-hours per liter (50% porosity), which is 10 times higher than that of a vanadium redox flow battery. The cell exhibits good electrochemical performance under a prolonged cycling test. Our prototype RFLB full cell paves the way toward the development of a new generation of flow batteries for large-scale energy storage. PMID:26702440
Material and Energy Flows in the Production of Cathode and Anode Materials for Lithium Ion Batteries
Dunn, Jennifer B.; James, Christine; Gaines, Linda; Gallagher, Kevin; Dai, Qiang; Kelly, Jarod C.
2015-09-01
The Greenhouse gases, Regulated Emissions and Energy use in Transportation (GREET) model has been expanded to include four new cathode materials that can be used in the analysis of battery-powered vehicles: lithium nickel cobalt manganese oxide (LiNi_{0.4}Co_{0.2}Mn_{0.4}O_{2} [NMC]), lithium iron phosphate (LiFePO_{4} [LFP]), lithium cobalt oxide (LiCoO_{2} [LCO]), and an advanced lithium cathode (0.5Li_{2}MnO_{3}∙0.5LiNi_{0.44}Co_{0.25}Mn_{0.31}O_{2} [LMR-NMC]). In GREET, these cathode materials are incorporated into batteries with graphite anodes. In the case of the LMR-NMC cathode, the anode is either graphite or a graphite-silicon blend. Lithium metal is also an emerging anode material. This report documents the material and energy flows of producing each of these cathode and anode materials from raw material extraction through the preparation stage. For some cathode materials, we considered solid state and hydrothermal preparation methods. Further, we used Argonne National Laboratory’s Battery Performance and Cost (BatPaC) model to determine battery composition (e.g., masses of cathode, anode, electrolyte, housing materials) when different cathode materials were used in the battery. Our analysis concluded that cobalt- and nickel-containing compounds are the most energy intensive to produce.
Free-surface interactions in the wake of an inverted cylindrical pendulum
NASA Astrophysics Data System (ADS)
Voorhees, A. V.; Wei, T.
2001-11-01
Vortex-induced vibration experiments of a low-mass circular cylinder were conducted in the Rutgers Free-Surface Water Tunnel. The focus of this study was to ascertain the extent of free-surface effects on the cylinder motion. A 2.54-cm diameter, 160 cm long circular cylinder was mounted vertically in the test section as an inverted pendulum. The bottom was attached to the tunnel floor by a leaf spring and the upper end protruded through the free surface. Strong vertical flows, associated with the Kármán vortex street, have been observed in the wake; such flows have not been found in transverse oscillator wakes that exhibit similar response data. In addition, these vertical flows are markedly altered near the free surface. The 3-D wake and free-surface interactions were examined in the synchronization regime at Re = 3800. In this regime, the amplitude response exhibits highly modulated, quasi-periodic beating. Temporally and spatially-resolved, phase-averaged DPIV measurements were made in horizontal and cross-stream planes. The oblique nature, vertical and streamwise, of the vortex shedding will be shown using composite time plots from both measurement planes. Changes in vortex shedding phase and formation length will be related to response amplitude modulations.
Free Surface and Flapping Foil Interactions
NASA Astrophysics Data System (ADS)
Ananthakrishnan, Palaniswamy
2014-11-01
Flapping foils for station-keeping of a near-surface body in a current is analyzed using a finite-difference method based on boundary-fitted coordinates. The foils are hinge-connected to the aft of the body and subject to pitch oscillation. Results are obtained for a range of Strouhal number, Froude number, unsteady frequency parameter τ, Reynolds number and the depth of foil submergence. Results show that at low Strouhal number (St < 0 . 1) and sub-critical unsteady parameter τ < 0 . 25 , the flapping generates drag instead of thrust. At high Strouhal number and super-critical value of the unsteady parameter (τ > 0 . 25) flapping generates high thrust with low efficiency. Thrust and efficiency are found to decrease with decreasing submergence depth of the foil. At the critical τ = 0 . 25 and shallow submergence of the foil, the standing wave generated above the foil continues to grow until breaking; both the thrust and efficiency of the foil are reduced at the critical τ. The necessary conditions for optimal thrust generation by a flapping foil underneath the free surface are found to be (i) Strouhal number in the range from 0.25 to 0.35, (ii) unsteady parameter τ > 0 . 25 and (iii) the maximum angle of attack less than 15° for the flat-plate foil. Supported by the US Office of Naval Research through the Naval Engineering Education Center (NEEC) Consortium of the University of Michigan, Ann Arbor.
46 CFR 42.20-10 - Free surface.
Code of Federal Regulations, 2014 CFR
2014-10-01
... 46 Shipping 2 2014-10-01 2014-10-01 false Free surface. 42.20-10 Section 42.20-10 Shipping COAST... Freeboards § 42.20-10 Free surface. When doing the calculations required in §§ 42.20-6(a), 42.20-7(a) and 42.20-8(a), the effect of free surface of the following liquids must be included: (a) For each type...
46 CFR 42.20-10 - Free surface.
Code of Federal Regulations, 2011 CFR
2011-10-01
... 46 Shipping 2 2011-10-01 2011-10-01 false Free surface. 42.20-10 Section 42.20-10 Shipping COAST... Freeboards § 42.20-10 Free surface. When doing the calculations required in §§ 42.20-6(a), 42.20-7(a) and 42.20-8(a), the effect of free surface of the following liquids must be included: (a) For each type...
46 CFR 42.20-10 - Free surface.
Code of Federal Regulations, 2010 CFR
2010-10-01
... 46 Shipping 2 2010-10-01 2010-10-01 false Free surface. 42.20-10 Section 42.20-10 Shipping COAST... Freeboards § 42.20-10 Free surface. When doing the calculations required in §§ 42.20-6(a), 42.20-7(a) and 42.20-8(a), the effect of free surface of the following liquids must be included: (a) For each type...
46 CFR 42.20-10 - Free surface.
Code of Federal Regulations, 2013 CFR
2013-10-01
... 46 Shipping 2 2013-10-01 2013-10-01 false Free surface. 42.20-10 Section 42.20-10 Shipping COAST... Freeboards § 42.20-10 Free surface. When doing the calculations required in §§ 42.20-6(a), 42.20-7(a) and 42.20-8(a), the effect of free surface of the following liquids must be included: (a) For each type...
46 CFR 42.20-10 - Free surface.
Code of Federal Regulations, 2012 CFR
2012-10-01
... 46 Shipping 2 2012-10-01 2012-10-01 false Free surface. 42.20-10 Section 42.20-10 Shipping COAST... Freeboards § 42.20-10 Free surface. When doing the calculations required in §§ 42.20-6(a), 42.20-7(a) and 42.20-8(a), the effect of free surface of the following liquids must be included: (a) For each type...
Development of liquid-lithium film jet-flow for the target of (7)Li(p,n)(7)Be reactions for BNCT.
Kobayashi, Tooru; Miura, Kuniaki; Hayashizaki, Noriyosu; Aritomi, Masanori
2014-06-01
A feasibility study on liquid lithium target in the form of a flowing film was performed to evaluate its potential use as a neutron generation target of (7)Li(p,n)(7)Be reaction in BNCT. The target is a windowless-type flowing film on a concave wall. Its configuration was adapted for a proton beam which is 30mm in diameter and with energy and current of up to 3MeV and 20mA, respectively. The flowing film of liquid lithium was 0.6mm in thickness, 50mm in width and 50mm in length. The shapes of the nozzle and concave back wall, which create a stable flowing film jet, were decided based on water experiments. A lithium hydrodynamic experiment was performed to observe the stability of liquid lithium flow behavior. The flowing film of liquid lithium was found to be feasible at temperatures below the liquid lithium boiling saturation of 342°C at the surface pressure of 1×10(-3)Pa. Using a proto-type liquid lithium-circulating loop for BNCT, the stability of the film flow was confirmed for velocities up to 30m/s at 220°C and 250°C in vacuum at a pressure lower than 10(-3) Pa. It is expected that for practical use, a flowing liquid lithium target of a windowless type can solve the problem of radiation damage and target cooling. PMID:24412425
On the interaction of a vertical shear layer with a free surface
NASA Astrophysics Data System (ADS)
Dabiri, Dana
2003-04-01
New experiments have been conducted using a combined free-surface gradient detector (FSGD) and digital particle image velocimetry (DPIV) technique to study the interaction between a vertical shear layer, created by a surface-piercing splitter plate, and a free surface. The emphasis of this study is on understanding aspects of the interaction between the free-surface deformation (FSD) and the near-surface turbulence through the correlations between the elevation and the vorticity fields, and the spectral behaviour of the near-surface pressure. The Reynolds number of the present study, based on visual thickness and the velocity average of the two streams, is 12 100. Mean results for the velocity and vorticity fields show that self-similarity is achieved. Instantaneous data sets show that at the free surface, vortex tubes within the main rollers connect normally with the free surface as is evidenced by strong vorticity as well as the strong deformations at the free surface. The instantaneous data sets also show that the streamwise vortices near the braid regions, while weaker than those seen in the main rollers, also reconnect with the free surface. Statistical analyses show that the FSD is strongly correlated with the near-surface vorticity field, as the correlation coefficients are quite high (˜ 0.7-0.8). The pressure spectrum slope within the shear layer near the surface is found to be -10/3, compared with the analytically derived value of -11/3 (George et al. 1984) found for a shear layer in unbounded flow.
Self-induced oscillation of free surface coupled with the axisymmetric jet and structure
Someya, Satoshi; Okamoto, Koji; Madarame, Haruki
1996-08-01
New self-induced free-surface oscillations were discovered, which were caused by the interaction between a jet, a free surface and a structure. A submerged upward round jet was injected into a cylindrical tank from an inlet nozzle at the bottom center. The jet impinged on a cylindrical rod (UIS) which was set just beneath the free surface of the tank. Then, the jet turned its direction along UIS and reached the free surface making a swell around UIS. Under a certain condition of jet velocity and UIS depth, several kinds of self-induced free surface oscillation were observed, which were roughly classified into two modes with their surface shapes. One had no diametric node and two nodal circles, (0,2), and the other had one node in each direction, (1,1). The (0,2) mode oscillation was divided into two types with their flow pattern. In the (0,2) mode with outward surface velocity, (0,2){sub OUT}, its oscillating frequency decreased with increasing jet velocity. The swell of the free surface characterized (0,2){sub OUT} mode oscillation which was different from the theoretical (0,2) mode sloshing. In the (0,2) mode with inward surface velocity, (0,2){sub IN}, its frequency was independent of jet velocity. Its growth mechanism was thought to be similar to that of self-induced sloshing in the previous study. The (1,1) mode oscillation was also divided into two types with the appearance of the swell around UIS. The (1,1) mode with swelling of the free surface, (1,1){sub {alpha}}, was found to be the same as Jet-Flutter in a cylindrical tank without UIS. This is important to evaluate performance of LMFBR.
Pan, Huilin; Wei, Xiaoliang; Henderson, Wesley A.; Shao, Yuyan; Chen, Junzheng; Bhattacharya, Priyanka; Xiao, Jie; Liu, Jun
2015-04-27
Lithium sulfur (Li-S) redox flow battery (RFB) is a promising candidate for high energy large-scale energy storage application due to good solubility of long-chain polysulfide species and low cost of sulfur. In this report, recent progress and new concepts for Li-S redox flow batteries are discussed with an emphasis on the fundamental understanding and control of lithium polysulfide chemistry to enable the development of liquid phase Li-S redox flow prototype cells. These differ significantly from conventional static Li-S batteries targeting for vehicle electrification. A high solubility of the different lithium polysulfides generated at different depths of discharge and states of charge is required for a flow battery in order to take full advantage of the multiple electron transitions between elemental sulfur and Li2S. A new DMSO-based electrolyte is proposed for Li-S redox flow batteries, which not only enables the high solubility of lithium polysulfide species, especially for the short-chain species, but also results in excellent cycling with a high Coulombic efficiency. The challenges and opportunities for the Li-S redox flow concept have also been discussed in depth.
Hydrodynamic coefficients of an oscillating ellipsoid moving in the free surface
NASA Astrophysics Data System (ADS)
Inglis, R. B.; Price, W. G.
1980-10-01
The frequency dependent heave, pitch, sway, and yaw hydrodynamic coefficients associated with an oscillating ellipsoid traveling with forward speed in the free surface are evaluated from a three-dimensional potential flow analysis. The free-surface boundary condition in the mathematical model either includes the influence of forward speed or is simplified to the equivalent zero speed case. This variation produces a velocity potential which is either frequency and speed dependent or just frequency dependent. The influence of forward speed on all the hydrodynamic coefficients is discussed.
Computational Fluid Dynamics Modeling of a Lithium/Thionyl Chloride Battery with Electrolyte Flow
Gu, W.B.; Jungst, Rudolph G.; Nagasubramanian, Ganesan; Wang, C.Y.; Weidner, John.
1999-06-11
A two-dimensional model is developed to simulate discharge of a lithium/thionyl chloride primary battery. The model accounts for not only transport of species and charge, but also the electrode porosity variations and the electrolyte flow induced by the volume reduction caused by electrochemical reactions. Numerical simulations are performed using a finite volume method of computational fluid dynamics. The predicted discharge curves for various temperatures are compared to the experimental data with excellent agreement. Moreover, the simulation results. in conjunction with computer visualization and animation techniques, confirm that cell utilization in the temperature and current range of interest is limited by pore plugging or clogging of the front side of the cathode as a result of LiCl precipitation. The detailed two-dimensional flow simulation also shows that the electrolyte is replenished from the cell header predominantly through the separator into the front of the cathode during most parts of the discharge, especially for higher cell temperatures.
A numerical procedure for transient free surface seepage through fracture networks
NASA Astrophysics Data System (ADS)
Jiang, Qinghui; Ye, Zuyang; Zhou, Chuangbing
2014-11-01
A parabolic variational inequality (PVI) formulation is presented for the transient free surface seepage problem defined for a whole fracture network. Because the seepage faces are specified as Signorini-type conditions, the PVI formulation can effectively eliminate the singularity of spillpoints that evolve with time. By introducing a continuous penalty function to replace the original Heaviside function, a finite element procedure based on the PVI formulation is developed to predict the transient free surface response in the fracture network. The effects of the penalty parameter on the solution precision are analyzed. A relative error formula for evaluating the flow losses at steady state caused by the penalty parameter is obtained. To validate the proposed method, three typical examples are solved. The solutions for the first example are compared with the experimental results. The results from the last two examples further demonstrate that the orientation, extent and density of fractures significantly affect the free surface seepage behavior in the fracture network.
Apple Snail: a Bio Cleaner of the Water Free Surface.
NASA Astrophysics Data System (ADS)
Bassiri, Golnaz
2005-11-01
Oil spills from tankers represent a threat for shorelines and marine life. Despite continuing research, there has been little change in the fundamental technology for dealing with oil spills. An experimental investigation of the feeding strategy of Apple snails from the water free surface, called surface film feeding, is being studied motivated by the need to develop new techniques to recover oil spills. To feed on floating food (usually a thin layer of microorganisms), the apple snail forms a funnel with its foot and pulls the free surface toward the funnel. High speed imaging and particle image velocimetry were used in the present investigation to measure the free surface motion and to investigate the mechanism used by the apple snails to pull the free surface. The results suggest that the snail pulls the free surface via the wavy motion of the muscles in its funnel.
Numerical Simulation of Floating Bodies in Extreme Free Surface Waves
NASA Astrophysics Data System (ADS)
Hu, Zheng Zheng; Causon, Derek; Mingham, Clive; Qiang, Ling
2010-05-01
and efficient. Firstly, extreme design wave conditions are generated in an empty NWT and compared with physical experiments as a precursor to calculations to investigate the survivability of the Bobber device operating in a challenging wave climate. Secondly, we consider a bench-mark test case involving in a first order regular wave maker acting on a fixed cylinder and Pelamis. Finally, a floating Bobber has been simulated under extreme wave conditions. These results will be reported at the meeting. Causon D.M., Ingram D.M., Mingham C.G., Yang G. Pearson R.V. (2000). Calculation of shallow water flows using a Cartesian cut cell approach. Advances in Water resources, 23: 545-562. Causon D.M., Ingram D.M., Mingham C.G. (2000). A Cartesian cut cell method for shallow water flows with moving boundaries. Advances in Water resources, 24: 899-911. Dalzell J.F. 1999 A note on finite depth second-order wave-wave interactions. Appl. Ocean Res. 21, 105-111. Ning D.Z., Zang J., Liu S.X. Eatock Taylor R. Teng B. & Taylor P.H. 2009 Free surface and wave kinematics for nonlinear focused wave groups. J. Ocean Engineering. Accepted. Hu Z.Z., Causon D.M., Mingham C.M. and Qian L.(2009). Numerical wave tank study of a wave energy converter in heave. Proceedlings 19th ISOPE conference, Osaka, Japan Qian L., Causon D.M. & Mingham C.G., Ingram D.M. 2006 A free-surface capturing method for two fluid flows with moving bodies. Proc. Roy. Soc. London, Vol. A 462 21-42.
Low-Reynolds-number rising of a bubble near a free surface at vanishing Bond number
NASA Astrophysics Data System (ADS)
Guémas, Marine; Sellier, Antoine; Pigeonneau, Franck
2016-06-01
This work considers a nearly spherical bubble and a nearly flat free surface interacting under buoyancy at vanishing Bond number Bo. For each perturbed surface, the deviation from the unperturbed shape is asymptotically obtained at leading order on Bo. The task appeals to the normal traction exerted on the unperturbed surface by the Stokes flow due to a spherical bubble translating toward a flat free surface. The free surface problem is then found to be well-posed and to admit a solution in closed form when gravity is still present in the linear differential equation governing the perturbed profile through a term proportional to Bo. In contrast, the bubble problem amazingly turns out to be over-determined. It however becomes well-posed if the requirement of horizontal tangent planes at the perturbed bubble north and south poles is discarded or if the term proportional to Bo is omitted. Both previous approaches turn out to predict for a small Bond number, quite close solutions except in the very vicinity of the bubble poles. The numerical solution of the proposed asymptotic analysis shows in the overlapping range Bo = O ( 0.1 ) and for both the bubble and the free surface perturbed shapes, a good agreement with a quite different boundary element approach developed in Pigeonneau and Sellier ["Low-Reynolds-number gravity-driven migration and deformation of bubbles near a free surface," Phys. Fluids 23, 092102 (2011)]. It also provides approximated bubble and free surface shapes whose sensitivity to the bubble location is examined.
Dynamics of free surface perturbations along an annular viscous film
NASA Astrophysics Data System (ADS)
Smolka, Linda B.; North, Justin; Guerra, Bree K.
2008-03-01
It is known that the free surface of an axisymmetric viscous film flowing down the outside of a thin vertical fiber under the influence of gravity becomes unstable to interfacial perturbations. We present an experimental study using fluids with different densities, surface tensions, and viscosities to investigate the growth and dynamics of these interfacial perturbations and to test the assumptions made by previous authors. We find that the initial perturbation growth is exponential, followed by a slower phase as the amplitude and wavelength saturate in size. Measurements of the perturbation growth for experiments conducted at low and moderate Reynolds numbers are compared to theoretical predictions developed from linear stability theory. Excellent agreement is found between predictions from a long-wave Stokes flow model [Craster and Matar, J. Fluid Mech. 553, 85 (2006)] and data, while fair to excellent agreement (depending on fiber size) is found between predictions from a moderate-Reynolds-number model [Sisoev , Chem. Eng. Sci. 61, 7279 (2006)] and data. Furthermore, we find that a known transition in the longer-time perturbation dynamics from unsteady to steady behavior at a critical flow rate Qc is correlated with a transition in the rate at which perturbations naturally form along the fiber. For Q
NASA Astrophysics Data System (ADS)
Lovecchio, Salvatore; Zonta, Francesco; Soldati, Alfredo
2014-10-01
The dispersion of floaters, small organic particles lighter than water, on the free surface of an open turbulent channel flow subject to thermal stratification is studied by Direct Numerical Simulation (DNS) of turbulence and Lagrangian Particle Tracking (LPT). Constant heat flux is maintained at the free surface of the channel, the bottom wall is adiabatic and the turbulent flow is driven by a pressure gradient. This archetypal flow setup mimics an environmentally plausible situation which can be found in terrestrial water bodies. The free surface turbulence characteristic of such flows has a strong influence on the distribution of the floaters: the objective of this work is to study the effect of different regimes of stable stratification on the surface distribution of floaters. The distribution of the floaters can possibly influence the transfer of chemical species across the water/atmosphere interface. Our results show that the modification of turbulence due to the thermal stratification strongly influences the settling velocity of floaters in the bulk of the flow. At the surface, stratification effects are also observed on the clustering of the floaters: the filamentary patterns of floaters observed in unstratified turbulence are progressively lost as thermal stratification increases, and the distribution of the floaters remains roughly two-dimensional.
Free surface height deformation upstream of an array of model marine hydrokinetic devices
NASA Astrophysics Data System (ADS)
Riley, D. R.; Beninati, M. L.; Volpe, M. A.; Krane, M.
2011-12-01
The changes in flow velocity and free surface height upstream of an array of model marine hydrokinetic (MHK) devices are measured in small-scale flume. These models consist of a stationary perforated plate mounted to a cylinder. The perforated plates mimic the blockage ratio and head loss effects of rotating blades on the flow, and the cylinders represent typical support structures for MHK devices. Experimental results are used to for two purposes. First, they will be used validate numerical simulations being conducted in parallel. Second, they will be used to determine a spatial arrangement for a field testing array that optimizes energy extraction. The operating hypothesis is that the proper spatial arrangement of a MHK array can modify the incident flow so as to increase the fluid energy incident upon some of the MHK devices in the array, increasing the energy produced by the array. The testing is conducted in the hydraulic flume facility (32 ft long, 4 ft wide, and 1.25 ft deep) in the Environmental Fluid Mechanics and Hydraulics Laboratory (EFM&HL) at Bucknell University. Within the facility a small-scale testing platform which consists of a nozzle insert is used to accelerate the flow into a test section (2 ft wide, 2.5 ft long, and 0.75 ft deep), where the cylinder arrays are located. The model MHK devices sizes are scaled using the test section floor turbulent boundary layer thickness. Flow field measurements are performed using an Acoustic Doppler Velocimeter. Mapping of the free surface is accomplished with an ultrasonic liquid level finder. Both devices are positioned using a programmable gantry system that has access to the entire test section. For each array configuration tested, both the free surface height and flow field upstream of the model MHK array are accurately mapped. These measurements are combined with the free surface height measurements to estimate the total fluid energy upstream of the array.
Yu, Mingzhe; McCulloch, William D; Beauchamp, Damian R; Huang, Zhongjie; Ren, Xiaodi; Wu, Yiying
2015-07-01
Integrating both photoelectric-conversion and energy-storage functions into one device allows for the more efficient solar energy usage. Here we demonstrate the concept of an aqueous lithium-iodine (Li-I) solar flow battery (SFB) by incorporation of a built-in dye-sensitized TiO2 photoelectrode in a Li-I redox flow battery via linkage of an I3(-)/I(-) based catholyte, for the simultaneous conversion and storage of solar energy. During the photoassisted charging process, I(-) ions are photoelectrochemically oxidized to I3(-), harvesting solar energy and storing it as chemical energy. The Li-I SFB can be charged at a voltage of 2.90 V under 1 sun AM 1.5 illumination, which is lower than its discharging voltage of 3.30 V. The charging voltage reduction translates to energy savings of close to 20% compared to conventional Li-I batteries. This concept also serves as a guiding design that can be extended to other metal-redox flow battery systems. PMID:26102317
Undulatory swimming on a free surface
NASA Astrophysics Data System (ADS)
Godoy-Diana, Ramiro; Ramananarivo, Sophie; Gann, Olivia; Thiria, Benjamin
2012-02-01
A wide variety of swimmers in nature use body undulations to generate a propulsive force, in part owing to the relative insensitivity of the principle of undulatory swimming to the value of the Reynolds number Re=UL/ν, which measures the relative importance of viscous and inertial forces in the flow considered (U and L being the typical speed and length of the animal, and ν the kinematic viscosity of the surrounding fluid). Here we study a flexible filament forced to oscillate by imposing a harmonic motion to one of its extremities (using magnetic interactions) and propelling itself at the surface of a water tank. This experiment serves as a canonical model for studying the interactions between an elastic structure undergoing complex deformations and the surrounding fluid.
Variational space-time (dis)continuous Galerkin method for nonlinear free surface water waves
NASA Astrophysics Data System (ADS)
Gagarina, E.; Ambati, V. R.; van der Vegt, J. J. W.; Bokhove, O.
2014-10-01
A new variational finite element method is developed for nonlinear free surface gravity water waves using the potential flow approximation. This method also handles waves generated by a wave maker. Its formulation stems from Miles' variational principle for water waves together with a finite element discretization that is continuous in space and discontinuous in time. One novel feature of this variational finite element approach is that the free surface evolution is variationally dependent on the mesh deformation vis-à-vis the mesh deformation being geometrically dependent on free surface evolution. Another key feature is the use of a variational (dis)continuous Galerkin finite element discretization in time. Moreover, in the absence of a wave maker, it is shown to be equivalent to the second order symplectic Störmer-Verlet time stepping scheme for the free-surface degrees of freedom. These key features add to the stability of the numerical method. Finally, the resulting numerical scheme is verified against nonlinear analytical solutions with long time simulations and validated against experimental measurements of driven wave solutions in a wave basin of the Maritime Research Institute Netherlands.
NASA Astrophysics Data System (ADS)
Zhang, Yin; Huang, Hu-Lin; Zhou, Xiao-Ming; Zhu, Gui-Ping; Zou, Yong
2016-04-01
Floating zone technique is a crucible-free process for growth of high quality single crystals. Unstable thermocapillary convection is a typical phenomenon during the process under microgravity. Therefore, it is very important to investigate the instability of thermocapillary convection in liquid bridges with deformable free-surface under microgravity. In this works, the Volume of Fluid (VOF) method is employed to track the free-surface movement. The results are presented as the behavior of flow structure and temperature distribution of the molten zone. The impact of Marangoni number ( Ma) is also investigated on free-surface deformation as well as the instability of thermocapillary convection. The free-surface exhibits a noticeable axisymmetric (but it is non-centrosymmetric) and elliptical shape along the circumferential direction. This specific surface shape presents a typical narrow `neck-shaped' structure with convex at two ends of the zone and concave at the mid-plane along the axial direction. At both θ = 0° and θ = 90°, the deformation ratio ξ increases rapidly with Ma at first, and then increases slowly. Moreover, the hydrothermal wave number m and the instability of thermocapillary convection increase with Ma.
An experimental investigation of spanwise vortices interacting with solid and free surfaces
NASA Astrophysics Data System (ADS)
Donnelly, Martin J.
Coherent vortices are generated in flow fields due to flow interaction with sharp solid surfaces. Such vortices generate significant disturbances in the flow and affect its further development. In this dissertation attention is focused on the interaction of vortices with solid or free liquid/air surfaces. We examine vortices with their axis parallel or normal to the surface. Three main cases were examined: the interaction of a vortex pair propagating towards a solid boundary, the interaction of spanwise vortices in a turbulent boundary layer, and finally the interaction of spanwise vortices with a flat-plate wake and a free liquid surface. These problems hold significance in several engineering applications, including investigations into trailing wing tip vortices and their interaction with the ground, vortical effects on the development of turbulent boundary layers and free surface signatures and their detection in ship/submarine wakes. Data are acquired with a laser Doppler velocimetry system (LDV) and with Particle-Image Velocimetry (PIV), using a high-speed digital video camera. The LDV system measures two components of velocity along appropriately chosen planes. Grids of data were acquired for different pitch rates of a disturbing flap that generates vortices. Phase-averaged vorticity and turbulence level contours are estimated and presented. It is found that vortices with diameter the order of the boundary layer quickly diffuse and disappear while their turbulent kinetic energy spreads uniformly across the entire boundary layer. Larger vortices have a considerably longer life span and in turn feed more vorticity into the boundary layer. Trailing edge vortices are generated in a water tunnel by sharp hinged motions of a flap. These vortices are allowed to reconnect with the free surface and mix with a turbulent free shear layer. The flow is conditionally sampled via frame grabbing of free surface shadowgraphs. It is found that the vortex core bends away from
NASA Astrophysics Data System (ADS)
Xu, Wenyu; Christenson, Michael; Fiflis, Peter; Curreli, Davide; Andruczyk, Daniel; Ruzic, David
2013-10-01
The application of liquid metal, especially liquid lithium has become an important topic for plasma facing component (PFC) design. A liquid PFC can effectively eliminate the erosion and thermal stress problems compared to the solid PFC while transferring heat and prolong the lifetime limit of the PFCs. A liquid lithium surface can also suppress the hydrogen isotopes recycling and getter the impurities in fusion reactors. The Lithium/metal infused trench (LiMIT) concept successfully proved that the thermoelectric effect can be utilized to drive liquid lithium flow within horizontally placed metallic open trenches in transverse magnetic field. A limiter based on this concept was tested in HT-7 and gave out positive results. However a broader application of this concept may require the trench be tilted or even placed vertically, for which strong capillary force caused by narrow trenches may be the solution. A new LiMIT design with very narrow trenches have been manufactured and tested in University of Illinois and related results will be presented. Based on this idea new limiters are designed for EAST and LTX and scheduled experiments on both devices will be discussed. This project is supported by DOE/ALPS contract: DEFG02- 99ER54515.
Zhao, Yu; Ding, Yu; Li, Yutao; Peng, Lele; Byon, Hye Ryung; Goodenough, John B; Yu, Guihua
2015-11-21
Electrical energy storage system such as secondary batteries is the principle power source for portable electronics, electric vehicles and stationary energy storage. As an emerging battery technology, Li-redox flow batteries inherit the advantageous features of modular design of conventional redox flow batteries and high voltage and energy efficiency of Li-ion batteries, showing great promise as efficient electrical energy storage system in transportation, commercial, and residential applications. The chemistry of lithium redox flow batteries with aqueous or non-aqueous electrolyte enables widened electrochemical potential window thus may provide much greater energy density and efficiency than conventional redox flow batteries based on proton chemistry. This Review summarizes the design rationale, fundamentals and characterization of Li-redox flow batteries from a chemistry and material perspective, with particular emphasis on the new chemistries and materials. The latest advances and associated challenges/opportunities are comprehensively discussed. PMID:26265165
Interfacial sliding near a free surface in a fibrous or layered composite during thermal cycling
NASA Technical Reports Server (NTRS)
Cox, B. N.
1990-01-01
This paper presents a simple shear lag model of interfacial sliding at a free surface in a layered or continuous fiber composite. The interface is characterized by a critical interfacial shear stress, tau0, which might represent the critical stress for frictional sliding at a weakly bonded interface, or the shear flow stress of a thin ductile interface layer at a well bonded interface. The history, during heating and cooling, of the relative normal displacement of the reinforcing inclusions and the matrix on a free surface cut normal to the inclusions is calculated and shown to depend on both the absolute value and the temperature dependence of tau0, as well as on the magnitudes of the bulk residual stresses. The variety of possible displacement histories suggests that they are a rich source of information about tau0 and the residual stresses.
NASA Astrophysics Data System (ADS)
Pottier, Basile; Talini, Laurence; Frétigny, Christian
We investigate the properties of liquids confined at nanometric scales from a solid wall with a new noninvasive technique. The optical technique used consists of measuring the height of fluctuations of the free surface, using the reflection of a laser beam on that surface. We hence measure the spontaneous thermal fluctuations of the free surfaces of liquids to probe their hydrodynamic boundary condition at a solid wall. The surface fluctuations of a silicon oil film could be described with a no-slip boundary condition for film thicknesses down to 20 nm. Oppositely, a 4 nm negative slip length had to be introduced to describe the behavior of n-hexadecane, consistently with previous surface force apparatus data on the same system. Our results demonstrate that at vanishing flow, a nanometric solid-like layer close to the wall may exist according to the nature of the liquid. currently at Laboratoire de Physique ENS Lyon.
High-flux neutron source based on a liquid-lithium target
Halfon, S.; Feinberg, G.; Paul, M.; Arenshtam, A.; Berkovits, D.; Kijel, D.; Nagler, A.; Eliyahu, I.; Silverman, I.
2013-04-19
A prototype compact Liquid Lithium Target (LiLiT), able to constitute an accelerator-based intense neutron source, was built. The neutron source is intended for nuclear astrophysical research, boron neutron capture therapy (BNCT) in hospitals and material studies for fusion reactors. The LiLiT setup is presently being commissioned at Soreq Nuclear research Center (SNRC). The lithium target will produce neutrons through the {sup 7}Li(p,n){sup 7}Be reaction and it will overcome the major problem of removing the thermal power generated by a high-intensity proton beam, necessary for intense neutron flux for the above applications. The liquid-lithium loop of LiLiT is designed to generate a stable lithium jet at high velocity on a concave supporting wall with free surface toward the incident proton beam (up to 10 kW). During off-line tests, liquid lithium was flown through the loop and generated a stable jet at velocity higher than 5 m/s on the concave supporting wall. The target is now under extensive test program using a high-power electron-gun. Up to 2 kW electron beam was applied on the lithium flow at velocity of 4 m/s without any flow instabilities or excessive evaporation. High-intensity proton beam irradiation will take place at SARAF (Soreq Applied Research Accelerator Facility) superconducting linear accelerator currently in commissioning at SNRC.
High-flux neutron source based on a liquid-lithium target
NASA Astrophysics Data System (ADS)
Halfon, S.; Feinberg, G.; Paul, M.; Arenshtam, A.; Berkovits, D.; Kijel, D.; Nagler, A.; Eliyahu, I.; Silverman, I.
2013-04-01
A prototype compact Liquid Lithium Target (LiLiT), able to constitute an accelerator-based intense neutron source, was built. The neutron source is intended for nuclear astrophysical research, boron neutron capture therapy (BNCT) in hospitals and material studies for fusion reactors. The LiLiT setup is presently being commissioned at Soreq Nuclear research Center (SNRC). The lithium target will produce neutrons through the 7Li(p,n)7Be reaction and it will overcome the major problem of removing the thermal power generated by a high-intensity proton beam, necessary for intense neutron flux for the above applications. The liquid-lithium loop of LiLiT is designed to generate a stable lithium jet at high velocity on a concave supporting wall with free surface toward the incident proton beam (up to 10 kW). During off-line tests, liquid lithium was flown through the loop and generated a stable jet at velocity higher than 5 m/s on the concave supporting wall. The target is now under extensive test program using a high-power electron-gun. Up to 2 kW electron beam was applied on the lithium flow at velocity of 4 m/s without any flow instabilities or excessive evaporation. High-intensity proton beam irradiation will take place at SARAF (Soreq Applied Research Accelerator Facility) superconducting linear accelerator currently in commissioning at SNRC.
NASA Astrophysics Data System (ADS)
Spitans, Sergejs; Baake, Egbert; Nacke, Bernard; Jakovics, Andris
2016-02-01
By means of external coupling between electromagnetic (EM) problem in ANSYS and hydrodynamic problem in FLUENT, a numerical model for the liquid metal free surface flow in an alternate EM field has been developed and verified in the first part of the article. Volume of Fluid ( VOF) algorithm has been used for tracking of free surface. In this work, improved performance of the model is presented. General validation of the VOF algorithm is performed by comparison of the calculated free oscillations of the liquid column to its analytical solution. The 3D/ VOF calculation of coupled EM field and free surface flow with Large Eddy Simulation turbulence description for the first time is applied for modeling of conventional EM levitation. Calculation results are compared with 2D/ VOF and 3D/ VOF models that use less precise k- ɛ and k- ω SST turbulence formulations. Obtained time-averaged droplet shapes are used for single-phase flow calculations with different turbulence models and free-slip/no-slip velocity conditions at the fixed free surface for validation of the flow. Meanwhile, series of levitation melting experiments are performed for verification of the simulated droplet shapes. In conclusion, parameter impact on the fully developed flow and the levitated droplet shape is discussed.
Front instability and energy of the free surface
NASA Astrophysics Data System (ADS)
Beltrame, Philippe
2014-05-01
In recent years, there has been a proliferation of research devoted to the formation of preferential flow paths occurring without macroscopic heterogeneity of the porous media. DiCarlo (2013) points out the connection between "overshoot" and the front instability. Extension of the standard Richards equation is required to capture this phenomenon. In most of the improvements, interfacial phenomena as the triple line at the front are considered. For instance, velocity dependent contact angle (Wang et al., 2013) or contact angle hysteresis (Rätz and Schweizer, 2012) allow to simulate successfully the instability. Another approach proposed by Cueto-Felgueroso and Juanes (2009) introduces a macroscopic surface tension related to the existence of the water/air interface. As previously, the simulation of an advancing front displays physical looking fingering displacements. The goal of this contribution is to better understand the role of the different surface energies in the emergence of the front instability. We propose a model involving both the macroscopic surface tension and the soil wettability. This latter allows to define a contact angle and possibly hysteresis using heterogeneous wettability (Beltrame et al., 2011). Therefore, we employ the phase field approach developed by Felgueroso and Juanes, 2009 to which we add a free energy term corresponding to the wettability: a disjoining or conjoining pressure resulting from effective molecular interactions between the substrate and the free surface (DeGennes, 1985). The difference with the classical suction pressure is the hydrophobic behavior for ultra-thin film (small water saturation). Such a water repellency was recently estimated in the soil (Diamantopoulos et al. 2013). Stability analysis of an advancing front in an uniform porous media shows that macroscopic surface tension and wettability may independently produce the instability growth. In contrast, for a front stopping when reaching the layers interface of
Bubble Bursting at a Free Surface in a Closed Domain
NASA Astrophysics Data System (ADS)
Liu, Nian-Nian; Zhang, Shuai; Wang, Shi-Ping
2016-06-01
When a charge explodes underwater near a free surface, a bubble would be generated and the surface pushed up very high. Experiments have shown that the motion of the spike lags a lot behind the bubble motion. Many studies only focus on the nonlinear interaction between the bubble and free surface while the water waves afterward is mainly studied based on the linear theory. The nonlinear motion of the water wave after the bubble pulsation is seldom studied. In this study, we concerns the interaction between underwater explosion generated bubble and a free surface and its bursting at a free surface in a closed domain. Suppose that the fluid outside the bubble is incompressible, non-viscous and irrotational and the velocity potential satisfies the Laplace equation. Boundary integral method is used to solve the Laplace equation for the velocity potential. The bubble content is described by an adiabatic law. The whole process of the bubble motion and subsequently the water wave propagation will be simulated in this paper. Particular attention will be focused on the phenomenon of water wave propagation in a closed domain.
NASA Astrophysics Data System (ADS)
Sobolev, V. Ya.; Sinyavsky, V. V.
2014-12-01
Sub-megawatt space NPP use lithium as a coolant and niobium alloy as a structural material. In order to refine the lithium-niobium technology of the material and design engineering, lithium-niobium loops were worked out in RSC Energia, and they were tested at a working temperature of lithium equal to 1070-1300 K. In order to measure the lithium flow and pressure, special gauges were developed, which made possible the calibration and checkout of the loops without their dismantling. The paper describes the architecture of the electromagnetic flowmeter and the electromagnetic vibrating-wire pressure transducer (gauge) for lithium coolant in the nuclear power plant cooling systems. The operating principles of these meters are presented. Flowmeters have been developed for channel diameters ranging from 10 to 100 mm, which are capable of measuring lithium flows in the range of 0.1 to 30 L/s with the error of 3% for design calibration and 1% for volume graduation. The temperature error of the pressure transducers does not exceed 0.4% per 100 K; the nonlinearity and hysteresis of the calibration curve do not exceed 0.3 and 0.4%, respectively. The transducer applications are illustrated by the examples of results obtained from tests on the NPP module mockup and heat pipes of a radiation cooler.
Passive manipulation of free-surface instability by deformable solid bilayers.
Sahu, Shivam; Shankar, V
2016-07-01
This study deals with the elastohydrodynamic coupling that occurs in the flow of a liquid layer down an inclined plane lined with a deformable solid bilayer and its consequences on the stability of the free surface of the liquid layer. The fluid is Newtonian and incompressible, while the linear elastic constitutive relation has been considered for the deformable solid bilayer, and the densities of the fluid and the two solids are kept equal. A temporal linear stability analysis is carried out for this coupled solid-fluid system. A long-wave asymptotic analysis is employed to obtain an analytical expression for the complex wavespeed in the low wave-number regime, and a numerical shooting method is used to solve the coupled set of governing differential equations in order to obtain the stability criterion for arbitrary values of the wave number. In a previous work on plane Couette flow past an elastic bilayer, Neelmegam et al. [Phys. Rev. E 90, 043004 (2014)PLEEE81539-375510.1103/PhysRevE.90.043004] showed that the instability of the flow can be significantly influenced by the nature of the solid layer, which is adjacent to the liquid layer. In stark contrast, for free-surface flow past a bilayer, our long-wave asymptotic analysis demonstrates that the stability of the free-surface mode is insensitive to the nature of the solid adjacent to the liquid layer. Instead, it is the effective shear modulus of the bilayer G_{eff} (given by H/G_{eff}=H_{1}/G_{1}+H_{2}/G_{2}, where H=H_{1}+H_{2} is the total thickness of the solid bilayer, H_{1} and H_{2} are the thicknesses of the two solid layers, and G_{1} and G_{2} are the shear moduli of the two solid layers) that determines the stability of the free surface in the long-wave limit. We show that for a given Reynolds number, the free-surface instability is stabilized when G_{eff} decreases below a critical value. At finite wave numbers, our numerical solution indicates that additional instabilities at the free surface and
Passive manipulation of free-surface instability by deformable solid bilayers
NASA Astrophysics Data System (ADS)
Sahu, Shivam; Shankar, V.
2016-07-01
This study deals with the elastohydrodynamic coupling that occurs in the flow of a liquid layer down an inclined plane lined with a deformable solid bilayer and its consequences on the stability of the free surface of the liquid layer. The fluid is Newtonian and incompressible, while the linear elastic constitutive relation has been considered for the deformable solid bilayer, and the densities of the fluid and the two solids are kept equal. A temporal linear stability analysis is carried out for this coupled solid-fluid system. A long-wave asymptotic analysis is employed to obtain an analytical expression for the complex wavespeed in the low wave-number regime, and a numerical shooting method is used to solve the coupled set of governing differential equations in order to obtain the stability criterion for arbitrary values of the wave number. In a previous work on plane Couette flow past an elastic bilayer, Neelmegam et al. [Phys. Rev. E 90, 043004 (2014), 10.1103/PhysRevE.90.043004] showed that the instability of the flow can be significantly influenced by the nature of the solid layer, which is adjacent to the liquid layer. In stark contrast, for free-surface flow past a bilayer, our long-wave asymptotic analysis demonstrates that the stability of the free-surface mode is insensitive to the nature of the solid adjacent to the liquid layer. Instead, it is the effective shear modulus of the bilayer Geff (given by H /Geff=H1/G1+H2/G2 , where H =H1+H2 is the total thickness of the solid bilayer, H1 and H2 are the thicknesses of the two solid layers, and G1 and G2 are the shear moduli of the two solid layers) that determines the stability of the free surface in the long-wave limit. We show that for a given Reynolds number, the free-surface instability is stabilized when Geff decreases below a critical value. At finite wave numbers, our numerical solution indicates that additional instabilities at the free surface and the liquid-solid interface can be induced by wall
Experimental investigation of a free-surface turbulent jet with Coanda effect
NASA Astrophysics Data System (ADS)
Miozzi, M.; Lalli, F.; Romano, G. P.
2010-07-01
The deviation of a jet from the straight direction due to the presence of a lateral wall is investigated from the experimental point of view. This flow condition is known as Coanda jet (from the Romanian aerodynamicist Henry Marie Coanda who discovered and applied it at the beginning of XXth century) or offset jet. The objective of the work is to detail the underlying mechanisms of such a phenomenon aiming to use it as a flow control method at polluted river flows mouth. To do this, a large laboratory free-surface tank with an incoming channel has been set up and velocity field measurements are performed by Optical Flow methods (namely Feature Tracking). Preliminary tests on the well-known free jet configuration without any marine structure ( i.e. lateral wall) are performed to allow comparison with free jet scaling and self-similar solutions. The presence of the free-surface gives rise to centerline velocity decay which is lower than in free unbounded plane or circular jets due to the vertically limited ambient fluid entrainment. In the second part of the paper, the effect of a lateral wall on the jet configuration is examined by placing it at different lateral distances from the jet outlet. The resulting velocity fields clearly show an inclined Coanda jet with details which seems to depend on the lateral wall distance itself. The analysis of self-similarity along the inclined jet direction reveals that for wall distances larger than 5 jet widths this dependence almost disappears.
Free surface oscillations of liquid hydrogen in microgravity conditions
NASA Astrophysics Data System (ADS)
Schmitt, Sebastian; Dreyer, Michael E.
2015-12-01
Experiments were performed to investigate the reorientation behavior and axial sloshing of liquid parahydrogen in a partly filled right circular cylinder due to a step reduction of gravity. Different temperature gradients along the cylinder wall in vertical direction were imposed to examine the influence of a wall superheat on the free surface. Experiments were conducted in the drop tower at the University of Bremen which provides a microgravity time of 4.7 s and a compensated gravity environment of 10-6g0 (acceleration due to gravity). The thermal preparation of the experiments allowed to create defined wall temperature gradients and a stratified or homogenous liquid temperature distribution. Several sensors along the cylinder wall and in the vapor region monitored the temperature. The pressure inside the experiment was recorded and visual access was enabled by an endoscope. The experiments showed that the wall superheat has an influence on the free surface as well as on the temperature and pressure evolutions.
Investigating Free-surface, MHD Instabilities in Liquid Metals
NASA Astrophysics Data System (ADS)
Brumfiel, Geoff; Ji, Hantao; Zweben, Stewart
1999-11-01
The addition of the Lorentz force (j × B) to MHD fluids, such as liquid metals, could introduce new instabilities to free-surface motions. An experiment is under development at PPPL that will contribute to the understanding of these instabilities. This experiment is designed to explore how instabilities form and propagate on the surface of liquid metals. Metals with low melting points and reasonable conductive properties (such as Ga) are melted in an eight inch, square Pyrex container. The container is placed in a large magnetic field (up to 5 kG) and a perpendicular current is set up in the metal. Waves are generated using a simple wave driving apparatus. A one dimensional, diode camera is used to monitor wave propagation perpendicular to the B-field. Initial results will be presented and discussed. This experiment will provide information that could eventually be used to better control free-surface motions in liquid Li walls in fusion reactors.
Reduced Viscosity of Free Surface in Entangled Polymer Melt Films
NASA Astrophysics Data System (ADS)
Koga, Tad; Li, C.; Endoh, M.; Koo, J.; Rafailovich, M.; Narayanan, S.; Lee, D.; Lurio, L.; Sinha, S.
2010-03-01
The dynamics of polymer chains near the surface of a melt and within thin films remains a subject of inquiry along with the nature of the glass transition in these systems. By embedding ``dilute'' gold nanoparticles in single polystyrene thin films as ``markers'', we could probe the local viscosity of the free surface at temperatures far above the glass transition temperature (Tg). The technique used was X-ray photon correlation spectroscopy with resonance-enhanced X-ray scattering. The results clearly showed the viscosity was about 30 % lower than the rest of the film. We found that this reduction is strongly associated with chain entanglements at the free surface rather than the reduction in Tg.
NASA Astrophysics Data System (ADS)
Pottier, B.; Frétigny, C.; Talini, L.
2015-06-01
We investigate the properties of nanometric liquid films with a new noninvasive technique. We measure the spontaneous thermal fluctuations of the free surfaces of liquids to probe their hydrodynamic boundary condition at a solid wall. The surface fluctuations of a silicon oil film could be described with a no-slip boundary condition for film thicknesses down to 20 nm. Oppositely, a 4 nm negative slip length had to be introduced to describe the behavior of n -hexadecane, consistently with previous surface force apparatus data on the same system. Our results demonstrate that at vanishing flow a nanometric solidlike layer close to the wall may exist according to the nature of the liquid.
Large-eddy simulation of free-surface decaying turbulence with dynamic subgrid-scale models
NASA Astrophysics Data System (ADS)
Salvetti, M. V.; Zang, Y.; Street, R. L.; Banerjee, S.
1997-08-01
This paper describes large-eddy simulations of decaying turbulence in an open channel, using different dynamic subgrade-scale models, viz. the dynamic model of Germano et al. [Phys. Fluids A 3, 1790 (1991)] (DSM), the dynamic mixed model in Zang et al. [Phys. Fluids A 5, 3186 (1993)] (DMM), and the dynamic two-parameter model of Salvetti and Banerjee [Phys. Fluids 7, 2831 (1995)] (DTM). These models are incorporated in a finite-volume solver of the Navier-Stokes equations. A direct numerical simulation of this flow conducted by Pan and Banerjee [Phys. Fluids 7, 1649 (1995)] showed that near the free surface turbulence has a quasi-two-dimensional behavior. Moreover, the quasi-two-dimensional region increases in thickness with the decay time, although the structure remains three-dimensional in the central regions of the flow. The results of the large-eddy simulations show that both the DMM and the DTM are able to reproduce the features of the decay process observed in the direct simulation and to handle the anisotropic nature of the flow. Nevertheless, the addition of the second model coefficient in the DTM improves the agreement with the direct simulation. When the DSM is used, significant discrepancies are observed between the large-eddy and the direct simulations during the decay process at the free surface.
NASA Astrophysics Data System (ADS)
Olmez, O.; Ozbulut, M.; Yildiz, M.; Goren, O.
2016-06-01
The present study investigates the vortical and nonlinear effects in the roll motion of a 2-D body with square cross-sections by using Smoothed Particle Hydrodynamics (SPH). A 2-D rigid body with square cross-section is taken into account for the benchmark study and subjected to the oscillatory roll motion with a given angular frequency. The governing equations are continuity equation and Euler's equation with artificial viscosity term. Weakly Compressible SPH (WCSPH) scheme is employed for the discretization of the governing equations. Velocities of the fluid particles are updated by means of XSPH+Artificial Particle Displacement (VXSPH+APD) algorithm. In this method only the free surface fluid particles are subjected to VXSPH algorithm while the APD algorithm is employed for the fully populated flow regions. The hybrid usage of numerical treatment keeps free surface particles together by creating an artificial surface tension on the free surface. VXSPH+APD is a proven numerical treatment to provide the most accurate results for this type of free surface flows (Ozbulut et al. 2014). The results of the present study are compared with those of the experimental studies as well as with those of the numerical methods obtained from the current literature.
NASA Astrophysics Data System (ADS)
Adamski, Samantha; Aliseda, Alberto
2012-11-01
Marine Hydrokinetic (MHK) turbines are a growing area of research in the renewable energy field because tidal currents are a highly predictable clean energy source. The presence of a free surface may influence the flow around the turbine and in the wake, critically affecting turbine performance and environmental effects through modification of wake physical variables. The characteristic Froude number that control these processes is still a matter of controversy, with the channel depth and turbine's depth, blade tip depth and diameter as potential candidates for a length scale used in literature. We use the Volume of Fluid model to track the free surface dynamics in a RANS simulation with a BEMT model of the turbine to understand the physics of the wake-free surface interactions. Pressure and flow rate boundary conditions for channel's inlet, outlet and air side have been tested in an effort to determine the optimum set of simulation conditions for MHK turbines in rivers or estuaries. Stability and accuracy in terms of power extraction and kinetic and potential energy budgets are considered. The goal of this research is to determine, quantitatively in non dimensional parameter space, the limit between negligible and significant free surface effects on MHK turbine analysis. Supported by DOE through the National Northwest Marine Renewable Energy Center.
NASA Astrophysics Data System (ADS)
Sarada Sree, Atchutuni; Tanaji, Kamble; Poulami, Chakraborty; Fotedar, R. K.; Rajendra Kumar, E.; Suri, A. K.; Platacis, E.; Ziks, A.; Bucenieks, I.; Poznjaks, A.; Shisko, A.
2014-08-01
To study the corrosion of P-91 (9% chromium and 1% molybdenum) material with lead-lithium (Pb-Li) eutectic, two experiments were carried out in a forced convection loop, at eutectic temperature of 550 °C. The first experiment was carried out at a velocity of 15 cm s-1 for 1000 h and the second experiment, at a velocity of 30 cm s-1 for 2700 h. In both the experiments, P-91 sample coupons were exposed to Pb-Li flow in the presence and absence of magnetic field. Samples were analyzed using an optical microscope, scanning electron microscope (SEM) and electron probe micro-analyzer (EPMA). Micro-Vickers hardness testing was also carried out. Dissolution of elements into liquid metal is the main corrosion mechanism. Iron and chromium were selectively getting leached out from the near-surface region (˜4 µm) in the first experiment and molybdenum and manganese were also found leaching from a greater depth in the second experiment. The samples kept in the magnetic field showed a higher corrosion rate (˜320 µm/year) as compared with the corrosion rate (˜200 µm/year) of the samples kept in non-magnetic field regions. Hardness of the exposed samples was lower than the unexposed samples in both the experiments. Hardness was found to be low in the near-surface region for all the samples in both the experiments.
Free Surface Wave Interaction with a Horizontal Cylinder
NASA Astrophysics Data System (ADS)
Oshkai, P.; Rockwell, D.
1999-10-01
Classes of vortex formation from a horizontal cylinder adjacent to an undulating free-surface wave are characterized using high-image-density particle image velocimetry. Instantaneous representations of the velocity field, streamline topology and vorticity patterns yield insight into the origin of unsteady loading of the cylinder. For sufficiently deep submergence of the cylinder, the orbital nature of the wave motion results in multiple sites of vortex development, i.e., onset of vorticity concentrations, along the surface of the cylinder, followed by distinctive types of shedding from the cylinder. All of these concentrations of vorticity then exhibit orbital motion about the cylinder. Their contributions to the instantaneous values of the force coefficients are assessed by calculating moments of vorticity. It is shown that large contributions to the moments and their rate of change with time can occur for those vorticity concentrations having relatively small amplitude orbital trajectories. In a limiting case, collision with the surface of the cylinder can occur. Such vortex-cylinder interactions exhibit abrupt changes in the streamline topology during the wave cycle, including abrupt switching of the location of saddle points in the wave. The effect of nominal depth of submergence of the cylinder is characterized in terms of the time history of patterns of vorticity generated from the cylinder and the free surface. Generally speaking, generic types of vorticity concentrations are formed from the cylinder during the cycle of the wave motion for all values of submergence. The proximity of the free surface, however, can exert a remarkable influence on the initial formation, the eventual strength, and the subsequent motion of concentrations of vorticity. For sufficiently shallow submergence, large-scale vortex formation from the upper surface of the cylinder is inhibited and, in contrast, that from the lower surface of the cylinder is intensified. Moreover
Micromechanics of slip bands on a free surface
NASA Technical Reports Server (NTRS)
Lin, S. R.; Lin, T. H.
1976-01-01
A micromechanics analysis for the formation and propagation of slip bands on the free surface of a polycrystal under monotonic loading is presented. For the growth of slip bands, the analysis satisfies the conditions of both equilibrium and displacement continuity, as well as the relation between slip and the resolved shear stress throughout the polycrystal. Numerical calculations show how the microstress field causes the concentration of plastic deformation in discrete sliding bands and give results which are in good qualitative agreement with known slip band observations on aluminum single crystals.
The motion of elliptic cylinder under free surface
NASA Astrophysics Data System (ADS)
Kostikov, V. K.; Makarenko, N. I.
2016-06-01
A problem on generation of unsteady nonlinear waves on the surface of an infinitely deep ideal fluid due to the motion of a submerged elliptical cylinder is considered. It is supposed that the cylinder can rotate in addition to translational two-dimensional motion. The initial formulation of the problem is reduced to an integrodifferential system of equations for the functions defining the free surface shape, the normal and tangential components of velocity on the free boundary. The small-time asymptotics of the solution is constructed in the case of the cylinder that moves with a constant acceleration from rest.
Study on Applicability of Numerical Simulation to Evaluation of Gas Entrainment From Free Surface
Kei Ito; Takaaki Sakai; Hiroyuki Ohshima
2006-07-01
An onset condition of gas entrainment (GE) due to free surface vortex has been studied to establish a design of fast breeder reactor with higher coolant velocity than conventional designs, because the GE might cause the reactor operation instability and therefore should be avoided. The onset condition of the GE has been investigated experimentally and theoretically, however, dependency of the vortex type GE on local geometry configuration of each experimental system and local velocity distribution has prevented researchers from formulating the universal onset condition of the vortex type GE. A real scale test is considered as an accurate method to evaluate the occurrence of the vortex type GE, but the real scale test is generally expensive and not useful in the design study of large and complicated FBR systems, because frequent displacement of inner equipments accompanied by the design change is difficult in the real scale test. Numerical simulation seems to be promising method as an alternative to the real scale test. In this research, to evaluate the applicability of the numerical simulation to the design work, numerical simulations were conducted on the basic experimental system of the vortex type GE. This basic experiment consisted of rectangular flow channel and two important equipments for vortex type GE in the channel, i.e. vortex generation and suction equipments. Generated vortex grew rapidly interacting with the suction flow and the grown vortex formed a free surface dent (gas core). When the tip of the gas core or the bubbles detached from the tip of the gas core reached the suction mouth, the gas was entrained to the suction tube. The results of numerical simulation under the experimental conditions were compared to the experiment in terms of velocity distributions and free surface shape. As a result, the numerical simulation showed qualitatively good agreement with experimental data. The numerical simulation results were similar to the experimental
Mathematical and computational studies of equilibrium capillary free surfaces
NASA Technical Reports Server (NTRS)
Albright, N.; Chen, N. F.; Concus, P.; Finn, R.
1977-01-01
The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.
Plankton dynamics in thermally-stratified free-surface turbulence
NASA Astrophysics Data System (ADS)
Lovecchio, Salvatore; Soldati, Alfredo
2015-11-01
Thermal stratification induced by solar heating near the ocean-atmosphere interface influences the transfer fluxes of heat, momentum and chemical species across the interface. Due to thermal stratification, a region of large temperature gradients (thermocline) may form with strong consequences for the marine ecosystem. In particular, the thermocline is believed to prevent phytoplankton from reaching the well-lit surface layer, where they can grow through the process of photosynthesis. In this paper, we use a DNS-based Eulerian-Lagrangian approach to examine the role of stratification on phytoplankton dynamics in thermally-stratified free-surface turbulence. We focus on gyrotactic self-propelled phytoplankton cells, considering different stratification levels (quantified by the Richardson number) and different gyro tactic re-orientation times. We show that the modulation of turbulent fluctuations induced by stable stratification has a strong effect on the orientation and distribution of phytoplankton, possibly leading to trapping of some species within the thermocline. Specifically, we observe the appearance of a depletion layer just below the free-surface as stratification increases, accompanied by a reduction in the vertical stability of phytoplankton cells.
Free surface BCP self-assembly process characterization with CDSEM
NASA Astrophysics Data System (ADS)
Levi, Shimon; Weinberg, Yakov; Adan, Ofer; Klinov, Michael; Argoud, Maxime; Claveau, Guillaume; Tiron, Raluca
2016-03-01
A simple and common practice to evaluate Block copolymers (BCP) self-assembly performances, is on a free surface wafer. With no guiding pattern the BCP designed to form line space pattern for example, spontaneously rearranges to form a random fingerprint type of a pattern. The nature of the rearrangement is dictated by the physical properties of the BCP moieties, wafer surface treatment and the self-assembly process parameters. Traditional CDSEM metrology algorithms are designed to measure pattern with predefined structure, like linespace or oval via holes. Measurement of pattern with expected geometry can reduce measurement uncertainty. Fingerprint type of structure explored in this dissertation, poses a challenge for CD-SEM measurement uncertainty and offers an opportunity to explore 2D metrology capabilities. To measure this fingerprints we developed a new metrology approach that combines image segmentation and edge detection to measure 2D pattern with arbitrary rearrangement. The segmentation approach enabled to quantify the quality of the BCP material and process, detecting 2D attributes such as: CD and CDU at one axis, and number of intersections, length and number of PS fragments, etched PMMA spaces and donut shapes numbers on the second axis. In this paper we propose a 2D metrology to measure arbitrary BCP pattern on a free surface wafer. We demonstrate experimental results demonstrating precision data, and characterization of PS-b-PMMA BCP, intrinsic period L0 = 38nm (Arkema), processed at different bake time and temperatures.
Magnetic gauge for free surface velocities due to rock blasts
NASA Astrophysics Data System (ADS)
Ashuach, Yecheskel; Gissis, Itai; Avinadav, Chen
2013-06-01
We developed a simple magnetic gauge for measuring free surface velocities of rock materials in the range of 0.1-20 m/s. The gauge consists of two elements: a NdFeB magnet and a pick-up coil. The coil is attached to the free surface at the point of interest. The magnet is placed a few centimeters away from the coil on its central axis, intact from the rock. Rock surface movement due to blast loading induces current in the coil due to change of the magnetic flux. The coil velocity is deduced from the measured current using a computational code. The gauge was tested and validated in a set of free-falling experiments. We present velocity measurements from various blast experiments in limestone and reinforced concrete, using both the magnetic gauge and a Doppler interferometer. The results obtained from the two measurement techniques were in good agreement during a few milliseconds. The magnetic gauge is cheap and very simple to operate, and therefore favorable for mapping the velocity distribution at multiple points of interest on the surface.
Review of free-surface MHD experiments and modeling.
Molokov, S.; Reed, C. B.
2000-06-02
This review paper was prepared to survey the present status of analytical and experimental work in the area of free surface MHD and thus provide a well informed starting point for further work by the Advanced Limiter-diverter Plasma-facing Systems (ALPS) program. ALPS were initiated to evaluate the potential for improved performance and lifetime for plasma-facing systems. The main goal of the program is to demonstrate the advantages of advanced limiter/diverter systems over conventional systems in terms of power density capability, component lifetime, and power conversion efficiency, while providing for safe operation and minimizing impurity concerns for the plasma. Most of the work to date has been applied to free surface liquids. A multi-disciplinary team from several institutions has been organized to address the key issues associated with these systems. The main performance goals for advanced limiters and diverters are a peak heat flux of >50 MW/m{sup 2}, elimination of a lifetime limit for erosion, and the ability to extract useful heat at high power conversion efficiency ({approximately}40%). The evaluation of various options is being conducted through a combination of laboratory experiments, modeling of key processes, and conceptual design studies.
46 CFR 170.290 - Free surface correction for damage stability calculations.
Code of Federal Regulations, 2011 CFR
2011-10-01
... 46 Shipping 7 2011-10-01 2011-10-01 false Free surface correction for damage stability...) SUBDIVISION AND STABILITY STABILITY REQUIREMENTS FOR ALL INSPECTED VESSELS Free Surface § 170.290 Free surface correction for damage stability calculations. (a) When doing the damage stability calculations required...
Self-induced sloshing caused by an upward round jet impinging on the free surface
Madarame, Haruki; Okamoto, Koji; Iida, Masao
1995-11-01
The demand for compact design of FBR made it necessary to study the interaction between high speed coolant flow and the free surface. Self-induced sloshing was found to occur with an upward round jet making a surface swell at the impinging point. The conditions under which the sloshing occurred were investigated experimentally together with the sloshing modes and frequencies. Without exception the induced mode had only one diametrical node, which was accompanied by the lateral motion of the swell. The mode was replaced by lower one with increasing the inlet-surface distance and the velocity, the reason of which was explained by the following assumption; the oscillation of water columns on both sides of the jet connected by a virtual cross flow through the jet, i.e., jet-flutter, was coupled with the sloshing and supplied energy to it.
Dunn, Jennifer B.; Gaines, Linda; Barnes, Matthew; Sullivan, John L.; Wang, Michael
2014-01-01
This document contains material and energy flows for lithium-ion batteries with an active cathode material of lithium manganese oxide (LiMn₂O₄). These data are incorporated into Argonne National Laboratory’s Greenhouse gases, Regulated Emissions, and Energy use in Transportation (GREET) model, replacing previous data for lithium-ion batteries that are based on a nickel/cobalt/manganese (Ni/Co/Mn) cathode chemistry. To identify and determine the mass of lithium-ion battery components, we modeled batteries with LiMn₂O₄ as the cathode material using Argonne’s Battery Performance and Cost (BatPaC) model for hybrid electric vehicles, plug-in hybrid electric vehicles, and electric vehicles. As input for GREET, we developed new or updated data for the cathode material and the following materials that are included in its supply chain: soda ash, lime, petroleum-derived ethanol, lithium brine, and lithium carbonate. Also as input to GREET, we calculated new emission factors for equipment (kilns, dryers, and calciners) that were not previously included in the model and developed new material and energy flows for the battery electrolyte, binder, and binder solvent. Finally, we revised the data included in GREET for graphite (the anode active material), battery electronics, and battery assembly. For the first time, we incorporated energy and material flows for battery recycling into GREET, considering four battery recycling processes: pyrometallurgical, hydrometallurgical, intermediate physical, and direct physical. Opportunities for future research include considering alternative battery chemistries and battery packaging. As battery assembly and recycling technologies develop, staying up to date with them will be critical to understanding the energy, materials, and emissions burdens associated with batteries.
Dunn, J.B.; Gaines, L.; Barnes, M.; Wang, M.; Sullivan, J.
2012-06-21
This document contains material and energy flows for lithium-ion batteries with an active cathode material of lithium manganese oxide (LiMn{sub 2}O{sub 4}). These data are incorporated into Argonne National Laboratory's Greenhouse gases, Regulated Emissions, and Energy use in Transportation (GREET) model, replacing previous data for lithium-ion batteries that are based on a nickel/cobalt/manganese (Ni/Co/Mn) cathode chemistry. To identify and determine the mass of lithium-ion battery components, we modeled batteries with LiMn{sub 2}O{sub 4} as the cathode material using Argonne's Battery Performance and Cost (BatPaC) model for hybrid electric vehicles, plug-in hybrid electric vehicles, and electric vehicles. As input for GREET, we developed new or updated data for the cathode material and the following materials that are included in its supply chain: soda ash, lime, petroleum-derived ethanol, lithium brine, and lithium carbonate. Also as input to GREET, we calculated new emission factors for equipment (kilns, dryers, and calciners) that were not previously included in the model and developed new material and energy flows for the battery electrolyte, binder, and binder solvent. Finally, we revised the data included in GREET for graphite (the anode active material), battery electronics, and battery assembly. For the first time, we incorporated energy and material flows for battery recycling into GREET, considering four battery recycling processes: pyrometallurgical, hydrometallurgical, intermediate physical, and direct physical. Opportunities for future research include considering alternative battery chemistries and battery packaging. As battery assembly and recycling technologies develop, staying up to date with them will be critical to understanding the energy, materials, and emissions burdens associated with batteries.
Free surface microfluidic/SERS for detection of gas-phase DNT
NASA Astrophysics Data System (ADS)
Piorek, Brian D.; Lee, Seung Joon; Judy, Nick; Meinhart, Carl D.; Moskovits, Martin; Fountain, Augustus; Christesen, Steven; Guicheteau, Jason
2010-04-01
A novel microfluidic/SERS platform has been developed for real time sensing of 2,4-DNT. The fundamental research is being conducted at UCSB, commercialized by SpectraFluidics, and validated at ECBC. The system leverages phenomena at multiple length scales, ranging from tens of micrometers to several nanometers. The key enabling technology is a newly developed invention termed Free-Surface Fluidics (FSF), where one or more fluidic surfaces are confined by surface tension forces, and exposed to the surrounding atmosphere. The free-surface fluidic architecture is combined with surface-enhanced Raman spectroscopy (SERS) for detection of 2,4-DNT. Once 2,4-DNT analyte molecules are absorbed into the flow, they can interact with gold or silver colloidal particles. This architecture allows for analysis and deterministic control of SERS 'hot spot' aggregation, which can increase Raman scattering signal strength by up to 10 orders in magnitude. We have successfully measured DNT vapor at concentrations as low as ~1 ppb. This sensitivity value is confirmed by orthogonal measurements using GC-mass spectroscopy at ECBC.
Modeling Outburst Flooding as a Turbulent Hydraulic Fracture Parallel to a Nearby Free Surface
NASA Astrophysics Data System (ADS)
Tsai, V. C.; Rice, J. R.
2009-12-01
Meltwater generated at the surface and base of glaciers and ice sheets is known to have a large impact on how ice masses behave dynamically, but much is still unknown about the physical processes responsible for how this meltwater drains out of the glacier. For example, little attention has been paid to short-timescale processes like turbulent hydraulic fracture, which is likely an important mechanism by which drainage channels initially form when water pressures are high. In recent work (Tsai and Rice [Fall AGU, 2008; JGR subm., 2009]), we have constructed a model of this turbulent hydraulic fracture process in which over-pressurized water is assumed to flow turbulently through a crack, leading to crack growth. However, one important limitation of this prior work is that it only strictly applies in the limit of short crack length 2L compared to glacier height H, whereas relevant observations of supraglacial lake drainage, jokulhlaups and sub-glacial lake-to-lake transport episodes do not fall in this regime. Here, we improve somewhat upon this model by explicitly accounting for a nearby free surface. We accomplish this by applying the approach of Erdogan et al. [Meth. Anal. Sol. Crack Prob., 1973] to numerically calculate elastic displacements consistent with crack pressure distribution for a crack near a free surface, and use these results as before to simultaneously satisfy the governing fluid, elastic and fracture equations. Our results are analogous to the zero fracture toughness results of Zhang et al. [Int. J. Numer. Anal. Meth. Geomech., 2005], but applied to the case of turbulent flow rather than laminar flow of a Newtonian viscous fluid. Our new results clarify the importance of the free surface and potentially explain discrepancies between our previous modeling results and observations of supraglacial lake drainage by Das et al. [Science, 2008]. However, the numerical challenges increase as 2L becomes comparable to or much larger than H. We hope to
Modeling Outburst Flooding as a Turbulent Hydraulic Fracture Parallel to a Nearby Free Surface
NASA Astrophysics Data System (ADS)
Tsai, Victor; Rice, James
2010-05-01
Meltwater generated at the surface and base of glaciers and ice sheets is known to have a large impact on how ice masses behave dynamically, but much is still unknown about the physical processes responsible for how this meltwater drains out of the glacier. For example, little attention has been paid to short-timescale processes like turbulent hydraulic fracture, which is likely an important mechanism by which drainage channels initially form when water pressures are high. In recent work (Tsai and Rice [Fall AGU, 2008; JGR subm., 2009]), we have constructed a model of this turbulent hydraulic fracture process in which over-pressurized water is assumed to flow turbulently through a crack, leading to crack growth. However, one important limitation of this prior work is that it only strictly applies in the limit of short crack length, 2L, compared to glacier height, H, whereas relevant observations of supraglacial lake drainage, jokulhlaups and sub-glacial lake-to-lake transport episodes do not fall in this regime. Here, we improve somewhat upon this model by explicitly accounting for a nearby free surface. We accomplish this by applying the approach of Erdogan et al. [Meth. Anal. Sol. Crack Prob., 1973] to numerically calculate elastic displacements consistent with crack pressure distribution for a crack near a free surface, and use these results as before to simultaneously satisfy the governing fluid, elastic and fracture equations. Our results are analogous to the zero fracture toughness results of Zhang et al. [Int. J. Numer. Anal. Meth. Geomech., 2005], but applied to the case of turbulent flow rather than laminar flow of a Newtonian viscous fluid. Our new results clarify the importance of the free surface and potentially explain discrepancies between our previous modeling results and observations of supraglacial lake drainage by Das et al. [Science, 2008]. However, the numerical challenges increase as 2L becomes comparable to or much larger than H. We hope to
Recent Advances in Applicability of TEMHD Driven Liquid Lithium as a Fusion Relevant PFC
NASA Astrophysics Data System (ADS)
Szott, Matthew; Fiflis, Peter; Kalathiparambil, Kishor; Ruzic, David N.
2015-11-01
Liquid lithium displays increasing promise as a replacement to solid plasma facing components (PFC) in fusion device applications. Liquid PFCs reduce erosion and thermal stress damage, prolonging device lifetime, while lithium has been shown to decrease edge recycling, reduce impurities, and enhance plasma performance. The Liquid Metal Infused Trench (LiMIT) concept developed at UIUC successfully demonstrates horizontal and vertical thermoelectric magnetohydrodynamic (TEMHD) flow of liquid lithium through metal trenches for use as a PFC. Installed in the HT-7 tokamak and at the Magnum-PSI linear plasma device, the system performed effectively in fusion relevant conditions. In high heat flux tests, trench dry-out was observed, which exposes solid trench material due to higher TEMHD force on the area with the highest heat flux. A 3D free surface fluid model of dry-out and experimental tests conducted to mitigate the detrimental effect are described. The final designs for the upcoming test of LiMIT as a limiter for the EAST tokamak are discussed, along with velocity characteristics of steady-state TEMHD driven flow through the LiMIT system inclined up to 180 degrees from horizontal, which is necessary for broad applicability of a liquid lithium PFC system.
Pottier, B; Frétigny, C; Talini, L
2015-06-01
We investigate the properties of nanometric liquid films with a new noninvasive technique. We measure the spontaneous thermal fluctuations of the free surfaces of liquids to probe their hydrodynamic boundary condition at a solid wall. The surface fluctuations of a silicon oil film could be described with a no-slip boundary condition for film thicknesses down to 20 nm. Oppositely, a 4 nm negative slip length had to be introduced to describe the behavior of n-hexadecane, consistently with previous surface force apparatus data on the same system. Our results demonstrate that at vanishing flow a nanometric solidlike layer close to the wall may exist according to the nature of the liquid. PMID:26196646
The threshold origin of solutocapillary Marangoni convection on a bounded free surface
NASA Astrophysics Data System (ADS)
Viviani, Antonio; Denisova, Maria; Kostarev, Konstantin; Zuev, Andrew
2014-09-01
The Marangoni convection is the key mechanism responsible for the motion of non-isothermal fluid with a free surface under microgravity conditions. It is generally assumed that the surface of a Newtonian fluid begins to move at any arbitrary small shear stress. However, the experiments show that under actual conditions the capillary fluid flow can be initiated by a certain (threshold) stress. For subthreshold shear stresses generated by the surface forces or volumetric flows, the fluid surface remains immovable. Such a behavior of the surface is caused by the existence of an adsorption layer formed by uncontrolled surface-active impurities contained in the fluid. The composition and amount of impurities absorbed on the surface are dependent on the degree of purification of the fluid and the value of its surface tension. The condition of a threshold initiation of the capillary motion has been determined experimentally for the case of development of the solutal Marangoni convection after placing a droplet of the surfactant liquor on the interface of water, aqueous surfactant solutions of lower concentrations and some organic fluids. A group of monohydric alcohols and carboxylic acids were used as the surfactants. It has been found that the value of the threshold difference in surfactant concentration rapidly increases with a decrease in fluid initial concentration and the characteristic dimension of the free surface. Due to the latter fact the threshold effects are of crucial importance for a variety of problems in microfluidics. It has been found that a threshold origin of a capillary motion is also peculiar to the liquid-liquid interface. The relationships between the threshold Marangoni number and the surface activity of a surfactant have been determined.
Glasses crystallize rapidly at free surfaces by growing crystals upward
Sun, Ye; Zhu, Lei; Kearns, Kenneth L.; Ediger, Mark D.; Yu, Lian
2011-01-01
The crystallization of glasses and amorphous solids is studied in many fields to understand the stability of amorphous materials, the fabrication of glass ceramics, and the mechanism of biomineralization. Recent studies have found that crystal growth in organic glasses can be orders of magnitude faster at the free surface than in the interior, a phenomenon potentially important for understanding glass crystallization in general. Current explanations differ for surface-enhanced crystal growth, including released tension and enhanced mobility at glass surfaces. We report here a feature of the phenomenon relevant for elucidating its mechanism: Despite their higher densities, surface crystals rise substantially above the glass surface as they grow laterally, without penetrating deep into the bulk. For indomethacin (IMC), an organic glass able to grow surface crystals in two polymorphs (α and γ), the growth front can be hundreds of nanometers above the glass surface. The process of surface crystal growth, meanwhile, is unperturbed by eliminating bulk material deeper than some threshold depth (ca. 300 nm for α IMC and less than 180 nm for γ IMC). As a growth strategy, the upward-lateral growth of surface crystals increases the system’s surface energy, but can effectively take advantage of surface mobility and circumvent slow growth in the bulk. PMID:21444775
Label-free surface plasmon sensing towards cancer diagnostics
NASA Astrophysics Data System (ADS)
Sankaranarayanan, Goutham
The main objective of this thesis is to develop a conventional, home-built SPR bio-sensor to demonstrate bio-sensing applications. This emphasizes the understanding of basic concepts of Surface Plasmon Resonance and various interrogation techniques. Intensity Modulation was opted to perform the label-free SPR bio-sensing experiments due to its cost-efficient and compact setup. Later, label-free surface plasmon sensing was carried out to study and understand the bio-molecular interactions between (1). BSA and Anti BSA molecules and (2). Exosome/Liposome on thin metal (Au) films. Exosomes are cell-derived vesicles present in bodily fluids like blood, saliva, urine, epididymal fluid containing miRNAs, RNA, proteins, etc., at stable quantities during normal health conditions. The exosomes comprise varied constituents based on their cell origin from where they are secreted and is specific to that particular origin. However an exacerbated release is observed during tumor or cancer conditions. This increased level of exosomes present in the sample, can be detected using the SPR bio-sensor demonstrated in this thesis and effective thickness of adsorption on Au surface can be estimated. Also, chemically synthesized liposome particles were studied to determine if they can generate an equivalent sensor response to that of exosomes to consider them as an alternate. Finally a 10ppb Mercury (Hg) sensing was performed as part of Environment Monitoring application and results have been tabulated and compared.
Spherical Demons: Fast Diffeomorphic Landmark-Free Surface Registration
Yeo, B.T. Thomas; Sabuncu, Mert R.; Vercauteren, Tom; Ayache, Nicholas; Fischl, Bruce; Golland, Polina
2010-01-01
We present the Spherical Demons algorithm for registering two spherical images. By exploiting spherical vector spline interpolation theory, we show that a large class of regularizors for the modified Demons objective function can be efficiently approximated on the sphere using iterative smoothing. Based on one parameter subgroups of diffeomorphisms, the resulting registration is diffeomorphic and fast. The Spherical Demons algorithm can also be modified to register a given spherical image to a probabilistic atlas. We demonstrate two variants of the algorithm corresponding to warping the atlas or warping the subject. Registration of a cortical surface mesh to an atlas mesh, both with more than 160k nodes requires less than 5 minutes when warping the atlas and less than 3 minutes when warping the subject on a Xeon 3.2GHz single processor machine. This is comparable to the fastest non-diffeomorphic landmark-free surface registration algorithms. Furthermore, the accuracy of our method compares favorably to the popular FreeSurfer registration algorithm. We validate the technique in two different applications that use registration to transfer segmentation labels onto a new image: (1) parcellation of in-vivo cortical surfaces and (2) Brodmann area localization in ex-vivo cortical surfaces. PMID:19709963
Fabrication of Adhesive Lenses Using Free Surface Shaping
NASA Astrophysics Data System (ADS)
Hoheisel, D.; Kelb, C.; Wall, M.; Roth, B.; Rissing, L.
2013-09-01
Two approaches for fabricating polymer lenses are presented in this paper. Both are based on filling circular holes with UV curing adhesives. Initially, the viscous adhesive material creates a liquid and spherical free surface due to its own surface tension. This shape is then preserved by curing with UV-hardening light. For the first approach, the holes are generated in a 4 inch Si-wafer by deep reactive ion etching (DRIE) and for the second, a polydimethylsiloxane (PDMS) mould is manufactured. Three types of UV-curing adhesives are investigated (NOA 61, NOA 88 and NEA 121 by Norland Products). Preliminary to the determination of the lens curvature, a contact angle goniometer is used for taking side view images of the lenses. The radius of curvature is then extracted via image processing with the software MATLABÂ®. Furthermore, the surface roughness of the PDMS mould and the generated lenses is measured with a white light interferometer to characterize the casting process. The resolution power of the generated lenses is evaluated by measurement of their point spread functions (psf) and modulation transfer functions (mtf), respectively.