SEARCH · Engineering Papers

Results for “unbounded domain”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

50 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1997

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗

NASA NTRS2003

Computational Methods for Jet Noise Simulation

The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Goodrich, John W.↗

NASA NTRS1989

Interaction dynamics of multiple mobile robots with simple navigation strategies

The global dynamic behavior of multiple interacting autonomous mobile robots with simple navigation strategies is studied. Here, the effective spatial domain of each robot is taken to be a closed ball about its mass center. It is assumed that each robot has a specified cone of visibility such that interaction with other robots takes place only when they enter its visibility cone. Based on a particle model for the robots, various simple homing and collision-avoidance navigation strategies are derived. Then, an analysis of the dynamical behavior of the interacting robots in unbounded spatial domains is made. The article concludes with the results of computer simulations studies of two or more interacting robots.

Wang, P. K. C.↗

NASA NTRS1989

Blow-up of unsteady two-dimensional Euler and Navier-Stokes solutions having stagnation-point form

The time-dependent form of the classic, two-dimensional stagnation-point solution of the Navier-Stokes equations is considered. If the viscosity is zero, a class of solutions of the initial-value problem can be found in closed form using Lagrangian coordinates. These solutions exhibit singular behavior in finite time, because of the infinite domain and unbounded initial vorticity. Thus, the blow-up found by Stuart in three dimensions using the stagnation-point form, also occurs in two. The singularity vanishes under a discrete, finite-dimensional 'point vortex' approximation, but is recovered as the number of vortices tends to infinity. We find that a small positive viscosity does not arrest the breakdown, but does strongly alter its form. Similar results are summarized for certain Boussinesq stratified flows.

Childress, S.↗

NASA NTRS1997

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗

NASA NTRS1994

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗

NASA NTRS1977

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗

NASA NTRS1995

Numerical Stability In Hyperbolic Boundary-Value Problems

Technical memorandum discusses stability of numerical solutions involving semidiscrete approximations to hyperbolic partial differential equations in initial-and-boundary-value problems. Topic of practical significance because hyperbolic partial differential equations arise in mathematical modeling of waves and blasts. Solutions often needed over restricted regions of unbounded spaces. Outer boundaries artificial, introduced only to limit domains of numerical solutions. Conditions at such artificial boundaries cause numerical instabilities that degrade computed solutions.

Warming, Robert F.↗

NASA NTRS2004

Signal Processing Variables for Optimization of Flaw Detection in Composites Using Ultrasonic Guided Wave Scanning

This study analyzes the effect of signal processing variables on the ability of the ultrasonic guided wave scan method at NASA Glenn Research Center to distinguish various flaw conditions in ceramic matrix composites samples. In the ultrasonic guided wave scan method, several time- and frequency-domain parameters are calculated from the ultrasonic guided wave signal at each scan location to form images. The parameters include power spectral density, centroid mean time, total energy (zeroth moment), centroid frequency, and ultrasonic decay rate. A number of signal processing variables are available to the user when calculating these parameters. These signal processing variables include 1) the time portion of the time-domain waveform processed, 2) integration type for the properties requiring integrations, 3) bounded versus unbounded integrations, 4) power spectral density window type, 5) and the number of time segments chosen if using the short-time fourier transform to calculate ultrasonic decay rate. Flaw conditions examined included delamination, cracking, and density variation.

Roth, Don J.↗

NASA NTRS2007

A Fast Method of Deriving the Kirchhoff Formula for Moving Surfaces

The Kirchhoff formula for a moving surface is very useful in many wave propagation problems, particularly in the prediction of noise from rotating machinery. Several publications in the last two decades have presented derivations of the Kirchhoff formula for moving surfaces in both time and frequency domains. Here we present a method originally developed by Farassat and Myers in time domain that is both simple and direct. It is based on generalized function theory and the useful concept of imbedding the problem in the unbounded three-dimensional space. We derive an inhomogeneous wave equation with the source terms that involve Dirac delta functions with their supports on the moving data surface. This wave equation is then solved using the simple free space Green's function of the wave equation resulting in the Kirchhoff formula. The algebraic manipulations are minimal and simple. We do not need the Green's theorem in four dimensions and there is no ambiguity in the interpretation of any terms in the final formulas. Furthermore, this method also gives the simplest derivation of the classical Kirchhoff formula which has a fairly lengthy derivation in physics and applied mathematics books. The Farassat-Myers method can be used easily in frequency domain.

Farassat, F.↗

NASA NTRS1977

Three-dimensional ion sound turbulence

A fast cold large diameter electron beam is injected into a uniform background plasma. The injected beam current is balanced by a field-aligned return current which gives rise to an ion acoustic instability. The plasma is essentially unbounded and collisionless on the time scale of the pulsed experiment. The fluctuations are analyzed by probe techniques in real time and space as well as in the frequency and wavenumber domains. Strong density fluctuations are observed in the entire ion acoustic spectrum reaching peak amplitudes at the low frequency end. The propagation of phase coherent test waves in the current carrying plasma is investigated. For the highly turbulent regime strong damping is observed while at lower drift velocities direction is observed. Oblique propagation of two strong test waves at different frequencies indicate the role of nonlinear ion Landau damping in absorption of wave energy.

Stenzel, R. L.↗

NASA NTRS1981

Far field conditions for compressible flows

The numerical treatment of exterior regions requires a method to convert the problem to one in a bounded region. One method is to truncate the unbounded region at some finite, artificial surface. This creates a finite computational region at the expense of imposing boundary conditions at the artificial boundary. If the problem has wave-like solutions near infinity, then these boundary conditions must simulate the radiation of energy out of the computational domain and towards infinity. Incorrect specification of these radiation boundary conditions can cause spurious reflected waves to be generated at the artificial boundary. Radiation boundary conditions which have the following properties were constructed: to accurately simulate the radiation of energy out of the computational domain; to accomplish the simulation with an accuracy which improves as fast as feasible as the artificial surface is moved outward; and to accelerate convergence to the steady state (by minimizing spurious reflections).

Bayliss, A.↗

NASA NTRS2020

Cross-field Streaming Instabilities in the Lower Hybrid Frequency Near the Front Pole Regions of a Magnetically Shielded Hall Thruster

The mechanism that led to the erosion of the inner front pole cover in wear tests of the Hall-effect Rocket with Magnetic Shielding (HERMeS) at 300 V and 12.5 kW has been elusive. Several classical mechanisms have been considered yet none has explained conclusively the source of the discrepancy between wear measurements and simulations. In this article we consider a new mechanism. The electrostatic dispersion relation for an unbounded, homogeneous plasma in a fixed magnetic field is solved numerically in the near-plume region of this thruster. The plasma conditions have been obtained from Hall2De (r-z) multifluid-PIC simulations. We find that the modified two-stream instability can be excited in most regions of the domain, with the fastest growth approaching the fluid-limit value of ωLH/2 near the magnetic poles where ωpe/ωce~1. In this region the relative drift between magnetized electrons and beam ions perpendicular to the magnetic field can exceed 3uTi, while that between counter-streaming beam and cathode ions can exceed 4uTi. The latter suggests that the ion-ion cross-field lower hybrid instability may also be active here. Unlike the classical ion-acoustic instability, lower hybrid instabilities can heat ions to temperatures that are comparable to the electron temperature. Such anomalous heating would enhance the average energy of these ions and, in turn, the erosion of the nearby surfaces they sputter. An idealized model that scales the anomalous ion heating with ωLH has been incorporated in Hall2De and the simulations now yield significantly improved agreement with the measured wear rates and IVDFs, as reported in our companion paper.

Ortega, Alejandro Lopez↗

NASA NTRS2018

Far-Field Lorenz-Mie Scattering in an Absorbing Host Medium: Theoretical Formalism and FORTRAN Program

In this paper we make practical use of the recently developed first-principles approach to electromagnetic scattering by particles immersed in an unbounded absorbing host medium. Specifically, we introduce an actual computational tool for the calculation of pertinent far-field optical observables in the context of the classical Lorenzâ€"Mie theory. The paper summarizes the relevant theoretical formalism, explains various aspects of the corresponding numerical algorithm, specifies the input and output parameters of a FORTRAN program available at https://www.giss.nasa.gov/staff/mmishchenko/Lorenz-Mie.html, and tabulates benchmark results useful for testing purposes. This public-domain FORTRAN program enables one to solve the following two important problems: (i) simulate theoretically the reading of a remote well-collimated radiometer measuring electromagnetic scattering by an individual spherical particle or a small random group of spherical particles; and (ii) compute the single-scattering parameters that enter the vector radiative transfer equation derived directly from the Maxwell equations.

Far-field electromagnetic scattering; Absorbing ho↗