Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two dimensional flow”

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

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1955

On the Kernel function of the integral equation relating the lift and downwash distributions of oscillating finite wings in subsonic flow

This report treats the Kernel function of an integral equation that relates a known prescribed downwash distribution to an unknown lift distribution for a harmonically oscillating finite wing in compressible subsonic flow. The Kernel function is reduced to a form that can be accurately evaluated by separating the Kernel function into two parts: a part in which the singularities are isolated and analytically expressed and a nonsingular part which may be tabulated. The form of the Kernel function for the sonic case (Mach number 1) is treated separately. In addition, results for the special cases of Mach number of 0 (incompressible case) and frequency of 0 (steady case) are given. The derivation of the integral equation which involves this Kernel function is reproduced as an appendix. Another appendix gives the reduction of the form of the Kernel function obtained herein for the three-dimensional case to a known result of Possio for two-dimensional flow. A third appendix contains some remarks on the evaluation of the Kernel function, and a fourth appendix presents an alternate form of expression for the Kernel function.

Watkins, Charles E↗
NASA NTRS1994

Effect of Reynolds number and turbulence on airfoil aerodynamics at -90-degree incidence

A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Stremel, Paul M.↗
NASA NTRS1985

Flight-measured laminar boundary-layer transition phenomena including stability theory analysis

Flight experiments were conducted on a single-engine turboprop aircraft fitted with a 92-in-chord, 3-ft-span natural laminar flow glove at glove section lift coefficients from 0.15 to 1.10. The boundary-layer transition measurement methods used included sublimating chemicals and surface hot-film sensors. Transition occurred downstream of the minimum pressure point. Hot-film sensors provided a well-defined indication of laminar, laminar-separation, transitional, and turbulent boundary layers. Theoretical calculations of the boundary-layer parameters provided close agreement between the predicted laminar-separation point and the measured transition location. Tollmien-Schlichting (T-S) wave growth n-factors between 15 and 17 were calculated at the predicted point of laminar separation. These results suggest that for many practical airplane cruise conditions, laminar separation (as opposed to T-S instability) is the major cause of transition in predominantly two-dimensional flows.

Obara, C. J.↗
NASA NTRS1985

On the use of wall functions as boundary conditions for two-dimensional separated compressible flows

A new and improved wall function method for compressible turbulent flows has been developed and tested. This method is applicable to attached and separated flows, to both high- and low-Reynolds number flows, and to flows with adiabatic and nonadiabatic surfaces. This wall function method has been applied to the Launder-Spalding k-epsilon two-equation model of turbulence. The tests consist of comparisons of calculated and experimental results for: (1) an axisymmetrical transonic shock-wave/boundary-wave interaction flow at low Reynolds number in an adiabatic tube, (2) an axisymmetrical high-Reynolds number transonic flow over a nonadiabatic bump, and (3) a two-dimensional supersonic high-Reynolds number flow on a nonadiabatic deflected flap. Each of these experiments had significant regions of flow separation. The calculations are performed with an implicit algorithm that solves the Reynolds-averaged Navier-Stokes equations. It is shown that the results obtained agree very well with the data for the complex compressible flows tested.

Viegas, J. R.↗
DOE OSTI2019

Probabilistic Forecast of Single–Phase Flow in Porous Media With Uncertain Properties

Uncertainty about geologic makeup and properties of the subsurface renders inadequate a unique quantitative prediction of flow and transport. Instead, multiple alternative scenarios have to be explored within the probabilistic framework, typically by means of Monte Carlo simulations (MCS). These can be computationally expensive, and often prohibitively so, especially when the goal is to compute the tails of a distribution, i.e., probabilities of rare events, which are necessary for risk assessment and decision making under uncertainty. We deploy the method of distributions to derive a deterministic equation for the cumulative distribution function (CDF) of hydraulic head in an aquifer with uncertain (random) hydraulic conductivity. The CDF equation relies on a self-consistent closure approximation, which ensures that the resulting CDF of hydraulic head has the same mean and variance as those computed with either statistical moment equation (the approach used herein) or MCS. We conduct a series of numerical experiments dealing with steady-state two-dimensional flow driven by either a natural hydraulic head gradient or a pumping well. Furthermore, these experiments reveal that the CDF method remains accurate and robust for highly heterogeneous formations with the variance of log conductivity as large as five. For the same accuracy, it is also up to four orders of magnitude faster than MCS in computing hydraulic head with a required degree of confidence (probability).

42 ENGINEERING↗
DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗
NASA NTRSDate not supplied

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗
DOE OSTI2017

Lubricated immersed boundary method in two dimensions

Many biological examples of fluid-structure interaction, including the transit of red blood cells through the narrow slits in the spleen and the intracellular trafficking of vesicles into dendritic spines, involve the near-contact of elastic structures separated by thin layers of fluid. Motivated by such problems, we introduce an immersed boundary method that uses elements of lubrication theory to resolve thin fluid layers between immersed boundaries. We demonstrate 2nd-order accurate convergence for simple two-dimensional flows with known exact solutions to showcase the increased accuracy of this method compared to the standard immersed boundary method. Motivated by the phenomenon of wall-induced migration, we apply the lubricated immersed boundary method to simulate an elastic vesicle near a wall in shear flow. We also simulate the dynamics of a vesicle traveling through a narrow channel and observe the ability of the lubricated method to capture the vesicle motion on relatively coarse fluid grids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Consensus Equilibrium for Subsurface Delineation

Heterogeneity and insufficient site characterization limit our knowledge of the subsurface. Inversion techniques, which minimize the mismatch between observations and model predictions, have become an essential tool of subsurface characterization. Most optimization-based approaches fail to incorporate various implicit priors and capture the geological complexity. We overcome these limitations by deploying the plug-and-play and consensus equilibrium (CE) strategies, which provide a flexible framework for image reconstruction. Our CE methodology for spatial delineation of geologic formations consists of an image denoiser and a variational auto-encoder (deep learning-based emulator). The former ameliorates the reconstruction noise, yielding well-defined geological structures; its mathematical equivalence with the proximal operator allows the deployment of advanced denoisers (e.g., CNN-based denoiser) that do not correspond to a regularization objective. The latter defines a geology prior that imposes a geological constraint, for example, continuity and shape of geological features, onto the reconstructed image. Here, we conduct a series of numerical experiments dealing with transient two-dimensional flow driven by a pumping well and natural hydraulic head gradient. They demonstrate the CE framework's ability to delineate, both probabilistically and deterministically, complex subsurface environments with sufficient quality.

42 ENGINEERING↗
DOE OSTI2019

Non-Intrusive Inference Reduced Order Model for Fluids Using Deep Multistep Neural Network

In this effort we propose a data-driven learning framework for reduced order modeling of fluid dynamics. Designing accurate and efficient reduced order models for nonlinear fluid dynamic problems is challenging for many practical engineering applications. Classical projection-based model reduction methods generate reduced systems by projecting full-order differential operators into low-dimensional subspaces. However, these techniques usually lead to severe instabilities in the presence of highly nonlinear dynamics, which dramatically deteriorates the accuracy of the reduced-order models. In contrast, our new framework exploits linear multistep networks, based on implicit Adams–Moulton schemes, to construct the reduced system. The advantage is that the method optimally approximates the full order model in the low-dimensional space with a given supervised learning task. Moreover, our approach is non-intrusive, such that it can be applied to other complex nonlinear dynamical systems with sophisticated legacy codes. We demonstrate the performance of our method through the numerical simulation of a two-dimensional flow past a circular cylinder with Reynolds number Re = 100. The results reveal that the new data-driven model is significantly more accurate than standard projection-based approaches.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1963

Some Effects of Sweep and Aspect Ratio on the Transonic Flutter Characteristics of a Series of Thin Cantilever Wings Having a Taper Ratio of 0.6

An investigation of the flutter characteristics of a series of thin cantilever wings having taper ratios of 0.6 was conducted in the Langley transonic blowdown tunnel at Mach numbers between 0.76 and 1.42. The angle of sweepback was varied from 0 degrees to 60 degrees on wings of aspect ratio 4, and the aspect ratio was varied from 2.4 to 6.4 on wings with 45 degrees of sweepback. The results are presented as ratios between the experimental flutter speeds and the reference flutter speeds calculated on the basis of incompressible two-dimensional flow. These ratios, designated the flutter-speed ratios, are given as functions of Mach number for the various wings. The flutter-speed ratios were characterized, in most cases, by values near 1.0 at subsonic speeds with large increases in the speed ratios in the range of supersonic speeds investigated. Increasing the sweep effected increases in the flutter-speed ratios between 0 degrees and 30 degrees followed by progressive reductions of the speed ratios to nearly 1.0 as the sweep was increased from 30 degrees to 60 degrees. Reducing the aspect ratio from 6.4 to 2.4 resulted in progressively larger values of the flutter-speed ratios throughout the Mach number range investigated.

CANTILEVER WING↗
NASA NTRS1971

Time dependent shear stress and temperature distribution over an insulated flat plate moving at hypersonic speed.

The laminar two-dimensional flow over a stepwise accelerated flat plate moving with hypersonic speed at zero angle of attack is analysed. The governing equations in the self-similar form are linearized and solved numerically for small times. The solutions obtained are the deviations of the velocity and the temperature profiles from those of steady state. The presented results may be used to find the first order boundary layer induced pressure on the plate.

Rodkiewicz, C. M.↗
NASA NTRS1974

Two-demensional analysis of heat and mass transfer in porous media using the strongly implicit procedure

Numerical results of the heat and mass transfer in a porous matrix are presented. The coupled, nonlinear partial differential equations describing this physical phenomenon are solved in finite difference form for two dimensions, using a new iterative technique (the strongly implicit procedure). The influence of the external environment conditions (heating and pressure) is shown to produce two-dimensional flow in the porous matrix. Typical fluid and solid temperature distributions in the porous matrix and internal pressure distributions are presented.

Curry, D. M.↗