Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

STARS: A general-purpose finite element computer program for analysis of engineering structures

STARS (Structural Analysis Routines) is primarily an interactive, graphics-oriented, finite-element computer program for analyzing the static, stability, free vibration, and dynamic responses of damped and undamped structures, including rotating systems. The element library consists of one-dimensional (1-D) line elements, two-dimensional (2-D) triangular and quadrilateral shell elements, and three-dimensional (3-D) tetrahedral and hexahedral solid elements. These elements enable the solution of structural problems that include truss, beam, space frame, plane, plate, shell, and solid structures, or any combination thereof. Zero, finite, and interdependent deflection boundary conditions can be implemented by the program. The associated dynamic response analysis capability provides for initial deformation and velocity inputs, whereas the transient excitation may be either forces or accelerations. An effective in-core or out-of-core solution strategy is automatically employed by the program, depending on the size of the problem. Data input may be at random within a data set, and the program offers certain automatic data-generation features. Input data are formatted as an optimal combination of free and fixed formats. Interactive graphics capabilities enable convenient display of nodal deformations, mode shapes, and element stresses.

Gupta, K. K.↗
NASA NTRS1975

Computer analysis of multicircuit shells of revolution by the field method

The field method, presented previously for the solution of even-order linear boundary value problems defined on one-dimensional open branch domains, is extended to boundary value problems defined on one-dimensional domains containing circuits. This method converts the boundary value problem into two successive numerically stable initial value problems, which may be solved by standard forward integration techniques. In addition, a new method for the treatment of singular boundary conditions is presented. This method, which amounts to a partial interchange of the roles of force and displacement variables, is problem independent with respect to both accuracy and speed of execution. This method was implemented in a computer program to calculate the static response of ring stiffened orthotropic multicircuit shells of revolution to asymmetric loads. Solutions are presented for sample problems which illustrate the accuracy and efficiency of the method.

Cohen, G. A.↗
NASA NTRS1974

COMOC: Three dimensional boundary region variant, programmer's manual

The three-dimensional boundary region variant of the COMOC computer program system solves the partial differential equation system governing certain three-dimensional flows of a viscous, heat conducting, multiple-species, compressible fluid including combustion. The solution is established in physical variables, using a finite element algorithm for the boundary value portion of the problem description in combination with an explicit marching technique for the initial value character. The computational lattice may be arbitrarily nonregular, and boundary condition constraints are readily applied. The theoretical foundation of the algorithm, a detailed description on the construction and operation of the program, and instructions on utilization of the many features of the code are presented.

Orzechowski, J. A.↗
NASA NTRS1987

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗
NASA NTRS1986

Turbulent relaxation processes in magnetohydrodynamics

Competing processes previously called 'selective decay' and 'dynamic alignment' are studied numerically for two-dimensional magnetohydrodynamic turbulence. In selective decay, the energy decays relatively to mean-square vector potential, and in dynamic alignment, the energy decays relatively to cross helicity. In the former case, the kinetic (fluid) energy decays to zero and the magnetic energy occupies the largest scales allowed by the boundary conditions. In the latter case, the velocity field and magnetic field become aligned and energetically equipartitioned. An extensive study of the initial value problem, with viscous and resistive dissipation, indicates that four distinct regimes of behavior are possible: a magnetically dominated regime, a velocity dominated regime, a dynamic alignment dominated regime, and a transition regime separating the others in parameter space. An analytical variational problem predicts several features seen in the computations, including geometrical alignment of velocity and magnetic fields, constraints on global quantities in the long-time limit, and in some cases, condensation to long wavelength modes.

Ting, Antonio C.↗
NASA NTRS1991

In-phase and out-of-phase axial-torsional fatigue behavior of Haynes 188 at 760 C

Isothermal, in-phase and out-of-phase axial-torsional fatigue experiments have been conducted at 760 C on uniform gage section, thin-walled tubular specimens of a wrought cobalt-base superalloy, Haynes 188. Test-control and data acquisition were accomplished with a minicomputer. Fatigue lives of the in- and out-of-phase axial-torsional fatigue tests have been estimated with four different multiaxial fatigue life prediction models that were developed primarly for predicting axial-torsional fatigue lives at room temperature. The models investigated were: (1) the von Mises equivalent strain range; (2) the Modified Multiaxiality Factor Approach; (3) the Modified Smith-Watson-Topper Parameter; and (4) the critical shear plane method of Fatemi, Socie, and Kurath. In general, life predictions by the von Mises equivalent strain range model were within a factor of 2 for a majority of the tests and the predictions by the Modified Multiaxiality Factor Approach were within a factor of 2, while predictions of the Modified Smith-Watson-Topper Parameter and of the critical shear plane method of Fatemi, Socie, and Kurath were unconservative and conservative, respectively, by up to factors of 4. In some of the specimens tested under combined axial-torsional loading conditions, fatigue cracks initiated near extensometer indentations. Two design modifications have been proposed to the thin-walled tubular specimen to overcome this problem.

Kalluri, Sreeramesh↗
NASA NTRS1994

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗
NASA NTRSDate not supplied

Comprehensive Analysis Aeroelastic Stability Predictions for the Mars Sample Recovery Helicopters and Ingenuity

A rotorcraft comprehensive analysis was used to generate aeroelastic stability predictions for both Ingenuity and the Mars Sample Recovery Helicopters. The analysis was performed multiple times for each rotor system with sequential increases to the level of modeling complexity. Initial cases used an independent blade analysis, trim inflow distribution, no airframe degrees of freedom, and a hover condition (a time-invariant problem). The final model used all blades, dynamic inflow, airframe degrees of freedom, and a Floquet theory approach to the flutter analysis. The approach of sequentially increasing modeling complexity was employed to ensure the causes of any observed instabilities could be identified. For each case, stability results were determined across a range of air densities, tip Mach numbers, and rotor collective pitch settings. For the Sample Recovery Helicopter, a forward flight condition analysis with the highest level of modeling complexity was also performed. The aeroelastic stability predictions showed both the Ingenuity and Sample Recovery Helicopter rotor systems to be stable for all analyzed conditions. Parameter excursions of structural properties were performed to explore conditions resulting in aeroelastic instabilities.

Aeroelastic↗
NASA NTRS1997

Three-Dimensional Viscous Alternating Direction Implicit Algorithm and Strategies for Shape Optimization

A gradient-based shape optimization based on quasi-analytical sensitivities has been extended for practical three-dimensional aerodynamic applications. The flow analysis has been rendered by a fully implicit, finite-volume formulation of the Euler and Thin-Layer Navier-Stokes (TLNS) equations. Initially, the viscous laminar flow analysis for a wing has been compared with an independent computational fluid dynamics (CFD) code which has been extensively validated. The new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4 with coarse- and fine-grid based computations performed with Euler and TLNS equations. The influence of the initial constraints on the geometry and aerodynamics of the optimized shape has been explored. Various final shapes generated for an identical initial problem formulation but with different optimization path options (coarse or fine grid, Euler or TLNS), have been aerodynamically evaluated via a common fine-grid TLNS-based analysis. The initial constraint conditions show significant bearing on the optimization results. Also, the results demonstrate that to produce an aerodynamically efficient design, it is imperative to include the viscous physics in the optimization procedure with the proper resolution. Based upon the present results, to better utilize the scarce computational resources, it is recommended that, a number of viscous coarse grid cases using either a preconditioned bi-conjugate gradient (PbCG) or an alternating-direction-implicit (ADI) method, should initially be employed to improve the optimization problem definition, the design space and initial shape. Optimized shapes should subsequently be analyzed using a high fidelity (viscous with fine-grid resolution) flow analysis to evaluate their true performance potential. Finally, a viscous fine-grid-based shape optimization should be conducted, using an ADI method, to accurately obtain the final optimized shape.

Pandya, Mohagna J.↗
NASA NTRS1980

Space Shuttle stability and control flight test techniques

A unique approach for obtaining vehicle aerodynamic characteristics during entry has been developed for the Space Shuttle. This is due to the high cost of Shuttle testing, the need to open constraints for operational flights, and the fact that all flight regimes are flown starting with the first flight. Because of uncertainties associated with predicted aerodynamic coefficients, nine flight conditions have been identified at which control problems could occur. A detailed test plan has been developed for testing at these conditions and is presented. Due to limited testing, precise computer initiated maneuvers are implemented. These maneuvers are designed to optimize the vehicle motion for determining aerodynamic coefficients. Special sensors and atmospheric measurements are required to provide stability and control flight data during an entire entry. The techniques employed in data reduction are proven programs developed and used at NASA/DFRC.

Cooke, D. R.↗
NASA NTRS1980

Computer code simulations of the formation of Meteor Crater, Arizona - Calculations MC-1 and MC-2

It has been widely accepted that hypervelocity impact processes play a major role in the evolution of the terrestrial planets and satellites. In connection with the development of quantitative methods for the description of impact cratering, it was found that the results provided by two-dimensional finite difference, computer codes is greatly improved when initial impact conditions can be defined and when the numerical results can be tested against field and laboratory data. In order to address this problem, a numerical code study of the formation of Meteor (Barringer) Crater, Arizona, has been undertaken. A description is presented of the major results from the first two code calculations, MC-1 and MC-2, that have been completed for Meteor Crater. Both calculations used an iron meteorite with a kinetic energy of 3.8 Megatons. Calculation MC-1 had an impact velocity of 25 km/sec and MC-2 had an impact velocity of 15 km/sec.

Roddy, D. J.↗
NASA NTRS1982

Utilization of normal mode initial conditions for detecting errors in the dynamics part of primitive equation global models

When a global atmospheric basic state has constant angular velocity and its temperature varies with altitude only, there exist normal mode solutions to the linearized global primitive equations. The use of these normal modes, which have known behavior in time, is superior to the use of the Rossby-Haurwitz wave as initial conditions for detecting errors in the dynamics part of primitive equation global models. With these initial conditions, integration through only one time step is sufficient to detect many formulation and coding errors. Other tests are still required for detecting problems of nonlinear instability and conservation of integral properties, however.

Chao, W. C.↗
NASA NTRS1991

A multi-point inverse airfoil design method based on conformal mapping

An exact method of multipoint inverse airfoil design for incompressible flow is presented. Multipoint design is handled by dividing the airfoil into a number of desired segments. For each segment the velocity distribution is prescribed together with an angle of attack at which the prescribed velocity distribution is to be achieved. In this manner, multipoint design objectives can be taken into account in the initial specification of the velocity distribution. In order for the multipoint inverse airfoil design problem to be well posed, three integral constraints and several conditions arise which must be satisfied. Further restrictions are imposed if the airfoil is to have a specified pitching moment, thickness ratio, or other constraints. The system of equations is solved partly as a linear system and partly through multidimensional Newton iteration. Since the velocity distribution is prescribed about the circle in the angular coordinate, specification of the velocity in terms of arc length is handled through the multidimensional Newton iteration as well. The current formulation sets the stage for a more general multipoint inverse airfoil design method in which it will be possible to specify the velocity distribution, some boundary-layer development, or the surface geometry along a segment.

Selig, Michael S.↗
NASA NTRS1993

Instabilities of two-dimensional inviscid compressible vortices

An investigation of the linear stability and subsequent nonlinear evolution and acoustic radiation of a planar inviscid compressible vortex is presented. The effects of the entropy gradient are investigated, and for the particular entropy profile chosen, the positive average entropy in the vortex core is destabilizing, while the opposite is true for the negative average entropy gradient. Finite-difference methods are used to study the linear initial value problem. These methods are found to be capable of accurately computing the frequencies and weak growth rates of the normal modes. When the initial condition consists of random perturbations, the long-time behavior is found to correspond to the most unstable normal mode in all cases. The numerical solution of the Euler equations is used to study the nonlinear development of an elliptic-mode perturbation.

Chan, W. M.↗
NASA NTRS1997

Hypersonic Vehicle Trajectory Optimization and Control

Two classes of neural networks have been developed for the study of hypersonic vehicle trajectory optimization and control. The first one is called an 'adaptive critic'. The uniqueness and main features of this approach are that: (1) they need no external training; (2) they allow variability of initial conditions; and (3) they can serve as feedback control. This is used to solve a 'free final time' two-point boundary value problem that maximizes the mass at the rocket burn-out while satisfying the pre-specified burn-out conditions in velocity, flightpath angle, and altitude. The second neural network is a recurrent network. An interesting feature of this network formulation is that when its inputs are the coefficients of the dynamics and control matrices, the network outputs are the Kalman sequences (with a quadratic cost function); the same network is also used for identifying the coefficients of the dynamics and control matrices. Consequently, we can use it to control a system whose parameters are uncertain. Numerical results are presented which illustrate the potential of these methods.

Balakrishnan, S. N.↗
NASA NTRS2006

Using Correlation to Compute Better Probability Estimates in Plan Graphs

Plan graphs are commonly used in planning to help compute heuristic "distance" estimates between states and goals. A few authors have also attempted to use plan graphs in probabilistic planning to compute estimates of the probability that propositions can be achieved and actions can be performed. This is done by propagating probability information forward through the plan graph from the initial conditions through each possible action to the action effects, and hence to the propositions at the next layer of the plan graph. The problem with these calculations is that they make very strong independence assumptions - in particular, they usually assume that the preconditions for each action are independent of each other. This can lead to gross overestimates in probability when the plans for those preconditions interfere with each other. It can also lead to gross underestimates of probability when there is synergy between the plans for two or more preconditions. In this paper we introduce a notion of the binary correlation between two propositions and actions within a plan graph, show how to propagate this information within a plan graph, and show how this improves probability estimates for planning. This notion of correlation can be thought of as a continuous generalization of the notion of mutual exclusion (mutex) often used in plan graphs. At one extreme (correlation=0) two propositions or actions are completely mutex. With correlation = 1, two propositions or actions are independent, and with correlation > 1, two propositions or actions are synergistic. Intermediate values can and do occur indicating different degrees to which propositions and action interfere or are synergistic. We compare this approach with another recent approach by Bryce that computes probability estimates using Monte Carlo simulation of possible worlds in plan graphs.

Bryce, Daniel↗
NASA NTRS1989

Study of the atmospheric general circulation as an initial value problem with the NCAR Community Climate Model

The essential elements of the transition of a model atmosphere to the equilibrium state from two extreme rest states, a cold (240 K) state and a warm (310 K) state, were examined using the NCAR Community Climate Model. It is shown that the same equilibrium state was reached in the two spinup conditions, which was also the same that was reported for this model by Williamson and Williamson (1984), who started from still a different initial condition. The spinup from the cold initial condition was faster than that from the warm condition, although the final stages of the approach to equilibrium were very simular in both cases. It was also found that the spinup of the barotropic flow lagged behind that of the baroclinic flow and the thermal field, and the model tropics spinup was faster than that of the mid-latitudes.

Chen, Tsing-Chang↗
NASA NTRS1972

Performance Analysis and Design Synthesis (PADS) computer program. Volume 2: Program description, part 2

The QL module of the Performance Analysis and Design Synthesis (PADS) computer program is described. Execution of this module is initiated when and if subroutine PADSI calls subroutine GROPE. Subroutine GROPE controls the high level logical flow of the QL module. The purpose of the module is to determine a trajectory that satisfies the necessary variational conditions for optimal performance. The module achieves this by solving a nonlinear multi-point boundary value problem. The numerical method employed is described. It is an iterative technique that converges quadratically when it does converge. The three basic steps of the module are: (1) initialization, (2) iteration, and (3) culmination. For Volume 1 see N73-13199.

Source record↗