Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Three-dimensional, transonic rotor flow field reconstructed from holographic interferogram data

Holographic interferometry and computer-assisted tomography (CAT) are used to determine the transonic flow field of a model rotor blade in hover. A pulsed ruby laser records 40 interferograms with a 61 cm-diam view field near the model rotor-blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting fringe-order functions, the data are transferred to a CAT code. The CAT code then calculates pressure coefficients in several planes above the blade surface. The values from the holography-CAT method compare favorably with previously obtained numerical computations and laser velocimeter measurements at most locations near the blade tip. The results demonstrate the technique's potential for three-dimensional transonic rotor flow studies.

Kittleson, J. K.↗
NASA NTRS1985

Reconstruction of a 3-dimensional transonic rotor flow field from holographic interferogram data

Holographic interferometry and computer-assisted tomography (CAT) are used to determine the transonic velocity field of a model rotor blade in hover. A pulsed ruby laser recorded 40 interferograms with a 2-ft-diam view field near the model rotor-blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting fringe-order functions, the data are transferred to a CAT code. The CAT code then calculates the perturbation velocity in seeral planes above the blade surface. The values from the holography-CAT method compare favorably with previously obtained numerical computations in most locations near the blade tip. The results demonstrate the technique's potential for three-dimensional transonic rotor flow studies.

Yu, Y. H.↗
NASA NTRS1985

Reconstruction of a three-dimensional, transonic rotor flow field from holographic interferogram data

Holographic interferometry and computer-assisted tomography (CAT) are used to determine the transonic velocity field of a model rotor blade in hover. A pulsed ruby laser recorded 40 interferograms with a 2-ft-diam view field near the model rotor-blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting fringe-order functions, the data are transferred to a CAT code. The CAT code then calculates the perturbation velocity in several planes above the blade surface. The value from the holography-CAT method compare favorably with previously obtained numerical computations in most locations near the blade tip. The results demonstrate the technique's potential for three-dimensional transonic rotor flow studies.

Yu, Y. H.↗
NASA NTRS1993

Normal modes of the world's oceans: A numerical investigation using Proudman functions

The numerical modeling of the normal modes of the global oceans is addressed. The results of such modeling could be expected to serve as a guide in the analysis of observations and measurements intended to detect these modes. The numerical computation of normal modes of the global oceans is a field in which several investigations have obtained results during the past 15 years. The results seem to be model-dependent to an unsatisfactory extent. Some modeling areas, such as higher resolution of the bathymetry, inclusion of self-attraction and loading, the role of the Arctic Ocean, and systematic testing by means of diagnostic models are addressed. The results show that the present state of the art is such that a final solution to the normal mode problem still lies in the future. The numerical experiments show where some of the difficulties are and give some insight as to how to proceed in the future.

Sanchez, Braulio V.↗
NASA NTRS1987

Transonic airfoil computation using the integral equation with and without embedded Euler domains

Two transonic computational schemes which are based on the Integral Equation Formulation of the full potential equation were presented. The first scheme is a Shock Capturing-Shock Fitting (SCSF) scheme which uses the full potential equation throughout with the exception of the shock wave where the Rankine-Hugoniot relations are used to cross and fit the shock. The second scheme is an Integral Equation with Embedded Euler (IEEE) scheme which uses the full potential equation with an embedded region where the Euler equations are used. The two schemes are applied to several transonic airfoil flows and the results were compared with numerous computational results and experimental domains with fine grids. The SCSF-scheme is restricted to flows with weak shock, while the IEEE-scheme can handle strong shocks. Currently, the IEEE scheme is applied to other transonic flows with strong shocks as well as to unsteady pitching oscillations.

Kandil, Osama A.↗
NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗
NASA NTRS1987

Computation of turbulent flows using an extended k-epsilon turbulence closure model

An extended kappa-epsilon turbulence model is proposed and tested with successful results. An improved transport equation for the rate of dissipation of the turbulent kinetic energy, epsilon, is proposed. The proposed model gives more effective response to the energy production rate than does the standard kappa-epsilon turbulence model. An extra time scale of the production range is included in the dissipation rate equation. This enables the present model to perform equally well for several turbulent flows with different characteristics, e.g., plane and axisymmetric jets, turbulent boundary layer flow, turbulent flow over a backward-facing step, and a confined turbulent swirling flow. A second-order accurate finite difference boundary layer code and a nearly second-order accurate finite difference elliptic flow solver are used for the present numerical computations.

Chen, Y.-S.↗
NASA NTRS2006

Stability of the Wave Bearing on an Elastic Support

Numerical computation predicts that an elastic support can substantially improve the stability of the wave bearing if the dynamic stiffness and damping of this support are in a specific range of values. To experimentally validate this prediction, the housing of a gas bearing was mounted on elastic O-rings and the threshold of sub-synchronous whirl motion was experimentally observed when the bearing runs unloaded with a rotating speed up to 30,000 RPM. The O-ring system was also dynamically characterized by measuring its stiffness and damping at various frequencies up to 500 Hz. Good correlation exists between the experimental data and numerical prediction.

Dimofte, Florin↗
NASA NTRS1951

General method and thermodynamic tables for computation of equilibrium composition and temperature of chemical reactions

A rapidly convergent successive approximation process is described that simultaneously determines both composition and temperature resulting from a chemical reaction. This method is suitable for use with any set of reactants over the complete range of mixture ratios as long as the products of reaction are ideal gases. An approximate treatment of limited amounts of liquids and solids is also included. This method is particularly suited to problems having a large number of products of reaction and to problems that require determination of such properties as specific heat or velocity of sound of a dissociating mixture. The method presented is applicable to a wide variety of problems that include (1) combustion at constant pressure or volume; and (2) isentropic expansion to an assigned pressure, temperature, or Mach number. Tables of thermodynamic functions needed with this method are included for 42 substances for convenience in numerical computations.

Huff, Vearl N↗
NASA NTRS1993

On acoustic radiation from a vibrating panel

An experimental and numerical study of the response and radiation from an acoustically loaded aircraft panel is presented. In the experiment, the panel is excited by a normally incident, harmonic wave. Measurements of the panel response and the resulting transmitted pressure are made in both the near- and acoustic far-fields. The numerical computations model the experiment and in particular account for the full coupling between the panel and the surrounding three dimensional acoustic fluid. The results demonstrate that for a sufficiently high excitation level, the panel response becomes nonlinear. The nonlinearity is characterized by the appearance of harmonics and subharmonics in the power spectral densities of the panel motion and consequently in the resulting acoustic radiation. The primary characteristic of the far-field acoustic radiation is the increase in harmonic content relative to the fundamental and subharmonics with increasing distance from the panel. This is shown to be a result of linear and weakly nonlinear wave propagation effects. The numerical results show that the radiated far-field pressure is strongly dependent on the position of the measurement point with respect to the panel center. The experimental and numerical results are in good qualitative agreement.

Frendi, Abdelkader↗
NASA NTRS2018

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗
NASA NTRSDate not supplied

Zemax Modeling of a Self-aligned Focusing Schlieren System with Simulated Compressible Flow

A self-aligned focusing schlieren system is simulated in the ANSYS Zemax OpticStudio software environment. The model uses a custom scatter-defining dynamic link library on the retroreflective background in order to scatter rays in the direction of the incident ray rather than the specular ray. Simulated compressible flow refractive index fields are imported into Zemax by use of the grid gradient object and a workaround to the maximum point limit is introduced. The Ronchi ruling is built as a high resolution bitmap image in MATLAB and set in Zemax as a slide object allowing the bitmap image to either transmit or terminate rays at each line pair. The Zemax simulations are compared to computational numerical schlieren and experimental conventional, path-integrated schlieren results from a hypersonic stagnation point injection flow configuration and an experimental image of a shock wave from a laser-induced breakdown. In the stagnation point injection simulations, the Zemax schlieren images capture complex flow features such as the Mach disk, shear layers, contact surface, and bow shock well, but contains long, curved features (striations) between the model and the bow shock that are not present in the synthetic schlieren image from the CFD or the experimental image. The image from Zemax of the spark-induced shock wave captures the shape and curvature of the shock but appears blurry due to grid discretization in the computational data.

Nicholas A Mejia↗
NASA NTRS2006

Beam Heating of Samples: Modeling and Verification

Energy absorbed from the X-ray beam by the sample requires cooling by forced convection (i.e. cryostream) to minimize temperature increase and the damage caused to the sample by the X-ray heating. In this presentation we will first review the current theoretical models and recent studies in the literature, which predict the sample temperature rise for a given set of beam parameters. It should be noted that a common weakness of these previous studies is that none of them provide actual experimental confirmation. This situation is now remedied in our investigation where the problem of x-ray sample heating is taken up once more. We have theoretically investigated, and at the same time, in addition to the numerical computations, performed experiments to validate the predictions. We have modeled, analyzed and experimentally tested the temperature rise of a 1 mm diameter glass sphere (sample surrogate) exposed to an intense synchrotron X-ray beam, while it is being cooled in a uniform flow of nitrogen gas. The heat transfer, including external convection and internal heat conduction was theoretically modeled using CFD to predict the temperature variation in the sphere during cooling and while it was subjected to an undulator (ID sector 19) X-ray beam at the APS. The surface temperature of the sphere during the X-ray beam heating was measured using the infrared camera measurement technique described in a previous talk. The temperatures from the numerical predictions and experimental measurements are compared and discussed. Additional results are reported for the two different sphere sizes and for two different supporting pin orientations.

Kazmierczak, Michael↗
NASA NTRS1976

A complete model of turbulence

A set of constitutive equations suitable for a priori computation of turbulent shear flows has been developed. Since no properties of a given turbulent flow need be known in advance in order to obtain a solution, the equations comprise a complete model of turbulence. Perturbation analysis shows that the model predicts a composite five-layer structure for an incompressible turbulent boundary layer, viz, a defect layer, a law-of-the-wall layer, a viscous sublayer, a near-surface roughness layer, and a viscous superlayer at the boundary-layer edge. Analysis of the defect layer demonstrates the key improvement of the model over its predecessor, the Saffman-Wilcox two-equation model of turbulence. Examination of model-predicted sublayer structure yields model-parameter boundary conditions appropriate for surfaces with roughness and mass injection. Results of numerical computations of compressible and incompressible equilibrium boundary layers show that, for such flows, the model is as accurate as mixing-length theory. Applications to transitional boundary layers and to nonequilibrium relaxation of a boundary layer passing from a rough to a smooth surface indicate that the model's applicability extends far beyond that of mixing-length theory's.

Wilcox, D. C.↗
NASA NTRS2001

Fluid Aspects of Solar Wind Disturbances Driven by Coronal Mass Ejections

Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.

Gosling, J. T.↗
NASA NTRS2018

Transition Delay via Vortex Generators in a Hypersonic Boundary Layer at Flight Conditions

The potential of realizable, stationary streaks undergoing non-modal growth to stabilize a hypersonic boundary-layer flow and, subsequently, delay the laminar-turbulent transition onset, is studied via numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach number of 5.30, freestream unit Reynolds number equal to 13.42 x 10(exp 6)/m, and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the test article. This paper investigates flow modifications induced by wall-mounted vortex generators (VGs), followed by an analysis of the modal instability of the perturbed, streaky boundary-layer flow. Results are presented both for a single array of VGs that was designed on the basis of optimal growth theory and for a VG configuration involving two separate arrays with opposite orientations that ware designed to provide staged control of flow instabilities while simultaneously reducing the amplification of streak instabilities resulting from the control devices. Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment, which did not include any control devices, correlated with an amplification factor of N = 14.7 for the planar Mack modes. If one assumes that the transition N -factor is not affected by the introduction of the VGs, then the control configurations based on a single array of VGs and two separate arrays would result in a transition delay of 17% and 40%, respectively. These findings suggest a passive flow control s to induce streaks that would delay transition in hypersonic boundary dominated by Mack-mode instabilities.

Incompressible boundary layer↗