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,009 records · Page 56

The lunar conductivity profile and the nonuniqueness of electromagnetic data inversion.

A review of the theory for the electromagnetic functional used to date to determine the lunar conductivity profile from spectral analyses of lunar magnetometer data. The 'hard' boundary condition used by Sonett et al. (1971) and others appears to be a good approximation for the sunlit lunar hemisphere. The use of the spectral data in conjunction with the functional to find a least squares conductivity profile is examined from the point of view of the nonuniqueness of nonlinear estimation. An unsuccessful search was carried out to find models that were both geophysically reasonable and consistent with the data. It is asserted that the failure to find these models is due to the inadequacy of the functional. Of the models generated, those that best fit the data all have a conductivity peak at a depth of 240 km similar to that of the model discussed by Sonett et al. However, all of these models are quite distinct elsewhere in the profile. It is shown how such models are dependent on both the nature of the functional as well as the initial guess in the least squares procedure. A correlation analysis shows that the only resolvable features from this type of model are a low conductivity crust, a high conductivity peak of limited radial extent, and a low conductivity zone beneath the peak.

Phillips, R. J.↗

Synthesis of multiple shaped beam antenna patterns

Results are presented of research into the problem of finding an excitation of a given antenna such that the desired radiation pattern is approximated to within acceptable limits. This is to be done in such a fashion that boundary conditions involving hardware limitations may be inserted into the problem. The intended application is synthesis of multiple shaped beam antennas. Since this is perhaps the most difficult synthesis problem an antenna engineer is likely to encounter, the approach taken was to include as a by-product capability for synthesizing simpler patterns. The synthesis technique has been almost totally computerized. The class of antennas which may be synthesized with the computer program are those which may be represented as planar (continuous or discrete) current distributions. The technique is not limited in this sense and could indeed by extended to include, for example, the synthesis of conformal arrays or current distributions on the surface of reflectors. The antenna types which the program is set up to synthesize are: line source, rectangular aperture, circular aperture, linear array, rectangular array, and arbitrary planar array.

Stutzman, W. L.↗

Interior radiances in optically deep absorbing media. 3: Scattering from Haze L

The interior radiances are calculated within an optically deep absorbing medium scattering according to the Haze L phase function. The dependence on the solar zenith angle, the single scattering albedo, and the optical depth within the medium is calculated by the matrix operator method. The development of the asymptotic angular distribution of the radiance in the diffusion region is illustrated through a number of examples; it depends only on the single scattering albedo and on the phase function for single scattering. The exact values of the radiance in the diffusion region are compared with values calculated from the approximate equations proposed by Van de Hulst. The variation of the radiance near the lower boundary of an optically thick medium is illustrated with examples. The attenuation length is calculated for various single scattering albedos and compared with the corresponding values for Rayleigh scattering. The ratio of the upward to the downward flux is found to be remarkably constant within the medium. The heating rate is calculated and found to have a maximum value at an optical depth of two within a Haze L layer when the sun is at the zenith.

Kattawar, G. W.↗

Supersonic laminar base pressure, heat-transfer, and upstream influence correlations for small steps

Some results of the application of a unified theory of small disturbances in nonuniform high-speed boundary layers to the small-step problem for the case of supersonic laminar nonadiabatic flow are described. Comparisons with experimental data are given which show good agreement in a variety of physical features, including useful new base pressure and heat-transfer correlation relations. Using a small perturbation approach, an analytical theory is developed that predicts the steady-state pressure, temperature, and heat-transfer disturbance fields including suction or blowing through the step. Particular emphasis is placed on the important effects of the highly nonuniform flow across the boundary layer, lateral pressure gradients, and the upstream influence associated with viscous-inviscid interaction. Approximate analytical expressions for the wall pressure, skin friction, and heat-transfer perturbations and upstream influence distance are derived by a Fourier transformation approach.

Inger, G. R.↗

Interior radiances in optically deep absorbing media. III Scattering from Haze L

The interior radiances are calculated within an optically deep absorbing medium scattering according to the Haze L phase function. The dependence on the solar zenith angle, the single scattering albedo, and the optical depth within the medium is calculated by the matrix operator method. The development of the asymptotic angular distribution of the radiance in the diffusion region is illustrated through a number of examples; it depends only on the single scattering albedo and on the phase function for single scattering. The exact values of the radiance in the diffusion region are compared with values calculated from the approximate equations proposed by Van de Hulst. The variation of the radiance near the lower boundary of an optically thick medium is illustrated with examples. The attenuation length is calculated for various single scattering albedos and compared with the corresponding values for Rayleigh scattering. The ratio of the upward to the downward flux is found to be remarkably constant within the medium.

Kattawar, G. W.↗

Crack-opening displacements in center-crack, compact, and crack-line wedge-loaded specimens

The theoretical crack-opening displacements for center-crack, compact, and crack-line wedge-loaded specimens (reported in the ASTM Proposed Recommended Practice for R-Curve Determination (1974)) disagree with experimental measurements in the literature. The disagreement is a result of using approximate specimen configurations and load representation to obtain the theoretical displacements. An improved method of boundary collocation is presented which was used to obtain the theoretical displacements in these three specimen types; the actual specimen configurations and more accurate load representation were used. In the analysis of crack-opening displacements in the compact and crack-line wedge-loaded specimens, the effects of the pin-loaded holes were also included. The theoretical calculations agree with the experimental measurements reported in the literature. Also examined are accurate polynomial expressions for crack-opening displacements in both compact and crack-line wedge-loaded specimens.

Newman, J. C., Jr.↗

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗

Approximate convective heating equations for hypersonic flows

Laminar and turbulent heating-rate equations appropriate for engineering predictions of the convective heating rates about blunt reentry spacecraft at hypersonic conditions are developed. The approximate methods are applicable to both nonreacting and reacting gas mixtures for either constant or variable-entropy edge conditions. A procedure which accounts for variable-entropy effects and is not based on mass balancing is presented. Results of the approximate heating methods are in good agreement with existing experimental results as well as boundary-layer and viscous-shock-layer solutions.

Zoby, E. V.↗

Three-dimensional Euler solutions for long-duct nacelles

A three-dimensional Euler-equation computational technique has been developed to solve for the transonic flow past flow-through nacelles. The technique employs an approximately-factored alternating-direction implicit numerical algorithm and a radiation treatment of the outflow boundary. Studies are presented which show that the radiation treatment gives better numerical convergence than the condition of specifying the pressure at the outflow boundary. Calculations made with the technique are presented for a long-duct turbofan engine nacelle at a Mach number of 0.80 and angles of attack of 0 deg and 4 deg. Good agreement is shown between the computational results and wind-tunnel data. Problem areas are identified and recommendations are made for further numerical studies.

Compton, W. B., III↗

Single scattering from nonspherical Chebyshev particles: A compendium of calculations

A large set of exact calculations of the scattering from a class of nonspherical particles known as Chebyshev particles' has been performed. Phase function and degree of polarization in random orientation, and parallel and perpendicular intensities in fixed orientations, are plotted for a variety of particles shapes and sizes. The intention is to furnish a data base against which both experimental data, and the predictions of approximate methods, can be tested. The calculations are performed with the widely-used Extended Boundary Condition Method. An extensive discussion of this method is given, including much material that is not easily available elsewhere (especially the analysis of its convergence properties). An extensive review is also given of all extant methods for nonspherical scattering calculations, as well as of the available pool of experimental data.

Wiscombe, W. J.↗

The construction and use of divergence free vector expansions for incompressible fluid flow calculations

For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.

Mhuiris, N. M. G.↗

Three-dimensional computation of mixing of transverse injector in a ducted supersonic airstream

Numerical solutions of the three-dimensional mass-averaged elliptic Navier-Stokes equations, including species transport, are obtained for nonreacting, turbulent, mixing flow fields for the case of transverse sonic injection of a secondary gas into a supersonic airstream through a circular orifice injector. Results are presented for flow through a constant area duct and through a duct with a rearward-facing step upstream of the injector. The equations are numerically integrated using MacCormack's explicit method and turbulence is included using the Baldwin-Lomax algebraic eddy viscosity model. In the species transport and energy equations, diffusion coefficients based on Fick's law and an assumption of unit Lewis number are applied. The computations were performed on a CDC-VPS-32 (extended version of Cyber-205) using a grid consisting of approximately 200,000 points. The computed results are compared with experimentally observed penetration and spreading boundaries for an injected gas at two dynamic pressure ratios. Three-dimensional flow field structures are dipicted in terms of static pressure, mass fractions of species and velocity vectors.

Uenishi, K.↗

BRYNTRN: A baryon transport model

The development of an interaction data base and a numerical solution to the transport of baryons through an arbitrary shield material based on a straight ahead approximation of the Boltzmann equation are described. The code is most accurate for continuous energy boundary values, but gives reasonable results for discrete spectra at the boundary using even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H2O). The resulting computer code is self-contained, efficient and ready to use. The code requires only a very small fraction of the computer resources required for Monte Carlo codes.

Wilson, John W.↗

Planetary circulations in the presence of transient and self-induced heating

The linearized global primitive equation model of Salby and Garcia (1987) was modified to study tropospheric circulations and their interaction with tropical convection. The vertical resolution was increased to approximately 1.5 km, with the effective vertical domain extending the approximately 40 km (higher altitudes in the calculation are used to ensure the radiation upper boundary condition is satisfied). A formal wave conditional instability of the second kind (CISK) parameterization was introduced. That parameterization allows low-frequency disturbances to interact with convection by organizing surface moisture. Boundary layer friction has two important effects in these calculations: it exaggerates surface convergence near the equator by driving the flow out of geostrophic balance; and it introduces a phase shift between the convergence pattern and the circulation, which drives the heating out of quadrature with the temperature field and allows a positive feedback to take place. The observational complement of this theoretical work involves Global Cloud Imagery (GCI) composited from six satellites in the International Satellite Cloud Climatology Project (ISCCP) which were simultaneously observing the Earth's cloud field. Having global coverage, horizontal resolution of half a degree, and produced synoptically every 3 hours, the GCI affords an unprecendented view of global convection. A full year (1983 to 1984) of GCI has now been successfully created. Among other phenomena, that imagery reveals a regular progression of easterly waves in the convective pattern of the Western Pacific. Those disturbances have odd symmetry about the equator and organize convection into a pattern with similar symmetry. In time-mean maps of cloud cover, these propagating disturbances produce a split Inter-Tropical Convergence Zone (ITCZ) in the western Pacific.

Salby, Murry L.↗

Experimental study of boundary layer transition with elevated freestream turbulence on a heated flat plate

A detailed investigation to document momentum and thermal development of boundary layers undergoing natural transition on a heated flat plate was performed. Experimental results of both overall and conditionally sampled characteristics of laminar, transitional, and low Reynolds number turbulent boundary layers are presented. Measurements were acquired in a low-speed, closed-loop wind tunnel with a freestream velocity of 100 ft/s and zero pressure gradient over a range of freestream turbulence intensities (TI) from 0.4 to 6 percent. The distributions of skin friction, heat transfer rate and Reynolds shear stress were all consistent with previously published data. Reynolds analogy factors for R(sub theta) is less than 2300 were found to be well predicted by laminar and turbulent correlations which accounted for an unheated starting length. The measured laminar value of Reynolds analogy factor was as much as 53 percent higher than the Pr(sup -2/3). A small dependence of turbulent results on TI was observed. Conditional sampling performed in the transitional boundary layer indicated the existence of a near-wall drop in intermittency, pronounced at certain low intermittencies, which is consistent with the cross-sectional shape of turbulent spots observed by others. Non-turbulent intervals were observed to possess large magnitudes of near-wall unsteadiness and turbulent intervals had peak values as much as 50 percent higher than were measured at fully turbulent stations. Non-turbulent and turbulent profiles in transitional boundary layers cannot be simply treated as Blasius and fully turbulent profiles, respectively. The boundary layer spectra indicate predicted selective amplification of T-S waves for TI is approximately 0.4 percent. However, for TI is approximately 0.8 and 1.1 percent, T-S waves are localized very near the wall and do not play a dominant role in transition process.

Sohn, Ki-Hyeon↗

Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

Deanna, R. G.↗

Winds from accretion disks - Ultraviolet line formation in cataclysmic variables

Winds from accretion disks in cataclysmic variable stars are ubiquitous. Observations by IUE reveal P Cygni-shaped profiles of high-ionization lines which are attributed to these winds. We have studied the formation of UV emission lines in cataclysmic variables by constructing kinematical models of biconical rotating outflows from disks around white dwarfs. The photoionization in the wind is calculated taking into account the radiation fields of the disk, the boundary layer, and the white dwarf. The 3D radiative transfer is solved in the Sobolev approximation. Effects on the line shapes of varying basic physical parameters of the wind are shown explicitly. We identify and map the resonant scattering regions in the wind which have strongly biconical character regardless of the assumed velocity and radiation fields. Rotation at the base of the wind introduces a radial shear which decreases the line optical depth and reduces the line core intensity. We find that it is possible to reproduce the observed P Cygni line shapes and make some predictions to be verified in high-resolution observations.

Shlosman, Isaac↗

A discussion of cone and flat-plate Reynolds numbers for equal ratios of the laminar shear to the shear caused by small velocity fluctuations in a laminar boundary layer

By use of the linear theory of boundary-layer stability and Schlichting's formula for the maximum amplification of a disturbance, an approximate relation is derived between the Reynolds number on a cone and the Reynolds number on a flat plate for equal closeness to transition. The indication is that the ratio of the cone Reynolds number for transition, based on the distance to the cone apex, to the plate Reynolds number for transition, based on the distance to the leading edge, is not in general equal to 3, as has been suggested by other investigators, but varies from 3 when transition occurs at the minimum critical Reynolds number to unity when transition occurs at a large multiple of the critical Reynolds number.

FLOW, TURBULENT↗