Modified thermal analyzer digital computer program (TAP)
Computer program for solving n-dimensional transient or steady state heat flow problems by creating electrical analogy of problem and solving by finite difference method
SEARCH · Engineering Papers
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.
Computer program for solving n-dimensional transient or steady state heat flow problems by creating electrical analogy of problem and solving by finite difference method
Natural convection effects in phase change thermal control devices were studied. A mathematical model was developed to evaluate natural convection effects in a phase change test cell undergoing solidification. Although natural convection effects are minimized in flight spacecraft, all phase change devices are ground tested. The mathematical approach to the problem was to first develop a transient two-dimensional conduction heat transfer model for the solidification of a normal paraffin of finite geometry. Next, a transient two-dimensional model was developed for the solidification of the same paraffin by a combined conduction-natural-convection heat transfer model. Throughout the study, n-hexadecane (n-C16H34) was used as the phase-change material in both the theoretical and the experimental work. The models were based on the transient two-dimensional finite difference solutions of the energy, continuity, and momentum equations.
Effects of inelastic materials models on computed stress-strain solutions for thermally loaded structures were studied by performing nonlinear (elastoplastic creep) and elastic structural analyses on a prismatic, double edge wedge specimen of IN 100 alloy that was subjected to thermal cycling in fluidized beds. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic kinematic, and combined plus transient creep) were exercised for the problem by using the MARC nonlinear, finite element computer program. Maximum total strain ranges computed from the elastic and nonlinear analyses agreed within 5 percent. Mean cyclic stresses, inelastic strain ranges, and inelastic work were significantly affected by the choice of inelastic constitutive model. The computing time per cycle for the nonlinear analyses was more than five times that required for the elastic analysis.
This effort is in support of the development of the expert system of computer codes to predict the loads on select structural components of a space propulsion engine. The development will be based primarily on the space shuttle main engine (SSME) test data base. Because of random variations of the many different sources of the loadings on the selected structural components and transients, a probabilistic approach to the problems was adopted. The goal of this task is to characterize all of the individual sources of loads at critical structural locations, such as the turbine blades, the transfer ducts, and liquid oxygen posts, using state-of-the-art probabilistic methods with varying levels of sophistication. The second phase of this work is the development of a composite load model based on a probabilistic synthesis of the individual load model previously developed. This model will be based on the stochastic combination of the load variables and not on the physical process for the combination of the individual loads to the composite load seen by the selected structural components.
Differences between retrieved temperatures from lidar and atmospheric model temperatures are calculated for four problems: (1) the effect of a transient thin aerosol layer at any altitude on temperature retrieved by a single-wavelength lidar; (2) the effect of residual aerosols on temperature retrieved by a two-wavelength lidar; (3) the errors due to splicing two backscatter signals which are obtained separately; and (4) the effect of multiple scattering. The results show that temperatures can be retrieved for + or - 3 K if Rayleigh backscatter can be measured for + or - 1.5 percent. The LITE lidar in a 240 km orbit aboard the Shuttle should be able to make a + or -3 K measurements from 10-40 km with vertical resolutions of 1 km and horizontal resolutions of 300 km.
The structural analysis methods research has several goals. One goal is to develop analysis methods that are general. This goal of generality leads naturally to finite-element methods, but the research will also include other structural analysis methods. Another goal is that the methods be amenable to error analysis; that is, given a physical problem and a mathematical model of that problem, an analyst would like to know the probable error in predicting a given response quantity. The ultimate objective is to specify the error tolerances and to use automated logic to adjust the mathematical model or solution strategy to obtain that accuracy. A third goal is to develop structural analysis methods that can exploit parallel processing computers. The structural analysis methods research will focus initially on three types of problems: local/global nonlinear stress analysis, nonlinear transient dynamics, and tire modeling.
Three problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. The first problem is the calculation of resonance within the continuous spectrum of the Blasius boundary layer. The second is calculation of the growth of Tollmien-Schlichting waves that are a direct result of disturbances that only lie outside of the boundary layer. And, the third problem is the calculation of non-parallel effects. Together, these problems represent a unified approach to the study of freestream disturbance effects that could lead to transition. Solutions to the temporal, initial-value problem with an inhomogeneous forcing term imposed upon the flow is sought. By solving a series of problems, it is shown that: A transient disturbance lying completely outside of the boundary layer can lead to the growth of an unstable Tollmien-Schlichting wave. A resonance with the continuous spectrum leads to strong amplification that may provide a mechanism for bypass transition once nonlinear effects are considered. A disturbance with a very weak unstable Tollmien-Schlichting wave can lead to a much stronger Tollmien-Schlichting wave downstream, if the original disturbance has a significant portion of its energy in the continuum modes.
Existing tendon-driven fingers have applied force control through independent tension controllers on each tendon, i.e. in the tendon-space. The coupled kinematics of the tendons, however, cause such controllers to exhibit a transient coupling in their response. This problem can be resolved by alternatively framing the controllers in the joint-space of the manipulator. This work presents a joint-space torque control law that demonstrates both a decoupled and significantly faster response than an equivalent tendon-space formulation. The law also demonstrates greater speed and robustness than comparable PI controllers. In addition, a tension distribution algorithm is presented here to allocate forces from the joints to the tendons. It allocates the tensions so that they satisfy both an upper and lower bound, and it does so without requiring linear programming or open-ended iterations. The control law and tension distribution algorithm are implemented on the robotic hand of Robonaut-2.
An effort to apply the reduced basis method to nonlinear transient thermal analysis is described. The method combines the classical Rayleigh-Ritz and modal superposition techniques with contemporary finite element methods to retain modeling versatility as the degrees of freedom in a problem are reduced. The essence of the method is to use a few thermal modes from eigenvalue analyses as basis vectors to represent the temperature response for a given thermal problem similar to the use of vibration modes to represent displacements in a dynamic response problem. Approximate temperature distributions were obtained using the reduced basis method for a small section of the Shuttle Orbiter lower wing undergoing reentry heating. Good agreement was obtained between the reduced basis method solutions and full system solutions with reductions in the degrees of freedom of up to a factor of four. The good agreement indicates the reduced basis method has the potential for significant reduction in computing effort for thermal analysis; however, considerable work remains to determine techniques for selecting the type and number of basis vectors needed for approximate solutions to more complex transient thermal problems.
A solution of the linear initial value problem for the model of Eady with the inclusion of Ekman damping is presented. This model exhibits large transient growth of perturbations for synoptic cyclone spatial scales and a realistic value of the vertical turbulent viscosity coefficient despite the fact that all normal modes are exponentially decaying. Similar results are found for the model of Charney, implying that exponential instability cannot, in general, serve to explain the occurrence of cyclone scale disturbances in midlatiudes. Rather these are seen to arise additionally and perhaps predominantly from the release of mean flow potential energy by favorably configured initial perturbations. The Petterssen criterion for midlatitude cyclogenesis results naturally from this development as does its extension to the formation of subtropical monsoon depressions. Implications for the maintenance of midlatitude temperature gradients are discussed.
The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. This proposed methodology would be readily adaptable to existing computer programs for structural thermal analysis.
The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.
The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally-useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.
There has been, recently, a revival of the stability problem of accretion disks. Much of this renewed interest is due to recent observational data on transient soft X-ray novae, which are low-mass X-ray binaries. It is widely believed that nonsteady mass transfer from the secondary onto the compact primary, through an accretion disk, is the reason for the observed spectacular events in the form of often repetitive outbursts, with recurrence times ranging from 1 to 60 yr and duration time on the scale of months. Though not having reached yet a consensus about the nature of the mechanism that regulates the mass transfer, the disk thermal instability model seems to be favored by the fact that the rise in the hard X-ray luminosity is prior to the rise in the soft X-ray luminosity, while the mass transfer instability model seems to be hindered by the fact that the luminosity during quiescence is unable to trigger the thermal instability. However, it should be stressed that, remarkably, the X-ray light curves of these X-ray novae all show overall exponential decays, a feature quite difficult to reproduce in the framework of the viscous disk model, which yields powerlike luminosity decay. Taking into account this observational constraint, we have studied the temporal evolution of perturbations in the accretion rate, under the assumption that alpha is radial and parameter dependent. The chosen dependence is such that the model can reproduce limit cycle behavior (the system is locally unstable but globally stable). However, the kind of dependence we are looking for in alpha does not allow us to use the usual Shakura and Sunyaev procedure in the sense that we no longer can obtain a linearized continuity equation without explicit dependence on the accretion rate. This is so because now we cannot eliminate the accretion rate by using the angular momentum conservation equation.
State and parameter identifications based on a form of the maximum likelihood method are applied to the problem of extracting linear perturbation models, including rotor dynamic inflow effects, from transient blade flapping measurements. The estimation method is first studied in computer simulations and then applied to cyclic pitch stirring transients generated with a four-bladed rotor model operating in hovering trim conditions. The analytical perturbation models extracted from the transient test results are compared with transient and frequency response tests not used in the state and parameter identification. The identified analytical perturbation model is also compared with a simple theory. The method that is applicable both to small scale and full scale dynamic rotor testing is being extended to perturbations from forward flight trim conditions.
The development of a 4-bladed model rotor is reported that can be excited with a simple eccentric mechanism in progressing and regressing modes with either harmonic or transient inputs. Parameter identification methods were applied to the problem of extracting parameters for linear perturbation models, including rotor dynamic inflow effects, from the measured blade flapping responses to transient pitch stirring excitations. These perturbation models were then used to predict blade flapping response to other pitch stirring transient inputs, and rotor wake and blade flapping responses to harmonic inputs. The viability and utility of using parameter identification methods for extracting the perturbation models from transients are demonstrated through these combined analytical and experimental studies.
Characterization of the three-dimensional structure of solar transients using incomplete plane of sky data is a difficult problem whose solutions have potential for societal benefit in terms of space weather applications. In this paper transients are characterized in three dimensions by means of conic coronal mass ejection (CME) approximation. A novel method for the automatic determination of cone model parameters from observed halo CMEs is introduced. The method uses both standard image processing techniques to extract the CME mass from white-light coronagraph images and a novel inversion routine providing the final cone parameters. A bootstrap technique is used to provide model parameter distributions. When combined with heliospheric modeling, the cone model parameter distributions will provide direct means for ensemble predictions of transient propagation in the heliosphere. An initial validation of the automatic method is carried by comparison to manually determined cone model parameters. It is shown using 14 halo CME events that there is reasonable agreement, especially between the heliocentric locations of the cones derived with the two methods. It is argued that both the heliocentric locations and the opening half-angles of the automatically determined cones may be more realistic than those obtained from the manual analysis
Experimental investigation of the problems associated with restarting hybrid rocket motors (i.e., motors wherein a liquid or gaseous oxidizer is injected into the port of a solid fuel grain with subsequent mixing and combustion of the oxidizer and fuel) following a brief period of extinguishment. The results include the finding that the ignition delay on restart is decreased because less energy is absorbed by the fuel before the surface reaches the ignition point.