Search NASA⌕ Search

SEARCH · Search NASA

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Coupling methods between finite element-based Boussinesq-type wave and particle-based free-surface flow models [Coupling methods between finite element-based Boussinesq-typewave and particle-based free-surface flow models]

We develop one-way coupling methods between a Boussinesq-type wave model based on the discontinuous Galerkin finite element method and a free-surface flow model based on a mesh-free particle method to strike a balance between accuracy and computational cost. In our proposed model, computation of the wave model in the global domain is conducted first, and the nonconstant velocity profiles in the vertical direction are reproduced by using its results. Computation of the free-surface flow is performed in a local domain included within the global domain with interface boundaries that move along the reproduced velocity field in a Lagrangian fashion. To represent the moving interfaces, we used a polygon wall boundary model for mesh-free particle methods. Verification and validation tests of our proposed model are performed, and results obtained by the model are compared with theoretical values and experimental results to show its accuracy and applicability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

High-resolution simulation of free-surface flow and tracer retention over streambeds with ripples

This study presents a novel high-resolution simulation of free-surface flow and tracer retention over a streambed with ripples based on varying ripple morphologies, surface hydraulics and the transport of a tracer pulse from surface water to surface dead zone. For the simulations, the computational fluid dynamics (CFD) model OpenFOAM was used to solve the three-dimensional Navier-Stokes equations in combination with an implemented transport equation. Pressure gradients at the streambed were used to account for hyporheic exchange, assuming water flow from high pressure zones to low pressure zones. Flow velocities, ripple sizes and spacing showed to significantly affect these pressure gradients, but also the transport of a passive tracer at the streambed, which was not investigated so far. Due to the velocity field, large parts of the tracer mass were transported alongside the main stream above the ripples. Tracer mass reaching the space between the ripples was temporarily retained due to low velocities and recirculations. It was shown that the retention is depending on the ripple size and space between the ripples as well as on the flow velocity. Decreasing ripple sizes and higher flow velocities lead to a smaller tracer retention. Furthermore we showed that the ripple length to height ratio controls the generation of recirculation zones which affect the residence time of the tracer significantly. Ripple spacing leads to temporarily higher tracer concentration at the streambed, but smaller tracer retention. We conclude that the impact of the streambed morphology on the hydraulics in combination with tracer retention should be addressed for a comprehensive understanding of compound movement, exchange and transformation within the hyporheic zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2017

Using the Multilayer Free-Surface Flow Model to Solve Wave Problems

A method is presented for changing over from a single-layer shallow-water model to a multilayer model with hydrostatic pressure profile and, then, to a multilayer model with nonhydrostatic pressure profile. The method does not require complex procedures for solving the discrete Poisson’s equation and features high computation efficiency. The results of validating the algorithm against experimental data critical for the numerical dissipation of the numerical scheme are presented. Examples are considered.

42 ENGINEERING↗
DOE OSTI2020

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↗
DOE OSTI2015

Casting Process Modeling and Validation

Computational simulations of gravity casting for metallic U-Pu-Zr nuclear fuel rods were performed with the free surface flow solver Flow3d. Two different casting apparatus designs were explored for studying the effects on the final cast rod profiles. Experimental data for the casting problem obtained from Idaho National Laboratory (INL) was validated with numerical simulations on one of the designs considered. The volume of the rods cast, length of the rods, and temperature profiles were predicted through the transient casting process and were found to match closely with experimental data. A sensitivity analysis with vent locations, pressure, mesh refinement, and effects of surface tension at the interface were explored. This work serves as a validation benchmark for the casting problem and provides guidelines when performing simulation exercises for modeling casting. This paper also provides the mathematical framework to enabling modeling of the casting process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Interfacial gauge methods for incompressible fluid dynamics

Designing numerical methods for incompressible fluid flow involving moving interfaces, for example, in the computational modeling of bubble dynamics, swimming organisms, or surface waves, presents challenges due to the coupling of interfacial forces with incompressibility constraints. A class of methods, denoted interfacial gauge methods, is introduced for computing solutions to the corresponding incompressible Navier-Stokes equations. These methods use a type of "gauge freedom" to reduce the numerical coupling between fluid velocity, pressure, and interface position, allowing high-order accurate numerical methods to be developed more easily. Making use of an implicit mesh discontinuous Galerkin framework, developed in tandem with this work, high-order results are demonstrated, including surface tension dynamics in which fluid velocity, pressure, and interface geometry are computed with fourth-order spatial accuracy in the maximum norm. Applications are demonstrated with two-phase fluid flow displaying fine-scaled capillary wave dynamics, rigid body fluid-structure interaction, and a fluid-jet free surface flow problem exhibiting vortex shedding induced by a type of Plateau-Rayleigh instability. The developed methods can be generalized to other types of interfacial flow and facilitate precise computation of complex fluid interface phenomena.

42 ENGINEERING↗
DOE OSTI2015

Plethora of transitions during breakup of liquid filaments

Significance Fluid flows, governed by nonlinear equations, permit formation of singularities. Often, singularities are artifacts of neglecting physical effects. However, free-surface flows exhibit observable singularities including filament pinch-off. As filaments thin, slightly (highly) viscous filaments are expected from theory to transition from an inertial (viscous) regime where viscosity (density) is negligible to an inertial–viscous regime where viscous and inertial effects are important. Previous works show this transition either does not occur or occurs for filament radii well below theoretical predictions. We demonstrate that thinning filaments unexpectedly pass through a number of intermediate transient regimes, thereby delaying onset of the final regime. The findings raise the question if similar dynamical transitions arise in problems that are not necessarily hydrodynamic in nature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity-Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three-Dimensional Smoothed Particle Hydrodynamics Simulations

In this work, we study unsaturated gravity-driven free-surface flows on a synthetic surface intersected by a horizontal fracture and the importance of flow modes (droplet; rivulet) on the partitioning behavior at the fracture intersection. We present laboratory experiments, three-dimensional smoothed particle hydrodynamics (SPH) simulations using a heavily parallelized code, and an analytical solution. The flow-rate-dependent mode switching from droplets to rivulets is reproduced by the SPH model, and the transition ranges agree with the laboratory experiments. We show that flow modes heavily influence the bypass behavior. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. This behavior is demonstrated for a multi-inlet laboratory setup where the inlet- specific flow rate is chosen so that either a droplet or rivulet flow persists. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture in- flow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), the flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer. For the case of rivulet flow, we provide an analytical solution that demonstrates the existence of classical Washburn flow within the horizontal fracture.

unsaturated flow, fractured media, smoothed partic↗
DOE OSTI2016

Stabilising falling liquid film flows using feedback control

Falling liquid films become unstable due to inertial effects when the fluid layer is sufficiently thick or the slope sufficiently steep. This free surface flow of a single fluid layer has industrial applications including coating and heat transfer, which benefit from smooth and wavy interfaces, respectively. Here, we discuss how the dynamics of the system are altered by feedback controls based on observations of the interface height, and supplied to the system via the perpendicular injection and suction of fluid through the wall. In this study, we model the system using both Benney and weighted-residual models that account for the fluid injection through the wall. We find that feedback using injection and suction is a remarkably effective control mechanism: the controls can be used to drive the system towards arbitrary steady states and travelling waves, and the qualitative effects are independent of the details of the flow modelling. Furthermore, we show that the system can still be successfully controlled when the feedback is applied via a set of localised actuators and only a small number of system observations are available, and that this is possible using both static (where the controls are based on only the most recent set of observations) and dynamic (where the controls are based on an approximation of the system which evolves over time) control schemes. This study thus provides a solid theoretical foundation for future experimental realisations of the active feedback control of falling liquid films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗