Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate boundary”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1954

Some Observations of Shock-induced Turbulent Separation on Supersonic Diffusers

A survey of experimental data at supersonic speed indicated that shock-induced separation of a turbulent boundary layer will result for Mach numbers of approximately 1.33 or greater when a theoretical stream static-pressure-rise ratio of approximately 1.89 occurs across a shock interacting with the boundary layer. The significance of this tentative criterion for turbulent boundary-layer separation is discussed with respect to the design of supersonic diffusers.

BOUNDARY-LAYER CONTROL--INTERNAL AERODYNAMICS↗
NASA NTRS1977

High coronal structure of high velocity solar wind stream sources

It is shown analytically that the transition from a high-speed stream source to the ambient coronal conditions is quite rapid in longitude in the high corona. This sharp eastern coronal boundary for the solar wind stream sources is strongly suggested by the solar wind 'dwells' which appear in plots of solar wind velocity against constant-radial-velocity-approximation source longitudes. The possibility of a systematic velocity-dependent effect in the constant-radial-velocity approximation, which would cause this boundary to appear sharper than it is, is investigated. A velocity-dependent interplanetary propagation effect or a velocity-dependent 'source altitude' are two possible sources of such a systematic effect. It is shown that, for at least some dwells, significant interplanetary effects are not likely. The variation of the Alfvenic critical radius in solar wind dwells is calculated, showing that the high-velocity stream originates from a significantly lower altitude than the ambient solar wind.

Nolte, J. T.↗
NASA NTRS1980

Spectral methods for problems in complex geometries

The properties of spectral methods are surveyed and their extension to solve problems in complex geometries is developed. A new iteration procedure is introduced to solve efficiently the full matrix equations resulting from spectral approximations to nonconstant coefficient boundary-value problems in complex geometries. It is shown that the work required to solve these spectral equations exceeds that of solving the lowest-order finite-difference approximation to the same problem by only O(N log N).

Orszag, S. A.↗
NASA NTRS2014

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗
NASA NTRS1994

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

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.

Hesthaven, J. S.↗
NASA NTRS2017

Stereo Particle Image Velocimetry Measurements of Transition Downstream of a Forward-Facing Step in a Swept-Wing Boundary Layer

Stereo particle image velocimetry measurements were performed downstream of a forward-facing step in a stationary-crossflow dominated flow. Three different step heights were studied with the same leading-edge roughness configuration to determine the effect of the step on the evolution of the stationary-crossflow. Above the critical step height, which is approximately 68% of the boundary-layer thickness at the step, the step caused a significant increase in the growth of the stationary crossflow. For the largest step height studied (68%), premature transition occurred shortly downstream of the step. The stationary crossflow amplitude only reached approximately 7% of U(sub e) in this case, which suggests that transition does not occur via the high-frequency secondary instabilities typically associated with stationary crossflow transition. The next largest step of 60% delta still caused a significant impact on the growth of the stationary crossflow downstream of the step, but the amplitude eventually returned to that of the baseline case, and the transition front remained the same. The smallest step height (56%) only caused a small increase in the stationary crossflow amplitude and no change in the transition front. A final case was studied in which the roughness on the leading edge of the model was enhanced for the lowest step height case to determine the impact of the stationary crossflow amplitude on transition. The stationary crossflow amplitude was increased by approximately four times, which resulted in premature transition for this step height. However, some notable differences were observed in the behavior of the stationary crossflow mode, which indicate that the interaction mechanism which results in the increased growth of the stationary crossflow downstream of the step may be different in this case compared to the larger step heights.

Eppink, Jenna L.↗
NASA NTRS1993

Kinematics of the New Zealand plate boundary: Relative motion by GPS across networks of 1000 km and 50 km spacing

The NASA/DOSE 'Kinematics of the New Zealand Plate Boundary' experiment is a four-year cooperative Global Positioning System (GPS) experiment involving 6 universities and institutions in New Zealand and the United States. The investigation covers two scales, the first on the scale of plates (approximately 1000 km) and the second is on the scale of the plate boundary zone (approximately 50 km). In the first portion of the experiment, phase A, the objective is to make direct measurements of tectonic plate motion between the Australian and Pacific plates using GPS in order to determine the Euler vector of this plate pair. The phase A portion of this experiment was initiated in December 1992 with the first-epoch baseline measurements on the large scale network. The network will be resurveyed two years later to obtain velocities. The stations which were observed for phase A are shown and listed. Additional regional stations which will be used for this study are listed and are part of either CIGNET or other global tracking networks. The phase A portion of the experiment is primarily the responsibility of the UNAVCO investigators. Therefore, this report concentrates on phase A. The first year of NASA funding for phase A included only support for the field work. Processing and analysis will take place with the second year of funding. The second part of the experiemnt measured relative motion between the Australian and Pacific plates across the pate boundary zone between Hokitika and Christchurch on the South Island of New Zealand. The extent and rate of deformation will be determined by comparisons with historical, conventional surveys and by repeated GPS measurements to be made in two years. This activity was the emphasis of the LDGO portion of the study. An ancillary experiment, phase C, concentrated on plate boundary deformation in the vicinity of Wellington and was done as part of training during the early portion of the field campaign. Details of the objectives of the field investigations are given in the appendix. An overview of the 1992 GPS field program is also given in the appendix.

Meertens, Charles M.↗
NASA NTRS1959

Space-Time Correlations and Spectra of Wall Pressure in a Turbulent Boundary Layer

Measurements of the statistical properties of the fluctuating wall pressure produced by a subsonic turbulent boundary layer are described. The measurements provide additional information about the structure of the turbulent boundary layer; they are applicable to the problems of boundary-layer induced noise inside an airplane fuselage and to the generation of waves-on water. The spectrum of the wall pressure is presented in dimensionless form. The ratio of the root-mean-square wall pressure to the free-stream dynamic pressure is found to be a constant square root of bar P(sup 2)/q(sub infinity) = 0.006 independent of Mach number and Reynolds number. In addition, space- time correlation measurements in the stream direction show that pressure fluctuations whose scale is greater than or equal to 0.3 times the boundary-layer thickness are convected with the convection speed U(sub c) = 0.82U(sub infinity) where U(infinity) is the free-stream velocity and have lost their identity in a distance approximately equal to 10 boundary-layer thicknesses.

Willmarth, W. W.↗
NASA NTRS1958

The Turbulent Boundary Layer on a Rough Curvilinear Surface

A number of semiempirical approximate methods exist for determining the characteristics of the turbulent boundary layer on a curvilinear surface. At present, among these methods, the one proposed by L. G. Loitsianskii is given frequent practical application. This method is sufficiently effective and permits, in the case of wing profiles with technically smooth surfaces, calculating the basic characteristics of the boundary layer and the values of the overall drag with an accuracy which suffices for practical purposes. The idea of making use of the basic integral momentum equation ((d delta(sup xx))/dx) + ((V' delta(sup xx))/V) (2 + H) = (tau(sub 0))/(rho V(exp 2)) proves to be fruitful also for the solution of the problems in the determination of the characteristics of the turbulent boundary layer on a rough surface.

Droblenkov, V. F.↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

A filament of energetic particles near the high-latitude dawn magnetopause

The Geotail satelite detected a filament of tailward-streaming energetic particles spatially separated from the boundary layer of energetic particles at the high-latitude dawn magnetopause at a downstream distance of approximately 80 R(sub E) on October 27, 1992. During this event, the composition and charge states of energetic ions at energies above approximately 10 keV show significant intermix of ions from solar wind and ionospheric sources. Detailed analysis leads to the deduction that the filament was moving southward towards the neutral sheet at an average speed of approximately 80 km/s, implying an average duskward electric field of approximately 1 mV/m. Its north-south dimension was approximately 1 R(sub E) and it was associated with an earthward directed field-aligned current of approximately 5 mA/m. The filament was separated from the energetic particle boundary layer straddling the magnetopause by approximately 0.8 R(sub E) and was inferred to be detached from the boundary layer at downstream distance beyond approximately 70 R(sub E) in the distant tail.

Lui, A. T. Y.↗
DOE OSTI2023

Quantifying unitary flow efficiency and entanglement for many-body localization

Holography tells us that there should be a mapping between a boundary Hamiltonian and bulk degrees of freedom. We probe the bulk geometry of the Wegner Wilson Flow (WWF) in the context of many-body localization by addressing efficiency and bulk entanglement growth measures through approximating upper bounds on the boundary entanglement entropy. Here, we connect these upper bounds to the Fubini-Study metric and clarify how a central quantity, the information fluctuation complexity, distinguishes bulk unitary rotation from entanglement production. We also give a short new proof of the small incremental entangling theorem in the absence of ancillas, achieving a dimension independent, universal factor of c = 2.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1952

Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System

The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.

Mager, Arthur↗
NASA NTRS2003

In-Flight Capability for Evaluating Skin-Friction Gages and Other Near-Wall Flow Sensors

An 8-in.-square boundary-layer sensor panel has been developed for in-flight evaluation of skin-friction gages and other near-wall flow sensors on the NASA Dryden Flight Research Center F-15B/Flight Test Fixture (FTF). Instrumentation on the sensor panel includes a boundary-layer rake, temperature sensors, static pressure taps, and a Preston tube. Space is also available for skin-friction gages or other near-wall flow sensors. Pretest analysis of previous F-15B/FTF flight data has identified flight conditions suitable for evaluating skin-friction gages. At subsonic Mach numbers, the boundary layer over the sensor panel closely approximates the two-dimensional (2D), law-of-the-wall turbulent boundary layer, and skin-friction estimates from the Preston tube and the rake (using the Clauser plot method) can be used to evaluate skin-friction gages. At supersonic Mach numbers, the boundary layer over the sensor panel becomes complex, and other means of measuring skin friction are needed to evaluate the accuracy of new skin-friction gages. Results from the flight test of a new rubber-damped skin-friction gage confirm that at subsonic Mach numbers, nearly 2D, law-of-the-wall turbulent boundary layers exist over the sensor panel. Sensor panel data also show that this new skin-friction gage prototype does not work in flight.

Bui, Trong T.↗
NASA NTRS1986

Estimation of coefficients and boundary parameters in hyperbolic systems

Semi-discrete Galerkin approximation schemes in connection with inverse problems for the estimation of spatially varying coefficients and boundary condition parameters in second order hyperbolic systems typical of those arising in 1-D surface seismic problems are considered. Spline based algorithms are proposed for which theoretical convergence results along with a representative sample of numerical findings are given.

Banks, H. T.↗
NASA NTRS1987

A general form for the dissipation length scale in turbulent shear flows

It has been found that, for a wide range of turbulent wall-bounded shear flows with mean velocity profile U(y), the scale length L-epsilon determining the dissipation is approximately described in terms of distance from the wall (y), the mean shear (dU/dy), and the variance of the normal component of turbulence (bar u-squared), by the formula: L-epsilon to the -1 power is approximately equal to (A sub B / y) + A sub S (dU/dy / square root of bar u-squared) where L-epsilon = epsilon/(bar u-squared) to the 3/2 power. To match with shear-free boundary layers, A sub B is approximately 0.27, and to match with the log layer, A sub S is approximately 0.46. The shear flows tested here were: boundary layers over a flat plate, sink flow, oscillatory flow, and channel flow. The use of the square root of bar u-squared as a velocity scale minimizes the effects of Reynolds number. However, the formula fails within a distance of order L-epsilon for the regions where dU/dy = 0.

Hunt, J. C. R.↗
NASA NTRS1972

Calculation of heat transfer on shuttle type configurations including the effects of variable entropy at boundary layer edge

A relatively simple method is presented for including the effect of variable entropy at the boundary-layer edge in a heat transfer method developed previously. For each inviscid surface streamline an approximate shockwave shape is calculated using a modified form of Maslen's method for inviscid axisymmetric flows. The entropy for the streamline at the edge of the boundary layer is determined by equating the mass flux through the shock wave to that inside the boundary layer. Approximations used in this technique allow the heating rates along each inviscid surface streamline to be calculated independent of the other streamlines. The shock standoff distances computed by the present method are found to compare well with those computed by Maslen's asymmetric method. Heating rates are presented for blunted circular and elliptical cones and a typical space shuttle orbiter at angles of attack. Variable entropy effects are found to increase heating rates downstream of the nose significantly higher than those computed using normal-shock entropy, and turbulent heating rates increased more than laminar rates. Effects of Reynolds number and angles of attack are also shown.

Dejarnette, F. R.↗