Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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 181 records · Page 10

The Ozone Mapping and Profiler Suite (OMPS) Limb Profiler (LP) Version 1 Aerosol Extinction Retrieval Algorithm: Theoretical Basis

The theoretical basis of the Ozone Mapping and Profiler Suite (OMPS) Limb Profiler (LP) Version 1 aerosol extinction retrieval algorithm is presented. The algorithm uses an assumed bimodal lognormal aerosol size distribution to retrieve aerosol extinction profiles at 675 nm from OMPS LP radiance measurements. A first-guess aerosol extinction profile is updated by iteration using the Chahine nonlinear relaxation method, based on comparisons between the measured radiance profile at 675 nm and the radiance profile calculated by the Gauss Seidel limb-scattering (GSLS) radiative transfer model for a spherical-shell atmosphere. This algorithm is discussed in the context of previous limb-scattering aerosol extinction retrieval algorithms, and the most significant error sources are enumerated. The retrieval algorithm is limited primarily by uncertainty about the aerosol phase function. Horizontal variations in aerosol extinction, which violate the spherical-shell atmosphere assumed in the version 1 algorithm, may also limit the quality of the retrieved aerosol extinction profiles significantly.

Loughman, Robert↗

Algorithm for Identifying Erroneous Rain-Gauge Readings

An algorithm analyzes rain-gauge data to identify statistical outliers that could be deemed to be erroneous readings. Heretofore, analyses of this type have been performed in burdensome manual procedures that have involved subjective judgements. Sometimes, the analyses have included computational assistance for detecting values falling outside of arbitrary limits. The analyses have been performed without statistically valid knowledge of the spatial and temporal variations of precipitation within rain events. In contrast, the present algorithm makes it possible to automate such an analysis, makes the analysis objective, takes account of the spatial distribution of rain gauges in conjunction with the statistical nature of spatial variations in rainfall readings, and minimizes the use of arbitrary criteria. The algorithm implements an iterative process that involves nonparametric statistics.

Rickman, Doug↗

Effects of Ice Particle Size Vertical Inhomogeneity on the Passive Remote Sensing of Ice Clouds

The solar reflectance bi-spectral (SRBS) and infrared split-window (IRSpW) methods are two of the most popular techniques for passive ice cloud property retrievals from multispectral imagers. Ice clouds are usually assumed to be vertically homogeneous in global operational algorithms based on these methods, although significant vertical variations of ice particle size are typically observed in ice clouds. In this Study we investigate uncertainties in retrieved optical thickness, effective particle size, and ice water path introduced by a homogeneous cloud assumption in both the SRBS and IRSpW methods, and focus on whether the assumption can lead to significant discrepancies between the two methods. The study simulates the upwelling spectral radiance associated with vertically structured clouds and passes the results through representative SRBS and IRSpW retrieval algorithms. Cloud optical thickness is limited to values for which IRSpW retrievals are possible (optical thickness less than about 7). When the ice cloud is optically thin and yet has a significant ice particle size vertical variation, it is found that both methods tend to underestimate the effective radius and ice water path. The reason for the underestimation is the nonlinear dependence of ice particle scattering properties (extinction and single scattering albedo) on the effective radius. Because the nonlinearity effect is stronger in the IRSpW than the SRBS method, the IRSpW-based IWP tends to be smaller than the SRBS counterpart. When the ice cloud is moderately optically thick, the IRSpW method is relatively insensitive to cloud vertical structure and effective radius retrieval is weighted toward smaller ice particle size, while the weighting function makes the SRBS method more sensitive to the ice particle size in the upper portion of the cloud. As a result, when ice particle size increases monotonically toward cloud base, the two methods are in qualitative agreement; in the event that ice particle size decreases toward cloud base, the effective radius and ice water path retrievals based on the SRBS method are substantially larger than those from the IRSpW. The main findings of this Study Suggest that the homogenous cloud assumption can affect the SRBS and IRSpW methods to different extents and, consequently, can lead to significantly different retrievals. Therefore caution should be taken when comparing and combining the ice cloud property retrievals from these two methods.

Zhang, Zhibo↗

Enhancements to Program LAURA for computation of three-dimensional hypersonic flow

Changes to Program Laura (Langley Aerothermodynamic Upwind Relaxation Algorithm) are presented which enhance both stability and accuracy of the algorithm. A discussion of iteration/sweeping strategies and their relation to computer architectures is included to best exploit the capabilities of serial, vector, and parallel processor machines. Test cases for Mach 10 perfect gas flow and Mach 32 real gas flow in chemical nonequilibrium over a blunt, raked elliptic cone using the thin-layer Navier-Stokes equations are presented in order to demonstrate the current improved capabilities. Algorithm changes include the use of volume averaging, application of a symmetric total variation diminishing (TVD) scheme, and stronger interaction between the grid/shock alignment routine and the relaxation algorithm. Good comparisons with heat transfer and pitching moment data at three different angles of attack for the Mach 10 tests serve to further validate the present algorithm. Parameters are defined which control the coupling of the specie continuity equations with the solution of the mixture conservation equations. A discussion of the consequences involved in the choice of strong versus weak coupling is presented, and a sample nonequilibrium calculation on a fine grid over a full scale model of the Aeroassist Flight Experiment (AFE) demonstrates current capabilities.

Gnoffo, Peter A.↗

Monte Carlo investigation of transient acoustic fields in partially or completely bounded medium

A simple repetitive calculation was used to investigate what happens to the field in terms of the signal paths of disturbances originating from the energy source. The computation allowed the field to be reconstructed as a function of space and time on a statistical basis. The suggested Monte Carlo method is in response to the need for a numerical method to supplement analytical methods of solution which are only valid when the boundaries have simple shapes, rather than for a medium that is bounded. For the analysis, a suitable model was created from which was developed an algorithm for the estimation of acoustic pressure variations in the region under investigation. The validity of the technique was demonstrated by analysis of simple physical models with the aid of a digital computer. The Monte Carlo method is applicable to a medium which is homogeneous and is enclosed by either rectangular or curved boundaries.

Thanedar, B. D.↗

Performance limits for optimal microburst encounter

An effort has been made to ascertain the envelope-edges for uneventful aircraft penetrations of microburst windshears on the basis of optimal aircraft control strategies. Over 1100 such trajectories have been computed for contemporary airliners and general aviation aircraft, in the case of idealized microbursts, using a successive quadratic program trajectory optimization algorithm able to directly handle inequality constraints. Variations of optimal performance with microburst type, intensity, length scale, and location, define performance limits; these limits fall into short, intermediate, and long microburst length scale regimes. The ability to safely transit a microburst also varies strongly with microburst location.

Psiaki, Mark L.↗

COBE experience with filter QUEST

A gyro based filter variation on the standard QUEST attitude determination algorithm is applied to the Cosmic Background Explorer (COBE). Filter QUEST is found to be three times as fast as the batch estimator and slightly more accurate than regular QUEST. Perhaps more important than its speed or accuracy is the fact that Filter QUEST can provide real time attitude solutions when regular QUEST cannot, due to lack of observability. Filter QUEST is also easy to use and adjust for the proper memory length. Suitable applications for Filter QUEST include coarse and real time attitude determination.

Filla, O.↗

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon↗

A comparison of CO(J = 1 yields 0) and CO(J = 2 yields 1) emission in the Milky Way molecular ring

We have carried out a CO (J = 2 right arrow 1) survey of the Scutum Arm region of the Milky Way molecular ring. Our goals are to compare CO (J = 2 right arrow 1) maps of individual Galactic clouds with the large-scale CO (J = 2 right arrow 1) emission from the Galactic plane, and to predict the CO (J = 2 right arrow 1) appearance of a Galactic cloud ensemble in an external galaxy. The angular resolution and spatial coverage of our survey are compatible with the existing CO (J = 1 right arrow 0) survey of this region by Sanders et al., which we use for comparison. We identify 34 molecular clouds in our map region; their relationships between size and line width and between virial mass and luminosity are consistent with the relationships seen in CO (J = 1 right arrow 0) emission. However, we note that previous studies have shown considerable variation in these relationships; we attribute much of this variation to differences in their cloud definition algorithm and statistical method. We find that the median ratio of integrated emission in the two lowest CO transitions for the clouds in our sample is I(2 right arrow 1)/ I(1 right arrow 0) = 0.69, implying that the typical emitting region in the line of sight contains cold gas that is not actively star-forming. Our conclusion that the molecular ring emission is not dominated by star-forming regions is consistent with other large-scale studies of the Milky Way. Our distributions of cloud size and temperature also imply that relatively massive molecular clouds that lac k star formation, such as Maddalena's cloud, are not rare in the inner Galaxy.

Chiar, Jean E.↗

Tiltrotor Vibration Reduction Through Higher Harmonic Control

The results of a joint NASA/Army/Bell Helicopter Textron wind-tunnel test to assess the potential of higher harmonic control (HHC) for reducing vibrations in tiltrotor aircraft operating in the airplane mode of flight, and to evaluate the effectiveness of a Bell-developed HHC algorithm called MAVSS (Multipoint Adaptive Vibration Suppression System) are presented. The test was conducted in the Langley Transonic Dynamics Tunnel using an unpowered 1/5-scale semispan aeroelastic model of the V-22 which was modified to incorporate an HHC system employing both the rotor swashplate and the wing flaperon. The effectiveness of the swashplate and the flaperon acting either singly or in combination in reducing IP and 3P wing vibrations over a wide range of tunnel airspeeds and rotor rotational speeds was demonstrated. The MAVSS algorithm was found to be robust to variations in tunnel airspeed and rotor speed, requiring only occasional on-line recalculations of the system transfer matrix. HHC had only a small (usually beneficial) effect on blade loads but increased pitch link loads by 25%. No degradation in aeroelastic stability was noted for any of the conditions tested.

Nixon, Mark W.↗

Geometric Representations of Condition Queries on Three-Dimensional Vector Fields

Condition queries on distributed data ask where particular conditions are satisfied. It is possible to represent condition queries as geometric objects by plotting field data in various spaces derived from the data, and by selecting loci within these derived spaces which signify the desired conditions. Rather simple geometric partitions of derived spaces can represent complex condition queries because much complexity can be encapsulated in the derived space mapping itself A geometric view of condition queries provides a useful conceptual unification, allowing one to intuitively understand many existing vector field feature detection algorithms -- and to design new ones -- as variations on a common theme. A geometric representation of condition queries also provides a simple and coherent basis for computer implementation, reducing a wide variety of existing and potential vector field feature detection techniques to a few simple geometric operations.

Henze, Chris↗

Utility of BRDF Models for Estimating Optimal View Angles in Classification of Remotely Sensed Images

Statistical classification of remotely sensed images attempts to discriminate between surface cover types on the basis of the spectral response recorded by a sensor. It is well known that surfaces reflect incident radiation as a function of wavelength producing a spectral signature specific to the material under investigation. Multispectral and hyperspectral sensors sample the spectral response over tens and even hundreds of wavelength bands to capture the variation of spectral response with wavelength. Classification algorithms then exploit these differences in spectral response to distinguish between materials of interest. Sensors of this type, however, collect detailed spectral information from one direction (usually nadir); consequently, do not consider the directional nature of reflectance potentially detectable at different sensor view angles. Improvements in sensor technology have resulted in remote sensing platforms capable of detecting reflected energy across wavelengths (spectral signatures) and from multiple view angles (angular signatures) in the fore and aft directions. Sensors of this type include: the moderate resolution imaging spectroradiometer (MODIS), the multiangle imaging spectroradiometer (MISR), and the airborne solid-state array spectroradiometer (ASAS). A goal of this paper, then, is to explore the utility of Bidirectional Reflectance Distribution Function (BRDF) models in the selection of optimal view angles for the classification of remotely sensed images by employing a strategy of searching for the maximum difference between surface BRDFs. After a brief discussion of directional reflect ante in Section 2, attention is directed to the Beard-Maxwell BRDF model and its use in predicting the bidirectional reflectance of a surface. The selection of optimal viewing angles is addressed in Section 3, followed by conclusions and future work in Section 4.

Valdez, P. F.↗

Higher Harmonic Control for Tiltrotor Vibration Reduction

The results of a joint NASA/Army/Bell Helicopter Textron wind-tunnel test to assess the potential of higher harmonic control (HHC) for reducing vibrations in tiltrotor aircraft operating in the airplane mode of flight, and to evaluate the effectiveness of a Bell-developed HHC algorithm called MAVSS (Multipoint Adaptive Vibration Suppression System) are presented. The test was conducted in the Langley Transonic Dynamics Tunnel using an unpowered 1/5- scale semispan aeroelastic model of the V-22 which was modified to incorporate an HHC system employing both the rotor swashplate and the wing flaperon. The effectiveness of the swashplate and the flaperon acting either singly or in combination in reducing 1P and 3P wing vibrations over a wide range of tunnel airspeeds and rotor rotational speeds was demonstrated. The MAVSS algorithm was found to be robust to variations in tunnel airspeed and rotor speed, requiring only occasion-al on-line recalculations of the system transfer matrix.

Nixon, Mark W.↗

Control algorithms for aerobraking in the Martian atmosphere

The Analytic Predictor Corrector (APC) and Energy Controller (EC) atmospheric guidance concepts were adapted to control an interplanetary vehicle aerobraking in the Martian atmosphere. Changes are made to the APC to improve its robustness to density variations. These changes include adaptation of a new exit phase algorithm, an adaptive transition velocity to initiate the exit phase, refinement of the reference dynamic pressure calculation and two improved density estimation techniques. The modified controller with the hybrid density estimation technique is called the Mars Hybrid Predictor Corrector (MHPC), while the modified controller with a polynomial density estimator is called the Mars Predictor Corrector (MPC). A Lyapunov Steepest Descent Controller (LSDC) is adapted to control the vehicle. The LSDC lacked robustness, so a Lyapunov tracking exit phase algorithm is developed to guide the vehicle along a reference trajectory. This algorithm, when using the hybrid density estimation technique to define the reference path, is called the Lyapunov Hybrid Tracking Controller (LHTC). With the polynomial density estimator used to define the reference trajectory, the algorithm is called the Lyapunov Tracking Controller (LTC). These four new controllers are tested using a six degree of freedom computer simulation to evaluate their robustness. The MHPC, MPC, LHTC, and LTC show dramatic improvements in robustness over the APC and EC.

Ward, Donald T.↗

Global and Local Helioseismic Studies of Solar Convection Zone Dynamics Using SOI-MDI on SOHO

Our joint collaborative analyses of global mode data to characterize the solar differential rotation (e.g. Thompson et al. 1996, Schou et al. 1998), and most recently to detect and analyze temporal variations in angular velocity Omega profiles both within the convection zone and in the deeper radiative interior (e.g. Howe et al 2000a,b; Toomre et al. 2000), have led to a series of fascinating discoveries. These should be pursued further as the solar cycle continues. The physical deductions being made from these studies have been greatly strengthened by utilizing both SOI-MDI and GONG data in order to have two independent observational realizations of Doppler images spanning a five-year interval, using two separate procedures to determine global mode splittings, and then analyzing those splitting data sets using both RLS and SOLA inversion procedures. There are considerable subtleties in the effects of instrumental response functions and calibrations, sensitivity of peak finding algorithms and their mode leakage estimates, and stochastic variations in mode amplitudes that can all contribute to apparent changes in the Omega profiles being inferred from sequences of helioseismic data. We have come to understand the implications of many of these calibration and analysis steps, greatly aided by frequent multi-week collaborative working sessions in our Helioseismic Analysis Facility (HAF) at JILA involving many members of the SO1 dynamics and inversion team, including most of our Co-Is during the summer months when we hold intensive working sessions. Considerable further focused attention is required in a collaborative setting on such global mode issues as we continue studying the changing sun.

Toomre, Juri↗

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗

Demand-Capacity Balancing Algorithms for Urban Air Mobility Operations

This paper proposes new Demand/Capacity Balancing (DCB) algorithms that resolve imbalances at enroute waypoints, such as crossing, merging, and UAM corridor entry or exit waypoints, in addition to vertiports; we refer to this algorithm as DCB-Waypoint, still using pre-departure delay as the sole resolution mechanism. Like DCB-Vertiport, DCB-Waypoint takes one flight at a time and resolves imbalances one waypoint at a time, starting from the origin vertiport followed by the sequence of constrained waypoints and, finally, the destination vertiport. The next advancement assigns airborne delays, in addition to pre-departure delays at vertiports. This algorithm, called DCB-Airborne-Delays, uses information on flight speeds to ensure that the assigned delays are feasible, meaning that variations in flight speeds are within feasible aircraft speed ranges. The proposed algorithms are being implemented in a new fast-time simulation tool developed especially for simulating UAM operations. The full paper will provide additional details of the new DCB algorithms that have been developed, as well as the results from simulations using these algorithms. Finally, the DCB algorithms will be compared against one another to obtain insights and recommendations for the future development of more advanced DCB algorithms.

urban air mobility, demand-capacity balancing↗