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.
A method is described for 3D-gas dynamics computer simulation in regions of complicated shape by means of nonadjusted rectangular grids providing unified treatment of various problems. Some test problem computation results are given.
Optimum sensitivity is defined as the derivative of the optimum design with respect to some problem parameter, P. The problem parameter is usually fixed during optimization, but may be changed later. Thus, optimum sensitivity is used to estimate the effect of changes in loads, materials or constraint bounds on the design without expensive re-optimization. Here, the general topic of optimum sensitivity is discussed, available methods identified, examples given, and the difficulties encountered in calculating this information in nonlinear constrained optimization are identified.
The development of efficient techniques for calculating sensitivity derivatives is studied. The objective is to present a computational procedure for calculating sensitivity derivatives as part of performing structural reanalysis for large-scale problems. The scope is limited to framed type structures. Both linear static analysis and free-vibration eigenvalue problems are considered.
The formulation and solution of inverse problems for the estimation of parameters which describe damping and other dynamic properties in distributed models for the vibration of flexible structures is considered. Motivated by a slewing beam experiment, the identification of a nonlinear velocity dependent term which models air drag damping in the Euler-Bernoulli equation is investigated. Galerkin techniques are used to generate finite dimensional approximations. Convergence estimates and numerical results are given. The modeling of, and related inverse problems for the dynamics of a high pressure hose line feeding a gas thruster actuator at the tip of a cantilevered beam are then considered. Approximation and convergence are discussed and numerical results involving experimental data are presented.
The first rotordynamics workshop proceedings (NASA CP-2133, 1980) emphasized a feeling of uncertainty in predicting the stability of characteristics of high-performance turbomachinery. In the second workshop proceedings (NASA CP-2250, 1982) these uncertainities were reduced through programs established to systematically resolve problems, with emphasis on experimental validiation of the forces that influence rotordynamics. In third proceedings (NASA CP-2338, 1984) many programs for predicting or measuring forces and force coefficients in high-performance turbomachinery produced results. Data became available for designing new machines with enhanced stability characteristics or for upgrading existing machines. The present workshop proceedings illustrates a continued trend toward a more unified view of rotordynamic instability problems and several encouraging new analytical developments.
A description of the finite element implementation of Robinson's unified viscoplastic model into the General Purpose Finite Element Program (MARC) is presented. To demonstrate its application, the implementation is applied to some uniaxial and multiaxial problems. A comparison of the results for the multiaxial problem of a thick internally pressurized cylinder, obtained using the finite element implementation and an analytical solution, is also presented. The excellent agreement obtained confirms the correct finite element implementation of Robinson's model.
A new algorithm for solving subset regression problems is described. The algorithm performs a QR decomposition with a new column-pivoting strategy, which permits subset selection directly from the originally defined regression parameters. This, in combination with a number of extensions of the new technique, makes the method a very flexible tool for analyzing subset regression problems in which the parameters have a physical meaning.
The problem of finding extraterrestrial civilizations and establishing contact with them is directly related to the problem of their evolution. Possible patterns in this evolution and the stages in the evolution of extraterrestrial civilizations are examined.
A method is proposed for the spectral element simulation of incompressible flow. This method constitutes in a well-posed optimal approximation of the steady Stokes problem with no spurious modes in the pressure. The resulting method is analyzed, and numerical results are presented for a model problem.
Element-by-element approximate factorization, implicit-explicit and adaptive implicit-explicit approximation procedures are presented for the finite-element formulations of large-scale fluid dynamics problems. The element-by-element approximation scheme totally eliminates the need for formation, storage and inversion of large global matrices. Implicit-explicit schemes, which are approximations to implicit schemes, substantially reduce the computational burden associated with large global matrices. In the adaptive implicit-explicit scheme, the implicit elements are selected dynamically based on element level stability and accuracy considerations. This scheme provides implicit refinement where it is needed. The methods are applied to various problems governed by the convection-diffusion and incompressible Navier-Stokes equations. In all cases studied, the results obtained are indistinguishable from those obtained by the implicit formulations.
The existence of a compensator which simultaneously renders a given r-tuple of multiinput-multioutput p x m linear dynamical systems internally stable is investigated. In particular, a set of simultaneously stabilizable r-tuples of plants is parametrized, and it is shown that, provided r = max(m,p) or less, the above set is semialgebraic and dense in the space Sigma of r-tuples of plants. An extension of the classical pole placement and stabilization problems is considered, and the simultaneous partial pole placement problem is investigated.
Adaptive finite element methods for steady and unsteady flow problems in two-dimensional domains are described. Details of a data management scheme are given that provide for the rapid implementation of various CFD algorithms on changing unstructured meshes. The results of several numerical experiments on subsonic and supersonic flow problems are discussed.
The inverse problem of satellite temperature profile retrieval is well known to be ill-posed. This means that not only is a vertical temperature profile solution not unique, but that two solutions can be very different from each other. A set of atmosphere-like, and true atmospheric examples of significantly dissimilar inverse solutions, were sought and found, using an 11-channel simulated HIRS sounding radiometer. Using the Riemann-Lebesgue Lemma for guidance, it is shown that simultaneous, numerical solutions of an atmospheric character may differ by as much as 10 K between 10-1000 mb. However, an empirical search for dissimilar solutions in the natural atmosphere reveals an extremely low probability of finding two significantly different RAOBs which produce radiance measurements whose differences cannot be resolved by the satellite radiometer. The empirical results are used to derive a first estimate of the limits of retrievability, analogous to the limits of predictability derivable from the ill-posed nature of the numerical weather prediction problem.
The paper surveys the last ten years of activity of the INME Swansea, dealing with problems of convection dominated flow. The basic explicit/implicit characteristic Galerkin process and its application to adaptive mesh refinement used in the solution of realistic problems is focused on.
The partitioning of a problem on a domain with unequal work estimates in different subdomains is considered in a way that balances the work load across multiple processors. Such a problem arises for example in solving partial differential equations using an adaptive method that places extra grid points in certain subregions of the domain. A binary decomposition of the domain is used to partition it into rectangles requiring equal computational effort. The communication costs of mapping this partitioning onto different microprocessors: a mesh-connected array, a tree machine and a hypercube is then studied. The communication cost expressions can be used to determine the optimal depth of the above partitioning.
This paper deals with nearly singular, possibly indefinite problems for which the usual multigrid solvers converge very slowly or even diverge. The main difficulty is related to some badly approximated smooth functions which correspond to eigenfunctions with nearly zero eigenvalues. A correction to the usual coarse-grid equations is derived, both in the correction scheme and in the full approximation scheme. The performance of the new algorithm using this correction is essentially as that of usual multigrid for definite problems.
Given a discretization stencil, partitioning the problem domain is an important first step for the efficient solution of partial differential equations on multiple processor systems. Partitions are derived that minimize interprocessor communication when the number of processors is known a priori and each domain partition is assigned to a different processor. This partitioning technique uses the stencil structure to select appropriate partition shapes. For square problem domains, it is shown that non-standard partitions (e.g., hexagons) are frequently preferable to the standard square partitions for a variety of commonly used stencils. This investigation is concluded with a formalization of the relationship between partition shape, stencil structure, and architecture, allowing selection of optimal partitions for a variety of parallel systems.
The formulation and recent applications of the finite-difference time-domain (FD-TD) method for the numerical modeling of electromagnetic scattering and interaction problems are considered. It is shown that improvements in FD-TD modeling concepts and software implementation often make it a preferable choice for structures which cannot be easily treated by conventional integral equations and asymptotic approaches. Recent FD-TD modeling validations in research areas including coupling to wires and wire bundles in free space and cavities, scattering from surfaces in relativistic motion, inverse scattering, and radiation condition theory, are reviewed. Finally, the advantages and disadvantages of FD-TD, and guidelines concerning when FD-TD should and should not be used in high-frequency electromagnetic modeling problems, are summarized.