SEARCH · Search NASA
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.
An explicit closed-form solution is obtained for the compensator transfer function for the compensator designed by Balakrishnan (1991) for stability enhancement of flexible structures with collocated controllers. The solution is positive real but nonrational and yields high stability and robustness.
Explicit solution for Lagrange multipliers in optimal thrust programming and steering problem
A derivation of an explicit solution to the two point boundary-value problem of exoatmospheric guidance and trajectory optimization is presented. Fixed initial conditions and continuous burn, multistage thrusting are assumed. Any number of end conditions from one to six (throttling is required in the case of six) can be satisfied in an explicit and practically optimal manner. The explicit equations converge for off nominal conditions such as engine failure, abort, target switch, etc. The self starting, predictor/corrector solution involves no Newton-Rhapson iterations, numerical integration, or first guess values, and converges rapidly if physically possible. A form of this algorithm has been chosen for onboard guidance, as well as real time and preflight ground targeting and trajectory shaping for the NASA Space Shuttle Program.
An explicit solution of the spectral radiation leaving an arbitrary point on the wall of a spherical cavity with diffuse reflectivity is obtained. A general measurement equation is obtained that describes the output of a sensor measuring the sphere output within a given field of view and with specified angular and spectral responses. The results are applied to the Earth Radiation Budget Experiment (ERBE) integrating sphere. It is shown that by using appropriate interpretation and processing, a high-accuracy shortwave calibration of the ERBE sensors can be achieved.
We present a class of compensators in explicit form (not requiring numerical computer calculations) for stabilizing flexible structures with collocated rate sensors. They are based on the explicit solution, valid for both Continuum and FEM Models, of the LQG problem for minimizing mean square rate. They are robust with respect to system stability (will not destabilize modes even with mismatch of parameters), can be instrumented in state space form suitable for digital controllers, and can be specified directly from the structure modes and mode 'signature' (displacement vectors at sensor locations). Some simulation results are presented for the NASA LaRC Phase-Zero Evolutionary Model - a modal Trust model with 86 modes - showing damping ratios attainable as a function of compensator design parameters and complexity.
The matrix of the normal equations of a central overlapping-plate solution with general models for frame parameters and the star parameters is inverted explicitly in terms of the inverses of the matrices in the isolated systems for the frame and star parameters. This makes it possible to calculate from explicit formulas the numerical estimates of the star and frame parameters which would follow from a rigorous overlapping-plate solution. Also, the covariance matrix of the adjustment parameters are calculated.
Studies aimed at an increase in the efficiency of calculating transient temperature fields in complex aerospace vehicle structures are reported. The advantages and disadvantages of explicit and implicit algorithms are discussed and a promising set of implicit algorithms with variable time steps, known as GEARIB, is described. Test problems, used for evaluating and comparing various algorithms, are discussed and finite element models of the configurations are described. These problems include a coarse model of the Space Shuttle wing, an insulated frame tst article, a metallic panel for a thermal protection system, and detailed models of sections of the Space Shuttle wing. Results generally indicate a preference for implicit over explicit algorithms for transient structural heat transfer problems when the governing equations are stiff (typical of many practical problems such as insulated metal structures). The effects on algorithm performance of different models of an insulated cylinder are demonstrated. The stiffness of the problem is highly sensitive to modeling details and careful modeling can reduce the stiffness of the equations to the extent that explicit methods may become the best choice. Preliminary applications of a mixed implicit-explicit algorithm and operator splitting techniques for speeding up the solution of the algebraic equations are also described.
Optimal thrust programming of single stage rocket vehicle for maximum payoff at end conditions using Lagrange multipliers for singular arc
A second-order time-accurate and spatially factored algorithm was used in a finite difference scheme for the numerical solution of the time-dependent, incompressible, two dimensional Navier-Stokes equations in conservation-law form using vorticity and stream function variables. The systems of equations are solved at each time step by an iterative technique. Numerical results were obtained for a circular cylinder at a Reynolds number of 15, and an NACA 0012 airfoil at zero angle of attack at Reynolds numbers of 10 to the third and 10 to the fourth powers. The results are in agreement with another numerical technique, and the computing time required to obtain the steady state solution at the Reynolds number of 10 to the 4th power was 49.7 sec on CDC 7600 computer using a 65 x 84 computational grind.
A Fredholm integral equation is solved in closed form by reducing it to a standard boundary value problem for a sectionally analytic function. The procedure maybe extended to more general classes of functions whose Fourier transforms only exist along certain lines in the complex k-plane. Examples are used to illustrate the results.
Modern military aircraft transparency systems, windshields and canopies, are complex systems which must meet a large and rapidly growing number of requirements. Many of these transparency system requirements are conflicting, presenting difficult balances which must be achieved. One example of a challenging requirements balance or trade is shaping for stealth versus aircrew vision. The large number of requirements involved may be grouped in a variety of areas including man-machine interface; structural integration with the airframe; combat hazards; environmental exposures; and supportability. Some individual requirements by themselves pose very difficult, severely nonlinear analysis problems. One such complex problem is that associated with the dynamic structural response resulting from high energy bird impact. An improved analytical capability for soft-body impact simulation was developed.
A general mathematical description of the response behavior of thin-skin pneumatic envelopes and many other membrane and cable structures produces under-constrained systems that pose severe difficulties to analysis. These systems are mobile, and the general mathematical description exposes the mobility. Yet the response behavior of special under-constrained structures under special loadings can be accurately predicted using a constrained mathematical description. The static response behavior of systems that are infinitesimally mobile, such as a non-slack membrane subtended from a rigid or elastic boundary frame, can be easily analyzed using such general mathematical description as afforded by the non-linear, finite element method using an implicit solution scheme if the incremental uploading is guided through a suitable path. Similarly, if such structures are assembled with structural lack of fit that provides suitable self-stress, then dynamic response behavior can be predicted by the non-linear, finite element method and an implicit solution scheme. An explicit solution scheme is available for evolution problems. Such scheme can be used via the method of dynamic relaxation to obtain the solution to a static problem. In some sense, pneumatic envelopes and many other compliant structures can be said to have destiny under a specified loading system. What that means to the analyst is that what happens on the evolution path of the solution is irrelevant as long as equilibrium is achieved at destiny under full load and that the equilibrium is stable in the vicinity of that load. The purpose of this paper is to alert practitioners to the fact that non-standard procedures in finite element analysis are useful and can be legitimate although they burden their users with the requirement to use special caution. Some interesting findings that are useful to the US Scientific Balloon Program and that could not be obtained without non-standard techniques are presented.
An optimal trajectory planning problem for a single-link, flexible joint manipulator is studied. A global feedback-linearization is first applied to formulate the nonlinear inequality-constrained optimization problem in a suitable way. Then, an exact and explicit structural formula for the optimal solution of the problem is derived and the solution is shown to be unique. It turns out that the optimal trajectory planning and control can be done off-line, so that the proposed method is applicable to both theoretical analysis and real time tele-robotics control engineering.
Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.
The explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid moving in any actual atmosphere was published several years ago. The solution yields the theoretical relation of l, the distance flown by the meteoroid in its trajectory, as a function of time, t, assuming that the height, h, is a known function of l. The photographic records of meteors and fireballs are coded by time marks, using a rotating shutter or a similar device to break the moving image. Time is, thus, the independent variable and for each time mark on a meteoroid trajector, the observed distance along the trajectory, l sub obs, as well as the double- or multiple- station photographs of the same meteoroid. Applying this solution to all available Prairie Network (PN) fireball-records, we recognized that the majority of them gave good solutions with standard deviations somewhat bigger than the intrinsic geometrical precision of the data. We also noticed that, on an average, previous methods of evaluation of the meteoroid velocities (interpolation polynomials, numerical differenciation of the observed l sub obs) used up to only several tens of percent of the intrinsic precision of the PN observational data. When residuals of these solutions, i.e. l sub obs - l sub com, were represented as a function of time for about 75 percent of solutions. The remaining 25 percent of residuals showed systematic changes with time exceeding one standard deviation. We tried to explain these systematic time course of residuals by using different meteoroids first computed theoretically and then analyzed by the same model as the natural PN fireballs were. The conclusion of these model computations is that systematic time changes of residuals in the nonfragmenting model exceeding one standard of deviation are caused by sudden gross fragmentation at one or more trajectory points. Thus, we generalized the explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid by allowing for one or more points, where sudden gross fragmentation can occur. Using this generalized solution, the distances along the meteoroid trajectory can be computed for any choice of input parameters and compared with the observed distances flown by the meteoroid. For the most precise and long fireball trajectories, the least-squares solution can, thus, yield the initial velocities, the ablation coefficients, the dynamical masses, the positions of gross-fragmentation points, and the terminal mass. At a gross-fragmentation point, the ratio of the main mass to all the remaining fragments can be compared with the dynamic mass determined from our gross-fragmentation model and, thus, the meteoroid bulk density can be evaluated. We applied the gross-fragmentation model to sever PN fireballs showing time changes of residuals, and we recognized that, in all these cases, the new computed bulk densities of meteoroids resulted higher in comparison with the meteoroid densities determined with the non-gross-fragmentation model. Other aspects of the study are discussed.
An algorithm for solving the nonlinear stationary Navier-Stokes problem is developed. Explicit error estimates are given. This mathematical technique is potentially adaptable to the separation problem.
The bifurcation diagram associated with the logistic equation v sup n+1 = av sup n (1-v sup n) is by now well known, as is its equivalence to solving the ordinary differential equation u prime = alpha u (1-u) by the explicit Euler difference scheme. It has also been noted by Iserles that other popular difference schemes may not only exhibit period doubling and chaotic phenomena but also possess spurious fixed points. Runge-Kutta schemes applied to both the equation u prime = alpha u (1-u) and the cubic equation u prime = alpha u (1-u)(b-u) were studied computationally and analytically and their behavior was contrasted with the explicit Euler scheme. Their spurious fixed points and periodic orbits were noted. In particular, it was observed that these may appear below the linearized stability limits of the scheme and, consequently, computation may lead to erroneous results.
The objective here is to describe the lightning hazards to buildings and their internal environments using advanced formulations of Maxwell's Equations. The method described is the Three Dimensional Finite Difference Time Domain Solution. It can be used to solve for the lightning interaction with such structures in three dimensions with the inclusion of a considerable amount of detail. Special techniques were developed for including wire, plumbing, and rebar into the model. Some buildings have provisions for lightning protection in the form of air terminals connected to a ground counterpoise system. It is shown that fields and currents within these structures can be significantly high during a lightning strike. Time lapse video presentations were made showing the electric and magnetic field distributions on selected cross sections of the buildings during a simulated lightning strike.