Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical stability”

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 1,117 records · Page 62

Simulation of Heavy Lift Airship dynamics over large ranges of incidence and speed

A nonlinear, multibody, six-degrees-of-freedom digital simulation has been developed to study generic Heavy Lift Airship (HLA) dynamics and control. The basic aerodynamic functions developed to model the hull, tail, and rotor loads continuously over all incidence ranges are reviewed and applied to a Quadrotor HLA with a low fineness ratio hull and a small vee-tail. Trim calculations for a test vehicle suggest control power deficiencies in crosswind stationkeeping for the unloaded vehicle. Gust responses show the importance of correctly calculating loads due to accelerated relative motion of air and hull. Numerically linearized dynamics for the test vehicle show the existence of a divergent yaw mode, and an oscillatory pitch mode whose stability characteristics are sensitive to flight speed. A considerable improvement in the vehicle's stability and response results from a simple multi-axis closed-loop control system operating on the rotors and propeller blades.

Tischler, M. B.↗

Some lessons from NACA/NASA aerodynamic studies following World War II

An historical account is presented of the new departures in aerodynamic research conducted by NACA, and subsequently NASA, as a result of novel aircraft technologies and operational regimes encountered in the course of the Second World War. The invention and initial development of the turbojet engine furnished the basis for a new speed/altitude regime in which numerous aerodynamic design problems arose. These included compressibility effects near the speed of sound, with attendant lift/drag efficiency reductions and longitudinal stability enhancements that were accompanied by a directional stability reduction. Major research initiatives were mounted in the investigation of swept, delta, trapezoidal and variable sweep wing configurations, sometimes conducted through flight testing of the 'X-series' aircraft. Attention is also given to the development of the first generation of supersonic fighter aircraft.

Spearman, M. L.↗

Evaluation of 2 1-D cloud models for the analysis of VAS soundings

Evaluation of the satellite Visual Infrared Spin Scan Radiometer Atmospheric Sounder (VISSR) has begun to document several of its critical shortcomings as far as numerical cloud models are concerned: excessive smoothing of thermal inversions; imprecise measurement of boundary layer moisture; and tendency to exaggerate atmospheric stability. The sensitivity of 1-D cloud models to their required inputs is stressed with special attention to those parameters obtained from atmospheric soundings taken by the VAS or rawinsonde. In addition to performing model experiments using temperature and moisture profiles having the general characteristics of VAS soundings, standard input sensitivity tests were made and 1-D model performance was compared with observations and the results of a 2-D model experiment using AVE/VAS data (Atmospheric Variability Experiment). Although very encouraging, the results are not sufficient to make any specific conclusions. In general, the VAS soundings are likely to be inadequate to provide the cloud base (and subcloud layer) information needed for inputs to current cumulus models. Above cloud base, the tendency to exaggerate the stability of the atmosphere requires solution before meaningful model experiments are run.

Emmitt, G. D.↗

Radiometric calibration analysis of SIR-B imagery

The second flight of the NASA Spaceborne Imaging Radar (SIR-B) collected nearly 8 hours of digital synthetic-aperture radar (SAR) data of the earth's surface. This paper analyzes the performance of the SAR to determine the accuracy of the relative radiometric calibration of these data. Procedures are presented for deriving accuracy estimates as well as limitations under which these estimates are valid. The impact on calibration of an RF breakdown in the primary antenna feed system cable is evaluated. An analysis of the platform stability based on the SAR echo data is also presented in conjunction with its resultant effect on the calibration accuracy. Finally, numerical error bounds are derived with guidelines on their utilization.

Wall, Stephen D.↗

Theoretical motions of hydrofoil systems

Results are presented of an investigation that has been undertaken to develop theoretical methods of treating the motions of hydrofoil systems and to determine some of the important parameters. Variations of parameters include three distributions of area between the hydrofoils, two rates of change of downwash angle with angle of attack, three depths of immersion, two dihedral angles, two rates of change of lift with immersion, three longitudinal hydrofoil spacings, two radii of gyration in pitching, and various horizontal and vertical locations of the center of gravity. Graphs are presented to show locations of the center of gravity for stable motion, values of the stability roots, and motions following the sudden application of a vertical force or a pitching moment to the hydrofoil system for numerous sets of values of the parameters.

Imlay, Frederick H↗

Bifurcations and sensitivity in parametric nonlinear programming

The parametric nonlinear programming problem is that of determining the behavior of solution(s) as a parameter or vector of parameters alpha belonging to R(sup r) varies over a region of interest for the problem: Minimize over x the set f(x, alpha):h(x, alpha) = 0, g(x, alpha) is greater than or equal to 0, where f:R(sup (n+r)) approaches R, h:R(sup (n+r)) approaches R(sup q) and g:R(sup (n+r)) approaches R(sup p) are assumed to be at least twice continuously differentiable. Some of these parameters may be fixed but not known precisely and others may be varied to enhance the performance of the system. In both cases a fundamentally important problem in the investigation of global sensitivity of the system is to determine the stability boundaries of the regions in parameter space which define regions of qualitatively similar solutions. The objective is to explain how numerical continuation and bifurcation techniques can be used to investigate the parametric nonlinear programming problem in a global sense. Thus, first the problem is converted to a closed system of parameterized nonlinear equations whose solution set contains all local minimizers of the original problem. This system, which will be represented as F(z,alpha) = O, will include all Karush-Kuhn-Tucker and Fritz John points, both feasible and infeasible solutions, and relative minima, maxima, and saddle points of the problem. The local existence and uniqueness of a solution path (z(alpha), alpha) of this system as well as the solution type persist as long as a singularity in the Jacobian D(sub z)F(z,alpha) is not encountered. Thus the nonsingularity of this Jacobian is characterized in terms of conditions on the problem itself. Then, a class of efficient predictor-corrector continuation procedures for tracing solution paths of the system F(z,alpha) = O which are tailored specifically to the parametric programming problem are described. Finally, these procedures and the obtained information are illustrated within the context of design optimization.

Lundberg, Bruce N.↗

On the Long Period Luni-Solar Effect in the Motion of an Artificial Satellite

Two systems of formulas are presented for the determination of the long period perturbations caused by the Sun and the Moon in the motion of an artificial satellite. The first system can be used to determine the lunar effect for all satellites. The second method is more convenient for finding the lunar effect for close satellites and the solar effect for all satellites. Knowledge of these effects is essential for determining the stability of the satellite orbit. The basic equations of both systems are arranged in a form which permits the use of numerical integration. The two theories are more accurate and more adaptable to the use of electronic machines than the analytical developments obtained previously.

Musen, Peter↗

Dynamics and Control of Articulated Anisotropic Timoshenko Beams

The paper illustrates the use of continuum models in control design for stabilizing flexible structures. A 6-DOF anisotropic Timoshenko beam with discrete nodes where lumped masses or actuators are located provides a sufficiently rich model to be of interest for mathematical theory as well as practical application. We develop concepts and tools to help answer engineering questions without having to resort to ad hoc heuristic ("physical") arguments or faith. In this sense the paper is more mathematically oriented than engineering papers and vice versa at the same time. For instance we make precise time-domain solutions using the theory of semigroups of operators rather than formal "inverse Laplace transforms." We show that the modes arise as eigenvalues of the generator of the semigroup, which are then related to the eigenvalues of the stiffness operator. With the feedback control, the modes are no longer orthogonal and the question naturally arises as to whether there is still a modal expansion. Here we prove that the eigenfunctions yield a biorthogonal Riesz basis and indicate the corresponding expansion. We prove mathematically that the number of eigenvalues is nonfinite, based on the theory of zeros of entire functions. We make precise the notion of asymptotic modes and indicate how to calculate them. Although limited by space, we do consider the root locus problem and show for instance that the damping at first increases as the control gain increases but starts to decrease at a critical value, and goes to zero as the gain increases without bound. The undamped oscillatory modes remain oscillatory and the rigid-body modes go over into deadbeat modes. The Timoshenko model dynamics are translated into a canonical wave equation in a Hilbert space. The solution is shown to require the use of an "energy" norm which is no more than the total energy: potential plus kinetic. We show that, under an appropriate extension of the notion of controllability, rate feedback with a collocated sensor can stabilize the structure in the sense that all modes are damped and the energy decays to zero. An example, non-numeric, is worked out in some detail illustrating the concepts and theory developed.

Balakrishnan, A. V.↗

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗

Considerations and Architectures for Inter-Satellite Communications in Distributed Spacecraft Systems

This paper will identify the important characteristics and requirements necessary for inter-satellite communications in distributed spacecraft systems and present analysis results focusing on architectural and protocol comparisons. Emerging spacecraft systems plan to deploy multiple satellites in various "distributed" configurations ranging from close proximity formation flying to widely separated constellations. Distributed spacecraft configurations provide advantages for science exploration and operations since many activities useful for missions may be better served by distributing them between spacecraft. For example, many scientific observations can be enhanced through spatially separated platforms, such as for deep space interferometry. operating multiple distributed spacecraft as a mission requires coordination that may be best provided through inter-satellite communications. For example, several future distributed spacecraft systems envision autonomous operations requiring relative navigational calculations and coordinated attitude and position corrections. To conduct these operations, data must be exchanged between spacecraft. Direct cross-links between satellites provides an efficient and practical method for transferring data and commands. Unlike existing "bent-pipe" relay networks supporting space missions, no standard or widely-used method exists for cross-link communications. Consequently, to support these future missions, the characteristics necessary for inter-satellite communications need to be examined. At first glance, all of the missions look extremely different. Some missions call for tens to hundreds of nano-satellites in constant communications in close proximity to each other. Other missions call for a handful of satellites communicating very slowly over thousands to hundreds of thousands of kilometers. The paper will first classify distributed spacecraft missions to help guide the evaluation and definition of cross-link architectures and approaches. Based on this general classification, the paper will examine general physical layer parameters, such as frequency bands and data rates, necessary to support the missions. The paper will also identify classes of communication architectures that may be employed, ranging from fully distributed to centralized topologies. Numerous factors, such as number of spacecraft, must be evaluated when attempting to pick a communications architecture. Also important is the stability of the formation from a communications standpoint. For example, do all of the spacecraft require equal bandwidth and are spacecraft allowed to enter and leave a formation? The type of science mission being attempted may also heavily influence the communications architecture. In addition, the paper will assess various parameters and characteristics typically associated with the data link layer. The paper will analyze the performance of various multiple access techniques given the operational scenario, requirements, and communication topologies envisioned for missions. This assessment will also include a survey of existing standards and their applicability for distributed spacecraft systems. An important consideration includes the interoperability of the lower layers (physical and data link) examined in this paper with the higher layer protocols(network) envisioned for future space internetworking. Finally, the paper will define a suggested path, including preliminary recommendations, for defining and developing a standard for intersatellite communications based on the classes of distributed spacecraft missions and analysis results.

Edwards, Bernard↗

On the Origin of the Saturnian Satellite System: Did Iapetus Form In-Situ?

Current models of planet and satellite formation are marred by our lack of understanding regarding the turbulent state of accretion disks. According to the Rayleigh criterion, Keplerian disks are hydrodynamically stable. Indeed, it has been argued that a carefully designed Taylor Couette experiment shows stability in the case of positive radial gradients in specific angular momentum even for high Reynolds numbers [1], in agreement with numerical simulations which consistently show turbulence decay [2]. Other possible sources of turbulence may fail due to low ionization, may decay as the optical depth decreases due to dust coagulation, may involve unrealistic boundary conditions, or result in limited transport. The difficulty stems not only from the degree of turbulence, but also from the kind of turbulence, and whether it may be characterized by an alpha parameter.

Mosqueira, I.↗

Quasi-Static Analysis of LaRC THUNDER Actuators

An analytic approach is developed to predict the shape and displacement with voltage in the quasi-static limit of LaRC Thunder Actuators. The problem is treated with classical lamination theory and Von Karman non-linear analysis. In the case of classical lamination theory exact analytic solutions are found. It is shown that classical lamination theory is insufficient to describe the physical situation for large actuators but is sufficient for very small actuators. Numerical results are presented for the non-linear analysis and compared with experimental measurements. Snap-through behavior, bifurcation, and stability are presented and discussed.

Campbell, Joel F.↗

Quasi-Static Analysis of Round LaRC THUNDER Actuators

An analytic approach is developed to predict the shape and displacement with voltage in the quasi-static limit of round LaRC Thunder Actuators. The problem is treated with classical lamination theory and Von Karman non-linear analysis. In the case of classical lamination theory exact analytic solutions are found. It is shown that classical lamination theory is insufficient to describe the physical situation for large actuators but is sufficient for very small actuators. Numerical results are presented for the non-linear analysis and compared with experimental measurements. Snap-through behavior, bifurcation, and stability are presented and discussed.

Campbell, Joel F.↗

Rotation of a Moonless Earth

We numerically explore the obliquity (axial tilt) variations of a hypothetical moonless Earth. Previous work has shown that the Earth's Moon stabilizes Earth's obliquity such that it remains within a narrow range, between 22.1 deg and 24.5 deg. Without lunar influence, a frequency-map analysis by Laskar et al. showed that the obliquity could vary between 0 deg. and 85 deg. This has left an impression in the astrobiology community that a large moon is necessary to maintain a habitable climate on an Earth-like planet. Using a modified version of the orbital integrator mercury, we calculate the obliquity evolution for moonless Earths with various initial conditions for up to 4 Gyr. We find that while obliquity varies significantly more than that of the actual Earth over 100,000 year timescales, the obliquity remains within a constrained range, typically 20-25 deg. in extent, for timescales of hundreds of millions of years. None of our Solar System integrations in which planetary orbits behave in a typical manner show obliquity accessing more than 65% of the full range allowed by frequency-map analysis. The obliquities of moonless Earths that rotate in the retrograde direction are more stable than those of pro-grade rotators. The total obliquity range explored for moonless Earths with rotation periods shorter than 12 h is much less than that for slower-rotating moonless Earths. A large moon thus does not seem to be needed to stabilize the obliquity of an Earth-like planet on timescales relevant to the development of advanced life.

orbital integrator mercury↗

Midlatitude Ice-Rich Ground on Mars: An Important Target for Science and In Situ Resource Utilization on Human Missions

The region of ROI is characterized by proven presence of near surface ground ice and numerous periglacial features. Midlatitude ground ice on Mars is of significant scientific interest for understanding the history and evolution of ice stability on Mars, the impact that changes in insolation produced by variations in Mars’ orbital parameters has on the regions climate, and could provide human exploration with a reliable and plentiful in situ resource. For both science and exploration, assessing the astrobiological potential of the ice is important in terms of (1) understanding the potential for life on Mars and (2) evaluating the presence of possible biohazards in advance of human exploration. Heldmann et al. (2014) studied locations on Mars in the Amazonis Planitia region where near surface ground ice was exposed by new impact craters (Byrne et al. 2009). The study examined whether sites in this region were suitable for human exploration including reviewing the evidence for midlatitude ground ice, discussing the possible explanations for its occurrence, assessing its potential habitability for modern life, and evaluating the resource potential. They systematically analyzed remote-sensing data sets to identify a viable landing site. Five sites where ground ice was exposed were examined with HiRise imaging and were classified according to (1) presence of polygons as a proxy for subsurface ice, (2) presence and abundance of rough topographic obstacles (e.g., large cracks, cliffs, uneven topography), (3) rock density, (4) presence and abundance of large boulders, and (5) presence of craters. A suitable landing site was found having ground ice at only 0.15m depth, and no landing site hazards within a 25 km landing ellipse. This paper presents results of that study and examines the relevance of this ROI to the workshop goals.

Human Missions↗

Experimental Characterization of the Dimensional Stability of Deployable Composite Booms During Stowage

Thin-ply carbon fiber reinforced polymer materials offer the opportunity to create long thin-shell booms with high strength to weight ratios. These booms can be rolled up and stowed in small volumes to be later deployed. These features make them ideal for use in gossamer solar sails as well as other deployable space structures. It has been observed that stowing the booms in rolled-up configurations can cause time-dependent deformations. A comprehensive understanding of the detrimental effects of stowing the coiled booms is necessary to implement these structures in spaceflight missions. During stowage inside the spacecraft, the boom is subjected to a constant deformation/strain that causes it to relax over time, measurable by a decrease in stress. Using large deformation bending tests, the complete fold-stow-unfold-recover cycle of thin-laminate candidates for the rollable booms is characterized. Further bending and stowage testing is done at the boom level. The results from this test campaign will be used to evaluate the dimensional stability of deployable booms made from thin-ply composites that exhibit visco-elasto-plastic behavior and calibrate and validate numerical finite element models that include visco-elasto-plastic material models.

Experimental Characterization↗

Linear instability of the supersonic wake behind a flat plate aligned with a uniform stream

A theoretical and numerical study is presented of the growth of linear disturbances in the laminar compressible supersonic wake behind a flat plate aligned with a uniform stream. The approach adopted is that of classical linear stability theory and covers both two- and three-dimensional disturbances. The basic flow is described, and the linear stability problem is formulated that leads to the compressible Rayleigh equation. The stability of short-wave disturbances is then analyzed in the very near wake where the compressible analog of the Goldstein double structure provides an asymptotic description for the unperturbed flow. The behavior of long waves is then considered. The waves are not as long as the plate length but are long compared to the boundary-layer thickness. Numerical solutions of the Rayleigh equation are discussed. Numerical results are then presented for a range of streamwise stations, Mach numbers, and wave propagation angles.

Papageorgiou, D. T.↗

Modeling, model simplification and stability robustness with aeroelastic vehicles

The linearization and simplification of a nonlinear, literal model for flexible aircraft is highlighted. Areas of model fidelity that are critical if the model is to be used for control system synthesis are developed, and several simplification techniques that can deliver the necessary model fidelity are discussed. These techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle mass properties, vibrational characteristics, and rigid-body and aeroelastic stability derivatives, thus leading to the underlying causes for critical dynamic characteristics.

Schmidt, David K.↗