Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 757 records · Page 42

A numerical parametric study of a scramjet inlet in a Mach 6 arc heated test facility

A numerical study of the NASA Langley Arc Heated Scramjet Test Facility (AHSTF) is reported. The facility test flow through the Mach 6 quasi-two-dimensional nozzle was computed, and the flow through the Langley parametric inlet was calculated using different regions of the computed facility nozzle exit flow as inflow, exploring in more detail the effect of boundary layer ingestion on inlet performance. The potential advantage of using a quasi-two-dimensional nozzle in the Langley AHSTF versus the square cross-sectioned nozzle currently used for Scramujet engine tests is explored in terms of resulting inlet performance.

Sekar, B.↗

Three-dimensional zonal grids about arbitrary shapes by Poisson's equation

A method for generating 3-D finite difference grids about or within arbitrary shapes is presented. The 3-D Poisson equations are solved numerically, with values for the inhomogeneous terms found automatically by the algorithm. Those inhomogeneous terms have the effect near boundaries of reducing cell skewness and imposing arbitrary cell height. The method allows the region of interest to be divided into zones (blocks), allowing the method to be applicable to almost any physical domain. A FORTRAN program called 3DGRAPE has been written to implement the algorithm. Lastly, a method for redistributing grid points along lines normal to boundaries will be described.

Sorenson, Reese L.↗

A numerical study of the effects of reverse sweep on a scramjet inlet performance

A comparative numerical study of performance parameters of a similar and an opposite sweep sidewall compression inlet is made. The focus of the study is the investigation of the impact of alternate backward-forward sweep on the compression sidewalls as opposed to back-ward sweep on all the sidewalls. Two equivalent scramjet inlet configurations are designed for this purpose. These inlets have the same wetted areas of compression and expansion and same height and width; but in one inlet all the compression surface are swept back (similar sweep inlet) whereas in the other inlet, alternate surfaces are swept backward and forward (opposite sweep inlet). The cowl closure in both cases begins at the start of the throat region. A three-dimensional Navier-Stokes code is used to calculate the flow through these inlets. Results of these calculations are used to compare the two designs for their performance and flow quality. Effects of boundary- layer ingestion on the performance and overall flow features are also investigated.

Kumar, Ajay↗

Numerical simulation of transverse fuel injection

A review of recent work at NASA Langley Research Center to compare the predictions of transverse fuel injector flow fields and mixing performance with experimental results is presented. Various cold (non-reactive) mixing studies were selected for code calibration which include the effects of boundary layer thickness and injection angle for sonic hydrogen injection into supersonic air. Angled injection of helium is also included. This study was performed using both the three-dimensional elliptic and the parabolized Navier-Stokes (PNS) versions of SPARK. Axial solution planes were passed from PNS to elliptic and elliptic to PNS in order to efficiently generate solutions. The PNS version is used both upstream and far downstream of the injector where the flow can be considered parabolic in nature. The comparisons are used to identify experimental deficiencies and computational procedures to improve agreement.

Mao, Marlon↗

A simple, analytical, axisymmetric microburst model for downdraft estimation

A simple analytical microburst model was developed for use in estimating vertical winds from horizontal wind measurements. It is an axisymmetric, steady state model that uses shaping functions to satisfy the mass continuity equation and simulate boundary layer effects. The model is defined through four model variables: the radius and altitude of the maximum horizontal wind, a shaping function variable, and a scale factor. The model closely agrees with a high fidelity analytical model and measured data, particularily in the radial direction and at lower altitudes. At higher altitudes, the model tends to overestimate the wind magnitude relative to the measured data.

Vicroy, Dan D.↗

Results from computational analysis of a mixed compression supersonic inlet

A numerical study was performed to simulate the critical flow through a supersonic inlet. This flowfield has many phenomena such as shock waves, strong viscous effects, turbulent boundary layer development, boundary layer separations and mass flow suction through the walls (bleed). The computational tools used in this study were two full Navier-Stokes (FNS) codes. The supersonic inlet that was analyzed in this study is the Variable Diameter Centerbody (VDC), inlet. This inlet is a candidate concept for the next generation supersonic transport. Application of the code to the inlet geometry involved effort in generating an efficient grid geometry and specifying boundary conditions, particularly in the bleed region and at the outflow boundary. Results for a critical inlet operation compare favorably to Method of Characteristics predictions and experimental data.

Saunders, J. D.↗

Results from computational analysis of a mixed compression supersonic inlet

A numerical study was performed to simulate the critical flow through a supersonic inlet. This flow field has many phenomena such as shock waves, strong viscous effects, turbulent boundary layer development, boundary layer separations, and mass flow suction through the walls, (bleed). The computational tools used were two full Navier-Stokes (FNS) codes. The supersonic inlet that was analyzed is the Variable Diameter Centerbody, (VDC), inlet. This inlet is a candidate concept for the next generation supersonic involved effort in generating an efficient grid geometry and specifying boundary conditions, particularly in the bleed region and at the outflow boundary. Results for a critical inlet operation compare favorably to Method of Characteristics predictions and experimental data.

Saunders, J. D.↗

Analytical study of nozzle performance for nuclear thermal rockets

Nuclear propulsion has been identified as one of the key technologies needed for human exploration of the Moon and Mars. The Nuclear Thermal Rocket (NTR) uses a nuclear reactor to heat hydrogen to a high temperature followed by expansion through a conventional convergent-divergent nozzle. A parametric study of NTR nozzles was performed using the Rocket Engine Design Expert System (REDES) at the NASA Lewis Research Center. The REDES used the JANNAF standard rigorous methodology to determine nozzle performance over a range of chamber temperatures, chamber pressures, thrust levels, and different nozzle configurations. A design condition was set by fixing the propulsion system exit radius at five meters and throat radius was varied to achieve a target thrust level. An adiabatic wall was assumed for the nozzle, and its length was assumed to be 80 percent of a 15 degree cone. The results conclude that although the performance of the NTR, based on infinite reaction rates, looks promising at low chamber pressures, finite rate chemical reactions will cause the actual performance to be considerably lower. Parameters which have a major influence on the delivered specific impulse value include the chamber temperature and the chamber pressures in the high thrust domain. Other parameters, such as 2-D and boundary layer effects, kinetic rates, and number of nozzles, affect the deliverable performance of an NTR nozzle to a lesser degree. For a single nozzle, maximum performance of 930 seconds and 1030 seconds occur at chamber temperatures of 2700 and 3100 K, respectively.

Davidian, Kenneth O.↗

The void spectrum in two-dimensional numerical simulations of gravitational clustering

An algorithm for deriving a spectrum of void sizes from two-dimensional high-resolution numerical simulations of gravitational clustering is tested, and it is verified that it produces the correct results where those results can be anticipated. The method is used to study the growth of voids as clustering proceeds. It is found that the most stable indicator of the characteristic void 'size' in the simulations is the mean fractional area covered by voids of diameter d, in a density field smoothed at its correlation length. Very accurate scaling behavior is found in power-law numerical models as they evolve. Eventually, this scaling breaks down as the nonlinearity reaches larger scales. It is shown that this breakdown is a manifestation of the undesirable effect of boundary conditions on simulations, even with the very large dynamic range possible here. A simple criterion is suggested for deciding when simulations with modest large-scale power may systematically underestimate the frequency of larger voids.

Kauffmann, Guinevere↗

Finite element developments for two dimensional multiple-interface phase change problems

Finite element developments for multiple phase change problems in two-dimensional models are presented for the first time. The enthalpy method is used to simulate latent heat release in conjunction with fixed grid techniques. An unconditionally stable implicit method is used for the time integration. The effects of boundary conditions and the different phase regions on the multiple phase front developments are examined for numerous examples. Discussions and conclusions are appropriately addressed.

Ouyang, Tianhong↗

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

Experimental verification of the asymtotic modal analysis method as applied to a rectangular acoustic cavity excited by structural vibration

An experiment was performed on a rigid wall rectangular acoustic cavity driven by a flexible plate mounted in a quarter of one end wall and excited by white noise. The experiment was designed so that the assumptions of Asymptotic Modal Analysis (AMA) were satisfied for certain bandwidths and center frequencies. Measurements of sound pressure levels at points along the boundaries and incrementally into tbe interior were taken. These were compared with the theoretical results predicted with AMA, and found to be in good agreement, particularly for moderate (1/3 octave) bandwidths and sufficiently high center frequencies. Sound pressure level measurements were also taken well into the cavity interior at various points along the 5 totally rigid walls. The AMA theory, including boundary intensification effects, was shown to be accurate provided the assumption of large number of acoustic modes is satisfied, and variables such as power spectra of the wall acceleration, frequency, and damping are slowly varying in the frequency of bandwidth.

Peretti, L. F.↗

A nonlinear theory for spinning anisotropic beams using restrained warping functions

A geometrically nonlinear theory is developed for spinning anisotropic beams having arbitrary cross sections. An assumed displacement field is developed using the standard 3D kinematics relations to describe the global beam behavior supplemented with an additional field that represents the local deformation within the cross section and warping out of the cross section plane. It is assumed that the magnitude of this additional field is directly proportional to the local stress resultants. In order to take into account the effects of boundary conditions, a restraining function is introduced. This function plays the role of reducing the amount of free warping deformation throughout the field due to the restraint of the cross section(s) at the end(s) of the beam, e.g., in the case of a cantilever beam. Using a developed ordering scheme, the nonlinear strains are calculated to the third order. The FEM is developed using the weak form variational formulation. Preliminary interesting numerical results have been obtained that indicate the role of the restraining function in the case of a cantilever beam with circular cross section. These results are for the cases of a tip displacement (static) and free vibration studies for both isotropic and anisotropic materials with varied fiber orientations.

Ie, C. A.↗

Transition effects on compressible dynamic stall of transiently pitching airfoils

Experimental results and analysis of the effects of boundary layer tripping on dynamic stall of a transiently pitching airfoil are presented. At low Mach numbers, the tripped airfoil exhibits qualitative similarity with the behavior of the untripped airfoil. However, the local supersonic flow at Mach numbers greater than 0.3 is significantly modified by the trip leading to vastly different shock/boundary layer interactions, dynamic stall onset and vortex formation angles. The formation of the laminar separation bubble is found to have a favorable influence in delaying dynamic stall on the untripped airfoil flow. In both Mach number regimes, the tripped flow actually stalls at slightly lower angles of attack. Further experimentation with three trips on an oscillating airfoil showed that the dynamic stall process is very sensitive to the state of the turbulence in the boundary layer. This sensitivity points to a need for finer turbulence modeling techniques for use in dynamic stall computations.

Wilder, M. C.↗