SEARCH · Engineering Papers
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.
Numerical simulations with periodic boundary conditions are widely used in cosmology. These have a multiply connected topology known as a three-torus. Such nontrivial topologies for the actual universe may have arisen in the Big Bang. A two-dimensional numerical model with a twisted topology, sometimes a Klein bottle, is shown as well as the fact that local properties of the model are not dependent on topology.
The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and non-zero inter-blade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions have been coded for the inlet and exit boundary conditions. Sample unsteady solutions have been performed for an oscillating cascade and compared to experimental data. Also, test cases were fun for a flat plate cascade to compare with an unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with non-zero inter-blade phase angles are in good agreement with experimental data and small-perturbation theory. The zero degree inter-blade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus non-reflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions. This code is expected to be used as a tool for reviewing simpler models that do not include the full nonlinear aerodynamics or as a final check for designs against flutter in turbomachinery.
The Burgers equation with a small viscosity term, initial and periodic boundary conditions is resolved using a spatial approximation constructed from an orthonormal basis of wavelets. The algorithm is directly derived from the notions of multiresolution analysis and tree algorithms. Before the numerical algorithm is described these notions are first recalled. The method uses extensively the localization properties of the wavelets in the physical and Fourier spaces. Moreover, the authors take advantage of the fact that the involved linear operators have constant coefficients. Finally, the algorithm can be considered as a time marching version of the tree algorithm. The most important point is that an adaptive version of the algorithm exists: it allows one to reduce in a significant way the number of degrees of freedom required for a good computation of the solution. Numerical results and description of the different elements of the algorithm are provided in combination with different mathematical comments on the method and some comparison with more classical numerical algorithms.
Atkins' (1959) basic theory of third sound in thick films of superfluid He-4 is extended here to include a heat source for exciting third sound. A one-dimensional model obeying periodic boundary conditions is developed which can represent a variety of experimentally important configurations. General solutions for the model equations are found which give a complete description of third-sound waves in terms of displacement of the film surface, superfluid velocity, and temperature change as functions of space and time. Algebraic expressions for quantities that occur in general formulas are found for special excitation signals. The structure of third sound is calculated for the first time and the results are found to be inconsistent with previous findings. A direct method for accurately measuring the vaporization coefficient is also presented.
Report describes computer-simulation of instability and transition to turbulence in plane Poisseuille flow in rotating channel. Motions of fluid governed by Navier-Stokes equations of incompressible flow, using periodic boundary conditions in streamwise and spanwise directions. Principal parameters affecting flow are Reynolds and rotation numbers.
The KIVA-2 code, originally developed to solve computational fluid dynamics problems in internal combustion engines, has been developed to solve rocket propulsion type flows. The objective of the work was to develop a code such that both liquid and solid particle motion could be simulated for arbitrary geometry and high speed as well as low speed reacting flows. Modification to the original code include: incorporating independently specific supersonic and subsonic inflows and outflows; symmetric as well as periodic boundary conditions; and the capability to use generalized single or multi-specie thermodynamic data and transport coefficients allowing the user to specify arbitrary wall temperature/heat flux distributions. This code has been shown to successfully solve rocket propulsion flows as well as flows with entrained particles for several different rocket nozzles.
High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth order central differences through Fast Fourier Transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large-scale features, such as the total circulation around the roll-up region, are adequately resolved.
The evolution of two neighboring nonsymmetric arcades with footpoint shearing profiles which can change in time is simulated in order to investigate the quasi-static evolution of a double arcade system. It is found that rapid forced reconnection can take place repeatedly, forming a series of plasmoids. A two-arcade system laterally confined by surrounding fluxes can develop instability with nonsymmetric footpoint shear, but the instability is significantly lower in comparison with the case in which the system is confined by a 'hard box' or periodic boundary conditions. The fast plasmoid formation process due to asymmetry switching is not due to the presence of simulation walls. The behavior of double arcades with asymmetry switching is described using equilibrium bifurcation diagrams with footpoint shear and degree of asymmetry, as two independent control parameters.
The program to obtain intergrid communications - Pegasus - was enhanced to make better use of computing resources. Periodic block tridiagonal and penta-diagonal diagonal routines in OVERFLOW were modified to use a better algorithm to speed up the calculation for grids with periodic boundary conditions. Several programs were added to collar grid tools and a user friendly shell script was developed to help users generate collar grids. User interface for HYPGEN was modified to cope with the changes in HYPGEN. ET/SRB attach hardware grids were added to the computational model for the space shuttle and is currently incorporated into the refined shuttle model jointly developed at Johnson Space Center and Ames Research Center. Flow simulation for the integrated space shuttle vehicle at flight Reynolds number was carried out and compared with flight data as well as the earlier simulation for wind tunnel Reynolds number.
Two fully-discrete finite-difference schemes for wave propagation problems are presented, a maximum-order scheme and an optimized (or spectral-like) scheme. Both combine a seven-point spatial operator and an explicit six-stage time-march method. The maximum-order operator is fifth-order in space and is sixth-order in time for a linear problem with periodic boundary conditions. The phase and amplitude errors of the schemes obtained using Fourier analysis are given and compared with a second-order and a fourth-order method. Numerical experiments are presented which demonstrate the usefulness of the schemes for a range of problems. For some problems, the optimized scheme leads to a reduction in global error compared to the maximum-order scheme with no additional computational expense.
The spectral method provides an efficient algorithm for solving the 3D incompressible Navier-Stokes equations in periodic boundaries. Most people, so far, have used vectorized machines, such as the CRAY-2, to implement fast Fourier transformations and time integrations in the spectral calculations. In this paper, new results are presented using the spectral calculations on the Connection Machine-2 with a parallel algorithm. The large memory of the Connection Machine-2 and the parallel algorithm allows, of the first time, to implement a 512-cubed mesh resolution for high Reynolds number flows. The computational speed of the present code is about 30 percent faster than the fastest CRAY-2 simulations with four processors. Parallel machines, such as the Connection Machine-2, will possibly provide new computational power for understanding the intermittency and cascade mechanism in fluid turbulence.
The extreme eigenvalues of the Chebyshev pseudospectral differentiation operator are O(N exp 2) where N is the number of grid points. As a result of this, the allowable time step in an explicit time marching algorithm is O(N exp -2) which, in many cases, is much below the time step dictated by the physics of the partial differential equation. A new set of interpolating points is introduced such that the eigenvalues of the differentiation operator are O(N) and the allowable time step is O(N exp -1). The properties of the new algorithm are similar to those of the Fourier method. The new algorithm also provides a highly accurate solution for non-periodic boundary value problems.
Studying the topology of large-scale structure as a way to better understand initial conditions has become more widespread in recent years. Studying topology of simulations (which have periodic boundary conditions) in redshift space produces results compatible with the real topological characteristics of the simulation. Thus we expect we can extract useful information from redshift surveys. However, with nonperiodic boundary conditions, the use of smoothing must result in the loss of information at survey boundaries. In this paper, we test different methods of smoothing samples with nonperiodic boundary conditions to see which most efficiently preserves the topological features of the real distribution. We find that a smoothing method which (unlike most previous published analysis) sums only over cells inside the survey volume produces the best results among the schemes tested.
Flexural stiffness properties of a textile composite beam are obtained from a finite-element model of the unit cell. Three linearly independent deformations, namely, pure extension, pure bending and pure shear, are applied to the unit cell. The top and bottom surfaces of the beam are assumed to be traction free. Periodic boundary conditions on the lateral boundaries of the unit cell are enforced by multi-point constraint elements. From the forces acting on the unit cell, the flexural stiffness coefficients of the composite beam are obtained. The difficulties in determining the transverse shear stiffness are discussed, and a modified approach is presented. The methods are first verified by applying them to isotropic and bimaterial beams for which the results are known, and then illustrated for a simple plain-weave textile composite.
Contains results of calculations made to determine the neutral oscillatory stability boundaries, period and time to damp of the oscillatory mode, and motions following disturbances. The calculations were made for Mach numbers of 0.75 and 2.3 at an attitude of 35,000 feet and for the landing condition at sea level.
The differentiation matrix for a Daubechies-based wavelet basis is constructed and superconvergence is proven. That is, it will be proven that under the assumption of periodic boundary conditions that the differentiation matrix is accurate of order 2M, even though the approximation subspace can represent exactly only polynomials up to degree M-1, where M is the number of vanishing moments of the associated wavelet. It is illustrated that Daubechies-based wavelet methods are equivalent to finite difference methods with grid refinement in regions of the domain where small-scale structure is present.
A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.
A plasma simulation has been developed to model interactions between inertial Alfven waves and double layers and to investigate their relative contributions at auroral particle acceleration. We use a novel one-dimensional particle-in-cell code, with periodic boundary conditions, to model the nonlinear excitation of current-driven weak double layers via free energy supplied by an inertial Alfven wave. Analysis of the simulation output shows that double layers are not the agent primarily responsible for electron acceleration. Rather, the inertial Alfven wave accelerates groups of electrons into a steepening beam as it encounters them. As the beam electrons reenter the main distribution, decelerated by anomalous resistive effects, they are replaced by electrons farther downstream. Hence, the particles do not free-stream over the length of the channel. Furthermore, this wave action persists even when the system is linearly stable to ion-acoustic modes, precluding the possibility that this behavior is brought about by the formation of ion-acoustic double layers.