Straight velocity boundaries in the lattice Boltzmann method

NASA Astrophysics Data System (ADS)

Latt, Jonas; Chopard, Bastien; Malaspinas, Orestis; Deville, Michel; Michler, Andreas

2008-05-01

Various ways of implementing boundary conditions for the numerical solution of the Navier-Stokes equations by a lattice Boltzmann method are discussed. Five commonly adopted approaches are reviewed, analyzed, and compared, including local and nonlocal methods. The discussion is restricted to velocity Dirichlet boundary conditions, and to straight on-lattice boundaries which are aligned with the horizontal and vertical lattice directions. The boundary conditions are first inspected analytically by applying systematically the results of a multiscale analysis to boundary nodes. This procedure makes it possible to compare boundary conditions on an equal footing, although they were originally derived from very different principles. It is concluded that all five boundary conditions exhibit second-order accuracy, consistent with the accuracy of the lattice Boltzmann method. The five methods are then compared numerically for accuracy and stability through benchmarks of two-dimensional and three-dimensional flows. None of the methods is found to be throughout superior to the others. Instead, the choice of a best boundary condition depends on the flow geometry, and on the desired trade-off between accuracy and stability. From the findings of the benchmarks, the boundary conditions can be classified into two major groups. The first group comprehends boundary conditions that preserve the information streaming from the bulk into boundary nodes and complete the missing information through closure relations. Boundary conditions in this group are found to be exceptionally accurate at low Reynolds number. Boundary conditions of the second group replace all variables on boundary nodes by new values. They exhibit generally much better numerical stability and are therefore dedicated for use in high Reynolds number flows.

Computation of viscous flows over airfoils, including separation, with a coupling approach

NASA Technical Reports Server (NTRS)

Leballeur, J. C.

1983-01-01

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Definition of boundary and initial conditions in the analysis of saturated ground-water flow systems - An introduction

USGS Publications Warehouse

Franke, O. Lehn; Reilly, Thomas E.; Bennett, Gordon D.

1987-01-01

Accurate definition of boundary and initial conditions is an essential part of conceptualizing and modeling ground-water flow systems. This report describes the properties of the seven most common boundary conditions encountered in ground-water systems and discusses major aspects of their application. It also discusses the significance and specification of initial conditions and evaluates some common errors in applying this concept to ground-water-system models. An appendix is included that discusses what the solution of a differential equation represents and how the solution relates to the boundary conditions defining the specific problem. This report considers only boundary conditions that apply to saturated ground-water systems.

Use of an Accurate DNS Particulate Flow Method to Supply and Validate Boundary Conditions for the MFIX Code

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhi-Gang Feng

2012-05-31

The simulation of particulate flows for industrial applications often requires the use of two-fluid models, where the solid particles are considered as a separate continuous phase. One of the underlining uncertainties in the use of the two-fluid models in multiphase computations comes from the boundary condition of the solid phase. Typically, the gas or liquid fluid boundary condition at a solid wall is the so called no-slip condition, which has been widely accepted to be valid for single-phase fluid dynamics provided that the Knudsen number is low. However, the boundary condition for the solid phase is not well understood. Themore » no-slip condition at a solid boundary is not a valid assumption for the solid phase. Instead, several researchers advocate a slip condition as a more appropriate boundary condition. However, the question on the selection of an exact slip length or a slip velocity coefficient is still unanswered. Experimental or numerical simulation data are needed in order to determinate the slip boundary condition that is applicable to a two-fluid model. The goal of this project is to improve the performance and accuracy of the boundary conditions used in two-fluid models such as the MFIX code, which is frequently used in multiphase flow simulations. The specific objectives of the project are to use first principles embedded in a validated Direct Numerical Simulation particulate flow numerical program, which uses the Immersed Boundary method (DNS-IB) and the Direct Forcing scheme in order to establish, modify and validate needed energy and momentum boundary conditions for the MFIX code. To achieve these objectives, we have developed a highly efficient DNS code and conducted numerical simulations to investigate the particle-wall and particle-particle interactions in particulate flows. Most of our research findings have been reported in major conferences and archived journals, which are listed in Section 7 of this report. In this report, we will present a brief description of these results.« less

Impact of the Diurnal Cycle of the Atmospheric Boundary Layer on Wind-Turbine Wakes: A Numerical Modelling Study

NASA Astrophysics Data System (ADS)

Englberger, Antonia; Dörnbrack, Andreas

2018-03-01

The wake characteristics of a wind turbine for different regimes occurring throughout the diurnal cycle are investigated systematically by means of large-eddy simulation. Idealized diurnal cycle simulations of the atmospheric boundary layer are performed with the geophysical flow solver EULAG over both homogeneous and heterogeneous terrain. Under homogeneous conditions, the diurnal cycle significantly affects the low-level wind shear and atmospheric turbulence. A strong vertical wind shear and veering with height occur in the nocturnal stable boundary layer and in the morning boundary layer, whereas atmospheric turbulence is much larger in the convective boundary layer and in the evening boundary layer. The increased shear under heterogeneous conditions changes these wind characteristics, counteracting the formation of the night-time Ekman spiral. The convective, stable, evening, and morning regimes of the atmospheric boundary layer over a homogeneous surface as well as the convective and stable regimes over a heterogeneous surface are used to study the flow in a wind-turbine wake. Synchronized turbulent inflow data from the idealized atmospheric boundary-layer simulations with periodic horizontal boundary conditions are applied to the wind-turbine simulations with open streamwise boundary conditions. The resulting wake is strongly influenced by the stability of the atmosphere. In both cases, the flow in the wake recovers more rapidly under convective conditions during the day than under stable conditions at night. The simulated wakes produced for the night-time situation completely differ between heterogeneous and homogeneous surface conditions. The wake characteristics of the transitional periods are influenced by the flow regime prior to the transition. Furthermore, there are different wake deflections over the height of the rotor, which reflect the incoming wind direction.

Time dependent inflow-outflow boundary conditions for 2D acoustic systems

NASA Technical Reports Server (NTRS)

Watson, Willie R.; Myers, Michael K.

1989-01-01

An analysis of the number and form of the required inflow-outflow boundary conditions for the full two-dimensional time-dependent nonlinear acoustic system in subsonic mean flow is performed. The explicit predictor-corrector method of MacCormack (1969) is used. The methodology is tested on both uniform and sheared mean flows with plane and nonplanar sources. Results show that the acoustic system requires three physical boundary conditions on the inflow and one on the outflow boundary. The most natural choice for the inflow boundary conditions is judged to be a specification of the vorticity, the normal acoustic impedance, and a pressure gradient-density gradient relationship normal to the boundary. Specification of the acoustic pressure at the outflow boundary along with these inflow boundary conditions is found to give consistent reliable results. A set of boundary conditions developed earlier, which were intended to be nonreflecting is tested using the current method and is shown to yield unstable results for nonplanar acoustic waves.

On Analysis of Stationary Viscous Incompressible Flow Through a Radial Blade Machine

NASA Astrophysics Data System (ADS)

Neustupa, Tomáš

2010-09-01

The paper is concerned with the analysis of the two dimensional model of incompressible, viscous, stationary flow through a radial blade machine. This type of turbine is sometimes called Kaplan's turbine. In the technical area the use is either to force some regular characteristic to the flow of the medium going through the turbine (flow of melted iron, air conditioning) or to gain some energy from the flowing medium (water). The inflow and outflow part of boundary are in general a concentric circles. The larger one represents an inflow part of boundary the smaller one the outflow part of boundary. Between them are regularly spaced the blades of the machine. We study the existence of the weak solution in the case of nonlinear boundary condition of the "do-nothing" type. The model is interesting for study the behavior of the flow when the boundary is formed by mutually disjoint and separated parts.

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

NASA Technical Reports Server (NTRS)

Baysal, Oktay; Kaushik, Dinesh K.; Idres, Moumen

1997-01-01

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

An immersed boundary-simplified sphere function-based gas kinetic scheme for simulation of 3D incompressible flows

NASA Astrophysics Data System (ADS)

Yang, L. M.; Shu, C.; Yang, W. M.; Wang, Y.; Wu, J.

2017-08-01

In this work, an immersed boundary-simplified sphere function-based gas kinetic scheme (SGKS) is presented for the simulation of 3D incompressible flows with curved and moving boundaries. At first, the SGKS [Yang et al., "A three-dimensional explicit sphere function-based gas-kinetic flux solver for simulation of inviscid compressible flows," J. Comput. Phys. 295, 322 (2015) and Yang et al., "Development of discrete gas kinetic scheme for simulation of 3D viscous incompressible and compressible flows," J. Comput. Phys. 319, 129 (2016)], which is often applied for the simulation of compressible flows, is simplified to improve the computational efficiency for the simulation of incompressible flows. In the original SGKS, the integral domain along the spherical surface for computing conservative variables and numerical fluxes is usually not symmetric at the cell interface. This leads the expression of numerical fluxes at the cell interface to be relatively complicated. For incompressible flows, the sphere at the cell interface can be approximately considered to be symmetric as shown in this work. Besides that, the energy equation is usually not needed for the simulation of incompressible isothermal flows. With all these simplifications, the simple and explicit formulations for the conservative variables and numerical fluxes at the cell interface can be obtained. Second, to effectively implement the no-slip boundary condition for fluid flow problems with complex geometry as well as moving boundary, the implicit boundary condition-enforced immersed boundary method [Wu and Shu, "Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications," J. Comput. Phys. 228, 1963 (2009)] is introduced into the simplified SGKS. That is, the flow field is solved by the simplified SGKS without considering the presence of an immersed body and the no-slip boundary condition is implemented by the immersed boundary method. The accuracy and efficiency of the present scheme are validated by simulating the decaying vortex flow, flow past a stationary and rotating sphere, flow past a stationary torus, and flows over dragonfly flight.

An Improved Treatment of External Boundary for Three-Dimensional Flow Computations

NASA Technical Reports Server (NTRS)

Tsynkov, Semyon V.; Vatsa, Veer N.

1997-01-01

We present an innovative numerical approach for setting highly accurate nonlocal boundary conditions at the external computational boundaries when calculating three-dimensional compressible viscous flows over finite bodies. The approach is based on application of the difference potentials method by V. S. Ryaben'kii and extends our previous technique developed for the two-dimensional case. The new boundary conditions methodology has been successfully combined with the NASA-developed code TLNS3D and used for the analysis of wing-shaped configurations in subsonic and transonic flow regimes. As demonstrated by the computational experiments, the improved external boundary conditions allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable speedup of convergence of the multigrid iterations.

An outflow boundary condition for aeroacoustic computations

NASA Technical Reports Server (NTRS)

Hayder, M. Ehtesham; Hagstrom, Thomas

1995-01-01

A formulation of boundary condition for flows with small disturbances is presented. The authors test their methodology in an axisymmetric jet flow calculation, using both the Navier-Stokes and Euler equations. Solutions in the far field are assumed to be oscillatory. If the oscillatory disturbances are small, the growth of the solution variables can be predicted by linear theory. Eigenfunctions of the linear theory are used explicitly in the formulation of the boundary conditions. This guarantees correct solutions at the boundary in the limit where the predictions of linear theory are valid.

An experimental investigation of wall boundary layer transition Reynolds numbers in an expansion tube

NASA Technical Reports Server (NTRS)

Weilmuenster, K. J.

1974-01-01

Experimental measurements of boundary-layer transition in an expansion-tube test-gas flow are presented along with radial distributions of pitot pressure. An integral method for calculating constant Reynolds number lines for an expansion-tube flow is introduced. Comparison of experimental data and constant Reynolds number calculations has shown that for given conditions, wall boundary-layer transition occurs at a constant Reynolds number in an expansion-tube flow. Operating conditions in the expansion tube were chosen so that the effects of test-gas nonequilibrium on boundary-layer transition could be studied.

Time-dependent and outflow boundary conditions for Dissipative Particle Dynamics

PubMed Central

Lei, Huan; Fedosov, Dmitry A.; Karniadakis, George Em

2011-01-01

We propose a simple method to impose both no-slip boundary conditions at fluid-wall interfaces and at outflow boundaries in fully developed regions for Dissipative Particle Dynamics (DPD) fluid systems. The procedure to enforce the no-slip condition is based on a velocity-dependent shear force, which is a generalized force to represent the presence of the solid-wall particles and to maintain locally thermodynamic consistency. We show that this method can be implemented in both steady and time-dependent fluid systems and compare the DPD results with the continuum limit (Navier-Stokes) results. We also develop a force-adaptive method to impose the outflow boundary conditions for fully developed flow with unspecified outflow velocity profile or pressure value. We study flows over the backward-facing step and in idealized arterial bifurcations using a combination of the two new boundary methods with different flow rates. Finally, we explore the applicability of the outflow method in time-dependent flow systems. The outflow boundary method works well for systems with Womersley number of O(1), i.e., when the pressure and flowrate at the outflow are approximately in-phase. PMID:21499548

A quiet tunnel investigation of hypersonic boundary-layer stability over a cooled, flared cone

NASA Technical Reports Server (NTRS)

Blanchard, Alan E.; Selby, Gregory V.; Wilkinson, Stephen P.

1996-01-01

A flared-cone model under adiabatic and cooled-wall conditions was placed in a calibrated, low-disturbance Mach 6 flow and the stability of the boundary layer was investigated using a prototype constant-voltage anemometer. The results were compared with linear-stability theory predictions and good agreement was found in the prediction of second-mode frequencies and growth. In addition, the same 'N = 10' criterion used to predict boundary-layer transition in subsonic, transonic, and supersonic flows under low freestream noise conditions was found to be applicable for the hypersonic flow regime as well. Under cooled-wall conditions, a unique set of spectral data was acquired that documents the linear, nonlinear, and breakdown regions associated with the transition of hypersonic flow under low-noise conditions.

Development of a Flow Field for Testing a Boundary-Layer-Ingesting Propulsor

NASA Technical Reports Server (NTRS)

Hirt, Stefanie M.; Arend, David J.; Wolter, John D.

2017-01-01

The test section of the 8- by 6-Foot Supersonic Wind Tunnel at NASA Glenn Research Center was modified to produce the test conditions for a boundary-layer-ingesting propulsor. A test was conducted to measure the flow properties in the modified test section before the propulsor was installed. Measured boundary layer and freestream conditions were compared to results from computational fluid dynamics simulations of the external surface for the reference vehicle. Testing showed that the desired freestream conditions and boundary layer thickness could be achieved; however, some non-uniformity of the freestream conditions, particularly the total temperature, were observed.

Lattice Boltzmann method for simulating the viscous flow in large distensible blood vessels

NASA Astrophysics Data System (ADS)

Fang, Haiping; Wang, Zuowei; Lin, Zhifang; Liu, Muren

2002-05-01

A lattice Boltzmann method for simulating the viscous flow in large distensible blood vessels is presented by introducing a boundary condition for elastic and moving boundaries. The mass conservation for the boundary condition is tested in detail. The viscous flow in elastic vessels is simulated with a pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady flow agree with the analytical prediction to very high accuracy, and the simulation results for pulsatile flow are comparable with those of the aortic flows observed experimentally. The model is expected to find many applications for studying blood flows in large distensible arteries, especially in those suffering from atherosclerosis, stenosis, aneurysm, etc.

Choice of no-slip curved boundary condition for lattice Boltzmann simulations of high-Reynolds-number flows.

PubMed

Sanjeevi, Sathish K P; Zarghami, Ahad; Padding, Johan T

2018-04-01

Various curved no-slip boundary conditions available in literature improve the accuracy of lattice Boltzmann simulations compared to the traditional staircase approximation of curved geometries. Usually, the required unknown distribution functions emerging from the solid nodes are computed based on the known distribution functions using interpolation or extrapolation schemes. On using such curved boundary schemes, there will be mass loss or gain at each time step during the simulations, especially apparent at high Reynolds numbers, which is called mass leakage. Such an issue becomes severe in periodic flows, where the mass leakage accumulation would affect the computed flow fields over time. In this paper, we examine mass leakage of the most well-known curved boundary treatments for high-Reynolds-number flows. Apart from the existing schemes, we also test different forced mass conservation schemes and a constant density scheme. The capability of each scheme is investigated and, finally, recommendations for choosing a proper boundary condition scheme are given for stable and accurate simulations.

Choice of no-slip curved boundary condition for lattice Boltzmann simulations of high-Reynolds-number flows

NASA Astrophysics Data System (ADS)

Sanjeevi, Sathish K. P.; Zarghami, Ahad; Padding, Johan T.

2018-04-01

Various curved no-slip boundary conditions available in literature improve the accuracy of lattice Boltzmann simulations compared to the traditional staircase approximation of curved geometries. Usually, the required unknown distribution functions emerging from the solid nodes are computed based on the known distribution functions using interpolation or extrapolation schemes. On using such curved boundary schemes, there will be mass loss or gain at each time step during the simulations, especially apparent at high Reynolds numbers, which is called mass leakage. Such an issue becomes severe in periodic flows, where the mass leakage accumulation would affect the computed flow fields over time. In this paper, we examine mass leakage of the most well-known curved boundary treatments for high-Reynolds-number flows. Apart from the existing schemes, we also test different forced mass conservation schemes and a constant density scheme. The capability of each scheme is investigated and, finally, recommendations for choosing a proper boundary condition scheme are given for stable and accurate simulations.

A non-hydrostatic flat-bottom ocean model entirely based on Fourier expansion

NASA Astrophysics Data System (ADS)

Wirth, A.

2005-01-01

We show how to implement free-slip and no-slip boundary conditions in a three dimensional Boussinesq flat-bottom ocean model based on Fourier expansion. Our method is inspired by the immersed or virtual boundary technique in which the effect of boundaries on the flow field is modeled by a virtual force field. Our method, however, explicitly depletes the velocity on the boundary induced by the pressure, while at the same time respecting the incompressibility of the flow field. Spurious spatial oscillations remain at a negligible level in the simulated flow field when using our technique and no filtering of the flow field is necessary. We furthermore show that by using the method presented here the residual velocities at the boundaries are easily reduced to a negligible value. This stands in contradistinction to previous calculations using the immersed or virtual boundary technique. The efficiency is demonstrated by simulating a Rayleigh impulsive flow, for which the time evolution of the simulated flow is compared to an analytic solution, and a three dimensional Boussinesq simulation of ocean convection. The second instance is taken form a well studied oceanographic context: A free slip boundary condition is applied on the upper surface, the modeled sea surface, and a no-slip boundary condition to the lower boundary, the modeled ocean floor. Convergence properties of the method are investigated by solving a two dimensional stationary problem at different spatial resolutions. The work presented here is restricted to a flat ocean floor. Extensions of our method to ocean models with a realistic topography are discussed.

A zonal method for modeling powered-lift aircraft flow fields

NASA Technical Reports Server (NTRS)

Roberts, D. W.

1989-01-01

A zonal method for modeling powered-lift aircraft flow fields is based on the coupling of a three-dimensional Navier-Stokes code to a potential flow code. By minimizing the extent of the viscous Navier-Stokes zones the zonal method can be a cost effective flow analysis tool. The successful coupling of the zonal solutions provides the viscous/inviscid interations that are necessary to achieve convergent and unique overall solutions. The feasibility of coupling the two vastly different codes is demonstrated. The interzone boundaries were overlapped to facilitate the passing of boundary condition information between the codes. Routines were developed to extract the normal velocity boundary conditions for the potential flow zone from the viscous zone solution. Similarly, the velocity vector direction along with the total conditions were obtained from the potential flow solution to provide boundary conditions for the Navier-Stokes solution. Studies were conducted to determine the influence of the overlap of the interzone boundaries and the convergence of the zonal solutions on the convergence of the overall solution. The zonal method was applied to a jet impingement problem to model the suckdown effect that results from the entrainment of the inviscid zone flow by the viscous zone jet. The resultant potential flow solution created a lower pressure on the base of the vehicle which produces the suckdown load. The feasibility of the zonal method was demonstrated. By enhancing the Navier-Stokes code for powered-lift flow fields and optimizing the convergence of the coupled analysis a practical flow analysis tool will result.

Understanding the impact of insulating and conducting endplate boundary conditions on turbulence in CSDX through nonlocal simulations

DOE PAGES

Vaezi, P.; Holland, C.; Thakur, S. C.; ...

2017-04-01

The Controlled Shear Decorrelation Experiment (CSDX) linear plasma device provides a unique platform for investigating the underlying physics of self-regulating drift-wave turbulence/zonal flow dynamics. A minimal model of 3D drift-reduced nonlocal cold ion fluid equations which evolves density, vorticity, and electron temperature fluctuations, with proper sheath boundary conditions, is used to simulate dynamics of the turbulence in CSDX and its response to changes in parallel boundary conditions. These simulations are then carried out using the BOUndary Turbulence (BOUT++) framework and use equilibrium electron density and temperature profiles taken from experimental measurements. The results show that density gradient-driven drift-waves are themore » dominant instability in CSDX. However, the choice of insulating or conducting endplate boundary conditions affects the linear growth rates and energy balance of the system due to the absence or addition of Kelvin-Helmholtz modes generated by the sheath-driven equilibrium E × B shear and sheath-driven temperature gradient instability. Moreover, nonlinear simulation results show that the boundary conditions impact the turbulence structure and zonal flow formation, resulting in less broadband (more quasi-coherent) turbulence and weaker zonal flow in conducting boundary condition case. These results are qualitatively consistent with earlier experimental observations.« less

Immersed Boundary Models for Quantifying Flow-Induced Mechanical Stimuli on Stem Cells Seeded on 3D Scaffolds in Perfusion Bioreactors.

PubMed

Guyot, Yann; Smeets, Bart; Odenthal, Tim; Subramani, Ramesh; Luyten, Frank P; Ramon, Herman; Papantoniou, Ioannis; Geris, Liesbet

2016-09-01

Perfusion bioreactors regulate flow conditions in order to provide cells with oxygen, nutrients and flow-associated mechanical stimuli. Locally, these flow conditions can vary depending on the scaffold geometry, cellular confluency and amount of extra cellular matrix deposition. In this study, a novel application of the immersed boundary method was introduced in order to represent a detailed deformable cell attached to a 3D scaffold inside a perfusion bioreactor and exposed to microscopic flow. The immersed boundary model permits the prediction of mechanical effects of the local flow conditions on the cell. Incorporating stiffness values measured with atomic force microscopy and micro-flow boundary conditions obtained from computational fluid dynamics simulations on the entire scaffold, we compared cell deformation, cortical tension, normal and shear pressure between different cell shapes and locations. We observed a large effect of the precise cell location on the local shear stress and we predicted flow-induced cortical tensions in the order of 5 pN/μm, at the lower end of the range reported in literature. The proposed method provides an interesting tool to study perfusion bioreactors processes down to the level of the individual cell's micro-environment, which can further aid in the achievement of robust bioprocess control for regenerative medicine applications.

Boundary Condition Study for the Juncture Flow Experiment in the NASA Langley 14x22-Foot Subsonic Wind Tunnel

NASA Technical Reports Server (NTRS)

Rumsey, C. L.; Carlson, J.-R.; Hannon, J. A.; Jenkins, L. N.; Bartram, S. M.; Pulliam, T. H.; Lee, H. C.

2017-01-01

Because future wind tunnel tests associated with the NASA Juncture Flow project are being designed for the purpose of CFD validation, considerable effort is going into the characterization of the wind tunnel boundary conditions, particularly at inflow. This is important not only because wind tunnel flowfield nonuniformities can play a role in integrated testing uncertainties, but also because the better the boundary conditions are known, the better CFD can accurately represent the experiment. This paper describes recent investigative wind tunnel tests involving two methods to measure and characterize the oncoming flow in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The features of each method, as well as some of their pros and cons, are highlighted. Boundary conditions and modeling tactics currently used by CFD for empty-tunnel simulations are also described, and some results using three different CFD codes are shown. Preliminary CFD parametric studies associated with the Juncture Flow model are summarized, to determine sensitivities of the flow near the wing-body juncture region of the model to a variety of modeling decisions.

New boundary conditions for fluid interaction with hydrophobic surface

NASA Astrophysics Data System (ADS)

Pochylý, František; Fialová, Simona; Havlásek, Michal

2018-06-01

Solution of both laminar and turbulent flow with consideration of hydrophobic surface is based on the original Navier assumption that the shear stress on the hydrophobic surface is directly proportional to the slipping velocity. In the previous work a laminar flow analysis with different boundary conditions was performed. The shear stress value on the tube walls directly depends on the pressure gradient. In the solution of the turbulent flow by the k-ɛ model, the occurrence of the fluctuation components of velocity on the hydrophobic surface is considered. The fluctuation components of the velocity affect the size of the adhesive forces. We assume that the boundary condition for ɛ depending on the velocity gradients will not need to be changed. When the liquid slips over the surface, non-zero fluctuation velocity components occur in the turbulent flow. These determine the non-zero value of the turbulent kinetic energy K. In addition, the fluctuation velocity components also influence the value of the adhesive forces, so it is necessary to include these in the formulation of new boundary conditions for turbulent flow on the hydrophobic surface.

Elasto visco-plastic flow with special attention to boundary conditions

NASA Technical Reports Server (NTRS)

Shimazaki, Y.; Thompson, E. G.

1981-01-01

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Assignment of boundary conditions in embedded ground water flow models

USGS Publications Warehouse

Leake, S.A.

1998-01-01

Many small-scale ground water models are too small to incorporate distant aquifer boundaries. If a larger-scale model exists for the area of interest, flow and head values can be specified for boundaries in the smaller-scale model using values from the larger-scale model. Flow components along rows and columns of a large-scale block-centered finite-difference model can be interpolated to compute horizontal flow across any segment of a perimeter of a small-scale model. Head at cell centers of the larger-scale model can be interpolated to compute head at points on a model perimeter. Simple linear interpolation is proposed for horizontal interpolation of horizontal-flow components. Bilinear interpolation is proposed for horizontal interpolation of head values. The methods of interpolation provided satisfactory boundary conditions in tests using models of hypothetical aquifers.Many small-scale ground water models are too small to incorporate distant aquifer boundaries. If a larger-scale model exists for the area of interest, flow and head values can be specified for boundaries in the smaller-scale model using values from the larger-scale model. Flow components along rows and columns of a large-scale block-centered finite-difference model can be interpolated to compute horizontal flow across any segment of a perimeter of a small-scale model. Head at cell centers of the larger.scale model can be interpolated to compute head at points on a model perimeter. Simple linear interpolation is proposed for horizontal interpolation of horizontal-flow components. Bilinear interpolation is proposed for horizontal interpolation of head values. The methods of interpolation provided satisfactory boundary conditions in tests using models of hypothetical aquifers.

Laminar supersonic flow over a backstep - A numerical solution at higher Reynolds numbers

NASA Technical Reports Server (NTRS)

Kronzon, Y.; Rom, J.; Seginer, A.

1976-01-01

The Allen-Cheng solution of the flow over a backward facing step is extended to Reynolds numbers up to 16,000 and to inflow boundary-layer height ratios as low as 0.1 by moving the downstream boundary into the recompression region and by smoothing the resulting errors. The boundary conditions in the supersonic outer flow and the downstream boundary conditions in the wake are determined by an extrapolation procedure. Computational results are compared with relevant experimental data. Fair agreement is found between the calculated base pressures and the experimental values, whereas agreement between heat transfer rates appears to be qualitative only.

Numerical analysis of thermal creep flow in curved channels for designing a prototype of Knudsen micropump

NASA Astrophysics Data System (ADS)

Leontidis, V.; Brandner, J. J.; Baldas, L.; Colin, S.

2012-05-01

The possibility to generate a gas flow inside a channel just by imposing a tangential temperature gradient along the walls without the existence of an initial pressure difference is well known. The gas must be under rarefied conditions, meaning that the system must operate between the slip and the free molecular flow regimes, either at low pressure or/and at micro/nano-scale dimensions. This phenomenon is at the basis of the operation principle of Knudsen pumps, which are actually compressors without any moving parts. Nowadays, gas flows in the slip flow regime through microchannels can be modeled using commercial Computational Fluid Dynamics softwares, because in this regime the compressible Navier-Stokes equations with appropriate boundary conditions are still valid. A simulation procedure has been developed for the modeling of thermal creep flow using ANSYS Fluent®. The implementation of the boundary conditions is achieved by developing User Defined Functions (UDFs) by means of C++ routines. The complete first order velocity slip boundary condition, including the thermal creep effects due to the axial temperature gradient and the effect of the wall curvature, and the temperature jump boundary condition are applied. The developed simulation tool is used for the preliminary design of Knudsen micropumps consisting of a sequence of curved and straight channels.

Cylindrical Couette flow of a rarefied gas: Effect of a boundary condition on the inverted velocity profile.

PubMed

Kosuge, Shingo

2015-07-01

The cylindrical Couette flow of a rarefied gas between a rotating inner cylinder and a stationary outer cylinder is investigated under the following two kinds of kinetic boundary conditions. One is the modified Maxwell-type boundary condition proposed by Dadzie and Méolans [J. Math. Phys. 45, 1804 (2004)] and the other is the Cercignani-Lampis condition, both of which have separate accommodation coefficients associated with the molecular velocity component normal to the boundary and with the tangential component. An asymptotic analysis of the Boltzmann equation for small Knudsen numbers and a numerical analysis of the Bhatnagar-Gross-Krook model equation for a wide range of the Knudsen number are performed to clarify the effect of each accommodation coefficient as well as of the boundary condition itself on the behavior of the gas, especially on the flow-velocity profile. As a result, the velocity-slip and temperature-jump conditions corresponding to the above kinetic boundary conditions are derived, which are necessary for the fluid-dynamic description of the problem for small Knudsen numbers. The parameter range for the onset of the velocity inversion phenomenon, which is related mainly to the decrease in the tangential momentum accommodation, is also obtained.

Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime: Application to plane boundaries

NASA Astrophysics Data System (ADS)

Silva, Goncalo; Semiao, Viriato

2017-07-01

The first nonequilibrium effect experienced by gaseous flows in contact with solid surfaces is the slip-flow regime. While the classical hydrodynamic description holds valid in bulk, at boundaries the fluid-wall interactions must consider slip. In comparison to the standard no-slip Dirichlet condition, the case of slip formulates as a Robin-type condition for the fluid tangential velocity. This makes its numerical modeling a challenging task, particularly in complex geometries. In this work, this issue is handled with the lattice Boltzmann method (LBM), motivated by the similarities between the closure relations of the reflection-type boundary schemes equipping the LBM equation and the slip velocity condition established by slip-flow theory. Based on this analogy, we derive, as central result, the structure of the LBM boundary closure relation that is consistent with the second-order slip velocity condition, applicable to planar walls. Subsequently, three tasks are performed. First, we clarify the limitations of existing slip velocity LBM schemes, based on discrete analogs of kinetic theory fluid-wall interaction models. Second, we present improved slip velocity LBM boundary schemes, constructed directly at discrete level, by extending the multireflection framework to the slip-flow regime. Here, two classes of slip velocity LBM boundary schemes are considered: (i) linear slip schemes, which are local but retain some calibration requirements and/or operation limitations, (ii) parabolic slip schemes, which use a two-point implementation but guarantee the consistent prescription of the intended slip velocity condition, at arbitrary plane wall discretizations, further dispensing any numerical calibration procedure. Third and final, we verify the improvements of our proposed slip velocity LBM boundary schemes against existing ones. The numerical tests evaluate the ability of the slip schemes to exactly accommodate the steady Poiseuille channel flow solution, over distinct wall slippage conditions, namely, no-slip, first-order slip, and second-order slip. The modeling of channel walls is discussed at both lattice-aligned and non-mesh-aligned configurations: the first case illustrates the numerical slip due to the incorrect modeling of slippage coefficients, whereas the second case adds the effect of spurious boundary layers created by the deficient accommodation of bulk solution. Finally, the slip-flow solutions predicted by LBM schemes are further evaluated for the Knudsen's paradox problem. As conclusion, this work establishes the parabolic accuracy of slip velocity schemes as the necessary condition for the consistent LBM modeling of the slip-flow regime.

Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime: Application to plane boundaries.

PubMed

Silva, Goncalo; Semiao, Viriato

2017-07-01

The first nonequilibrium effect experienced by gaseous flows in contact with solid surfaces is the slip-flow regime. While the classical hydrodynamic description holds valid in bulk, at boundaries the fluid-wall interactions must consider slip. In comparison to the standard no-slip Dirichlet condition, the case of slip formulates as a Robin-type condition for the fluid tangential velocity. This makes its numerical modeling a challenging task, particularly in complex geometries. In this work, this issue is handled with the lattice Boltzmann method (LBM), motivated by the similarities between the closure relations of the reflection-type boundary schemes equipping the LBM equation and the slip velocity condition established by slip-flow theory. Based on this analogy, we derive, as central result, the structure of the LBM boundary closure relation that is consistent with the second-order slip velocity condition, applicable to planar walls. Subsequently, three tasks are performed. First, we clarify the limitations of existing slip velocity LBM schemes, based on discrete analogs of kinetic theory fluid-wall interaction models. Second, we present improved slip velocity LBM boundary schemes, constructed directly at discrete level, by extending the multireflection framework to the slip-flow regime. Here, two classes of slip velocity LBM boundary schemes are considered: (i) linear slip schemes, which are local but retain some calibration requirements and/or operation limitations, (ii) parabolic slip schemes, which use a two-point implementation but guarantee the consistent prescription of the intended slip velocity condition, at arbitrary plane wall discretizations, further dispensing any numerical calibration procedure. Third and final, we verify the improvements of our proposed slip velocity LBM boundary schemes against existing ones. The numerical tests evaluate the ability of the slip schemes to exactly accommodate the steady Poiseuille channel flow solution, over distinct wall slippage conditions, namely, no-slip, first-order slip, and second-order slip. The modeling of channel walls is discussed at both lattice-aligned and non-mesh-aligned configurations: the first case illustrates the numerical slip due to the incorrect modeling of slippage coefficients, whereas the second case adds the effect of spurious boundary layers created by the deficient accommodation of bulk solution. Finally, the slip-flow solutions predicted by LBM schemes are further evaluated for the Knudsen's paradox problem. As conclusion, this work establishes the parabolic accuracy of slip velocity schemes as the necessary condition for the consistent LBM modeling of the slip-flow regime.

Rarefaction effects in microchannel gas flow driven by rhythmic wall contractions

NASA Astrophysics Data System (ADS)

Chatterjee, Krishnashis; Staples, Anne; Department of Biomedical Engineering; Mechanics, Virginia Tech Collaboration

2015-11-01

Current state of the art microfluidic devices employ precise and timely operation of a complex arrangement of micropumps and valves for fluid transport. A much more novel flow transport mechanism is found in entomological respiratory systems, which involve rhythmic wall contractions for driving the fluid flow. The practical viability of using this technique in future microfluidic devices has been studied earlier. The present study investigates the incorporation of rarefaction effects in the above model of microscale gas flow by including slip boundary conditions. The Navier Stokes equations for gas flow in rectangular microchannel are solved analytically with microscale and lubrication theory assumptions. First order slip boundary conditions are incorporated to account for the rarefaction effects. The dependence of fluid velocities and pressure gradient on the slip boundary conditions is studied. Time averaged unidirectional fluid flow rates are plotted for different phase lags between the contractions, with and without slip in order to obtain an optimum range under different conditions.

Experimental evaluation of heat transfer on a 1030:1 area ratio rocket nozzle

NASA Technical Reports Server (NTRS)

Kacynski, Kenneth J.; Pavli, Albert J.; Smith, Tamara A.

1987-01-01

A 1030:1 carbon steel, heat-sink nozzle was tested. The test conditions included a nominal chamber pressure of 2413 kN/sq m and a mixture ratio range of 2.78 to 5.49. The propellants were gaseous oxygen and gaseous hydrogen. Outer wall temperature measurements were used to calculate the inner wall temperature and the heat flux and heat rate to the nozzle at specified axial locations. The experimental heat fluxes were compared to those predicted by the Two-Dimensional Kinetics (TDK) computer model analysis program. When laminar boundary layer flow was assumed in the analysis, the predicted values were within 15 percent of the experimental values for the area ratios of 20 to 975. However, when turbulent boundary layer conditions were assumed, the predicted values were approximately 120 percent higher than the experimental values. A study was performed to determine if the conditions within the nozzle could sustain a laminar boundary layer. Using the flow properties predicted by TDK, the momentum-thickness Reynolds number was calculated, and the point of transition to turbulent flow was predicted. The predicted transition point was within 0.5 inches of the nozzle throat. Calculations of the acceleration parameter were then made to determine if the flow conditions could produce relaminarization of the boundary layer. It was determined that if the boundary layer flow was inclined to transition to turbulent, the acceleration conditions within the nozzle would tend to suppress turbulence and keep the flow laminar-like.

Experimental evaluation of heat transfer on a 1030:1 area ratio rocket nozzle

NASA Technical Reports Server (NTRS)

Kacynski, Kenneth J.; Pavli, Albert J.; Smith, Tamara A.

1987-01-01

A 1030:1 carbon steel, heat-sink nozzle was tested. The test conditions included a nominal chamber pressure of 2413 kN/sq m and a mixture ratio range of 2.78 to 5.49. The propellants were gaseous oxygen and gaseous hydrogen. Outer wall temperature measurements were used to calculate the inner wall temperature and the heat flux and heat rate to the nozzle at specified axial locations. The experimental heat fluxes were compared to those predicted by the Two-Dimensional Kinetics (TDK) computer model analysis program. When laminar boundary layer flow was assumed in the analysis, the predicted values were within 15% of the experimental values for the area ratios of 20 to 975. However, when turbulent boundary layer conditions were assumed, the predicted values were approximately 120% higher than the experimental values. A study was performed to determine if the conditions within the nozzle could sustain a laminar boundary layer. Using the flow properties predicted by TDK, the momentum-thickness Reynolds number was calculated, and the point of transition to turbulent flow was predicted. The predicted transition point was within 0.5 inches of the nozzle throat. Calculations of the acceleration parameter were then made to determine if the flow conditions could produce relaminarization of the boundary layer. It was determined that if the boundary layer flow was inclined to transition to turbulent, the acceleration conditions within the nozzle would tend to suppress turbulence and keep the flow laminar-like.

High order accurate solutions of viscous problems

NASA Technical Reports Server (NTRS)

Hayder, M. Ehtesham; Turkel, Eli

1993-01-01

We consider a fourth order extension to MacCormack's scheme. The original extension was fourth order only for the inviscid terms but was second order for the viscous terms. We show how to modify the viscous terms so that the scheme is uniformly fourth order in the spatial derivatives. Applications are given to some boundary layer flows. In addition, for applications to shear flows the effect of the outflow boundary conditions are very important. We compare the accuracy of several of these different boundary conditions for both boundary layer and shear flows. Stretching at the outflow usually increases the oscillations in the numerical solution but the addition of a filtered sponge layer (with or without stretching) reduces such oscillations. The oscillations are generated by insufficient resolution of the shear layer. When the shear layer is sufficiently resolved then oscillations are not generated and there is less of a need for a nonreflecting boundary condition.

Global Artificial Boundary Conditions for Computation of External Flow Problems with Propulsive Jets

NASA Technical Reports Server (NTRS)

Tsynkov, Semyon; Abarbanel, Saul; Nordstrom, Jan; Ryabenkii, Viktor; Vatsa, Veer

1998-01-01

We propose new global artificial boundary conditions (ABC's) for computation of flows with propulsive jets. The algorithm is based on application of the difference potentials method (DPM). Previously, similar boundary conditions have been implemented for calculation of external compressible viscous flows around finite bodies. The proposed modification substantially extends the applicability range of the DPM-based algorithm. In the paper, we present the general formulation of the problem, describe our numerical methodology, and discuss the corresponding computational results. The particular configuration that we analyze is a slender three-dimensional body with boat-tail geometry and supersonic jet exhaust in a subsonic external flow under zero angle of attack. Similarly to the results obtained earlier for the flows around airfoils and wings, current results for the jet flow case corroborate the superiority of the DPM-based ABC's over standard local methodologies from the standpoints of accuracy, overall numerical performance, and robustness.

The internal boundary layer — A review

NASA Astrophysics Data System (ADS)

Garratt, J. R.

1990-03-01

A review is given of relevant work on the internal boundary layer (IBL) associated with: (i) Small-scale flow in neutral conditions across an abrupt change in surface roughness, (ii) Small-scale flow in non-neutral conditions across an abrupt change in surface roughness, temperature or heat/moisture flux, (iii) Mesoscale flow, with emphasis on flow across the coastline for both convective and stably stratified conditions. The major theme in all cases is on the downstream, modified profile form (wind and temperature), and on the growth relations for IBL depth.

A comparative study of various inflow boundary conditions and turbulence models for wind turbine wake predictions

NASA Astrophysics Data System (ADS)

Tian, Lin-Lin; Zhao, Ning; Song, Yi-Lei; Zhu, Chun-Ling

2018-05-01

This work is devoted to perform systematic sensitivity analysis of different turbulence models and various inflow boundary conditions in predicting the wake flow behind a horizontal axis wind turbine represented by an actuator disc (AD). The tested turbulence models are the standard k-𝜀 model and the Reynolds Stress Model (RSM). A single wind turbine immersed in both uniform flows and in modeled atmospheric boundary layer (ABL) flows is studied. Simulation results are validated against the field experimental data in terms of wake velocity and turbulence intensity.

A critical assessment of viscous models of trench topography and corner flow

NASA Technical Reports Server (NTRS)

Zhang, J.; Hager, B. H.; Raefsky, A.

1984-01-01

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vaezi, P.; Holland, C.; Thakur, S. C.

The Controlled Shear Decorrelation Experiment (CSDX) linear plasma device provides a unique platform for investigating the underlying physics of self-regulating drift-wave turbulence/zonal flow dynamics. A minimal model of 3D drift-reduced nonlocal cold ion fluid equations which evolves density, vorticity, and electron temperature fluctuations, with proper sheath boundary conditions, is used to simulate dynamics of the turbulence in CSDX and its response to changes in parallel boundary conditions. These simulations are then carried out using the BOUndary Turbulence (BOUT++) framework and use equilibrium electron density and temperature profiles taken from experimental measurements. The results show that density gradient-driven drift-waves are themore » dominant instability in CSDX. However, the choice of insulating or conducting endplate boundary conditions affects the linear growth rates and energy balance of the system due to the absence or addition of Kelvin-Helmholtz modes generated by the sheath-driven equilibrium E × B shear and sheath-driven temperature gradient instability. Moreover, nonlinear simulation results show that the boundary conditions impact the turbulence structure and zonal flow formation, resulting in less broadband (more quasi-coherent) turbulence and weaker zonal flow in conducting boundary condition case. These results are qualitatively consistent with earlier experimental observations.« less

A Comparative Experimental Study of Fixed Temperature and Fixed Heat Flux Boundary Conditions in Turbulent Thermal Convection

NASA Astrophysics Data System (ADS)

Huang, Shi-Di; Wang, Fei; Xi, Heng-Dong; Xia, Ke-Qing

2014-11-01

We report an experimental study of the influences of thermal boundary condition in turbulent thermal convection. Two configurations were examined: one was fixed heat flux at the bottom boundary and fixed temperature at the top (HC cells); the other was fixed temperature at both boundaries (CC cells). It is found that the flow strength in the CC cells is on average 9% larger than that in the HC ones, which could be understood as change in plume emission ability under different boundary conditions. It is further found, rather surprisingly, that flow reversals of the large-scale circulation occur more frequently in the CC cell, despite a stronger large-scale flow and more uniform temperature distribution over the boundaries. These findings provide new insights into turbulent thermal convection and should stimulate further studies, especially experimental ones. This work is supported by the Hong Kong Research Grants Council under Grant No. CUHK 403712.

The respective roles of bulk friction and slip velocity during a granular mass flow

NASA Astrophysics Data System (ADS)

Staron, Lydie

2016-04-01

Catastrophic granular mass flows form an important natural hazard. Mitigation has motivated numerous studies on the properties of natural granular flows, and in particular, their ability to travel long distances away from the release point. The mobility of granular flows is commonly characterised through the definition of rheological properties and effective friction. Yet, it is widely accepted that the description in term of effective friction may include various lubrication effects, softening at the base of the flow and large slip velocities being a most likely one. In this case, flow bulk properties may obliterate the flow boundary conditions. In this contribution, we investigate how disentangling bulk properties from boundary conditions may improve our understanding of the flow. Using discrete simulations, we induce increasing slip velocities in different flow configurations. We show that increased mobility may be achieved without changing bulk properties. The results are interpreted in terms of a Robin-Navier slip condition and implemented in a continuum Navier-Stokes solver. We quantify the respective role of rheological bulk properties and boundary conditions in the general behaviour of a transient mass flow. We show that omitting the description of boundary conditions leads to misinterpretation of the flow properties. The outcome is discussed in terms of models reliability. References P.-Y. Lagrée et al, The granular column collapse as a continuum: validity of a two-dimensional Navier-Stokes model with the mu(I) rheology, J. Fluid Mech. 686, 378-408 (2011) L. Staron and E. Lajeunesse, Understanding how the volume affects the mobility of dry debris flows, Geophys. Res. Lett. 36, L12402 (2009) L. Staron, Mobility of long-runout rock flows: a discrete numerical investigation, Geophys. J. Int. 172, 455-463 (2008)

Boundary layer effects on liners for aircraft engines

NASA Astrophysics Data System (ADS)

Gabard, Gwénaël

2016-10-01

The performance of acoustic treatments installed on aircraft engines is strongly influenced by the boundary layer of the grazing flow on the surface of the liner. The parametric study presented in this paper illustrates the extent of this effect and identifies when it is significant. The acoustic modes of a circular duct with flow are calculated using a finite difference method. The parameters are representative of the flow conditions, liners and sound fields found in current turbofan engines. Both the intake and bypass ducts are considered. Results show that there is a complex interplay between the boundary layer thickness, the direction of propagation and the liner impedance and that the boundary layer can have a strong impact on liner performance for typical configurations (including changes of the order of 30 dB on the attenuation of modes associated with tonal fan noise). A modified impedance condition including the effect of a small but finite boundary layer thickness is considered and compared to the standard Myers condition based on an infinitely thin boundary layer. We show how this impedance condition can be implemented in a mode calculation method by introducing auxiliary variables. This condition is able to capture the trends associated with the boundary layer effects and in most cases provides improved predictions of liner performance.

Computation of the shock-wave boundary layer interaction with flow separation

NASA Technical Reports Server (NTRS)

Ardonceau, P.; Alziary, T.; Aymer, D.

1980-01-01

The boundary layer concept is used to describe the flow near the wall. The external flow is approximated by a pressure displacement relationship (tangent wedge in linearized supersonic flow). The boundary layer equations are solved in finite difference form and the question of the presence and unicity of the solution is considered for the direct problem (assumed pressure) or converse problem (assumed displacement thickness, friction ratio). The coupling algorithm presented implicitly processes the downstream boundary condition necessary to correctly define the interacting boundary layer problem. The algorithm uses a Newton linearization technique to provide a fast convergence.

Efficient Fluid Dynamic Design Optimization Using Cartesian Grids

NASA Technical Reports Server (NTRS)

Dadone, A.; Grossman, B.; Sellers, Bill (Technical Monitor)

2004-01-01

This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.

Calculation of Multistage Turbomachinery Using Steady Characteristic Boundary Conditions

NASA Technical Reports Server (NTRS)

Chima, Rodrick V.

1998-01-01

A multiblock Navier-Stokes analysis code for turbomachinery has been modified to allow analysis of multistage turbomachines. A steady averaging-plane approach was used to pass information between blade rows. Characteristic boundary conditions written in terms of perturbations about the mean flow from the neighboring blade row were used to allow close spacing between the blade rows without forcing the flow to be axisymmetric. In this report the multiblock code is described briefly and the characteristic boundary conditions and the averaging-plane implementation are described in detail. Two approaches for averaging the flow properties are also described. A two-dimensional turbine stator case was used to compare the characteristic boundary conditions with standard axisymmetric boundary conditions. Differences were apparent but small in this low-speed case. The two-stage fuel turbine used on the space shuttle main engines was then analyzed using a three-dimensional averaging-plane approach. Computed surface pressure distributions on the stator blades and endwalls and computed distributions of blade surface heat transfer coefficient on three blades showed very good agreement with experimental data from two tests.

A physical approach to the numerical treatment of boundaries in gas dynamics

NASA Technical Reports Server (NTRS)

Moretti, G.

1981-01-01

Two types of boundaries are considered: rigid walls, and artificial (open) boundaries which were arbitrarily drawn somewhere across a wider flow field. A set of partial differential equations (typically, the Euler equations) has an infinite number of solutions, each one defined by a set of initial and boundary conditions. The initial conditions remaining the same, any change in the boundary conditions will produce a new solution. To pose the problem well, a necessary and sufficient number of boundary conditions are prescribed.

User's guide to the Variably Saturated Flow (VSF) process to MODFLOW

USGS Publications Warehouse

Thoms, R. Brad; Johnson, Richard L.; Healy, Richard W.

2006-01-01

A new process for simulating three-dimensional (3-D) variably saturated flow (VSF) using Richards' equation has been added to the 3-D modular finite-difference ground-water model MODFLOW. Five new packages are presented here as part of the VSF Process--the Richards' Equation Flow (REF1) Package, the Seepage Face (SPF1) Package, the Surface Ponding (PND1) Package, the Surface Evaporation (SEV1) Package, and the Root Zone Evapotranspiration (RZE1) Package. Additionally, a new Adaptive Time-Stepping (ATS1) Package is presented for use by both the Ground-Water Flow (GWF) Process and VSF. The VSF Process allows simulation of flow in unsaturated media above the ground-water zone and facilitates modeling of ground-water/surface-water interactions. Model performance is evaluated by comparison to an analytical solution for one-dimensional (1-D) constant-head infiltration (Dirichlet boundary condition), field experimental data for a 1-D constant-head infiltration, laboratory experimental data for two-dimensional (2-D) constant-flux infiltration (Neumann boundary condition), laboratory experimental data for 2-D transient drainage through a seepage face, and numerical model results (VS2DT) of a 2-D flow-path simulation using realistic surface boundary conditions. A hypothetical 3-D example case also is presented to demonstrate the new capability using periodic boundary conditions (for example, daily precipitation) and varied surface topography over a larger spatial scale (0.133 square kilometer). The new model capabilities retain the modular structure of the MODFLOW code and preserve MODFLOW's existing capabilities as well as compatibility with commercial pre-/post-processors. The overall success of the VSF Process in simulating mixed boundary conditions and variable soil types demonstrates its utility for future hydrologic investigations. This report presents a new flow package implementing the governing equations for variably saturated ground-water flow, four new boundary condition packages unique to unsaturated flow, the Adaptive Time-Stepping Package for use with both the GWF Process and the new VSF Process, detailed descriptions of the input and output files for each package, and six simulation examples verifying model performance.

Laser Doppler velocimeter measurements of boundary layer velocity and turbulent intensities in Mach 2.5 flow

NASA Technical Reports Server (NTRS)

Sewell, Jesse; Chew, Larry

1994-01-01

In recent years, the interest in developing a high-speed civil transport has increased. This has led to an increase in research activity on compressible supersonic flows, in particular the boundary layer. The structure of subsonic boundary layers has been extensively documented using conditional sampling techniques which exploit the knowledge of both u and v velocities. Researchers using these techniques have been able to explore some of the complex three-dimensional motions which are responsible for Reynolds stress production and transport in the boundary layer. As interest in turbulent structure has grown to include supersonic flows, a need for simultaneous multicomponent velocity measurements in these flows has developed. The success of conditional analysis in determining the characteristics of coherent motions and structures in the boundary layer relies on accurate, simultaneous measurement of two instantaneous velocity components.

Unsteady-flow-field predictions for oscillating cascades

NASA Technical Reports Server (NTRS)

Huff, Dennis L.

1991-01-01

The unsteady flow field around an oscillating cascade of flat plates with zero stagger was studied by using a time marching Euler code. This case had an exact solution based on linear theory and served as a model problem for studying pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step size was shown for a moderate reduced frequency. Results show that an approximate nonreflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer to the airfoils than when reflective boundaries are used. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates captures the unsteady flow field better than when uniform grids are used as long as the 'Courant Friedrichs Levy' (CFL) number is less than 1 for a sufficient portion of the grid. Finally, a solution based on an optimization of grid, CFL number, and boundary conditions shows good agreement with linear theory.

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

NASA Technical Reports Server (NTRS)

Hu, Fang; Pizzo, Michelle E.; Nark, Douglas M.

2017-01-01

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Well-posedness of the free boundary problem in compressible elastodynamics

NASA Astrophysics Data System (ADS)

Trakhinin, Yuri

2018-02-01

We study the free boundary problem for the flow of a compressible isentropic inviscid elastic fluid. At the free boundary moving with the velocity of the fluid particles the columns of the deformation gradient are tangent to the boundary and the pressure vanishes outside the flow domain. We prove the local-in-time existence of a unique smooth solution of the free boundary problem provided that among three columns of the deformation gradient there are two which are non-collinear vectors at each point of the initial free boundary. If this non-collinearity condition fails, the local-in-time existence is proved under the classical Rayleigh-Taylor sign condition satisfied at the first moment. By constructing an Hadamard-type ill-posedness example for the frozen coefficients linearized problem we show that the simultaneous failure of the non-collinearity condition and the Rayleigh-Taylor sign condition leads to Rayleigh-Taylor instability.

Flow boundary conditions for chain-end adsorbing polymer blends.

PubMed

Zhou, Xin; Andrienko, Denis; Delle Site, Luigi; Kremer, Kurt

2005-09-08

Using the phenol-terminated polycarbonate blend as an example, we demonstrate that the hydrodynamic boundary conditions for a flow of an adsorbing polymer melt are extremely sensitive to the structure of the epitaxial layer. Under shear, the adsorbed parts (chain ends) of the polymer melt move along the equipotential lines of the surface potential whereas the adsorbed additives serve as the surface defects. In response to the increase of the number of the adsorbed additives the surface layer becomes thinner and solidifies. This results in a gradual transition from the slip to the no-slip boundary condition for the melt flow, with a nonmonotonic dependence of the slip length on the surface concentration of the adsorbed ends.

Helmholtz decomposition revisited: Vorticity generation and trailing edge condition. I - Incompressible flows

NASA Technical Reports Server (NTRS)

Morino, L.

1986-01-01

Using the decomposition for the infinite-space, the issue of the nonuniqueness of the Helmholtz decomposition for the problem of the three-dimensional unsteady incompressible flow around a body is considered. A representation for the velocity that is valid for both the fluid region and the region inside the boundary surface is employed, and the motion of the boundary is described as the limiting case of a sequence of impulsive accelerations. At each instant of velocity discontinuity, vorticity is shown to be generated by the boundary condition on the normal component of the velocity, for both inviscid and viscous flows. In viscous flows, the vorticity is shown to diffuse into the surroundings, and the no-slip conditions are automatically satisfied. A trailing edge condition must be satisfied for the solution to the Euler equations to be the limit of the solution of the Navier-Stokes equations.

Application of a low order panel method to complex three-dimensional internal flow problems

NASA Technical Reports Server (NTRS)

Ashby, D. L.; Sandlin, D. R.

1986-01-01

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

Near Continuum Velocity and Temperature Coupled Compressible Boundary Layer Flow over a Flat Plate

NASA Astrophysics Data System (ADS)

He, Xin; Cai, Chunpei

2017-04-01

The problem of a compressible gas flows over a flat plate with the velocity-slip and temperature-jump boundary conditions are being studied. The standard single- shooting method is applied to obtain the exact solutions for velocity and temperature profiles when the momentum and energy equations are weakly coupled. A double-shooting method is applied if these two equations are closely coupled. If the temperature affects the velocity directly, more significant velocity slip happens at locations closer to the plate's leading edge, and inflections on the velocity profiles appear, indicating flows may become unstable. As a consequence, the temperature-jump and velocity-slip boundary conditions may trigger earlier flow transitions from a laminar to a turbulent flow state.

Open-ocean boundary conditions from interior data: Local and remote forcing of Massachusetts Bay

USGS Publications Warehouse

Bogden, P.S.; Malanotte-Rizzoli, P.; Signell, R.

1996-01-01

Massachusetts and Cape Cod Bays form a semienclosed coastal basin that opens onto the much larger Gulf of Maine. Subtidal circulation in the bay is driven by local winds and remotely driven flows from the gulf. The local-wind forced flow is estimated with a regional shallow water model driven by wind measurements. The model uses a gravity wave radiation condition along the open-ocean boundary. Results compare reasonably well with observed currents near the coast. In some offshore regions however, modeled flows are an order of magnitude less energetic than the data. Strong flows are observed even during periods of weak local wind forcing. Poor model-data comparisons are attributable, at least in part, to open-ocean boundary conditions that neglect the effects of remote forcing. Velocity measurements from within Massachusetts Bay are used to estimate the remotely forced component of the flow. The data are combined with shallow water dynamics in an inverse-model formulation that follows the theory of Bennett and McIntosh [1982], who considered tides. We extend their analysis to consider the subtidal response to transient forcing. The inverse model adjusts the a priori open-ocean boundary condition, thereby minimizing a combined measure of model-data misfit and boundary condition adjustment. A "consistency criterion" determines the optimal trade-off between the two. The criterion is based on a measure of plausibility for the inverse solution. The "consistent" inverse solution reproduces 56% of the average squared variation in the data. The local-wind-driven flow alone accounts for half of the model skill. The other half is attributable to remotely forced flows from the Gulf of Maine. The unexplained 44% comes from measurement errors and model errors that are not accounted for in the analysis. 

Numerical and Experimental Studies of the Natural Convection Flow Within a Horizontal Cylinder Subjected to a Uniformly Cold Wall Boundary Condition. Ph.D. Thesis - Va. Poly. Inst. and State Univ.

NASA Technical Reports Server (NTRS)

Stewart, R. B.

1972-01-01

Numberical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time dependent natural convection flow within a horizontal cylinder. The early time flow development and wall heat transfer is obtained after imposing a uniformly cold wall boundary condition on the cylinder. Solutions are also obtained for the case of a time varying cold wall boundary condition. Windware explicit differ-encing is used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first order accuracy is maintained in time and space. The results encompass a range of Grashof numbers from 8.34 times 10,000 to 7 times 10 to the 7th power which is within the laminar flow regime for gravitationally driven fluid flows. Experiments within a small scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer and also the decay of wall heat transfer with time.

An Experimental Investigation of Wall-Cooling Effects on Hypersonic Boundary-Layer Stability in a Quiet Wind Tunnel

NASA Technical Reports Server (NTRS)

Blanchard, Alan E.; Selby, Gregory V.

1996-01-01

One of the primary reasons for developing quiet tunnels is for the investigation of high-speed boundary-layer stability and transition phenomena without the transition-promoting effects of acoustic radiation from tunnel walls. In this experiment, a flared-cone model under adiabatic- and cooled-wall conditions was placed in a calibrated, 'quiet' Mach 6 flow and the stability of the boundary layer was investigated using a prototype constant-voltage anemometer. The results were compared with linear-stability theory predictions and good agreement was found in the prediction of second-mode frequencies and growth. In addition, the same 'N=10' criterion used to predict boundary-layer transition in subsonic, transonic, and supersonic flows was found to be applicable for the hypersonic flow regime as well. Under cooled-wall conditions, a unique set of continuous spectra data was acquired that documents the linear, nonlinear, and breakdown regions associated with the transition of hypersonic flow under low-noise conditions.

Industrial application of ultrasound based in-line rheometry: From stationary to pulsating pipe flow of chocolate suspension in precrystallization process

NASA Astrophysics Data System (ADS)

Ouriev, Boris; Windhab, Erich; Braun, Peter; Birkhofer, Beat

2004-10-01

In-line visualization and on-line characterization of nontransparent fluids becomes an important subject for process development in food and nonfood industries. In our work, a noninvasive Doppler ultrasound-based technique is introduced. Such a technique is applied for investigation of nonstationary flow in the chocolate precrystallization process. Unstable flow conditions were induced by abrupt flow interruption and were followed up by strong flow pulsations in the piping system. While relying on available process information, such as absolute pressures and temperatures, no analyses of flow conditions or characterization of suspension properties could possibly be done. It is obvious that chocolate flow properties are sensitive to flow boundary conditions. Therefore, it becomes essential to perform reliable structure state monitoring and particularly in application to nonstationary flow processes. Such flow instabilities in chocolate processing can often lead to failed product quality with interruption of the mainstream production. As will be discussed, a combination of flow velocity profiles, on-line fit into flow profiles, and pressure difference measurement are sufficient for reliable analyses of fluid properties and flow boundary conditions as well as monitoring of the flow state. Analyses of the flow state and flow properties of chocolate suspension are based on on-line measurement of one-dimensional velocity profiles across the flow channel and their on-line characterization with the power-law model. Conclusions about flow boundary conditions were drawn from a calculated velocity standard mean deviation, the parameters of power-law fit into velocity profiles, and volumetric flow rate information.

Numerical simulations of the flow in the HYPULSE expansion tube

NASA Technical Reports Server (NTRS)

Wilson, Gregory J.; Sussman, Myles A.; Bakos, Robert J.

1995-01-01

Axisymmetric numerical simulations with finite-rate chemistry are presented for two operating conditions in the HYPULSE expansion tube. The operating gas for these two cases is nitrogen and the computations are compared to experimental data. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low static pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a static pressure of 38 kPa and is characterized by a turbulent boundary layer. For both cases, the time-varying test gas pressure predicted by the simulations is in good agreement with experimental data. The computations are also found to be in good agreement with Mirels' correlations for shock tube flow. It is shown that the nonuniformity of the test gas observed in the HYPULSE expansion tube is strongly linked to the boundary layer thickness. The turbulent flow investigated has a larger boundary layer and greater test gas nonuniformity. In order to investigate possibilities of improving expansion tube flow quality by reducing the boundary layer thickness, parametric studies showing the effect of density and turbulent transition point on the test conditions are also presented. Although an increase in the expansion tube operating pressure level would reduce the boundary layer thickness, the simulations indicate that the reduction would be less than what is predicted by flat plate boundary layer correlations.

Constant-concentration boundary condition: Lessons from the HYDROCOIN variable-density groundwater benchmark problem

USGS Publications Warehouse

Konikow, Leonard F.; Sanford, W.E.; Campbell, P.J.

1997-01-01

In a solute-transport model, if a constant-concentration boundary condition is applied at a node in an active flow field, a solute flux can occur by both advective and dispersive processes. The potential for advective release is demonstrated by reexamining the Hydrologic Code Intercomparison (HYDROCOIN) project case 5 problem, which represents a salt dome overlain by a shallow groundwater system. The resulting flow field includes significant salinity and fluid density variations. Several independent teams simulated this problem using finite difference or finite element numerical models. We applied a method-of-characteristics model (MOCDENSE). The previous numerical implementations by HYDROCOIN teams of a constant-concentration boundary to represent salt release by lateral dispersion only (as stipulated in the original problem definition) was flawed because this boundary condition allows the release of salt into the flow field by both dispersion and advection. When the constant-concentration boundary is modified to allow salt release by dispersion only, significantly less salt is released into the flow field. The calculated brine distribution for case 5 depends very little on which numerical model is used, as long as the selected model is solving the proper equations. Instead, the accuracy of the solution depends strongly on the proper conceptualization of the problem, including the detailed design of the constant-concentration boundary condition. The importance and sensitivity to the manner of specification of this boundary does not appear to have been recognized previously in the analysis of this problem.

An Evaluation of a Phase-Lag Boundary Condition for Francis Hydroturbine Simulations Using a Pressure-Based Solver

NASA Astrophysics Data System (ADS)

Wouden, Alex; Cimbala, John; Lewis, Bryan

2014-11-01

While the periodic boundary condition is useful for handling rotational symmetry in many axisymmetric geometries, its application fails for analysis of rotor-stator interaction (RSI) in multi-stage turbomachinery flow. The inadequacy arises from the underlying geometry where the blade counts per row differ, since the blade counts are crafted to deter the destructive harmonic forces of synchronous blade passing. Therefore, to achieve the computational advantage of modeling a single blade passage per row while preserving the integrity of the RSI, a phase-lag boundary condition is adapted to OpenFOAM® software's incompressible pressure-based solver. The phase-lag construct is accomplished through restating the implicit periodic boundary condition as a constant boundary condition that is updated at each time step with phase-shifted data from the coupled cells adjacent to the boundary. Its effectiveness is demonstrated using a typical Francis hydroturbine modeled as single- and double-passages with phase-lag boundary conditions. The evaluation of the phase-lag condition is based on the correspondence of the overall computational performance and the calculated flow parameters of the phase-lag simulations with those of a baseline full-wheel simulation. Funded in part by DOE Award Number: DE-EE0002667.

Large eddy simulation for atmospheric boundary layer flow over flat and complex terrains

NASA Astrophysics Data System (ADS)

Han, Yi; Stoellinger, Michael; Naughton, Jonathan

2016-09-01

In this work, we present Large Eddy Simulation (LES) results of atmospheric boundary layer (ABL) flow over complex terrain with neutral stratification using the OpenFOAM-based simulator for on/offshore wind farm applications (SOWFA). The complete work flow to investigate the LES for the ABL over real complex terrain is described including meteorological-tower data analysis, mesh generation and case set-up. New boundary conditions for the lateral and top boundaries are developed and validated to allow inflow and outflow as required in complex terrain simulations. The turbulent inflow data for the terrain simulation is generated using a precursor simulation of a flat and neutral ABL. Conditionally averaged met-tower data is used to specify the conditions for the flat precursor simulation and is also used for comparison with the simulation results of the terrain LES. A qualitative analysis of the simulation results reveals boundary layer separation and recirculation downstream of a prominent ridge that runs across the simulation domain. Comparisons of mean wind speed, standard deviation and direction between the computed results and the conditionally averaged tower data show a reasonable agreement.

Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body

NASA Astrophysics Data System (ADS)

Brykina, I. G.; Rogov, B. V.; Semenov, I. L.; Tirskiy, G. A.

2011-05-01

Super- and hypersonic rarefied gas flow over blunt bodies is investigated by using asymptotically correct viscous shock layer (VSL) model with effective boundary conditions and thin viscous shock layer model. Correct shock and wall conditions for VSL are proposed with taking into account terms due to the curvature which are significant at low Reynolds number. These conditions improve original Davis's VSL model [1]. Numerical calculation of Krook equation [2] is carried out to verify continuum results. Continuum numerical and asymptotic solutions are compared with kinetic solution, free-molecule flow solution and with DSMC solutions [3, 4, 5] over a wide range of free-stream Knudsen number Kn∞. It is shown that taking into account terms with shock and surface curvatures have a pronounced effect on skin friction and heat-transfer in transitional flow regime. Using the asymptotically correct VSL model with effective boundary conditions significantly extends the range of its applicability to higher Kn∞ numbers.

Turbulence Modeling Validation, Testing, and Development

NASA Technical Reports Server (NTRS)

Bardina, J. E.; Huang, P. G.; Coakley, T. J.

1997-01-01

The primary objective of this work is to provide accurate numerical solutions for selected flow fields and to compare and evaluate the performance of selected turbulence models with experimental results. Four popular turbulence models have been tested and validated against experimental data often turbulent flows. The models are: (1) the two-equation k-epsilon model of Wilcox, (2) the two-equation k-epsilon model of Launder and Sharma, (3) the two-equation k-omega/k-epsilon SST model of Menter, and (4) the one-equation model of Spalart and Allmaras. The flows investigated are five free shear flows consisting of a mixing layer, a round jet, a plane jet, a plane wake, and a compressible mixing layer; and five boundary layer flows consisting of an incompressible flat plate, a Mach 5 adiabatic flat plate, a separated boundary layer, an axisymmetric shock-wave/boundary layer interaction, and an RAE 2822 transonic airfoil. The experimental data for these flows are well established and have been extensively used in model developments. The results are shown in the following four sections: Part A describes the equations of motion and boundary conditions; Part B describes the model equations, constants, parameters, boundary conditions, and numerical implementation; and Parts C and D describe the experimental data and the performance of the models in the free-shear flows and the boundary layer flows, respectively.

Thermal Convection in a Creeping Solid With Melting/Freezing Interfaces at Either or Both Boundaries

NASA Astrophysics Data System (ADS)

Labrosse, S.; Morison, A.; Deguen, R.; Alboussiere, T.; Tackley, P. J.; Agrusta, R.

2017-12-01

Thermal convection in the solid mantles of the Earth, other terrestrial planets and icy satellites sets in while it is still crystallising from a liquid layer (see abstract by Morison et al, this conference). The existence of an ocean (water or magma) either or both below and above the solid mantle modifies the conditions applying at the boundary since matter can flow through it by changing phase. Adapting the boundary conditions developed for the dynamics of the inner core by Deguen et al (GJI 2013) to the plane layer and the spherical shell, we solve the linear stability problem and obtain weakly non-linear solutions as well as direct numerical solutions in both geometries, with a liquid-solid phase change at either or both boundaries. The phase change boundary condition is controlled by a dimensionless number, Φ , which when small, allows easy flow through the boundary while the classical non-penetrating boundary condition is recovered for large values. If both boundaries have a phase change, the preferred wavelength of the flow is large, i.e. λ ∝Φ -1/2 in a plane layer and degree 1 in a spherical shell, and the critical Rayleigh number is of order Φ . The heat transfer efficiency, as measured by the dependence of the Nusselt number on the Rayleigh number also increases indefinitely for decreasing values of Φ . If only one boundary has a phase change condition, the critical wavelength is increased by about a factor 2 and the critical Rayleigh number is decreased by about a factor 4. The dynamics is controlled entirely by the boundary layer opposite to the phase change interface and the geometry of the flow. This model provides a natural explanation for the emergence of degree 1 convection in thin ice layers and implies a style of early mantle dynamics on Earth very different from what is classically envisioned.

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

NASA Technical Reports Server (NTRS)

Garai, Anirban; Murman, Scott M.; Madavan, Nateri K.

2016-01-01

The numerical simulation of many aerodynamic non-periodic flows of practical interest involves discretized computational domains that often must be artificially truncated. Appropriate boundary conditions are required at these truncated domain boundaries, and ideally, these boundary conditions should be perfectly "absorbing" or "nonreflecting" so that they do not contaminate the flow field in the interior of the domain. The proper specification of these boundaries is critical to the stability, accuracy, convergence, and quality of the numerical solution, and has been the topic of considerable research. The need for accurate boundary specification has been underscored in recent years with efforts to apply higher-fidelity methods (DNS, LES) in conjunction with high-order low-dissipation numerical schemes to realistic flow configurations. One of the most popular choices for specifying these boundaries is the characteristics-based boundary condition where the linearized flow field at the boundaries are decomposed into characteristic waves using either one-dimensional Riemann or other multi-dimensional Riemann approximations. The values of incoming characteristics are then suitably modified. The incoming characteristics are specified at the in flow boundaries, and at the out flow boundaries the variation of the incoming characteristic is zeroed out to ensure no reflection. This, however, makes the problem ill-posed requiring the use of an ad-hoc parameter to allow small reflections that make the solution stable. Generally speaking, such boundary conditions work reasonably well when the characteristic flow direction is normal to the boundary, but reflects spurious energy otherwise. An alternative to the characteristic-based boundary condition is to add additional "buffer" regions to the main computational domain near the artificial boundaries, and solve a different set of equations in the buffer region in order to minimize acoustic reflections. One approach that has been used involves modeling the pressure fluctuations as acoustic waves propagating in the far-field relative to a single noise-source inside the buffer region. This approach treats vorticity-induced pressure fluctuations the same as acoustic waves. Another popular approach, often referred to as the "sponge layer," attempts to dampen the flow perturbations by introducing artificial dissipation in the buffer region. Although the artificial dissipation removes all perturbations inside the sponge layer, incoming waves are still reflected from the interface boundary between the computational domain and the sponge layer. The effect of these refkections can be somewhat mitigated by appropriately selecting the artificial dissipation strength and the extent of the sponge layer. One of the most promising variants on the buffer region approach is the Perfectly Matched Layer (PML) technique. The PML technique mitigates spurious reflections from boundaries and interfaces by dampening the perturbation modes inside the buffer region such that their eigenfunctions remain unchanged. The technique was first developed by Berenger for application to problems involving electromagnetic wave propagation. It was later extended to the linearized Euler, Euler and Navier-Stokes equations by Hu and his coauthors. The PML technique ensures the no-reflection property for all waves, irrespective of incidence angle, wavelength, and propagation direction. Although the technique requires the solution of a set of auxiliary equations, the computational overhead is easily justified since it allows smaller domain sizes and can provide better accuracy, stability, and convergence of the numerical solution. In this paper, the PML technique is developed in the context of a high-order spectral-element Discontinuous Galerkin (DG) method. The technique is compared to other approaches to treating the in flow and out flow boundary, such as those based on using characteristic boundary conditions and sponge layers. The superiority of the current PML technique over other approaches is demonstrated for a range of test cases, viz., acoustic pulse propagation, convective vortex, shear layer flow, and low-pressure turbine cascade flow. The paper is structured as follows. We first derive the PML equations from the non{linear Euler equations. A short description of the higher-order DG method used is then described. Preliminary results for the four test cases considered are then presented and discussed. Details regarding current work that will be included in the final paper are also provided.

A Boundary Condition for Simulation of Flow Over Porous Surfaces

NASA Technical Reports Server (NTRS)

Frink, Neal T.; Bonhaus, Daryl L.; Vatsa, Veer N.; Bauer, Steven X. S.; Tinetti, Ana F.

2001-01-01

A new boundary condition is presented.for simulating the flow over passively porous surfaces. The model builds on the prior work of R.H. Bush to eliminate the need for constructing grid within an underlying plenum, thereby simplifying the numerical modeling of passively porous flow control systems and reducing computation cost. Code experts.for two structured-grid.flow solvers, TLNS3D and CFL3D. and one unstructured solver, USM3Dns, collaborated with an experimental porosity expert to develop the model and implement it into their respective codes. Results presented,for the three codes on a slender forebody with circumferential porosity and a wing with leading-edge porosity demonstrate a good agreement with experimental data and a remarkable ability to predict the aggregate aerodynamic effects of surface porosity with a simple boundary condition.

A complex-lamellar description of boundary layer transition

NASA Astrophysics Data System (ADS)

Kolla, Maureen Louise

Flow transition is important, in both practical and phenomenological terms. However, there is currently no method for identifying the spatial locations associated with transition, such as the start and end of intermittency. The concept of flow stability and experimental correlations have been used, however, flow stability only identifies the location where disturbances begin to grow in the laminar flow and experimental correlations can only give approximations as measuring the start and end of intermittency is difficult. Therefore, the focus of this work is to construct a method to identify the start and end of intermittency, for a natural boundary layer transition and a separated flow transition. We obtain these locations by deriving a complex-lamellar description of the velocity field that exists between a fully laminar and fully turbulent boundary condition. Mathematically, this complex-lamellar decomposition, which is constructed from the classical Darwin-Lighthill-Hawthorne drift function and the transport of enstrophy, describes the flow that exists between the fully laminar Pohlhausen equations and Prandtl's fully turbulent one seventh power law. We approximate the difference in enstrophy density between the boundary conditions using a power series. The slope of the power series is scaled by using the shape of the universal intermittency distribution within the intermittency region. We solve the complex-lamellar decomposition of the velocity field along with the slope of the difference in enstrophy density function to determine the location of the laminar and turbulent boundary conditions. Then from the difference in enstrophy density function we calculate the start and end of intermittency. We perform this calculation on a natural boundary layer transition over a flat plate for zero pressure gradient flow and for separated shear flow over a separation bubble. We compare these results to existing experimental results and verify the accuracy of our transition model.

Computer program for design of two-dimensional supersonic turbine rotor blades with boundary-layer correction

NASA Technical Reports Server (NTRS)

Goldman, L. J.; Scullin, V. J.

1971-01-01

A FORTRAN 4 computer program for the design of two-dimensional supersonic rotor blade sections corrected for boundary-layer displacement thickness is presented. The ideal rotor is designed by the method of characteristics to produce vortex flow within the blade passage. The boundary-layer parameters are calculated by Cohen and Reshotoko's method for laminar flow and Sasman and Cresci's method for turbulent flow. The program input consists essentially of the blade surface Mach number distribution and total flow conditions. The primary output is the corrected blade profile and the boundary-layer parameters.

Theoretical and Experimental Studies of the Transonic Flow Field and Associated Boundary Conditions near a Longitudinally-Slotted Wind-Tunnel Wall. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Everhart, Joel Lee

1988-01-01

A theoretical examination of the slotted-wall flow field is conducted to determine the appropriate wall pressure drop (or boundary condition) equation. This analysis improves the understanding of the fluid physics of these types of flow fields and helps in evaluating the uncertainties and limitations existing in previous mathematical developments. It is shown that the resulting slotted-wall boundary condition contains contributions from the airfoil-induced streamline curvature and the non-linear, quadratic, slot crossflow in addition to an often neglected linear term which results from viscous shearing in the slot. Existing and newly acquired experimental data are examined in the light of this formulation and theoretical developments.

Multi-resolution MPS method

NASA Astrophysics Data System (ADS)

Tanaka, Masayuki; Cardoso, Rui; Bahai, Hamid

2018-04-01

In this work, the Moving Particle Semi-implicit (MPS) method is enhanced for multi-resolution problems with different resolutions at different parts of the domain utilising a particle splitting algorithm for the finer resolution and a particle merging algorithm for the coarser resolution. The Least Square MPS (LSMPS) method is used for higher stability and accuracy. Novel boundary conditions are developed for the treatment of wall and pressure boundaries for the Multi-Resolution LSMPS method. A wall is represented by polygons for effective simulations of fluid flows with complex wall geometries and the pressure boundary condition allows arbitrary inflow and outflow, making the method easier to be used in flow simulations of channel flows. By conducting simulations of channel flows and free surface flows, the accuracy of the proposed method was verified.

In-Flight Boundary-Layer Transition of a Large Flat Plate at Supersonic Speeds

NASA Technical Reports Server (NTRS)

Banks, D. W.; Frederick, M. A.; Tracy, R. R.; Matisheck, J. R.; Vanecek, N. D.

2012-01-01

A flight experiment was conducted to investigate the pressure distribution, local-flow conditions, and boundary-layer transition characteristics on a large flat plate in flight at supersonic speeds up to Mach 2.00. The tests used a NASA testbed aircraft with a bottom centerline mounted test fixture. The primary objective of the test was to characterize the local flow field in preparation for future tests of a high Reynolds number natural laminar flow test article. A second objective was to determine the boundary-layer transition characteristics on the flat plate and the effectiveness of using a simplified surface coating. Boundary-layer transition was captured in both analog and digital formats using an onboard infrared imaging system. Surface pressures were measured on the surface of the flat plate. Flow field measurements near the leading edge of the test fixture revealed the local flow characteristics including downwash, sidewash, and local Mach number. Results also indicated that the simplified surface coating did not provide sufficient insulation from the metallic structure, which likely had a substantial effect on boundary-layer transition compared with that of an adiabatic surface. Cold wall conditions were predominant during the acceleration to maximum Mach number, and warm wall conditions were evident during the subsequent deceleration.

Numerical Treatment of Stokes Solvent Flow and Solute-Solvent Interfacial Dynamics for Nonpolar Molecules.

PubMed

Sun, Hui; Zhou, Shenggao; Moore, David K; Cheng, Li-Tien; Li, Bo

2016-05-01

We design and implement numerical methods for the incompressible Stokes solvent flow and solute-solvent interface motion for nonpolar molecules in aqueous solvent. The balance of viscous force, surface tension, and van der Waals type dispersive force leads to a traction boundary condition on the solute-solvent interface. To allow the change of solute volume, we design special numerical boundary conditions on the boundary of a computational domain through a consistency condition. We use a finite difference ghost fluid scheme to discretize the Stokes equation with such boundary conditions. The method is tested to have a second-order accuracy. We combine this ghost fluid method with the level-set method to simulate the motion of the solute-solvent interface that is governed by the solvent fluid velocity. Numerical examples show that our method can predict accurately the blow up time for a test example of curvature flow and reproduce the polymodal (e.g., dry and wet) states of hydration of some simple model molecular systems.

Numerical Treatment of Stokes Solvent Flow and Solute-Solvent Interfacial Dynamics for Nonpolar Molecules

PubMed Central

Sun, Hui; Zhou, Shenggao; Moore, David K.; Cheng, Li-Tien; Li, Bo

2015-01-01

We design and implement numerical methods for the incompressible Stokes solvent flow and solute-solvent interface motion for nonpolar molecules in aqueous solvent. The balance of viscous force, surface tension, and van der Waals type dispersive force leads to a traction boundary condition on the solute-solvent interface. To allow the change of solute volume, we design special numerical boundary conditions on the boundary of a computational domain through a consistency condition. We use a finite difference ghost fluid scheme to discretize the Stokes equation with such boundary conditions. The method is tested to have a second-order accuracy. We combine this ghost fluid method with the level-set method to simulate the motion of the solute-solvent interface that is governed by the solvent fluid velocity. Numerical examples show that our method can predict accurately the blow up time for a test example of curvature flow and reproduce the polymodal (e.g., dry and wet) states of hydration of some simple model molecular systems. PMID:27365866

Intermittent Behavior of the Separated Boundary Layer along the Suction Surface of a Low Pressure Turbine Blade under Periodic Unsteady Flow Conditions

NASA Technical Reports Server (NTRS)

Oeztuerk, B; Schobeiri, M. T.; Ashpis, David E.

2005-01-01

The paper experimentally and theoretically studies the effects of periodic unsteady wake flow and aerodynamic characteristics on boundary layer development, separation and re-attachment along the suction surface of a low pressure turbine blade. The experiments were carried out at Reynolds number of 110,000 (based on suction surface length and exit velocity). For one steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, intermittency behaviors were experimentally and theoretically investigated. The current investigation attempts to extend the intermittency unsteady boundary layer transition model developed in previously to the LPT cases, where separation occurs on the suction surface at a low Reynolds number. The results of the unsteady boundary layer measurements and the intermittency analysis were presented in the ensemble-averaged and contour plot forms. The analysis of the boundary layer experimental data with the flow separation, confirms the universal character of the relative intermittency function which is described by a Gausssian function.

Sediment Vertical Flux in Unsteady Sheet Flows

NASA Astrophysics Data System (ADS)

Hsu, T.; Jenkins, J. T.; Liu, P. L.

2002-12-01

In models for sediment suspension, two different boundary conditions have been employed at the sediment bed. Either the sediment concentration is given or the vertical flux of sediment is specified. The specification of the latter is usually called the pick-up function. Recently, several developments towards a better understanding of the sediment bed boundary condition have been reported. Nielson et al (Coastal Engineering 2002, 45, p61-68) have indicated a better performance using the sediment vertical flux as the bed boundary condition in comparisons with experimental data. Also, Drake and Calantoni (Journal of Geophysical Research 2001, 106, C9, p19859-19868) have suggested that in the nearshore environment with its various unsteady flow conditions, the appropriate sediment boundary conditions of a large-scale morphology model must consider both the magnitude the free stream velocity and the acceleration of the flow. In this research, a small-scale sheet flow model based on the two-phase theory is implemented to further study these issues. Averaged two-phase continuum equations are presented for concentrated flows of sediment that are driven by strong, fully developed, unsteady turbulent shear flows over a mobile bed. The particle inter-granular stress is modeled using collisional granular flow theory and a two-equation closure for the fluid turbulence is adopted. In the context of the two-phase theory, sediment is transported through the sediment vertical velocity. Using the fully developed sediment phase continuity equation, it can be shown that the vertical velocity of the sediment must vanish when the flow reaches a steady state. In other words, in fully developed conditions, it is the unsteadiness of the flow that induces the vertical motion of the sediment and that changes the sediment concentration profile. Therefore, implementing a boundary condition based on sediment vertical flux is consistent with both the two-phase theory and with the observation that the flow acceleration is an important parameter. In this paper, the vertical flux of sediment is studied under various combinations of free stream velocity, acceleration, and sediment material properties using the two-phase sheet flow model. Some interesting features of sediment dynamics within the sheet, such as time history of sediment vertical velocity, collisional and turbulent suspension mechanisms are presented.

Tidal Boundary Conditions in SEAWAT

USGS Publications Warehouse

Mulligan, Ann E.; Langevin, Christian; Post, Vincent E.A.

2011-01-01

SEAWAT, a U.S. Geological Survey groundwater flow and transport code, is increasingly used to model the effects of tidal motion on coastal aquifers. Different options are available to simulate tidal boundaries but no guidelines exist nor have comparisons been made to identify the most effective approach. We test seven methods to simulate a sloping beach and a tidal flat. The ocean is represented in one of the three ways: directly using a high hydraulic conductivity (high-K) zone and indirect simulation via specified head boundaries using either the General Head Boundary (GHB) or the new Periodic Boundary Condition (PBC) package. All beach models simulate similar water fluxes across the upland boundary and across the sediment-water interface although the ratio of intertidal to subtidal flow is different at low tide. Simulating a seepage face results in larger intertidal fluxes and influences near-shore heads and salinity. Major differences in flow occur in the tidal flat simulations. Because SEAWAT does not simulate unsaturated flow the water table only rises via flow through the saturated zone. This results in delayed propagation of the rising tidal signal inland. Inundation of the tidal flat is delayed as is flow into the aquifer across the flat. This is severe in the high-K and PBC models but mild in the GHB models. Results indicate that any of the tidal boundary options are fine if the ocean-aquifer interface is steep. However, as the slope of that interface decreases, the high-K and PBC approaches perform poorly and the GHB boundary is preferable.

The effects of boundary conditions on the steady-state response of three hypothetical ground-water systems; results and implications of numerical experiments

USGS Publications Warehouse

Franke, O. Lehn; Reilly, Thomas E.

1987-01-01

The most critical and difficult aspect of defining a groundwater system or problem for conceptual analysis or numerical simulation is the selection of boundary conditions . This report demonstrates the effects of different boundary conditions on the steady-state response of otherwise similar ground-water systems to a pumping stress. Three series of numerical experiments illustrate the behavior of three hypothetical groundwater systems that are rectangular sand prisms with the same dimensions but with different combinations of constant-head, specified-head, no-flow, and constant-flux boundary conditions. In the first series of numerical experiments, the heads and flows in all three systems are identical, as are the hydraulic conductivity and system geometry . However, when the systems are subjected to an equal stress by a pumping well in the third series, each differs significantly in its response . The highest heads (smallest drawdowns) and flows occur in the systems most constrained by constant- or specified-head boundaries. These and other observations described herein are important in steady-state calibration, which is an integral part of simulating many ground-water systems. Because the effects of boundary conditions on model response often become evident only when the system is stressed, a close match between the potential distribution in the model and that in the unstressed natural system does not guarantee that the model boundary conditions correctly represent those in the natural system . In conclusion, the boundary conditions that are selected for simulation of a ground-water system are fundamentally important to groundwater systems analysis and warrant continual reevaluation and modification as investigation proceeds and new information and understanding are acquired.

Solute transport with time-variable flow paths during upward and downward flux in a heterogeneous unsaturated porous medium

NASA Astrophysics Data System (ADS)

Cremer, Clemens; Neuweiler, Insa; Bechtold, Michel; Vanderborght, Jan

2014-05-01

To acquire knowledge of solute transport through the unsaturated zone in the shallow subsurface is decisive to assess groundwater quality, nutrient cycling or to plan remediation strategies. The shallow subsurface is characterized by structural heterogeneity and strongly influenced by atmospheric conditions. This leads to changing flow directions, strong temporal changes in saturation and heterogeneous water fluxes during infiltration and evaporation events. Recent studies (e.g. Lehmann and Or, 2009; Bechtold et al.,2011) demonstrated the importance of lateral flow and solute transport during evaporation conditions (upward flux). The heterogeneous structure in these studies was constructed using two types of sand with strong material contrasts and arranged in parallel with a vertical orientation. Lateral transport and redistribution of solute from coarse to fine media was observed deeper in the soil column and from fine to coarse close to the soil surface. However, if boundary conditions are reversed due to precipitation, the flow field is not necessarily reversed in the same manner, resulting in entirely different transport patterns for downward and upward flow. Therefore, considering net-flow rates alone is misleading when describing transport under those conditions. In this contribution we analyze transport of a solute in the shallow subsurface to assess effects resulting from the temporal change of heterogeneous soil structures due to dynamic flow conditions. Two-dimensional numerical simulations of unsaturated flow and transport are conducted using a coupled finite volume and random walk particle tracking algorithm to quantify solute transport and leaching rates. Following previous studies (Lehmann and Or, 2009; Bechtold et al., 2011), the chosen domain is composed of two materials, coarse and fine sand, arranged in parallel with a vertical orientation. Hence, one sharp interface of strong material heterogeneity is induced. During evaporation both sands are assumed to stay under liquid-flow dominated evaporation conditions ("stage 1"). Simulations considering dynamic (infiltration-evaporation) and steady (solely infiltration) boundary conditions are carried out. The influence of dynamic boundary conditions (intensity and duration of precipitation and evaporation events) is examined in a multitude of simulations. If flow rates smaller than the saturated hydraulic conductivity of both materials are chosen to be applied as boundary condition, simulation results indicate that the flow field within the domain is exactly reversed. However, if applied flow rates exceed the saturated hydraulic conductivity of one material, the flow field is not just reversed, but different flow paths during downward and upward flow are observed. Results show the tendency of faster solute leaching under dynamic boundary conditions compared to steady infiltration conditions with the same net-infiltration rate. We use a double domain transport method as an upscaled model to reproduce vertically averaged concentration profiles with net flux only and compare the model parameters for information about flow dynamics and soil heterogeneity.

Adaptive wing and flow control technology

NASA Astrophysics Data System (ADS)

Stanewsky, E.

2001-10-01

The development of the boundary layer and the interaction of the boundary layer with the outer “inviscid” flow field, exacerbated at high speed by the occurrence of shock waves, essentially determine the performance boundaries of high-speed flight. Furthermore, flight and freestream conditions may change considerably during an aircraft mission while the aircraft itself is only designed for multiple but fixed design points thus impairing overall performance. Consequently, flow and boundary layer control and adaptive wing technology may have revolutionary new benefits for take-off, landing and cruise operating conditions for many aircraft by enabling real-time effective geometry optimization relative to the flight conditions. In this paper we will consider various conventional and novel means of boundary layer and flow control applied to moderate-to-large aspect ratio wings, delta wings and bodies with the specific objectives of drag reduction, lift enhancement, separation suppression and the improvement of air-vehicle control effectiveness. In addition, adaptive wing concepts of varying complexity and corresponding aerodynamic performance gains will be discussed, also giving some examples of possible structural realizations. Furthermore, penalties associated with the implementation of control and adaptation mechanisms into actual aircraft will be addressed. Note that the present contribution is rather application oriented.

DESIGN OF TWO-DIMENSIONAL SUPERSONIC TURBINE ROTOR BLADES WITH BOUNDARY-LAYER CORRECTION

NASA Technical Reports Server (NTRS)

Goldman, L. J.

1994-01-01

A computer program has been developed for the design of supersonic rotor blades where losses are accounted for by correcting the ideal blade geometry for boundary layer displacement thickness. The ideal blade passage is designed by the method of characteristics and is based on establishing vortex flow within the passage. Boundary-layer parameters (displacement and momentum thicknesses) are calculated for the ideal passage, and the final blade geometry is obtained by adding the displacement thicknesses to the ideal nozzle coordinates. The boundary-layer parameters are also used to calculate the aftermixing conditions downstream of the rotor blades assuming the flow mixes to a uniform state. The computer program input consists essentially of the rotor inlet and outlet Mach numbers, upper- and lower-surface Mach numbers, inlet flow angle, specific heat ratio, and total flow conditions. The program gas properties are set up for air. Additional gases require changes to be made to the program. The computer output consists of the corrected rotor blade coordinates, the principal boundary-layer parameters, and the aftermixing conditions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 7094. This program was developed in 1971.

Group solution for unsteady free-convection flow from a vertical moving plate subjected to constant heat flux

NASA Astrophysics Data System (ADS)

Kassem, M.

2006-03-01

The problem of heat and mass transfer in an unsteady free-convection flow over a continuous moving vertical sheet in an ambient fluid is investigated for constant heat flux using the group theoretical method. The nonlinear coupled partial differential equation governing the flow and the boundary conditions are transformed to a system of ordinary differential equations with appropriate boundary conditions. The obtained ordinary differential equations are solved numerically using the shooting method. The effect of Prandlt number on the velocity and temperature of the boundary-layer is plotted in curves. A comparison with previous work is presented.

Multicomponent Gas Diffusion and an Appropriate Momentum Boundary Condition

NASA Technical Reports Server (NTRS)

Noever, David A.

1994-01-01

Multicomponent gas diffusion is reviewed with particular emphasis on gas flows near solid boundaries-the so-called Kramers-Kistemaker effect. The aim is to derive an appropriate momentum boundary condition which governs many gaseous species diffusing together. The many species' generalization of the traditional single gas condition, either as slip or stick (no-slip), is not obvious, particularly for technologically important cases of lower gas pressures and very dissimilar molecular weight gases. No convincing theoretical case exists for why two gases should interact with solid boundaries equally but in opposite flow directions, such that the total gas flow exactly vanishes. ln this way, the multicomponent no-slip boundary requires careful treatment The approaches discussed here generally adopt a microscopic model for gas-solid contact. The method has the advantage that the mathematics remain tractable and hence experimentally testable. Two new proposals are put forward, the first building in some molecular collision physics, the second drawing on a detailed view of surface diffusion which does not unphysically extrapolate bulk gas properties to govern the adsorbed molecules. The outcome is a better accounting of previously anomalous experiments. Models predict novel slip conditions appearing even for the case of equal molecular weight components. These approaches become particularly significant in view of a conceptual contradiction found to arise in previous derivations of the appropriate boundary conditions. The analogous case of three gases, one of which is uniformly distributed and hence non-diffusing, presents a further refinement which gives unexpected flow reversals near solid boundaries. This case is investigated alone and for aggregating gas species near their condensation point. In addition to predicting new physics, this investigation carries practical implications for controlling vapor diffusion in the growth of crystals used in medical diagnosis (e.g. mercuric iodide) and semiconductors.

Flow analysis for the nacelle of an advanced ducted propeller at high angle-of-attack and at cruise with boundary layer control

NASA Technical Reports Server (NTRS)

Hwang, D. P.; Boldman, D. R.; Hughes, C. E.

1994-01-01

An axisymmetric panel code and a three dimensional Navier-Stokes code (used as an inviscid Euler code) were verified for low speed, high angle of attack flow conditions. A three dimensional Navier-Stokes code (used as an inviscid code), and an axisymmetric Navier-Stokes code (used as both viscous and inviscid code) were also assessed for high Mach number cruise conditions. The boundary layer calculations were made by using the results from the panel code or Euler calculation. The panel method can predict the internal surface pressure distributions very well if no shock exists. However, only Euler and Navier-Stokes calculations can provide a good prediction of the surface static pressure distribution including the pressure rise across the shock. Because of the high CPU time required for a three dimensional Navier-Stokes calculation, only the axisymmetric Navier-Stokes calculation was considered at cruise conditions. The use of suction and tangential blowing boundary layer control to eliminate the flow separation on the internal surface was demonstrated for low free stream Mach number and high angle of attack cases. The calculation also shows that transition from laminar flow to turbulent flow on the external cowl surface can be delayed by using suction boundary layer control at cruise flow conditions. The results were compared with experimental data where possible.

Numerical method for predicting flow characteristics and performance of nonaxisymmetric nozzles, theory

NASA Technical Reports Server (NTRS)

Thomas, P. D.

1979-01-01

The theoretical foundation and formulation of a numerical method for predicting the viscous flowfield in and about isolated three dimensional nozzles of geometrically complex configuration are presented. High Reynolds number turbulent flows are of primary interest for any combination of subsonic, transonic, and supersonic flow conditions inside or outside the nozzle. An alternating-direction implicit (ADI) numerical technique is employed to integrate the unsteady Navier-Stokes equations until an asymptotic steady-state solution is reached. Boundary conditions are computed with an implicit technique compatible with the ADI technique employed at interior points of the flow region. The equations are formulated and solved in a boundary-conforming curvilinear coordinate system. The curvilinear coordinate system and computational grid is generated numerically as the solution to an elliptic boundary value problem. A method is developed that automatically adjusts the elliptic system so that the interior grid spacing is controlled directly by the a priori selection of the grid spacing on the boundaries of the flow region.

Study of the Effect of Wind Speed on Evaporation from Soil Through Integrated Modeling of Atmospheric Boundary Layer and Shallow Subsurface

NASA Astrophysics Data System (ADS)

Smits, K. M.; Davarzani, H.; Illangasekare, T. H.

2012-12-01

The study of the interaction between the land and atmosphere is paramount to our understanding of many emerging problems to include climate change and the movement of green house gases such as possible leaking of sequestered CO2. Soil moisture distribution in the shallow subsurface becomes a critical factor in these problems. The heat and mass flux in the form of soil evaporation across the land surface couples the atmospheric boundary layer to the shallow subsurface. The coupling between land and the atmosphere leads to highly dynamic interactions between the porous media properties, transport processes and boundary conditions, resulting in dynamic evaporative behavior. However, the coupling at the land-atmospheric interface is rarely considered in most current models and their validation for practical applications. This is due to the complexity of the problem in field scenarios and the scarcity of field or laboratory data capable of testing and refining coupled energy and mass transfer theories. In most efforts to compute evaporation from soil, only indirect coupling is provided to characterize the interaction between non-isothermal multiphase flows under realistic atmospheric conditions even though heat and mass flux are controlled by the coupled dynamics of the land and the atmospheric boundary layer. In earlier drying modeling concepts, imposing evaporation flux (kinetic of relative humidity) and temperature as surface boundary condition is often needed. With the goal of improving our understanding of the land/atmospheric coupling, we developed a model based on the coupling of Navier-Stokes free flow and Darcy flow in porous medium. The model consists of the coupled equations of mass conservation for the liquid phase (water) and gas phase (water vapor and air) in porous medium with gas phase (water vapor and air) in free flow domain under non-isothermal, non-equilibrium conditions. The boundary conditions at the porous medium-free flow medium interface include dynamical, thermal and solutal equilibriums, and using the Beavers-Joseph slip boundary condition. What is unique about this model is that the evaporation rate and soil surface temperature conditions come directly from the model output. In order to experimentally validate the numerical results, we developed and used a unique two dimensional wind tunnel placed above a soil tank equipped with a network of different sensors. A series of experiments under varying boundary conditions were performed. Precision data under well-controlled transient heat and wind boundary conditions was generated. Results from numerical simulations were compared with experimental data. Results demonstrate that the coupling concept can predict the different stages of the drying process in porous media with good accuracy. Increasing the wind speed increases the first stage evaporation rate and decreases the transition time at low velocity values; then, at high values of wind speed the evaporation rate becomes less dependent of flow in free fluid. In the opposite, the impact of the wind speed on the second stage evaporation (diffusion dominant stage) is not significant. The proposed theoretical model can be used to predict the evaporation process where a porous medium flow is coupled to a free flow for different practical applications.

A VERSATILE SHARP INTERFACE IMMERSED BOUNDARY METHOD FOR INCOMPRESSIBLE FLOWS WITH COMPLEX BOUNDARIES

PubMed Central

Mittal, R.; Dong, H.; Bozkurttas, M.; Najjar, F.M.; Vargas, A.; von Loebbecke, A.

2010-01-01

A sharp interface immersed boundary method for simulating incompressible viscous flow past three-dimensional immersed bodies is described. The method employs a multi-dimensional ghost-cell methodology to satisfy the boundary conditions on the immersed boundary and the method is designed to handle highly complex three-dimensional, stationary, moving and/or deforming bodies. The complex immersed surfaces are represented by grids consisting of unstructured triangular elements; while the flow is computed on non-uniform Cartesian grids. The paper describes the salient features of the methodology with special emphasis on the immersed boundary treatment for stationary and moving boundaries. Simulations of a number of canonical two- and three-dimensional flows are used to verify the accuracy and fidelity of the solver over a range of Reynolds numbers. Flow past suddenly accelerated bodies are used to validate the solver for moving boundary problems. Finally two cases inspired from biology with highly complex three-dimensional bodies are simulated in order to demonstrate the versatility of the method. PMID:20216919

Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to Robin boundary condition

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre M.

2014-02-15

Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to variousmore » forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.« less

CFD determination of flow perturbation boundary conditions for seal rotordynamic modeling

NASA Astrophysics Data System (ADS)

Venkatesan, Ganesh

2002-09-01

A new approach has been developed and utilized to determine the flow field perturbations (i.e. disturbance due to rotor eccentricity and/or motion) upstream of and within a non-contacting seal. The results are proposed for use with bulk-flow perturbation and CFD-perturbation seal rotordynamic models, as well as in fully 3-D CFD models, to specify approximate boundary conditions for the first-order variables at the computational domain inlet. The perturbation quantities were evaluated by subtracting the numerical flow field solutions corresponding to the concentric rotor position from that for an eccentric rotor position. The disturbance pressure quantities predicted from the numerical solutions were validated by comparing with previous pressure measurements. A parametric study was performed to understand the influence of upstream chamber height, seal clearance, shaft speed, whirl speed, zeroth-order streamwise and swirl velocities, and downstream pressure on the distribution of the first-order quantities in the upstream chamber, seal inlet and seal exit regions. Radially bulk-averaged first-order quantities were evaluated in the upstream chamber, as well as at the seal inlet and exit. The results were finally presented in the form of generalized dimensionless boundary condition correlations so that they can be applied to seal rotordynamic models over a wide range of operating conditions and geometries. To examine the effect of the proposed, approximate first-order boundary conditions on the solutions of the fully 3-D CFD rotordynamic models, the first-order boundary condition correlations for the upstream chamber were used to adjust the circumferential distribution of domain inlet values. The benefit of the boundary condition expressions was assessed for two previously measured test cases, one for a gas seal and the other for a liquid seal. For the gas seal case, a significant improvement in the prediction of the cross-coupled stiffness, when including the proposed first-order inlet boundary values, was found. In the case of liquid seals the tangential impedance values obtained with boundary condition adjustments showed a very slight improvement for a range of whirl speeds over those obtained without them. The radial impedance values obtained with the new adjustments showed a significant improvement over those obtained without them.

Artificial Boundary Conditions for Computation of Oscillating External Flows

NASA Technical Reports Server (NTRS)

Tsynkov, S. V.

1996-01-01

In this paper, we propose a new technique for the numerical treatment of external flow problems with oscillatory behavior of the solution in time. Specifically, we consider the case of unbounded compressible viscous plane flow past a finite body (airfoil). Oscillations of the flow in time may be caused by the time-periodic injection of fluid into the boundary layer, which in accordance with experimental data, may essentially increase the performance of the airfoil. To conduct the actual computations, we have to somehow restrict the original unbounded domain, that is, to introduce an artificial (external) boundary and to further consider only a finite computational domain. Consequently, we will need to formulate some artificial boundary conditions (ABC's) at the introduced external boundary. The ABC's we are aiming to obtain must meet a fundamental requirement. One should be able to uniquely complement the solution calculated inside the finite computational domain to its infinite exterior so that the original problem is solved within the desired accuracy. Our construction of such ABC's for oscillating flows is based on an essential assumption: the Navier-Stokes equations can be linearized in the far field against the free-stream back- ground. To actually compute the ABC's, we represent the far-field solution as a Fourier series in time and then apply the Difference Potentials Method (DPM) of V. S. Ryaben'kii. This paper contains a general theoretical description of the algorithm for setting the DPM-based ABC's for time-periodic external flows. Based on our experience in implementing analogous ABC's for steady-state problems (a simpler case), we expect that these boundary conditions will become an effective tool for constructing robust numerical methods to calculate oscillatory flows.

In-Flight Boundary-Layer Transition on a Large Flat Plate at Supersonic Speeds

NASA Technical Reports Server (NTRS)

Banks, Daniel W.; Fredericks, Michael Alan; Tracy, Richard R.; Matisheck, Jason R.; Vanecek, Neal D.

2012-01-01

A flight experiment was conducted to investigate the pressure distribution, local flow conditions, and boundary-layer transition characteristics on a large flat plate in flight at supersonic speeds up to Mach 2.0. The primary objective of the test was to characterize the local flow field in preparation for future tests of a high Reynolds number natural laminar flow test article. The tests used a F-15B testbed aircraft with a bottom centerline mounted test fixture. A second objective was to determine the boundary-layer transition characteristics on the flat plate and the effectiveness of using a simplified surface coating for future laminar flow flight tests employing infrared thermography. Boundary-layer transition was captured using an onboard infrared imaging system. The infrared imagery was captured in both analog and digital formats. Surface pressures were measured with electronically scanned pressure modules connected to 60 surface-mounted pressure orifices. The local flow field was measured with five 5-hole conical probes mounted near the leading edge of the test fixture. Flow field measurements revealed the local flow characteristics including downwash, sidewash, and local Mach number. Results also indicated that the simplified surface coating did not provide sufficient insulation from the metallic structure, which likely had a substantial effect on boundary-layer transition compared with that of an adiabatic surface. Cold wall conditions were predominant during the acceleration to maximum Mach number, and warm wall conditions were evident during the subsequent deceleration. The infrared imaging system was able to capture shock wave impingement on the surface of the flat plate in addition to indicating laminar-to-turbulent boundary-layer transition.

Intrawellbore kinematic and frictional losses in a horizontal well in a bounded confined aquifer

NASA Astrophysics Data System (ADS)

Wang, Quanrong; Zhan, Hongbin

2017-01-01

Horizontal drilling has become an appealing technology for water resource exploration or aquifer remediation in recent decades, due to decreasing operational cost and many technical advantages over vertical wells. However, many previous studies on flow into horizontal wells were based on the Uniform Flux Boundary Condition (UFBC), which does not reflect the physical processes of flow inside the well accurately. In this study, we investigated transient flow into a horizontal well in an anisotropic confined aquifer laterally bounded by two constant-head boundaries. Three types of boundary conditions were employed to treat the horizontal well, including UFBC, Uniform-Head Boundary Condition (UHBC), and Mixed-Type Boundary Condition (MTBC). The MTBC model considered both kinematic and frictional effects inside the horizontal well, in which the kinematic effect referred to the accelerational and fluid-inflow effects. A new solution of UFBC was derived by superimposing the point sink/source solutions along the axis of a horizontal well with a uniform flux distribution. New solutions of UHBC and MTBC were obtained by a hybrid analytical-numerical method, and an iterative method was proposed to determine the well discretization required for achieving sufficiently accurate results. This study showed that the differences among the UFBC, UHBC, and MTBC solutions were obvious near the well screen, decreased with distance from the well, and became negligible near the constant-head boundary. The relationship between the flow rate and the drawdown was nonlinear for the MTBC solution, while it was linear for the UFBC and UHBC solutions.

Numerical solutions of the Navier-Stokes equations for the supersonic laminar flow over a two-dimensional compression corner

NASA Technical Reports Server (NTRS)

Carter, J. E.

1972-01-01

Numerical solutions have been obtained for the supersonic, laminar flow over a two-dimensional compression corner. These solutions were obtained as steady-state solutions to the unsteady Navier-Stokes equations using the finite difference method of Brailovskaya, which has second-order accuracy in the spatial coordinates. Good agreement was obtained between the computed results and wall pressure distributions measured experimentally for Mach numbers of 4 and 6.06, and respective Reynolds numbers, based on free-stream conditions and the distance from the leading edge to the corner. In those calculations, as well as in others, sufficient resolution was obtained to show the streamline pattern in the separation bubble. Upstream boundary conditions to the compression corner flow were provided by numerically solving the unsteady Navier-Stokes equations for the flat plate flow field, beginning at the leading edge. The compression corner flow field was enclosed by a computational boundary with the unknown boundary conditions supplied by extrapolation from internally computed points.

Nested large-eddy simulations of nighttime shear-instability waves and transient warming in a steep valley

NASA Astrophysics Data System (ADS)

Zhou, Bowen; Chow, Fotini

2012-11-01

This numerical study investigates the nighttime flow dynamics in a steep valley. The Owens Valley in California is highly complex, and represents a challenging terrain for large-eddy simulations (LES). To ensure a faithful representation of the nighttime atmospheric boundary layer (ABL), realistic external boundary conditions are provided through grid nesting. The model obtains initial and lateral boundary conditions from reanalysis data, and bottom boundary conditions from a land-surface model. We demonstrate the ability to extend a mesoscale model to LES resolutions through a systematic grid-nesting framework, achieving accurate simulations of the stable ABL over complex terrain. Nighttime cold-air flow was channeled through a gap on the valley sidewall. The resulting katabatic current induced a cross-valley flow. Directional shear against the down-valley flow in the lower layers of the valley led to breaking Kelvin-Helmholtz waves at the interface, which is captured only on the LES grid. Later that night, the flow transitioned from down-slope to down-valley near the western sidewall, leading to a transient warming episode. Simulation results are verified against field observations and reveal good spatial and temporal precision. Supported by NSF grant ATM-0645784.

Numerical Simulation of Dynamic Contact Angles and Contact Lines in Multiphase Flows using Level Set Method

NASA Astrophysics Data System (ADS)

Pendota, Premchand

Many physical phenomena and industrial applications involve multiphase fluid flows and hence it is of high importance to be able to simulate various aspects of these flows accurately. The Dynamic Contact Angles (DCA) and the contact lines at the wall boundaries are a couple of such important aspects. In the past few decades, many mathematical models were developed for predicting the contact angles of the inter-face with the wall boundary under various flow conditions. These models are used to incorporate the physics of DCA and contact line motion in numerical simulations using various interface capturing/tracking techniques. In the current thesis, a simple approach to incorporate the static and dynamic contact angle boundary conditions using the level set method is developed and implemented in multiphase CFD codes, LIT (Level set Interface Tracking) (Herrmann (2008)) and NGA (flow solver) (Desjardins et al (2008)). Various DCA models and associated boundary conditions are reviewed. In addition, numerical aspects such as the occurrence of a stress singularity at the contact lines and grid convergence of macroscopic interface shape are dealt with in the context of the level set approach.

A Cartesian-based embedded geometry technique with adaptive high-order finite differences for compressible flow around complex geometries

NASA Astrophysics Data System (ADS)

Uddin, H.; Kramer, R. M. J.; Pantano, C.

2014-04-01

An immersed boundary methodology to solve the compressible Navier-Stokes equations around complex geometries in Cartesian fluid dynamics solvers is described. The objective of the new approach is to enable smooth reconstruction of pressure and viscous stresses around the embedded objects without spurious numerical artifacts. A standard level set represents the boundary of the object and defines a fictitious domain into which the flow fields are smoothly extended. Boundary conditions on the surface are enforced by an approach inspired by analytic continuation. Each fluid field is extended independently, constrained only by the boundary condition associated with that field. Unlike most existing methods, no jump conditions or explicit derivation of them from the boundary conditions are required in this approach. Numerical stiffness that arises when the fluid-solid interface is close to grid points of the mesh is addressed by preconditioning. In addition, the embedded geometry technique is coupled with a stable high-order adaptive discretization that is enabled around the object boundary to enhance resolution. The stencils used to transition the order of accuracy of the discretization are derived using the summation-by-parts technique that ensures stability. Applications to shock reflections, shock-ramp interactions, and supersonic and low-Mach number flows over two- and three-dimensional geometries are presented.

BLSTA: A boundary layer code for stability analysis

NASA Technical Reports Server (NTRS)

Wie, Yong-Sun

1992-01-01

A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.

CFD Modeling of a CFB Riser Using Improved Inlet Boundary Conditions

NASA Astrophysics Data System (ADS)

Peng, B. T.; Zhang, C.; Zhu, J. X.; Qi, X. B.

2010-03-01

A computational fluid dynamics (CFD) model based on Eulerian-Eulerian approach coupled with granular kinetics theory was adopted to investigate the hydrodynamics and flow structures in a circulating fluidized bed (CFB) riser column. A new approach to specify the inlet boundary conditions was proposed in this study to simulate gas-solids flow in CFB risers more accurately. Simulation results were compared with the experimental data, and good agreement between the numerical results and experimental data was observed under different operating conditions, which indicates the effectiveness and accuracy of the CFD model with the proposed inlet boundary conditions. The results also illustrate a clear core annulus structure in the CFB riser under all operating conditions both experimentally and numerically.

Air flow in the boundary layer near a plate

NASA Technical Reports Server (NTRS)

Dryden, Hugh L

1937-01-01

The published data on the distribution of speed near a thin flat plate with sharp leading edge placed parallel to the flow (skin friction plate) are reviewed and the results of some additional measurements are described. The purpose of the experiments was to study the basic phenomena of boundary-layer flow under simple conditions.

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

NASA Astrophysics Data System (ADS)

Zhang, Yu-Dong; Xu, Ai-Guo; Zhang, Guang-Cai; Chen, Zhi-Hua

2018-01-01

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model (DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation. It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model. Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model. To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified. Support of National Natural Science Foundation of China under Grant Nos. 11475028, 11772064, and 11502117 Science Challenge Project under Grant Nos. JCKY2016212A501 and TZ2016002

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

NASA Technical Reports Server (NTRS)

Tsynkov, Semyon V.

1997-01-01

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

A theoretical method for the analysis and design of axisymmetric bodies. [flow distribution and incompressible fluids

NASA Technical Reports Server (NTRS)

Beatty, T. D.

1975-01-01

A theoretical method is presented for the computation of the flow field about an axisymmetric body operating in a viscous, incompressible fluid. A potential flow method was used to determine the inviscid flow field and to yield the boundary conditions for the boundary layer solutions. Boundary layer effects in the forces of displacement thickness and empirically modeled separation streamlines are accounted for in subsequent potential flow solutions. This procedure is repeated until the solutions converge. An empirical method was used to determine base drag allowing configuration drag to be computed.

Design criterion for the heat-transfer coefficient in opposing flow, mixed convention heat transfer in a vertical tube

DOE Office of Scientific and Technical Information (OSTI.GOV)

Joye, D.D.

1996-07-01

Mixed convection heat transfer in a vertical tube with opposing flow (downflow heating) was studied experimentally for Reynolds numbers ranging from about 1,000 to 30,000 at constant Grashof numbers ranging about 1{1/2} orders of magnitude under constant wall temperature (CWT) conditions. Three correlations developed for opposing mixed convection flows in vertical conduits predicted the data reasonably well, except near and into the asymptote region for which these equations were not designed. A critical Reynolds number is developed here, above which these equations can be used for design purposes regardless of the boundary condition. Below Re{sub crit}, the correlations, the asymptotemore » equation should be used for the CWT boundary condition, which is more prevalent in process situations than the uniform heat flux (UHF) boundary condition.« less

A Global Interpolation Function (GIF) boundary element code for viscous flows

NASA Technical Reports Server (NTRS)

Reddy, D. R.; Lafe, O.; Cheng, A. H-D.

1995-01-01

Using global interpolation functions (GIF's), boundary element solutions are obtained for two- and three-dimensional viscous flows. The solution is obtained in the form of a boundary integral plus a series of global basis functions. The unknown coefficients of the GIF's are determined to ensure the satisfaction of the governing equations at selected collocation points. The values of the coefficients involved in the boundary integral equations are determined by enforcing the boundary conditions. Both primitive variable and vorticity-velocity formulations are examined.

The response of a laminar boundary layer in supersonic flow to small amplitude progressive waves

NASA Technical Reports Server (NTRS)

Duck, Peter W.

1989-01-01

The effect of a small amplitude progressive wave on the laminar boundary layer on a semi-infinite flat plate, due to a uniform supersonic freestream flow, is considered. The perturbation to the flow divides into two streamwise zones. In the first, relatively close to the leading edge of the plate, on a transverse scale comparable to the boundary layer thickness, the perturbation flow is described by a form of the unsteady linearized compressible boundary layer equations. In the freestream, this component of flow is governed by the wave equation, the solution of which provides the outer velocity conditions for the boundary layer. This system is solved numerically, and also the asymptotic structure in the far downstream limit is studied. This reveals a breakdown and a subsequent second streamwise zone, where the flow disturbance is predominantly inviscid. The two zones are shown to match in a proper asymptotic sense.

Wave dynamics in an extended macroscopic traffic flow model with periodic boundaries

NASA Astrophysics Data System (ADS)

Wang, Yu-Qing; Chu, Xing-Jian; Zhou, Chao-Fan; Yan, Bo-Wen; Jia, Bin; Fang, Chen-Hao

2018-06-01

Motivated by the previous traffic flow model considering the real-time traffic state, a modified macroscopic traffic flow model is established. The periodic boundary condition is applied to the car-following model. Besides, the traffic state factor R is defined in order to correct the real traffic conditions in a more reasonable way. It is a key step that we introduce the relaxation time as a density-dependent function and provide corresponding evolvement of traffic flow. Three different typical initial densities, namely the high density, the medium one and the low one, are intensively investigated. It can be found that the hysteresis loop exists in the proposed periodic-boundary system. Furthermore, the linear and nonlinear stability analyses are performed in order to test the robustness of the system.

An experimental study of the turbulent boundary layer on a transport wing in subsonic and transonic flow

NASA Technical Reports Server (NTRS)

Spaid, Frank W.; Roos, Frederick W.; Hicks, Raymond M.

1990-01-01

The upper surface boundary layer on a transport wing model was extensively surveyed with miniature yaw probes at a subsonic and a transonic cruise condition. Additional data were obtained at a second transonic test condition, for which a separated region was present at mid-semispan, aft of mid-chord. Significant variation in flow direction with distance from the surface was observed near the trailing edge except at the wing root and tip. The data collected at the transonic cruise condition show boundary layer growth associated with shock wave/boundary layer interaction, followed by recovery of the boundary layer downstream of the shock. Measurements of fluctuating surface pressure and wingtip acceleration were also obtained. The influence of flow field unsteadiness on the boundary layer data is discussed. Comparisons among the data and predictions from a variety of computational methods are presented. The computed predictions are in reasonable agreement with the experimental data in the outboard regions where 3-D effects are moderate and adverse pressure gradients are mild. In the more highly loaded mid-span region near the trailing edge, displacement thickness growth was significantly underpredicted, except when unrealistically severe adverse pressure gradients associated with inviscid calculations were used to perform boundary layer calculations.

Effects of inlet boundary conditions, on the computed flow in the Turbine-99 draft tube, using OpenFOAM and CFX

NASA Astrophysics Data System (ADS)

Nilsson, H.; Cervantes, M. J.

2012-11-01

The flow in the Turbine-99 Kaplan draft tube was thoroughly investigated at three workshops (1999, 2001, 2005), which aimed at determining the state of the art of draft tube simulations. The flow is challenging due to the different flow phenomena appearing simultaneously such as unsteadiness, separation, swirl, turbulence, and a strong adverse pressure gradient. The geometry and the experimentally determined inlet boundary conditions were provided to the Turbine-99 workshop participants. At the final workshop, angular resolved inlet velocity boundary conditions were provided. The rotating non-axi-symmetry of the inlet flow due to the runner blades was thus included. The effect of the rotating angular resolution was however not fully investigated at that workshop. The first purpose of this work is to further investigate this effect. Several different inlet boundary conditions are applied - the angular resolved experimental data distributed at the Turbine-99 workshop, the angular resolved results of a runner simulation with interpolated values using different resolution in the tangential and radial directions, and an axi-symmetric variant of the same numerical data. The second purpose of this work is to compare the results from the OpenFOAM and CFX CFD codes, using as similar settings as possible. The present results suggest that the experimental angular inlet boundary conditions proposed to the workshop are not adequate to simulate accurately the flow in the T-99 draft tube. The reason for this is that the experimental phase-averaged data has some important differences compared to the previously measured time-averaged data. Using the interpolated data from the runner simulation as inlet boundary condition however gives good results as long as the resolution of that data is sufficient. It is shown that the difference between the results using the angular-resolved and the corresponding symmetric inlet data is very small, suggesting that the importance of the angular resolution is small. The results from OpenFOAM and CFX are very similar as long as the inlet data resolution is fine enough. CFX seems to be more sensitive to that resolution.

Simulation of turbulent shear flows at Stanford and NASA-Ames - What can we do and what have we learned?

NASA Technical Reports Server (NTRS)

Reynolds, W. C.

1983-01-01

The capabilities and limitations of large eddy simulation (LES) and full turbulence simulation (FTS) are outlined. It is pointed out that LES, although limited at the present time by the need for periodic boundary conditions, produces large-scale flow behavior in general agreement with experiments. What is more, FTS computations produce small-scale behavior that is consistent with available experiments. The importance of the development work being done on the National Aerodynamic Simulator is emphasized. Studies at present are limited to situations in which periodic boundary conditions can be applied on boundaries of the computational domain where the flow is turbulent.

Step scaling and the Yang-Mills gradient flow

NASA Astrophysics Data System (ADS)

Lüscher, Martin

2014-06-01

The use of the Yang-Mills gradient flow in step-scaling studies of lattice QCD is expected to lead to results of unprecedented precision. Step scaling is usually based on the Schrödinger functional, where time ranges over an interval [0 , T] and all fields satisfy Dirichlet boundary conditions at time 0 and T. In these calculations, potentially important sources of systematic errors are boundary lattice effects and the infamous topology-freezing problem. The latter is here shown to be absent if Neumann instead of Dirichlet boundary conditions are imposed on the gauge field at time 0. Moreover, the expectation values of gauge-invariant local fields at positive flow time (and of other well localized observables) that reside in the center of the space-time volume are found to be largely insensitive to the boundary lattice effects.

Calibration of Axisymmetric and Quasi-1D Solvers for High Enthalpy Nozzles

NASA Technical Reports Server (NTRS)

Papadopoulos, P. E.; Gochberg, L. A.; Tokarcik-Polsky, S.; Venkatapathy, E.; Deiwert, G. S.; Edwards, Thomas A. (Technical Monitor)

1994-01-01

The proposed paper will present a numerical investigation of the flow characteristics and boundary layer development in the nozzles of high enthalpy shock tunnel facilities used for hypersonic propulsion testing. The computed flow will be validated against existing experimental data. Pitot pressure data obtained at the entrance of the test cabin will be used to validate the numerical simulations. It is necessary to accurately model the facility nozzles in order to characterize the test article flow conditions. Initially the axisymmetric nozzle flow will be computed using a Navier Stokes solver for a range of reservoir conditions. The calculated solutions will be compared and calibrated against available experimental data from the DLR HEG piston-driven shock tunnel and the 16-inch shock tunnel at NASA Ames Research Center. The Reynolds number is assumed to be high enough at the throat that the boundary layer flow is assumed turbulent at this point downstream. The real gas affects will be examined. In high Mach number facilities the boundary layer is thick. Attempts will be made to correlate the boundary layer displacement thickness. The displacement thickness correlation will be used to calibrate the quasi-1D codes NENZF and LSENS in order to provide fast and efficient tools of characterizing the facility nozzles. The calibrated quasi-1D codes will be implemented to study the effects of chemistry and the flow condition variations at the test section due to small variations in the driver gas conditions.

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Silva, Goncalo, E-mail: goncalo.nuno.silva@gmail.com; Talon, Laurent, E-mail: talon@fast.u-psud.fr; Ginzburg, Irina, E-mail: irina.ginzburg@irstea.fr

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies themore » effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.« less

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

NASA Astrophysics Data System (ADS)

Silva, Goncalo; Talon, Laurent; Ginzburg, Irina

2017-04-01

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes-Brinkman-Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes-Brinkman-Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

Unified solver for fluid dynamics and aeroacoustics in isentropic gas flows

NASA Astrophysics Data System (ADS)

Pont, Arnau; Codina, Ramon; Baiges, Joan; Guasch, Oriol

2018-06-01

The high computational cost of solving numerically the fully compressible Navier-Stokes equations, together with the poor performance of most numerical formulations for compressible flow in the low Mach number regime, has led to the necessity for more affordable numerical models for Computational Aeroacoustics. For low Mach number subsonic flows with neither shocks nor thermal coupling, both flow dynamics and wave propagation can be considered isentropic. Therefore, a joint isentropic formulation for flow and aeroacoustics can be devised which avoids the need for segregating flow and acoustic scales. Under these assumptions density and pressure fluctuations are directly proportional, and a two field velocity-pressure compressible formulation can be derived as an extension of an incompressible solver. Moreover, the linear system of equations which arises from the proposed isentropic formulation is better conditioned than the homologous incompressible one due to the presence of a pressure time derivative. Similarly to other compressible formulations the prescription of boundary conditions will have to deal with the backscattering of acoustic waves. In this sense, a separated imposition of boundary conditions for flow and acoustic scales which allows the evacuation of waves through Dirichlet boundaries without using any tailored damping model will be presented.

Developments in flow visualization methods for flight research

NASA Technical Reports Server (NTRS)

Holmes, Bruce J.; Obara, Clifford J.; Manuel, Gregory S.; Lee, Cynthia C.

1990-01-01

With the introduction of modern airplanes utilizing laminar flow, flow visualization has become an important diagnostic tool in determining aerodynamic characteristics such as surface flow direction and boundary-layer state. A refinement of the sublimating chemical technique has been developed to define both the boundary-layer transition location and the transition mode. In response to the need for flow visualization at subsonic and transonic speeds and altitudes above 20,000 feet, the liquid crystal technique has been developed. A third flow visualization technique that has been used is infrared imaging, which offers non-intrusive testing over a wide range of test conditions. A review of these flow visualization methods and recent flight results is presented for a variety of modern aircraft and flight conditions.

Absorbing boundary conditions for second-order hyperbolic equations

NASA Technical Reports Server (NTRS)

Jiang, Hong; Wong, Yau Shu

1989-01-01

A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.

On flows of viscoelastic fluids under threshold-slip boundary conditions

NASA Astrophysics Data System (ADS)

Baranovskii, E. S.

2018-03-01

We investigate a boundary-value problem for the steady isothermal flow of an incompressible viscoelastic fluid of Oldroyd type in a 3D bounded domain with impermeable walls. We use the Fujita threshold-slip boundary condition. This condition states that the fluid can slip along a solid surface when the shear stresses reach a certain critical value; otherwise the slipping velocity is zero. Assuming that the flow domain is not rotationally symmetric, we prove an existence theorem for the corresponding slip problem in the framework of weak solutions. The proof uses methods for solving variational inequalities with pseudo-monotone operators and convex functionals, the method of introduction of auxiliary viscosity, as well as a passage-to-limit procedure based on energy estimates of approximate solutions, Korn’s inequality, and compactness arguments. Also, some properties and estimates of weak solutions are established.

Study of the effect of wind speed on evaporation from soil through integrated modeling of atmospheric boundary layer and shallow subsurface

NASA Astrophysics Data System (ADS)

Davarzani, Hossein; Smits, Kathleen; Tolene, Ryan; Illangasekare, Tissa

2013-04-01

The study of the interaction between the land and atmosphere is paramount to our understanding of many emerging problems to include climate change, the movement of green house gases such as possible leaking of sequestered CO2 and the accurate detection of buried objects such as landmines. Soil moisture distribution in the shallow subsurface becomes a critical factor in all these problems. The heat and mass flux in the form of soil evaporation across the land surface couples the atmospheric boundary layer to the shallow subsurface. The coupling between land and the atmosphere leads to highly dynamic interactions between the porous media properties, transport processes and boundary conditions, resulting in dynamic evaporative behavior. However, the coupling at the land-atmospheric interface is rarely considered in most current models and their validation for practical applications. This is due to the complexity of the problem in field scenarios and the scarcity of field or laboratory data capable of testing and refining coupled energy and mass transfer theories. In most efforts to compute evaporation from soil, only indirect coupling is provided to characterize the interaction between non-isothermal multiphase flows under realistic atmospheric conditions even though heat and mass flux are controlled by the coupled dynamics of the land and the atmospheric boundary layer. In earlier drying modeling concepts, imposing evaporation flux (kinetic of relative humidity) and temperature as surface boundary condition is often needed. With the goal of improving our understanding of the land/atmospheric coupling, we developed a model based on the coupling of Navier-Stokes free flow and Darcy flow in porous medium. The model consists of the coupled equations of mass conservation for the liquid phase (water) and gas phase (water vapor and air) in porous medium with gas phase (water vapor and air) in free flow domain under non-isothermal, non-equilibrium conditions. The boundary conditions at the porous medium-free flow medium interface include dynamical, thermal and solutal equilibriums, and using the Beavers-Joseph slip boundary condition. What is unique about this model is that the evaporation rate and soil surface temperature conditions come directly from the model output. In order to experimentally validate the numerical results, we developed and used a unique two dimensional wind tunnel placed above a soil tank equipped with a network of different sensors. A series of experiments under varying boundary conditions, using a test sand for which the hydraulic and thermal properties were well characterized, were performed. Precision data for soil moisture, soil and air temperature and relative humidity, and also wind velocity under well-controlled transient heat and wind boundary conditions was generated. Results from numerical simulations were compared with experimental data. Results demonstrate that the coupling concept can predict the different stages of the drying process in porous media with good accuracy. Increasing the wind speed increases the first stage evaporation rate and decreases the transition time at low velocity values; then, at high values of wind speed the evaporation rate becomes less dependent of flow in free fluid. In the opposite, the impact of the wind speed on the second stage evaporation (diffusion dominant stage) is not significant. The proposed theoretical model can be used to predict the evaporation process where a porous medium flow is coupled to a free flow for different practical applications.

Inverse boundary-layer theory and comparison with experiment

NASA Technical Reports Server (NTRS)

Carter, J. E.

1978-01-01

Inverse boundary layer computational procedures, which permit nonsingular solutions at separation and reattachment, are presented. In the first technique, which is for incompressible flow, the displacement thickness is prescribed; in the second technique, for compressible flow, a perturbation mass flow is the prescribed condition. The pressure is deduced implicitly along with the solution in each of these techniques. Laminar and turbulent computations, which are typical of separated flow, are presented and comparisons are made with experimental data. In both inverse procedures, finite difference techniques are used along with Newton iteration. The resulting procedure is no more complicated than conventional boundary layer computations. These separated boundary layer techniques appear to be well suited for complete viscous-inviscid interaction computations.

An Implementation of Hydrostatic Boundary Conditions for Variable Density Lattice Boltzmann Methods

NASA Astrophysics Data System (ADS)

Bardsley, K. J.; Thorne, D. T.; Lee, J. S.; Sukop, M. C.

2006-12-01

Lattice Boltzmann Methods (LBMs) have been under development for the last two decades and have become another capable numerical method for simulating fluid flow. Recent advances in lattice Boltzmann applications involve simulation of density-dependent fluid flow in closed (Dixit and Babu, 2006; D'Orazio et al., 2004) or periodic (Guo and Zhao, 2005) domains. However, standard pressure boundary conditions (BCs) are incompatible with concentration-dependent density flow simulations that use a body force for gravity. An implementation of hydrostatic BCs for use under these conditions is proposed here. The basis of this new implementation is an additional term in the pressure BC. It is derived to account for the incorporation of gravity as a body force and the effect of varying concentration in the fluid. The hydrostatic BC expands the potential of density-dependent LBM to simulate domains with boundaries other than the closed or periodic boundaries that have appeared in previous literature on LBM simulations. With this new implementation, LBM will be able to simulate complex concentration-dependent density flows, such as salt water intrusion in the classic Henry and Henry-Hilleke problems. This is demonstrated using various examples, beginning with a closed box system, and ending with a system containing two solid walls, one velocity boundary and one pressure boundary, as in the Henry problem. References Dixit, H. N., V. Babu, (2006), Simulation of high Rayleigh number natural convection in a square cavity using the lattice Boltzmann method, Int. J. Heat Mass Transfer, 49, 727-739. D'Orazio, A., M. Corcione, G.P. Celata, (2004), Application to natural convection enclosed flows of a lattice Boltzmann BGK model coupled with a general purpose thermal boundary conditions, Int. J. Thermal Sci., 43, 575-586. Gou, Z., T.S. Zhao, (2005), Lattice Boltzmann simulation of natural convection with temperature-dependant viscosity in a porous cavity, Numerical Heat Transfer, Part B, 47, 157-177.

Analytical method for predicting the pressure distribution about a nacelle at transonic speeds

NASA Technical Reports Server (NTRS)

Keith, J. S.; Ferguson, D. R.; Merkle, C. L.; Heck, P. H.; Lahti, D. J.

1973-01-01

The formulation and development of a computer analysis for the calculation of streamlines and pressure distributions around two-dimensional (planar and axisymmetric) isolated nacelles at transonic speeds are described. The computerized flow field analysis is designed to predict the transonic flow around long and short high-bypass-ratio fan duct nacelles with inlet flows and with exhaust flows having appropriate aerothermodynamic properties. The flow field boundaries are located as far upstream and downstream as necessary to obtain minimum disturbances at the boundary. The far-field lateral flow field boundary is analytically defined to exactly represent free-flight conditions or solid wind tunnel wall effects. The inviscid solution technique is based on a Streamtube Curvature Analysis. The computer program utilizes an automatic grid refinement procedure and solves the flow field equations with a matrix relaxation technique. The boundary layer displacement effects and the onset of turbulent separation are included, based on the compressible turbulent boundary layer solution method of Stratford and Beavers and on the turbulent separation prediction method of Stratford.

Pressure wave propagation studies for oscillating cascades

NASA Technical Reports Server (NTRS)

Huff, Dennis L.

1992-01-01

The unsteady flow field around an oscillating cascade of flat plates is studied using a time marching Euler code. Exact solutions based on linear theory serve as model problems to study pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step is demonstrated. Results show that an approximate non-reflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer than cases using reflective boundary conditions. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates does a better job of capturing the unsteady flow field than cases using uniform grids as long as the CFL number is less than one for a sufficient portion of the grid. Results for various stagger angles and oscillation frequencies show good agreement with linear theory as long as the grid is properly resolved.

Friction-term response to boundary-condition type in flow models

USGS Publications Warehouse

Schaffranek, R.W.; Lai, C.

1996-01-01

The friction-slope term in the unsteady open-channel flow equations is examined using two numerical models based on different formulations of the governing equations and employing different solution methods. The purposes of the study are to analyze, evaluate, and demonstrate the behavior of the term in a set of controlled numerical experiments using varied types and combinations of boundary conditions. Results of numerical experiments illustrate that a given model can respond inconsistently for the identical resistance-coefficient value under different types and combinations of boundary conditions. Findings also demonstrate that two models employing different dependent variables and solution methods can respond similarly for the identical resistance-coefficient value under similar types and combinations of boundary conditions. Discussion of qualitative considerations and quantitative experimental results provides insight into the proper treatment, evaluation, and significance of the friction-slope term, thereby offering practical guidelines for model implementation and calibration.

Scale effect of slip boundary condition at solid–liquid interface

PubMed Central

Nagayama, Gyoko; Matsumoto, Takenori; Fukushima, Kohei; Tsuruta, Takaharu

2017-01-01

Rapid advances in microelectromechanical systems have stimulated the development of compact devices, which require effective cooling technologies (e.g., microchannel cooling). However, the inconsistencies between experimental and classical theoretical predictions for the liquid flow in microchannel remain unclarified. Given the larger surface/volume ratio of microchannel, the surface effects increase as channel scale decreases. Here we show the scale effect of the boundary condition at the solid–liquid interface on single-phase convective heat transfer characteristics in microchannels. We demonstrate that the deviation from classical theory with a reduction in hydraulic diameters is due to the breakdown of the continuum solid–liquid boundary condition. The forced convective heat transfer characteristics of single-phase laminar flow in a parallel-plate microchannel are investigated. Using the theoretical Poiseuille and Nusselt numbers derived under the slip boundary condition at the solid–liquid interface, we estimate the slip length and thermal slip length at the interface. PMID:28256536

Calculations of unsteady turbulent boundary layers with flow reversal

NASA Technical Reports Server (NTRS)

Nash, J. F.; Patel, V. C.

1975-01-01

The results are presented of a series of computational experiments aimed at studying the characteristics of time-dependent turbulent boundary layers with embedded reversed-flow regions. A calculation method developed earlier was extended to boundary layers with reversed flows for this purpose. The calculations were performed for an idealized family of external velocity distributions, and covered a range of degrees of unsteadiness. The results confirmed those of previous studies in demonstrating that the point of flow reversal is nonsingular in a time-dependent boundary layer. A singularity was observed to develop downstream of reversal, under certain conditions, accompanied by the breakdown of the boundary-layer approximations. A tentative hypothesis was advanced in an attempt to predict the appearance of the singularity, and is shown to be consistent with the calculated results.

MHD Free Convective Boundary Layer Flow of a Nanofluid past a Flat Vertical Plate with Newtonian Heating Boundary Condition

PubMed Central

Uddin, Mohammed J.; Khan, Waqar A.; Ismail, Ahmed I.

2012-01-01

Steady two dimensional MHD laminar free convective boundary layer flows of an electrically conducting Newtonian nanofluid over a solid stationary vertical plate in a quiescent fluid taking into account the Newtonian heating boundary condition is investigated numerically. A magnetic field can be used to control the motion of an electrically conducting fluid in micro/nano scale systems used for transportation of fluid. The transport equations along with the boundary conditions are first converted into dimensionless form and then using linear group of transformations, the similarity governing equations are developed. The transformed equations are solved numerically using the Runge-Kutta-Fehlberg fourth-fifth order method with shooting technique. The effects of different controlling parameters, namely, Lewis number, Prandtl number, buoyancy ratio, thermophoresis, Brownian motion, magnetic field and Newtonian heating on the flow and heat transfer are investigated. The numerical results for the dimensionless axial velocity, temperature and nanoparticle volume fraction as well as the reduced Nusselt and Sherwood number have been presented graphically and discussed. It is found that the rate of heat and mass transfer increase as Newtonian heating parameter increases. The dimensionless velocity and temperature distributions increase with the increase of Newtonian heating parameter. The results of the reduced heat transfer rate is compared for convective heating boundary condition and found an excellent agreement. PMID:23166688

Sound-turbulence interaction in transonic boundary layers

NASA Astrophysics Data System (ADS)

Lelostec, Ludovic; Scalo, Carlo; Lele, Sanjiva

2014-11-01

Acoustic wave scattering in a transonic boundary layer is investigated through a novel approach. Instead of simulating directly the interaction of an incoming oblique acoustic wave with a turbulent boundary layer, suitable Dirichlet conditions are imposed at the wall to reproduce only the reflected wave resulting from the interaction of the incident wave with the boundary layer. The method is first validated using the laminar boundary layer profiles in a parallel flow approximation. For this scattering problem an exact inviscid solution can be found in the frequency domain which requires numerical solution of an ODE. The Dirichlet conditions are imposed in a high-fidelity unstructured compressible flow solver for Large Eddy Simulation (LES), CharLESx. The acoustic field of the reflected wave is then solved and the interaction between the boundary layer and sound scattering can be studied.

Simulation of ground-water flow in coastal Georgia and adjacent parts of South Carolina and Florida-predevelopment, 1980, and 2000

USGS Publications Warehouse

Payne, Dorothy F.; Rumman, Malek Abu; Clarke, John S.

2005-01-01

A digital model was developed to simulate steady-state ground-water flow in a 42,155-square-mile area of coastal Georgia and adjacent parts of South Carolina and Florida. The model was developed to (1) understand and refine the conceptual model of regional ground-water flow, (2) serve as a framework for the development of digital subregional ground-water flow and solute-transport models, and (3) serve as a tool for future evaluations of hypothetical pumping scenarios used to facilitate water management in the coastal area. Single-density ground-water flow was simulated using the U.S. Geological Survey finite-difference code MODFLOW-2000 for mean-annual conditions during predevelopment (pre?1900) and the years 1980 and 2000. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units. A combination of boundary conditions was applied, including a general-head boundary condition on the top active cells of the model and a time-variable fixed-head boundary condition along part of the southern lateral boundary. Simulated heads for 1980 and 2000 conditions indicate a good match to observed values, based on a plus-or-minus 10-foot (ft) calibration target and calibration statistics. The root-mean square of residual water levels for the Upper Floridan aquifer was 13.0 ft for the 1980 calibration and 9.94 ft for the 2000 calibration. Some spatial patterns of residuals were indicated for the 1980 and 2000 simulations, and are likely a result of model-grid cell size and insufficiently detailed hydraulic-property and pumpage data in some areas. Simulated potentiometric surfaces for predevelopment, 1980, and 2000 conditions all show major flow system features that are indicated by estimated peotentiometric maps. During 1980?2000, simulated water levels at the centers of pumping at Savannah and Brunswick rose more than 20 ft and 8 ft, respectively, in response to decreased pumping. Simulated drawdown exceeded 10 ft in the Upper Floridan aquifer across much of the western half of the model area, with drawdown exceeding 20 ft along parts of the western, northern, and southern boundaries where irrigation pumping increased during this period. From predevelopment to 2000 conditions, the simulated water budget showed an increase in inflow from, and decrease in outflow to, the general-head boundaries, and a reversal from net seaward flow to net landward flow across the coastline. Simulated changes in recharge and discharge distribution from predevelopment to 2000 conditions showed an increase in extent and magnitude of net recharge cells in the northern part of the model area, and a decrease in discharge or change to recharge in cells containing major streams and beneath major pumping centers. The model is relatively sensitive to pumping and the controlling head at the fixed-head boundary and less sensitive to the distribution of aquifer properties in general. Model limitations include: (1) its spatial scale and discretization, (2) the extent to which data are available to physically define the flow system, (3) the type of boundary conditions and controlling parameters used, (4) uncertainty in the distribution of pumping, and (5) uncertainty in field-scale hydraulic properties. The model could be improved with more accurate estimates of ground-water pumpage and better characterization of recharge and discharge.

A non-local computational boundary condition for duct acoustics

NASA Technical Reports Server (NTRS)

Zorumski, William E.; Watson, Willie R.; Hodge, Steve L.

1994-01-01

A non-local boundary condition is formulated for acoustic waves in ducts without flow. The ducts are two dimensional with constant area, but with variable impedance wall lining. Extension of the formulation to three dimensional and variable area ducts is straightforward in principle, but requires significantly more computation. The boundary condition simulates a nonreflecting wave field in an infinite duct. It is implemented by a constant matrix operator which is applied at the boundary of the computational domain. An efficient computational solution scheme is developed which allows calculations for high frequencies and long duct lengths. This computational solution utilizes the boundary condition to limit the computational space while preserving the radiation boundary condition. The boundary condition is tested for several sources. It is demonstrated that the boundary condition can be applied close to the sound sources, rendering the computational domain small. Computational solutions with the new non-local boundary condition are shown to be consistent with the known solutions for nonreflecting wavefields in an infinite uniform duct.

A computer program for calculating laminar and turbulent boundary layers for two-dimensional time-dependent flows

NASA Technical Reports Server (NTRS)

Cebeci, T.; Carr, L. W.

1978-01-01

A computer program is described which provides solutions of two dimensional equations appropriate to laminar and turbulent boundary layers for boundary conditions with an external flow which fluctuates in magnitude. The program is based on the numerical solution of the governing boundary layer equations by an efficient two point finite difference method. An eddy viscosity formulation was used to model the Reynolds shear stress term. The main features of the method are briefly described and instructions for the computer program with a listing are provided. Sample calculations to demonstrate its usage and capabilities for laminar and turbulent unsteady boundary layers with an external flow which fluctuated in magnitude are presented.

Generalized Forchheimer Flows of Isentropic Gases

NASA Astrophysics Data System (ADS)

Celik, Emine; Hoang, Luan; Kieu, Thinh

2018-03-01

We consider generalized Forchheimer flows of either isentropic gases or slightly compressible fluids in porous media. By using Muskat's and Ward's general form of the Forchheimer equations, we describe the fluid dynamics by a doubly nonlinear parabolic equation for the appropriately defined pseudo-pressure. The volumetric flux boundary condition is converted to a time-dependent Robin-type boundary condition for this pseudo-pressure. We study the corresponding initial boundary value problem, and estimate the L^∞ and W^{1,2-a} (with 0

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).

Topology optimization of thermal fluid flows with an adjoint Lattice Boltzmann Method

NASA Astrophysics Data System (ADS)

Dugast, Florian; Favennec, Yann; Josset, Christophe; Fan, Yilin; Luo, Lingai

2018-07-01

This paper presents an adjoint Lattice Boltzmann Method (LBM) coupled with the Level-Set Method (LSM) for topology optimization of thermal fluid flows. The adjoint-state formulation implies discrete velocity directions in order to take into account the LBM boundary conditions. These boundary conditions are introduced at the beginning of the adjoint-state method as the LBM residuals, so that the adjoint-state boundary conditions can appear directly during the adjoint-state equation formulation. The proposed method is tested with 3 numerical examples concerning thermal fluid flows, but with different objectives: minimization of the mean temperature in the domain, maximization of the heat evacuated by the fluid, and maximization of the heat exchange with heated solid parts. This latter example, treated in several articles, is used to validate our method. In these optimization problems, a limitation of the maximal pressure drop and of the porosity (number of fluid elements) is also applied. The obtained results demonstrate that the method is robust and effective for solving topology optimization of thermal fluid flows.

Complex Wall Boundary Conditions for Modeling Combustion in Catalytic Channels

NASA Astrophysics Data System (ADS)

Zhu, Huayang; Jackson, Gregory

2000-11-01

Monolith catalytic reactors for exothermic oxidation are being used in automobile exhaust clean-up and ultra-low emissions combustion systems. The reactors present a unique coupling between mass, heat, and momentum transport in a channel flow configuration. The use of porous catalytic coatings along the channel wall presents a complex boundary condition when modeled with the two-dimensional channel flow. This current work presents a 2-D transient model for predicting the performance of catalytic combustion systems for methane oxidation on Pd catalysts. The model solves the 2-D compressible transport equations for momentum, species, and energy, which are solved with a porous washcoat model for the wall boundary conditions. A time-splitting algorithm is used to separate the stiff chemical reactions from the convective/diffusive equations for the channel flow. A detailed surface chemistry mechanism is incorporated for the catalytic wall model and is used to predict transient ignition and steady-state conversion of CH4-air flows in the catalytic reactor.

On the definition of dominant force regimes for flow boiling heat transfer by using single mini-tubes

NASA Astrophysics Data System (ADS)

Baba, Soumei; Sawada, Kenichiro; Kubota, Chisato; Kawanami, Osamu; Asano, Hitoshi; Inoue, Koichi; Ohta, Haruhiko

Recent increase in the size of space platforms requires the management of larger amount of waste heat under high heat flux conditions and the transportation of it along a long distance to the radiator. Flow boiling applied to the thermal management system in space attracts much attention as promising means to realize high-performance heat transfer and transport because of large latent heat of vaporization. In microgravity two-phase flow phenomena are quite different from those under 1-g condition because buoyancy effects are significantly reduced and surface tension becomes dominant. By the similar reason, flow boiling characteristics in mini channels are not the same as those in channels of normal sizes. In the present stage, however, the boundary between the regimes of body force dominated and of surface tension dominated is not clear. The design of space thermal devices, operated under the conditions where no effect of gravity is expected, will improve the reliability of their ground tests, provided that the boundaries of dominant force regimes are clarified quantitatively in advance. In flow boiling in mini channels or in parallel channels, back flow could be occurred because of rapid growth of bubbles in a confined space, resulting flow rate fluctuation. Flow boiling heat transfer characteristics in mini channels can be changed considerably by the existence of inlet flow rate fluctuation. It is important to pay attention to experimental accuracy and to use a single circular mini-tube to compare heat transfer characteristics with those of normal size tubes. In the present paper, effects of tube orientations, i.e. vertical upward flow, vertical downward flow and horizontal flow, on flow boiling heat transfer characteristics is investigated for FC72 flowing in single mini-tubes with inner diameters of 0.13 and 0.51 mm to establish a reliable dominant force regime map. If the regime map is described by using dimensionless groups of Bond, Weber and Froude numbers, the boundary of dominant forces of inertia and body force is examined by using the mini-tube of the larger diameter at constant Bond number Bo = 0.51, and the boundary of inertia and surface tension by using the mini-tube of smaller diameter at Bo = 0.033. The influence of inertia is varied by the change of vapor quality, i.e. ratio of vapor mass flow rate to the total, under constant mass velocities, where the velocity of liquid-vapor mixtures is increased with increasing vapor quality. For the tube diameter of 0.51 mm, under low inertia conditions at Froude number Fr < 5, heat transfer coefficients were influenced by the tube orientation, while the heat transfer coefficients were almost independent of the orientation for Fr > 5. The results indicated that the boundary between the body force dominated and the inertia force dominated regimes was given by Froude number as Fr ˜ 5. On the other hand, for tube diameter of 0.13 mm, almost no effect of tube = orientation was observed for all combinations of mass velocity and vapor quality, and heat transfer coefficients were independent of vapor quality under low inertia conditions at Weber number We < 5, and vice versa. The results implied the boundary between the surface tension dominated and the inertia force dominated regimes was represented by We ˜ 5. = In addition, by the reflection of both results on the two-dimensional regime map, the boundary between the surface tension dominated and the body force dominated regimes was approx-imately evaluated as Bo ˜ 0.25 from the crossing point of two boundary lines. This value = located in the range of 0.033 < Bo < 0.51 is consistent with the boundaries between the sur-face tension dominated and the body force dominated regimes classified for the smaller and larger mini-tubes, respectively, under low inertia conditions.

Flow separation of currents in shallow water

USGS Publications Warehouse

Signell, Richard P.

1989-01-01

Flow separation of currents in shallow coastal areas is investigated using a boundary layer model for two-dimensional (depth-averaged) tidal flow past an elliptic headland. If the shoaling region near the coast is narrow compared to the scale of the headland, bottom friction causes the flow to separate just downstream of the point where the pressure gradient switches from favoring to adverse. As long as the shoaling region at the coast is well resolved, the inclusion of eddy viscosity and a no-slip boundary condition have no effect on this result. An approximate analytic solution for the pressure gradient along the boundary is obtained by assuming the flow away from the immediate vicinity of the boundary is irrotational. On the basis of the pressure gradient obtained from the irrotational flow solution, flow separation is a strong function of the headland aspect ratio, an equivalent Reynolds number, and a Keulegan-Carpenter number.

Stokesian dynamics of pill-shaped Janus particles with stick and slip boundary conditions

NASA Astrophysics Data System (ADS)

Sun, Qiang; Klaseboer, Evert; Khoo, Boo Cheong; Chan, Derek Y. C.

2013-04-01

We study the forces and torques experienced by pill-shaped Janus particles of different aspect ratios where half of the surface obeys the no-slip boundary condition and the other half obeys the Navier slip condition of varying slip lengths. Using a recently developed boundary integral formulation whereby the traditional singular behavior of this approach is removed analytically, we quantify the strength of the forces and torques experienced by such particles in a uniform flow field in the Stokes regime. Depending on the aspect ratio and the slip length, the force transverse to the flow direction can change sign. This is a novel property unique to the Janus nature of the particles.

Nonequilibrium viscous flow over Jovian entry probes at high altitudes

NASA Technical Reports Server (NTRS)

Kumar, A.; Szema, K. Y.; Tiwari, S. N.

1979-01-01

The viscous chemical nonequilibrium flow around a Jovian entry body is investigated at high altitudes using two different methods. First method is only for the stagnation region and integrates the full Navier-Stokes equations from the body surface to the freestream. The second method uses viscous shock layer equations between the body surface and the shock. Due to low Reynolds numbers, both methods use surface slip boundary conditions and the second method also uses shock slip boundary conditions. The results of the two methods are compared at the stagnation point. It is found that the entire shock layer is under chemical nonequilibrium at higher altitudes and that the slip boundary conditions are important at these altitudes.

Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices

NASA Technical Reports Server (NTRS)

Nance, Robert P.; Hash, David B.; Hassan, H. A.

1997-01-01

A study is made of the issues surrounding prediction of microchannel flows using the direct simulation Monte Carlo method. This investigation includes the introduction and use of new inflow and outflow boundary conditions suitable for subsonic flows. A series of test simulations for a moderate-size microchannel indicates that a high degree of grid under-resolution in the streamwise direction may be tolerated without loss of accuracy. In addition, the results demonstrate the importance of physically correct boundary conditions, as well as possibilities for reducing the time associated with the transient phase of a simulation. These results imply that simulations of longer ducts may be more feasible than previously envisioned.

A level-set method for two-phase flows with moving contact line and insoluble surfactant

NASA Astrophysics Data System (ADS)

Xu, Jian-Jun; Ren, Weiqing

2014-04-01

A level-set method for two-phase flows with moving contact line and insoluble surfactant is presented. The mathematical model consists of the Navier-Stokes equation for the flow field, a convection-diffusion equation for the surfactant concentration, together with the Navier boundary condition and a condition for the dynamic contact angle derived by Ren et al. (2010) [37]. The numerical method is based on the level-set continuum surface force method for two-phase flows with surfactant developed by Xu et al. (2012) [54] with some cautious treatment for the boundary conditions. The numerical method consists of three components: a flow solver for the velocity field, a solver for the surfactant concentration, and a solver for the level-set function. In the flow solver, the surface force is dealt with using the continuum surface force model. The unbalanced Young stress at the moving contact line is incorporated into the Navier boundary condition. A convergence study of the numerical method and a parametric study are presented. The influence of surfactant on the dynamics of the moving contact line is illustrated using examples. The capability of the level-set method to handle complex geometries is demonstrated by simulating a pendant drop detaching from a wall under gravity.

Thin-layer approximation and algebraic model for separated turbulent flows

NASA Technical Reports Server (NTRS)

Baldwin, B.; Lomax, H.

1978-01-01

An algebraic turbulence model for two- and three-dimensional separated flows is specified that avoids the necessity for finding the edge of the boundary layer. Properties of the model are determined and comparisons made with experiment for an incident shock on a flat plate, separated flow over a compression corner, and transonic flow over an airfoil. Separation and reattachment points from numerical Navier-Stokes solutions agree with experiment within one boundary-layer thickness. Use of law-of-the-wall boundary conditions does not alter the predictions significantly. Applications of the model to other cases are contained in companion papers.

Numerical Investigation of Flow in an Over-Expanded Nozzle with Porous Surfaces

NASA Technical Reports Server (NTRS)

Elmiligui, Alaa; Abdol-Hamid, K. S.; Hunter, Craig A.

2005-01-01

A new porous condition has been implemented in the PAB3D solver for simulating the flow over porous surfaces. The newly-added boundary condition is utilized to compute the flow field of a non-axisymmetric, convergent-divergent nozzle incorporating porous cavities for shock-boundary layer interaction control. The nozzle has an expansion ratio (exit area/throat area) of 1.797 and a design nozzle pressure ratio of 8.78. The flow fields for a baseline nozzle (no porosity) and for a nozzle with porous surfaces (10% porosity ratio) are computed for NPR varying from 2.01 to 9.54. Computational model results indicate that the over-expanded nozzle flow was dominated by shock-induced boundary-layer separation. Porous configurations were capable of controlling off-design separation in the nozzle by encouraging stable separation of the exhaust flow. Computational simulation results, wall centerline pressure, mach contours, and thrust efficiency ratio are presented and discussed. Computed results are in excellent agreement with experimental data.

Numerical Investigation of Flow in an Over-expanded Nozzle with Porous Surfaces

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.; Elmilingui, Alaa A.; Hunter, Craig A.

2006-01-01

A new porous condition has been implemented in the PAB3D solver for simulating the flow over porous surfaces. The newly-added boundary condition is utilized to compute the flow field of a non-axisymmetric, convergent-divergent nozzle incorporating porous cavities for shock-boundary layer interaction control. The nozzle has an expansion ratio (exit area/throat area) of 1.797 and a design nozzle pressure ratio of 8.78. The flow fields for a baseline nozzle (no porosity) and for a nozzle with porous surfaces (10% porosity ratio) are computed for NPR varying from 2.01 to 9.54. Computational model results indicate that the over-expanded nozzle flow is dominated by shock-induced boundary-layer separation. Porous configurations are capable of controlling off-design separation in the nozzle by encouraging stable separation of the exhaust flow. Computational simulation results, wall centerline pressure, mach contours, and thrust efficiency ratio are presented and discussed. Computed results are in excellent agreement with experimental data.

Patient-specific computational modeling of blood flow in the pulmonary arterial circulation.

PubMed

Kheyfets, Vitaly O; Rios, Lourdes; Smith, Triston; Schroeder, Theodore; Mueller, Jeffrey; Murali, Srinivas; Lasorda, David; Zikos, Anthony; Spotti, Jennifer; Reilly, John J; Finol, Ender A

2015-07-01

Computational fluid dynamics (CFD) modeling of the pulmonary vasculature has the potential to reveal continuum metrics associated with the hemodynamic stress acting on the vascular endothelium. It is widely accepted that the endothelium responds to flow-induced stress by releasing vasoactive substances that can dilate and constrict blood vessels locally. The objectives of this study are to examine the extent of patient specificity required to obtain a significant association of CFD output metrics and clinical measures in models of the pulmonary arterial circulation, and to evaluate the potential correlation of wall shear stress (WSS) with established metrics indicative of right ventricular (RV) afterload in pulmonary hypertension (PH). Right Heart Catheterization (RHC) hemodynamic data and contrast-enhanced computed tomography (CT) imaging were retrospectively acquired for 10 PH patients and processed to simulate blood flow in the pulmonary arteries. While conducting CFD modeling of the reconstructed patient-specific vasculatures, we experimented with three different outflow boundary conditions to investigate the potential for using computationally derived spatially averaged wall shear stress (SAWSS) as a metric of RV afterload. SAWSS was correlated with both pulmonary vascular resistance (PVR) (R(2)=0.77, P<0.05) and arterial compliance (C) (R(2)=0.63, P<0.05), but the extent of the correlation was affected by the degree of patient specificity incorporated in the fluid flow boundary conditions. We found that decreasing the distal PVR alters the flow distribution and changes the local velocity profile in the distal vessels, thereby increasing the local WSS. Nevertheless, implementing generic outflow boundary conditions still resulted in statistically significant SAWSS correlations with respect to both metrics of RV afterload, suggesting that the CFD model could be executed without the need for complex outflow boundary conditions that require invasively obtained patient-specific data. A preliminary study investigating the relationship between outlet diameter and flow distribution in the pulmonary tree offers a potential computationally inexpensive alternative to pressure based outflow boundary conditions. Copyright © 2015 Elsevier Ireland Ltd. All rights reserved.

Effect of inlet conditions for numerical modelling of the urban boundary layer

NASA Astrophysics Data System (ADS)

Gnatowska, Renata

2018-01-01

The paper presents the numerical results obtained with the use of the ANSYS FLUENT commercial code for analysing the flow structure around two rectangular inline surface-mounted bluff bodies immersed in a boundary layer. The effects of the inflow boundary layer for the accuracy of the numerical modelling of the flow field around a simple system of objects are described. The analysis was performed for two concepts. In the former case, the inlet velocity profile was defined using the power law, whereas the kinetic and dissipation energy was defined from the equations according to Richards and Hoxey [1]. In the latter case, the inlet conditions were calculated for the flow over the rough area composed of the rectangular components.

Similarity solutions for unsteady free-convection flow from a continuous moving vertical surface

NASA Astrophysics Data System (ADS)

Abd-El-Malek, Mina B.; Kassem, Magda M.; Mekky, Mohammad L.

2004-03-01

The transformation group theoretic approach is applied to present an analysis of the problem of unsteady free convection flow over a continuous moving vertical sheet in an ambient fluid. The thermal boundary layer induced within a vertical semi-infinite layer of Boussinseq fluid by a constant heated bounding plate. The application of two-parameter groups reduces the number of independent variables by two, and consequently the system of governing partial differential equations with the boundary conditions reduces to a system of ordinary differential equations with appropriate boundary conditions. The obtained ordinary differential equations are solved analytically for the temperature and numerically for the velocity using the shooting method. Effect of Prandtl number on the thermal boundary-layer and velocity boundary-layer are studied and plotted in curves.

Effects of different boundary conditions on the simulation of groundwater flow in a multi-layered coastal aquifer system (Taranto Gulf, southern Italy)

NASA Astrophysics Data System (ADS)

De Filippis, Giovanna; Foglia, Laura; Giudici, Mauro; Mehl, Steffen; Margiotta, Stefano; Negri, Sergio L.

2017-11-01

The evaluation of the accuracy or reasonableness of numerical models of groundwater flow is a complex task, due to the uncertainties in hydrodynamic properties and boundary conditions and the scarcity of good-quality field data. To assess model reliability, different calibration techniques are joined to evaluate the effects of different kinds of boundary conditions on the groundwater flow in a coastal multi-layered aquifer in southern Italy. In particular, both direct and indirect approaches for inverse modeling were joined through the calibration of one of the most uncertain parameters, namely the hydraulic conductivity of the karst deep hydrostratigraphic unit. The methodology proposed here, and applied to a real case study, confirmed that the selection of boundary conditions is among the most critical and difficult aspects of the characterization of a groundwater system for conceptual analysis or numerical simulation. The practical tests conducted in this study show that incorrect specification of boundary conditions prevents an acceptable match between the model response to the hydraulic stresses and the behavior of the natural system. Such effects have a negative impact on the applicability of numerical modeling to simulate groundwater dynamics in complex hydrogeological situations. This is particularly important for management of the aquifer system investigated in this work, which represents the only available freshwater resource of the study area, and is threatened by overexploitation and saltwater intrusion.

Pressure driven laminar flow of a power-law fluid in a T-channel

NASA Astrophysics Data System (ADS)

Dyakova, O. A.; Frolov, O. Yu

2017-10-01

Planar flow of a non-Newtonian fluid in a T-channel is investigated. The viscosity is determined by the Ostwald-de Waele power law. Motion of the fluid is caused by pressure drop given in boundary sections of the T-channel. On the solid walls, the no slip boundary condition is used. The problem is numerically solved with using a finite difference method based on the SIMPLE procedure. As a result of this study, characteristic flow regimes have been found. Influence of main parameters on the flow pattern has been demonstrated. Criteria dependences describing basic characteristics of the flow under conditions of the present work have been shown.

A direct-inverse method for transonic and separated flows about airfoils

NASA Technical Reports Server (NTRS)

Carlson, Leland A.

1990-01-01

A direct-inverse technique and computer program called TAMSEP that can be used for the analysis of the flow about airfoils at subsonic and low transonic freestream velocities is presented. The method is based upon a direct-inverse nonconservative full potential inviscid method, a Thwaites laminar boundary layer technique, and the Barnwell turbulent momentum integral scheme; and it is formulated using Cartesian coordinates. Since the method utilizes inverse boundary conditions in regions of separated flow, it is suitable for predicting the flow field about airfoils having trailing edge separated flow under high lift conditions. Comparisons with experimental data indicate that the method should be a useful tool for applied aerodynamic analyses.

Effect of wall pattern configurations on Stokes flow through a microchannel with superhydrophobic slip

NASA Astrophysics Data System (ADS)

Mak, H. M.; Ng, C. O.

2010-11-01

The present work aims to study low-Reynolds-number flow through a microchannel with superhydrophobic surfaces, which contain a periodic array of parallel ribs on the upper and lower walls. Mimicking impregnation, the liquid is allowed to penetrate the grooves between the ribs which are filled with an inviscid gas. The array of ribs and grooves gives a heterogeneous wall boundary condition to the channel flow, with partial-slip boundary condition on the solid surface and no-shear boundary condition on the liquid-gas interface. Using the method of eigenfunction expansions and domain decomposition, semi-analytical models are developed for four configurations. Two of them are for longitudinal flow and the others are for transverse flow. For each flow orientation, in-phase and out-phase alignments of ribs between the upper and lower walls are analyzed. The effect of the phase alignments of ribs is appreciable when the channel height is sufficiently small. In-phase alignment gives rise to a larger effective slip length in longitudinal flow. On the contrary, out-phase alignment will yield a larger effective slip length in transverse flow. This work was supported by the Research Grants Council of the Hong Kong Special Administrative Region, China, through Project HKU 7156/09E.

Front propagation in a regular vortex lattice: Dependence on the vortex structure.

PubMed

Beauvier, E; Bodea, S; Pocheau, A

2017-11-01

We investigate the dependence on the vortex structure of the propagation of fronts in stirred flows. For this, we consider a regular set of vortices whose structure is changed by varying both their boundary conditions and their aspect ratios. These configurations are investigated experimentally in autocatalytic solutions stirred by electroconvective flows and numerically from kinematic simulations based on the determination of the dominant Fourier mode of the vortex stream function in each of them. For free lateral boundary conditions, i.e., in an extended vortex lattice, it is found that both the flow structure and the front propagation negligibly depend on vortex aspect ratios. For rigid lateral boundary conditions, i.e., in a vortex chain, vortices involve a slight dependence on their aspect ratios which surprisingly yields a noticeable decrease of the enhancement of front velocity by flow advection. These different behaviors reveal a sensitivity of the mean front velocity on the flow subscales. It emphasizes the intrinsic multiscale nature of front propagation in stirred flows and the need to take into account not only the intensity of vortex flows but also their inner structure to determine front propagation at a large scale. Differences between experiments and simulations suggest the occurrence of secondary flows in vortex chains at large velocity and large aspect ratios.

An Immersed Boundary-Lattice Boltzmann Method for Simulating Particulate Flows

NASA Astrophysics Data System (ADS)

Zhang, Baili; Cheng, Ming; Lou, Jing

2013-11-01

A two-dimensional momentum exchange-based immersed boundary-lattice Boltzmann method developed by X.D. Niu et al. (2006) has been extended in three-dimensions for solving fluid-particles interaction problems. This method combines the most desirable features of the lattice Boltzmann method and the immersed boundary method by using a regular Eulerian mesh for the flow domain and a Lagrangian mesh for the moving particles in the flow field. The non-slip boundary conditions for the fluid and the particles are enforced by adding a force density term into the lattice Boltzmann equation, and the forcing term is simply calculated by the momentum exchange of the boundary particle density distribution functions, which are interpolated by the Lagrangian polynomials from the underlying Eulerian mesh. This method preserves the advantages of lattice Boltzmann method in tracking a group of particles and, at the same time, provides an alternative approach to treat solid-fluid boundary conditions. Numerical validations show that the present method is very accurate and efficient. The present method will be further developed to simulate more complex problems with particle deformation, particle-bubble and particle-droplet interactions.

Numerical Study of Outlet Boundary Conditions for Unsteady Turbulent Internal Flows Using the NCC

NASA Technical Reports Server (NTRS)

Liu, Nan-Suey; Shih, Tsan-Hsing

2009-01-01

This paper presents the results of studies on the outlet boundary conditions for turbulent internal flow simulations. Several outlet boundary conditions have been investigated by applying the National Combustion Code (NCC) to the configuration of a LM6000 single injector flame tube. First of all, very large eddy simulations (VLES) have been performed using the partially resolved numerical simulation (PRNS) approach, in which both the nonlinear and linear dynamic subscale models were employed. Secondly, unsteady Reynolds averaged Navier- Stokes (URANS) simulations have also been performed for the same configuration to investigate the effects of different outlet boundary conditions in the context of URANS. Thirdly, the possible role of the initial condition is inspected by using three different initial flow fields for both the PRNS/VLES simulation and the URANS simulation. The same grid is used for all the simulations and the number of mesh element is about 0.5 million. The main purpose of this study is to examine the long-time behavior of the solution as determined by the imposed outlet boundary conditions. For a particular simulation to be considered as successful under the given initial and boundary conditions, the solution must be sustainable in a physically meaningful manner over a sufficiently long period of time. The commonly used outlet boundary condition for steady Reynolds averaged Navier-Stokes (RANS) simulation is a fixed pressure at the outlet with all the other dependent variables being extrapolated from the interior. The results of the present study suggest that this is also workable for the URANS simulation of the LM6000 injector flame tube. However, it does not work for the PRNS/VLES simulation due to the unphysical reflections of the pressure disturbances at the outlet boundary. This undesirable situation can be practically alleviated by applying a simple unsteady convection equation for the pressure disturbances at the outlet boundary. The numerical results presented in this paper suggest that this unsteady convection of pressure disturbances at the outlet works very well for all the unsteady simulations (both PRNS/VLES and URANS) of the LM6000 single injector flame tube.

Sediment transport under wave groups: Relative importance between nonlinear waveshape and nonlinear boundary layer streaming

USGS Publications Warehouse

Yu, X.; Hsu, T.-J.; Hanes, D.M.

2010-01-01

Sediment transport under nonlinear waves in a predominately sheet flow condition is investigated using a two-phase model. Specifically, we study the relative importance between the nonlinear waveshape and nonlinear boundary layer streaming on cross-shore sand transport. Terms in the governing equations because of the nonlinear boundary layer process are included in this one-dimensional vertical (1DV) model by simplifying the two-dimensional vertical (2DV) ensemble-averaged two-phase equations with the assumption that waves propagate without changing their form. The model is first driven by measured time series of near-bed flow velocity because of a wave group during the SISTEX99 large wave flume experiment and validated with the measured sand concentration in the sheet flow layer. Additional studies are then carried out by including and excluding the nonlinear boundary layer terms. It is found that for the grain diameter (0.24 mm) and high-velocity skewness wave condition considered here, nonlinear waveshape (e.g., skewness) is the dominant mechanism causing net onshore transport and nonlinear boundary layer streaming effect only causes an additional 36% onshore transport. However, for conditions of relatively low-wave skewness and a stronger offshore directed current, nonlinear boundary layer streaming plays a more critical role in determining the net transport. Numerical experiments further suggest that the nonlinear boundary layer streaming effect becomes increasingly important for finer grain. When the numerical model is driven by measured near-bed flow velocity in a more realistic surf zone setting, model results suggest nonlinear boundary layer processes may nearly double the onshore transport purely because of nonlinear waveshape. Copyright 2010 by the American Geophysical Union.

Influence of electrical boundary conditions on molecular dynamics simulations of ionic liquid electrosprays.

PubMed

Borner, Arnaud; Wang, Pengxiang; Levin, Deborah A

2014-12-01

Molecular dynamics (MD) simulations are coupled to solutions of Poisson's equation to study the effects of the electrical boundary conditions on the emission modes of an electrospray thruster fed with an ionic liquid. A comparison of a new tip boundary condition with an analytical model based on a semihyperboloidal shape offers good agreement, although the analytical model overestimates the maximum value of the tangential electric field since it does not take into account the space charge that reduces the field at the liquid surface. It is found that a constant electric field model gives similar agreement to the more rigorous and computationally expensive tip boundary condition at lower flow rates. However, at higher mass flow rates the constant electric field produces extruded particles with higher Coulomb energy per ion, consistent with droplet formation. Furthermore, the MD simulations show that ion emission sites differ based on the boundary condition and snapshots offer an explanation as to why some boundary condition models will predict emission in a purely ionic mode, whereas others suggest a mixed ion-droplet regime. Finally, specific impulses and thrusts are compared for the different models and are found to vary up to 30% due to differences in the average charge to mass ratio.

Influence of electrical boundary conditions on molecular dynamics simulations of ionic liquid electrosprays

NASA Astrophysics Data System (ADS)

Borner, Arnaud; Wang, Pengxiang; Levin, Deborah A.

2014-12-01

Molecular dynamics (MD) simulations are coupled to solutions of Poisson's equation to study the effects of the electrical boundary conditions on the emission modes of an electrospray thruster fed with an ionic liquid. A comparison of a new tip boundary condition with an analytical model based on a semihyperboloidal shape offers good agreement, although the analytical model overestimates the maximum value of the tangential electric field since it does not take into account the space charge that reduces the field at the liquid surface. It is found that a constant electric field model gives similar agreement to the more rigorous and computationally expensive tip boundary condition at lower flow rates. However, at higher mass flow rates the constant electric field produces extruded particles with higher Coulomb energy per ion, consistent with droplet formation. Furthermore, the MD simulations show that ion emission sites differ based on the boundary condition and snapshots offer an explanation as to why some boundary condition models will predict emission in a purely ionic mode, whereas others suggest a mixed ion-droplet regime. Finally, specific impulses and thrusts are compared for the different models and are found to vary up to 30% due to differences in the average charge to mass ratio.

Simulation of Couette flow using conventional Burnett equations with modified slip boundary conditions

NASA Astrophysics Data System (ADS)

Liu, Hualin; Zhao, Wenwen; Chen, Weifang

2016-11-01

Gas or liquid flow through small channels has become more and more popular due to the micro-electro-mechanical systems (MEMS) fabrication technologies such as micro-motors, electrostatic comb-drive, micro-chromatographs, micro-actuators, micro-turbines and micro-pumps, etc. The flow conditions in and around these systems are always recognized as typical transitional regimes. Under these conditions, the mean free path of gas molecules approaches the characteristic scale of the micro-devices itself, and due to the little collisions the heat and momentum cannot equilibrate between the wall and fluids quickly. Couette flow is a simple and critical model in fluid dynamics which focuses on the mechanism of the heat transfer in shear-driven micro-cavities or micro-channels. Despite numerous work on the numerical solutions of the Couette flow, how to propose stable and accurate slip boundary conditions in rarefied flow conditions still remains to be elucidated. In this paper, converged solutions for steady-state micro Couette flows are obtained by using conventional Burnett equations with a set of modified slip boundary conditions. Instead of using the physical variables at the wall, the modified slip conditions use the variables at the edge of the Knudsen layer based on a physically plausible assumption in literature that Knudsen layer has a thickness only in the order of a mean free path and molecules are likely to travel without collision in this layer. Numerical results for non-dimensional wall shear stress and heat flux are compared with those of the DSMC solutions. Although there are not much improvement in the accuracy by using this modified slip conditions, the modified conditions perform much better than the unmodified slip conditions for numerical stabilization. All results show that the set of conventional Burnett equations with second order modified conditions are proved to be an appropriate model for the micro-Couette flows.

Finite-Difference Solutions for Compressible Laminar Boundary-Layer Flows of a Dusty Gas over a Semi-Infinite Flat Plate.

DTIC Science & Technology

1986-08-01

AD-A174 952 FINITE - DIFFERENCE SOLUTIONS FOR CONPRESSIBLE LANINAR 1/2 BOUNDARY-LAYER FLOUS (U) TORONTO UNIV DOWNSVIEW (ONTARIO) INST FOR AEROSPACE...dilute dusty gas over a semi-infinite flat plate. Details are given of the impliit finite , difference schemes as well as the boundary conditions... FINITE - DIFFERENCE SOLUTIONS FOR COMPRESSIBLE LAMINAR BOUNDARY-LAYER FLOWS OF A DUSTY GAS OVER A SEMI-INFINITE FLAT PLATE by B. Y. Wang and I. I

Boundary layer transition observations on a body of revolution with surface heating and cooling in water

NASA Astrophysics Data System (ADS)

Arakeri, V. H.

1980-04-01

Boundary layer flow visualization in water with surface heat transfer was carried out on a body of revolution which had the predicted possibility of laminar separation under isothermal conditions. Flow visualization was by in-line holographic technique. Boundary layer stabilization, including elimination of laminar separation, was observed to take place on surface heating. Conversely, boundary layer destabilization was observed on surface cooling. These findings are consistent with the theoretical predictions of Wazzan et al. (1970).

Scaling Relations and Self-Similarity of 3-Dimensional Reynolds-Averaged Navier-Stokes Equations.

PubMed

Ercan, Ali; Kavvas, M Levent

2017-07-25

Scaling conditions to achieve self-similar solutions of 3-Dimensional (3D) Reynolds-Averaged Navier-Stokes Equations, as an initial and boundary value problem, are obtained by utilizing Lie Group of Point Scaling Transformations. By means of an open-source Navier-Stokes solver and the derived self-similarity conditions, we demonstrated self-similarity within the time variation of flow dynamics for a rigid-lid cavity problem under both up-scaled and down-scaled domains. The strength of the proposed approach lies in its ability to consider the underlying flow dynamics through not only from the governing equations under consideration but also from the initial and boundary conditions, hence allowing to obtain perfect self-similarity in different time and space scales. The proposed methodology can be a valuable tool in obtaining self-similar flow dynamics under preferred level of detail, which can be represented by initial and boundary value problems under specific assumptions.

Boundary condition computational procedures for inviscid, supersonic steady flow field calculations

NASA Technical Reports Server (NTRS)

Abbett, M. J.

1971-01-01

Results are given of a comparative study of numerical procedures for computing solid wall boundary points in supersonic inviscid flow calculatons. Twenty five different calculation procedures were tested on two sample problems: a simple expansion wave and a simple compression (two-dimensional steady flow). A simple calculation procedure was developed. The merits and shortcomings of the various procedures are discussed, along with complications for three-dimensional and time-dependent flows.

Knowledge Sharing and Organizational Change in Higher Education

ERIC Educational Resources Information Center

Dee, Jay; Leisyte, Liudvika

2017-01-01

Purpose: Organizational learning in higher education institutions depends upon the ability of managers and academics to maintain a flow of knowledge across the structural boundaries of the university. This paper aims to understand the boundary conditions that foster or impede the flow of knowledge during organizational change at a large public…

The boundary layer as a means of controlling the flow of liquids and gases

NASA Technical Reports Server (NTRS)

Schrenk, Oskar

1930-01-01

According to one of the main propositions of the boundary layer theory the scarcely noticeable boundary layer may, under certain conditions, have a decisive influence on the form of the external flow by causing it to separate from the wing surface. These phenomena are known to be caused by a kind of stagnation of the boundary layer at the point of separation. The present report deals with similar phenomena. It is important to note that usually the cause (external interference) directly affects only the layer close to the wall, while its indirect effect extends to a large portion of the external flow.

Zero-G experiments in two-phase fluids flow regimes

NASA Technical Reports Server (NTRS)

Heppner, D. B.; King, C. D.; Littles, J. W.

1975-01-01

The two-phase flows studied were liquid and gas mixtures in a straight flow channel of circular cross-section. Boundaries between flow regimes have been defined for normogravity on coordinates of gas quality and total mass velocity; and, when combined with boundary expressions having a Froude number term, an analytical model was derived predicting boundary shifts with changes in gravity level. Experiments with air and water were performed, first in the normogravity environment of a ground laboratory and then in 'zero gravity' aboard a KC-135 aircraft flying parabolic trajectories. Data reduction confirmed regime boundary shifts in the direction predicted, although the magnitude was a little less than predicted. Pressure drop measurements showed significant increases for the low gravity condition.

Hypersonic Boundary Layer Instability Over a Corner

NASA Technical Reports Server (NTRS)

Balakumar, Ponnampalam; Zhao, Hong-Wu; McClinton, Charles (Technical Monitor)

2001-01-01

A boundary-layer transition study over a compression corner was conducted under a hypersonic flow condition. Due to the discontinuities in boundary layer flow, the full Navier-Stokes equations were solved to simulate the development of disturbance in the boundary layer. A linear stability analysis and PSE method were used to get the initial disturbance for parallel and non-parallel flow respectively. A 2-D code was developed to solve the full Navier-stokes by using WENO(weighted essentially non-oscillating) scheme. The given numerical results show the evolution of the linear disturbance for the most amplified disturbance in supersonic and hypersonic flow over a compression ramp. The nonlinear computations also determined the minimal amplitudes necessary to cause transition at a designed location.

Effects of the bottom boundary condition in numerical investigations of dense water cascading on a slope

NASA Astrophysics Data System (ADS)

Berntsen, Jarle; Alendal, Guttorm; Avlesen, Helge; Thiem, Øyvind

2018-05-01

The flow of dense water along continental slopes is considered. There is a large literature on the topic based on observations and laboratory experiments. In addition, there are many analytical and numerical studies of dense water flows. In particular, there is a sequence of numerical investigations using the dynamics of overflow mixing and entrainment (DOME) setup. In these papers, the sensitivity of the solutions to numerical parameters such as grid size and numerical viscosity coefficients and to the choices of methods and models is investigated. In earlier DOME studies, three different bottom boundary conditions and a range of vertical grid sizes are applied. In other parts of the literature on numerical studies of oceanic gravity currents, there are statements that appear to contradict choices made on bottom boundary conditions in some of the DOME papers. In the present study, we therefore address the effects of the bottom boundary condition and vertical resolution in numerical investigations of dense water cascading on a slope. The main finding of the present paper is that it is feasible to capture the bottom Ekman layer dynamics adequately and cost efficiently by using a terrain-following model system using a quadratic drag law with a drag coefficient computed to give near-bottom velocity profiles in agreement with the logarithmic law of the wall. Many studies of dense water flows are performed with a quadratic bottom drag law and a constant drag coefficient. It is shown that when using this bottom boundary condition, Ekman drainage will not be adequately represented. In other studies of gravity flow, a no-slip bottom boundary condition is applied. With no-slip and a very fine resolution near the seabed, the solutions are essentially equal to the solutions obtained with a quadratic drag law and a drag coefficient computed to produce velocity profiles matching the logarithmic law of the wall. However, with coarser resolution near the seabed, there may be a substantial artificial blocking effect when using no-slip.

Off-Body Boundary-Layer Measurement Techniques Development for Supersonic Low-Disturbance Flows

NASA Technical Reports Server (NTRS)

Owens, Lewis R.; Kegerise, Michael A.; Wilkinson, Stephen P.

2011-01-01

Investigations were performed to develop accurate boundary-layer measurement techniques in a Mach 3.5 laminar boundary layer on a 7 half-angle cone at 0 angle of attack. A discussion of the measurement challenges is presented as well as how each was addressed. A computational study was performed to minimize the probe aerodynamic interference effects resulting in improved pitot and hot-wire probe designs. Probe calibration and positioning processes were also developed with the goal of reducing the measurement uncertainties from 10% levels to less than 5% levels. Efforts were made to define the experimental boundary conditions for the cone flow so comparisons could be made with a set of companion computational simulations. The development status of the mean and dynamic boundary-layer flow measurements for a nominally sharp cone in a low-disturbance supersonic flow is presented.

Numerical analysis of hypersonic turbulent film cooling flows

NASA Technical Reports Server (NTRS)

Chen, Y. S.; Chen, C. P.; Wei, H.

1992-01-01

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

A closed-form analytical model for predicting 3D boundary layer displacement thickness for the validation of viscous flow solvers

NASA Astrophysics Data System (ADS)

Kumar, V. R. Sanal; Sankar, Vigneshwaran; Chandrasekaran, Nichith; Saravanan, Vignesh; Natarajan, Vishnu; Padmanabhan, Sathyan; Sukumaran, Ajith; Mani, Sivabalan; Rameshkumar, Tharikaa; Nagaraju Doddi, Hema Sai; Vysaprasad, Krithika; Sharan, Sharad; Murugesh, Pavithra; Shankar, S. Ganesh; Nejaamtheen, Mohammed Niyasdeen; Baskaran, Roshan Vignesh; Rahman Mohamed Rafic, Sulthan Ariff; Harisrinivasan, Ukeshkumar; Srinivasan, Vivek

2018-02-01

A closed-form analytical model is developed for estimating the 3D boundary-layer-displacement thickness of an internal flow system at the Sanal flow choking condition for adiabatic flows obeying the physics of compressible viscous fluids. At this unique condition the boundary-layer blockage induced fluid-throat choking and the adiabatic wall-friction persuaded flow choking occur at a single sonic-fluid-throat location. The beauty and novelty of this model is that without missing the flow physics we could predict the exact boundary-layer blockage of both 2D and 3D cases at the sonic-fluid-throat from the known values of the inlet Mach number, the adiabatic index of the gas and the inlet port diameter of the internal flow system. We found that the 3D blockage factor is 47.33 % lower than the 2D blockage factor with air as the working fluid. We concluded that the exact prediction of the boundary-layer-displacement thickness at the sonic-fluid-throat provides a means to correctly pinpoint the causes of errors of the viscous flow solvers. The methodology presented herein with state-of-the-art will play pivotal roles in future physical and biological sciences for a credible verification, calibration and validation of various viscous flow solvers for high-fidelity 2D/3D numerical simulations of real-world flows. Furthermore, our closed-form analytical model will be useful for the solid and hybrid rocket designers for the grain-port-geometry optimization of new generation single-stage-to-orbit dual-thrust-motors with the highest promising propellant loading density within the given envelope without manifestation of the Sanal flow choking leading to possible shock waves causing catastrophic failures.

Real Gas Scale Effects on Hypersonic Laminar Boundary-Layer Parameters Including Effects of Entropy-Layer Swallowing

DTIC Science & Technology

1975-12-01

crossed the essentially normal portion of the bow shock is swallowed by the boundary layer. The flow along the edge of the boundary layer on the aft...portions hf the body will then have passed through an oblique part of the bow shock and will be in a different state than had it passed through a normal...determination of the local edge flow conditions may be improvedby taking into con- sideration the inclination of the bow shock where the local flow stream- line

Reynolds number scaling of the influence of boundary layers on the global behavior of laboratory quasi-Keplerian flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Edlund, E. M.; Ji, H.

2015-10-06

Here, we present fluid velocity measurements in a modified Taylor-Couette device operated in the quasi-Keplerian regime, where it is observed that nearly ideal flows exhibit self-similarity under scaling of the Reynolds number. In contrast, nonideal flows show progressive departure from ideal Couette as the Reynolds number is increased. We present a model that describes the observed departures from ideal Couette rotation as a function of the fluxes of angular momentum across the boundaries, capturing the dependence on Reynolds number and boundary conditions.

Reynolds number scaling of the influence of boundary layers on the global behavior of laboratory quasi-Keplerian flows.

PubMed

Edlund, E M; Ji, H

2015-10-01

We present fluid velocity measurements in a modified Taylor-Couette device operated in the quasi-Keplerian regime, where it is observed that nearly ideal flows exhibit self-similarity under scaling of the Reynolds number. In contrast, nonideal flows show progressive departure from ideal Couette as the Reynolds number is increased. We present a model that describes the observed departures from ideal Couette rotation as a function of the fluxes of angular momentum across the boundaries, capturing the dependence on Reynolds number and boundary conditions.

Boundary-Layer-Ingesting Inlet Flow Control

NASA Technical Reports Server (NTRS)

Owens, Lewis R.; Allan, Brian G.; Gorton, Susan A.

2008-01-01

An experimental study was conducted to provide the first demonstration of an active flow control system for a flush-mounted inlet with significant boundary-layer-ingestion in transonic flow conditions. The effectiveness of the flow control in reducing the circumferential distortion at the engine fan-face location was assessed using a 2.5%-scale model of a boundary-layer-ingesting offset diffusing inlet. The inlet was flush mounted to the tunnel wall and ingested a large boundary layer with a boundary-layer-to-inlet height ratio of 35%. Different jet distribution patterns and jet mass flow rates were used in the inlet to control distortion. A vane configuration was also tested. Finally a hybrid vane/jet configuration was tested leveraging strengths of both types of devices. Measurements were made of the onset boundary layer, the duct surface static pressures, and the mass flow rates through the duct and the flow control actuators. The distortion and pressure recovery were measured at the aerodynamic interface plane. The data show that control jets and vanes reduce circumferential distortion to acceptable levels. The point-design vane configuration produced higher distortion levels at off-design settings. The hybrid vane/jet flow control configuration reduced the off-design distortion levels to acceptable ones and used less than 0.5% of the inlet mass flow to supply the jets.

A new boundary scheme for simulation of gas flow in kerogen pores with considering surface diffusion effect

NASA Astrophysics Data System (ADS)

Wang, Lingquan; Zeng, Zhong; Zhang, Liangqi; Qiao, Long; Zhang, Yi; Lu, Yiyu

2018-04-01

Navier-Stokes (NS) equations with no-slip boundary conditions fail to realistically describe micro-flows with considering nanoscale phenomena. Particularly, in kerogen pores, slip-flow and surface diffusion are important. In this study, we propose a new slip boundary scheme for the lattice Boltzmann (LB) method through the non-equilibrium extrapolation scheme to simulate the slip-flow considering surface diffusion effect. Meanwhile, the second-order slip velocity can be taken into account. The predicted characteristics in a two-dimensional micro-flow, including slip-velocity, velocity distribution along the flow direction with/without surface diffusion are present. The results in this study are compared with available analytical and reference results, and good agreements are achieved.

Modeling Plasma Turbulence and Flows in LAPD using BOUT++

NASA Astrophysics Data System (ADS)

Friedman, B.; Carter, T. A.; Schaffner, D.; Popovich, P.; Umansky, M. V.; Dudson, B.

2010-11-01

A Braginskii fluid model of plasma turbulence in the BOUT code has recently been applied to LAPD at UCLA [1]. While these initial simulations with a reduced model and periodic axial boundary conditions have shown good agreement with measurements (e.g. power spectrum, correlation lengths), these simulations have lacked physics essential for modeling self-consistent, quantitatively correct flows. In particular, the model did not contain parallel plasma flow induced by sheath boundary conditions, and the axisymmetric radial electric field was not consistent with experiment. This work addresses these issues by extending the simulation model in the BOUT++ code [2], a more advanced version of BOUT. Specifically, end-plate sheath boundary conditions are added, as well as equations to evolve electron temperature and parallel ion velocity. Finally, various techniques are used to attempt to match the experimental electric potential profile, including fixing an equilibrium profile, fixing the radial boundaries, and adding an angular momentum source. [4pt] [1] Popovich et al., http://arxiv.org/abs/1005.2418 (2010).[0pt] [2] Dudson et al., Computer Physics Communications 180 (2009).

The Perfectly Matched Layer absorbing boundary for fluid-structure interactions using the Immersed Finite Element Method.

PubMed

Yang, Jubiao; Yu, Feimi; Krane, Michael; Zhang, Lucy T

2018-01-01

In this work, a non-reflective boundary condition, the Perfectly Matched Layer (PML) technique, is adapted and implemented in a fluid-structure interaction numerical framework to demonstrate that proper boundary conditions are not only necessary to capture correct wave propagations in a flow field, but also its interacted solid behavior and responses. While most research on the topics of the non-reflective boundary conditions are focused on fluids, little effort has been done in a fluid-structure interaction setting. In this study, the effectiveness of the PML is closely examined in both pure fluid and fluid-structure interaction settings upon incorporating the PML algorithm in a fully-coupled fluid-structure interaction framework, the Immersed Finite Element Method. The performance of the PML boundary condition is evaluated and compared to reference solutions with a variety of benchmark test cases including known and expected solutions of aeroacoustic wave propagation as well as vortex shedding and advection. The application of the PML in numerical simulations of fluid-structure interaction is then investigated to demonstrate the efficacy and necessity of such boundary treatment in order to capture the correct solid deformation and flow field without the requirement of a significantly large computational domain.

Measurements in a Transitioning Cone Boundary Layer at Freestream Mach 3.5

NASA Technical Reports Server (NTRS)

King, Rudolph A.; Chou, Amanda; Balakumar, Ponnampalam; Owens, Lewis R.; Kegerise, Michael A.

2016-01-01

An experimental study was conducted in the Supersonic Low-Disturbance Tunnel to investigate naturally-occurring instabilities in a supersonic boundary layer on a 7 deg half- angle cone. All tests were conducted with a nominal freestream Mach number of M(sub infinity) = 3:5, total temperature of T(sub 0) = 299:8 K, and unit Reynolds numbers of Re(sub infinity) x 10(exp -6) = 9:89, 13.85, 21.77, and 25.73 m(exp -1). Instability measurements were acquired under noisy- ow and quiet- ow conditions. Measurements were made to document the freestream and the boundary-layer edge environment, to document the cone baseline flow, and to establish the stability characteristics of the transitioning flow. Pitot pressure and hot-wire boundary- layer measurements were obtained using a model-integrated traverse system. All hot- wire results were single-point measurements and were acquired with a sensor calibrated to mass ux. For the noisy-flow conditions, excellent agreement for the growth rates and mode shapes was achieved between the measured results and linear stability theory (LST). The corresponding N factor at transition from LST is N 3:9. The stability measurements for the quiet-flow conditions were limited to the aft end of the cone. The most unstable first-mode instabilities as predicted by LST were successfully measured, but this unstable first mode was not the dominant instability measured in the boundary layer. Instead, the dominant instabilities were found to be the less-amplified, low-frequency disturbances predicted by linear stability theory, and these instabilities grew according to linear theory. These low-frequency unstable disturbances were initiated by freestream acoustic disturbances through a receptivity process that is believed to occur near the branch I locations of the cone. Under quiet-flow conditions, the boundary layer remained laminar up to the last measurement station for the largest Re1, implying a transition N factor of N greater than 8:5.

Hydrodynamic boundary condition of water on hydrophobic surfaces.

PubMed

Schaeffel, David; Yordanov, Stoyan; Schmelzeisen, Marcus; Yamamoto, Tetsuya; Kappl, Michael; Schmitz, Roman; Dünweg, Burkhard; Butt, Hans-Jürgen; Koynov, Kaloian

2013-05-01

By combining total internal reflection fluorescence cross-correlation spectroscopy with Brownian dynamics simulations, we were able to measure the hydrodynamic boundary condition of water flowing over a smooth solid surface with exceptional accuracy. We analyzed the flow of aqueous electrolytes over glass coated with a layer of poly(dimethylsiloxane) (advancing contact angle Θ = 108°) or perfluorosilane (Θ = 113°). Within an error of better than 10 nm the slip length was indistinguishable from zero on all surfaces.

Shock Tunnel Operation and Correlation of Boundary Layer Transition on a Cone in Hypervelocity Flow

DTIC Science & Technology

2013-07-01

conditions from the ideal reflected-shock pressure to measured reservoir pressure using an isentropic expan- sion. Furthermore, the 1-D nozzle computation...does not account for boundary layer growth within the nozzle , off-design operation conditions that lead to flow nonuni- formity, or vibration...translation nonequilibrium and freezing within the nozzle , which is significant for the N2 cases. For the uncertainties that can be quantified, we have combined

Discontinuous Galerkin Methods for Turbulence Simulation

NASA Technical Reports Server (NTRS)

Collis, S. Scott

2002-01-01

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Improving a complex finite-difference ground water flow model through the use of an analytic element screening model

USGS Publications Warehouse

Hunt, R.J.; Anderson, M.P.; Kelson, V.A.

1998-01-01

This paper demonstrates that analytic element models have potential as powerful screening tools that can facilitate or improve calibration of more complicated finite-difference and finite-element models. We demonstrate how a two-dimensional analytic element model was used to identify errors in a complex three-dimensional finite-difference model caused by incorrect specification of boundary conditions. An improved finite-difference model was developed using boundary conditions developed from a far-field analytic element model. Calibration of a revised finite-difference model was achieved using fewer zones of hydraulic conductivity and lake bed conductance than the original finite-difference model. Calibration statistics were also improved in that simulated base-flows were much closer to measured values. The improved calibration is due mainly to improved specification of the boundary conditions made possible by first solving the far-field problem with an analytic element model.This paper demonstrates that analytic element models have potential as powerful screening tools that can facilitate or improve calibration of more complicated finite-difference and finite-element models. We demonstrate how a two-dimensional analytic element model was used to identify errors in a complex three-dimensional finite-difference model caused by incorrect specification of boundary conditions. An improved finite-difference model was developed using boundary conditions developed from a far-field analytic element model. Calibration of a revised finite-difference model was achieved using fewer zones of hydraulic conductivity and lake bed conductance than the original finite-difference model. Calibration statistics were also improved in that simulated base-flows were much closer to measured values. The improved calibration is due mainly to improved specification of the boundary conditions made possible by first solving the far-field problem with an analytic element model.

Evaluation of Wall Boundary Conditions for Impedance Eduction Using a Dual-Source Method

NASA Technical Reports Server (NTRS)

Watson, W. R.; Jones, M. G.

2012-01-01

The accuracy of the Ingard-Myers boundary condition and a recently proposed modified Ingard-Myers boundary condition is evaluated for use in impedance eduction under the assumption of uniform mean flow. The evaluation is performed at three centerline Mach numbers, using data acquired in a grazing flow impedance tube, using both upstream and downstream propagating sound sources, and on a database of test liners for which the expected behavior of the impedance spectra is known. The test liners are a hard-wall insert consisting of 12.6 mm thick aluminum, a linear liner without a facesheet consisting of a number of small diameter but long cylindrical channels embedded in a ceramic material, and two conventional nonlinear liners consisting of a perforated facesheet bonded to a honeycomb core. The study is restricted to a frequency range for which only plane waves are cut on in the hard-wall sections of the flow impedance tube. The metrics used to evaluate each boundary condition are 1) how well it educes the same impedance for upstream and downstream propagating sources, and 2) how well it predicts the expected behavior of the impedance spectra over the Mach number range. The primary conclusions of the study are that the same impedance is educed for upstream and downstream propagating sources except at the highest Mach number, that an effective impedance based on both the upstream and downstream measurements is more accurate than an impedance based on the upstream or downstream data alone, and that the Ingard-Myers boundary condition with an effective impedance produces results similar to that achieved with the modified Ingard-Myers boundary condition.

Numerical modeling of Stokes flows over a superhydrophobic surface containing gas bubbles

NASA Astrophysics Data System (ADS)

Ageev, A. I.; Golubkina, I. V.; Osiptsov, A. N.

2017-10-01

This paper continues the numerical modeling of Stokes flows near cavities of a superhydrophobic surface, occupied by gas bubbles, based on the Boundary Element Method (BEM). The aim of the present study is to estimate the friction reduction (pressure drop) in a microchannel with a bottom superhydrophobic surface, the texture of which is formed by a periodic system of striped rectangular microcavities containing compressible gas bubbles. The model proposed takes into account the streamwise variation of the bubble shift into the cavities, caused by the longitudinal pressure gradient in the channel flow. The solution for the macroscopic (averaged) flow in the microchannel, constructed using an effective slip boundary condition on the superhydrophobic bottom wall, is matched with the solution of the Stokes problem at the microscale of a single cavity containing a gas bubble. The 2D Stokes problems of fluid flow over single cavities containing curved phase interfaces with the condition of zero shear stress are reduced to the boundary integral equations which are solved using the BEM method.

Which key properties controls the preferential transport in the vadose zone under transient hydrological conditions

NASA Astrophysics Data System (ADS)

Groh, J.; Vanderborght, J.; Puetz, T.; Gerke, H. H.; Rupp, H.; Wollschlaeger, U.; Stumpp, C.; Priesack, E.; Vereecken, H.

2015-12-01

Understanding water flow and solute transport in the unsaturated zone is of great importance for an appropriate land use management strategy. The quantification and prediction of water and solute fluxes through the vadose zone can help to improve management practices in order to limit potential risk on our fresh water resources. Water related solute transport and residence time is strongly affected by preferential flow paths in the soil. Water flow in soils depends on soil properties and site factors (climate or experiment conditions, land use) and are therefore important factors to understand preferential solute transport in the unsaturated zone. However our understanding and knowledge of which on-site properties or conditions define and enhance preferential flow and transport is still poor and mostly limited onto laboratory experimental conditions (small column length and steady state boundary conditions). Within the TERENO SOILCan lysimeter network, which was designed to study the effects of climate change on soil functions, a bromide tracer was applied on 62 lysimeter at eight different test sites between Dec. 2013 and Jan. 2014. The TERENO SOILCan infrastructure offers the unique possibility to study the occurrence of preferential flow and transport of various soil types under different natural transient hydrological conditions and land use (crop, bare and grassland) at eight TERENO SOILCan observatories. Working with lysimeter replicates at each observatory allows defining the spatial variability of preferential transport and flow. Additionally lysimeters in the network were transferred within and between observatories in order to subject them to different rainfall and temperature regimes and enable us to relate the soil type susceptibility of preferential flow and transport not only to site specific physical and land use properties, but also to different transient boundary conditions. Comparison and statistical analysis between preferential flow indicators 5% arrival time and potential key soil properties, site factors and boundary conditions will be presented in order to identify key properties which control the preferential transport in the vadose zone under transient hydrological conditions.

Velocity boundary conditions for vorticity formulations of the incompressible Navier-Stokes equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kempka, S.N.; Strickland, J.H.; Glass, M.W.

1995-04-01

formulation to satisfy velocity boundary conditions for the vorticity form of the incompressible, viscous fluid momentum equations is presented. The tangential and normal components of the velocity boundary condition are satisfied simultaneously by creating vorticity adjacent to boundaries. The newly created vorticity is determined using a kinematical formulation which is a generalization of Helmholtz` decomposition of a vector field. Though it has not been generally recognized, these formulations resolve the over-specification issue associated with creating voracity to satisfy velocity boundary conditions. The generalized decomposition has not been widely used, apparently due to a lack of a useful physical interpretation. Anmore » analysis is presented which shows that the generalized decomposition has a relatively simple physical interpretation which facilitates its numerical implementation. The implementation of the generalized decomposition is discussed in detail. As an example the flow in a two-dimensional lid-driven cavity is simulated. The solution technique is based on a Lagrangian transport algorithm in the hydrocode ALEGRA. ALEGRA`s Lagrangian transport algorithm has been modified to solve the vorticity transport equation and the generalized decomposition, thus providing a new, accurate method to simulate incompressible flows. This numerical implementation and the new boundary condition formulation allow vorticity-based formulations to be used in a wider range of engineering problems.« less

Further investigation of a finite difference procedure for analyzing the transonic flow about harmonically oscillating airfoils and wings

NASA Technical Reports Server (NTRS)

Weatherill, W. H.; Ehlers, F. E.; Yip, E.; Sebastian, J. D.

1980-01-01

Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. The steady velocity potential is obtained first from the well-known nonlinear equation for steady transonic flow. The unsteady velocity potential is then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. An out-of-core direct solution procedure was developed and applied to two-dimensional sections. Results are presented for a section of vanishing thickness in subsonic flow and an NACA 64A006 airfoil in supersonic flow. Good correlation is obtained in the first case at values of Mach number and reduced frequency of direct interest in flutter analyses. Reasonable results are obtained in the second case. Comparisons of two-dimensional finite difference solutions with exact analytic solutions indicate that the accuracy of the difference solution is dependent on the boundary conditions used on the outer boundaries. Homogeneous boundary conditions on the mesh edges that yield complex eigenvalues give the most accurate finite difference solutions. The plane outgoing wave boundary conditions meet these requirements.

An integral transform approach for a mixed boundary problem involving a flowing partially penetrating well with infinitesimal well skin

NASA Astrophysics Data System (ADS)

Chang, Chien-Chieh; Chen, Chia-Shyun

2002-06-01

A flowing partially penetrating well with infinitesimal well skin is a mixed boundary because a Cauchy condition is prescribed along the screen length and a Neumann condition of no flux is stipulated over the remaining unscreened part. An analytical approach based on the integral transform technique is developed to determine the Laplace domain solution for such a mixed boundary problem in a confined aquifer of finite thickness. First, the mixed boundary is changed into a homogeneous Neumann boundary by substituting the Cauchy condition with a Neumann condition in terms of well bore flux that varies along the screen length and is time dependent. Despite the well bore flux being unknown a priori, the modified model containing this homogeneous Neumann boundary can be solved with the Laplace and the finite Fourier cosine transforms. To determine well bore flux, screen length is discretized into a finite number of segments, to which the Cauchy condition is reinstated. This reinstatement also restores the relation between the original model and the solutions obtained. For a given time, the numerical inversion of the Laplace domain solution yields the drawdown distributions, well bore flux, and the well discharge. This analytical approach provides an alternative for dealing with the mixed boundary problems, especially when aquifer thickness is assumed to be finite.

Considerations on the moving contact-line singularity, with application to frictional drag on a slender drop

NASA Technical Reports Server (NTRS)

Durbin, P. A.

1988-01-01

It has previously been shown that the no-slip boundary conditions leads to a singularity at a moving contact line and that this presumes some form of slip. Present considerations on the energetics of slip due to shear stress lead to a yield stress boundary condition. A model for the distortion of the liquid state near solid boundaries gives a physical basis for this boundary condition. The yield stress condition is illustrated by an analysis of a slender drop rolling down an incline. That analysis provides a formula for the frictional drag resisting the drop movement. With the present boundary condition, the length of the slip region becomes a property of the fluid flow.

Computational Analysis of Human Blood Flow

NASA Astrophysics Data System (ADS)

Panta, Yogendra; Marie, Hazel; Harvey, Mark

2009-11-01

Fluid flow modeling with commercially available computational fluid dynamics (CFD) software is widely used to visualize and predict physical phenomena related to various biological systems. In this presentation, a typical human aorta model was analyzed assuming the blood flow as laminar with complaint cardiac muscle wall boundaries. FLUENT, a commercially available finite volume software, coupled with Solidworks, a modeling software, was employed for the preprocessing, simulation and postprocessing of all the models.The analysis mainly consists of a fluid-dynamics analysis including a calculation of the velocity field and pressure distribution in the blood and a mechanical analysis of the deformation of the tissue and artery in terms of wall shear stress. A number of other models e.g. T branches, angle shaped were previously analyzed and compared their results for consistency for similar boundary conditions. The velocities, pressures and wall shear stress distributions achieved in all models were as expected given the similar boundary conditions. The three dimensional time dependent analysis of blood flow accounting the effect of body forces with a complaint boundary was also performed.

Dependence of energy characteristics of ascending swirling air flow on velocity of vertical blowing

NASA Astrophysics Data System (ADS)

Volkov, R. E.; Obukhov, A. G.; Kutrunov, V. N.

2018-05-01

In the model of a compressible continuous medium, for the complete Navier-Stokes system of equations, an initial boundary problem is proposed that corresponds to the conducted and planned experiments and describes complex three-dimensional flows of a viscous compressible heat-conducting gas in ascending swirling flows that are initiated by a vertical cold blowing. Using parallelization methods, three-dimensional nonstationary flows of a polytropic viscous compressible heat-conducting gas are constructed numerically in different scaled ascending swirling flows under the condition when gravity and Coriolis forces act. With the help of explicit difference schemes and the proposed initial boundary conditions, approximate solutions of the complete system of Navier-Stokes equations are constructed as well as the velocity and energy characteristics of three-dimensional nonstationary gas flows in ascending swirling flows are determined.

Numerical study of the effects of Planetary Boundary Layer structure on the pollutant dispersion within built-up areas.

PubMed

Miao, Yucong; Liu, Shuhua; Zheng, Yijia; Wang, Shu; Liu, Zhenxin; Zhang, Bihui

2015-06-01

The effects of different Planetary Boundary Layer (PBL) structures on pollutant dispersion processes within two idealized street canyon configurations and a realistic urban area were numerically examined by a Computational Fluid Dynamics (CFD) model. The boundary conditions of different PBL structures/conditions were provided by simulations of the Weather Researching and Forecasting model. The simulated results of the idealized 2D and 3D street canyon experiments showed that the increment of PBL instability favored the downward transport of momentum from the upper flow above the roof to the pedestrian level within the street canyon. As a result, the flow and turbulent fields within the street canyon under the more unstable PBL condition are stronger. Therefore, more pollutants within the street canyon would be removed by the stronger advection and turbulent diffusion processes under the unstable PBL condition. On the contrary, more pollutants would be concentrated in the street canyon under the stable PBL condition. In addition, the simulations of the realistic building cluster experiments showed that the density of buildings was a crucial factor determining the dynamic effects of the PBL structure on the flow patterns. The momentum field within a denser building configuration was mostly transported from the upper flow, and was more sensitive to the PBL structures than that of the sparser building configuration. Finally, it was recommended to use the Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme, which can explicitly output the needed turbulent variables, to provide the boundary conditions to the CFD simulation. Copyright © 2015. Published by Elsevier B.V.

Two-phase simulations of the full load surge in Francis turbines

NASA Astrophysics Data System (ADS)

Wack, J.; Riedelbauch, S.

2016-11-01

At off-design conditions, Francis turbines experience cavitation which may reduce the power output and can cause severe damage in the machine. Certain conditions can cause self-excited oscillations of the vortex rope in the draft tube at full load operating point. For the presented work, two-phase simulations are carried out at model scale on a domain ranging from the inlet of the spiral case to the outlet of the draft tube. At different locations, wall pressure measurements are available and compared to the simulation results. Furthermore, the dynamics of the cavity volume in the draft tube cone and at the trailing edge of the runner blades are investigated by comparing with high speed visualization. To account for the selfexcited behaviour, proper boundary conditions need to be set. In this work, the focus lies on the treatment of the boundary condition at the inlet. In the first step, the dynamic behaviour of the cavity regions is investigated using a constant mass flow. Thereafter, oscillations of the total pressure and mass flow rate are prescribed using various frequencies and amplitudes. This methodology enables to examine the response of the cavity dynamics due to different excitations. It can be observed that setting a constant mass flow boundary condition is not suitable to account for the self-excited behaviour. Prescribing the total pressure has the result that the frequency of the vapour volume oscillation is the same as the frequency of the excitation signal. Contrary to that, for an excitation with a mass flow boundary condition, the response of the system is not equal to the excitation.

A direct-inverse method for transonic and separated flows about airfoils

NASA Technical Reports Server (NTRS)

Carlson, K. D.

1985-01-01

A direct-inverse technique and computer program called TAMSEP that can be sued for the analysis of the flow about airfoils at subsonic and low transonic freestream velocities is presented. The method is based upon a direct-inverse nonconservative full potential inviscid method, a Thwaites laminar boundary layer technique, and the Barnwell turbulent momentum integral scheme; and it is formulated using Cartesian coordinates. Since the method utilizes inverse boundary conditions in regions of separated flow, it is suitable for predicing the flowfield about airfoils having trailing edge separated flow under high lift conditions. Comparisons with experimental data indicate that the method should be a useful tool for applied aerodynamic analyses.

A numerical method for the solution of internal pipe/channel flows in laminar or turbulent motion

NASA Astrophysics Data System (ADS)

Lourenco, L.; Essers, J. A.

1981-11-01

A computer program which is useful in the solution of problems of internal turbulent or laminar flow without recirculation is described. The flow is treated in terms of parabolic boundary layer differential equations. The eddy diffusivity concept is used to model turbulent stresses. Two turbulent models are available: the Prandtl mixing length model and the Nee-Kovasznay model for the effective viscosity. Fluid is considered incompressible, but little program modification is needed to treat compressible flows. Initial conditions are prescribed as well as the boundary conditions. The differencing scheme employed is fully implicit for the dependent variables. This allows the use of relatively large forward steps without stability problems.

Low-to-moderate Reynolds number swirling flow in an annular channel with a rotating end wall.

PubMed

Davoust, Laurent; Achard, Jean-Luc; Drazek, Laurent

2015-02-01

This paper presents a new method for solving analytically the axisymmetric swirling flow generated in a finite annular channel from a rotating end wall, with no-slip boundary conditions along stationary side walls and a slip condition along the free surface opposite the rotating floor. In this case, the end-driven swirling flow can be described from the coupling between an azimuthal shear flow and a two-dimensional meridional flow driven by the centrifugal force along the rotating floor. A regular asymptotic expansion based on a small but finite Reynolds number is used to calculate centrifugation-induced first-order correction to the azimuthal Stokes flow obtained as the solution at leading order. For solving the first-order problem, the use of an integral boundary condition for the vorticity is found to be a convenient way to attribute boundary conditions in excess for the stream function to the vorticity. The annular geometry is characterized by both vertical and horizontal aspect ratios, whose respective influences on flow patterns are investigated. The vertical aspect ratio is found to involve nontrivial changes in flow patterns essentially due to the role of corner eddies located on the left and right sides of the rotating floor. The present analytical method can be ultimately extended to cylindrical geometries, irrespective of the surface opposite the rotating floor: a wall or a free surface. It can also serve as an analytical tool for monitoring confined rotating flows in applications related to surface viscosimetry or crystal growth from the melt.

Locomotion of bacteria in liquid flow and the boundary layer effect on bacterial attachment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Chao, E-mail: zhangchao@cqu.edu.cn; Institute of Engineering Thermophysics, Chongqing University, Chongqing 400030; Liao, Qiang, E-mail: lqzx@cqu.edu.cn

The formation of biofilm greatly affects the performance of biological reactors, which highly depends on bacterial swimming and attachment that usually takes place in liquid flow. Therefore, bacterial swimming and attachment on flat and circular surfaces with the consideration of flow was studied experimentally. Besides, a mathematical model comprehensively combining bacterial swimming and motion with flow is proposed for the simulation of bacterial locomotion and attachment in flow. Both experimental and theoretical results revealed that attached bacteria density increases with decreasing boundary layer thickness on both flat and circular surfaces, the consequence of which is inherently related to the competitionmore » between bacterial swimming and the non-slip motion with flow evaluated by the Péclet number. In the boundary layer, where the Péclet number is relatively higher, bacterial locomotion mainly depends on bacterial swimming. Thinner boundary layer promotes bacterial swimming towards the surface, leading to higher attachment density. To enhance the performance of biofilm reactors, it is effective to reduce the boundary layer thickness on desired surfaces. - Highlights: • Study of bacterial locomotion in flow as an early stage in biofilm formation. • Mathematical model combining bacterial swimming and the motion with flow. • Boundary layer plays a key role in bacterial attachment under flow condition. • The competition between bacterial swimming and the motion with flow is evaluated.« less

Calculation of steady and unsteady transonic flow using a Cartesian mesh and gridless boundary conditions with application to aeroelasticity

NASA Astrophysics Data System (ADS)

Kirshman, David

A numerical method for the solution of inviscid compressible flow using an array of embedded Cartesian meshes in conjunction with gridless surface boundary conditions is developed. The gridless boundary treatment is implemented by means of a least squares fitting of the conserved flux variables using a cloud of nodes in the vicinity of the surface geometry. The method allows for accurate treatment of the surface boundary conditions using a grid resolution an order of magnitude coarser than required of typical Cartesian approaches. Additionally, the method does not suffer from issues associated with thin body geometry or extremely fine cut cells near the body. Unlike some methods that consider a gridless (or "meshless") treatment throughout the entire domain, multi-grid acceleration can be effectively incorporated and issues associated with global conservation are alleviated. The "gridless" surface boundary condition provides for efficient and simple problem set up since definition of the body geometry is generated independently from the field mesh, and automatically incorporated into the field discretization of the domain. The applicability of the method is first demonstrated for steady flow of single and multi-element airfoil configurations. Using this method, comparisons with traditional body-fitted grid simulations reveal that steady flow solutions can be obtained accurately with minimal effort associated with grid generation. The method is then extended to unsteady flow predictions. In this application, flow field simulations for the prescribed oscillation of an airfoil indicate excellent agreement with experimental data. Furthermore, it is shown that the phase lag associated with shock oscillation is accurately predicted without the need for a deformable mesh. Lastly, the method is applied to the prediction of transonic flutter using a two-dimensional wing model, in which comparisons with moving mesh simulations yield nearly identical results. As a result, applicability of the method to transient and vibrating fluid-structure interaction problems is established in which the requirement for a deformable mesh is eliminated.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bühler, Stefan; Obrist, Dominik; Kleiser, Leonhard

We investigate numerically the effects of nozzle-exit flow conditions on the jet-flow development and the near-field sound at a diameter-based Reynolds number of Re{sub D} = 18 100 and Mach number Ma = 0.9. Our computational setup features the inclusion of a cylindrical nozzle which allows to establish a physical nozzle-exit flow and therefore well-defined initial jet-flow conditions. Within the nozzle, the flow is modeled by a potential flow core and a laminar, transitional, or developing turbulent boundary layer. The goal is to document and to compare the effects of the different jet inflows on the jet flow development and themore » sound radiation. For laminar and transitional boundary layers, transition to turbulence in the jet shear layer is governed by the development of Kelvin-Helmholtz instabilities. With the turbulent nozzle boundary layer, the jet flow development is characterized by a rapid changeover to a turbulent free shear layer within about one nozzle diameter. Sound pressure levels are strongly enhanced for laminar and transitional exit conditions compared to the turbulent case. However, a frequency and frequency-wavenumber analysis of the near-field pressure indicates that the dominant sound radiation characteristics remain largely unaffected. By applying a recently developed scaling procedure, we obtain a close match of the scaled near-field sound spectra for all nozzle-exit turbulence levels and also a reasonable agreement with experimental far-field data.« less

Assessment of applications of transport models on regional scale solute transport

NASA Astrophysics Data System (ADS)

Guo, Z.; Fogg, G. E.; Henri, C.; Pauloo, R.

2017-12-01

Regional scale transport models are needed to support the long-term evaluation of groundwater quality and to develop management strategies aiming to prevent serious groundwater degradation. The purpose of this study is to evaluate the capacity of previously-developed upscaling approaches to accurately describe main solute transport processes including the capture of late-time tails under changing boundary conditions. Advective-dispersive contaminant transport in a 3D heterogeneous domain was simulated and used as a reference solution. Equivalent transport under homogeneous flow conditions were then evaluated applying the Multi-Rate Mass Transfer (MRMT) model. The random walk particle tracking method was used for both heterogeneous and homogeneous-MRMT scenarios under steady state and transient conditions. The results indicate that the MRMT model can capture the tails satisfactorily for plume transported with ambient steady-state flow field. However, when boundary conditions change, the mass transfer model calibrated for transport under steady-state conditions cannot accurately reproduce the tailing effect observed for the heterogeneous scenario. The deteriorating impact of transient boundary conditions on the upscaled model is more significant for regions where flow fields are dramatically affected, highlighting the poor applicability of the MRMT approach for complex field settings. Accurately simulating mass in both mobile and immobile zones is critical to represent the transport process under transient flow conditions and will be the future focus of our study.

On an Asymptotically Consistent Unsteady Interacting Boundary Layer

NASA Technical Reports Server (NTRS)

Bartels, Robert E.

2007-01-01

This paper develops the asymptotic matching of an unsteady compressible boundary layer to an inviscid flow. Of particular importance is the velocity injection or transpiration boundary condition derived by this theory. It is found that in general the transpiration will contain a slope of the displacement thickness and a time derivative of a density integral. The conditions under which the second term may be neglected, and its consistency with the established results of interacting boundary layer are discussed.

Vorticity interaction effects on blunt bodies. [hypersonic viscous shock layers

NASA Technical Reports Server (NTRS)

Anderson, E. C.; Wilcox, D. C.

1977-01-01

Numerical solutions of the viscous shock layer equations governing laminar and turbulent flows of a perfect gas and radiating and nonradiating mixtures of perfect gases in chemical equilibrium are presented for hypersonic flow over spherically blunted cones and hyperboloids. Turbulent properties are described in terms of the classical mixing length. Results are compared with boundary layer and inviscid flowfield solutions; agreement with inviscid flowfield data is satisfactory. Agreement with boundary layer solutions is good except in regions of strong vorticity interaction; in these flow regions, the viscous shock layer solutions appear to be more satisfactory than the boundary layer solutions. Boundary conditions suitable for hypersonic viscous shock layers are devised for an advanced turbulence theory.

CFD study of the flow pattern in an ultrasonic horn reactor: Introducing a realistic vibrating boundary condition.

PubMed

Rahimi, Masoud; Movahedirad, Salman; Shahhosseini, Shahrokh

2017-03-01

Recently, great attention has been paid to predict the acoustic streaming field distribution inside the sonoreactors, induced by high-power ultrasonic wave generator. The focus of this paper is to model an ultrasonic vibrating horn and study the induced flow pattern with a newly developed moving boundary condition. The numerical simulation utilizes the modified cavitation model along with the "mixture" model for turbulent flow (RNG, k-ε), and a moving boundary condition with an oscillating parabolic-logarithmic profile, applied to the horn tip. This moving-boundary provides the situation in which the center of the horn tip vibrates stronger than that of the peripheral regions. The velocity field obtained by computational fluid dynamic was in a reasonably good agreement with the PIV results. The moving boundary model is more accurate since it better approximates the movement of the horn tip in the ultrasonic assisted process. From an optimizing point of view, the model with the new moving boundary is more suitable than the conventional models for design purposes because the displacement magnitude of the horn tip is the only fitting parameter. After developing and validating the numerical model, the model was utilized to predict various quantities such as cavitation zone, pressure field and stream function that are not experimentally feasible to measure. Copyright © 2016 Elsevier B.V. All rights reserved.

Simulation of confined magnetohydrodynamic flows with Dirichlet boundary conditions using a pseudo-spectral method with volume penalization

DOE Office of Scientific and Technical Information (OSTI.GOV)

Morales, Jorge A.; Leroy, Matthieu; Bos, Wouter J.T.

A volume penalization approach to simulate magnetohydrodynamic (MHD) flows in confined domains is presented. Here the incompressible visco-resistive MHD equations are solved using parallel pseudo-spectral solvers in Cartesian geometries. The volume penalization technique is an immersed boundary method which is characterized by a high flexibility for the geometry of the considered flow. In the present case, it allows to use other than periodic boundary conditions in a Fourier pseudo-spectral approach. The numerical method is validated and its convergence is assessed for two- and three-dimensional hydrodynamic (HD) and MHD flows, by comparing the numerical results with results from literature and analyticalmore » solutions. The test cases considered are two-dimensional Taylor–Couette flow, the z-pinch configuration, three dimensional Orszag–Tang flow, Ohmic-decay in a periodic cylinder, three-dimensional Taylor–Couette flow with and without axial magnetic field and three-dimensional Hartmann-instabilities in a cylinder with an imposed helical magnetic field. Finally, we present a magnetohydrodynamic flow simulation in toroidal geometry with non-symmetric cross section and imposing a helical magnetic field to illustrate the potential of the method.« less

A new ghost-node method for linking different models and initial investigations of heterogeneity and nonmatching grids

USGS Publications Warehouse

Dickinson, J.E.; James, S.C.; Mehl, S.; Hill, M.C.; Leake, S.A.; Zyvoloski, G.A.; Faunt, C.C.; Eddebbarh, A.-A.

2007-01-01

A flexible, robust method for linking parent (regional-scale) and child (local-scale) grids of locally refined models that use different numerical methods is developed based on a new, iterative ghost-node method. Tests are presented for two-dimensional and three-dimensional pumped systems that are homogeneous or that have simple heterogeneity. The parent and child grids are simulated using the block-centered finite-difference MODFLOW and control-volume finite-element FEHM models, respectively. The models are solved iteratively through head-dependent (child model) and specified-flow (parent model) boundary conditions. Boundary conditions for models with nonmatching grids or zones of different hydraulic conductivity are derived and tested against heads and flows from analytical or globally-refined models. Results indicate that for homogeneous two- and three-dimensional models with matched grids (integer number of child cells per parent cell), the new method is nearly as accurate as the coupling of two MODFLOW models using the shared-node method and, surprisingly, errors are slightly lower for nonmatching grids (noninteger number of child cells per parent cell). For heterogeneous three-dimensional systems, this paper compares two methods for each of the two sets of boundary conditions: external heads at head-dependent boundary conditions for the child model are calculated using bilinear interpolation or a Darcy-weighted interpolation; specified-flow boundary conditions for the parent model are calculated using model-grid or hydrogeologic-unit hydraulic conductivities. Results suggest that significantly more accurate heads and flows are produced when both Darcy-weighted interpolation and hydrogeologic-unit hydraulic conductivities are used, while the other methods produce larger errors at the boundary between the regional and local models. The tests suggest that, if posed correctly, the ghost-node method performs well. Additional testing is needed for highly heterogeneous systems. ?? 2007 Elsevier Ltd. All rights reserved.

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2014-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2015-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

A physical-based gas-surface interaction model for rarefied gas flow simulation

NASA Astrophysics Data System (ADS)

Liang, Tengfei; Li, Qi; Ye, Wenjing

2018-01-01

Empirical gas-surface interaction models, such as the Maxwell model and the Cercignani-Lampis model, are widely used as the boundary condition in rarefied gas flow simulations. The accuracy of these models in the prediction of macroscopic behavior of rarefied gas flows is less satisfactory in some cases especially the highly non-equilibrium ones. Molecular dynamics simulation can accurately resolve the gas-surface interaction process at atomic scale, and hence can predict accurate macroscopic behavior. They are however too computationally expensive to be applied in real problems. In this work, a statistical physical-based gas-surface interaction model, which complies with the basic relations of boundary condition, is developed based on the framework of the washboard model. In virtue of its physical basis, this new model is capable of capturing some important relations/trends for which the classic empirical models fail to model correctly. As such, the new model is much more accurate than the classic models, and in the meantime is more efficient than MD simulations. Therefore, it can serve as a more accurate and efficient boundary condition for rarefied gas flow simulations.

Evaluation of general non-reflecting boundary conditions for industrial CFD applications

NASA Astrophysics Data System (ADS)

Basara, Branislav; Frolov, Sergei; Lidskii, Boris; Posvyanskii, Vladimir

2007-11-01

The importance of having proper boundary conditions for the calculation domain is a known issue in Computational Fluid Dynamics (CFD). In many situations, it is very difficult to define a correct boundary condition. The flow may enter and leave the computational domain at the same time and at the same boundary. In such circumstances, it is important that numerical implementation of boundary conditions enforces certain physical constraints leading to correct results which then ensures a better convergence rate. The aim of this paper is to evaluate recently proposed non-reflecting boundary conditions (Frolov et al., 2001, Advances in Chemical Propulsion) on industrial CFD applications. Derivation of the local non-reflecting boundary conditions at the open boundary is based on finding the solution of linearized Euler equations vanishing at infinity for both incompressible and compressible formulations. This is implemented into the in-house CFD package AVL FIRE and some numerical details will be presented as well. The key applications in this paper are from automotive industry, e.g. an external car aerodynamics, an intake port, etc. The results will show benefits of using effective non-reflecting boundary conditions.

A High-Order Immersed Boundary Method for Acoustic Wave Scattering and Low-Mach Number Flow-Induced Sound in Complex Geometries

PubMed Central

Seo, Jung Hee; Mittal, Rajat

2010-01-01

A new sharp-interface immersed boundary method based approach for the computation of low-Mach number flow-induced sound around complex geometries is described. The underlying approach is based on a hydrodynamic/acoustic splitting technique where the incompressible flow is first computed using a second-order accurate immersed boundary solver. This is followed by the computation of sound using the linearized perturbed compressible equations (LPCE). The primary contribution of the current work is the development of a versatile, high-order accurate immersed boundary method for solving the LPCE in complex domains. This new method applies the boundary condition on the immersed boundary to a high-order by combining the ghost-cell approach with a weighted least-squares error method based on a high-order approximating polynomial. The method is validated for canonical acoustic wave scattering and flow-induced noise problems. Applications of this technique to relatively complex cases of practical interest are also presented. PMID:21318129

Experimental analysis of the flow in a two stage axial compressor at off-design conditions

NASA Astrophysics Data System (ADS)

Massardo, Aristide; Satta, Antonio

1987-05-01

The experimental analysis of the flow that develops in a two-stage axial flow compressor at off-design conditions is presented. The measurements are performed upstream, between, and downstream of the four blade rows of the compressor. The analysis shows the off-design effects on the local conditions of the flow field. Low-energy flow zones are identified, and the development of annulus-boundary-layer, secondary, and tip-clearance flows is shown. The tip-clearance flows are also present in the stator rows with various outlying conditions (stationary or rotating hub).

Entrainment of bed material by Earth-surface mass flows: review and reformulation of depth-integrated theory

USGS Publications Warehouse

Iverson, Richard M.; Chaojun Ouyang,

2015-01-01

Earth-surface mass flows such as debris flows, rock avalanches, and dam-break floods can grow greatly in size and destructive potential by entraining bed material they encounter. Increasing use of depth-integrated mass- and momentum-conservation equations to model these erosive flows motivates a review of the underlying theory. Our review indicates that many existing models apply depth-integrated conservation principles incorrectly, leading to spurious inferences about the role of mass and momentum exchanges at flow-bed boundaries. Model discrepancies can be rectified by analyzing conservation of mass and momentum in a two-layer system consisting of a moving upper layer and static lower layer. Our analysis shows that erosion or deposition rates at the interface between layers must in general satisfy three jump conditions. These conditions impose constraints on valid erosion formulas, and they help determine the correct forms of depth-integrated conservation equations. Two of the three jump conditions are closely analogous to Rankine-Hugoniot conditions that describe the behavior of shocks in compressible gasses, and the third jump condition describes shear traction discontinuities that necessarily exist across eroding boundaries. Grain-fluid mixtures commonly behave as compressible materials as they undergo entrainment, because changes in bulk density occur as the mixtures mobilize and merge with an overriding flow. If no bulk density change occurs, then only the shear-traction jump condition applies. Even for this special case, however, accurate formulation of depth-integrated momentum equations requires a clear distinction between boundary shear tractions that exist in the presence or absence of bed erosion.

Investigation of pumping mechanism for non-Newtonian blood flow with AC electrothermal forces in a microchannel by hybrid boundary element method and immersed boundary-lattice Boltzmann method.

PubMed

Ren, Qinlong

2018-02-10

Efficient pumping of blood flow in a microfluidic device is essential for rapid detection of bacterial bloodstream infections (BSI) using alternating current (AC) electrokinetics. Compared with AC electro-osmosis (ACEO) phenomenon, the advantage of AC electrothermal (ACET) mechanism is its capability of pumping biofluids with high electrical conductivities at a relatively high AC voltage frequency. In the current work, the microfluidic pumping of non-Newtonian blood flow using ACET forces is investigated in detail by modeling its multi-physics process with hybrid boundary element method (BEM) and immersed boundary-lattice Boltzmann method (IB-LBM). The Carreau-Yasuda model is used to simulate the realistic rheological behavior of blood flow. The ACET pumping efficiency of blood flow is studied in terms of different AC voltage magnitudes and frequencies, thermal boundary conditions of electrodes, electrode configurations, channel height, and the channel length per electrode pair. Besides, the effect of rheological behavior on the blood flow velocity is theoretically analyzed by comparing with the Newtonian fluid flow using scaling law analysis under the same physical conditions. The results indicate that the rheological behavior of blood flow and its frequency-dependent dielectric property make the pumping phenomenon of blood flow different from that of the common Newtonian aqueous solutions. It is also demonstrated that using a thermally insulated electrode could enhance the pumping efficiency dramatically. Besides, the results conclude that increasing the AC voltage magnitude is a more economical pumping approach than adding the number of electrodes with the same energy consumption when the Joule heating effect is acceptable. © 2018 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

The `Henry Problem' of `density-driven' groundwater flow versus Tothian `groundwater flow systems' with variable density: A review of the influential Biscayne aquifer data.

NASA Astrophysics Data System (ADS)

Weyer, K. U.

2017-12-01

Coastal groundwater flow investigations at the Biscayne Bay, south of Miami, Florida, gave rise to the concept of density-driven flow of seawater into coastal aquifers creating a saltwater wedge. Within that wedge, convection-driven return flow of seawater and a dispersion zone were assumed by Cooper et al. (1964) to be the cause of the Biscayne aquifer `sea water wedge'. This conclusion was based on the chloride distribution within the aquifer and on an analytical model concept assuming convection flow within a confined aquifer without taking non-chemical field data into consideration. This concept was later labelled the `Henry Problem', which any numerical variable density flow program must be able to simulate to be considered acceptable. Both, `density-driven flow' and Tothian `groundwater flow systems' (with or without variable density conditions) are driven by gravitation. The difference between the two are the boundary conditions. 'Density-driven flow' occurs under hydrostatic boundary conditions while Tothian `groundwater flow systems' occur under hydrodynamic boundary conditions. Revisiting the Cooper et al. (1964) publication with its record of piezometric field data (heads) showed that the so-called sea water wedge has been caused by discharging deep saline groundwater driven by gravitational flow and not by denser sea water. Density driven flow of seawater into the aquifer was not found reflected in the head measurements for low and high tide conditions which had been taken contemporaneously with the chloride measurements. These head measurements had not been included in the flow interpretation. The very same head measurements indicated a clear dividing line between shallow local fresh groundwater flow and saline deep groundwater flow without the existence of a dispersion zone or a convection cell. The Biscayne situation emphasizes the need for any chemical interpretation of flow pattern to be supported by head data as energy indicators of flow fields. At the Biscayne site density-driven flow of seawater did and does not exist. Instead this site and the Florida coast line in general are the end points of local fresh and regional saline groundwater flow systems driven by gravity forces and not by density differences.

Heat and mass transfer boundary conditions at the surface of a heated sessile droplet

NASA Astrophysics Data System (ADS)

Ljung, Anna-Lena; Lundström, T. Staffan

2017-12-01

This work numerically investigates how the boundary conditions of a heated sessile water droplet should be defined in order to include effects of both ambient and internal flow. Significance of water vapor, Marangoni convection, separate simulations of the external and internal flow, and influence of contact angle throughout drying is studied. The quasi-steady simulations are carried out with Computational Fluid Dynamics and conduction, natural convection and Marangoni convection are accounted for inside the droplet. For the studied conditions, a noticeable effect of buoyancy due to evaporation is observed. Hence, the inclusion of moisture increases the maximum velocities in the external flow. Marangoni convection will, in its turn, increase the velocity within the droplet with up to three orders of magnitude. Results furthermore show that the internal and ambient flow can be simulated separately for the conditions studied, and the accuracy is improved if the internal temperature gradient is low, e.g. if Marangoni convection is present. Simultaneous simulations of the domains are however preferred at high plate temperatures if both internal and external flows are dominated by buoyancy and natural convection. The importance of a spatially resolved heat and mass transfer boundary condition is, in its turn, increased if the internal velocity is small or if there is a large variation of the transfer coefficients at the surface. Finally, the results indicate that when the internal convective heat transport is small, a rather constant evaporation rate may be obtained throughout the drying at certain conditions.

Infrared Images of Boundary Layer Transition on the D8 Transport Configuration in the LaRC 14- by 22-Foot Subsonic Tunnel

NASA Technical Reports Server (NTRS)

Mason, Michelle L.; Gatlin, Gregory M.

2015-01-01

Grit, trip tape, or trip dots are routinely applied on the leading-edge regions of the fuselage, wings, tails or nacelles of wind tunnel models to trip the flow from laminar to turbulent. The thickness of the model's boundary layer is calculated for nominal conditions in the wind tunnel test to determine the effective size of the trip dots, but the flow over the model may not transition as intended for runs with different flow conditions. Temperature gradients measured with an infrared camera can be used to detect laminar to turbulent boundary layer transition on a wind tunnel model. This non-intrusive technique was used in the NASA Langley 14- by 22-Foot Subsonic Tunnel to visualize the behavior of the flow over a D8 transport configuration model. As the flow through the wind tunnel either increased to or decreased from the run conditions, a sufficient temperature difference existed between the air and the model to visualize the transition location (due to different heat transfer rates through the laminar and the turbulent boundary layers) for several runs in this test. Transition phenomena were visible without active temperature control in the atmospheric wind tunnel, whether the air was cooler than the model or vice-versa. However, when the temperature of the model relative to the air was purposely changed, the ability to detect transition in the infrared images was enhanced. Flow characteristics such as a wing root horseshoe vortex or the presence of fore-body vortical flows also were observed in the infrared images. The images of flow features obtained for this study demonstrate the usefulness of current infrared technology in subsonic wind tunnel tests.

Boundary Layer Transition in the Leading Edge Region of a Swept Cylinder in High Speed Flow

NASA Technical Reports Server (NTRS)

Coleman, Colin P.

1998-01-01

Experiments were conducted on a 76 degree swept cylinder to establish the behavior of the attachment line transition process in a low-disturbance level, Mach number 1.6 flow. For a near adiabatic wall condition, the attachment-line boundary layer remained laminar up to the highest attainable Reynolds number. The attachment-line boundary layer transition under the influence of trip wires depended on wind tunnel disturbance level, and a transition onset condition for this flow is established. Internal heating raised the surface temperature of the attachment line to induce boundary layer instabilities. This was demonstrated experimentally for the first time and the frequencies of the most amplified disturbances were determined over a range of temperature settings. Results were in excellent agreement to those predicted by a linear stability code, and provide the first experimental verification of theory. Transition onset along the heated attachment line at an R-bar of 800 under quiet tunnel conditions was found to correlate with an N factor of 13.2. Increased tunnel disturbance levels caused the transition onset to occur at lower cylinder surface temperatures and was found to correlate with an approximate N factor of 1 1.9, so demonstrating that the attachment-line boundary layer is receptive to increases in the tunnel disturbance level.

Shock wave boundary layer interaction on suction side of compressor profile in single passage test section

NASA Astrophysics Data System (ADS)

Flaszynski, Pawel; Doerffer, Piotr; Szwaba, Ryszard; Kaczynski, Piotr; Piotrowicz, Michal

2015-11-01

The shock wave boundary layer interaction on the suction side of transonic compressor blade is one of the main objectives of TFAST project (Transition Location Effect on Shock Wave Boundary Layer Interaction). In order to investigate the flow structure on the suction side of a profile, a design of a generic test section in linear transonic wind tunnel was proposed. The experimental and numerical results for the flow structure investigations are shown for the flow conditions as the existing ones on the suction side of the compressor profile. Near the sidewalls the suction slots are applied for the corner flow structure control. It allows to control the Axial Velocity Density Ratio (AVDR), important parameter for compressor cascade investigations. Numerical results for Explicit Algebraic Reynolds Stress Model with transition modeling are compared with oil flow visualization, schlieren and Pressure Sensitive Paint. Boundary layer transition location is detected by Temperature Sensitive Paint.

Boundary layers and global stability of laboratory quasi-Keplerian flow

NASA Astrophysics Data System (ADS)

Edlund, E. M.; Ji, H.

2013-11-01

Studies in the HTX device at PPPL, a modified Taylor-Couette experiment, have demonstrated a robust stability of astrophysically relevant, quasi-Keplerian flows. Independent rings on the axial boundary can be used to fine tune the rotation profile, allowing ideal Couette rotation to be achieved over nearly the entire radial gap. Fluctuation levels in these flows are observed to be at nearly the noise floor of the laser Doppler velocimetry (LDV) diagnostic, in agreement with prior studies under similar conditions. Deviations from optimal operating parameters illustrate the importance of centrifugally unstable boundary layers in Taylor-Couette devices of the classical configuration where the axial boundaries rotate with the outer cylinder. The global stability of nearly ideal-Couette flows, with implications for astrophysical systems, will be discussed in light of the global stability of these flows with respect to externally applied perturbations of large magnitude.

Improved Boundary Conditions for Cell-centered Difference Schemes

NASA Technical Reports Server (NTRS)

VanderWijngaart, Rob F.; Klopfer, Goetz H.; Chancellor, Marisa K. (Technical Monitor)

1997-01-01

Cell-centered finite-volume (CCFV) schemes have certain attractive properties for the solution of the equations governing compressible fluid flow. Among others, they provide a natural vehicle for specifying flux conditions at the boundaries of the physical domain. Unfortunately, they lead to slow convergence for numerical programs utilizing them. In this report a method for investigating and improving the convergence of CCFV schemes is presented, which focuses on the effect of the numerical boundary conditions. The key to the method is the computation of the spectral radius of the iteration matrix of the entire demoralized system of equations, not just of the interior point scheme or the boundary conditions.

Assessment of a 3-D boundary layer code to predict heat transfer and flow losses in a turbine

NASA Technical Reports Server (NTRS)

Anderson, O. L.

1984-01-01

Zonal concepts are utilized to delineate regions of application of three-dimensional boundary layer (DBL) theory. The zonal approach requires three distinct analyses. A modified version of the 3-DBL code named TABLET is used to analyze the boundary layer flow. This modified code solves the finite difference form of the compressible 3-DBL equations in a nonorthogonal surface coordinate system which includes coriolis forces produced by coordinate rotation. These equations are solved using an efficient, implicit, fully coupled finite difference procedure. The nonorthogonal surface coordinate system is calculated using a general analysis based on the transfinite mapping of Gordon which is valid for any arbitrary surface. Experimental data is used to determine the boundary layer edge conditions. The boundary layer edge conditions are determined by integrating the boundary layer edge equations, which are the Euler equations at the edge of the boundary layer, using the known experimental wall pressure distribution. Starting solutions along the inflow boundaries are estimated by solving the appropriate limiting form of the 3-DBL equations.

A nonperturbing boundary-layer transition detection

NASA Astrophysics Data System (ADS)

Ohare, J. E.

1985-01-01

A laser interferometer technique is being applied to the characterization of boundary-layer conditions on models in supersonic and hypersonic wind tunnels in the von Karman Facility at Arnold Engineering Development Center (AEDC). The Boundary-Layer Transition Detector (BLTD), based on lateral interferometry, is applicable for determining the turbulence frequency spectrum of boundary layers in compressible flow. The turbulence, in terms of air density fluctuations, is detected by monitoring interferometric fringe phase shifts (in real time) formed by one beam which passes through the boundary layer and a reference beam which is outside the boundary layer. This technique is nonintrusive to the flow field unlike other commonly used methods such as pitot tube probing and hot-wire anemometry. Model boundary-layer data are presented at Mach 8 and compared with data recorded using other methods during boundary-layer transition from laminar to turbulent flow. Spectra from the BLTD reveal the presence of a high-frequency peak during transition, which is characteristic of spectra obtained with hot wires. The BLTD is described along with operational requirements and limitations.

A Nonperturbing Boundary-Layer Transition Detector

NASA Astrophysics Data System (ADS)

O'Hare, J. E.

1986-01-01

A laser interferometer technique is being applied to the characterization of boundary-layer conditions on models in supersonic and hypersonic wind tunnels in the von Kaman Facility at Arnold Engineering Development Center (AEDC). The Boundary-Layer Transition Detector (BLTD), based on lateral interferometry, is applicable for determining the turbulence frequency spectrum of boundary layers in compressible flow. The turbulence, in terms of air density fluctuations, is detected by monitoring interferometric fringe phase shifts (in real time) formed by one beam which passes through the boundary layer and a reference beam which is outside the boundary layer. This technique is nonintrusive to the flow field unlike other commonly used methods such as pitot tube probing and hot-wire anemometry. Model boundary-layer data are presented at Mach 8 and compared with data recorded using other methods during boundary-layer transition from laminar to turbulent flow. Spectra from the BLTD reveal the presence of a high-frequency peak during transition, which is characteristic of spectra obtained with hot wires. The BLTD is described along with operational requirements and limitations.

Re-Innovating Recycling for Turbulent Boundary Layer Simulations

NASA Astrophysics Data System (ADS)

Ruan, Joseph; Blanquart, Guillaume

2017-11-01

Historically, turbulent boundary layers along a flat plate have been expensive to simulate numerically, in part due to the difficulty of initializing the inflow with ``realistic'' turbulence, but also due to boundary layer growth. The former has been resolved in several ways, primarily dedicating a region of at least 10 boundary layer thicknesses in width to rescale and recycle flow or by extending the region far enough downstream to allow a laminar flow to develop into turbulence. Both of these methods are relatively costly. We propose a new method to remove the need for an inflow region, thus reducing computational costs significantly. Leveraging the scale similarity of the mean flow profiles, we introduce a coordinate transformation so that the boundary layer problem can be solved as a parallel flow problem with additional source terms. The solutions in the new coordinate system are statistically homogeneous in the downstream direction and so the problem can be solved with periodic boundary conditions. The present study shows the stability of this method, its implementation and its validation for a few laminar and turbulent boundary layer cases.

A Lagrangian Analysis of a Developing and Non-Developing Disturbance Observed During the PREDICT Experiment

DTIC Science & Technology

2012-12-03

paper provides an introduction of Lagrangian techniques for locating flow boundaries that encompass regions of recirculation in time- dependent flows...the low- to mid- level embryonic vortex from adverse conditions, while the 1The glossary on NOAA’s Hurricane Research Division’s web - site uses...wave or disturbance. This paper provides an introduction of Lagrangian techniques for locating flow boundaries that encompass regions of recirculation

The NASA Lewis Research Center Internal Fluid Mechanics Facility

NASA Technical Reports Server (NTRS)

Porro, A. R.; Hingst, W. R.; Wasserbauer, C. A.; Andrews, T. B.

1991-01-01

An experimental facility specifically designed to investigate internal fluid duct flows is described. It is built in a modular fashion so that a variety of internal flow test hardware can be installed in the facility with minimal facility reconfiguration. The facility and test hardware interfaces are discussed along with design constraints of future test hardware. The plenum flow conditioning approach is also detailed. Available instrumentation and data acquisition capabilities are discussed. The incoming flow quality was documented over the current facility operating range. The incoming flow produces well behaved turbulent boundary layers with a uniform core. For the calibration duct used, the boundary layers approached 10 percent of the duct radius. Freestream turbulence levels at the various operating conditions varied from 0.64 to 0.69 percent of the average freestream velocity.

Numerical calculation of the parameters of the efflux from a helium dewar used for cooling of heat shields in a satellite

NASA Technical Reports Server (NTRS)

Brendley, K.; Chato, J. C.

1982-01-01

The parameters of the efflux from a helium dewar in space were numerically calculated. The flow was modeled as a one dimensional compressible ideal gas with variable properties. The primary boundary conditions are flow with friction and flow with heat transfer and friction. Two PASCAL programs were developed to calculate the efflux parameters: EFFLUZD and EFFLUXM. EFFLUXD calculates the minimum mass flow for the given shield temperatures and shield heat inputs. It then calculates the pipe lengths, diameter, and fluid parameters which satisfy all boundary conditions. Since the diameter returned by EFFLUXD is only rarely of nominal size, EFFLUXM calculates the mass flow and shield heat exchange for given pipe lengths, diameter, and shield temperatures.

Darcy-Forchheimer Three-Dimensional Flow of Williamson Nanofluid over a Convectively Heated Nonlinear Stretching Surface

NASA Astrophysics Data System (ADS)

Hayat, Tasawar; Aziz, Arsalan; Muhammad, Taseer; Alsaedi, Ahmed

2017-09-01

The present study elaborates three-dimensional flow of Williamson nanoliquid over a nonlinear stretchable surface. Fluid flow obeys Darcy-Forchheimer porous medium. A bidirectional nonlinear stretching surface generates the flow. Convective surface condition of heat transfer is taken into consideration. Further the zero nanoparticles mass flux condition is imposed at the boundary. Effects of thermophoresis and Brownian diffusion are considered. Assumption of boundary layer has been employed in the problem formulation. Convergent series solutions for the nonlinear governing system are established through the optimal homotopy analysis method (OHAM). Graphs have been sketched in order to analyze that how the velocity, temperature and concentration distributions are affected by distinct emerging flow parameters. Skin friction coefficients and local Nusselt number are also computed and discussed.

On the Physics of Flow Separation Along a Low Pressure Turbine Blade Under Unsteady Flow Conditions

NASA Technical Reports Server (NTRS)

Schobeiri, Meinhard T.; Ozturk, Burak; Ashpis, David E.

2005-01-01

The present study, which is the first of a series of investigations dealing with specific issues of low pressure turbine (LPT) boundary layer aerodynamics, is aimed at providing detailed unsteady boundary flow information to understand the underlying physics of the inception, onset, and extent of the separation zone. A detailed experimental study on the behavior of the separation zone on the suction surface of a highly loaded LPT-blade under periodic unsteady wake flow is presented. Experimental investigations were performed at Texas A&M Turbomachinery Performance and Flow Research Laboratory using a large-scale unsteady turbine cascade research facility with an integrated wake generator and test section unit. To account for a high flow deflection of LPT-cascades at design and off-design operating points, the entire wake generator and test section unit including the traversing system is designed to allow a precise angle adjustment of the cascade relative to the incoming flow. This is done by a hydraulic platform, which simultaneously lifts and rotates the wake generator and test section unit. The unit is then attached to the tunnel exit nozzle with an angular accuracy of better than 0.05 , which is measured electronically. Utilizing a Reynolds number of 110,000 based on the blade suction surface length and the exit velocity, one steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, wake velocities and turbulence intensities are investigated using hot-wire anemometry. In addition to the unsteady boundary layer measurements, blade surface pressure measurements were performed at Re=50,000, 75,000, 100,000, and 125,000 at one steady and two periodic unsteady inlet flow conditions. Detailed unsteady boundary layer measurement identifies the onset and extent of the separation zone as well as its behavior under unsteady wake flow. The results presented in ensemble-averaged and contour plot forms contribute to understanding the physics of the separation phenomenon under periodic unsteady wake flow. Several physical mechanisms are discussed.

Effect of Boundary Conditions on Numerically Simulated Tornado-like Vortices.

NASA Astrophysics Data System (ADS)

Smith, David R.

1987-02-01

The boundary conditions for Rotunno's numerical model which simulates tornado-like vortices are examined. In particular, the lateral boundary condition for tangential velocity and the upper boundary condition for radial and tangential velocities are considered to determine if they have any significant impact on vortex development.The choice of the lateral boundary condition did not appear to have any real effect on the development of the vortex over the range of swirl ratios studied (0.87-2.61).The upper boundary conditions attempt to simulate both the presence and absence of the flow-straightening baffle. The boundary condition corresponding to the baffle in place produced a distinct boundary layer in the u and v field and very strong upflow and downflow within the vortex core. When this condition is removed, there is both radial and tangential motion throughout the domain and a reduction of the vertical velocity. At small swirl ratio (S = 0.87) this boundary condition has a profound impact on the narrow vortex, producing changes in the pressure field that intensifies the vortex. At higher swirl ratio the vortex is apparently broad enough to better adjust to the changes of the upper boundary condition and, thus, experiences little change in the development of the vortex.

Influence of spatial variation of phenomenological parameters on the modeling of boundary conditions for flows with dynamic wetting

NASA Astrophysics Data System (ADS)

Hizumi, Yuka; Omori, Takeshi; Yamaguchi, Yasutaka; Kajisima, Takeo

2014-11-01

For reliable prediction of multiphase flows in micro- and nano-scales, continuum models are expected to account for small scale physics near the contact line (CL) region. Some existing works (for example the series of papers by the group of Qian and Ren) have been successful in deriving continuum models and corresponding boundary conditions which reproduce well the molecular dynamics (MD) simulation results. Their studies, however, did not fully address the issue of adsorption layer especially in the CL region, and it is still not clear if general conclusion can be deduced from their results. In the present study we investigate in detail the local viscosity and the corresponding stress tensor formulation in the solid-liquid interface and in the CL region of immiscible two-phase Couette flows by means of MD simulation. The application limit of the generalized Navier boundary condition and the continuum model with uniform viscosity is addressed by systematic coarse-graining of sampling bins.

Investigation of blood flow in the external carotid artery and its branches with a new 0D peripheral model.

PubMed

Ohhara, Yoshihito; Oshima, Marie; Iwai, Toshinori; Kitajima, Hiroaki; Yajima, Yasuharu; Mitsudo, Kenji; Krdy, Absy; Tohnai, Iwai

2016-02-04

Patient-specific modelling in clinical studies requires a realistic simulation to be performed within a reasonable computational time. The aim of this study was to develop simple but realistic outflow boundary conditions for patient-specific blood flow simulation which can be used to clarify the distribution of the anticancer agent in intra-arterial chemotherapy for oral cancer. In this study, the boundary conditions are expressed as a zero dimension (0D) resistance model of the peripheral vessel network based on the fractal characteristics of branching arteries combined with knowledge of the circulatory system and the energy minimization principle. This resistance model was applied to four patient-specific blood flow simulations at the region where the common carotid artery bifurcates into the internal and external carotid arteries. Results of these simulations with the proposed boundary conditions were compared with the results of ultrasound measurements for the same patients. The pressure was found to be within the physiological range. The difference in velocity in the superficial temporal artery results in an error of 5.21 ± 0.78 % between the numerical results and the measurement data. The proposed outflow boundary conditions, therefore, constitute a simple resistance-based model and can be used for performing accurate simulations with commercial fluid dynamics software.

A stable penalty method for the compressible Navier-Stokes equations. 1: Open boundary conditions

NASA Technical Reports Server (NTRS)

Hesthaven, J. S.; Gottlieb, D.

1994-01-01

The purpose of this paper is to present asymptotically stable open boundary conditions for the numerical approximation of the compressible Navier-Stokes equations in three spatial dimensions. The treatment uses the conservation form of the Navier-Stokes equations and utilizes linearization and localization at the boundaries based on these variables. The proposed boundary conditions are applied through a penalty procedure, thus ensuring correct behavior of the scheme as the Reynolds number tends to infinity. The versatility of this method is demonstrated for the problem of a compressible flow past a circular cylinder.

Boundary streaming with Navier boundary condition.

PubMed

Xie, Jin-Han; Vanneste, Jacques

2014-06-01

In microfluidic applications involving high-frequency acoustic waves over a solid boundary, the Stokes boundary-layer thickness δ is so small that some non-negligible slip may occur at the fluid-solid interface. This paper assesses the impact of this slip by revisiting the classical problem of steady acoustic streaming over a flat boundary, replacing the no-slip boundary condition with the Navier condition u|_{y=0}=L_{s}∂_{y}u|_{y=0}, where u is the velocity tangent to the boundary y=0, and the parameter L_{s} is the slip length. A general expression is obtained for the streaming velocity across the boundary layer as a function of the dimensionless parameter L_{s}/δ. The limit outside the boundary layer provides an effective slip velocity satisfied by the interior mean flow. Particularizing to traveling and standing waves shows that the boundary slip respectively increases and decreases the streaming velocity.

Higher modes of the Orr-Sommerfeld problem for boundary layer flows

NASA Technical Reports Server (NTRS)

Lakin, W. D.; Grosch, C. E.

1983-01-01

The discrete spectrum of the Orr-Sommerfeld problem of hydrodynamic stability for boundary layer flows in semi-infinite regions is examined. Related questions concerning the continuous spectrum are also addressed. Emphasis is placed on the stability problem for the Blasius boundary layer profile. A general theoretical result is given which proves that the discrete spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 0,0) has only a finite number of discrete modes when U(y) has derivatives of all orders. Details are given of a highly accurate numerical technique based on collocation with splines for the calculation of stability characteristics. The technique includes replacement of 'outer' boundary conditions by asymptotic forms based on the proper large parameter in the stability problem. Implementation of the asymptotic boundary conditions is such that there is no need to make apriori distinctions between subcases of the discrete spectrum or between the discrete and continuous spectrums. Typical calculations for the usual Blasius problem are presented.

Dust transportation in bounday layers on complex areas

NASA Astrophysics Data System (ADS)

Karelsky, Kirill; Petrosyan, Arakel

2017-04-01

This presentation is aimed at creating and realization of new physical model of impurity transfer (solid particles and heavy gases) in areas with non-flat and/or nonstationary boundaries. The main idea of suggested method is to use non-viscous equations for solid particles transport modeling in the vicinity of complex boundary. In viscous atmosphere with as small as one likes coefficient of molecular viscosity, the non-slip boundary condition on solid surface must be observed. This postulates the reduction of velocity to zero at a solid surface. It is unconditionally in this case Prandtle hypothesis must be observed: for rather wide range of conditions in the surface neighboring layers energy dissipation of atmosphere flows is comparable by magnitude with manifestation of inertia forces. That is why according to Prandtle hypothesis in atmosphere movement characterizing by a high Reynolds number the boundary layer is forming near a planet surface, within which the required transition from zero velocities at the surface to magnitudes at the external boundary of the layer that are quite close to ones in ideal atmosphere flow. In that layer fast velocity gradients cause viscous effects to be comparable in magnitude with inertia forces influence. For conditions considered essential changes of hydrodynamic fields near solid boundary caused not only by nonslip condition but also by a various relief of surface: mountains, street canyons, individual buildings. Transport of solid particles, their ascent and precipitation also result in dramatic changes of meteorological fields. As dynamic processes of solid particles transfer accompanying the flow past of complex relief surface by wind flows is of our main interest we are to use equations of non-viscous hydrodynamic. We should put up with on the one hand idea of high wind gradients in the boundary layer and on the other hand disregard of molecular viscosity in two-phase atmosphere equations. We deal with describing high field gradients with the aid of scheme viscosity of numerical algorithm used to model near-surface phenomena. This idea is implemented in the model of ideal gas equations with variable equation of state describing particulates transportation within boundary layer with obstacles.

Chemical nonequilibrium Navier-Stokes solutions for hypersonic flow over an ablating graphite nosetip

NASA Technical Reports Server (NTRS)

Chen, Y. K.; Henline, W. D.

1993-01-01

The general boundary conditions including mass and energy balances of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver for predictions of the flow field, surface temperature, and surface ablation rates over re-entry space vehicles with ablating Thermal Protection Systems (TPS). The Navier-Stokes solver with general surface thermochemistry boundary conditions can predict more realistic solutions and provide useful information for the design of TPS. A test case with a proposed hypersonic test vehicle configuration and associated free stream conditions was developed. Solutions with various surface boundary conditions were obtained, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined. The solutions of the GASP code with complete ablating surface conditions were compared with those of the ASC code. The direction of future work is also discussed.

Computer simulations of rapid granular flows of spheres interacting with a flat, frictional boundary

DOE Office of Scientific and Technical Information (OSTI.GOV)

Louge, M.Y.

This paper employs computer simulations to test the theory of Jenkins [J. Applied Mech. [bold 59], 120 (1992)] for the interaction between a rapid granular flow of spheres and a flat, frictional wall. This paper examines the boundary conditions that relate the shear stress and energy flux at the wall to the normal stress, slip velocity, and fluctuation energy, and to the parameters that characterize a collision. It is found that while the theory captures the trends of the boundary conditions at low friction, it does not anticipate their behavior at large friction. A critical evaluation of Jenkins' assumptions suggestsmore » where his theory may be improved.« less

Subauroral polarization stream on the outer boundary of the ring current during an energetic ion injection event

NASA Astrophysics Data System (ADS)

Yuan, Zhigang; Qiao, Zheng; Li, Haimeng; Huang, Shiyong; Wang, Dedong; Yu, Xiongdong; Yu, Tao

2017-04-01

Subauroral polarization stream (SAPS) electric field can play an important role in the coupling between the inner magnetosphere and ionosphere; however, the production mechanism of SAPS has not been yet solved. During an energetic ion injection event on 26 March 2004, at latitudes lower than the equatorward boundaries of precipitating plasma sheet electrons and ions, the Defense Meteorological Satellite Program (DMSP) F13 satellite simultaneously observed a strong SAPS with the peak velocity of 1294 m/s and downward flowing field-aligned currents (FACs). Conjugate observations of DMSP F13 and NOAA 15 satellites have shown that FACs flowing into the ionosphere just lie in the outer boundary of the ring current (RC). The downward flowing FACs were observed in a region of positive latitudinal gradients of the ion energy density, implying that the downward flowing FACs are more likely linked to the azimuthal gradient than the radial gradient of the RC ion pressure. Our result demonstrates that RC ion pressure gradients on the outer boundary of the RC in the evening sector during energetic ion injection events can lead to downward flowing FACs so as to cause strong SAPS in condition of low ionospheric conductivities.