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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
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.↗
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.↗
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.↗
NASA NTRS1972

Decision surface estimate of nonlinear system stability domain by Lie series method.

The Lie series recursive algorithm for Zubov's partial differential equation is used to generate two sets of points, where one represents the exact asymptotic stability boundary of an equilibrium state of the nonlinear system under consideration and the other is interior to it. Based on these two sets of data as training samples of two classes, a decision hypersurface can be determined such that it is a close approximation of the asymptotic stability boundary.

Kormanik, J.↗
NASA NTRS2018

Quantifying Shortcomings in the 2-D Pc Calculation

The probability of collision, Pc, between two Earth-orbiting satellites can often, but not always, be approximated adequately using the "2D-Pc" formulation. The objective is to find a set of "boundary conditions" that ensure the 2D-Pc approximation be sufficiently accurate, so that it may be determined when computationally-intensive Brute Force Monte Carlo1 (BFMC) Pc estimates are required. We critically examine the assumptions used in the formulation of the 2D-Pc approximation and then formulate 2D-Pc boundary condition tests to check if these assumptions are satisfied adequately.

Analysis↗
NASA NTRS1989

Development of an integrated BEM approach for hot fluid structure interaction

The progress made toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-Orbit engine hot section components is reported. The convective viscous integral formulation was derived and implemented in the general purpose computer program GP-BEST. The new convective kernel functions, in turn, necessitated the development of refined integration techniques. As a result, however, since the physics of the problem is embedded in these kernels, boundary element solutions can now be obtained at very high Reynolds number. Flow around obstacles can be solved approximately with an efficient linearized boundary-only analysis or, more exactly, by including all of the nonlinearities present in the neighborhood of the obstacle. The other major accomplishment was the development of a comprehensive fluid-structure interaction capability within GP-BEST. This new facility is implemented in a completely general manner, so that quite arbitrary geometry, material properties and boundary conditions may be specified. Thus, a single analysis code (GP-BEST) can be used to run structures-only problems, fluids-only problems, or the combined fluid-structure problem. In all three cases, steady or transient conditions can be selected, with or without thermal effects. Nonlinear analyses can be solved via direct iteration or by employing a modified Newton-Raphson approach.

Dargush, G. F.↗
NASA NTRS1990

Acoustic waves superimposed on incompressible flows

The use of incompressible approximations in deriving solutions to the Lighthill wave equation was investigated for problems where an analytical solution could be found. A particular model problem involves the determination of the sound field of a spherical oscillating bubble in an ideal fluid. It is found that use of incompressible boundary conditions leads to good approximations in the important region of high acoustic wave number.

Hodge, Steve↗
NASA NTRS1978

Magnetic fields and flows between 1 AU and 0.3 AU during the primary mission of HELIOS 1

The recurrent flow and field patterns observed by HELIOS 1, and the relation between these patterns and coronal holes are discussed. Four types of recurrent patterns were observed: a large recurrent stream, a recurrent slow (quiet) flow, a rapidly evolving flow, and a recurrent compound stream. There recurrent streams were not stationary, for although the sources recurred at approximately the same longitudes on successive rotations, the shapes and latitudinal patterns changed from one rotation to the next. A type of magnetic field and plasma structure characterized by a low ion temperature and a high magnetic field intensity is described as well as the structures of stream boundaries between the sun at approximately 0.3 AU.

Burlaga, L. F.↗