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 649 records · Page 36

Robust Stability Analysis of the Space Launch System Control Design: A Singular Value Approach

Classical stability analysis consists of breaking the feedback loops one at a time and determining separately how much gain or phase variations would destabilize the stable nominal feedback system. For typical launch vehicle control design, classical control techniques are generally employed. In addition to stability margins, frequency domain Monte Carlo methods are used to evaluate the robustness of the design. However, such techniques were developed for Single-Input-Single-Output (SISO) systems and do not take into consideration the off-diagonal terms in the transfer function matrix of Multi-Input-Multi-Output (MIMO) systems. Robust stability analysis techniques such as H(sub infinity) and mu are applicable to MIMO systems but have not been adopted as standard practices within the launch vehicle controls community. This paper took advantage of a simple singular-value-based MIMO stability margin evaluation method based on work done by Mukhopadhyay and Newsom and applied it to the SLS high-fidelity dynamics model. The method computes a simultaneous multi-loop gain and phase margin that could be related back to classical margins. The results presented in this paper suggest that for the SLS system, traditional SISO stability margins are similar to the MIMO margins. This additional level of verification provides confidence in the robustness of the control design.

Pei, Jing↗

Orbiter altitude at ALT interface based on ALSES (Approach and Landing Shuttle Engineering Simulation) analyses. Mission planning, mission analysis and software formulation

The details of the orbiter altitude attainable at the approach and landing tests (ALT) interface determined by the approach and landing shuttle engineering simulation (manned) were documented. The analysis culminated in the verification of the trends observed in a similar analysis performed previously on the space vehicle dynamics simulation (unmanned). Altitude variations attributable to pilot steering variability ranged between 492 ft higher to 383 ft lower. The requirement for this parametric analysis is first elaborated. The specifications, assumptions, and analytical approach used to determine the orbiter altitude at the ALT interface are then presented, followed by the results of the analytical approach and the conclusions and recommendations.

Colwell, C. L.↗

A two-dimensional, time-dependent, near-earth magnetotail

A kinetic approach is presented for obtaining slowly time-dependent, self-consistent plasma, and magnetic field configurations for this region, using anisotropic particle distributions. At present, a 2D configuration is used. The time variation of the plasma/field configuration is assumed to be driven by a slowly varying, externally imposed electric field (representing, for example, a change in solar wind conditions). The numerical approach involves solving a 2D nonlinear Poisson equation for the magnetic vector potential at each time step. The local current density is calculated from particle velocity distributions under the assumption of adiabatic particle motion. The obtained equations of state are compared with those of double-adiabatic theory. Numerical examples of self-consistent responses of the near-earth plasma and field to an externally imposed electric field pulse are presented.

Whipple, Elden↗

Elevated temperature crack growth

Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Kim, K. S.↗

Remote Sensing of Crystal Shapes in Ice Clouds

Ice crystals in clouds exist in a virtually limitless variation of geometries. The most basic shapes of ice crystals are columnar or plate-like hexagonal prisms with aspect ratios determined by relative humidity and temperature. However, crystals in ice clouds generally display more complex structures owing to aggregation, riming and growth histories through varying temperature and humidity regimes. Crystal shape is relevant for cloud evolution as it affects microphysical properties such as fall speeds and aggregation efficiency. Furthermore, the scattering properties of ice crystals are affected by their general shape, as well as by microscopic features such as surface roughness, impurities and internal structure. To improve the representation of ice clouds in climate models, increased understanding of the global variation of crystal shape and how it relates to, e.g., location, cloud temperature and atmospheric state is crucial. Here, the remote sensing of ice crystal macroscale and microscale structure from airborne and space-based lidar depolarization observations and multi-directional measurements of total and polarized reflectances is reviewed. In addition, a brief overview is given of in situ and laboratory observations of ice crystal shape as well as the optical properties of ice crystals that serve as foundations for the remote sensing approaches. Lidar depolarization is generally found to increase with increasing cloud height and to vary with latitude. Although this variation is generally linked to the variation of ice crystal shape, the interpretation of the depolarization remains largely qualitative and more research is needed before quantitative conclusions about ice shape can be deduced. The angular variation of total and polarized reflectances of ice clouds has been analyzed by numerous studies in order to infer information about ice crystal shapes from them. From these studies it is apparent that pristine crystals with smooth surfaces are generally inconsistent with the data and thus crystal impurity, distortion or surface roughness is prevalent. However, conclusions about the dominating ice shapes are often inconclusive and contradictory and are highly dependent on the limited selection of shapes included in the investigations. Since ice crystal optical properties are mostly determined by the aspect ratios of the crystal components and their microscale structure, it is advised that remote sensing applications focus on the variation of these ice shape characteristics, rather than on the macroscale shape or habit. Recent studies use databases with nearly continuous ranges of crystal component aspect ratio and-or roughness levels to infer the variation of ice crystal shape from satellite and airborne remote sensing measurements. Here, the rationale and results of varying strategies for the remote sensing of ice crystal shape are reviewed. Observed systematic variations of ice crystal geometry with location, cloud height and atmospheric state suggested by the data are discussed. Finally, a prospective is given on the future of the remote sensing of ice cloud particle shapes.

clouds↗

Solar wind mass and momentum flux variations at 0.3 AU

In the past we have used electron Thomson scattering brightness observations, obtained with the zodiacal-light photometers on board the spacecraft Helios 1 and Helios 2, to study the global density structure of the quiet corona and inner heliosphere (greater than 17 solar radii). This was done by means of a comparison of synoptic maps based on these Thomson scattering observations and synoptic maps based on other solar/heliospheric data, such as IPS velocity, K-coronameter brightness and magnetic source surface data. In this paper we continue this approach by combining the Helios Thomson scattering maps (which provide density information) with IPS solar wind velocity maps to map out variations in mass and momentum flux of the solar wind as a function of latitude and phase of the solar cycle. The method used to construct the Helios and IPS synoptic maps emphasizes the global, persistent (as opposed to transient) structures, and thus can be viewed as approximating conditions in the quiet corona and inner heliosphere.

Hick, P.↗

A New Neural Network Approach Including First-Guess for Retrieval of Atmospheric Water Vapor, Cloud Liquid Water Path, Surface Temperature and Emissivities Over Land From Satellite Microwave Observations

The analysis of microwave observations over land to determine atmospheric and surface parameters is still limited due to the complexity of the inverse problem. Neural network techniques have already proved successful as the basis of efficient retrieval methods for non-linear cases, however, first-guess estimates, which are used in variational methods to avoid problems of solution non-uniqueness or other forms of solution irregularity, have up to now not been used with neural network methods. In this study, a neural network approach is developed that uses a first-guess. Conceptual bridges are established between the neural network and variational methods. The new neural method retrieves the surface skin temperature, the integrated water vapor content, the cloud liquid water path and the microwave surface emissivities between 19 and 85 GHz over land from SSM/I observations. The retrieval, in parallel, of all these quantities improves the results for consistency reasons. A data base to train the neural network is calculated with a radiative transfer model and a a global collection of coincident surface and atmospheric parameters extracted from the National Center for Environmental Prediction reanalysis, from the International Satellite Cloud Climatology Project data and from microwave emissivity atlases previously calculated. The results of the neural network inversion are very encouraging. The r.m.s. error of the surface temperature retrieval over the globe is 1.3 K in clear sky conditions and 1.6 K in cloudy scenes. Water vapor is retrieved with a r.m.s. error of 3.8 kg/sq m in clear conditions and 4.9 kg/sq m in cloudy situations. The r.m.s. error in cloud liquid water path is 0.08 kg/sq m . The surface emissivities are retrieved with an accuracy of better than 0.008 in clear conditions and 0.010 in cloudy conditions. Microwave land surface temperature retrieval presents a very attractive complement to the infrared estimates in cloudy areas: time record of land surface temperature will be produced.

Aires, F.↗

Computational Methods for Frictional Contact With Applications to the Space Shuttle Orbiter Nose-Gear Tire

A computational procedure is presented for the solution of frictional contact problems for aircraft tires. A Space Shuttle nose-gear tire is modeled using a two-dimensional laminated anisotropic shell theory which includes the effects of variations in material and geometric parameters, transverse-shear deformation, and geometric nonlinearities. Contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with both contact and friction conditions. The contact-friction algorithm is based on a modified Coulomb friction law. A modified two-field, mixed-variational principle is used to obtain elemental arrays. This modification consists of augmenting the functional of that principle by two terms: the Lagrange multiplier vector associated with normal and tangential node contact-load intensities and a regularization term that is quadratic in the Lagrange multiplier vector. These capabilities and computational features are incorporated into an in-house computer code. Experimental measurements were taken to define the response of the Space Shuttle nose-gear tire to inflation-pressure loads and to inflation-pressure loads combined with normal static loads against a rigid flat plate. These experimental results describe the meridional growth of the tire cross section caused by inflation loading, the static load-deflection characteristics of the tire, the geometry of the tire footprint under static loading conditions, and the normal and tangential load-intensity distributions in the tire footprint for the various static vertical loading conditions. Numerical results were obtained for the Space Shuttle nose-gear tire subjected to inflation pressure loads and combined inflation pressure and contact loads against a rigid flat plate. The experimental measurements and the numerical results are compared.

Tanner, John A.↗

Computational methods for frictional contact with applications to the Space Shuttle orbiter nose-gear tire: Comparisons of experimental measurements and analytical predictions

A computational procedure is presented for the solution of frictional contact problems for aircraft tires. A Space Shuttle nose-gear tire is modeled using a two-dimensional laminated anisotropic shell theory which includes the effects of variations in material and geometric parameters, transverse-shear deformation, and geometric nonlinearities. Contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with both contact and friction conditions. The contact-friction algorithm is based on a modified Coulomb friction law. A modified two-field, mixed-variational principle is used to obtain elemental arrays. This modification consists of augmenting the functional of that principle by two terms: the Lagrange multiplier vector associated with normal and tangential node contact-load intensities and a regularization term that is quadratic in the Lagrange multiplier vector. These capabilities and computational features are incorporated into an in-house computer code. Experimental measurements were taken to define the response of the Space Shuttle nose-gear tire to inflation-pressure loads and to inflation-pressure loads combined with normal static loads against a rigid flat plate. These experimental results describe the meridional growth of the tire cross section caused by inflation loading, the static load-deflection characteristics of the tire, the geometry of the tire footprint under static loading conditions, and the normal and tangential load-intensity distributions in the tire footprint for the various static vertical-loading conditions. Numerical results were obtained for the Space Shuttle nose-gear tire subjected to inflation pressure loads and combined inflation pressure and contact loads against a rigid flat plate. The experimental measurements and the numerical results are compared.

Tanner, John A.↗

Development of an Advanced Annular Combustor

The objective of the effort described in this report was to determine the structural durability of a full-scale advanced annular turbojet combustor using ASTM A-1 type fuel and operating at conditions typical of advanced supersonic aircraft. A full-scale annular combustor of the ram-induction type was fabricated and subjected to a 325-hour cyclic endurance test at conditions representative of operation in a Mach 3.0 aircraft. The combustor exhibited extensive cracking and scoop burning at the end of the test program. But these defects had no appreciable effect on combustor performance, as performance remained at a high level throughout the endurance program. Most performance goals were achieved with pressure loss values near 6% and 8%, and temperature rise variation ratio (deltaTVR) values near 1.25 and l.22 at takeoff and cruise conditions, respectively. Combustion efficiencies approached l004 and the exit radial temperature profiles were approximately as desired.

Rusnak, J. P.↗

Synthesis and dynamic characteristics of large structures with rotating substructures

The variational equations of motion of large structures with rotating substructures are derived by the substructure synthesis approach, whereby a discretization procedure akin to the Rayleigh-Ritz method is used to represent the elastic motion of every substructure by a suitable set of admissible functions. Using an inclusion principle for gyroscopic systems, the effects on the predicted system dynamic characteristics of truncating the number of admissible functions used for each substructure can be assessed. Criteria for rational selection of admissible functions are presented.

Meirovitch, L.↗

High temperature thermal energy storage in steel and sand

The technical and economic potential for high temperature (343 C, 650 F) thermal energy storage in hollow steel ingots, pipes embedded in concrete, and for pipes buried in sand was evaluated. Because it was determined that concrete would separate from pipes due to thermal stresses, concrete was replaced by sand, which is free from thermal stresses. Variations of the steel ingot concept were not cost effective compared to the sand-pipe approach, therefore, the sand-pipe thermal storage unit (TSU) was evaluated in depth to assess the approximate tube spacing requirements consistent with different system performance characteristics and also attendant system costs. For large TSUs which do not require fast response times, the sand-pipe approach offers attractive possibilities. A pipe diameter about 9 cm (3.5 in) and pipe spacing of approximately 25 cm (10 in), with sand filling the interspaces, appears appropriate. Such a TSU system designed for 8 hours charge/discharge cycle has an energy unit storage cost (CE) of $2.63/kWhr-t and a power unit storage cost (Cp) of $42/kW-t (in 1977 dollars).

Turner, R. H.↗

Daylight imagery of ocean surface waves for wave spectra

A simple second-order model is developed to describe the transformation between surface wave slope and radiance incident on an optical sensor. The model includes azimuth as well as elevation variations in the sky-radiance distribution and their effects on the imaging of wave slopes approaching the sensor azimuth from various directions. The model also includes upwelling radiance. Comparison between the exact transformation from wave slope to radiance and the second-order model suggests that they yield quantitatively similar results. The model is used in physical interpretations of how wave slopes are imaged and in indicating a way of optimizing the linearity and contrast of the transfer function relating wave slope to radiance for various sensor geometries, thereby minimizing sources of error. One-dimensional numerical simulations of the imaging of an analytically generated wave-slope profile demonstrate the utility of this technique in obtaining slope spectra.

Monaldo, F. M.↗

Models of the structure of the atmosphere of Venus from the surface to 100 kilometers altitude

From a critical comparison and synthesis of data from the four Pioneer Venus Probes, the Pioneer Venus Orbiter, and the Venera 10, 12, and 13 landers, models of the lower and middle atmosphere of Venus are derived. The models are consistent with the data sets within the measurement uncertainties and established variability of the atmosphere. The models represent the observed variations of state properties with latitude, and preserve the observed static stability. The rationale and the approach used to derive the models are discussed, and the remaining uncertainties are estimated.

Seiff, A.↗

A Design Study of Onboard Navigation and Guidance During Aerocapture at Mars

The navigation and guidance of a high lift-to-drag ratio sample return vehicle during aerocapture at Mars are investigated. Emphasis is placed on integrated systems design, with guidance algorithm synthesis and analysis based on vehicle state and atmospheric density uncertainty estimates provided by the navigation system. The latter utilizes a Kalman filter for state vector estimation, with useful update information obtained through radar altimeter measurements and density altitude measurements based on IMU-measured drag acceleration. A three-phase guidance algorithm, featuring constant bank numeric predictor/corrector atmospheric capture and exit phases and an extended constant altitude cruise phase, is developed to provide controlled capture and depletion of orbital energy, orbital plane control, and exit apoapsis control. Integrated navigation and guidance systems performance are analyzed using a four degree-of-freedom computer simulation. The simulation environment includes an atmospheric density model with spatially correlated perturbations to provide realistic variations over the vehicle trajectory. Navigation filter initial conditions for the analysis are based on planetary approach optical navigation results. Results from a selection of test cases are presented to give insight into systems performance.

Fuhry, Douglas Paul↗

Static multipole polarisabilities and second-order Stark shift in francium

The multipole polarizability of the ground state of francium is calculated by utilizing both the variational technique of Davison and the quantum defect theory underlying the Bates-Damgaard method. This approach is also shown to yield reasonable results for other alkali atoms. Second-order Stark shift for the ground state of francium is presented as a function of field strength for possible future experimental comparison.

Khan, F.↗

Assimilation of satellite surface wind speed data using the GLA analysis/forecast system

Research at the Goddard Laboratory for Atmospheres to evaluate methods to assign directions to both real and simulated surface wind speed data is discussed. Surface wind speed measurements are obtained from satellites, including Seasat, Nimbus-7, and Geosat. The methods include the interpolation of modeled forecast winds to the wind speed datum location, the use of surface pressures with a balance relation, and a variational analysis method. It is found that the best estimates may be obtained using a multipass approach to perform a surface analysis incorporating all conventional surface data.

Bloom, S. C.↗

Eigenvalue sensitivity analysis of planar frames with variable joint and support locations

Two sensitivity equations are derived in this study based upon the continuum approach for eigenvalue sensitivity analysis of planar frame structures with variable joint and support locations. A variational form of an eigenvalue equation is first derived in which all of the quantities are expressed in the local coordinate system attached to each member. Material derivative of this variational equation is then sought to account for changes in member's length and orientation resulting form the perturbation of joint and support locations. Finally, eigenvalue sensitivity equations are formulated in either domain quantities (by the domain method) or boundary quantities (by the boundary method). It is concluded that the sensitivity equation derived by the boundary method is more efficient in computation but less accurate than that of the domain method. Nevertheless, both of them in terms of computational efficiency are superior to the conventional direct differentiation method and the finite difference method.

Chuang, Ching H.↗