Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 1,135 records · Page 63

Numerical calculation of transonic axial turbomachinery flows

This paper presents a numerical method and the results of a computer program for solving an exact, three-dimensional, full-potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analysing flow fields about arbitrarily shaped wing-body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type-dependent differencing, successive line overrelaxation, and sequential boundary-conforming grid refinement. An artificial viscosity is added in fully conservative form; and an initial guess for the potential field is applied, as determined by a two-dimensional cascade analysis.

Dulikravich, D. S.↗

Numerical solutions of the Navier-Stokes equations for transonic afterbody flows

The time dependent Navier-Stokes equations in mass averaged variables are solved for transonic flow over axisymmetric boattail plume simulator configurations. Numerical solution of these equations is accomplished with the unsplit explict finite difference algorithm of MacCormack. A grid subcycling procedure and computer code vectorization are used to improve computational efficiency. The two layer algebraic turbulence models of Cebeci-Smith and Baldwin-Lomax are employed for investigating turbulence closure. Two relaxation models based on these baseline models are also considered. Results in the form of surface pressure distribution for three different circular arc boattails at two free stream Mach numbers are compared with experimental data. The pressures in the recirculating flow region for all separated cases are poorly predicted with the baseline turbulence models. Significant improvements in the predictions are usually obtained by using the relaxation models.

Swanson, R. C., Jr.↗

Penetrative convection

The current state of understanding of the most directly observable solar convection, the granulation and supergranulation is summarized. The body of work in which the complete time dependent Navier-Stokes equations and entropy transport equation are solved for a fully compressible atmosphere is considered. Relevant anelastic and incompressible calculations in two dimensions are also discussed.

Moore, D. R.↗

New discretization and solution techniques for incompressible viscous flow problems

Several topics arising in the finite element solution of the incompressible Navier-Stokes equations are considered. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. The role of artificial viscosity in viscous flow calculations is studied, emphasizing work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some modifications are mentioned.

Gunzburger, M. D.↗

New discretization and solution techniques for incompressible viscous flow problems

This paper considers several topics arising in the finite element solution of the incompressible Navier-Stokes equations. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. Following this, the role of artificial viscosity in viscous flow calculations is studied, emphasizing recent work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some recent modifications are mentioned.

Gunzburger, M. D.↗

Three-dimensional flow over a conical afterbody containing a centered propulsive jet - A numerical simulation

The supersonic flow field over a body of revolution incident to the free stream is simulated numerically on a large, array processor (the CDC Cyber 205). The configuration is composed of a cone-cylinder forebody followed by a conical afterbody from which emanates a centered, supersonic propulsive jet. The free-stream Mach number is 2, the jet-exit Mach number is 2.5, and the jet-to-free-stream static pressure ratio is 3. Both the external flow and the exhaust are ideal air at a common total temperature. The thin-layer approximation to the time-dependent, compressible, Reynolds-averaged Navier-Stokes equations are solved using an implicit finite-difference algorithm. The data base, of 5 million words, is structured in a 'pencil' format so that efficient use of the array processor can be realized. The computer code is completely vectorized to take advantage of the data structure. Turbulence closure is achieved using an empirical algebraic eddy-viscosity model. The configuration and flow conditions correspond to published experimental tests and the computed solutions are consistent with the experimental data.

Deiwert, G. S.↗

Analysis of three-dimensional transonic compressors

A method for computing the three-dimensional transonic flow around the blades of a compressor or of a propeller is given. The method is based on the use of the velocity potential, on the hypothesis that the flow is inviscid, irrotational and isentropic. The equation of the potential is solved in a transformed space such that the surface of the blade is mapped into a plane where the periodicity is implicit. This equation is in a nonconservative form and is solved with the help of a finite difference method using artificial time. A computer code is provided and some sample results are given in order to demonstrate the influence of three-dimensional effects and the blade's rotation.

Bourgeade, A.↗

Exact solution of three-dimensional transport problems using one-dimensional models

Several parameters of certain three-dimensional semiconductor devices including diodes, transistors, and solar cells can be determined without solving the actual boundary-value problem. The recombination current, transit time, and open-circuit voltage of planar diodes are emphasized here. The resulting analytical expressions enable determination of the surface recombination velocity of shallow planar diodes. The method involves introducing corresponding one-dimensional models having the same values of these parameters.

Misiakos, K.↗

Rapid computation of chemical equilibrium composition - An application to hydrocarbon combustion

A scheme for rapidly computing the chemical equilibrium composition of hydrocarbon combustion products is derived. A set of ten governing equations is reduced to a single equation that is solved by the Newton iteration method. Computation speeds are approximately 80 times faster than the often used free-energy minimization method. The general approach also has application to many other chemical systems.

Erickson, W. D.↗

Temporal features of the refilling of a plasmaspheric flux tube

The refilling of plasmaspheric flux tubes was studied by assuming that the protonosphere provides an ionospheric boundary where the H(+) density can be assumed; the supersonic flow in the flux tube is driven by the depletion of the plasma from the flux tube, while the base density and pressure in the protonosphere remain constant. The time-dependent continuity and momentum equations for the H(+) ions were solved. The electron gas was assumed to obey the Boltzmann law, and the proton gas was assumed to be isothermal. In agreement with the postulate of Banks et al. (1971), it was found that an important feature of the refilling is the formation of a shock pair at the equator; as the shocks propagate toward the ionosphere, the refilling occurs. Depending on the density at the ionospheric boundaries, a fair agreement was found between the refilling rates obtained for L = 6.6 and those from the GEOS 2 observations.

Singh, Nagendra↗

A technique for the prediction of airfoil flutter characteristics in separated flow

A solution procedure is described for determining the two-dimensional, one- or two-degree-of-freedom flutter characteristics of arbitrary airfoils at large angles of attack. The same procedure is used to predict stall flutter. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. Results for a number of special cases are presented to demonstrate the suitability of this scheme to predict flutter at large mean angles of attack. Some stall flutter applications are also presented.

Wu, Jiunn-Chi↗

Range and mission scheduling automation using combined AI and operations research techniques

Ground-based systems for Satellite Command, Control, and Communications (C3) operations require a method for planning, scheduling and assigning the range resources such as: antenna systems scattered around the world, communications systems, and personnel. The method must accommodate user priorities, last minute changes, maintenance requirements, and exceptions from nominal requirements. Described are computer programs which solve 24 hour scheduling problems, using heuristic algorithms and a real time interactive scheduling process.

Arbabi, Mansur↗

Application of Navier-Stokes analysis to stall flutter

A solution procedure was developed to investigate the two-dimensional, one- or two-dimensional flutter characteristics of arbitrary airfoils. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. The flutter characteristics of airfoils in subsonic speed at high angles of attack and airfoils in high subsonic and transonic speeds at low angles of attack are investigated. The stall flutter characteristics are also predicted using the same procedure.

Wu, J. C.↗

A computational study of heat transfer for laminar hypervelocity flows

Computational Fluid Dynamics (CFD) codes are routinely used to predict the flowfield and the heating environment around complex reentry configurations. At hypervelocities, where the velocity is greater than 3 km/sec, the AFWAL version of the blunt body code predicts the correct surface pressure distributions but underpredicts laminar wall heat fluxes. This study was performed to determine the reasons for the underprediction. The computer code chosen solves thin-layer Navier-Stokes equations in a time-asymptotic manner and assumes a constant isentropic exponent. Flowfields around a spherical configuration at various entry velocities are computed. The computed pressure distributions agree well with the tabulated, inviscid results of Lyubimov and Rusanov for entry velocities ranging from 0.6 to 5.92 km/sec. At hypervelocities, the calculated stagnation point heat transfer rates were lower by roughly fifty percent when compared to engineering correlations available in the literature. Good comparisons between heat transfer rates are obtained at hypervelocity entry conditions provided the CFD code is modified to include equilibrium air properties.

Balakrishnan, A.↗

Cracking of coated materials under transient thermal stresses

The crack problem for a relatively thin layer bonded to a very thick substrate under thermal shock conditions is considered. The effect of surface cooling rate is studied by assuming the temperature boundary condition to be a ramp function. Among the crack geometries considered are the edge crack in the coating layer, the broken layer, the edge crack going through the interface, the undercoat crack in the substrate and the embedded crack crossing the interface. The primary calculated quantity is the stress intensity factor at various singular points and the main variables are the relative sizes and locations of cracks, the time, and the duration of the cooling ramp. The problem is solved and rather extensive results are given for two material pairs, namely a stainless steel layer welded on a ferritic medium and a ceramic coating on a steel substrate.

Rizk, A. A.↗

Unsteady stagnation-point heat transfer during passage of a concentrated vortex

The unsteady boundary layer due to a single rectilinear vortex filament approaching a 2-D stagnation point is investigated. Assuming the vortex remains far from the surface, incompressible potential flow theory is used to determine the time dependent inviscid flow field. The unsteady boundary layer equations are solved by an alternating-direction-implicit finite-difference method. Two mechanisms which cause fluctuations in heat transfer are the unsteady velocity field in the boundary layer and secondly, the unsteady total temperature at the edge of the boundary layer. The relative importance of these mechanisms is dependent upon the total temperature fluctuations relative to the imposed temperature difference. As a vortex approaches a stagnation point it may be forced to one side of the stagnation line or the other, depending on its initial position. Results are presented for both of these cases.

Rigby, David L.↗

Solar coronal loop heating by cross-field wave transport

Solar coronal arches heated by turbulent ion-cyclotron waves may suffer significant cross-field transport by these waves. Nonlinear processes fix the wave-propagation speed at about a tenth of the ion thermal velocity, which seems sufficient to spread heat from a central core into a large cool surrounding cocoon. Waves heat cocoon ions both through classical ion-electron collisions and by turbulent stochastic ion motions. Plausible cocoon sizes set by wave damping are in roughly kilometers, although the wave-emitting core may be only 100 m wide. Detailed study of nonlinear stabilization and energy-deposition rates predicts that nearby regions can heat to values intermediate between the roughly electron volt foot-point temperatures and the about 100 eV core, which is heated by anomalous Ohmic losses. A volume of 100 times the core volume may be affected. This qualitative result may solve a persistent problem with current-driven coronal heating; that it affects only small volumes and provides no way to produce the extended warm structures perceptible to existing instruments.

Amendt, Peter↗