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 937 records · Page 52

Transition on a sharp cone at high enthalpy - New measurements in the shock tunnel T5 at GALCIT

An exploratory study of the laminar, transitional and turbulent boundary layer on a slender cone in hypervelocity flow was conducted in the high-enthalpy shock tunnel T5. A novel flow visualization technique using sodium seeding to increase the sensitivity of conventional interferometric techniques by resonant enhancement of the refractivity of the medium was developed to study its structure. The experiments were designed to cover a large range of specific reservoir enthalpy, ranging from the perfect-gas region to the range where significant oxygen and some nitrogen dissociation and recombination effects may be expected in the boundary layer. Surface heat flux measurements indicate that the transition Reynolds number is approximately the same as that observed in cold hypersonic flow in conventional continuous flow facilities, but significantly lower than that observed in 'quiet' tunnels. No significant effect of reservoir specific enthalpy on transition was observed. The heat flux distribution showed differences between the air and the nitrogen experiments in the laminar regime. These are consistent with the presence of recombination of oxygen at the wall. This view is supported by the fact that this effect is absent in the low-enthalpy region, where the heat flux distribution agrees with the similarity theory.

Germain, P.↗

A preliminary investigation of boundary-layer transition along a flat plate with adverse pressure gradient

Boundary-layer surveys were made throughout the transition region along a smooth flat plate placed in an airstream of practically zero turbulence and with an adverse pressure gradient. The boundary-layer Reynolds number at the laminar separation point was varied from 1,800 to 2,600. The test data, when considered in the light of certain theoretical deductions, indicated that transition probably began with separation of the laminar boundary layer. The extent of the transition region, defined as the distance from a calculated laminar separation point to the position of the first fully developed turbulent boundary-layer profile, could be expressed as a constant Reynolds number run of approximately 70,000. Some speculations are presented concerning the application of the foregoing concepts, after certain assumptions have been made, to the problem of the connection between transition on the upper surface of an airfoil at high angles of attack and the maximum lift.

Von Doenhoff, Albert E↗

Proton aurora and substorm intensifications

Ground based measurements from the CANOPUS array of meridian scanning photometers and precipitating ion and electron data from the DMSP F9 satellite show that the electron arc which brightens to initiate substorm intensifications is formed within a region of intense proton precipitation that is well equatorward (approximately four to six degrees) of the nightside open-closed field line boundary. The precipitating protons are from a population that is energized via earthward convection from the magnetotail into the dipolar region of the magnetosphere and may play an important role in the formation of the electron arcs leading to substorm intensifications on dipole-like field lines.

Samson, J. C.↗

Aerodynamic parameter estimation via Fourier modulating function techniques

Parameter estimation algorithms are developed in the frequency domain for systems modeled by input/output ordinary differential equations. The approach is based on Shinbrot's method of moment functionals utilizing Fourier based modulating functions. Assuming white measurement noises for linear multivariable system models, an adaptive weighted least squares algorithm is developed which approximates a maximum likelihood estimate and cannot be biased by unknown initial or boundary conditions in the data owing to a special property attending Shinbrot-type modulating functions. Application is made to perturbation equation modeling of the longitudinal and lateral dynamics of a high performance aircraft using flight-test data. Comparative studies are included which demonstrate potential advantages of the algorithm relative to some well established techniques for parameter identification. Deterministic least squares extensions of the approach are made to the frequency transfer function identification problem for linear systems and to the parameter identification problem for a class of nonlinear-time-varying differential system models.

Pearson, A. E.↗

Blended Wing Body Systems Studies: Boundary Layer Ingestion Inlets With Active Flow Control

A CFD analysis was performed on a Blended Wing Body (BWB) aircraft with advanced, turbofan engines analyzing various inlet configurations atop the aft end of the aircraft. The results are presented showing that the optimal design for best aircraft fuel efficiency would be a configuration with a partially buried engine, short offset diffuser using active flow control, and a D-shaped inlet duct that partially ingests the boundary layer air in flight. The CFD models showed that if active flow control technology can be satisfactorily developed, it might be able to control the inlet flow distortion to the engine fan face and reduce the powerplant performance losses to an acceptable level. The weight and surface area drag benefits of a partially submerged engine shows that it might offset the penalties of ingesting the low energy boundary layer air. The combined airplane performance of such a design might deliver approximately 5.5% better aircraft fuel efficiency over a conventionally designed, pod-mounted engine.

Geiselhart, Karl A.↗

Evaluation of Arctic Sea Ice Thickness Simulated by Arctic Ocean Model Intercomparison Project Models

Six Arctic Ocean Model Intercomparison Project model simulations are compared with estimates of sea ice thickness derived from pan-Arctic satellite freeboard measurements (2004-2008); airborne electromagnetic measurements (2001-2009); ice draft data from moored instruments in Fram Strait, the Greenland Sea, and the Beaufort Sea (1992-2008) and from submarines (1975-2000); and drill hole data from the Arctic basin, Laptev, and East Siberian marginal seas (1982-1986) and coastal stations (1998-2009). Despite an assessment of six models that differ in numerical methods, resolution, domain, forcing, and boundary conditions, the models generally overestimate the thickness of measured ice thinner than approximately 2 mand underestimate the thickness of ice measured thicker than about approximately 2m. In the regions of flat immobile landfast ice (shallow Siberian Seas with depths less than 25-30 m), the models generally overestimate both the total observed sea ice thickness and rates of September and October ice growth from observations by more than 4 times and more than one standard deviation, respectively. The models do not reproduce conditions of fast ice formation and growth. Instead, the modeled fast ice is replaced with pack ice which drifts, generating ridges of increasing ice thickness, in addition to thermodynamic ice growth. Considering all observational data sets, the better correlations and smaller differences from observations are from the Estimating the Circulation and Climate of the Ocean, Phase II and Pan-Arctic Ice Ocean Modeling and Assimilation System models.

sea ice↗

Aerodynamic Characterization of a Thin, High-Performance Airfoil for Use in Ground Fluids Testing

The FAA has worked with Transport Canada and others to develop allowance times for aircraft operating in ice-pellet precipitation. Wind-tunnel testing has been carried out to better understand the flowoff characteristics and resulting aerodynamic effects of anti-icing fluids contaminated with ice pellets using a thin, high-performance wing section at the National Research Council of Canada Propulsion and Icing Wind Tunnel. The objective of this paper is to characterize the aerodynamic behavior of this wing section in order to better understand the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination. Aerodynamic performance data, boundary-layer surveys and flow visualization were conducted at a Reynolds number of approximately 6.0×10(exp 6) and a Mach number of 0.12. The clean, baseline model exhibited leading-edge stall characteristics including a leading-edge laminar separation bubble and minimal or no separation on the trailing edge of the main element or flap. These results were consistent with expected 2-D aerodynamics and showed no anomalies that could adversely affect the evaluation of anti-icing fluids and ice-pellet contamination on the wing. Tests conducted with roughness and leading-edge flow disturbances helped to explain the aerodynamic impact of the anti-icing fluids and contamination. The stalling characteristics of the wing section with fluid and contamination appear to be driven at least partially by the effects of a secondary wave of fluid that forms near the leading edge as the wing is rotated in the simulated takeoff profile. These results have provided a much more complete understanding of the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination on this wing section. This is important since these results are used, in part, to develop the ice-pellet allowance times that are applicable to many different airplanes.

Broeren, Andy P.↗

Aerodynamic Characterization of a Thin, High-Performance Airfoil for Use in Ground Fluids Testing

The FAA has worked with Transport Canada and others to develop allowance times for aircraft operating in ice-pellet precipitation. Wind-tunnel testing has been carried out to better understand the flowoff characteristics and resulting aerodynamic effects of anti-icing fluids contaminated with ice pellets using a thin, high-performance wing section at the National Research Council of Canada Propulsion and Icing Wind Tunnel. The objective of this paper is to characterize the aerodynamic behavior of this wing section in order to better understand the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination. Aerodynamic performance data, boundary-layer surveys and flow visualization were conducted at a Reynolds number of approximately 6.0 x 10(exp 6) and a Mach number of 0.12. The clean, baseline model exhibited leading-edge stall characteristics including a leading-edge laminar separation bubble and minimal or no separation on the trailing edge of the main element or flap. These results were consistent with expected 2-D aerodynamics and showed no anomalies that could adversely affect the evaluation of anti-icing fluids and ice-pellet contamination on the wing. Tests conducted with roughness and leading-edge flow disturbances helped to explain the aerodynamic impact of the anti-icing fluids and contamination. The stalling characteristics of the wing section with fluid and contamination appear to be driven at least partially by the effects of a secondary wave of fluid that forms near the leading edge as the wing is rotated in the simulated takeoff profile. These results have provided a much more complete understanding of the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination on this wing section. This is important since these results are used, in part, to develop the ice-pellet allowance times that are applicable to many different airplanes.

Broeren, Andy P.↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low 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 pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jet flow control used to suppress boundary layer separation and increase the maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate streamwise direction. It also generates turbulent eddies at the oscillation frequency, which are typically large relative to the scales of boundary layer turbulence, and which augment mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from a nozzle downstream of the oscillator, and an array of sweep jets which suppresses boundary layer separation are performed. Simulation results are compared to data from a dedicated validation experiment of a single oscillator and its sweep jet, and from a wind tunnel test of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets. A critical step in the work is the development of realistic time-dependent sweep jet inflow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the computational fluid dynamics (CFD) solver Overow, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

The SMM Model as a Boundary Value Problem Using the Discrete Diffusion Equation

A generalized single step stepwise mutation model (SMM) is developed that takes into account an arbitrary initial state to a certain partial difference equation. This is solved in both the approximate continuum limit and the more exact discrete form. A time evolution model is developed for Y DNA or mtDNA that takes into account the reflective boundary modeling minimum microsatellite length and the original difference equation. A comparison is made between the more widely known continuum Gaussian model and a discrete model, which is based on modified Bessel functions of the first kind. A correction is made to the SMM model for the probability that two individuals are related that takes into account a reflecting boundary modeling minimum microsatellite length. This method is generalized to take into account the general n-step model and exact solutions are found. A new model is proposed for the step distribution.

Campbell, Joel↗

Accurate computation and continuation of homoclinic and heteroclinic orbits for singular perturbation problems

In earlier papers, Doedel and the authors have developed a numerical method and derived error estimates for the computation of branches of heteroclinic orbits for a system of autonomous ordinary differential equations in R(exp n). The idea of the method is to reduce a boundary value problem on the real line to a boundary value problem on a finite interval by using a local (linear or higher order) approximation of the stable and unstable manifolds. A practical limitation for the computation of homoclinic and heteroclinic orbits has been the difficulty in obtaining starting orbits. Typically these were obtained from a closed form solution or via a homotopy from a known solution. Here we consider extensions of our algorithm which allow us to obtain starting orbits on the continuation branch in a more systematic way as well as make the continuation algorithm more flexible. In applications, we use the continuation software package AUTO in combination with some initial value software. The examples considered include computation of homoclinic orbits in a singular perturbation problem and in a turbulent fluid boundary layer in the wall region problem.

Vaughan, William W.↗

Turbulent heat transfer prediction method for application to scramjet engines

An integral method for predicting boundary layer development in turbulent flow regions on two-dimensional or axisymmetric bodies was developed. The method has the capability of approximating nonequilibrium velocity profiles as well as the local surface friction in the presence of a pressure gradient. An approach was developed for the problem of predicting the heat transfer in a turbulent boundary layer in the presence of a high pressure gradient. The solution was derived with particular emphasis on its applicability to supersonic combustion; thus, the effects of real gas flows were included. The resulting integrodifferential boundary layer method permits the estimation of cooling reguirements for scramjet engines. Theoretical heat transfer results are compared with experimental combustor and noncombustor heat transfer data. The heat transfer method was used in the development of engine design concepts which will produce an engine with reduced cooling requirements. The Langley scramjet engine module was designed by utilizing these design concepts and this engine design is discussed along with its corresponding cooling requirements. The heat transfer method was also used to develop a combustor cooling correlation for a combustor whose local properties are computed one dimensionally by assuming a linear area variation and a given heat release schedule.

Pinckney, S. Z.↗

Heat-Transfer Measurements on a 5.5- Inch-Diameter Hemispherical Concave Nose in Free Flight at Mach Numbers up to 6.6

The aerodynamic heat transfer to a hemispherical concave nose has been measured in free flight at Mach numbers from 3.5 to 6.6 with corresponding Reynolds numbers based on nose diameter from 7.4 x 10(exp 6) to 14 x 10(exp 6). Over the test Mach number range the heating on the cup nose, expressed as a ratio to the theoretical stagnation-point heating on a hemisphere nose of the same diameter, varied from 0.05 to 0.13 at the stagnation point of the cup, was approximately 0.1 at other locations within 40 deg of the stagnation point, and varied from 0.6 to 0.8 just inside the lip where the highest heating rates occurred. At a Mach number of 5 the total heat input integrated over the surface of the cup nose including the lip was 0.55 times the theoretical value for a hemisphere nose with laminar boundary layer and 0.76 times that for a flat face. The heating at the stagnation point was approximately 1/5 as great as steady-flow tunnel results. Extremely high heating rates at the stagnation point (on the order of 30 times the stagnation-point values of the present test), which have occurred in conjunction with unsteady oscillatory flow around cup noses in wind-tunnel tests at Mach and Reynolds numbers within the present test range, were not observed.

Levine, Jack↗

A study of the effects of Lebu devices on turbulent boundary layer drag

Initial measurements of the changes in local skin friction, velocity profile shape, and turbulence structure which result from the placement of tandem plates parallel to the wall in the outer region of thick turbulent boundary layers were made. Using a tunnel with a .75 m x 1.2 m x 7.3 m test section, which diverged so as to keep the pressure gradient less than 2x1000/ft, on the test wall, a skin friction reduction of approximately 30% was measured at xi/h = 62. This relaxed to a reduction of approximately 16% at xi/h = 124 for h/delta M = .6. The c sub f measurements for both the normal and modified boundary layers were obtained by measuring the slope of the velocity profile within the linear sublayer. Visual results indicated a continued presence of strong large eddy structure downstream of the devices. Local skin friction reduction of 12% at xi/h = 62 was also obtained with the manipulators above the boundary layer at y/delta m = 1.1.

Falco, R. E.↗

Accuracy of approximations to the Navier-Stokes equations

The systems of truncated differential equations that have been proposed to reduce the complexity and large computational costs of solutions to the full Navier-Stokes equations are considered. These systems are computationally efficient and capture all the physically relevant behavior. The systems follow a certain hierarchy: (1) the classical boundary-layer equations with specified edge properties (usually the streamwise pressure distribution); (2) the coupled boundary-layer/inviscid equations; (3) the so-called thin-layer equations that discard streamwise diffusion; and (4) the Navier-Stokes equations. Consideration is given to each of these approximations applied to an incompressible, laminar-separating flow at low and moderate Reynolds numbers. It is pointed out that for any flow or region of flow for which viscous-inviscid interaction effects are small, classical boundary-layer equations will provide a satisfactory description of the viscous flow at a fraction of the computational cost of any higher approximations.

Murphy, J. D.↗

Delineation of the boundaries of a buried pre-glacial valley with LANDSAT-1 data

The continuity of a narrow meandering strip of Udoll (prairie) soils running east and west for approximately 40 miles across north central Indiana in an area predominantly of Udalfs (timber soils) was detected from LANDSAT-1 data taken on June 9, 1973. This data was processed through a clustering procedure and classified with resulting increased definition of the boundaries among soils grouped according to nine categories and vegetation to two categories of reflectance. This dark stretch of prairie soil is believed to have formed in the heavy textured, poorly drained glacial debris which filled a major pre-glacial tributary of the Teays River System. Ready identification and location of the valley has significance to soil survey and land classification people as a guide to soil classification and land use and to geologists as a guide to location of a potentially economically significant aquifer.

Peterson, J. B.↗

Warm core ring dynamics derived from satellite imagery

A reinterpretation of the life history of warm core ring 82-b is constructed from drifter trajectories, hydrographic profiles, and satellite-derived sea surface temperature fields. An analysis of the observations strongly suggests this ring is fundamentally a dipole structure for much of its existence. The dipole has unequal strength vortices, with the cyclone being the weaker component. Quantification of the dipole indicates the anticyclonic is oval shaped and rotates at approximately 7.7 deg per day, while the cyclone revolves around the anticyclone at about 14.5 deg per day, yielding an orbit period of approximately 25 days for the cyclone. The elusive part of the dipole explanation is in interpreting its surface signature as a function of time, as the cyclone appears to disappear when in proximity to the steep continental slope. Vortex modeling suggests that when the weaker cyclonic eddy is placed close to a boundary, the cyclone is drawn into a filament as it is advected through the gap between the anticyclone and the boundary. This means the cyclone is either in a filamentation state or undergoing the straining associated with filamentation approximately three fourths of this time, which is in complete agreement with the satellite imagery. In addition, the 82-B dipole configuration is sufficiently robust to survive continuous interaction with the continental slope and the onset of a Gulf Stream interaction. Although only warm core ring 82-B is analyzed in detail (and found to be a vortex pair), dipole configurations are found in the warm core ring 81-F time series as well.

Hooker, Stanford B.↗