Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solutions”

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,801 records · Page 100

A comparison of Euler and Navier-Stokes solutions for supersonic flow over a conical delta wing

The flow around a thin elliptic cone at Mach number 2 and angle of attack 10 deg is simulated numerically, comparing the results obtained using a fine grid and the Navier-Stokes (NS) equations with those obtained using both coarse and fine grids and the Euler (E) equations. On the coarse grid, the E results are found to reproduce roughly the NS predictions of leading-edge separation and primary and secondary vortices on the leeward surface, but the E solution on the fine NS grid (showing no separation and a crossflow-shock-induced vorticity) reveals this agreement to be spurious. The fine-grid E prediction is shown to be a valid solution of the E equations, indicating their unsuitability as approximators of viscous flow. The imposition of a Kutta condition on the E computation is found to be ineffective in correcting this discrepancy.

Newsome, R. W.↗

Solution method for a hovering helicopter rotor using the Euler equations

A method for the calculation of the flow field of a helicopter rotor blade in hover is presented. The approach uses a finite volume solution of the three dimensional Euler equations for the blade near field. In the Euler solver the velocity field is decomposed into two parts. One is the induced velocity of the vortex wake extending below the blade, found from a free wake calculation procedure. The other part is the unknown additional velocity field of the rotor blade. This approach eliminates numerical diffusion of the rolled up wake vorticity due to truncation error and artificial viscosity. Also, the effects of the far wake are included in the limited computational domain. Solutions are presented for an isolated wing and a model helicopter rotor and compared to experiment.

Roberts, T. W.↗

Efficient self-consistent viscous-inviscid solutions for unsteady transonic flow

An improved method is presented for coupling a boundary layer code with an unsteady inviscid transonic computer code in a quasi-steady fashion. At each fixed time step, the boundary layer and inviscid equations are successively solved until the process converges. An explicit coupling of the equations is described which greatly accelerates the convergence process. Computer times for converged viscous-inviscid solutions are about 1.8 times the comparable inviscid values. Comparison of the results obtained with experimental data on three airfoils are presented. These comparisons demonstrate that the explicitly coupled viscous-inviscid solutions can provide efficient predictions of pressure distributions and lift for unsteady two-dimensional transonic flows.

Howlett, J. T.↗

Recent advances in methods for numerical solution of O.D.E. initial value problems

In the mathematical modeling of physical systems, it is often necessary to solve an initial value problem (IVP), consisting of a system of ordinary differential equations (ODE). A typical program produces approximate solutions at certain mesh points. Almost all existing codes try to control the local truncation error, while the user is really interested in controlling the true or global error. The present investigation provides a review of recent advances regarding the solution of the IVP, giving particular attention to stiff systems. Stiff phenomena are customarily defined in terms of the eigenvalues of the Jacobian. There are, however, some difficulties connected with this approach. It is pointed out that an estimate of the Lipschitz constant proves to be a very practical way to determine the stiffness of a problem.

Bui, T. D.↗

On the development of hierarchical solution strategies for nonlinear finite element formulations

This paper develops a hierarchical type solution scheme which can handle the field equations associated with nonlinear finite element simulations. The overall procedure possesses various levels of application namely degree of freedom, nodal, elemental, substructural as well as global. In particular iteration, updating, assembly and solution control occurs at the various hierarchical levels. Due to the manner of formulation, the degree of matrix inversion depends on the size of the various hierarchical partitioned groups. In this context, degree of freedom partitioning requires no inversion. To benchmark the overall scheme, the results of several numerical examples are presented.

Padovan, J.↗

Approximate analytic solutions for singular non-linear oscillators

Mickens (1981, 1984) has considered analytic techniques for obtaining approximate solutions to one-dimensional nonlinear oscillatory systems x(double-dot) + x = lambda f(x, x/dot/, lambda) where lambda is a small positive parameter and f is a nonlinear polynomial function of its arguments. However, in certain cases there is interest in the analysis of physical systems for which the nonlinear function f(x, x/dot/, lambda) is singular for finite values of x or x(dot). The present investigation is concerned with the use of existing approximate analytic schemes to obtain solutions to singular nonlinear oscillatory differential equations.

Bota, K. B.↗

Multigrid techniques for the solution of the passive scalar advection-diffusion equation

The solution of elliptic passive scalar advection-diffusion equations is required in the analysis of many turbulent flow and convective heat transfer problems. The accuracy of the solution may be affected by the presence of regions containing large gradients of the dependent variables. The multigrid concept of local grid refinement is a method for improving the accuracy of the calculations in these problems. In combination with the multilevel acceleration techniques, an accurate and efficient computational procedure is developed. In addition, a robust implementation of the QUICK finite-difference scheme is described. Calculations of a test problem are presented to quantitatively demonstrate the advantages of the multilevel-multigrid method.

Phillips, R. E.↗

Closed-form solutions for a class of optimal quadratic tracking problems

Closed-form solutions are derived for a class of tracking problems including a linear optimal regulator and a prefilter for a time-invariant plant. The solutions for the prefilter equation and state trajectory coupled by the Riccati equation are exponentially related to the stability matrix of the plant. A computational procedure is presented in recursive form when the desired output state dynamics is assumed linear and time-invariant. Several examples are given for illustration.

Turner, J. D.↗

Refined numerical solution of the transonic flow past a wedge

A numerical procedure combining the ideas of solving a modified difference equation and of adaptive mesh refinement is introduced. The numerical solution on a fixed grid is improved by using better approximations of the truncation error computed from local subdomain grid refinements. This technique is used to obtain refined solutions of steady, inviscid, transonic flow past a wedge. The effects of truncation error on the pressure distribution, wave drag, sonic line, and shock position are investigated. By comparing the pressure drag on the wedge and wave drag due to the shocks, a supersonic-to-supersonic shock originating from the wedge shoulder is confirmed.

Liang, S.-M.↗

Growth of Solid Solution Crystals

The major objective of this program is to determine the conditions under which single crystals of solid solutions can be grown from the melt in a Bridgman configuration with a high degree of chemical homogeneity. The central aim is to assess the role of gravity in the growth process and to explore the possible advantages for growth in the absence of gravity. The alloy system being investigated is the solid solution semiconductor with x-values appropriate for infrared detector applications in Hg sub (1-x) Cd sub x Te the 8 to 14 micro m wavelength region. Both melt and Te-solvent growth are being considered. The study consists of an extensive ground-based experimental and theoretical research effort followed by flight experimentation where appropriate. Experimental facilities have been established for the purification, casting, and crystal growth of the alloy system. Facilities have been also established for the metallurgical, compositional, electric and optical characterization of the alloys. Crystals are being grown by the Bridgman-Stockbarger method and are analyzed by various experimental techniques to evaluate the effects of growth conditions on the longitudinal and radial compositional variations and defect densities in the crystals.

Lehoczky, S. L.↗

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Abolhassani, J. S.↗

Electrochemical control of pH in a hydroponic nutrient solution

The electrochemical pH control system described was found to provide a feasible alternative method of controlling nutrient solution pH for CELSS applications. The plants grown in nutrient solution in which the pH was controlled electrochemically showed no adverse effects. Further research into the design of a larger capacity electrode bridge for better control is indicated by the results of this experiment, and is currently under way.

Schwartzkopf, S. H.↗

Cyclic creep analysis from elastic finite-element solutions

A uniaxial approach was developed for calculating cyclic creep and stress relaxation at the critical location of a structure subjected to cyclic thermomechanical loading. This approach was incorporated into a simplified analytical procedure for predicting the stress-strain history at a crack initiation site for life prediction purposes. An elastic finite-element solution for the problem was used as input for the simplified procedure. The creep analysis includes a self-adaptive time incrementing scheme. Cumulative creep is the sum of the initial creep, the recovery from the stress relaxation and the incremental creep. The simplified analysis was exercised for four cases involving a benchmark notched plate problem. Comparisons were made with elastic-plastic-creep solutions for these cases using the MARC nonlinear finite-element computer code.

Kaufman, A.↗

Effect of triangular element orientation on finite element solutions of the Helmholtz equation

The Galerkin finite element solutions for the scalar homogeneous Helmholtz equation are presented for no reflection, hard wall, and potential relief exit terminations with a variety of triangular element orientations. For this group of problems, the correlation between the accuracy of the solution and the orientation of the linear triangle is examined. Nonsymmetric element patterns are found to give generally poor results in the model problems investigated, particularly for cases where standing waves exist. For a fixed number of vertical elements, the results showed that symmetric element patterns give much better agreement with corresponding exact analytical results. In laminated wave guide application, the symmetric pyramid pattern is convenient to use and is shown to give excellent results.

Baumeister, K. J.↗

Solution of viscous transonic flow over wings

Since the calculations of plane steady transonic flows conducted by Murman and Cole (1971), steady progress has been made with respect to the computation of inviscid transonic flows. It has been found that it is inadequate to consider practical wing design at transonic speeds without considering the effects of viscosity and turbulence. The present study has the objective to develop a zonal method for viscous transonic flow over three-dimensional (3-D) wings. The employed approach follows closely the viscous/inviscid interaction techniques discussed by Melnik et al. (1983) for viscous flow about airfoils. Attention is given to inviscid flow equations and boundary conditions, the solution of 3-D boundary layer and wake, an iterative solution to viscid-inviscid interaction analysis, and results obtained for a transonic cruise wing of transport type.

Chow, R. R.↗

On the solution topologies of polytropic winds

Steady polytropic wind flows are studied by examining their solution topologies. The problem of whether the degeneracy of the continuous and standing shocked solutions in the isothermal case persist in the more general case of a polytropic expanding atmosphere is addressed. The effects of departures from spherical symmetry and nonthermal momentum addition are considered within the context of the new class of magnetospheric models proposed by Low and Tsinganos (1985). The flow near the boundary of a polar coronal hole is considered for the specific case of a magnetosphere in which this boundary asymptotically approaches the equator at large heliocentric distances.

Bailyn, C.↗

Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies

The particular problems investigated in the present study arise from nonlinear transient heat conduction. One of two types of nonlinearities considered is related to a material temperature dependence which is frequently needed to accurately model behavior over the range of temperature of engineering interest. The second nonlinearity is introduced by radiation boundary conditions. The finite element equations arising from the solution of nonlinear transient heat conduction problems are formulated. The finite element matrix equations are temporally discretized, and a nonlinear iterative solution algorithm is proposed. Algorithms for solving the linear problem are discussed, taking into account the form of the matrix equations, Gaussian elimination, cost, and iterative techniques. Attention is also given to approximate factorization, implementational aspects, and numerical results.

Winget, J. M.↗

Refine your search to explore more results.