Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

A numerical study of fundamental shock noise mechanisms

The results of this thesis demonstrate that direct numerical simulation can predict sound generation in unsteady aerodynamic flows containing shock waves. Shock waves can be significant sources of sound in high speed jet flows, on helicopter blades, and in supersonic combustion inlets. Direct computation of sound permits the prediction of noise levels in the preliminary design stage and can be used as a tool to focus experimental studies, thereby reducing cost and increasing the probability of a successfully quiet product in less time. This thesis reveals and investigates two mechanisms fundamental to sound generation by shocked flows: shock motion and shock deformation. Shock motion is modeled by the interaction of a sound wave with a shock. During the interaction, the shock wave begins to move and the sound pressure is amplified as the wave passes through the shock. The numerical approach presented in this thesis is validated by the comparison of results obtained in a quasi-one dimensional simulation with linear theory. Analysis of the perturbation energy demonstrated for the first time that acoustic energy is generated by the interaction. Shock deformation is investigated by the numerical simulation of a ring vortex interacting with a shock. This interaction models the passage of turbulent structures through the shock wave. The simulation demonstrates that both acoustic waves and contact surfaces are generated downstream during the interaction. Analysis demonstrates that the acoustic wave spreads cylindrically, that the sound intensity is highly directional, and that the sound pressure level increases significantly with increasing shock strength. The effect of shock strength on sound pressure level is consistent with experimental observations of shock noise, indicating that the interaction of a ring vortex with a shock wave correctly models a dominant mechanism of shock noise generation.

Meadows, Kristine R.↗
NASA NTRS1994

Electron transport and energy degradation in the ionosphere: Evaluation of the numerical solution, comparison with laboratory experiments and auroral observations

Auroral electron transport calculations are a critical part of auroral models. We evaluate a numerical solution to the transport and energy degradation problem. The numerical solution is verified by reproducing simplified problems to which analytic solutions exist, internal self-consistency tests, comparison with laboratory experiments of electron beams penetrating a collision chamber, and by comparison with auroral observations, particularly the emission ratio of the N2 second positive to N2(+) first negative emissions. Our numerical solutions agree with range measurements in collision chambers. The calculated N(2)2P to N2(+)1N emission ratio is independent of the spectral characteristics of the incident electrons, and agrees with the value observed in aurora. Using different sets of energy loss cross sections and different functions to describe the energy distribution of secondary electrons that emerge from ionization collisions, we discuss the uncertainties of the solutions to the electron transport equation resulting from the uncertainties of these input parameters.

Lummerzheim, D.↗
NASA NTRS1992

A numerical study of three-dimensional surface tension driven convection with fre surface deformation

The steady three-dimensional thermocapillary motion with a deformable free surface is studied numerically in both normal and zero gravity environments. Flow configurations consist of a square cavity heated from the side. In the analysis, the free surface is allowed to deform and the grid distribution is adapted to the surface deformation. The divergence-free condition is satisfied by using a dual time-stepping approach in the numerical scheme. Convective flux derivatives are evaluated using a third-order accurate upwind-biased flux-split differencing technique. The numerical solutions at the midplane of the square cavity are compared with the results from two-dimensional calculations. In addition, numerial results for cases under zero and normal gravity conditions are compared. Significantly different flow structures and surface deformation have been observed. The comparison of calculated results will be compared with experimental data in the updated version of this paper.

Hsieh, Kwang-Chung↗
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 NTRS1995

Spurious Numerical Solutions Of Differential Equations

Paper presents detailed study of spurious steady-state numerical solutions of differential equations that contain nonlinear source terms. Main objectives of this study are (1) to investigate how well numerical steady-state solutions of model nonlinear reaction/convection boundary-value problem mimic true steady-state solutions and (2) to relate findings of this investigation to implications for interpretation of numerical results from computational-fluid-dynamics algorithms and computer codes used to simulate reacting flows.

Lafon, A.↗
NASA NTRS1996

Numerical Algorithms for Acoustic Integrals - The Devil is in the Details

The accurate prediction of the aeroacoustic field generated by aerospace vehicles or nonaerospace machinery is necessary for designers to control and reduce source noise. Powerful computational aeroacoustic methods, based on various acoustic analogies (primarily the Lighthill acoustic analogy) and Kirchhoff methods, have been developed for prediction of noise from complicated sources, such as rotating blades. Both methods ultimately predict the noise through a numerical evaluation of an integral formulation. In this paper, we consider three generic acoustic formulations and several numerical algorithms that have been used to compute the solutions to these formulations. Algorithms for retarded-time formulations are the most efficient and robust, but they are difficult to implement for supersonic-source motion. Collapsing-sphere and emission-surface formulations are good alternatives when supersonic-source motion is present, but the numerical implementations of these formulations are more computationally demanding. New algorithms - which utilize solution adaptation to provide a specified error level - are needed.

Brentner, Kenneth S.↗
NASA NTRS1996

A wavelet-optimized, very high order adaptive grid and order numerical method

Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.

Jameson, Leland↗
NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗
NASA NTRS1990

Direct Numerical Simulations of Turbulent Convection: Variable Gravity and Uniform Rotation - Part 1

Turbulent convection may have played a major role in determining the structure and evolution of the primordial solar nebula, but current, incomplete models of convection and turbulence give very different results and remain largely untested in the absence of detailed astronomical observations. Numerical simulations provide an 'experimental' database for comparison with these models, and, to this end, direct numerical simulations of turbulent convection were performed with modifications intended to mimic some of the unique physical features of thin accretion disks, such as the primordial solar nebula: internal heating, a gravitational acceleration that is linearly proportional to the distance from midplane of the nebula, and rapid rotation. Peclet numbers in the simulations are comparable to those in solar nebula models; Rossby numbers in the simulations are an order of magnitude larger than those in solar nebula models because of the unrealistically high Prandtl and low, Reynolds numbers required to resolve all scales of the convective flow. We find that, despite the loss of buoyancy at midplane, turbulent motions easily penetrate the midplane region with little loss of intensity, providing efficient transport of heat and turbulent kinetic energy throughout the interior. A simple mixing length model modified to include rotation is found to give convective heat fluxes for the interior flow in rough agreement with the numerical simulations. We discuss these preliminary results with regard to assumptions about heating distributions and convective heat fluxes made in standard solar nebula models. More definitive comparisons with solar nebula modelling will become possible when more realistic effects of shear, density, stratification, and compression are included.

Cabot, W.↗
NASA NTRS1997

Numerical modeling of HgCdTe solidification: Effects of phase diagram, double-diffusion convection and microgravity level

A numerical model of HgCdTe solidification was implemented using finite the element code FIDAP. Model verification was done using both experimental data and numerical test problems. The model was used to evaluate possible effects of double-diffusion convection in molten material, and microgravity level on concentration distribution in the solidified HgCdTe. Particular attention was paid to incorporation of HgCdTe phase diagram. It was found, that below a critical microgravity amplitude, the maximum convective velocity in the melt appears virtually independent on the microgravity vector orientation. Good agreement between predicted interface shape and an interface obtained experimentally by quenching was achieved. The results of numerical modeling are presented in the form of video film.

Bune, Andris V.↗
NASA NTRS1998

Microgravity Experiments and Numerical Modeling of Rotating Buoyant Convection in a Spherical Shell with Latitudinal Thermal Gradients

The results of numerical model calculations are compared with space laboratory experiments for rotating, thermally driven flow in a hemispherical shell, and various flow regimes are described. A spherically symmetric body force, analogous to gravity, is imposed in the radial direction through the use of a dielectric fluid and an electrostatic potential difference across the gap. The spherical boundaries are maintained at constant temperature profiles (with the pole and the inner sphere being the warmer surfaces), and the equatorial wall is an insulator. Typical parameter combinations result in highly nonlinear, but laminar, flow. For weak enough buoyant forcing, the flow is axisymmetric. It consists of a single meridional cell, rising in warm latitudes (in this case, near the pole) and sinking in cool latitudes, with prograde flow in the equatorial region and near the inner hemisphere for other latitudes, and retrograde motion elsewhere. For fast rotation, the first transition due to the instability of this simple flow is also axisymmetric, consisting of rings of convection in the polar region, which propagate poleward. The first nonaxisymmetric convection occurs at Rayleigh numbers which increase with rotation rate. The form of the convection near the transition also depends upon the rate of rotation. Selected flow patterns near the transition as well as those beyond it are studied numerically. For those cases where there exist laboratory experiments with which to compare, the numerical and experimental results agree very well.

Leslie, Fred W.↗
NASA NTRS1997

On the Universality of the Kolmogorov Constant in Numerical Simulations of Turbulence

Motivated by a recent survey of experimental data, we examine data on the Kolmogorov spectrum constant in numerical simulations of isotropic turbulence, using results both from previous studies and from new direct numerical simulations over a range of Reynolds numbers (up to 240 on the Taylor scale) at grid resolutions up to 512(exp 3). It is noted that in addition to k(exp -5/3) scaling, identification of a true inertial range requires spectral isotropy in the same wavenumber range. We found that a plateau in the compensated three-dimensional energy spectrum at k(eta) approx. = 0.1 - -0.2, commonly used to infer the Kolmogorov constant from the compensated three-dimensional energy spectrum, actually does not represent proper inertial range behavior. Rather, a proper, if still approximate, inertial range emerges at k(eta) approx. = 0.02 - 0.05 when R(sub lambda) increases beyond 140. The new simulations indicate proportionality constants C(sub 1) and C in the one- and three-dimensional energy spectra respectively about 0.60 and 1.62. If the turbulence were perfectly isotropic then use of isotropy relations in wavenumber space (C(sub 1) = 18/55 C) would imply that C(sub 1) approx. = 0.53 for C = 1.62, in excellent agreement with experiments. However the one- and three-dimensional estimates are not fully consistent, because of departures (due to numerical and statistical limitations) from isotropy of the computed spectra at low wavenumbers. The inertial scaling of structure functions in physical space is briefly addressed. Since DNS is still restricted to moderate Reynolds numbers, an accurate evaluation of the Kolmogorov constant is very difficult. We focus on providing new insights on the interpretation of Kolmogorov 1941 similarity in the DNS literature and do not consider issues pertaining to the refined similarity hypotheses of Kolmogorov (K62).

Yeung, P. K.↗
NASA NTRS1998

A Quantitative Investigation of Entrainment and Detrainment in Numerically Simulated Convective Clouds: Model Development - Pt. 1

A method is developed which uses numerical tracers to make accurate diagnoses of entraimnent and detrainment rates and of the properties of the entrained and detrained air in numerically simulated clouds. The numerical advection scheme is modified to make it nondispersive, as required by the use of the tracers. Tests of the new method are made, and an appropriate definition of clouds is selected. Distributions of mixing fractions in the model consistently show maximums at the end points, for nearly undilute environmental air or nearly undilute cloud air, with a uniform distribution between. The cumulonimbus clouds simulated here entrain air that had been substantially changed by the clouds, and detrained air that is not necessarily representative of the cloud air at the same level.

Cohen, Charles↗
NASA NTRS1999

Numerical Study of Buoyancy and Differential Diffusion Effects on the Structure and Dynamics of Triple Flames

Triple flames arise in a number of practical configurations where fuel and oxidizer are partially premixed, such as in the base of a lifted jet flame. Past experimental studies, theoretical analyses, and numerical modeling of triple flames suggested the potential role of triple flames in stabilizing turbulent flames and in promoting flame propagation. From recent numerical simulations of laminar triple flames, a strong influence of differential diffusion among species and heat on the triple flame structure has been gradually appreciated. This paper reports preliminary numerical results on the influence of gravity and differential diffusion effects on the structure and dynamics of triple flames with a one-step global irreversible chemistry model.

Chen, J. -Y.↗
NASA NTRS1999

Implicit and Multigrid Method for Ideal Multigrid Convergence: Direct Numerical Simulation of Separated Flow Around NACA 0012 Airfoil

Numerical investigation of flow separation over a NACA 0012 airfoil at large angles of attack has been carried out. The numerical calculation is performed by solving the full Navier-Stokes equations in generalized curvilinear coordinates. The second-order LU-SGS implicit scheme is applied for time integration. This scheme requires no tridiagonal inversion and is capable of being completely vectorized, provided the corresponding Jacobian matrices are properly selected. A fourth-order centered compact scheme is used for spatial derivatives. In order to reduce numerical oscillation, a sixth-order implicit filter is employed. Non-reflecting boundary conditions are imposed at the far-field and outlet boundaries to avoid possible non-physical wave reflection. Complex flow separation and vortex shedding phenomenon have been observed and discussed.

Liu, Chao-Qun↗
NASA NTRS2000

Numerical Simulation of Aircraft Trailing Vortices

The increase in air traffic is currently outpacing the development of new airport runways. This is leading to greater air traffic congestion, resulting in costly delays and cancellations. The National Aeronautics and Space Administration (NASA) under its Terminal Area Productivity (TAP) program is investigating new technologies that will allow increased airport capacity while maintaining the present standards for safety. As an element of this program, the Aircraft Vortex Spacing System (AVOSS) is being demonstrated in July 2000, at Dallas Ft-Worth Airport. This system allows reduced aircraft separations, thus increasing the arrival and departure rates, while insuring that wake vortices from a leading aircraft do not endanger trailing aircraft. The system uses predictions or wake vortex position and strength based on input from the current weather state. This prediction is accomplished by a semi-empirical model developed from theory, field observations, and relationships derived from numerical wake vortex simulations. Numerical experiments with a Large Eddy Simulation (LES) model are being conducted in order to provide guidance for the enhancement of these prediction algorithms. The LES Simulations of wake vortices are carried out with NASA's Terminal Area Simulation System (TASS). Previous wake vortex investigations with TASS are described. The primary objective of these numerical studies has been to quantify vortex transport and decay in relation to atmospheric variables. This paper summarizes many of the previous investigations with the TASS model and presents some new results regarding the onset of wake vortex decay.

Proctor, Fred H.↗
NASA NTRS1997

A Numerical Simulation of a Normal Sonic Jet into a Hypersonic Cross-Flow

This study involves numerical modeling of a normal sonic jet injection into a hypersonic cross-flow. The numerical code used for simulation is GASP (General Aerodynamic Simulation Program.) First the numerical predictions are compared with well established solutions for compressible laminar flow. Then comparisons are made with non-injection test case measurements of surface pressure distributions. Good agreement with the measurements is observed. Currently comparisons are underway with the injection case. All the experimental data were generated at the Southampton University Light Piston Isentropic Compression Tube.

Jeffries, Damon K.↗
NASA NTRS2000

Numerical Simulations of Homogeneous Turbulence Using Lagrangian-Averaged Navier-Stokes Equations

The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Mohseni, Kamran↗