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

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↗

Application of the Analogy Between Water Flow with a Free Surface and Two-dimensional Compressible Gas Flow

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.

Orlin, W James↗

Application of the Analogy Between Water Flow with a Free Surface and Two-Dimensional Compressible Gas Flow

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.

Orlin, W James↗

Three-Dimensional, Multi-Phase Unsteady Flow of Water Undergoing Multiple Flow Regimes (Free Surface, Droplet, and Fluid Film)

Upgrades to the Mobile Launcher Rainbird nozzles were implemented for the Artemis II launch to increase water flow rates and improve the mobile launcher deck water coverage while minimizing the water throw distance to avoid impacting the Space Launch System (SLS) rocket nozzles. This analysis describes the use of a three-dimensional, multi-phase, free surface transient computational fluid dynamics (CFD) model to simulate the water flowing out through the Rainbird nozzle and ensure compliance with current Artemis requirements. Flow rates, pressures, velocities, nozzle throw distances and coverage areas were analyzed to ensure that system requirements are met. The completed results were used to make necessary modifications to the Ignition Overpressure Protection and Sound Suppression (IOP/SS) system Rainbird design by reshaping the nozzle used in the system.

CFD↗

Implementation of nonreflective boundary condition at the outflow boundary

Numerical techniques for the solution of free surface flows about a disturbance on or near a free surface were developed. Any disturbance moving on or near a liquid surface generates surface waves that propagate away from the disturbance. The potential flow problems in which the evolution of free surface waves is the principal phenomenon of interest are considered. The distinct feature of these problems is that the position of the surface is also unknown and is to be determined as part of the solution. Two boundary conditions are required at the free surface since there are two unknown functions there, the free surface height and the potential function. The required boundary conditions are satisfied by the kinematic condition which states that no flow can cross the surface and the dynamic condition which matches the fluid pressure to the pressure acting on the surface through Bernoulli's equation.

Hall, D. R.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Numerical simulation of fluid flow and heat transfer in a thin liquid film over a rotating disk

The results of a numerical simulation of the flow field and associated heat transfer coefficient are presented for the free surface flow of a thin liquid film adjacent to a horizontal rotating disk. The computation has been performed for different flow rates and rotational velocities using a three-dimensional 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 interative procedure had to be used to ascertain its location. The computed film height agreed well with existing experimental measurements. The flow was 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 enhanced the heat transfer coefficient by a significant amount.

Rahman, M. M.↗

New concepts for Reynolds stress transport equation modeling of inhomogeneous flows

The ability to model turbulence near solid walls and other types of boundaries is important in predicting complex engineering flows. Most turbulence modeling has concentrated either on flows which are nearly homogeneous or isotropic, or on turbulent boundary layers. Boundary layer models usually rely very heavily on the presence of mean shear and the production of turbulence due to that mean shear. Most other turbulence models are based on the assumption of quasi-homogeneity. However, there are many situations of engineering interest which do not involve large shear rates and which are not quasi-homogeneous or isotropic. Shear-free turbulent boundary layers are the prototypical example of such flows, with practical situations being separation and reattachment, bluff body flow, high free-stream turbulence, and free surface flows. Although these situations are not as common as the variants of the flat plate turbulent boundary layer, they tend to be critical factors in complex engineering situations. The models developed are intended to extend classical quasi-homogeneous models into regions of large inhomogeneity. These models do not rely on the presence of mean shear or production, but are still applicable when those additional effects are included. Although the focus is on shear-free boundary layers as tests for these models, results for standard shearing boundary layers are also shown.

Perot, J. Blair↗

Modeling of Turbulent Flow in Electromagnetically Levitated Metal Droplets

This article details an effort to improve the understanding and prediction of turbulent flow inside a droplet of molten metal levitated in an electromagnetic field. It is shown that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is in the transitional regime between laminar and turbulent flow. Past research efforts have used laminar, enhanced viscosity, and k-epsilon turbulence models to describe the flow. The method highlighted in our study is the renormalization group (RNG) algorithm. We show that an accurate description of the turbulent eddy viscosity is critical in order to obtain realistic velocity fields, and that the turbulent eddy viscosity cannot be uniform in levitated droplets. The RNG method does not impose isotropic length or time scales on the flow field, thus allowing such nonuniform features to be captured. A number of other materials processing applications exhibit similarly complex flow characteristics, such as highly recirculating, transitional, and free surface flows, for which this modeling approach may prove useful.

Berry, S.↗

Dynamics of electromechanical flow structures.

Free-surface gravity flows and capillary wicking provide examples of flow structures with fluid partially ducted at free surfaces by external forces. Wall-less electromechanical flow structures are developed which have a similar nature, but with polarization forces providing the orientation at free surfaces. Like their mechanical counterparts, these have the ability to ingest liquid or expel vapor through their walls. The structures consist of electrodes running in the flow direction z with slowly varying cross sections in a plane transverse to the flow. A formulation is given of the long-wave nonlinear (principal mode) dynamics, with use made of energy functions to represent a broad class of possible mechanical and electrical structure geometries.

Jones, T. B., Jr.↗

Simulation of three-dimensional liquid sloshing flows using a strongly implicit calculation procedure

A coupled strongly implicit solution strategy for unsteady three-dimensional free surface flows has been developed based on an artificial compressibility formulation for the incompressible Navier-Stokes equations. A pseudotime term has been used in the continuity equation to permit time accurate calculations to be achieved. The scheme appears capable of tracking the free surface reasonably accurately inside a partially-filled spherical container undergoing a general rotating motion characteristic of that experienced by a spin-stabilized satellite. Five different free surface calculations have been presented. Some of the results exhibit an interesting Reynolds number dependent oscillatory behavior which is believed to be physical although no experimental results appears to be available for verification to date.

Chen, K.-H.↗

The flow of a thin liquid film on a stationary and rotating disk. II - Theoretical prediction

The existing theoretical models are improved and a systematic procedure to compute the free surface flow of a thin liquid film is suggested. The solutions for axisymmetric radial flow on a stationary horizontal disk and for the disk rotating around its axis are presented. The theoretical predictions are compared with the experimental data presented in Part I of this report. The analysis shows results for both supercritical and subcritical flows and the flow structure in the vicinity of a hydraulic jump which isolates these two flow types. The detailed flow structure in a hydraulic jump was computed and shown to contain regions of separation including a 'surface roller'. The effects of surface tension are found to be important near the outer edge of the disk where the fluid experiences a free fall. At other locations, the surface tension is negligible. For a rotating disk, the frictional resistance in the angular direction is found to be as important as that in the radial direction.

Rahman, M. M.↗

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗

2-D moving mesh modeling of lithium dryout in open surface liquid metal flow applications

Liquid lithium displays increasing promise as a replacement for solid plasma facing components (PFCs) in fusion device applications. The Liquid Metal Infused Trench (LiMIT) system, developed at the University of Illinois (UIUC), has demonstrated how thermoelectric magnetohydrodynamics (TEMHD) can be harnessed to drive liquid lithium flow in an open surface PFC. However, in the highest heat flux applications, large local acceleration is created via TEMHD, and the sudden increase in velocity can cause the liquid level to expose the underlying solid, eliminating the protective benefits of the lithium. In order to study potential mitigation strategies, a 2-D COMSOL Multiphysics model was developed using the moving mesh module to capture free surface flow. The model depicts the development of the dryout phenomenon for 2 test cases – slow (1 cm/s) and medium (10 cm/s) flow in 5 mm deep trenches – including the liquid level reduction under the high heat flux and the pileup of slower flow both upstream and downstream of the heat stripe. The effectiveness of trench shaping dryout mitigation strategies is examined. For the slow flow case, it is shown that a 1.8 mm ledge placed under the heat stripe will stop dryout, and for the medium flow case, a 2.7 mm ledge is required to mitigate the effect. This model can be used to identify strategies for increasing the viable heat load for open surface liquid lithium PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heat transfer to a thin liquid film with a free surface

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.

Rahman, M. M.↗

Numerical solution of flow problems using body-fitted coordinate systems

The paper deals with numerically generated boundary-fitted coordinate systems. This procedure eliminates the shape of the boundaries as a complicating factor and allows the flow about arbitrary boundaries to be treated essentially as easily as that about simple boundaries. The technique of boundary-fitted coordinate systems is based on a method of automatic numerical generation of a general curvilinear coordinate system having a coordinate line coincident with each boundary of a general multiconnected region involving any number of arbitrarily shaped boundaries. Once the curvilinear coordinate system is generated, any partial differential system of interest may be solved on the coordinate system by transforming the equations and solving the resulting system in finite-difference approximation on the rectangular transformed plane. Attention is given to the types of boundary-fitted coordinate systems, coordinate system control, operation of the coordinate codes, solution of partial differential equations, application to free-surface flow, and other applications of interest.

Thompson, J. F.↗