Search NASA⌕ Search

SEARCH · Search NASA

Results for “Smoothing 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 253 records · Page 14

Interface of an uncoupled boundary layer algorithm with an inviscid core flow algorithm for unsteady supersonic engine inlets

An uncoupled boundary layer algorithm was combined with an inviscid core flow algorithm to model flows within supersonic engine inlets. The inviscid flow algorithm that was used was the LArge Perturbation INlet Code (LAPIN). The boundary layer and inviscid core flow algorithms were formulated in different manners. The boundary layer algorithm was two dimensional and solved in nonconservation form, while the core flow algorithm was one dimensional and solved in conservation form. In order to interface the two codes, the following modifications were important. The coordinate system was set up to maintain the parabolic nature of the boundary layer algorithm while approaching the one dimensional core flow solution far from a wall. The pressure gradient used in the boundary layer equation was calculated using the core flow values and the boundary layer equations, so the boundary layer solution smoothly approached the core flow values far from the wall. Flaring was used for the advection terms perpendicular to the core flow to maintain the stability of the algorithm. With these modifications, the combined viscous/inviscid algorithm matched well with experimental observations of pressure distributions with a supersonic inlet.

Darling, Douglas↗

Interface of an uncoupled boundary layer algorithm with an inviscid core flow algorithm for unsteady supersonic engine inlets

An uncoupled boundary layer algorithm was combined with an inviscid core flow algorithm to model flows within supersonic engine inlets. The inviscid flow algorithm that was used was the LArge Perturbation INlet Code (LAPIN). The boundary layer and inviscid core flow algorithms were formulated in different manners. The boundary layer algorithm was two dimensional and solved in nonconservation form, while the core flow algorithm was one dimensional and solved in conservation form. In order to interface the two codes, the following modifications were important. The coordinate system was set up to maintain the parabolic nature of the boundary layer algorithm while approaching the one dimensional core flow solution far from a wall. The pressure gradient used in the boundary layer equation was calculated using the core flow values and the boundary layer equations, so the boundary layer solution smoothly approached the core flow values far from the wall. Flaring was used for the advection terms perpendicular to the core flow to maintain the stability of the algorithm. With these modifications, the combined viscous/inviscid algorithm matched well experimental observations of pressure distributions with a supersonic inlet.

Darling, Douglas↗

Optimization of satellite altimeter and wave height measurements

Two techniques for simultaneously estimating altitude, ocean wave height, and signal-to-noise ratio from the GEOS-C satellite altimeter data are described. One technique was based on maximum likelihood estimation, MLE, and the other on minimum mean square error estimation, MMSE. Performance was determined by comparing the variance and bias of each technique with the variance and bias of the smoothed output from the Geos altimeter tracker. Ocean wave height tracking performance for the MLE and MMSE algorithms was measured by comparing the variance and bias of the wave height estimates with that of the expression for the return waveform obtained by a fit to the average output of the 16 waveform sampling gates.

Dooley, R. P.↗

Transonic flow calculations using a flux vector splitting method for the Euler equations

A study of the flux vector splitting method of Steger and Warming for the solution of the time dependent Euler equations in strong conservation law form for arbitrary two-dimensional geometries is presented. The procedure employed here differs from that of Buning and Steger in that it uses a different algorithm and employs implicit boundary conditions. Moreover, the method, as implemented here, does not contain any explicit smoothing or any adjustable parameters. Calculations were carried out for an NACA 0012 airfoil at various Mach numbers and angles of attack, and cylinders. Steady symmetric solutions were obtained for the full cylinder at a freestream Mach number of .5 without imposing a symmetry condition. In general, good agreement with other methods was obtained.

Seaford, C. M.↗

LANDSAT instruments characterization

Work performed for the LANDSAT instrument characterization task in the areas of absolute radiometry, coherent noise analysis, and between-date smoothing is reported. Absolute radiometric calibration for LANDSAT-5 TM under ambient conditions was performed. The TM Radiometric Algorithms and Performance Program (TRAPP) was modified to create optional midscan data files and to match the TM Image Processing System (TIPS) algorithm for pulse determination. Several data reduction programs were developed, including a linear regression and its plotted result. A fast Fourier transformation study was conducted on the resequenced TM data. Subscenes of homogeneous water within scenes over Pensacola, Florida were used for testing the FFT on the resequenced data. Finally, a gain and pulse height stability study of LANDSAT 5 TM spectral bands was performed.

Lee, Y.↗

The MSFC/J70 orbital atmosphere model and the data bases for the MSFC solar activity prediction technique

The MSFC/J70 Orbital Atmospheric Density Model, a modified version of the Smithsonian Astrophysical Observatory Jacchia 1970 model is explained. The algorithms describing the MSFC/J70 model are included as well as listing of the computer program. The 13-month smoothed values of solar flux (F sub 10.7) and geomagnetic index (S sub p), which are required as inputs for the MSFC/J70 model, are also included and discussed.

Johnson, D. L.↗

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. 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.

Anderson, W. Kyle↗

Volume 2: Explicit, multistage upwind schemes for Euler and Navier-Stokes equations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using monotonic upstream schemes for conservation laws (MUSCL)-type differencing to obtain state variables at cell interface. Variable interpolations were written in the k-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator. Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-state predictor-corrector schemes, and multistage time-stepping schemes. The coefficients of the multistage time-stepping schemes have been modified successfully to achieve better performance with upwind differencing. A technique was developed to optimize the coefficients for good high-frequency damping at relatively high CFL numbers. Local time-stepping, implicit residual smoothing, and multigrid procedure were added to the explicit time stepping scheme to accelerate convergence to steady-state. The developed algorithm was implemented successfully in a multi-block code, which provides complete topological and geometric flexibility. The only requirement is C degree continuity of the grid across the block interface. The algorithm has been validated on a diverse set of three-dimensional test cases of increasing complexity. The cases studied were: (1) supersonic corner flow; (2) supersonic plume flow; (3) laminar and turbulent flow over a flat plate; (4) transonic flow over an ONERA M6 wing; and (5) unsteady flow of a compressible jet impinging on a ground plane (with and without cross flow). The emphasis of the test cases was validation of code, and assessment of performance, as well as demonstration of flexibility.

Elmiligui, Alaa↗

Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields

Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.

Tannehill, J. C.↗

The topology of large-scale structure. V - Two-dimensional topology of sky maps

A 2D algorithm is applied to observed sky maps and numerical simulations. It is found that when topology is studied on smoothing scales larger than the correlation length, the topology is approximately in agreement with the random phase formula for the 2D genus-threshold density relation, G2(nu) varies as nu(e) exp-nu-squared/2. Some samples show small 'meatball shifts' similar to those seen in corresponding 3D observational samples and similar to those produced by biasing in cold dark matter simulations. The observational results are thus consistent with the standard model in which the structure in the universe today has grown from small fluctuations caused by random quantum noise in the early universe.

Gott, J. R., III↗

Satellite laser station Helwan status 1992

The Satellite Laser Station Helwan has been operated jointly by the National Research Institute of Astronomy and Geophysics in Helwan, Egypt and the Czech Technical University in Prague, Czechslovakia. The station components have been carefully tuned to increase the systems overall stability and reliability critical for the remote location. The mount correction model based on the Gaussian smoothing has been implemented to simplify the blind satellite acquisition and tracking. The on-site normal points generation algorithm has been implemented, the station has been connected to the international information network. The ERS-1 satellite has been included into the tracking schedule. The station range capability has been verified by experimental Etalon 1 ranging by April 1992. The ranging precision of 2-3 centimeters is obtained when ranging to ERS-1, Starlette, and Lageos satellites.

Cech, M.↗

A Full Multi-Grid Method for the Solution of the Cell Vertex Finite Volume Cauchy-Riemann Equations

The system of inhomogeneous Cauchy-Riemann equations defined on a square domain and subject to Dirichlet boundary conditions is considered. This problem is discretised by using the cell vertex finite volume method on quadrilateral meshes. The resulting algebraic problem is overdetermined and the solution is defined in a least squares sense. By this approach a consistent algebraic problem is obtained which differs from the original one by O(h(exp 2)) perturbations of the right-hand side. A suitable cell-based convergent smoothing iteration is presented which is naturally linked to the least squares formulation. Hence, a standard multi-grid algorithm is reported which combines the given smoother and cell-based transfer operators. Some remarkable reduction properties of these operators are shown. A full multi-grid method is discussed which solves the discrete problem to the level of truncation error by employing one multi-grid cycle at each current level of discretisation. Experiments and applications of the full multi-grid scheme are presented.

Borzi, A.↗

Visualization and 3D Mesh Generation for ASSEMBLERS Range of Motion

In an effort to create a base and foundation for ASSEMBLERS with an informed layout, data representing the range of motion of varying heights of Stewart platforms were analyzed. 3D shapes were generated first to be approximately representative of these data in the form of rotationally symmetric shells, making use of 2D alpha shapes to create a border that did not account for asymmetry inherent in the true range of motion. Further work was performed to generate mesh files more fully representative of reach, creating 3D alpha shapes based directly on point clouds generated from uniformly random platform positions. These shapes were initially not representative of the entire range of motion of a platform stack due to the low likelihood of the most extreme positions being added to a given point cloud but were more representative of the asymmetry inherent in the possible positions of a Stewart platform as they no longer relied on a single 2D contour. The same algorithm for generating 3D alpha shapes was then applied to point clouds designed to represent extreme positions, generating larger smooth shells with six clear planes of symmetry. All shells were converted to STLs to aid in further modeling of range of motion for base layouts and in visualizing reachable area for varying stack heights, successfully creating closed shells representative of typical and extreme positions for platform stacks of varying heights.

Jonah DeGuire Vanke↗

An efficient and robust grid optimization algorithm

The development of an efficient and robust grid optimization is presented. This algorithm is developed by combining the best characteristics of algebraic, elliptic, and hyperbolic grid generation techniques. This development is based on the following observations and evaluations: (1) algebraic systems are fast and economical; (2) precise spacing control is always achieved; (3) grid generation by elliptic systems is always smooth; and (4) the hyperbolic system preserves the orthogonality at the solid boundary and the point distribution in the field. Computational examples representing practical internal flow configurations are presented to demonstrate the algorithm.

Soni, Bharat K.↗

Efficient computation of optimal low thrust gravity perturbed orbit transfers

We have developed a new method for solving low-thrust fuel-optimal orbit transfer problems in the vicinity of a large body (planet or asteroid), considering a high-fidelity spherical harmonic gravity model. The algorithm is formulated via the indirect optimization method, leading to a two-point boundary value problem (TPBVP). We make use of a hyperbolic tangent smoothing law for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and thus promote convergence. The TPBVP is solved using the method of particular solutions (MPS) shooting method and Picard-Chebyshev numerical integration. Application of Picard-Chebyshev integration affords an avenue for increased efficiency that is not available with step-by-step integrators. We demonstrate that computing the particular solutions with only a low-fidelity force model greatly increases the efficiency of the algorithm while ultimately achieving near machine precision accuracy. A salient feature of the MPS is that it is parallelizable, and thus further speedups are available. It is also shown that, for near Earth orbits and over a small number of en-route revolutions around the Earth, only the zonal perturbation terms are required in the costate equations to obtain a solution that is accurate to machine precision and optimal to engineering precision. The proposed framework can be used for trajectory design around small asteroids and also for orbit debris rendezvous and removal tasks.

Junkins, John L.↗

A global atlas of GEOS-3 significant waveheight data and comparison of the data with national buoy data

The accuracy of the GEOS-3 significant waveheight estimates compared with buoy measurements of significant waveheight were determined. A global atlas of the GEOS-3 significant waveheight estimates gathered is presented. The GEOS-3 significant waveheight estimation algorithm is derived by analyzing the return waveform characteristics of the altimeter. Convergence considerations are examined, the rationale for a smoothing technique is presented and the convergence characteristics of the smoothed estimate are discussed. The GEOS-3 data are selected for comparison with buoy measurements. The GEOS-3 significant waveheight estimates are assembled in the form of a global atlas of contour maps. Both high and low sea state contour maps are presented, and the data are displayed both by seasons and for the entire duration of the GEOS-3 mission.

Mcmillan, J. D.↗

A new approach to remove interorder background in high dispersion IUE images

A two dimensional Chebyshev interpolating scheme, followed by modeling of point spread functions as Voigt profiles, is proposed for final archiving IUESIPS to remove the interorder overlap from the smoothed background fluxes of high dispersion International Ultraviolet Explorer (IUE) images. Tests of the first and last stages of the algorithm suggest that it can give a considerably better zero point than the old IUESIPS algorithm which samples local interorder background.

Smith, Myron A.↗

Potential, velocity, and density fields from redshift-distance samples: Application - Cosmography within 6000 kilometers per second

A potential flow reconstruction algorithm has been applied to the real universe to reconstruct the three-dimensional potential, velocity, and mass density fields smoothed on large scales. The results are shown as maps of these fields, revealing the three-dimensional structure within 6000 km/s distance from the Local Group. The dominant structure is an extended deep potential well in the Hydra-Centaurus region, stretching across the Galactic plane toward Pavo, broadly confirming the Great Attractor (GA) model of Lynden-Bell et al. (1988). The Local Supercluster appears to be an extended ridge on the near flank of the GA, proceeding through the Virgo Southern Extension to the Virgo and Ursa Major clusters. The Virgo cluster and the Local Group are both falling toward the bottom of the GA potential well with peculiar velocities of 658 + or - 121 km/s and 565 + or - 125 km/s, respectively.

Bertschinger, Edmund↗