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 1,063 records · Page 59

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

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗

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

Evidence of a large-scale gradient in the solar wind velocity

Measurements of the solar wind velocity are compared at two widely separated locations, using plasma data obtained by Mariner 5, en route to Venus, and with the near-earth satellites Explorer 33, 34, and 35. A previous study of the propagation of interplanetary sector boundaries between Mariner and the earth had implied the existence of large scale velocity gradient which was interpreted as a latitude gradient of approximately 13 km/sec per degree of latitude. The results of the present investigation in which the earlier results were extended to the overall solar wind, without regard to the presence of sector boundaries, suggests that a latitude dependence of the solar wind velocity is the most plausible interpretation of the large-scale velocity gradient.

Rhodes, E. J., Jr.↗

Lectures series in computational fluid dynamics

The lecture notes cover the basic principles of computational fluid dynamics (CFD). They are oriented more toward practical applications than theory, and are intended to serve as a unified source for basic material in the CFD field as well as an introduction to more specialized topics in artificial viscosity and boundary conditions. Each chapter in the test is associated with a videotaped lecture. The basic properties of conservation laws, wave equations, and shock waves are described. The duality of the conservation law and wave representations is investigated, and shock waves are examined in some detail. Finite difference techniques are introduced for the solution of wave equations and conservation laws. Stability analysis for finite difference approximations are presented. A consistent description of artificial viscosity methods are provided. Finally, the problem of nonreflecting boundary conditions are treated.

Thompson, Kevin W.↗

Entropy and climate. I - ERBE observations of the entropy production of the earth

An approximate method for estimating the global distributions of the entropy fluxes flowing through the upper boundary of the climate system is introduced, and an estimate of the entropy exchange between the earth and space and the entropy production of the planet is provided. Entropy fluxes calculated from the Earth Radiation Budget Experiment measurements show how the long-wave entropy flux densities dominate the total entropy fluxes at all latitudes compared with the entropy flux densities associated with reflected sunlight, although the short-wave flux densities are important in the context of clear sky-cloudy sky net entropy flux differences. It is suggested that the entropy production of the planet is both constant for the 36 months of data considered and very near its maximum possible value. The mean value of this production is 0.68 x 10 exp 15 W/K, and the amplitude of the annual cycle is approximately 1 to 2 percent of this value.

Stephens, G. L.↗