Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Free surface flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records

Finite element simulation of temperature dependent free surface flows

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.

Engelman, M. S.↗

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite element solutions of free surface flows

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.

Zarda, P. R.↗

Continuation in a parameter - Experience with viscous and free surface flows

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.

Kheshgi, H. S.↗

Distribution of resistive body-force in curved free-surface flow

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.

Sivakumaran, N. S.↗

Computational and theoretical analysis of free surface flow in a thin liquid film under zero and normal gravity

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.

Faghri, Amir↗

Application of the methods of gas dynamics to water flows with free surface I : flows with no energy dissipation

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.

Preiswerk, Ernst↗

RIPPLE - A new model for incompressible flows with free surfaces

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.

Kothe, D. B.↗

An experimental study of free surface deformation in oscillatory thermocapillary flow

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.

Lin, Jianlian↗

Computing Incompressible Flows With Free Surfaces

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.

Kothe, D.↗

Project JOVE

The goal of this project is to investigate new areas of research pertaining to free surface-interface fluids mechanics and/or microgravity which have potential commercial applications. This paper presents an introduction to ferrohydrodynamics (FHD), and discusses some applications. Also, computational methods for solving free surface flow problems are presented in detail. Both have diverse applications in industry and in microgravity fluids applications. Three different modeling schemes for FHD flows are addressed and the governing equations, including Maxwell's equations, are introduced. In the area of computational modeling of free surface flows, both Eulerian and Lagrangian schemes are discussed. The state of the art in computational methods applied to free surface flows is elucidated. In particular, adaptive grids and re-zoning methods are discussed. Additional research results are addressed and copies of the publications produced under the JOVE Project are included.

Lyell, M. J.↗

Numerical simulation of fluid flow and heat transfer in a thin liquid film over a stationary and rotating disk and comparison with experimental data

In the first section, improvements in the theoretical model and computational procedure for the prediction of film height and heat-transfer coefficient of the free surface flow of a radially-spreading thin liquid film adjacent to a flat horizontal surface of finite extent are presented. Flows in the presence and absence of gravity are considered. Theoretical results are compared to available experimental data with good agreement. In the presence of gravity, a hydraulic jump is present, isolating the flow into two regimes: supercritical upstream from the jump and subcritical downstream of it. In this situation, the effects of surface tension are important near the outer edge of the disk where the fluid experiences a free fall. A region of flow separation is present just downstream of the jump. In the absence of gravity, no hydraulic jump or separated flow region is present. The variation of the heat-transfer coefficient for flows in the presence and absence of gravity are also presented. In the second section, the results of a numerical simulation of the flow field and associated heat transfer coefficients are presented for the free surface flow of a thin liquid film adjacent to a horizontal rotating disk. The computation was performed for different flow rates and rotational velocities using a 3-D boundary-fitted coordinate system. Since the geometry of the free surface is unknown and dependent on flow rate, rate of rotation, and other parameters, an iterative procedure had to be used to ascertain its location. The computed film height agreed well with existing experimental measurements. The flow is found to be dominated by inertia near the entrance and close to the free surface and dominated by centrifugal force at larger radii and adjacent to the disk. The rotation enhances the heat transfer coefficient by a significant amount.

Faghri, Amir↗

Numerical and experimental analysis of a thin liquid film on a rotating disk related to development of a spacecraft absorption cooling system

The numerical and experimental analysis of a thin liquid film on a rotating and a stationary disk related to the development of an absorber unit for a high capacity spacecraft absorption cooling system, is described. The creation of artificial gravity by the use of a centrifugal field was focused upon in this report. Areas covered include: (1) One-dimensional computation of thin liquid film flows; (2) Experimental measurement of film height and visualization of flow; (3) Two-dimensional computation of the free surface flow of a thin liquid film using a pressure optimization method; (4) Computation of heat transfer in two-dimensional thin film flow; (5) Development of a new computational methodology for the free surface flows using a permeable wall; (6) Analysis of fluid flow and heat transfer in a thin film in the presence and absence of gravity; and (7) Comparison of theoretical prediction and experimental data. The basic phenomena related to fluid flow and heat transfer on rotating systems reported here can also be applied to other areas of space systems.

Faghri, Amir↗

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗