Search NASA⌕ Search

SEARCH · Search NASA

Results for “Embedded 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Multiple grid problems on concurrent-processing computers

Three computer codes were studied which make use of concurrent processing computer architectures in computational fluid dynamics (CFD). The three parallel codes were tested on a two processor multiple-instruction/multiple-data (MIMD) facility at NASA Ames Research Center, and are suggested for efficient parallel computations. The first code is a well-known program which makes use of the Beam and Warming, implicit, approximate factored algorithm. This study demonstrates the parallelism found in a well-known scheme and it achieved speedups exceeding 1.9 on the two processor MIMD test facility. The second code studied made use of an embedded grid scheme which is used to solve problems having complex geometries. The particular application for this study considered an airfoil/flap geometry in an incompressible flow. The scheme eliminates some of the inherent difficulties found in adapting approximate factorization techniques onto MIMD machines and allows the use of chaotic relaxation and asynchronous iteration techniques. The third code studied is an application of overset grids to a supersonic blunt body problem. The code addresses the difficulties encountered when using embedded grids on a compressible, and therefore nonlinear, problem. The complex numerical boundary system associated with overset grids is discussed and several boundary schemes are suggested. A boundary scheme based on the method of characteristics achieved the best results.

Eberhardt, D. S.↗
NASA NTRS1988

Generation of dynamo magnetic fields in protoplanetary and other astrophysical accretion disks

A computational method for treating the generation of dynamo magnetic fields in astrophysical disks is presented. The numerical difficulty of handling the boundary condition at infinity in the cylindrical disk geometry is overcome by embedding the disk in a spherical computational space and matching the solutions to analytically tractable spherical functions in the surrounding space. The lowest lying dynamo normal modes for a 'thick' astrophysical disk are calculated. The generated modes found are all oscillatory and spatially localized. Tha potential implications of the results for the properties of dynamo magnetic fields in real astrophysical disks are discussed.

Stepinski, T. F.↗
NASA NTRS1998

The Aircraft Morphing Program

In the last decade smart technologies have become enablers that cut across traditional boundaries in materials science and engineering. Here we define smart to mean embedded actuation, sensing, and control logic in a tightly coupled feedback loop. While multiple successes have been achieved in the laboratory, we have yet to see the general applicability of smart devices to real aircraft systems. The NASA Aircraft Morphing program is an attempt to couple research across a wide range of disciplines to integrate smart technologies into high payoff aircraft applications. The program bridges research in seven individual disciplines and combines the effort into activities in three primary program thrusts. System studies are used to assess the highest- payoff program objectives, and specific research activities are defined to address the technologies required for development of smart aircraft systems. In this paper we address the overall program goals and programmatic structure, and discuss the challenges associated with bringing the technologies to fruition.

Wlezien, R. W.↗
NASA NTRS1992

A boundary element alternating method for two-dimensional mixed-mode fracture problems

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Raju, I. S.↗
NASA NTRS1989

Boundary value problems in magnetohydrodynamics (and fluid dynamics). I - Radiation boundary condition

In view of the difficulties that emerge in numerical simulations of time-dependent phenomena involving MHD waves due to waves' artificial reflection at numerical boundaries, where time-dependent magnetic fields outside those boundaries are unknown, a method previously developed for the determination of radiative boundary conditions is generalized for the treatment of isothermal MHD waves. The illustrative problem thus treated is that of the magnetic braking of a rigidly rotating disk which is embedded in a differentially-rotating fluid. The analytic results obtained are compared with those from numerical simulations with and without the present boundary treatment.

Vanajakshi, T. C.↗
NASA NTRS1984

Magnetic structure of the distant geotail from -60 to -220 earth radii - ISEE-3

ISEE-3 magnetic-field measurements in the region of the geomagnetic tail from -80 to -220 earth radii are reported and discussed. A well-ordered field structure is found, comprising two 7-8-nT lobes separated by a plasma sheet, an embedded neutral sheet with significant By fields, and an intermittent plasma-sheet boundary layer with 5-nT-amplitude (peak-to-peak) electromagnetic waves. The plasma-sheet Bz distribution changes from principally northern orientation near the earth to an approximately equal north-south distribution at 200-220 earth radii. These findings are considered to be in general agreement with magnetic-reconnection models of the magnetosphere, with reconnection either throughout the region observed (in tearing-mode or plasmoid-formation models) or at a constant (about 220-earth-radii) or variable (40-80 to 220-earth-radii) X line (in X-line models).

Tsurutani, B. T.↗
NASA NTRS1989

An Euler analysis of a High-Speed Civil Transport (HSCT) concept at Mach 3

A marching Euler solver, GEM3D, was used to predict the Mach 3 flow field for the wing and body of a High-Speed Civil Transport concept. The analysis focused on a typical cruise lift coefficient of 0.1 at alpha = 3 deg. The Euler solution indicated that embedded shocks formed on the upper surface of the inboard wing panel and at the leading-edge of the outboard wing panel, due to its supersonic leading edge condition. According to a simple static-pressure criterion, the embedded wing upper-surface shocks were sufficiently strong to separate a turbulent boundary layer. Comparison of aerodynamic coefficients from the Euler solver with those from linear theory shows that the linear theory estimates of lift and drag are optimistic, which would lead to optimistic estimates of cruise range.

Pittman, James L.↗
NASA NTRS2002

Finite-Difference Time Domain Solution of Light Scattering and Absorption by Particles in an Absorbing Medium

The three-dimensional (3-D) finite-difference time-domain (FDTD) technique has been extended to simulate light scattering and absorption by nonspherical particles embedded in an absorbing dielectric medium. A uniaxial perfectly matched layer (UPML) absorbing boundary condition is used to truncate the computational domain. When computing the single-scattering properties of a particle in an absorbing dielectric medium, we derive the single-scattering properties including scattering phase functions, extinction, and absorption efficiencies using a volume integration of the internal field. A Mie solution for light scattering and absorption by spherical particles in an absorbing medium is used to examine the accuracy of the 3-D UPML FDTD code. It is found that the errors in the extinction and absorption efficiencies from the 3-D UPML FDTD are less than similar to 2%. The errors in the scattering phase functions are typically less than similar to 5%. The errors in the asymmetry factors are less than similar to 0.l%. For light scattering by particles in free space, the accuracy of the 3-D UPML FDTD scheme is similar to a previous model.

Sun, W.↗
NASA NTRS2011

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗
NASA NTRS1980

Supersonic flow over three-dimensional ablated nosetips using an unsteady implicit numerical procedure

The three-dimensional supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical domain, containing multiple zones of embedded subsonic flow and separated flow regions into a rectangular computational volume at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing numerical solutions and experimental data for simple spheres which validates the program. Results are also presented for analytically defined indented bodies for both laminar and turbulent flow conditions that demonstrate the program's capability for computing such flows.

Kutler, P.↗
NASA NTRS1981

Numerical analysis of the scramjet-inlet flow field by using two-dimensional Navier-Stokes equations

A computer code was developed to solve the full two dimensional Navier-Stokes equations in a supersonic combustion ramjet (scramjet) inlet. In order to be able to consider a general inlet geometry with embedded bodies, a numerical coordinate transformation is used which generates a set of boundary-fitted curvilinear coordinates. The explicit finite difference algorithm of MacCormack is used to solve the governing equations. An algebraic, two-layer eddy-viscosity model is used for the turbulent flow. The code can analyze both inviscid and viscous flows with no strut, one strut, or multiple struts in the flow field. The application of the two dimensional analysis in the preliminary parametric design studies of a scramjet inlet is discussed. Detailed results are presented for one model problem and for several actual scramjet-inlet configurations.

Kumar, A.↗
NASA NTRS1982

Transonic potential flow and coordinate generation for bodies in a wind tunnel

An accurate method has been developed for computation of transonic potential flow about a 2-D lifting airfoil or an axisymmetric body in a wind tunnel. The computational mesh is nearly orthogonal everywhere and is generated by a sequence of Schwarz-Christoffel transformations and shearings to obtain an '0' grid near the body. A conservative finite-volume scheme for the full-potential equation and exact boundary conditions is used together with 'retarded' density to solve efficiently transonic flow with embedded shocks and large regions of supersonic flow. Supersonic free-stream flows are also solved with captured bow shocks and embedded subsonic regions.

Doria, M. L.↗
NASA NTRS1984

Calculating C-grids with fine and embedded mesh regions

A program for calculating a C-type mesh around airfoil like shapes is described. The Jameson/Caughey approach is used: a parabolic transformation coupled with a shearing transformation. The algebraic algorithm is capable of efficiently generating nearly orthogonal grids. A high degree of grid control is possible. The user may specify grid boundaries, number of grid lines, and location of (and spacing in) trailing edge and leading edge fine mesh areas. The capability of embedding fine mesh regions, for use with new adaptive grid techniques, is being developed. Grids generated by the program were used in Euler flow flow calculatons by W. Usab. Results superior to results calculated on previous O-type grids were obtained. Specifically, calculations converged faster using C-grids rather than 0-grids, total pressure loss spikes at the trailing edge of the airfoil were eliminated, and the Ni method converged with zero artificial smoothing for a subcritical case (resulting overall total pressure loss was then nearly zero). These improvements were attributed to higher grid orthogonality, especially at the trailing edge. The program itself is fairly straightforward. Roughly half of the 800 code lines are comment lines.

Loyd, B.↗
NASA NTRS1984

Numerical algorithm conjugating steady and transient, separated, compressible flow and a solid body having arbitrarily distributed heat sources

A numerical algorithm that is second-order accurate in time has been developed for the conjugated problem of a separated, compressible flow field and a conductive solid body. The full two-dimensional time-dependent Navier-Stokes equations are coupled with the time-dependent energy equation for the solid body and are solved simultaneously. using implicit algorithms. The energy equation for the solid body may include arbitrarily distributed heat sources. The algorithm has been exmined for the case of two-dimensional supersonic compression-corner interaction, with a heat source embedded in the wall in the vicinity of the separation bubble and the attached boundary layer. The effect of the heat source on the flow field is studied for steady and transient cases.

Degani, D.↗
NASA NTRS1992

CFD analysis of pump consortium impeller

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design impact in a productive manner. The main goal of this study is to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A Navier-Stokes flow solver, FDNS, embedded with the extended k-epsilon turbulence model and with appropriate moving interface boundary conditions, is developed to analyze turbulent flows in the turbomachinery devices. The FDNS code was benchmarked with its numerical predictions of the pump consortium inducer, and provides satisfactory results. In the present study, a CFD analysis of the pump consortium impeller will be conducted with the application of the FDNS code. The pump consortium impeller, with partial blades, is the new design concept of the advanced rocket engine.

Cheng, Gary C.↗
NASA NTRS1993

Birkeland currents in the plasma sheet

A search was conducted for the signatures of Birkeland currents in the Earth's magnetic tail, using observed values of B(sub x) and B(sub y) from large sets of spacecraft data. The data were binned by x and y for -10 greater than x(sub GSM) greater than -35 and absolute value of y(sub GSM) less than or equal to 20 R(sub E) (less than or equal to 30 R(sub E) for x(sub GSM) less than or equal to -25 R(sub E)) and in each bin their distribution in the (B(sub x), B(sub y)) plane was fitted by least squares to a piecewise linear function. That gave average x-y distributions of the flaring angle between B(sub xy) and the x direction, as well as that angle's variation across the thickness of the plasma sheet. Angles obtained in the central plasma sheet differed from those derived near the lobe boundary. That is the expected signature if earthward or tailward Birkeland current sheets are embedded in the plasma sheet, and from this dfiference we derived the dawn-dusk profiles of the tail Birkeland currents for several x(sub GSM) intervals. It was found that (1) the Birkeland currents have the sense of region 1 currents, when mapped to the ionosphere; (2) both the linear current density (kiloamperes/R(sub E)) and the net magnitude of the field-aligned currents decrease rapidly down the tail; (3) the total Birkeland current at x approximately equals -10 R(sub E) equals approximately equals 500-700 kA, which is approx. 30% of the net region 1 current observed at ionospheric altitudes, in agreement with model mapping results; and (4) the B(sub z) and B(sub y) components of the interplanetary magnetic field influence the distribution of Birkeland currents in the tail.

Tsyganenko, Nikolai A.↗
NASA NTRS1999

Automated Development of Accurate Algorithms and Efficient Codes for Computational Aeroacoustics

The simulation of sound generation and propagation in three space dimensions with realistic aircraft components is a very large time dependent computation with fine details. Simulations in open domains with embedded objects require accurate and robust algorithms for propagation, for artificial inflow and outflow boundaries, and for the definition of geometrically complex objects. The development, implementation, and validation of methods for solving these demanding problems is being done to support the NASA pillar goals for reducing aircraft noise levels. Our goal is to provide algorithms which are sufficiently accurate and efficient to produce usable results rapidly enough to allow design engineers to study the effects on sound levels of design changes in propulsion systems, and in the integration of propulsion systems with airframes. There is a lack of design tools for these purposes at this time. Our technical approach to this problem combines the development of new, algorithms with the use of Mathematica and Unix utilities to automate the algorithm development, code implementation, and validation. We use explicit methods to ensure effective implementation by domain decomposition for SPMD parallel computing. There are several orders of magnitude difference in the computational efficiencies of the algorithms which we have considered. We currently have new artificial inflow and outflow boundary conditions that are stable, accurate, and unobtrusive, with implementations that match the accuracy and efficiency of the propagation methods. The artificial numerical boundary treatments have been proven to have solutions which converge to the full open domain problems, so that the error from the boundary treatments can be driven as low as is required. The purpose of this paper is to briefly present a method for developing highly accurate algorithms for computational aeroacoustics, the use of computer automation in this process, and a brief survey of the algorithms that have resulted from this work. A review of computational aeroacoustics has recently been given by Lele.

Goodrich, John W.↗
NASA NTRS2018

Use of Simple Continuum Solutions in Finite Element Alternating Method for Fracture Problems

The performance of the finite element alternating (FEAM) method for two-dimensional crack problems is studied with respect to a polynomial pressure distribution fitted to the crack face stresses. The FEAM alternates between the analytical solution of crack in an infinite plate subjected to arbitrary polynomial distribution and a finite element solution of an uncracked body to satisfy the required boundary conditions in the crack problem. In this paper, the FEAM is applied to embedded crack and edge crack problems. For embedded crack problems, all of the constant, linear, and quadratic ( N=0,1, or 2, respectively) pressure distributions yield very accurate results with this algorithm with 4 to 5 iterations. The edge crack problems, on the other hand, require much higher order polynomials distributions (N=5 to 6) to yield accurate solutions. For slant edge crack problems, the mode-I stress-intensity factors have better accuracy than the mode-II stress-intensity factors for the same convergence tolerance.

Krishnamurthy, T.↗