Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Numerical solution for supersonic flow near the trailing edge of a flat plate

In the present study, problems of laminar and turbulent two-dimensional flow of a viscous compressible fluid near the trailing edge of a thin flat plate are considered. The complete set of Navier-Stokes equations is solved by the finite-difference method of MacCormack (MacCormack and Baldwin, 1975). It is an explicit, predictor-corrector, time-splitting method of second order acuracy. The computational mesh employed has sufficient resolution for all the characteristic lengths suggested by theory. In the laminar case, the present results are compared with the triple deck solution of Daniels (1974). This comparison indicates that the asymptotic triple deck theory for supersonic trailing edge flow is accurate within five percent for Reynolds numbers greater than 1000. In the turbulent case, the Prandtl-Van Driest-Clauser algebraic eddy viscosity model is used. The numerical results show that the region of upstream influence is approximately of the order of the boundary layer thickness. The solutions for skin-friction, pressure and wake center-line velocity are presented.

Degani, D.↗
NASA NTRS1979

The solution of the three-dimensional viscous-compressible Navier-Stokes equations on a vector computer

The development of a vectorized computer code for the solution of the three-dimensional viscous-compressible Navier-Stokes equations is described. The code is applied on the CDC STAR-100 vector computer which is capable of achieving high result rates when a high degree of parallelism is present in the computations. The computational technique is an explicit time-split MacCormack predictor-corrector algorithm. Since a large volume of data is processed and virtual memory utilized, a data management scheme based on interleaving is used. The program has been applied to obtain the solution of the laminar supersonic flow about a family of three-dimensional corners. The equations of motion are expressed in a generalized form relative to a uniform rectangular computational domain. The metric coefficient and boundary conditions must be supplied for the corresponding physical domain. For calculations with 30,000 grid points, a computational rate of 0.00015 seconds per grid point per time step is observed.

Smith, R. E.↗
NASA NTRS1992

Diffraction by a symmetric material junction simulated with generalized sheet transistion conditions

When diffraction by a material discontinuity in a thick dielectric/ferrite slab is treated for the dual integral equation approach (a variation of the Weiner-Hopf method), the resulting solution displays unknown constants, characteristic of the incompleteness associated with the generalized transition conditions employed in modeling a thick slab. The constants can be shown to depend on the geometry and properties of the discontinuity. Their explicit determination by enforcing field continuity across the slab junction is treated. The field internal to the slab is determined by first expressing the field as an expansion of discrete and continuous eigenmodes valid in the exterior and interior slab region. The expansion or mode coefficients are then identified by recasting the Weiner-Hopf solution in a suitable form, enabling the solution constants to be numerically computed by matching the field expansions to the left and right of the slab discontinuity at a discrete set of points.

Ricoy, Mark A.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. The governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models are described in detail.

Towne, Charles E.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗
NASA NTRS1962

BUCKLING OF INITIALLY IMPERFECT AXIALLY COMPRESSED CYLINDRICAL SHELLS

The effects of imperfections on the buckling strength of pressurized cylinders under axial compression are studied. The imperfection factor is considered in the finite-deflection compatibility equation, as well as in the equilibrium equation. A new relation to express the imperfection as an explicit exponential function of pressure and radius thickness ratio is proposed. A method for finding the exponential function is described, and the solution for the critical stress is found in a fairly simple form. The relation of decrease of stability with respect to magnitude of imperfections if found.

CYLINDRICAL SHELL↗
NASA NTRS1979

Space Shuttle ascent guidance, navigation, and control

The factors leading to the particular design of the Shuttle guidance, navigation and control software are discussed. The derivation of explicit guidance equations satisfying a wide range of different maneuver constraints and steering equations that create attitude steering errors from the guidance solutions is presented, as are navigation equations, and equations for identifying faulty instruments from sets of redundant instrument measurements.

R. L. Mchenry↗
NASA NTRS1979

Dynamic acoustics for the STAR-100

An algorithm is described to compute time dependent acoustic waves in a jet. The method differs from previous methods in that no harmonic time dependence is assumed, thus permitting the study of nonharmonic acoustical behavior. Large grids are required to resolve the acoustic waves. Since the problem is nonstiff, explicit high order schemes can be used. These have been adapted to the STAR-100 with great efficiencies and permitted the efficient solution of problems which would not be feasible on a scalar machine.

Bayliss, A.↗
NASA NTRS1981

Role of the doubly stochastic Neyman type-A and Thomas counting distributions in photon detection

The sums of Meyman type-A and Thomas random variables are shown to retain their form under the constant multiplication parameter constraint. The conditions under which the two random variables converge in distribution to the fixed multiplicative Poisson, and to the Gaussian, are presented. It is shown that the latter result is important, in that it provides a solution to likelihood-ratio detection, estimation and discrimination problems in the presence of many kinds of signal and noise. Among the explicit applications presented are: (1) the photon-counter scintillation detection of nuclear particles when particle flux is low, (2) the photon-counting detection of weak optical signals in the presence of ionizing radiation, and (3) the design of a star-scanner spacecraft guidance system adequate for hostile space environments.

Teich, M. C.↗
NASA NTRS1984

Dual approximations in optimal control

A dual approximation for the solution to an optimal control problem is analyzed. The differential equation is handled with a Lagrange multiplier while other constraints are treated explicitly. An algorithm for solving the dual problem is presented.

Hager, W. W.↗
NASA NTRS1987

Use of an implicit formulation based on quasilinearization for the aeroelastic response and stability of rotor blades in forward flight

This paper describes a new methodology for the formulation of the aeroelastic stability and response problem for helicopter rotor blades. The mathematical expressions for the aerodynamic loads need not be explicit functions of the blade displacement quantities. This methodology is combined with a finite element model of the blade, and a quasilinearization solution technique. The resulting computer program is used to study the behavior of blades with noncoincident elastic axis, aerodynamic centers, and centers of mass.

Celi, R.↗
NASA NTRS1991

Extension of multigrid methodology to supersonic/hypersonic 3-D viscous flows

A multigrid acceleration technique developed for solving 3-D Navier-Stokes equations for subsonic/transonic flows was extended to supersonic/hypersonic flows. An explicit multistage Runge-Kutta type of time stepping scheme is used as the basic algorithm in conjunction with the multigrid scheme. Solutions were obtained for a blunt conical frustum at Mach 6 to demonstrate the applicability of the multigrid scheme to high speed flows. Computations were performed for a generic High Speed Civil Transport configuration designed to cruise at Mach 3. These solutions show both the efficiency and accuracy of the present scheme for computing high speed viscous flows over configurations of practical interest.

Vatsa, Veer N.↗
NASA NTRS1992

Finite-element solution scheme for compressible viscous flow

Adaptivity offers the possibility of using FEM for problems that are otherwise intractable. The solution algorithm is combined with an adaptive approach in which the unstructured portion of the grid is adaptively regenerated as the solution proceeds; the structured portion of the mesh is modified accordingly, to ensure both compatibility along the interface between the implicit and explicit element groups and increased resolution of the viscous layer. The approach offers significant computational savings in 3D simulations.

Hassan, O.↗
NASA NTRS1994

Hypersonic vehicle model and control law development using H(infinity) and micron synthesis

The control system design for a Single Stage To Orbit (SSTO) air breathing vehicle will be central to a successful mission because a precise ascent trajectory will preserve narrow payload margins. The air breathing propulsion system requires the vehicle to fly roughly halfway around the Earth through atmospheric turbulence. The turbulence, the high sensitivity of the propulsion system to inlet flow conditions, the relatively large uncertainty of the parameters characterizing the vehicle, and continuous acceleration make the problem especially challenging. Adequate stability margins must be provided without sacrificing payload mass since payload margins are critical. Therefore, a multivariable control theory capable of explicitly including both uncertainty and performance is needed. The H(infinity) controller in general provides good robustness but can result in conservative solutions for practical problems involving structured uncertainty. Structured singular value mu framework for analysis and synthesis is potentially much less conservative and hence more appropriate for problems with tight margins. An SSTO control system requires: highly accurate tracking of velocity and altitude commands while limiting angle-of-attack oscillations, minimized control power usage, and a stabilized vehicle when atmospheric turbulence and system uncertainty are present. The controller designs using H(infinity) and mu-synthesis procedures were compared. An integrated flight/propulsion dynamic mathematical model of a conical accelerator vehicle was linearized as the vehicle accelerated through Mach 8. Vehicle acceleration through the selected flight condition gives rise to parametric variation that was modeled as a structured uncertainty. The mu-analysis approach was used in the frequency domain to conduct controller analysis and was confirmed by time history plots. Results demonstrate the inherent advantages of the mu framework for this class of problems.

Gregory, Irene M.↗
NASA NTRS1998

Optimality of a Fully Stressed Design

For a truss a fully stressed state is reached and when all its members are utilized to their full strength capacity. Historically, engineers considered such a design optimum. But recently this optimality has been questioned, especially since the weight of the structure is not explicitly used in fully stressed design calculations. This paper examines optimality of the full stressed design (FSD) with analytical and graphical illustrations. Solutions for a set of examples obtained by using the FSD method and optimization methods numerically confirm the optimality of the FSD. The FSD, which can be obtained with a small amount of calculation, can be extended to displacement constraints and to nontruss-type structures.

Patnaik, Surya N.↗
NASA NTRS2005

Assessment of an Explicit Algebraic Reynolds Stress Model

This study assesses an explicit algebraic Reynolds stress turbulence model in the in the three-dimensional Reynolds averaged Navier-Stokes (RANS) solver, ISAAC (Integrated Solution Algorithm for Arbitrary Con gurations). Additionally, it compares solutions for two select configurations between ISAAC and the RANS solver PAB3D. This study compares with either direct numerical simulation data, experimental data, or empirical models for several different geometries with compressible, separated, and high Reynolds number flows. In general, the turbulence model matched data or followed experimental trends well, and for the selected configurations, the computational results of ISAAC closely matched those of PAB3D using the same turbulence model.

Carlson, Jan-Renee↗
NASA NTRS2007

Estimating the Inertia Matrix of a Spacecraft

A paper presents a method of utilizing some flight data, aboard a spacecraft that includes reaction wheels for attitude control, to estimate the inertia matrix of the spacecraft. The required data are digitized samples of (1) the spacecraft attitude in an inertial reference frame as measured, for example, by use of a star tracker and (2) speeds of rotation of the reaction wheels, the moments of inertia of which are deemed to be known. Starting from the classical equations for conservation of angular momentum of a rigid body, the inertia-matrix-estimation problem is formulated as a constrained least-squares minimization problem with explicit bounds on the inertia matrix incorporated as linear matrix inequalities. The explicit bounds reflect physical bounds on the inertia matrix and reduce the volume of data that must be processed to obtain a solution. The resulting minimization problem is a semidefinite optimization problem that can be solved efficiently, with guaranteed convergence to the global optimum, by use of readily available algorithms. In a test case involving a model attitude platform rotating on an air bearing, it is shown that, relative to a prior method, the present method produces better estimates from few data.

Acikmese, Behcet↗
NASA NTRS1999

Solving Upwind-Biased Discretizations: Multigrid Solver Using Semicoarsening - 2

This paper studies a novel multigrid approach to the solution for a second order upwind biased discretization of the convection equation in two dimensions. This approach is based on semi-coarsening and well balanced explicit correction terms added to coarse-grid operators to maintain on coarse-grid the same cross-characteristic interaction as on the target (fine) grid. Colored relaxation schemes are used on all the levels allowing a very efficient parallel implementation. The results of the numerical tests can be summarized as follows: 1) The residual asymptotic convergence rate of the proposed V(0, 2) multigrid cycle is about 3 per cycle. This convergence rate far surpasses the theoretical limit (4/3) predicted for standard multigrid algorithms using full coarsening. The reported efficiency does not deteriorate with increasing the cycle, depth (number of levels) and/or refining the target-grid mesh spacing. 2) The full multi-grid algorithm (FMG) with two V(0, 2) cycles on the target grid and just one V(0, 2) cycle on all the coarse grids always provides an approximate solution with the algebraic error less than the discretization error. Estimates of the total work in the FMG algorithm are ranged between 18 and 30 minimal work units (depending on the target (discretizatioin). Thus, the overall efficiency of the FMG solver closely approaches (if does not achieve) the goal of the textbook multigrid efficiency. 3) A novel approach to deriving a discrete solution approximating the true continuous solution with a relative accuracy given in advance is developed. An adaptive multigrid algorithm (AMA) using comparison of the solutions on two successive target grids to estimate the accuracy of the current target-grid solution is defined. A desired relative accuracy is accepted as an input parameter. The final target grid on which this accuracy can be achieved is chosen automatically in the solution process. the actual relative accuracy of the discrete solution approximation obtained by AMA is always better than the required accuracy; the computational complexity of the AMA algorithm is (nearly) optimal (comparable with the complexity of the FMG algorithm applied to solve the problem on the optimally spaced target grid).

Diskin, Boris↗