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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Tungsten solution kinetics and amorphization of nickel in mechanically alloyed Ni-W alloys

The kinetics of solution of W, and the subsequent amorphization of Ni, in mechanically alloyed Ni-W alloys has been investigated. As W is a highly abrasive material in the energy intensive devices used for mechanical alloying, we studied the above reactions in different mills. One used hardened steel balls as the grinding media, and the other Al2O3. Abrasion is common to both mills, but Fe wear debris from the hardened steel enters into solution in the Ni rich phases whereas Al2O3 debris is present as small dispersoids. The kinetics of W solution and those of subsequent amorphization do not appear strongly affected by the Fe in solution or the Al2O3 dispersoid. Tungsten dissolves in crystalline Ni in amounts in excess of the equilibrium solubility during alloying. Amorphization of the Ni phase occurs if the W content in this phase exceeds ca. 28 at. pct.

Aning, A. O.↗
NASA NTRS1993

Three-dimensional Navier-Stokes calculations using solution-adapted grids

A three-dimensional solution-adaptive grid generation technique is presented. The adaptation technique redistributes grid points to improve the accuracy of a flow solution without increasing the number of grid points. It is applicable to structured grids with a multiblock topology. The method uses a numerical mapping and potential theory to modify the initial grid distribution based on the properties of the flow solution on the initial grid. The technique is demonstrated with two examples - a transonic finite wing and a supersonic blunt fin. The advantages are shown by comparing flow solutions on the adapted grids with those on the initial grids.

Henderson, T. L.↗
NASA NTRS1993

On the accuracy of limiters and convergence to steady state solutions

This paper addresses the practical problem of obtaining convergence to steady state solutions when limiters are used in conjunction with upwind schemes on unstructured grids. The base scheme forms a gradient and limits it by imposing monotonicity conditions in the reconstruction stage. It is shown by analysis in one dimension that such an approach leads to various schemes meeting TVD requirements in one dimension. It is further shown that these formally second order accurate schemes are less than second order accurate in practice because of the action of the limiter function in smooth regions of the solution. Modifications are proposed to the limiter that restore the second order accuracy. In multiple dimensions these schemes produce steady state solutions that are monotone and devoid of oscillations. However, convergence stalls after a few orders of reduction in the residual. With the modified limiter, on the other hand, it is shown that converged steady state solutions can be obtained.

Venkatakrishnan, V.↗
NASA NTRS1993

Application of the multigrid solution technique to hypersonic entry vehicles

A multigrid solution procedure has been incorporated in a version of the Langley Aerothermodynamic Upwind Relaxation Algorithm. The multigrid scheme is based on the Full Approximation Storage approach and uses Full Multigrid to obtain a well defined fine mesh starting solution. Predictions were obtained using standard transfer operators and a 'V-cycle' was used to control grid sequencing. Computed hypersonic flow solutions compared with experimental data for a 15 degree sphere cone, blended-wing body, and shuttle-like geometries are presented. It is shown that the algorithm accurately predicts heating rates, and when compared with the single grid algorithm computes solutions in one-third the computational time.

Greene, Francis A.↗
NASA NTRS1992

The processing of the Viking Orbiter range data and its contribution to Mars gravity solutions

The processing of Doppler data has been the primary method for deriving models of the Mars gravity field. Since the Mariner 9 and Viking spacecraft were placed in orbit about Mars, many models from degree and order 6 to degree and order 50 have been developed. However, during the Viking mission, some 26,000 range measurements to the two Viking Orbiters were also obtained. These data have not previously been used in the derivation of Mars gravity models. A portion of these range data have been processed simultaneously with the Doppler data. Normal equations were generated for both sets of data and were used to create two solutions complete to degree and order 30: a nominal solution including both the range and the Doppler data (MGM-R100), and another solution including only the Doppler data (MGM-R101). Tests with the covariances of these solutions, as well as with orbit overlap tests indicate that the interplanetary range data can be used to improve the modeling of the Mars gravity field.

Lemoine, Frank G.↗
NASA NTRS1993

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme-based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗
NASA NTRS1993

Pinching solutions of slender cylindrical jets

Simplified equations for slender jets are derived for a circular jet of one fluid flowing into an ambient second fluid, the flow being confined in a circular tank. Inviscid flows are studied which include both surface tension effects and Kelvin-Helmholtz instability. For slender jets a coupled nonlinear system of equations is found for the jet shape and the axial velocity jump across it. The equations can break down after a finite time and similarity solutions are constructed, and studied analytically and numerically. The break-ups found pertain to the jet pinching after a finite time, without violation of the slender jet ansatz. The system is conservative and admissible singular solutions are those which conserve the total energy, mass, and momentum. Such solutions are constructed analytically and numerically, and in the case of vortex sheets with no surface tension certain solutions are given in closed form.

Papageorgiou, Demetrios T.↗
NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗
NASA NTRS1993

Projection techniques for iterative solution of A(bar x) equals (bar b) with successive right-hand sides

Two projection techniques for computing approximate solutions to linear systems of the form A(bar x)(sup n) = (bar b)(sup n), for a sequence n = 1, 2, ..., e.g., such as arises from time discretization of a partial differential equation, are presented. The inexpensive approximate solutions can be used as initial guesses for iterative solution of the system, resulting in significantly reduced computational expense. Examples of two- and three-dimensional incompressible Navier-Stokes calculations are presented in which x represents the pressure, and A is a discrete Poisson operator. In flows containing significant dynamic activity, these projection techniques lead to as much as a two-fold reduction in solution time.

Fischer, Paul F.↗
NASA NTRS1994

Improved solution accuracy for Landsat-4 (TDRSS-user) orbit determination

This paper presents the results of a study to compare the orbit determination accuracy for a Tracking and Data Relay Satellite System (TDRSS) user spacecraft, Landsat-4, obtained using a Prototype Filter Smoother (PFS), with the accuracy of an established batch-least-squares system, the Goddard Trajectory Determination System (GTDS). The results of Landsat-4 orbit determination will provide useful experience for the Earth Observing System (EOS) series of satellites. The Landsat-4 ephemerides were estimated for the January 17-23, 1991, timeframe, during which intensive TDRSS tracking data for Landsat-4 were available. Independent assessments were made of the consistencies (overlap comparisons for the batch case and convariances for the sequential case) of solutions produced by the batch and sequential methods. The filtered and smoothed PFS orbit solutions were compared with the definitive GTDS orbit solutions for Landsat-4; the solution differences were generally less than 15 meters.

Oza, D. H.↗
NASA NTRS1994

Open solutions to distributed control in ground tracking stations

The advent of high speed local area networks has made it possible to interconnect small, powerful computers to function together as a single large computer. Today, distributed computer systems are the new paradigm for large scale computing systems. However, the communications provided by the local area network is only one part of the solution. The services and protocols used by the application programs to communicate across the network are as indispensable as the local area network. And the selection of services and protocols that do not match the system requirements will limit the capabilities, performance, and expansion of the system. Proprietary solutions are available but are usually limited to a select set of equipment. However, there are two solutions based on 'open' standards. The question that must be answered is 'which one is the best one for my job?' This paper examines a model for tracking stations and their requirements for interprocessor communications in the next century. The model and requirements are matched with the model and services provided by the five different software architectures and supporting protocol solutions. Several key services are examined in detail to determine which services and protocols most closely match the requirements for the tracking station environment. The study reveals that the protocols are tailored to the problem domains for which they were originally designed. Further, the study reveals that the process control model is the closest match to the tracking station model.

Heuser, William Randy↗
NASA NTRS1995

Validation of the SINDA/FLUINT code using several analytical solutions

The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Keller, John R.↗
NASA NTRS1995

Developmental problems and their solution for the Space Shuttle main engine alternate liquid oxygen high-pressure turbopump: Anomaly or failure investigation the key

The Space Shuttle main engine (SSME) alternate high-pressure liquid oxygen pump experienced synchronous vibration and ball bearing life problems that were program threatening. The success of the program hinged on the ability to solve these development problems. The design and solutions to these problems are engirded in the lessons learned and experiences from prior programs, technology programs, and the ability to properly conduct failure or anomaly investigations. The failure investigation determines the problem cause and is the basis for recommending design solutions. For a complex problem, a comprehensive solution requires that formal investigation procedures be used, including fault trees, resolution logic, and action items worked through a concurrent engineering-multidiscipline team. The normal tendency to use an intuitive, cut-and-try approach will usually prove to be costly, both in money and time and will reach a less than optimum, poorly understood answer. The SSME alternate high-pressure oxidizer turbopump development has had two complex problems critical to program success: (1) high synchronous vibrations and (2) excessive ball bearing wear. This paper will use these two problems as examples of this formal failure investigation approach. The results of the team's investigation provides insight into the complexity of the turbomachinery technical discipline interacting/sensitivities and the fine balance of competing investigations required to solve problems and guarantee program success. It is very important to the solution process that maximum use be made of the resources that both the contractor and Government can bring to the problem in a supporting and noncompeting way. There is no place for the not-invented-here attitude. The resources include, but are not limited to: (1) specially skilled professionals; (2) supporting technologies; (3) computational codes and capabilities; and (4) test and manufacturing facilities.

Ryan, R.↗
NASA NTRS1995

Solution of the three-dimensional Helmholtz equation with nonlocal boundary conditions

The Helmholtz equation is solved within a three-dimensional rectangular duct with a nonlocal radiation boundary condition at the duct exit plane. This condition accurately models the acoustic admittance at an arbitrarily-located computational boundary plane. A linear system of equations is constructed with second-order central differences for the Helmholtz operator and second-order backward differences for both local admittance conditions and the gradient term in the nonlocal radiation boundary condition. The resulting matrix equation is large, sparse, and non-Hermitian. The size and structure of the matrix makes direct solution techniques impractical; as a result, a nonstationary iterative technique is used for its solution. The theory behind the nonstationary technique is reviewed, and numerical results are presented for radiation from both a point source and a planar acoustic source. The solutions with the nonlocal boundary conditions are invariant to the location of the computational boundary, and the same nonlocal conditions are valid for all solutions. The nonlocal conditions thus provide a means of minimizing the size of three-dimensional computational domains.

Hodge, Steve L.↗
NASA NTRS1994

Unique determination of solutions to the Burnett equations

Recent success in applying the Burnett equations to the one-dimensional shock-structure problem has raised the issue of whether the full Burnett equations can be used replace the Navier-Stokes equations for solving boundary-value problems in rarefied gasdynamics. As is familiar from the classical rarefied gasdynamics literature, the Burnett equations, if not solved as a successive approximation to the Navier-Stokes equations for a small Knudsen number, would require more boundary conditions than those in the Navier-Stokes system, owing to the presence of the higher-order derivatives. In this paper, this issue is examined with concrete solution examples for the steady Couette flows, addressing specifically whether solutions to the full Burnett equations can be uniquely determined without adding more boundary conditions than those in the Navier-Stokes system. The analysis, supported by detailed numerical solutions, confirms that additional boundary conditions are needed as long as the Knudsen number is not identically zero, lest the solution to the Burnett equations is not unique.

Lee, C. J.↗
NASA NTRS1995

Dissipative MHD solutions for resonant Alfven waves in 1-dimensional magnetic flux tubes

The present paper extends the analysis by Sakurai, Goossens, and Hollweg (1991) on resonant Alfven waves in nonuniform magnetic flux tubes. It proves that the fundamental conservation law for resonant Alfven waves found in ideal MHD by Sakurai, Goossens, and Hollweg remains valid in dissipative MHD. This guarantees that the jump conditions of Sakurai, Goossens, and Hollweg, that connect the ideal MHD solutions for xi(sub r), and P' across the dissipative layer, are correct. In addition, the present paper replaces the complicated dissipative MHD solutions obtained by Sakurai, Goossens, and Hollweg for xi(sub r), and P' in terms of double integrals of Hankel functions of complex argument of order 1/3 with compact analytical solutions that allow a straight- forward mathematical and physical interpretation. Finally, it presents an analytical dissipative MHD solution for the component of the Lagrangian displacement in the magnetic surfaces perpen- dicular to the magnetic field lines xi(sub perpendicular) which enables us to determine the dominant dynamics of resonant Alfven waves in dissipative MHD.

Goossens, Marcel↗
NASA NTRS1995

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - 2

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Yee, H. C.↗
NASA NTRS1996

Asymptotic solution of the turbulent mixing layer for velocity ratio close to unity

The equations describing the first two terms of an asymptotic expansion of the solution of the planar turbulent mixing layer for values of the velocity ratio close to one are obtained. The first term of this expansion is the solution of the well-known time-evolving problem and the second, which includes the effects of the increase of the turbulence scales in the stream-wise direction, obeys a linear system of equations. Numerical solutions of these equations for a two-dimensional reacting mixing layer show that the correction to the time-evolving solution may explain the asymmetry of the entrainment and the differences in product generation observed in flip experiments.

Higuera, F. J.↗