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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Generation of auroral kilometric radiation and the structure of auroral acceleration region

Generation of auroral kilometric radiation (AKR) in the auroral acceleration region is studied. It is shown that auroral kilometric radiation can be generated by backscattered electrons trapped in the acceleration region via a cyclotron maser process. The parallel electric field in the acceleration region is required to be distributed over 1-2 earth radii. The observed AKR frequency spectrum can be used to estimate the altitude range of the auroral acceleration region. The altitudes of the lower and upper boundaries of the acceleration region determined from the AKR data are respectively approximately 2000 and 9000 km.

Lee, L. C.↗
NASA NTRS1982

Aerothermodynamic environment of a Titan aerocapture vehicle

The extent of convective and radiative heating for a Titan aerocapture vehicle is investigated. The flow in the shock layer is assumed to be axisymmetric, steady, viscous, and compressible. It is further assumed that the gas is in chemical and local thermodynamic equilibrium and tangent slab approximation is used for the radiative transport. The effect of the slip boundary conditions on the body surface and at the shock wave are included in the analysis of high-altitude entry conditions. The implicit finite difference techniques is used to solve the viscous shock-layer equations for a 45 degree sphere cone at zero angle of attack. Different compositions for the Titan atmosphere are assumed, and results are obtained for the entry conditions specified by the Jet Propulsion Laboratory.

Tiwari, S. N.↗
NASA NTRS1983

Sound radiation from a vibrating cylinder in flow

An alternative approach to the conventional manner of calculating net power output of a source by integrating the energy flux over a distant control surface is presented. The method gives a means of calculating source radiation by an integral only over the source region and eliminates the need for a possibly difficult flux integration. The radiator chosen here is a finite, open ended, cylindrical shell of negligible thickness in the presence of uniform axial flow. A homogeneous medium without external boundaries is assumed. An integral representation and, in the long wave length approximation, an analytical expression for the power radiated by the vibrating cylinder are obtained. The results show power radiation may be increased or decreased by the presence of flow, depending on the relationship between the wave number of the vibration and the size of the radiator. The result is compared with results for simpler geometries.

Smith, M.↗
NASA NTRS1984

Implicit treatment of the unsteady full potential equation in conservation form

An implicit, conservative treatment for the unsteady full potential equation in two-dimensions is presented. The method employs a local time linearization for density, and introduces flux biasing concepts based on sonic conditions for the generation of artificial viscosity to capture shocks without any overshoots. The boundary condition is treated implicitly using a splitting procedure consistent with the approximate factorization scheme. This allows for extremely large Courant numbers, even for nonorthogonal grid at the body. The method has application not only to unsteady problems, but also to generate the starting blunt body solution for a supersonic full potential marching code. Results are presented for flows over cylinders, spheres and airfoils. Comparisons are made with available Euler and full potential results, and are in excellent agreement.

Shankar, V.↗
NASA NTRS1985

Modelling Equilibrium and Fractional Crystallization in the System MgO-FeO-CaO-Al2O3-SiO2

A mathematical modelling technique for use in petrogenesis calculations in the system MgO-FeO-CaO-Al2O3-SiO2 is reported. Semiempirical phase boundary and elemental distribution information was combined with mass balance to compute approximate equilibrium crystallization paths for arbitrary system compositions. The calculation is applicable to a range of system compositions and fractionation calculations are possible. The goal of the calculation is the computation of the composition and quantity of each phase present as a function of the degree of solidification. The degree of solidification is parameterized by the heat released by the solidifying phases. The mathematical requirement for the solution of this problem is: (1) An equation constraining the composition of the magma for each solid phase in equilibrium with the liquidus phase, and (2) an equation for each solid phase and each component giving the distribution of that element between that phase and the magma.

Herbert, F.↗
NASA NTRS1985

Stability of interfaces with mesh refinement

A study is conducted of the stability of mesh refinement in space and time for several different interface equations and finite-difference approximations. First, a root condition which implies stability for the initial-boundary value problem for this type of interface is derived. From the root condition, the stability of several interface equations is proved, using the maximum principle. In some cases, the final verification steps can be done analytically; in other cases, a simple computer program has been written to check the condition for values of a parameter along the boundary of the unit circle. Using this method, stability for Lax-Wendroff with all the interface conditions considered, and for Leapfrog with interpolation interface conditions when the fine and coarse grids overlap is proved.

Berger, M. J.↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
NASA NTRS1987

Experimental cavity pressure distributions at supersonic speeds

An investigation was conducted to define pressure distributions for rectangular cavities over a range of free-stream Mach numbers and cavity dimensions. These pressure distributions together with schlieren photographs are used to define the critical values of cavity length-to-depth ratio that separate open type cavity flows from closed type cavity flows. For closed type cavity flow, the shear layer expands over the cavity leading edge and impinges on the cavity floor, whereas for open type cavity flow, the shear layer bridges the cavity. The tests were conducted by using a flat-plate model permitting the cavity length to be remotely varied from 0.5 to 12 in. Cavity depths and widths were varied from 0.5 to 2.5 in. The flat-plate boundary layer approaching the cavity was turbulent and had a thickness of approximately 0.2 in. at the cavity front face for the range of test Mach numbers from 1.5 to 2.86. Presented are a discussion of the results and a complete tabulation of the experimental data.

Stallings, Robert L., Jr.↗
NASA NTRS1987

Computational fluid dynamics applications at McDonnel Douglas

Representative examples are presented of applications and development of advanced Computational Fluid Dynamics (CFD) codes for aerodynamic design at the McDonnell Douglas Corporation (MDC). Transonic potential and Euler codes, interactively coupled with boundary layer computation, and solutions of slender-layer Navier-Stokes approximation are applied to aircraft wing/body calculations. An optimization procedure using evolution theory is described in the context of transonic wing design. Euler methods are presented for analysis of hypersonic configurations, and helicopter rotors in hover and forward flight. Several of these projects were accepted for access to the Numerical Aerodynamic Simulation (NAS) facility at the NASA-Ames Research Center.

Hakkinen, R. J.↗
NASA NTRS1987

Comparing the jerk with other global models of the geomagnetic field from 1960 to 1978

About 3300 satellite values of geomagnetic intensity and about 700 observatory values of annual mean magnetic vector components from 1960 to 1978 were fitted by three global models of the geomagnetic field B. Each model includes a spatially constant external field whose time dependence is a constant plus another constant times the Dst index, and each model accepts a time-independent station correction at each observatory. The time dependence of the internal Gauss coefficients is either cubic, quintic, or biquadratic (two independent quadratics, one before and one after January 1, 1970); and g1(0) also has an induced term proportional to the Dst index. The rms residual of the data fit is the same for the cubic and biquadratic models and insignificantly smaller for the quintic model. The quintic and biquadratic models have 1164 adjustable parameters, and the cubic has 1038. At a high level of significance the parameters of the best fitting biquadratic rule out a physical model for the magnetic impulse of 1969 in which the level surfaces of electrical conductivity in the lower mantle are approximately spherical, and the radial magnetic field at the core-mantle boundary goes from one quadratic time dependence to another in a year or less.

Backus, G. E.↗
NASA NTRS1989

Supercritical flow past a symmetrical bicircular arc airfoil

A numerical scheme is developed for computing steady supercritical flow about symmetrical airfoils, applying it to an ellipse for zero angle of attack. An algorithmic description of this new scheme is presented. Application to a symmetrical bicircular arc airfoil is also proposed. The flow field before the shock is region 1. For transonic flow, singularity can be avoided by integrating the resulting ordinary differential equations away from the body. Region 2 contains the shock which will be located by shock fitting techniques. The shock divides region 2 into supersonic and subsonic regions and there is no singularity problem in this case. The Method of Lines is used in this region and it is advantageous to integrate the resulting ordinary differential equation along the body for shock fitting. Coaxial coordinates have to be used for the bicircular arc airfoil so that boundary values on the airfoil body can be taken with one direction of the coaxial coordinates fixed. To avoid taking boundary values at + or - infinity in the coaxial co-ordinary system, approximate analytical representation of the flow field near the tips of the airfoil is proposed.

Holt, Maurice↗
NASA NTRS1988

Environmental effects of large impacts on the earth; relation to extinction mechanisms

Since Alvarez et al., discovered a worldwide approx. cm-thick layer of fine sediments laden with platinum group elements in approximately chondritic proportions exactly at the Cretaceous-Tertiary (C-T) boundary, and proposed bolide-impact as triggering mass extinctions, many have studied this hypothesis and the layer itself with its associated spherules and shocked quartz. At issue is whether the mass extinctions, and this horizon has an impact versus volcanic origin. A critical feature of the Alvarez hypothesis is the suggestion that the bolide or possibly a shower of objects delivered to the earth approx. 0.6 x 10 to the 18th power g of material which resulted in aerosol-sized ejecta such that global insolation was drastically reduced for significant periods. Such an event would lower temperatures on continents and halt photosynthesis in the upper 200 m of th eocean. The latter would strangle the marine food chain and thus produce the major marine faunal extinctions which mark the C-T boundary. Crucial issues examined include: What are the dynamics of atmospheric flow occurring upon impact of a large bolide with the earth; What is the size distributions of the very fine impact ejecta and how do these compare to the models of ejecta which are used to model the earth's radiative thermal balance. The flow field due to passage of a 10 km diameter bolide through an exponential atmosphere and the interaction of the gas flow and bolide with the solid ear was calculated. The CO2 released upon impact onto shallow marine carbonate sections was modeled and found that the mass of CO2 released exceeds the present 10 to the 18th power g CO2 budget of the earth's atmosphere by several times. Using the calculations of Kasting and Toon it was found that to compute the temperature rise of the earth's surface as a function of CO2 content, it was found that sudden and prolonged global increases are induced from impact of 20 to 50 km radius projectiles and propose that sudden terrestrial greenhouse-induced heating, not cooling, produced the highly variable extinctions seen at the C-T boundary.

Okeefe, John D.↗
NASA NTRS1987

Current distribution in a three-dimensional IC analyzed by a perturbation method. Part 1: A simple steady state theory

The steady state current distribution in a three dimensional integrated circuit is presented. A device physics approach, based on a perturbation method rather than an equivalent lumped circuit approach, is used. The perturbation method allows the various currents to be expressed in terms of elementary solutions which are solutions to very simple boundary value problems. A Simple Steady State Theory is the subtitle because the most obvious limitation of the present version of the analysis is that all depletion region boundary surfaces are treated as equipotential surfaces. This may be an adequate approximation in some applications but it is an obvious weakness in the theory when applied to latched states. Examples that illustrate the use of these analytical methods are not given because they will be presented in detail in the future.

Edmonds, Larry D.↗
NASA NTRS1988

Methods for analysis of cracks in three-dimensional solids

Various methods used for determining stress-intensity factors for cracked three-dimensional bodies are reviewed. The review is limited to the determination of mode-I stress-intensity factors. Some exact solutions for cracks in infinite solids are presented, and various approximate methods that have been used for the solution of the boundary-value problem of finite solids with cracks are addressed. The techniques used to extract the stress-intensity factors from these solutions are considered, and various methods used to calculate the stress-intensity factors for through-the-thickness cracks, semielliptical surface cracks, and quarter-elliptical corner cracks at holes are compared for the case of remote tensile loading.

Raju, I. S.↗
NASA NTRS1991

Symmetry plane model for turbulent flows with vortex generators

An approximate procedure is proposed for predicting the performance of counterrotating vortex-generator installations in incompressible flow. An inviscid calculation that includes the motion of the vortices is used to obtain crossflow velocities at the boundary-layer edge as a function of initial position, spacing, and strength of the vortices, and local values of the spanwise gradient are then folded into an integral turbulent-boundary layer procedure applied in the plane of symmetry. Special attention is paid to the consistency of the approximations and equations used. The two-dimensional aerodynamics of vortex generator installations on a NACA 0016 airfoil at angle-of-attack are estimated in this manner, and the results compared with experiments carried out with a 30-cm chord wing mounted in a 2.4 x 3.6-m cross-section wind tunnel and tested at chord Reynolds numbers of 0.7 and 1.4 x 10 to the 6th. Agreement in the separation location is found for these complex flows for a range of conditions.

Arnaud, Gilles L.↗
NASA NTRS1991

Evolution on large-scale plasma structures in comets: Kinematics and physics

The disconnection event (DE) consists of the periodic loss of a comet's entire plasma tail and the growth of a new one. This spectacular phenomenon is not understood. The strategy was to assemble a data base of specific events studied in detail, determine the solar wind conditions responsible for DEs, and develop a consistent physical model. Analysis is complete for the sequence of DEs that took place during 13 to 18 April 1986. The first DE correlates well with a sector boundary crossing for the comet and a group of DEs that occurred approximately a day later could have produced by polarity reversals seen in the IMP-8 data. Thus, these DEs are consistent with the frontside, magnetic reconnection mechanism.

Brandt, John C.↗
NASA NTRS1992

UV line diagnostics of accretion disk winds in cataclysmic variables

The IUE data base is used to analyze the UV line shapes of cataclysmic variables RW Sex, RW Tri, and V Sge. Observed lines are compared to synthetic line profiles computed using a model of rotating bi-conical winds from accretion disks. The wind model calculates the wind ionization structure self-consistently including photoionization from the disk and boundary layer and treats 3-D line radiation transfer in the Sobolev approximation. It is found that winds from accretion disks provide a good fit for reasonable parameters to the observed UV lines which include the P Cygni profiles for low inclination systems and pure emission at large inclination. Disk winds are preferable to spherical winds which originate on the white dwarf because they (1) require a much lower ratio of mass loss rate to accretion rate and are therefore more plausible energetically, (2) provide a natural source for a bi-conical distribution of mass outflow which produces strong scattering far above the disk leading to P Cygni profiles for low inclination systems, and pure line emission profiles at high inclination with the absence of eclipses in UV lines, and (3) produce rotation broadened pure emission lines at high inclination.

Vitello, Peter↗