Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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,135 records · Page 63

A Parallel Genetic Algorithm for Automated Electronic Circuit Design

Parallelized versions of genetic algorithms (GAs) are popular primarily for three reasons: the GA is an inherently parallel algorithm, typical GA applications are very compute intensive, and powerful computing platforms, especially Beowulf-style computing clusters, are becoming more affordable and easier to implement. In addition, the low communication bandwidth required allows the use of inexpensive networking hardware such as standard office ethernet. In this paper we describe a parallel GA and its use in automated high-level circuit design. Genetic algorithms are a type of trial-and-error search technique that are guided by principles of Darwinian evolution. Just as the genetic material of two living organisms can intermix to produce offspring that are better adapted to their environment, GAs expose genetic material, frequently strings of 1s and Os, to the forces of artificial evolution: selection, mutation, recombination, etc. GAs start with a pool of randomly-generated candidate solutions which are then tested and scored with respect to their utility. Solutions are then bred by probabilistically selecting high quality parents and recombining their genetic representations to produce offspring solutions. Offspring are typically subjected to a small amount of random mutation. After a pool of offspring is produced, this process iterates until a satisfactory solution is found or an iteration limit is reached. Genetic algorithms have been applied to a wide variety of problems in many fields, including chemistry, biology, and many engineering disciplines. There are many styles of parallelism used in implementing parallel GAs. One such method is called the master-slave or processor farm approach. In this technique, slave nodes are used solely to compute fitness evaluations (the most time consuming part). The master processor collects fitness scores from the nodes and performs the genetic operators (selection, reproduction, variation, etc.). Because of dependency issues in the GA, it is possible to have idle processors. However, as long as the load at each processing node is similar, the processors are kept busy nearly all of the time. In applying GAs to circuit design, a suitable genetic representation 'is that of a circuit-construction program. We discuss one such circuit-construction programming language and show how evolution can generate useful analog circuit designs. This language has the desirable property that virtually all sets of combinations of primitives result in valid circuit graphs. Our system allows circuit size (number of devices), circuit topology, and device values to be evolved. Using a parallel genetic algorithm and circuit simulation software, we present experimental results as applied to three analog filter and two amplifier design tasks. For example, a figure shows an 85 dB amplifier design evolved by our system, and another figure shows the performance of that circuit (gain and frequency response). In all tasks, our system is able to generate circuits that achieve the target specifications.

Long, Jason D.↗

Encounter with Asteroids Gaspra and Ida: Imaging Results

A progress report is given on imaging science results from the Galileo spacecraft flybys of the S-type main-belt asteroids Gaspra and Ida. Images were obtained in six colors of both asteroids, a few at resolutions better than 100 m. The overall sizes, shapes, and spin vectors of both objects are defined, and more detailed geological studies are possible for approximately half the surface of each (the hemisphere facing the spacecraft at closest approach). Each is elongated and highly irregular in profile, with heavily cratered surfaces. A number of grooves, other linear features, and surface blocks may also be associated with a recent history of impacts. Modest variations in surface albedo and color are probably due to differences in regolith development rather than compositional differences across the surfaces of each. Both the density and the size distributions of craters differ from Gaspra to Ida, in spite of its younger anticipated dynamical age, Ida is more density cratered. Ida also has a small satellite in a roughly 100-km orbit, allowing the possibility of determining the density of Ida with sufficient precision to constrain alternative models of its bulk composition.

Morrison, David↗

Satellite and Skin Layer Effects on the Accuracy of Sea Surface Temperature Measurements from the GOES Satellites

The latest Geostationary Operational Environmental Satellites (GOES) have facilitated significant improvements in our ability to measure sea surface temperature (SST) from geostationary satellites. Nonetheless, difficulties associated with sensor calibration and oceanic near-surface temperature gradients affect the accuracy of the measurements and our ability to estimate and interpret the diurnal cycle of the bulk SST. Overall, measurements of SST from the GOES Imagers on the GOES 8-10 satellites are shown to have very small bias (less than 0.02 K) and rms differences of between 0.6 and 0.9 K relative to buoy observations. Separate consideration of individual measurement times, however, demonstrates systematic bias variations of over 0.6 K with measurement hour. These bias variations significantly affect both the amplitude and shape of estimates of the diurnal SST cycle. Modeled estimates of the temperature difference across the oceanic cool skin and diurnal thermocline show that bias variations up to 0.3 K can result from variability in the near-surface layer. Oceanic near-surface layer and known "satellite midnight" calibration effects, however, explain only a portion of the observed bias variations, suggesting other possible calibration concerns. Methods of explicitly incorporating skin layer and diurnal thermocline effects in satellite bulk SST measurements were explored in an effort to further improve the measurement accuracy. While the approaches contain more complete physics, they do not yet significantly improve the accuracy of bulk SST measurements due to remaining uncertainties in the temperature difference across the near-surface layer.

Wick, Gary A.↗

Simulator Investigation of Pilot Aids for Helicopter Terminal Area Operations with One Engine Inoperative

Two recent piloted simulation experiments have investigated advanced display concepts applied to civil transport helicopter terminal area operations. Civil Category A helicopter operations apply to multi-engine helicopters wherein a safe recovery (land or fly out) is required in the event of a single engine failure. The investigation used the NASA Ames Research Center Vertical Motion Simulator, which has a full six degrees of freedom, to simulate the flight task as closely as possible. The goal of these experiments was to use advanced cockpit displays to improve flight safety and enhance the mission performance of Category A terminal area operations in confined areas. The first experiment investigated the use of military display formats to assist civil rotorcraft in performing a Category A takeoff in confined terminal areas. Specifically, it addressed how well a difficult hovering backup path could be followed using conventional instruments in comparison to panel mounted integrated displays. The hovering backup takeoff, which enables pilots to land back to the confined area pad in the event of an engine failure, was chosen since it is a difficult task to perform. Seven NASA and Army test pilots participated in the experiment. Evaluations, based on task performance and pilot workload, showed that an integrated display enabled the pilot to consistently achieve adequate or desired performance with reasonable pilot workload. Use of conventional instruments, however, frequently resulted in unacceptable performance (poor flight path tracking), higher pilot workload, and poor situational awareness. Although OEI landbacks were considered a visual task, the improved performance on the backup portion, in conjunction with increased situational awareness resulting from use of integrated displays, enabled the pilots to handle an engine failure and land back safely. In contrast, use of conventional instruments frequently led to excessive rates of sink at touchdown. A second simulation (in progress - July - August) is being conducted to investigate the use of advanced displays to perform vertical and short takeoffs and landings. One Engine Inoperative trajectories, which were optimized based on safety of flight restrictions, are utilized. Based on comments from the first experiment and further analytic development, appropriate fly out and approach guidance was added. Displays include conventional instruments with raw data, and the following integrated displays: multi-view and side-view hover displays based on the Apache Pilot Night Vision System, and variations of the pathway-in-the-sky displays with a flight-path-vector, a leader and flight director modifications. Panel mounted and head-up displays are being evaluated. Engine modifications have been incorporated to simulate 30 second and 2 minute contingency power ratings. Evaluations are based on task performance and pilot workload. NASA, Army, FAA, and industry test pilots participated. Details concerning the design, conduct, and the results of the experiment will be reported in the proposed paper.

Iseler, Laura↗

Evidence of Temporal Variation of Titan Atmospheric Density in 2005-2013

One major science objective of the Cassini mission is an investigation of Titan's atmosphere constituent abundances. Titan's atmospheric density is of interest not only to planetary scientists but also to mission design and mission control engineers. Knowledge of the dependency of Titan's atmospheric density with altitude is important because any unexpectedly high atmospheric density has the potential to tumble the spacecraft during a flyby. During low-altitude Titan flyby, thrusters are fired to counter the torque imparted on the spacecraft due to the Titan atmosphere. The denser the Titan's atmosphere is, the higher are the duty cycles of the thruster firings. Therefore thruster firing telemetry data could be used to estimate the atmospheric torque imparted on the spacecraft. Since the atmospheric torque imparted on the spacecraft is related to the Titan's atmospheric density, atmospheric densities are estimated accordingly. In 2005-2013, forty-three low-altitude Titan flybys were executed. The closest approach altitudes of these Titan flybys ranged from 878 to 1,074.8 km. Our density results are also compared with those reported by other investigation teams: Voyager-1 (in November 1980) and the Huygens Atmospheric Structure Instrument, HASI (in January 2005). From our results, we observe a temporal variation of the Titan atmospheric density in 2005-2013. The observed temporal variation is significant and it isn't due to the estimation uncertainty (5.8%, 1 sigma) of the density estimation methodology. Factors that contributed to this temporal variation have been conjectured but are largely unknown. The observed temporal variation will require synergetic analysis with measurements made by other Cassini science instruments and future years of laboratory and modeling efforts to solve. The estimated atmospheric density results are given in this paper help scientists to better understand and model the density structure of the Titan atmosphere.

Atmosphere↗

Rational approach for assumed stress finite elements

A new method for the formulation of hybrid elements by the Hellinger-Reissner principle is established by expanding the essential terms of the assumed stresses as complete polynomials in the natural coordinates of the element. The equilibrium conditions are imposed in a variational sense through the internal displacements which are also expanded in the natural co-ordinates. The resulting element possesses all the ideal qualities, i.e. it is invariant, it is less sensitive to geometric distortion, it contains a minimum number of stress parameters and it provides accurate stress calculations. For the formulation of a 4-node plane stress element, a small perturbation method is used to determine the equilibrium constraint equations. The element has been proved to be always rank sufficient.

Pian, T. H. H.↗

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather↗

The inverse gravimetric problem in gravity modelling

One of the main purposes of geodesy is to determine the gravity field of the Earth in the space outside its physical surface. This purpose can be pursued without any particular knowledge of the internal density even if the exact shape of the physical surface of the Earth is not known, though this seems to entangle the two domains, as it was in the old Stoke's theory before the appearance of Molodensky's approach. Nevertheless, even when large, dense and homogeneous data sets are available, it was always recognized that subtracting from the gravity field the effect of the outer layer of the masses (topographic effect) yields a much smoother field. This is obviously more important when a sparse data set is bad so that any smoothing of the gravity field helps in interpolating between the data without raising the modeling error, this approach is generally followed because it has become very cheap in terms of computing time since the appearance of spectral techniques. The mathematical description of the Inverse Gravimetric Problem (IGP) is dominated mainly by two principles, which in loose terms can be formulated as follows: the knowledge of the external gravity field determines mainly the lateral variations of the density; and the deeper the density anomaly giving rise to a gravity anomaly, the more improperly posed is the problem of recovering the former from the latter. The statistical relation between rho and n (and its inverse) is also investigated in its general form, proving that degree cross-covariances have to be introduced to describe the behavior of rho. The problem of the simultaneous estimate of a spherical anomalous potential and of the external, topographic masses is addressed criticizing the choice of the mixed collection approach.

Sanso, F.↗

The elemental abundance ratios of interstellar secondary and primary cosmic rays

We report new observations of abundances in the charge range (Z) between 2 and 10, which were obtained with a dE/dx-Cerenkov detector launched into a polar orbit on OGO-6 as part of the Caltech Solar and Galactic Cosmic Ray Experiment. Integral rigidity spectra of all the elements observed have shapes similar to that of the helium spectrum in the rigidity range of 2 to 14 GV, approaching a power law with exponent -1.6 above 8 GV. Calculations of interstellar propagation assuming a steady-state model and including the presence of interstellar helium and the effects of solar modulation predict a variation with rigidity of ratios such as Be-O and B/O, which is not observed. The data can be explained by assuming a rigidity-dependent confinement of cosmic rays within the Galaxy.

Brown, J. W.↗

Correlations of supersonic boundary-layer transition on cones including effects of large axial variations in wind-tunnel noise

Transition data on sharp tip cones in two pilot low disturbance wind tunnels at Mach numbers of 3.5 and 5 were correlated in terms of noise parameters with data from several conventional wind tunnels and with data from supersonic flight tests on a transition cone. The noise parameters were developed to account for the large axial variations of the free stream noise and the very high frequency noise spectra that occurred in the low disturbance tunnels for some test conditions. The noise could be varied in these tunnels from high levels, approaching those in conventional tunnels, to extremely low levels. The correlations indicated that transition in the low disturbance tunnels was dominated by the local stream noise that was incident on the cone boundary layer unstream of the neutral stability point. The correlation results also suggested that high frequency components of the low disturbance tunnel noise spectra had significant effects on transition when the noise was incident on the boundary layer both upstream and downstream of the neutral stability point.

Chen, F. J.↗

Dynamics of meteor streams

The overwhelming majority of meteor streams are generally assumed to be formed due to the decay of comets. The most effective process of the release of solid particles from a cometary nucleus is their ejection by sublimating gases when the comet approaches the Sun. The results of investigation of the Geminids and Quadrantids meteor stream evolution show that under the influence of planetary perturbations, the stream may originally be flat but then thicken depending on the variation range of orbital inclinations. Eventually, due to planetary perturbations, a meteor stream may take such a shape as to cause the start of several active showers at different solar longitudes.

Babadjanov, P. B.↗

Solar wind and coronal structure near sunspot minimum - Pioneer and SMM observations from 1985-1987

Changes in solar wind speed and magnetic polarity observed at the Pioneer spacecraft are discussed here in terms of the changing magnetic geometry implied by SMM coronagraph observations over the period 1985-1987. The pattern of recurrent solar wind streams, the long-term average speed, and the sector polarity of the interplanetary magnetic field all changed in a manner suggesting both a temporal variation, and a changing dependence on heliographic latitude. Coronal observations during this epoch show a systematic variation in coronal structure and the magnetic structure imposed on the expanding solar wind. These observations suggest interpretation of the solar wind speed variations in terms of the familiar model where the speed increases with distance from a nearly flat interplanetary current sheet, and where this current sheet becomes aligned with the solar equatorial plane as sunspot minimum approaches, but deviates rapidly from that orientation after minimum.

Mihalov, J. D.↗

High-temperature thermal storage systems for advanced solar receivers materials selections

Advanced space power systems that use solar energy and Brayton or Stirling heat engines require thermal energy storage (TES) systems to operate continuously through periods of shade. The receiver storage units, key elements in both Brayton and Stirling systems, are designed to use the latent heat of fusion of phase-change materials (PCMs). The power systems under current consideration for near-future National Aeronautics and Space Administration space missions require working fluid temperatures in the 1100 to 1400 K range. The PCMs under current investigation that gave liquid temperatures within this range are the fluoride family of salts. However, these salts have low thermal conductivity, which causes large temperature gradients in the storage systems. Improvements can be obtained, however, with the use of thermal conductivity enhancements or metallic PCMs. In fact, if suitable containment materials can be found, the use of metallic PCMs would virtually eliminate the orbit associated temperature variations in TES systems. The high thermal conductivity and generally low volume change on melting of germanium and alloys based on silicon make them attractive for storage of thermal energy in space power systems. An approach to solving the containment problem, involving both chemical and physical compatibility, preparation of NiSi/NiSi2, and initial results for containment of germanium and NiSi/NiSi2, are presented.

Wilson, D. F.↗

State-space self-tuner for on-line adaptive control

Dynamic systems, such as flight vehicles, satellites and space stations, operating in real environments, constantly face parameter and/or structural variations owing to nonlinear behavior of actuators, failure of sensors, changes in operating conditions, disturbances acting on the system, etc. In the past three decades, adaptive control has been shown to be effective in dealing with dynamic systems in the presence of parameter uncertainties, structural perturbations, random disturbances and environmental variations. Among the existing adaptive control methodologies, the state-space self-tuning control methods, initially proposed by us, are shown to be effective in designing advanced adaptive controllers for multivariable systems. In our approaches, we have embedded the standard Kalman state-estimation algorithm into an online parameter estimation algorithm. Thus, the advanced state-feedback controllers can be easily established for digital adaptive control of continuous-time stochastic multivariable systems. A state-space self-tuner for a general multivariable stochastic system has been developed and successfully applied to the space station for on-line adaptive control. Also, a technique for multistage design of an optimal momentum management controller for the space station has been developed and reported in. Moreover, we have successfully developed various digital redesign techniques which can convert a continuous-time controller to an equivalent digital controller. As a result, the expensive and unreliable continuous-time controller can be implemented using low-cost and high performance microprocessors. Recently, we have developed a new hybrid state-space self tuner using a new dual-rate sampling scheme for on-line adaptive control of continuous-time uncertain systems.

Shieh, L. S.↗

Sensitivity analysis of frictional contact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Karaoglan, I.↗

Dynamic sensitivity analysis of frictional contact/impact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Karaoglan, Levent↗

Catalogue of X-Ray Texture Data for Al-Cu-Li Alloy 1460, 2090, 2096 and 2195 Near-Net-Shape Extrusions, Sheet and Plate

The effect of crystallographic texture on the mechanical properties of near-net-shape extrusions is of major interest ff these products are to find application in launch vehicle or aircraft structures. The objective of this research was to produce a catalogue containing quantitative texture information for extruded product, sheet and plate. The material characterized was extracted from wide, integrally stiffened panels fabricated from the Al-Cu-Li alloys 1460, 2090, 2096 and 2195. The textural characteristics of sheet and plate products of the same alloys were determined for comparison purposes. The approach involved using X-ray diffraction to generate pole figures in combination with orientation distribution function analysis. The data were compiled as a function of location in the extruded cross-sections and the variation in the major deformation- and recrystallization-related texture components was identified.

Hales, Stephen J.↗

A Multifunctional Interface Method for Coupling Finite Element and Finite Difference Methods: Two-Dimensional Scalar-Field Problems

A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.

Ransom, Jonathan B.↗