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.

SourceResearch connectionFrom yearThrough year

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

At least 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

NASA NTRS1996

A numerical investigation of premixed combustion in wave rotors

Wave rotor cycles which utilize premixed combustion processes within the passages are examined numerically using a one-dimensional CFD-based simulation. Internal-combustion wave rotors are envisioned for use as pressure-gain combustors in gas turbine engines. The simulation methodology is described, including a presentation of the assumed governing equations for the flow and reaction in the channels, the numerical integration method used, and the modeling of external components such as recirculation ducts. A number of cycle simulations are then presented which illustrate both turbulent-deflagration and detonation modes of combustion. Estimates of performance and rotor wall temperatures for the various cycles are made, and the advantages and disadvantages of each are discussed.

Nalim, M. Razi↗

NASA NTRS1996

Theoretical and Numerical Investigation of Radiative Extinction of Diffusion Flames

The influence of soot radiation on diffusion flames was investigated using both analytical and numerical techniques. Soot generated in diffusion flames dominate the flame radiation over gaseous combustion products and can significantly lower the temperature of the flame. In low gravity situations there can be significant accumulation of soot and combustion products in the vicinity of the primary reaction zone owing to the absence of any convective buoyant flow. Such situations may result in substantial suppression of chemical activities in a flame, and the possibility of a radiative extinction may also be anticipated. The purpose of this work was to not only investigate the possibility of radiative extinction of a diffusion flame but also to qualitatively and quantitatively analyze the influence of soot radiation on a diffusion flame. In this study, first a hypothetical radiative loss profile of the form of a sech(sup 2) was assumed to influence a pure diffusion flame. It was observed that the reaction zone can, under certain circumstances, move through the radiative loss zone and locate itself on the fuel side of the loss zone contrary to our initial postulate. On increasing the intensity and/or width of the loss zone it was possible to extinguish the flame, and extinction plots were generated. In the presence of a convective flow, however, the movement of the temperature and reaction rate peaks indicated that the flame behavior is more complicated compared to a pure diffusional flame. A comprehensive model of soot formation, oxidation and radiation was used in a more involved analysis. The soot model of Syed, Stewart and Moss was used for soot nucleation and growth and the model of Nagle and Strickland-Constable was used for soot oxidation. The soot radiation was considered in the optically thin limit. An analysis of the flame structure revealed that the radiative loss term is countered both by the reaction term and the diffusion term. The essential balance for the soot volume fraction was found to be between the processes of soot convection and soot growth. Such a balance yielded to analytical treatment and the soot volume fraction could be expressed in the form of an integral. The integral was evaluated using two approximate methods and the results agreed very well with the numerical solutions for all cases examined.

Ray, Anjan↗

NASA NTRS1996

Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

Mickens, Ronald E.↗

NASA NTRS1996

Numerical Algorithms Based on Biorthogonal Wavelets

Wavelet bases are used to generate spaces of approximation for the resolution of bidimensional elliptic and parabolic problems. Under some specific hypotheses relating the properties of the wavelets to the order of the involved operators, it is shown that an approximate solution can be built. This approximation is then stable and converges towards the exact solution. It is designed such that fast algorithms involving biorthogonal multi resolution analyses can be used to resolve the corresponding numerical problems. Detailed algorithms are provided as well as the results of numerical tests on partial differential equations defined on the bidimensional torus.

Ponenti, Pj.↗

NASA NTRS1994

Numerical Simulation of the Instability of an Inviscid Liquid Jet in a Coflowing Compressible Airstream

The nonlinear interfacial instability of a liquid jet in a coflowing compressible airstream is studied numerically. A high-resolution scheme which has second-order accuracy in space and time is coupled with a Lagrangian marker particle algorithm to visualize the large-scale motion of the interfaces in compressible flow. A numerical algorithm based on an approximate equation of state of a compressible liquid is developed to allow this two-fluid system to be governed by the nonlinear unsteady Euler equations in conservative form. The initial growth of small disturbances given by the simulations agrees well with linear theory. The process of jet disruption in compressible flow is demonstrated to consist of the formation of liquid spikes, interweaving of the gas and liquid and stretching and detachment of the liquid main center core.

Li, H. -S.↗

NASA NTRS1997

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗

NASA NTRS1996

Numerical modelling of the circulation of superrotating atmospheres: Venus and Titan

Atmospheric superrotation is known to exist on Venus where the cloud cover, at 70 km above the surface, rotates about 60 times faster than the solid planet. Superrotation is strongly suspected on Titan on the basis of numerical studies and observations. The numerical simulation of superrotation is described. It is shown how the general circulation models, originally developed for terrestrial meteorology and climatology, can be used to infer the circulation in such atmospheres. The emphasis is on the case of Titan for which a specific general circulation model was developed.

Hourding, F.↗

NASA NTRS1997

Science Support for Space-Based Droplet Combustion: Drop Tower Experiments and Detailed Numerical Modeling

This program supports the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies. Experimental emphasis is on the study of simple alcohols (methanol, ethanol) and alkanes (n-heptane, n-decane) as fuels with time dependent measurements of drop size, flame-stand-off, liquid-phase composition, and finally, extinction. Experiments have included bench-scale studies at Princeton, studies in the 2.2 and 5.18 drop towers at NASA-LeRC, and both the Fiber Supported Droplet Combustion (FSDC-1, FSDC-2) and the free Droplet Combustion Experiment (DCE) studies aboard the shuttle. Test matrix and data interpretation are performed through spherically-symmetric, time-dependent numerical computations which embody detailed sub-models for physical and chemical processes. The computed burning rate, flame stand-off, and extinction diameter are compared with the respective measurements for each individual experiment. In particular, the data from FSDC-1 and subsequent space-based experiments provide the opportunity to compare all three types of data simultaneously with the computed parameters. Recent numerical efforts are extending the computational tools to consider time dependent, axisymmetric 2-dimensional reactive flow situations.

Marchese, Anthony J.↗

NASA NTRS1997

Numerical Optimization Using Computer Experiments

Engineering design optimization often gives rise to problems in which expensive objective functions are minimized by derivative-free methods. We propose a method for solving such problems that synthesizes ideas from the numerical optimization and computer experiment literatures. Our approach relies on kriging known function values to construct a sequence of surrogate models of the objective function that are used to guide a grid search for a minimizer. Results from numerical experiments on a standard test problem are presented.

Trosset, Michael W.↗

NASA NTRS1997

Numerical Model Studies of the Martian Mesoscale Circulations

The study objectives were to evaluate by numerical modeling various possible mesoscale circulation on Mars and related atmospheric boundary layer processes. The study was in collaboration with J. Tillman of the University of Washington (who supported the study observationally). Interaction has been made with J. Prusa of Iowa State University in numerical modeling investigation of dynamical effects of topographically-influenced flow. Modeling simulations included evaluations of surface physical characteristics on: (i) the Martian atmospheric boundary layer and (ii) their impact on thermally and dynamically forced mesoscale flows. Special model evaluations were made in support of selection of the Pathfinder landing sites. J. Tillman's finding of VL-2 inter-annual temperature difference was followed by model simulations attempting to point out the forcing for this feature. Publication of the results in the reviewed literature in pending upon completion of the manuscripts in preparation as indicated later.

Segal, Moti↗

NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗

NASA NTRS1994

A Numerical Comparison of Barrier and Modified Barrier Methods for Large-Scale Bound-Constrained Optimization

When a classical barrier method is applied to the solution of a nonlinear programming problem with inequality constraints, the Hessian matrix of the barrier function becomes increasingly ill-conditioned as the solution is approached. As a result, it may be desirable to consider alternative numerical algorithms. We compare the performance of two methods motivated by barrier functions. The first is a stabilized form of the classical barrier method, where a numerically stable approximation to the Newton direction is used when the barrier parameter is small. The second is a modified barrier method where a barrier function is applied to a shifted form of the problem, and the resulting barrier terms are scaled by estimates of the optimal Lagrange multipliers. The condition number of the Hessian matrix of the resulting modified barrier function remains bounded as the solution to the constrained optimization problem is approached. Both of these techniques can be used in the context of a truncated-Newton method, and hence can be applied to large problems, as well as on parallel computers. In this paper, both techniques are applied to problems with bound constraints and we compare their practical behavior.

Nash, Stephen G.↗

NASA NTRS1997

Numerical Computation of Sensitivities and the Adjoint Approach

We discuss the numerical computation of sensitivities via the adjoint approach in optimization problems governed by differential equations. We focus on the adjoint problem in its weak form. We show how one can avoid some of the problems with the adjoint approach, such as deriving suitable boundary conditions for the adjoint equation. We discuss the convergence of numerical approximations of the costate computed via the weak form of the adjoint problem and show the significance for the discrete adjoint problem.

Lewis, Robert Michael↗

NASA NTRS1998

Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft

This report compiles the various research activities conducted under the auspices of the NASA Grant NAG3-1026, "Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft" during the period of April 1989 to April 1994. The effort involved the development of multigrid based algorithms and computer programs for the calculation of the flow and temperature fields generated by Short Take-off and Vertical Landing (STOVL) aircraft, while hovering in ground proximity. Of particular importance has been the interaction of the exhaust jets with the head wind which gives rise to the hot gas ingestion process. The objective of new STOVL designs to reduce the temperature of the gases ingested into the engine. The present work describes a solution algorithm for the multi-dimensional elliptic partial-differential equations governing fluid flow and heat transfer in general curvilinear coordinates. The solution algorithm is based on the multigrid technique which obtains rapid convergence of the iterative numerical procedure for the discrete equations. Initial efforts were concerned with the solution of the Cartesian form of the equations. This algorithm was applied to a simulated STOVL configuration in rectangular coordinates. In the next phase of the work, a computer code for general curvilinear coordinates was constructed. This was applied to model STOVL geometries on curvilinear grids. The code was also validated in model problems. In all these efforts, the standard k-Epsilon model was used.

Vanka, S. P.↗

NASA NTRS1998

Numerical Simulation of the Generation of Axisymmetric Mode Jet Screech Tones

An imperfectly expanded supersonic jet, invariably, radiates both broadband noise and discrete frequency sound called screech tones. Screech tones are known to be generated by a feedback loop driven by the large scale instability waves of the jet flow. Inside the jet plume is a quasi-periodic shock cell structure. The interaction of the instability waves and the shock cell structure, as the former propagates through the latter, is responsible for the generation of the tones. Presently, there are formulas that can predict the tone frequency fairly accurately. However, there is no known way to predict the screech tone intensity. In this work, the screech phenomenon of an axisymmetric jet at low supersonic Mach number is reproduced by numerical simulation. The computed mean velocity profiles and the shock cell pressure distribution of the jet are found to be in good agreement with experimental measurements. The same is true with the simulated screech frequency. Calculated screech tone intensity and directivity at selected jet Mach number are reported in this paper. The present results demonstrate that numerical simulation using computational aeroacoustics methods offers not only a reliable way to determine the screech tone intensity and directivity but also an opportunity to study the physics and detailed mechanisms of the phenomenon by an entirely new approach.

Shen, Hao↗

NASA NTRS1960

The Numerical Calculation of Flow Past Conical Bodies Supporting Elliptic Conical Shock Waves at Finite Angles of Incidence

The inverse method, with the shock wave prescribed to be an elliptic cone at a finite angle of incidence, is applied to calculate numerically the supersonic perfect-gas flow past conical bodies not having axial symmetry. Two formulations of the problem are employed, one using a pair of stream functions and the other involving entropy and components of velocity. A number of solutions are presented, illustrating the numerical methods employed, and showing the effects of moderate variation of the initial parameters.

Briggs, Benjamin R.↗

NASA NTRS1998

Numerical Study of Rotating Turbulence with External Forcing

Direct numerical simulation at 256(exp 3) resolution have been carried out to study the response of isotropic turbulence to the concurrent effects of solid-body rotation and numerical forcing at the large scales. Because energy transfer to the smaller scales is weakened by rotation, energy input from forcing gradually builds up at the large scales, causing the overall kinetic energy to increase. At intermediate wavenumbers the energy spectrum undergoes a transition from a limited k(exp -5/3) inertial range to k(exp -2) scaling recently predicted in the literature. Although the Reynolds stress tensor remains approximately isotropic and three-components, evidence for anisotropy and quasi- two-dimensionality in length scales and spectra in different velocity components and directions is strong. The small scales are found to deviate from local isotropy, primarily as a result of anisotropic transfer to the high wavenumbers. To understand the spectral dynamics of this flow we study the detailed behavior of nonlinear triadic interactions in wavenumber space. Spectral transfer in the velocity component parallel to the axis of rotation is qualitatively similar to that in non-rotating turbulence; however the perpendicular component is characterized by a greatly suppressed energy cascade at high wavenumber and a local reverse transfer at the largest scales. The broader implications of this work are briefly addressed.

Yeung, P. K.↗

NASA NTRS1998

The Space-Time Conservation Element and Solution Element Method-A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: Numerical Simulation of Shock Waves and Contact Discontinuities - 2

Without resorting to special treatment for each individual test case, the 1D and 2D CE/SE shock-capturing schemes described previously (in Part I) are used to simulate flows involving phenomena such as shock waves, contact discontinuities, expansion waves and their interactions. Five 1D and six 2D problems are considered to examine the capability and robustness of these schemes. Despite their simple logical structures and low computational cost (for the 2D CE/SE shock-capturing scheme, the CPU time is about 2 micro-secs per mesh point per marching step on a Cray C90 machine), the numerical results, when compared with experimental data, exact solutions or numerical solutions by other methods, indicate that these schemes can accurately resolve shock and contact discontinuities consistently.

Wang, Xiao-Yen↗