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 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Numerical simulation of two-dimensional heat transfer in composite bodies with application to de-icing of aircraft components

Transient, numerical simulations of the de-icing of composite aircraft components by electrothermal heating were performed for a two dimensional rectangular geometry. The implicit Crank-Nicolson formulation was used to insure stability of the finite-difference heat conduction equations and the phase change in the ice layer was simulated using the Enthalpy method. The Gauss-Seidel point iterative method was used to solve the system of difference equations. Numerical solutions illustrating de-icer performance for various composite aircraft structures and environmental conditions are presented. Comparisons are made with previous studies. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Chao, D. F. K.↗

NASA NTRS1983

Analytical and Numerical Solution for a Solidifying Liquid Alloy Slab

Numerical and analytical solutions are presented for the temperature and concentration distributions during the solidification of a binary liquid alloy slab. The slab is taken to be of a finite depth but infinite in the horizontal direction. The solidification process is started by withdrawing a fixed amount of heat from the lower surface of the slab. The upper surface of the slab is subjected to both radiation and convective conditions. The solution gives the concentration and temperature profiles and the interface position as a function of time. Due to the smallness of the mass diffusion coefficient in the solid, the numerical solution method breaks down whenever the ratio of the diffusivities in the solid and the liquid falls below a certain value. An analytical method is developed which gives accurate solution for any value of the diffusivity ratio.

Antar, B. N.↗

NASA NTRS1983

Numerical solution of the Navier-Stokes equations for arbitrary 2-dimensional multi-element airfoils

Numerical solutions of the Navier-Stokes equations, with an algebraic turbulence model, for time-dependent two dimensional flow about multi-element airfoils were developed. Fundamental to these solutions was the use of numerically-generated boundary-conforming curvilinear coordinate systems to allow bodies of arbitrary shape to be treated. A general two dimensional grid generation code for multiple-body configuration was written as a part of this project and made available through the COSMIC code library.

Thompson, J. F.↗

NASA NTRS1984

Operational numerical weather prediction on the CYBER 205 at the National Meteorological Center

The Development Division of the National Meteorological Center (NMC), having the responsibility of maintaining and developing the numerical weather forecasting systems of the center, is discussed. Because of the mission of NMC data products must be produced reliably and on time twice daily free of surprises for forecasters. Personnel of Development Division are in a rather unique situation. They must develop new advanced techniques for numerical analysis and prediction utilizing current state-of-the-art techniques, and implement them in an operational fashion without damaging the operations of the center. With the computational speeds and resources now available from the CYBER 205, Development Division Personnel will be able to introduce advanced analysis and prediction techniques into the operational job suite without disrupting the daily schedule. The capabilities of the CYBER 205 are discussed.

Deaven, D.↗

NASA NTRS1984

Dynamics analysis of electrodynamic satellite tethers. Equations of motion and numerical solution algorithms for the tether

The equations of motion are developed for a perfectly flexible, inelastic tether with a satellite at its extremity. The tether is attached to a space vehicle in orbit. The tether is allowed to possess electrical conductivity. A numerical solution algorithm to provide the motion of the tether and satellite system is presented. The resulting differential equations can be solved by various existing standard numerical integration computer programs. The resulting differential equations allow the introduction of approximations that can lead to analytical, approximate general solutions. The differential equations allow more dynamical insight of the motion.

Nacozy, P. E.↗

NASA NTRS1983

Aerodynamic design using numerical optimization

The procedure of using numerical optimization methods coupled with computational fluid dynamic (CFD) codes for the development of an aerodynamic design is examined. Several approaches that replace wind tunnel tests, develop pressure distributions and derive designs, or fulfill preset design criteria are presented. The method of Aerodynamic Design by Numerical Optimization (ADNO) is described and illustrated with examples.

Murman, E. M.↗

NASA NTRS1984

Numerical simulation of boundary-layer transition

The transition to turbulence in boundary layers was investigated by direct numerical solution of the nonlinear, three-dimensional, incompressible Navier-Stokes equations in the half-infinite domain over a flat plate. Periodicity was imposed in the streamwise and spanwise directions. A body force was applied to approximate the effect of a nonparallel mean flow. The numerical method was spectra, based on Fourier series and Jacobi polynomials, and used divergence-free basis functions. Extremely rapid convergence was obtained when solving the linear Orr-Sommerfeld equation. The early nonlinear and three-dimensional stages of transition, in a boundary layer disturbed by a vibrating ribbon, were successfully simulated. Excellent qualitative agreement was observed with either experiments or weakly nonlinear theories. In particular, the breakdown pattern was staggered or nonstaggered depending on the disturbance amplitude.

Spalart, P. R.↗

NASA NTRS1983

Numerical approximations to nonlinear conservation laws with locally varying time and space grids

Numerical approximations to the initial value problem for nonlinear systems of conservation laws are considered. The considered system is said to be hyperbolic when all eigenvalues of every real linear combination of the Jacobian matrices are real. Solutions may develop discontinuities in finite time, even when the initial data are smooth. In the investigation, explicit finite difference methods which use locally varying time grids are considered. The global CFL restriction is replaced by a local restriction. The numerical flux function is studied from a finite volume viewpoint, and a differencing technique is developed at interface points between regions of distinct time increments.

Osher, S.↗

NASA NTRS1984

Comparison of an analytical approximation for particle motion in a current sheet with precise numerical calculations

Approximate analytic solutions exist for particle motion in a one-dimensional current sheet with a constant normal magnetic field component. These solutions are tested against precise numerical calculations, and a range of validity of the analytic solutions is inferred. For example, in the geomagnetic tail neutral sheet, for a dawn-dusk electric field of 0.1-1 mV/m, lobe field of 10-40 nT, and sheet thickness of 1000 km, the analytic solutions serve as a good predictor of particle motion when the normal magnetic field component is less than 3 or 4 nT. By using the analytic solutions, initial distribution functions are mapped into final (accelerated) distributions, and the analytic mappings are compared with numerical mappings.

Speiser, T. W.↗

NASA NTRS1984

A numerical weather prediction system designed to simulate atmospheric downburst phenomena

It is pointed out that an increase in the understanding of weather-related aircraft accidents can save hundreds of human lives and million of dollars. A better understanding regarding the interaction between aircraft operation and severe weather conditions can be obtained with the aid of flight simulator facilities. It is shown that numerical weather modeling is one of the most precise and cost-effective inputs for flight simulators in the long run. A comprehensive weather modeling system is being developed for the simulation of different scales of atmospheric phenomena. The modeling system utilizes two numerical weather models, including the Mesoscale Atmospheric Simulation system, and the Terminal Area Simulation System.

Chuang, S.↗

NASA NTRS1984

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem. Previously announced in STAR as N83-30109

Phillips, T. N.↗

NASA NTRS1984

Numerical simulations of the Ross Sea tides

Tidal currents below the floating Ross Ice shelf are reconstructed by using a numerical tidal model. They are predominantly diurnal, achieve maximum strength in regions near where the ice shelf runs aground, and are significantly enhanced by topographic Rossby wave propagation along the ice front. A comparison with observations of the vertical motion of the ice shelf surface indicates that the model reproduces the diurnal tidal characteristics within 20 percent. Similar agreement for the relatively weak semidiurnal tides was not obtained, and this calls attention to possible errors of the open boundary forcing obtained from global-ocean tidal simulations and to possible errors in mapping zones of ice shelf grounding. Air-sea contact below the ice shelf is eliminated by the thick ice cover. The dominant sub-ice-shelf circulation may thus be tidally induced. A preliminary assessment of sub-ice-shelf conditions based on the numerical tidal simulations suggests that (1) strong barotropic circulation is driven along the ice front and (2) tidal fronts may form in the sub-ice-shelf cavity where the water column is thin and where the buoyancy input is weak.

Macayeal, D. R.↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included. Previously announced in STAR as N82-28577

Moin, P.↗

NASA NTRS1982

Calculation of pressure statistics in turbulent free shear flows by direct numerical simulation

Accurate numerical simulations of turbulent free shear flows were performed in order to study the behavior of certain critical turbulence quantities which have been difficult to measure in the laboratory. The three-dimensional time-dependent nonlinear Navier-Stokes equations are solved by the direct numerical simulation technique without turbulence modeling or ad hoc closure assumptions. The pressure interaction terms were factored into diffusion and return-to-isotropy components according to two models. Simulations of the evolution of a turbulent wake with and without swirl present were performed at sufficiently high Reynolds numbers to insure that the flows were strongly nonlinear and turbulent, but low enough so that all scales of motion containing significant energy are adequately resolved. It is noted that the presence of swirl tends to enhance the amplitude of the rms pressure levels.

Metcalfe, R. W.↗

NASA NTRS1984

Interactive phenomena in supersonic jet mixing problems. I Phenomenology and numerical modeling techniques

The interactive phenomena that occur in supersonic jet mixing flowfields, and numerical modeling techniques developed to analyze such phenomena are discussed. A spatial marching procedure based on solving the parabolized Navier-Stokes jet mixing equations is presented. This procedure combines shock-capturing methodology for the analysis of supersonic mixing regions with pressure-split methodology for the analysis of subsonic mixing regions. The two regions are coupled at viscous sonic lines utilizing a viscous-characteristic coupling procedure. Specialized techniques for the treatment of jet boundary growth, strong discontinuties (Mach disks), and small embedded subsonic zones (behind Mach disks) are presented. Turbulent processes are represented by two-equation turbulence model formulations. In Part II of this article, numerical studies are presented for a variety of supersonic jet interactive phenomena.

Dash, S. M.↗

NASA NTRS1984

Numerical simulation of controlled flow tunnel for V/STOL testing

A 'smart' wind tunnel concept for V/STOL transition regime studies is presented and performance is projected with numerical simulations. The wind tunnel would provide a flow that faithfully represents free air flow and thereby permits accurate lift calculations for V/STOL aircraft. The design concept is based on free-air tunnels used for planar wings. Modifications are necessary to account for the highly deflected, vortical wake produced in low-speed take-off and landing situations. A three-step numerical procedure covering the potential flow of a high lift system in free-air, controlled flow data with a model in a wind tunnel and addition of the wall conditions to the simulation as a feedback loop to correct the closed-tunnel is described. The initial simulation is for three-dimensional full span jet flapped wings, and is derived from extended two-dimensional calculations. The results indicate that only partial control of the tunnel flow is sufficient to serve as a valid environment for V/STOl studies.

Joppa, R. G.↗

NASA NTRS1984

A numerical simulation of barotropic instability. I Wave-mean flow interaction

A numerical model is used to study the evolution of the barotropic point jet instability as it interacts with the mean flow. The linearized instability solution agrees well with the recent analytical solutions of Lindzen. Stabilization of the point jet instability occurs as the mean flow is modified by wave vorticity transport. Assuming stabilization occurs when the meridional gradient of the zonal mean vorticity is no longer negative, the maximum integrated wave enstrophy can be predicted. In addition, an estimate of the integrated wave enstrophy at steady state can be made by balancing the generation of vorticity against dissipation. These limits are found to be in good agreement with the numerical results.

Schoeberl, M. R.↗

NASA NTRS1984

Numerical optimization of circulation control airfoil at high subsonic speed

A numerical procedure for optimizing the design of the circulation control airfoil for use at high subsonic speeds is presented. The procedure consists of an optimization scheme coupled with a viscous potential flow analysis for the blowing jet. The desired airfoil is defined by a combination of three baseline shapes (cambered ellipse and cambered ellipse with drooped and spiraled trailing edges). The coefficients of these shapes are used as design variables in the optimization process. Under the constraints of lift augmentation and lift-to-drag ratios, the airfoil, optimized at free-stream Mach 0.54 and alpha = -2 degrees can be characterized as a cambered ellipse with a drooped trailing edge. Experimental tests support the performance improvement predicted by numerical optimization.

Tai, T. C.↗