Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximation algorithm”

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,225 records · Page 68

Nonlinear effect on modal data analysis method

How the presence of nonlinearities in structural test data can be detected when using modern linear modal data analysis methods is discussed. The extent to which linear algorithms can provide useful information on nonlinear systems was discussed. The modal data analysis method successfully identified the frequency components (but not damping) of the true solution of a nonlinear systems. It was found that the approximate modal parameters of the linear system can be identified by the modal method when low levels of nonlinearities are present. The modal method successfully identified the modal parameters of a linear mode in the presence of a nonlinear response.

Horta, L. G.↗

Multi-grid calculation of three-dimensional transonic potential flows

A multi-grid algorithm has been developed to speed the iterative convergence of calculations for the transonic potential flow past swept wings and wing-fuselage combinations. The method is based upon a fully-conservative, finite-volume approximation to the steady potential equation which is second-order accurate everywhere in the flow field except near shock waves. The multi-grid scheme is incorporated within the framework of an alternating successive-line-overrelaxation (SLOR) solver of the difference equations. Computed results confirm the second-order accuracy of the scheme, and demonstrate the effectiveness of the multi-grid procedure.

Caughey, D. A.↗

Dynamics analysis of electrodynamic satellite tethers. Equations of motion and numerical solution algorithms for the tether

The equations of motion are developed for a perfectly flexible, inelastic tether with a satellite at its extremity. The tether is attached to a space vehicle in orbit. The tether is allowed to possess electrical conductivity. A numerical solution algorithm to provide the motion of the tether and satellite system is presented. The resulting differential equations can be solved by various existing standard numerical integration computer programs. The resulting differential equations allow the introduction of approximations that can lead to analytical, approximate general solutions. The differential equations allow more dynamical insight of the motion.

Nacozy, P. E.↗

A fully Galerkin method for the recovery of stiffness and damping parameters in Euler-Bernoulli beam models

A fully Sinc-Galerkin method for recovering the spatially varying stiffness and damping parameters in Euler-Bernoulli beam models is presented. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which converges exponentially and is valid on the infinite time interval. Hence the method avoids the time-stepping which is characteristic of many of the forward schemes which are used in parameter recovery algorithms. Tikhonov regularization is used to stabilize the resulting inverse problem, and the L-curve method for determining an appropriate value of the regularization parameter is briefly discussed. Numerical examples are given which demonstrate the applicability of the method for both individual and simultaneous recovery of the material parameters.

Smith, R. C.↗

The formation of protostellar disks. I - 1 M(solar)

Hydrodynamical calculations of the collapse of an axisymmetric, rotating one solar mass protostellar cloud, including the effects of radiative transfer and radiative acceleration but without magnetic fields, are presented. The results include calculations of infrared protostellar spectra as a function of time and viewing angle. A numerical algorithm involving explicit nested grids is used to resolve the region of initial disk formation and at the same time to include the outer regions in the calculation. The central part of the protostar is modeled approximately. Initial conditions are systematically varied to investigate their influence on the evolution and final configuration of central star plus circumstellar disk. The initial state for the standard case is a centrally condensed molecular cloud core of one solar mass with a mean density of 8 x 10 exp -18 g/cu cm and a specific angular momentum at the outer edge of 7 x 10 exp 20 sq cm/s. The collapse is followed for 8 x 10 exp 4 yr, at which point 0.45 solar mass is contained in a rapidly rotating central object and most of the remainder in a surrounding equilibrium disk. The stability of this final structure is qualitatively analyzed.

Yorke, Harold W.↗

A six degree of freedom, plume-fuel optimal trajectory planner for spacecraft proximity operations using an A* node search

Spacecraft proximity operations are complicated by the fact that exhaust plume impingement from the reaction control jets of space vehicles can cause structural damage, contamination of sensitive arrays and instruments, or attitude misalignment during docking. The occurrence and effect of jet plume impingement can be reduced by planning approach trajectories with plume effects considered. An A* node search is used to find plume-fuel optimal trajectories through a discretized six dimensional attitude-translation space. A plume cost function which approximates jet plume isopressure envelopes is presented. The function is then applied to find relative costs for predictable 'trajectory altering' firings and unpredictable 'deadbanding' firings. Trajectory altering firings are calculated by running the spacecraft jet selection algorithm and summing the cost contribution from each jet fired. A 'deadbanding effects' function is defined and integrated to determine the potential for deadbanding impingement along candidate trajectories. Plume costs are weighed against fuel costs in finding the optimal solution. A* convergence speed is improved by solving approach trajectory problems in reverse time. Results are obtained on a high fidelity space shuttle/space station simulation. Trajectory following is accomplished by a six degree of freedom autopilot. Trajectories planned with, and without, plume costs are compared in terms of force applied to the target structure.

Jackson, Mark Charles↗

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗

A dual-low-frequency radar for sub-canopy and deep soil-moisture measurements

Measurements of deep and sub-canopy soil moisture are critical in understanding the global water and carbon energy cycle, but are not presently available on a synoptic basis. In this paper, we discuss a proposed spaceborne dual-frequency (UHF and VHF) radar that can provide global & these key measurements. This system is polarimetric and the low transmit frequencies chosen for their penetration abilities necessitate a large antenna that has an aperture of approximately 30m by 11m at VHF, and 30m by 3m at UHF. We describe the mission concept, overall system design and performance characteristics, and discuss ongoing tasks to prototype key system components, and verify the retrieval algorithms. We are also developing a tower-based prototype radar system. This system will, through field observations, demonstrate the scientific effectiveness of the measurement concept and provide critical data for algorithm development. We provide details of the ground experimentation including issues unique to operating at the low-frequencies chosen for these systems.

radar hazard detection↗

Quasi-Global Precipitation as Depicted in the GPCPV2.2 and TMPA V7

After a lengthy incubation period, the year 2012 saw the release of the Global Precipitation Climatology Project (GPCP) Version 2.2 monthly dataset and the TRMM Multi-satellite Precipitation Analysis (TMPA) Version 7. One primary feature of the new data sets is that DMSP SSMIS data are now used, which entailed a great deal of development work to overcome calibration issues. In addition, the GPCP V2.2 included a slight upgrade to the gauge analysis input datasets, particularly over China, while the TMPA V7 saw more-substantial upgrades: 1) The gauge analysis record in Version 6 used the (older) GPCP monitoring product through April 2005 and the CAMS analysis thereafter, which introduced an inhomogeneity. Version 7 uses the Version 6 GPCC Full analysis, switching to the Version 4 Monitoring analysis thereafter. 2) The inhomogeneously processed AMSU record in Version 6 is uniformly processed in Version 7. 3) The TMI and SSMI input data have been upgraded to the GPROF2010 algorithm. The global-change, water cycle, and other user communities are acutely interested in how these data sets compare, as consistency between differently processed, long-term, quasi-global data sets provides some assurance that the statistics computed from them provide a good representation of the atmosphere's behavior. Within resolution differences, the two data sets agree well over land as the gauge data (which tend to dominate the land results) are the same in both. Over ocean the results differ more because the satellite products used for calibration are based on very different algorithms and the dominant input data sets are different. The time series of tropical (30 N-S) ocean average precipitation shows that the TMPA V7 follows the TMI-PR Combined Product calibrator, although running approximately 5% higher on average. The GPCP and TMPA time series are fairly consistent, although the GPCP runs approximately 10% lower than the TMPA, and has a somewhat larger interannual variation. As well, the GPCP and TMPA interannual variations have an apparent phase shift, with GPCP running a few months later. Additional diagnostics will include mean maps and selected scatter plots.

Huffman, George J.↗

ARIEL: Autonomous Excavation Site Selection for Europa Lander Mission Concept

This paper presents ARIEL (Autonomous Ranking and Interrogation of Excavation Location), an autonomy system for selecting an excavation site on-board for NASA’s Europa Lander Mission Concept. Historically, excavation site selection has been performed by a lengthy ground-in-the-loop (GITL) process involving manual inspections, assessments, and decision making in past missions. However, as Europa Lander would have approximately 20 days of lifetime after the landing, many surface activities, including excavation site selection, must be autonomously performed on-board. This paper describes the overall system of ARIEL as well as its two major algorithmic components: vision-based candidate selection and smart interrogation, which estimates the physical properties of the icy surface through physical contact with the robotic arm’s endeffector. Preliminary results are presented using images from Earth analogue sites. The Europa Lander mission returned to the formulation phases in early 2019 while ARIEL was at an early stage of development. Described in this paper is a snapshot of ARIEL as of the project suspension. This paper also describes the remaining challenges to be solved, should the mission resume in the future.

Ono, Hiro↗

Computer Vision Dataset for Aircraft Taxi Operations

The development and democratization of computer vision algorithms are contingent on the availability of high-quality datasets. In this paper, we introduce a database of forward-facing videos from taxiing aircraft as well as the time-correlated flight data at approximately 1 to10 Hz containing aircraft state data and environmental conditions. The video data is sourced from the National Aeronautics and Space Administration Airborne Science Program archive and includes over 33 hours of 4k, 1080p, and 720p video from twenty-two airports around the world. This paper describes the method of the database construction and a brief analysis of its contents.

Ryan Horn↗

GPM Pre-Launch Algorithm Development for Physically-Based Falling Snow Retrievals

In this work we compare and correlate the long time series (Nov.-March) n~easurements of precipitation rate from the Parsivels and 2DVD to the passive (89, 150, 183+/-1, +/-3, +/-7 GHz) observations of NOAA's AMSU-B radiometer. There are approximately 5-8 AMSU-B overpass views of the CARE site a day. We separate the comparisons into categories of no precipitation, liquid rain and falling snow precipitation. Scatterplots between the Parsivel snowfall rates and AMSU-B brightness temperatures (TBs) did not show an exploitable relationship for retrievals. We further compared and contrasted brightness temperatures to other surface measurements such as temperature and relative humidity with equally unsatisfying results. We found that there are similar TBs (especially at 89 and 150 GHz) for cases with falling snow and for non-precipitating cases. The comparisons indicate that surface emissivity contributions to the satellite observed TB over land can add uncertainty in detecting and estimating falling snow. The newest results show that the cloud icc scattering signal in the AMSU-B data call be detected by computing clear air TBs based on CARE radiosonde data and a rough estimate of surface emissivity. That is the differences in computed TI3 and AMSU-B TB for precipitating and nonprecipitating cases are unique such that the precipitating versus ~lon-precipitating cases can be identified. These results require that the radiosonde releases are within an hour of the AMSU-B data and allow for three surface types: no snow on the ground, less than 5 cm snow on the ground, and greater than 5 cm on the ground (as given by ground station data). Forest fraction and measured emissivities were combined to calculate the surface emissivities. The above work and future work to incorporate knowledge about falling snow retrievals into the framework of the expected GPM Bayesian retrievals will be described during this presentation.

Gail Skofronick-Jackson↗

Extending "Deep Blue" Aerosol Retrieval Coverage to Cases of Absorbing Aerosols Above Clouds: Sensitivity Analysis and First Case Studies

Cases of absorbing aerosols above clouds (AACs), such as smoke or mineral dust, are omitted from most routinely processed space-based aerosol optical depth (AOD) data products, including those from the Moderate Resolution Imaging Spectroradiometer (MODIS). This study presents a sensitivity analysis and preliminary algorithm to retrieve above-cloud AOD and liquid cloud optical depth (COD) for AAC cases from MODIS or similar sensors, for incorporation into a future version of the "Deep Blue" AOD data product. Detailed retrieval simulations suggest that these sensors should be able to determine AAC AOD with a typical level of uncertainty approximately 25-50 percent (with lower uncertainties for more strongly absorbing aerosol types) and COD with an uncertainty approximately10-20 percent, if an appropriate aerosol optical model is known beforehand. Errors are larger, particularly if the aerosols are only weakly absorbing, if the aerosol optical properties are not known, and the appropriate model to use must also be retrieved. Actual retrieval errors are also compared to uncertainty envelopes obtained through the optimal estimation (OE) technique; OE-based uncertainties are found to be generally reasonable for COD but larger than actual retrieval errors for AOD, due in part to difficulties in quantifying the degree of spectral correlation of forward model error. The algorithm is also applied to two MODIS scenes (one smoke and one dust) for which near-coincident NASA Ames Airborne Tracking Sun photometer (AATS) data were available to use as a ground truth AOD data source, and found to be in good agreement, demonstrating the validity of the technique with real observations.

aerosols above clouds↗

The role of under-determined approximations in engineering and science application

There is currently a great deal of interest in using response surfaces in the optimization of aircraft performance. The objective function and/or constraint equations involved in these optimization problems may come from numerous disciplines such as structures, aerodynamics, environmental engineering, etc. In each of these disciplines, the mathematical complexity of the governing equations usually dictates that numerical results be obtained from large computer programs such as a finite element method program. Thus, when performing optimization studies, response surfaces are a convenient way of transferring information from the various disciplines to the optimization algorithm as opposed to bringing all the sundry computer programs together in a massive computer code. Response surfaces offer another advantage in the optimization of aircraft structures. A characteristic of these types of optimization problems is that evaluation of the objective function and response equations (referred to as a functional evaluation) can be very expensive in a computational sense. Because of the computational expense in obtaining functional evaluations, the present study was undertaken to investigate under-determinined approximations. An under-determined approximation is one in which there are fewer training pairs (pieces of information about a function) than there are undetermined parameters (coefficients or weights) associated with the approximation. Both polynomial approximations and neural net approximations were examined. Three main example problems were investigated: (1) a function of one design variable was considered; (2) a function of two design variables was considered; and (3) a 35 bar truss with 4 design variables was considered.

Carpenter, William C.↗

Hybird state of the tail mangetic configuration during steady convection events

Previous observations have shown that during periods of steady magnetospheric convection (SMC) a large amount of magnetic flux crosses the plasma sheet (corresponding to approximately 10 deg wide auroral oval at the nightside) and that the magnetic configuration in the midtail is relaxed (the curent sheet is thick and contains enhanced B(sub Z). These signatures are typical for the substorm recovery phase. Using near-geostationary magnetic field data, magnetic field modeling and a noval diagostic technique (isotropic boundary algorithm), we show that in the near-Earth tail the magnetic confirguration is very stretched during the SMC events. This stretching is caused by an intense, thin westward current. Because of the srongly depressed B(sub Z), there is a large radial gradient in the near-tail magetic field. These signatures have been peviously associated only with the substorm growth phase. Our results indicate that during the SMC periods the magnetic configuration is very peculiar, with co-existing thin near-Earth current sheet and thick midtail plasma sheet. The deep local minimum of the equatorial B(sub Z) that devleops at R approximately 12 R(sub E) is consistent with steady, adiabatic, Earthward convection in the midtail. These results impose contraints on the existing substorm theories, and call for an explanation of how such a stressed configuration can persist for such a long time without tail current disruptions that occur at the end of a substorm growth phase.

Sergeev, V. A.↗

The analysis of GEOS-3 altimeter data in the Tasman and Coral seas

A technique was developed for preprocessing GEOS-3 altimetry data to establish a model of the regional sea surface. The algorithms developed models for a 35,000,000 sq km area with an internal precision of + or - 1 m. There were discrepancies between the sea surface model so obtained and GEM6 based geoid profiles with wavelengths of approximately 2500 km and amplitudes of up to 5 m in this region. The amplitudes were smaller when compared with GEM10-based geoid determinations. However, the comparison of 14 pairs of overlapping passes in the region indicated altimeter resolution of the + or - 25 cm level if the wavelength corresponding to the Nyquist frequency were 30 km. The spectral analysis of such comparisons indicated the existence of significant signal strength in the discrepancies after least squares fitting, with wavelengths in excess of 200 km.

Mather, R. S.↗

Computer-aided analysis of nonlinear problems in transport phenomena

The paper describes algorithms for equilibrium and steady-state problems with coefficients in the expansions derived by the Galerkin weighted residual method and calculated from the resulting sets of nonlinear algebraic equations by the Newton-Raphson method. Initial approximations are obtained from nearby solutions by continuation techniques as parameters are varied. The Newton-Raphson technique is preferred because the Jacobian of the solution is useful for continuation, for analyzing the stability of solutions, for detecting bifurcation of solution families, and for computing asymptotic estimates of the effects on any solution of small changes in parameters, boundary conditions, and boundary shape.

Brown, R. A.↗

Estimation of coefficients and boundary parameters in hyperbolic systems

Semi-discrete Galerkin approximation schemes are considered in connection with inverse problems for the estimation of spatially varying coefficients and boundary condition parameters in second order hyperbolic systems typical of those arising in 1-D surface seismic problems. Spline based algorithms are proposed for which theoretical convergence results along with a representative sample of numerical findings are given.

Banks, H. T.↗