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.

SourceResearch connectionFrom yearThrough year

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

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗

NASA NTRS1984

Canonical solutions for unsteady flow fields

The initial value problem of one-dimensional gas-dynamics involving discontinuous, nonuniform initial data is discussed. Canonical solutions which are valid in a small x, t region aroung a discontinuity, and which include the first order effects of nonuniformities in the data, are derived explicitly. The theory is derived by considering a group of elementary piston problems. Solutions with a shock or with a centered expansion wave are worked out individually in order to relate initial flow properties and their gradients to the speed and acceleration of the discontinuity waves. They are then combined to represent the solution of a general initial value problem by regarding the piston path as a contact line. In addition, problems with chemical reaction are discussed in terms of elementary piston problems which involve strong detonation waves, Chapman-Jouguet detonation waves, and deflagration waves.

Liu, G. C.↗

NASA NTRS1984

Numerical solution for the interaction of shock wave with laminar boundary layer in two-dimensional flow on a flat plate

The finite difference computation method was investigated for solving problems of interaction between a shock wave and a laminar boundary layer, through solution of the complete Navier-Stokes equations. This method provided excellent solutions, was simple to perform and needed a relatively short solution time. A large number of runs for various flow conditions could be carried out from which the interaction characteristics and principal factors that influence interaction could be studied.

Landau, U.↗

NASA NTRS1984

Transonic solutions for a multielement airfoil using the full-potential equation

Transonic flow solutions are obtained over a multielement airfoil (augmentor-wing) using the full-potential equation. Solutions obtained for a subcritical case and a strong shock case show good quantitative agreement with experiment in regions not dominated by viscous effects. In those regions where viscous effects are dominant, the results are still in good qualitative agreement. For the strong shock case, Mach number and angle-of-attack corrections were necessary to match experimental coefficient of lift. Typical results from the transonic augmentor-wing Potential Code on the Cray-1S computer require about 10 sec of CPU time for a three-order-of-magnitude drop in the maximum residual. The speed with which solutions can be generated, and the associated low cost, will make this code a practical tool for the design aerodynamicist.

Flores, J.↗

NASA NTRS1983

A numerical model of acoustic choking. I - Shock free solutions

The phenomenon of acoustic choking in near sonic flows is investigated using the one-dimensional gasdynamic equations for an ideal gas. It is pointed out that this approach eliminates the need to make the classical small disturbance assumption. A finite difference scheme is elaborated to approximate these equations. Boundary conditions that correctly model the physics of the acoustics problems are chosen. The use of a fast Fourier transform routine allows the results to be compared with those obtained by using the conventional harmonic assumptions. Only solutions in which shocks have not fully developed are reported. Continuous solutions are sought, and these are compared with the Fubini solution and with certain results of Myers and Callegari (1977, 1978, 1979).

Walkington, N. J.↗

NASA NTRS1982

Computation of consistent boundary quantities in finite element thermal-fluid solutions

The consistent boundary quantity method for computing derived quantities from finite element nodal variable solutions is investigated. The method calculates consistent, continuous boundary surface quantities such as heat fluxes, flow velocities, and surface tractions from nodal variables such as temperatures, velocity potentials, and displacements. Consistent and lumped coefficient matrix solutions for such problems are compared. The consistent approach may produce more accurate boundary quantities, but spurious oscillations may be produced in the vicinity of discontinuities. The uncoupled computations of the lumped approach provide greater flexibility in dealing with discontinuities and provide increased computational efficiency. The consistent boundary quantity approach can be applied to solution boundaries other than those with Dirichlet boundary conditions, and provides more accurate results than the customary method of differentiation of interpolation polynomials.

Thornton, E. A.↗

NASA NTRS1984

Navier-Stokes solutions for two-dimensional subsonic base flow

Methods for determining the effects of mass injection from the trailing edge of a bluff body at low speeds and in transonic flow were numerically studied along with an unmodified blunt-based body to gain insight into the effects of vortex shedding on the base drag. The methodology used to obtain finite-difference solutions to the Navier-Stokes equations for subsonic compressible two-dimensional near-wake flows is presented. The effectiveness of an introduced outflow boundary condition which minimizes reflections back into the computational domain was demonstrated with the solution of a model vortex problem. Calculations of the near-wake flow past a circular cylinder were in excellent agreement with experimental data. Laminar-flow solutions for a blunt-based model with and without a base cavity and with mass injection into the wake agreed qualitatively with experimental observations. The drag reduction capability provided by such base modifications was demonstrated.

Rudy, D. H.↗

NASA NTRS1984

A finite difference solution for the propagation of sound in near sonic flows

An explicit time/space finite difference procedure is used to model the propagation of sound in a quasi one-dimensional duct containing high Mach number subsonic flow. Nonlinear acoustic equations are derived by perturbing the time-dependent Euler equations about a steady, compressible mean flow. The governing difference relations are based on a fourth-order, two-step (predictor-corrector) MacCormack scheme. The solution algorithm functions by switching on a time harmonic source and allowing the difference equations to iterate to a steady state. The principal effect of the non-linearities was to shift aocustical energy to higher harmonics. With increased source strength, wave steepening was observed. This phenomenon suggests that the acoustical response may approach a shock behavior at higher sound pressure level as the throat Mach number approaches unity. On a peak level basis, good agreement between the nonlinear finite difference and linear finite element solutions was observed, even through a peak sound pressure level of about 150 dB occurred in the throat region. Nonlinear steady state waveform solutions are shown to be in excellent agreement with a nonlinear asymptotic theory. Previously announced in STAR as N83-30167

Hariharan, S. I.↗

NASA NTRS1984

Analytic solutions for dual-spin spacecraft during platform motion

This paper presents analytic solutions to the dynamic equations of dual-spinners during platform slews. The rotor is assumed to be an axisymmetric balanced rigid body. The platform is an asymmetric rigid body with small static and dynamic imbalances. The solutions given in this paper extend a previous theory which provided an approximate solution only in the steady state.

Hayati, S.↗

NASA NTRS1984

Closed-form solutions for a class of optimal quadratic regulator problems with terminal constraints

Closed-form solutions are derived for coupled Riccati-like matrix differential equations describing the solution of a class of optimal finite time quadratic regulator problems with terminal constraints. Analytical solutions are obtained for the feedback gains and the closed-loop response trajectory. A computational procedure is presented which introduces new variables for efficient computation of the terminal control law. Two examples are given to illustrate the validity and usefulness of the theory.

Juang, J.-N.↗

NASA NTRS1984

Error analysis and correction of discrete solutions from finite element codes

Many structures are an assembly of individual shell components. Therefore, results for stresses and deflections from finite element solutions for each shell component should agree with the equations of shell theory. This paper examines the problem of applying shell theory to the error analysis and the correction of finite element results. The general approach to error analysis and correction is discussed first. Relaxation methods are suggested as one approach to correcting finite element results for all or parts of shell structures. Next, the problem of error analysis of plate structures is examined in more detail. The method of successive approximations is adapted to take discrete finite element solutions and to generate continuous approximate solutions for postbuckled plates. Preliminary numerical results are included.

Thurston, G. A.↗

NASA NTRS1984

A solution of Mx(double dot) + Cx(dot) + Kx = 0 applicable to the design of active dampers

A solution is presented for the equations of motion for the damped linear oscillator, Mx(double dot) + Cx(dot) + Kx = 0. The algorithm solves a transformed set of equations in terms of the modal variables of the undamped system and, at the same time, solves the adjoint equation of the transformed problem. The adjoint solution is normalized to give the inverse of the solution matrix of the transformed problem. The normalized inverse is useful in design for direct computation of sensitivity derivatives of damping ratios with respect to damping rates. The algorithm is programmed to reduce storage requirements by a factor of three-fourths compared to standard complex eigenvalue subroutines. A numerical example is included.

Thurston, G. A.↗

NASA NTRS1984

An analytic solution for the orbital perturbations of the Venus Radar Mapper due to gravitational harmonics

Hill's variational equations are solved analytically for the orbital perturbations of a spacecraft nominally in an elliptic orbit around a non-spherical body. The rotation of the central planet about its spin-axis is not considered in the analysis. The perturbations are restricted to the planetary gravitational harmonics only. An extremely simple algorithm is derived to transform the spherical harmonic potentials to the orbital coordinate system, and the resulting accelerations are shown to be simply trigonometric functions of the true anomaly. With the principal matrix solution for the differential equations of the adjoint system given in closed form, the orthogonality of the trigonometric functions makes it possible to obtain an analytic solution for the non-homogeneous problem, at intervals of 2 pi in true anomaly. The solution for orbital perturbations can be extended over several revolutions by applying well-known results from Floquet's theory. The technique is demonstrated with results presented on the spacecraft periapsis altitude for the forthcoming Venus Radar Mapper Mission.

Vijayaraghavan, A.↗

NASA NTRS1984

Solution of the wave equation for open surfaces involving a line integral over the edge

A simple mathematical model of a stationary source distribution for the supersonic-propeller noise-prediction formula of Farassat (1983) is developed to test the validity of the formula solutions. The conventional thickness source term is used in place of the Isom thickness formula; the relative importance of the line and surface integrals in the solutions is evaluated; and the numerical results are compared with those obtained with a conventional retarded-time solution in tables. Good agreement is obtained over elevation angles from 10 to 90 deg, and the line-integral contribution is found to be significant at all elevation angles and of the same order of magnitude as the surface-integral contribution at angles less than 30 deg. The amplitude-normalized directivity patterns for the four cases computed (x = 1.5 or 10; k = 5.0 or 50) are presented graphically.

Farassat, F.↗

NASA NTRS1984

Existence and non-uniqueness of similarity solutions of a boundary layer problem

A Blasius boundary value problem with inhomogeneous lower boundary conditions f(0) = 0 and f'(0) = - lambda with lambda strictly positive was considered. The Crocco variable formulation of this problem has a key term which changes sign in the interval of interest. It is shown that solutions of the boundary value problem do not exist for values of lambda larger than a positive critical value lambda. The existence of solutions is proven for 0 lambda lambda by considering an equivalent initial value problem. It is found however that for 0 lambda lambda, solutions of the boundary value problem are nonunique. Physically, this nonuniqueness is related to multiple values of the skin friction.

Hussaini, M. Y.↗

NASA NTRS1984

Second-order p-iterative solution of the Lambert/Gauss problem

An algorithm is presented for efficient p-iterative solution of the Lambert/Gauss orbit-determination problem using second-order Newton iteration. The algorithm is based on a universal transformation of Kepler's time-of-flight equation and approximate inverse solutions of this equation for short-way and long-way flight paths. The approximate solutions provide both good starting values for iteration and simplified computation of the second-order term in the iteration formula. Numerical results are presented which indicate that in many cases of practical significance (except those having collinear position vectors) the algorithm produces at least eight significant digits of accuracy with just two or three steps of iteration.

Boltz, F. W.↗

NASA NTRS1984

Finite difference solutions of heat conduction problems in multi-layered bodies with complex geometries

A numerical procedure is presented for analyzing a wide variety of heat conduction problems in multilayered bodies having complex geometry. The method is based on a finite difference solution of the heat conduction equation using a body fitted coordinate system transformation. Solution techniques are described for steady and transient problems with and without internal energy generation. Results are found to compare favorably with several well known solutions.

Masiulaniec, K. C.↗

NASA NTRS1985

Gravitational lensing effects of vacuum strings - Exact solutions

Exact interior and exterior solutions to Einstein's field equations are derived for vacuum strings. The exterior solution for a uniform density vacuum string corresponds to a conical space while the interior solution is that of a spherical cap. For Mu equals 0-1/4 the external metric is ds-squared = -dt-squared + dr-squared + (1-4 Mu)-squared r-squared dphi-squared + dz-squared, where Mu is the mass per unit length in the string in Planck masses per Planck length. A maximum mass per unit length for a string is 6.73 x 10 to the 27th g/cm. It is shown that strings cause temperature fluctuations in the cosmic microwave background and produce equal brightness double QSO images separated by up to several minutes of arc. Formulae for lensing probabilities, image splittings, and time delays are derived for strings in a realistic cosmological setting. String searches using ST, the VLA, and the COBE satellite are discussed.

Gott, J. R., III↗