Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Polarimetric thermal emission from rough surfaces

Recent theoretical works have suggested the potential of passive polarimetry in the remote sensing of geophysical media. It was shown that the third Stokes parameter U of the thermal emission may become larger for azimuthally asymmetric fields of observation. In order to investigate the potential applicability of passive polarimetry to the remote sensing of ocean surface, measurements of the polarimetric thermal emission from a sinusoidal water surface and a numerical study of the polarimetric thermal emission from randomly rough ocean surfaces were performed. Measurements of sinusoidal water surface thermal emission were performed using a sinusoidal water surface which was created by placing a thin sheet of fiberglass with a sinusoidal profile in two dimensions extended infinitely in the third dimension onto a water surface. The theory of thermal emission from a 'two-layer' periodic surface is derived and the exact solution is performed using both the extended boundary condition method (EBC) and the method of moments (MOM). The theoretical predictions are found to be in good agreement with the experimental results once the effects of the radiometer antenna pattern are included and the contribution of background noise to the measurements is modeled. The experimental results show that the U parameter indicates the direction of periodicity of the water surface and can approach values of up to 30 K for the surface observed. Next, a numerical study of polarimetric thermal emission from randomly rough surfaces was performed. A Monte Carlo technique utilizing an exact method for calculating thermal emission was chosen for the study to avoid any of the limitations of the commonly used approximate methods in rough surface scattering. In this Monte Carlo technique, a set of finite rough surface profiles in two dimensions with desired statistics was generated and extended periodically. The polarimetric thermal emission from each surface of the set was then calculated using both the EBC and the MOM and the results were averaged. The surface statistics chosen were intended to model a wind perturbed ocean surface in the X to Ku band microwave region. The results indicate that the U parameter is sensitive to the azimuthal angle between the surface periodicity and the looking angle and to the rms height of the surface, and that the U parameter is fairly insensitive to variations in polar angle, permittivity, surface power law spectrum, and surface spectrum high frequency cutoff. These properties give further strength to the idea of using the U parameter to detect wind direction over the ocean.

Johnson, J. T.↗
NASA NTRS1972

The determination of the satellite orbit of Mariner 9.

This paper presents a comprehensive analysis of the Mars orbital phase of the Mariner 9 trajectory as determined from Earth based radio data. Both the method and accuracy of the orbit determination process are reviewed. Analysis is presented to show the effects of Mars gravity model and node in the plane of the sky errors on the accuracy of orbit determination. In addition the long term evolution of the orbit from insertion to date is presented, and is decomposed into effects from the Mars gravity field, n-body perturbations, and solar radiation pressure. Since the orbit period is nearly commensurable with the Mars rotational period, the orbit experiences significant resonance perturbations. The primary perturbation is in-track with a maximum amplitude of 1000 km and a wavelength of 39 revolutions.

Born, G. H.↗
NASA NTRS1994

Heating of the nighttime D region by very low frequency transmitters

VLF signals propagating in the Earth-ionosphere waveguide are used to probe the heated nighttime D region over three U.S. Navy very low frequency (VLF,3-30 kHz) transmitters. Ionospheric cooling and heating are observed when a transmitter turns off and on in the course of normal operations. Heating by the 24.0-kHz NAA transmitter in Cutler, Maine, (1000 kW radiated power) was observed by this method in 41 of 52 off/on episodes during December 1992, increasing the amplitude and retarding the phase of the 21.4-kHz NSS probe wave propagating from Annapolis, Maryland, to Gander, Newfoundland, by as much as 0.84 dB and 5.3 deg, respectively. In 6 of these 41 episodes, the amplitude of the 28.5-kHz NAU probe wave propagating from Puerto Rico to Gander was also perturbed by as much as 0.29 dB. The latter observations were unexpected due to the greater than 770 km distance between NAA and the NAU-Gander great circle path. Heating by the NSS (21.4 kHz, 265 kW) and NLK (24.8 kHz, 850 kW) transmitters was observed serendipitously in data from earlier measurements of the amplitudes of VLF signals propagating in the Earth-ionosphere waveguide. A three-dimensional model of wave absorption and electron heating in a magnetized, weakly ionized plasma is used to calculate the extent nad shape of the collision frequency (i.e., electron temperature) enhancement above a VLF transmitter. The enhancements are annular, with a geomagnetic north-south asymmetry and a radius at the outer half-maximum of the collision frequency enhancement of about 150 km. Heating by the NAA transmitter is predicted to increase the nighttime D region electron temperature by as much as a factor of 3. The calculated changes in the D region conductivity are used in a three-dimensional model of propagation in the Earth-ionosphere wavelength to predict the effect of the heated patch on a subionospheric VLF probe wave. The range of predicted scattered field amplitudes is in general consistent with the observed signal perturbations. Discrepanices in the predictions are attributed to lack of knowledge of the D region electron density profile along the probe wave great circle paths.

Rodriguez, Juan V.↗
NASA NTRS2004

Alternative Attitude Commanding and Control for Precise Spacecraft Landing

A report proposes an alternative method of control for precision landing on a remote planet. In the traditional method, the attitude of a spacecraft is required to track a commanded translational acceleration vector, which is generated at each time step by solving a two-point boundary value problem. No requirement of continuity is imposed on the acceleration. The translational acceleration does not necessarily vary smoothly. Tracking of a non-smooth acceleration causes the vehicle attitude to exhibit undesirable transients and poor pointing stability behavior. In the alternative method, the two-point boundary value problem is not solved at each time step. A smooth reference position profile is computed. The profile is recomputed only when the control errors get sufficiently large. The nominal attitude is still required to track the smooth reference acceleration command. A steering logic is proposed that controls the position and velocity errors about the reference profile by perturbing the attitude slightly about the nominal attitude. The overall pointing behavior is therefore smooth, greatly reducing the degree of pointing instability.

Singh, Gurkirpal↗
NASA NTRSDate not supplied

Lapse-Rate Feeback: The Key to Reconciling TOA and Surface Attributions of Surface Warming

The impact of climate feedbacks on surface warming is conventionally evaluated using a decomposition of the top-of-atmosphere (TOA) energy budget. Alternatively, the climate feedback analysis can also be carried out using the surface energy budget. However, the two perspectives do not provide the same interpretation of process contributions to surface warming, particularly when executing a spatial analysis. The TOA energy budget is equal to the sum of the surface and atmospheric energy budgets. Using the CMIP5 RCP 8.5 model projections, we show that the major discrepancies between the surface and TOA climate feedback attributions of surface warming are due to non-negligible changes in the atmospheric energy budget that differ from their counterparts at the surface. Individual radiative and non-radiative processes cause vertically non-uniform energy flux perturbations, in response to an external forcing, that naturally lead to a vertically non-uniform temperature change. The TOA lapse-rate feedback is the manifestation of these multiple processes that produce a vertically non-uniform warming response such that it accounts for the asymmetry between the changes in the atmospheric and surface energy budgets. Using the climate feedback-response analysis method, we can decompose the lapse-rate feedback into contributions by individual processes. The negative lapse-rate feedback in the tropics is due to greater moist convection and condensational heating along with ocean heat storage. On the other hand, the positive lapse-rate feedback in polar regions is primarily due to the surface albedo and water vapor feedbacks. Combining the process contributions that are hidden within the lapse-rate feedback with their respective direct impacts on the TOA energy budget allows for a very consistent picture of process contributions to surface warming and its inter-model spread as that given by the surface energy budget approach. With the modified approach, both perspectives indicate water vapor feedback is the largest contributor in the tropics, while surface albedo feedback is the greatest contributor in polar regions. Dynamics and ocean heat storage are the main suppressors of surface warming, except over Antarctica. Both perspectives show there is large inter-model uncertainty in the contributions of water vapor, clouds, albedo, and dynamics plus ocean heat storage.

Sergio A Sejas↗
NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗
NASA NTRS1985

A remark about pointed bubbles

The polymer expansion is a formal algebraic identity between a partition function and logarithm in statistical physics problems. The expansion gives a systematic method to control the free energy or to establish exponential tree-graph decay of connected correlations. Here, the convergence properties of the polymer expansion are analyzed in connection with three practical examples, including: intersecting bonds in chemical polymer chains; a connected closed hypersurface built from the (d-1)-faces of the d-dimensional unit cubes; and the set of Feynamn diagrams in the perturbation series of the Euclidean field theory partition function Z. The example of connected polymer chains is generalized to apply to other lattice models, including n-state Ising models at high temperature; short range lattice gases at high temperature; and weak coupling lattice field and gauge theories.

Garabedian, P. R.↗
NASA NTRS1987

Stability and resonance in grooved-channel flows

Numerical simulation using the spectral element method is used to study the stability and resonant response of incompressible moderate-Reynolds-number flows in periodically-grooved channels. For Reynolds numbers below a critical value, the flow is shown to approach a stable steady-state, and the least stable modes resemble Tollmien-Schlichting channel waves. For Reynolds numbers greater than the critical value, self-sustained oscillations result. Oscillatory perturbation of the grooved-channel flow at the frequency of the least stable mode of the linearized system is found to result in subcritical resonant excitation as the critical Reynolds number is approached. The importance of nonhomogeneous geometry in the forced response of the flow is discussed.

Ghaddar, Nasreen K.↗
NASA NTRS1992

Methodology for the passive control of orbital inclination and mean local time to meet sun-synchronous orbit requirements

The mean local time (MLT) of equatorial crossing of a sun-synchronous Earth-observing spacecraft orbit drifts with inclination; therefore, in order to maintain the MLT, the inclination must be controlled. Inclination may be maintained actively by costly out-of-plane maneuvers or passively by using the perturbing forces due to the sun and moon. This paper examines the passive control approach using the Earth Observing System (EOS) as a basis for the discussion. Applications to Landsat and National Oceanic and Atmospheric Administration (NOAA) spacecraft are presented for comparison. This technique is especially beneficial to spacecraft lacking propulsion systems. The results indicate that passive inclination control appears to be the preferable maintenance method when spacecraft weight restrictions, operational considerations, and scientific requirements apply.

Folta, David↗
NASA NTRS2010

Transport and Use of a Centaur Second Stage in Space

As nations continue to explore space, the desire to reduce costs will continue to grow. As a method of cost reduction, transporting and/or use of launch system components as integral components of missions may become more commonplace in the future. There have been numerous scenarios written for using launch vehicle components (primarily space shuttle used external tanks) as part of flight missions or future habitats. Future studies for possible uses of launch vehicle upper stages might include asteroid diverter using gravity orbital perturbation, orbiting station component, raw material at an outpost, and kinetic impactor. The LCROSS (Lunar CRater Observation and Sensing Satellite) mission was conceived as a low-cost means of determining whether water exists at the polar regions of the moon. Manifested as a secondary payload with the LRO (Lunar Reconnaissance Orbiter) spacecraft aboard an Atlas V launch vehicle, LCROSS guided its spent Centaur Earth Departure Upper Stage (EDUS) into the lunar crater Cabeu's, as a kinetic impactor. This paper describes some of the challenges that the LCROSS project encountered in planning, designing, launching with and carrying the Centaur upper stage to the moon.

Strong, James M.↗
NASA NTRS1988

A perturbation technique for parameter identification in distributed structures

Structures are often characterized by parameters, such as mass and stiffness, that are spatially distributed. Parameter identification of distributed structures is subject to many of the difficulties involved in the modeling problem, and the choice of the model can greatly affect the results of the parameter identification process. Analogously to control spillover in the control of distributed-parameter systems, identification spillover is shown to exist as well and its effect is to degrade the parameter estimates. Moreover, as in modeling by the Rayleigh-Ritz method, it is shown that, for a Rayleigh-Ritz type identification algorithm, an inclusion principle exists in the identification of distributed-parameter systems as well, so that the identified natural frequencies approach the actual natural frequencies monotonically from above.

Meirovitch, L.↗
NASA NTRS1985

Use of ground radar to detect reentering debris

The velocity of the particles is required to identify the type of particles producing the ionization trails. A method of approximating the velocity of a meteor from radar data was developed. The method requires the time between the spacings of the Fresnel interference fringes, the range to the ionization trail, and the wavelength of the radar system. The orbital mechanics of the problem are evaluated, if the particles originate with the shuttle, the orbital mechanics will substantiate the relative position of the particles with the position of the shuttle. A program to determine spacecraft orbital decay due to perturbations is utilized for a preliminary evaluation of the orbital mechanics of the problem. Many assumptions concerning the size, shape, density, etc. of the particles are necessary for the preliminary evaluation. The results do not negate the possibility that the events observed by the radar are reentering particles originating from the shuttle.

Crews, J. L.↗
NASA NTRS1993

Linear, multivariable robust control with a mu perspective

The structured singular value is a linear algebra tool developed to study a particular class of matrix perturbation problems arising in robust feedback control of multivariable systems. These perturbations are called linear fractional, and are a natural way to model many types of uncertainty in linear systems, including state-space parameter uncertainty, multiplicative and additive unmodeled dynamics uncertainty, and coprime factor and gap metric uncertainty. The structured singular value theory provides a natural extension of classical SISO robustness measures and concepts to MIMO systems. The structured singular value analysis, coupled with approximate synthesis methods, make it possible to study the tradeoff between performance and uncertainty that occurs in all feedback systems. In MIMO systems, the complexity of the spatial interactions in the loop gains make it difficult to heuristically quantify the tradeoffs that must occur. This paper examines the role played by the structured singular value (and its computable bounds) in answering these questions, as well as its role in the general robust, multivariable control analysis and design problem.

Packard, Andy↗
NASA NTRS2004

Fracture and Stress Evolution on Europa: New Insights Into Fracture Interpretation and Ice Thickness Estimates Using Fracture Mechanics Analyses

The work completed during the funding period has provided many important insights into fracturing behavior in Europa's ice shell. It has been determined that fracturing through time is likely to have been controlled by the effects of nonsynchronous rotation stresses and that as much as 720 deg of said rotation may have occurred during the visible geologic history. It has been determined that there are at least two distinct styles of strike-slip faulting and that their mutual evolutionary styles are likely to have been different, with one involving a significant dilational component during shear motion. It has been determined that secondary fracturing in perturbed stress fields adjacent to older structures such as faults is a prevalent process on Europa. It has been determined that cycloidal ridges are likely to experience shear stresses along the existing segment portions as they propagate, which affects propagation direction and ultimately induces tailcracking at the segment tip than then initiates a new cycle of cycloid segment growth. Finally, it has been established that mechanical methods (e.g., flexure analysis) can be used to determine the elastic thickness of the ice shell, which, although probably only several km thick, is likely to be spatially variable, being thinner under bands but thicker under ridged plains terrain.

Kattenhorn, Simon↗
NASA NTRS1976

A new unified approach for analyzing wing-body-tail configurations with control surfaces

A general theory for steady and unsteady, subsonic and supersonic potential aerodynamics for complex configurations is presented. Special attention is given to the theoretical formulation and the corresponding numerical implementation for coplanar interfering surfaces. Applying the Green's function method to the equation of the velocity potential and discretizing the spatial problem by using the finite-element technique, yields a set of differential-delay equations in time relating the potential to the normal wash. For fully unsteady flow, the motion is assumed to consist of constant subsonic or supersonic speed for time t less than or equal to 0 (steady state) and of small perturbations around the steady state for time t greater than 0; the solution is obtained in the Laplace domain. From the potential, the aerodynamic pressure and the generalized forces are evaluated. The program SOUSSA (Steady, Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) is briefly described. Numerical results obtained with SOUSSA are presented.

Tseng, K.↗
NASA NTRS1980

Reformulation of Possio's kernel with application to unsteady wind tunnel interference

An efficient method for computing the Possio kernel has remained elusive up to the present time. In this paper the Possio is reformulated so that it can be computed accurately using existing high precision numerical quadrature techniques. Convergence to the correct values is demonstrated and optimization of the integration procedures is discussed. Since more general kernels such as those associated with unsteady flows in ventilated wind tunnels are analytic perturbations of the Possio free air kernel, a more accurate evaluation of their collocation matrices results with an exponential improvement in convergence. An application to predicting frequency response of an airfoil-trailing edge control system in a wind tunnel compared with that in free air is given showing strong interference effects.

Fromme, J. A.↗
NASA NTRS1982

Time-marching transonic flutter solutions including angle-of-attack effects

Transonic aeroelastic solutions based upon the transonic small perturbation potential equation were studied. Time-marching transient solutions of plunging and pitching airfoils were analyzed using a complex exponential modal identification technique, and seven alternative integration techniques for the structural equations were evaluated. The HYTRAN2 code was used to determine transonic flutter boundaries versus Mach number and angle-of-attack for NACA 64A010 and MBB A-3 airfoils. In the code, a monotone differencing method, which eliminates leading edge expansion shocks, is used to solve the potential equation. When the effect of static pitching moment upon the angle-of-attack is included, the MBB A-3 airfoil can have multiple flutter speeds at a given Mach number.

Edwards, J. W.↗
NASA NTRS1982

Time-marching transonic flutter solutions including angle-of-attack effects

Transonic aeroelastic solutions based upon the transonic small perturbation potential equation are studied. Time-marching transient solutions of plunging and pitching airfoils are analyzed using a complex exponential modal identification technique, and seven alternative integration techniques for the structural equations are evaluated. The HYTRAN2 code is used to determine transonic flutter boundaries versus Mach number and angle-of-attack for NACA 64A010 and MBB A-3 airfoils. In that code, a monotone differencing method, which eliminates leading edge expansion shocks, is used to solve the potential equation. When the effect of static pitching moment upon the angle-of-attack is included, the MBB A-3 airfoil can have multiple flutter speeds at a given Mach number.

Edwards, J. W.↗