Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 343 records · Page 19

Statistical Mechanics and Dynamics of the Outer Solar System.I. The Jupiter/Saturn Zone

We report on numerical simulations designed to understand how the solar system evolved through a winnowing of planetesimals accreeted from the early solar nebula. This sorting process is driven by the energy and angular momentum and continues to the present day. We reconsider the existence and importance of stable niches in the Jupiter/Saturn Zone using greatly improved numerical techniques based on high-order optimized multi-step integration schemes coupled to roundoff error minimizing methods.

numerical simulations planetesimals planet formati↗

Land-Atmosphere Interactions: Successes, Problems and Prospects

After two decades of active research, a much better understanding of the broader role of biospheric processes on the local climate has emerged. A surface-albedo increase, particularly in desert border regions of the subtropics (as well as the deforested tropical regions), leads to a net surface energy deficit, which in turn leads to a relative sinking and reduced rainfall. On the other hand, studies of the influence of altered ratios of evapotranspiration and sensible fluxes, in situations where the net solar income is unchanged, show that evapotranspiration is a more desirable flux for increased precipitation and vitality of the biosphere. Besides providing water vapor and convective available potential energy (CAPE) to the lower troposphere, evapotranspiration helps in building larger CAPE before "turning on" the moist-convection. Larger CAPE in the lower troposphere enables convection to reach into the deeper atmosphere thereby heating the upper troposphere; indeed, moist-convection is also accompanied by the evaporation of falling precipitation that cools and moistens the lower atmosphere. While convective, as opposed to stratiform, precipitation reduces the fractional cloud cover; it also allows more solar radiation to reach the surface thereby invigorating surface fluxes. These, together with moist convection and associated downdrafts help to maintain the characteristic upper temperature limit(s) of the moist-land as well as oceanic regions. Regardless of the above understanding, several important problems continue to hinder the accurate simulation of a realistic land atmosphere interaction in a numerical model (both GCM and/or Meso-scale models). Among the unsolved problems are parameterization of sub-grid scale land processes that include small-scale variability of soil moisture, snow-cover and snow-physics, the biodiversity of the biosphere, orography, local drainage characteristics under natural conditions, and surface flow over the natural terrain. A well-known non-linear response of surface fluxes to these variations makes the problem of parameterizing land-atmosphere interaction processes hard-to-address and simulate, particularly in a GCM. In our presentation, we will discuss how orographic, snow-cover, and water table interactions can be included into a Simple Biosphere Model such as SiB/SSiB. Figure I shows how, in the Russian region, spring snowmelt affects the soil moisture profile. Corresponding figure 2 shows how interaction with the water table decreases the natural evapotranspiration in the Sahel region simulation. While these simulations need better validation with data, the simulations reveal that surface processes are sensitive to these parameterizations. With these developments, we continue to advance our understanding of the interaction of land with the atmosphere aloft, but the intrinsic variability of the newer parameters, e. g., hydraulic properties of the soil, diminish the positive influences of these advances on the improved climate simulation with GCMs.

Sud, Y. C.↗

Thermal Conductivity Estimation from Transient Test Data with Embedded Thermocouples using Genetic Algorithm Optimization

Inverse heat transfer methodology previously developed to estimate thermal properties of high temperature fibrous insulation from embedded thermocouples was applied to the thermoplastic polymer, polyether-ether-ketone (PEEK). A small experimental setup was utilized to test the PEEK sample between temperatures of 300 K and 525 K at atmospheric pressure. Cylindrical plugs of PEEK with three thermocouples embedded at various depths were incorporated in the test sample. The experimental PEEK thermocouple data were used as the boundary and initial conditions of a one-dimensional numerical thermal model to predict the internal temperatures of the material. The thermal conductivity of PEEK was estimated with the Continuous Genetic Algorithm optimization technique by searching for the coefficients of a functional form of thermal conductivity that minimized the difference between the experimentally measured and model predicted internal temperature values. The thermal testing, one-dimensional numerical thermal model, optimization algorithm, analysis, and results are presented.

Thermal Properties↗

Study of thermal effects on nickel-cadmium batteries

Isothermal continuous flow calorimeter is designed to test a nickel-cadmium battery under numerous orbital conditions. This sensitive calorimeter collects cell data such as oxygen pressure and rate of heat generation, and calculates changes in enthalpy.

Foley, R. T.↗

Effects of inclination and eccentricity on optimal trajectories between earth and Venus

The true optimal transfers, including the effects of the inclination and eccentricity of the planets' orbits, between earth and Venus are presented as functions of the corresponding idealized Hohmann transfers. The method of determining the optimal transfers using the calculus of variations is presented. For every possible Hohmann window, specified as a continuous function of the longitude of perihelion of the Hohmann trajectory, the corresponding numerically exact optimal two-impulse transfers are given in graphical form. The cases for which the optimal two-impulse transfer is the absolute optimal, and those for which a three-impulse transfer provides the absolute optimal transfer are indicated. This information furnishes everything necessary for quick and accurate orbit calculations for preliminary Venus mission analysis. This makes it possible to use the actual optimal transfers for advanced planning in place of the standard Hohmann transfers.

Gravier, J.-P.↗

Mars - Paleostratigraphic restoration of buried surfaces in Tharsis montes

Volcanism in the Tharsis province of Mars occurred in several different areas and was generally continuous without large time intervals between eruptive episodes. Major lava flow units are numerous and extensive, but relatively thin. In many places, impact craters on buried surfaces project above younger flows that overlie them. A new application of crater dating methods has been developed to aid in the identification of these buried surfaces and to determine their lateral extent. The technique is especially adaptable to the Tharsis region where the stratigraphic succession of major flow units has been established by detailed geologic mapping. Knowledge of the overall stratigraphy allows correlations to be made between known and unknown surfaces by comparing their crater frequencies at diameters large enough to insure their recognition on the buried unit. The method has been applied to aid in the restoration of buried rock units and to construct a series of paleostratigraphic maps showing the sequence of major eruptive events in the Tharsis region.

Scott, D. H.↗

A numerical solution for two-dimensional Fredholm integral equations of the second kind with kernels of the logarithmic potential form

Two dimensional Fredholm integral equations with logarithmic potential kernels are numerically solved. The explicit consequence of these solutions to their true solutions is demonstrated. The results are based on a previous work in which numerical solutions were obtained for Fredholm integral equations of the second kind with continuous kernels.

Gabrielsen, R. E.↗

Elasticity Solution of an Adhesively Bonded Cover Plate of Various Geometries

The plane strain of adhesively bonded structures consisting of two different isotropic adherends is considered. By expressing the x-y components of the displacements in terms of Fourier integrals and using the corresponding boundary and continuity conditions, the integral equations for the general problem are obtained and solved numerically by applying Gauss-Chebyshev integration scheme. The shear and the normal stresses in the adhesive are calculated for various geometries and material properties for a stiffened plate under uniaxial tension. Numerical results involving the stress intensity factors and the strain energy release rate are presented. The closed-form expressions for the Fredholm kernels are provided to obtain the solution for an arbitrary geometry and material properties. For the general geometry, the contribution of the normal stress is quite significant, while for symmetric geometries, the shear stress is dominant, the normal stress vanishes if the adherends are of the same material and the same thickness.

Aksel, G. N.↗

Analysis of three-dimensional solidification interface shape

The determination of the shape of the solidification interface in a cooling melt by the thermal boundary conditions is investigated analytically, extending the analysis of Siegel (1982) and Siegel and Sosoka (1982) to the three-dimensional case. The solutions are generated by continuation into the solid region from conditions at the cooled wall, and some numerical results are presented in tables and graphs. These findings and the analytical approach employed are applicable to crystallization and metal-forming processes.

Siegel, R.↗

High spatial resolution passive microwave sounding systems

Two extensive series of flights aboard the ER-2 aircraft were conducted with the MIT 118 GHz imaging spectrometer together with a 53.6 GHz nadir channel and a TV camera record of the mission. Other microwave sensors, including a 183 GHz imaging spectrometer were flown simultaneously by other research groups. Work also continued on evaluating the impact of high-resolution passive microwave soundings upon numerical weather prediction models.

Staelin, D. H.↗

Iridium abundance measurements across bio-event horizons in the geological record

Geochemical studies have been performed on thousands of rock samples collected across bio-event horizons in the fossil record using INAA for about 40 common and trace elements and radiochemical isolation procedures for Os, Ir, Pt, and Au on selected samples. These studies were begun soon after the Alvarez team announced their discovery of the Cretaceous-Tertiary (K-T) Ir anomaly in marine rock sequences in Europe. With their encouragement the Authors searched for the anomaly in nearby continental (freshwater coal swamp) deposits. In collaboration with scientists from the U.S.G.S. in Denver, the anomaly was located and it was observed that a floral crisis occurred at the same stratigraphic position as the Ir spike. Further work in the Raton Basin has turned up numerous well-preserved K-T boundary sections. Although the Authors have continued to study the K-T boundary and provide geochemical measurements for other groups trying to precisely locate it, the primary effort was turned to examining the other bio-events in the Phanerozoic, especially to those that are older than the terminal Cretaceous. A list of horizons that were examined in collaboration with paleontologists and geologists is given. Results are also given and discussed.

Orth, C. J.↗

Neural computing for numeric-to-symbolic conversion in control systems

A type of neural network, the multilayer perceptron, is used to classify numeric data and assign appropriate symbols to various classes. This numeric-to-symbolic conversion results in a type of information extraction, which is similar to what is called data reduction in pattern recognition. The use of the neural network as a numeric-to-symbolic converter is introduced, its application in autonomous control is discussed, and several applications are studied. The perceptron is used as a numeric-to-symbolic converter for a discrete-event system controller supervising a continuous variable dynamic system. It is also shown how the perceptron can implement fault trees, which provide useful information (alarms) in a biological system and information for failure diagnosis and control purposes in an aircraft example.

Passino, Kevin M.↗

Manipulation strategies for massive space payloads

Control for the bracing strategy is being examined. It was concluded earlier that trajectory planning must be improved to best achieve the bracing motion. Very interesting results were achieved which enable the inverse dynamics of flexible arms to be calculated for linearized motion in a more efficient manner than previously published. The desired motion of the end point beginning at t=0 and ending at t=t sub f is used to calculate the required torque at the joint. The solution is separated into a causal function that is zero for t is less than 0 and an accusal function which is zero for t is greater than t sub f. A number of alternative end point trajectories were explored in terms of the peak torque required, the amount of anticipatory action, and other issues. The single link case is the immediate subject and an experimental verification of that case is being performed. Modeling with experimental verification of closed chain dynamics continues. Modeling effort has pointed out inaccuracies that result from the choice of numerical techniques used to incorporate the closed chain constraints when modeling our experimental prototype RALF (Robotic Arm Large and Flexible). Results were compared to TREETOPS, a multi body code. The experimental verification work is suggesting new ways to make comparisons with systems having structural linearity and joint and geometric nonlinearity. The generation of inertial forces was studied with a small arm that will damp the large arm's vibration.

Book, Wayne J.↗

Iterative LQG Controller Design Through Closed-Loop Identification

This paper presents an iterative Linear Quadratic Gaussian (LQG) controller design approach for a linear stochastic system with an uncertain open-loop model and unknown noise statistics. This approach consists of closed-loop identification and controller redesign cycles. In each cycle, the closed-loop identification method is used to identify an open-loop model and a steady-state Kalman filter gain from closed-loop input/output test data obtained by using a feedback LQG controller designed from the previous cycle. Then the identified open-loop model is used to redesign the state feedback. The state feedback and the identified Kalman filter gain are used to form an updated LQC controller for the next cycle. This iterative process continues until the updated controller converges. The proposed controller design is demonstrated by numerical simulations and experiments on a highly unstable large-gap magnetic suspension system.

Hsiao, Min-Hung↗

Transient Analysis Used to Study Thermal Radiation Effects in Single and Composite Semitransparent Layers

In a continuing in-house program at the NASA Lewis Research Center, analytical and numerical methods are being developed to apply radiative analysis to predict transient temperature distributions and heat flows in partially transmitting materials. Results have been obtained for a single plane layer, and a transient analysis is being developed for a two-layer composite where each layer has a different refractive index. Because the ceramic refractive indices are larger than one, internal reflections are produced at the surfaces and at the internal interface. Reflections tend to distribute energy within a layer, and this affects the transient temperature distributions.

Source record↗

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

Aero-Heating of Shallow Cavities in Hypersonic Freestream Flow

The purpose of these experiments and analysis was to augment the heating database and tools used for assessment of impact-induced shallow-cavity damage to the thermal protection system of the Space Shuttle Orbiter. The effect of length and depth on the local heating disturbance of rectangular cavities tested at hypersonic freestream conditions has been globally assessed using the two-color phosphor thermography method. These rapid-response experiments were conducted in the Langley 31-Inch Mach 10 Tunnel and were initiated immediately prior to the launch of STS-114, the initial flight in the Space Shuttle Return-To-Flight Program, and continued during the first week of the mission. Previously-designed and numerically-characterized blunted-nose baseline flat plates were used as the test surfaces. Three-dimensional computational predictions of the entire model geometry were used as a check on the design process and the two-dimensional flow assumptions used for the data analysis. The experimental boundary layer state conditions were inferred using the measured heating distributions on a no-cavity test article. Two test plates were developed, each containing 4 equally-spaced spanwise-distributed cavities. The first test plate contained cavities with a constant length-to-depth ratio of 8 with design point depth-to-boundary-layer-thickness ratios of 0.1, 0.2, 0.35, and 0.5. The second test plate contained cavities with a constant design point depth-to-boundary-layer-thickness ratio of 0.35 with length-to-depth ratios of 8, 12, 16, and 20. Cavity design parameters and the test condition matrix were established using the computational predictions. Preliminary results indicate that the floor-averaged Bump Factor (local heating rate nondimensionalized by upstream reference) at the tested conditions is approximately 0.3 with a standard deviation of 0.04 for laminar-in/laminar-out conditions when the cavity length-to-boundary-layer thickness is between 2.5 and 10 and for cavities in the depth-to-boundary-layer-thickness range of 0.3 to 0.8. Over this same range of conditions and parameters, preliminary results also indicate that the maximum Bump Factor on the cavity centerline falls between 2.0 and 2.75, as long as the cavity-exit conditions remain laminar. Cavities with length-to-boundary-layer-thickness ratio less than 2.5 can not be easily classified with this approach and require further analysis.

Everhart, Joel L.↗

Analysis of Separation Corridors for Visiting Vehicles from the International Space Station

The International Space Station (ISS) is a very dynamic vehicle with many operational constraints that affect its performance, operations, and vehicle lifetime. Most constraints are designed to alleviate various safety concerns that are a result of dynamic activities between the ISS and various Visiting Vehicles (VVs). One such constraint that has been in place for Russian Vehicle (RV) operations is the limitation placed on Solar Array (SA) positioning in order to prevent collisions during separation and subsequent relative motion of VVs. An unintended consequence of the SA constraint has been the impacts to the operational flexibility of the ISS resulting from the reduced power generation capability as well as from a reduction in the operational lifetime of various SA components. The purpose of this paper is to discuss the technique and the analysis that were applied in order to relax the SA constraints for RV undockings, thereby improving both the ISS operational flexibility and extending its lifetime for many years to come. This analysis focused on the effects of the dynamic motion that occur both prior to and following RV separations. The analysis involved a parametric approach in the conservative application of various initial conditions and assumptions. These included the use of the worst case minimum and maximum vehicle configurations, worst case initial attitudes and attitude rates, and the worst case docking port separation dynamics. Separations were calculated for multiple ISS docking ports, at varied deviations from the nominal undocking attitudes and included the use of two separate attitude control schemes: continuous free-drift and a post separation attitude hold. The analysis required numerical propagation of both the separation motion and the vehicle attitudes using 3-degree-of-freedom (DOF) relative motion equations coupled with rigid body rotational dynamics to generate a large set of separation trajectories.

Zaczek, Mariusz P.↗