SEARCH · Engineering Papers

Results for “Numerical techniques”

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,279 records · Page 71

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Numerical simulation of the flow about the F-18 HARV at high angle of attack

This research has been aimed at validating numerical methods for computing the flow about the complete F-18 HARV at alpha = 30 deg and alpha = 45 deg. At 30 deg angle of attack, the flow about the F-18 is dominated by the formation, and subsequent breakdown, of strong vortices over the wing leading-edge extensions (LEX). As the angle of attack is increased to alpha = 45 deg, the fuselage forebody of the F-18 contains significant laminar and transitional regions which are not present at alpha = 30 deg. Further, the flow over the LEX at alpha = 45 deg is dominated by an unsteady shedding in time, rather than strong coherent vortices. This complex physics, combined with the complex geometry of a full-aircraft configuration, provides a challenge for current computational fluid dynamics (CFD) techniques. The following sections present the numerical method and grid generation scheme that was used, a review of prior research done to numerically model the F-18 HARV, and a discussion of the current research. The current research is broken into three main topics; the effect of engine-inlet mass-flow rate on the F-18 vortex breakdown position, the results using a refined F-18 computational model to compute the flow at alpha = 30 deg and alpha = 45 deg, and research done using the simplified geometry of an ogive-cylinder configuration to investigate the physics of unsteady shear-layer shedding. The last section briefly summarizes the discussion.

Murman, Scott M.↗

NASA NTRS1988

Material parameter measurements at high temperatures

Alternate fixtures of techniques for the measurement of the constitutive material parameters at elevated temperatures are presented. The technique utilizes scattered field data from material coated cylinders between parallel plates or material coated hemispheres over a finite size groundplane. The data acquisition is centered around the HP 8510B Network Analyzer. The parameters are then found from a numerical search algorithm using the Newton-Ralphson technique with the measured and calculated fields from these canonical scatters. Numerical and experimental results are shown.

Dominek, A.↗

NASA NTRS1973

Use of Polarization Lidar for Investigation of Meteorological Formations

This paper presents the results of theoretical and experimental, investigations of depolarization characteristics of different meteorological formations. Experimental investigations are carried out with a monostatic lidar. The ruby laser radiation is polarized in a vertical plane. The radiation reflected is accepted by a lens system of 150 rom in diameter and a viewing angle of 12' and further it is divided by Wollaston prism into the components polarized orthogonally. In this case the principal plane of the prism is exposed parallel with the laser polarization plane. Investigations show the degree of radiation polarization, reflected from water clouds, to be changed within 1/0.7 (seldom up to 0.6) depending on their density. In most cases a signal reflected from the cloud leading edge is polarized completely. The time shift is observed between polarized and crosspolarized components of a Fignal, reflected from a cloud, depending on the density of a meteorological object. While penetrating into the cloud depth a degree of polarization decreases up to 0.8-0.7, and the character of this decrease is different for various types of clouds. For crystal clouds the shift between the components of the reflected signal is not observed and the magnitude of polarization degree amounts to 0.1/0.3 in comparison with water clouds. The polarization degree of radiation reflected by fog is not less than 0.6, and that in the rains of average intensity (about 5 mm/h) is always about 1. The authors have suggested an algorithm of numerical solution of nonstationary transfer equation in the vector form to forecast the influence of multiple scattering effects on polarization characteristics of the lidar light signal. The method of statistical simulation (Monte-Carlo technique) forms the basis of the algorithm. Numerical estimates obtained for a model of stratocumulus at lambda = 0.6943 microns under boundary conditions close to the conditions of natural experiment being discussed proved to be in a good agreement with the results of observation. Specifically, Fig. 1 shows the profiles of polarization (p) versus depth (L) of the following drop formations: fog (curve 1) with horizontal meteorological visibility of 4 km two stratocumulus at a height of 1100 m with the attenuation factors delta = 0.01m(exp -1) (curve 2) and delta = 0.05m(exp -1) (curve 3). Curve 3 shows the results of numerical estimates and the value of their statistical error.

Balin, Yu. S.↗

NASA NTRS1974

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗

NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 2: Solutions and numerical results

This paper discusses an incremental-iterative nonlinear solution technique for solving the nonlinear finite element equations of the superconducting state of a superconductor. The untreated equations are highly ill-conditioned and are impossible to solve within the typical 16-place double precision supplied by most computers. A combination of matrix scaling and mesh grading techniques is used to reduce the condition number of the tangent stiffness matrix and increase the accuracy of the current carrying boundary layer representation. Numerical results for a one-dimensional model of a time-independent superconductor treated by the Ginzburg-Landau model are presented and discussed. The computed solutions clearly display the Meissner effect of magnetic field expulsion from the central region of the superconductor. These results are compared to the physics of a low-viscosity fluid problem. From this analogy, a physical argument is advanced about the macroscopic behavior of superconductors.

Schuler, James J.↗

NASA NTRS1983

A survey of grid generation techniques in computational fluid dynamics

The use of numerically-generated body-conforming coordinate systems has become well-established in computational fluid dynamics, allowing general codes to be written which can treat boundaries of arbitrary shape. There is now a rapidly increasing interest in the dynamic coupling of the grid with the physical solution so that the grid lines continually move to achieve concentration in regions of strong variation of the physical solution and alignment with shocks, flame fronts, etc. Such adaptive coordinate systems have been developed using a variety of techniques for sensing and reacting to the developing regions of concentration and/or alignment. This paper reviews the techniques and applications of these dynamically adaptive systems and other systems.

Thompson, J. F.↗

NASA NTRS1991

Spectrum transformation for divergent iterations

Certain spectrum transformation techniques are described that can be used to transform a diverging iteration into a converging one. Two techniques are considered called spectrum scaling and spectrum enveloping and how to obtain the optimum values of the transformation parameters is discussed. Numerical examples are given to show how this technique can be used to transform diverging iterations into converging ones; this technique can also be used to accelerate the convergence of otherwise convergent iterations.

Gupta, Murli M.↗

NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗

NASA NTRS1973

Application of phase-change technique to thin sections with heating on both surfaces

A numerical digital computer program was developed to calculate the heat-transfer coefficients for both sides of a one-dimensional finite slab subject to the boundary conditions ascribed to the phase-change coating heat-transfer technique. In a typical tunnel test situation where a thin wing was exposed to heating on both sides, the data reduction procedures for a semi-infinite slab gave heat-transfer coefficients as much as 375 percent too high on the side with the lowest heating. The results from the one-dimensional finite-slab procedure are presented in the form of correction factors to the solution for a semi-infinite slab in terms of parameters normally used with the phase-change heat-transfer technique. These correlations are not restricted to slab thickness or thermophysical properties and are easily used to obtain accurate data on thin model sections.

Hunt, J. L.↗

NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗

NASA NTRS1990

Short-time Lyapunov exponent analysis

A new technique for analyzing complicated fluid flows in numerical simulations has been successfully tested. The analysis uses short time Lyapunov exponent contributions and the associated Lyapunov perturbation fields. A direct simulation of the Taylor-Couette flow just past the onset of chaos demonstrated that this new technique marks important times during the system evolution and identifies the important flow features at those times. This new technique will now be applied to a 'minimal' turbulent channel.

Vastano, J. A.↗

NASA NTRS1992

Processing infrared images of aircraft lapjoints

Techniques for processing IR images of aging aircraft lapjoint data are discussed. Attention is given to a technique for detecting disbonds in aircraft lapjoints which clearly delineates the disbonded region from the bonded regions. The technique is weak on unpainted aircraft skin surfaces, but can be overridden by using a self-adhering contact sheet. Neural network analysis on raw temperature data has been shown to be an effective tool for visualization of images. Numerical simulation results show the above processing technique to be an effective tool in delineating the disbonds.

Syed, Hazari↗

NASA NTRS1986

A multistage time-stepping scheme for the thin-layer Navier-Stokes equations

A finite-volume scheme for numerical integration of the Euler equations was extended to allow solution of the thin-layer Navier-Stokes equations in two and three dimensions. The extended algorithm, which is based on a class of four-stage Runge-Kutta time-stepping schemes, was made numerically efficient through the following convergence acceleration technique: (1) local time stepping, (2) enthalpy damping, and (3) residual smoothing. Also, the high degree of vectorization possible with the algorithm has yielded an efficient program for vector processors. The scheme was evaluated by solving laminar and turbulent flows. Numerical results have compared well with either theoretical or other numerical solutions and/or experimental data.

Swanson, R. C., Jr.↗

NASA NTRS2019

Thermal Property Estimation of Fibrous Insulation: Heat Transfer Modeling and the Continuous Genetic Algorithm

Thermal properties of high-temperature fibrous insulation materials were estimated from transient thermal tests using an inverse heat transfer technique. Transient temperature data from an experimental set up was collected and a simple, one-dimensional numerical model was constructed to replicate the temperatures within the test assembly. The Continuous Genetic Algorithm optimization technique in conjunction with a numerical thermal model and transient test data was used to estimate coefficients of a functional representation of the thermal property. The thermal properties, i.e., thermal conductivity and specific heat, of an alumina insulation felt were estimated over the temperature range of 300 K to 1700 K at various constant static pressures in nitrogen gas and compared with published data. The resulting thermal property estimates were within 10% of published values over the entire temperature range at various pressures. The methodology, application, and results are presented.

Frye, Elora↗

NASA NTRS1983

Asymptotic analysis of numerical wave propagation in finite difference equations

An asymptotic technique is developed for analyzing the propagation and dissipation of wave-like solutions to finite difference equations. It is shown that for each fixed complex frequency there are usually several wave solutions with different wavenumbers and the slowly varying amplitude of each satisfies an asymptotic amplitude equation which includes the effects of smoothly varying coefficients in the finite difference equations. The local group velocity appears in this equation as the velocity of convection of the amplitude. Asymptotic boundary conditions coupling the amplitudes of the different wave solutions are also derived. A wavepacket theory is developed which predicts the motion, and interaction at boundaries, of wavepackets, wave-like disturbances of finite length. Comparison with numerical experiments demonstrates the success and limitations of the theory. Finally an asymptotic global stability analysis is developed.

Giles, M.↗

NASA NTRS1969

Dynamic characteristics of a variable-mass flexible missile: Dynamics of a two-stage variable-mass flexible rocket

The dynamic characteristics of two-stage slender elastic body were investigated. The first stage, containing a solid-fuel rocket, possesses variable mass while the second stage, envisioned as a flexible case, contains packaged instruments of constant mass. The mathematical formulation was in terms of vector equations of motion transformed by a variational principle into sets of scalar differential equations in terms of generalized coordinates. Solutions to the complete equations were obtained numerically by means of finite difference techniques. The problem has been programmed in the FORTRAN 4 language and solved on an IBM 360/50 computer. Results for limited cases are presented showing the nature of the solutions.

Meirovitch, L.↗

NASA NTRS1972

Flow field predictions for a slab delta wing at incidence

Theoretical results are presented for the structure of the hypersonic flow field of a blunt slab delta wing at moderately high angle of attack. Special attention is devoted to the interaction between the boundary layer and the inviscid entropy layer. The results are compared with experimental data. The three-dimensional inviscid flow is computed numerically by a marching finite difference method. Attention is concentrated on the windward side of the delta wing, where detailed comparisons are made with the data for shock shape and surface pressure distributions. Surface streamlines are generated, and used in the boundary layer analysis. The three-dimensional laminar boundary layer is computed numerically using a specially-developed technique based on small cross-flow in streamline coordinates. In the rear sections of the wing the boundary layer decreases drastically in the spanwise direction, so that it is still submerged in the entropy layer at the centerline, but surpasses it near the leading edge. Predicted heat transfer distributions are compared with experimental data.

Conti, R. J.↗

NASA NTRS1971

Evaluation of two statistical models using the shock structure problem.

The accuracy of two statistical models for the collision integral of the Boltzmann equation has been evaluated by applying the models to the solution of the problem of shock structure in a monatomic gas and then comparing the theoretical results with available ex perimental data. The two models considered here are the Bhatnagar-Gross-Krook and ellipsoidal statistical models. The Mach number range covered is 1.59-10.7 and profiles for density and, where available, temperature are compared. The method of numerical solution is the discrete ordinate technique which looks quite promising for application to more complicated models. The results indicate that the ellipsoidal statistical model, which gives a correct value for the Prandtl number, gives accurate results for a low Mach number shock. However, the accuracy degenerates as the Mach number increases. The Bhatnagar-Gross-Krook model gives poorer agreement with experimental data in all cases examined.

Giddens, D. P.↗