Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convex approximation”

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 55 records · Page 3

The Effect of Curvature on the Transition from Laminar to Turbulent Boundary Layer

In the flow over the upper surface of a wing, a discrepancy between the predicted and actual point of transition from laminar to turbulent boundary layer was found. This effect may be due to the comparatively small radius of curvature of the upper surface of the wing. Tests were undertaken to investigate this effect. As far as the authors know, the present investigation is the first to show that curvature has a pronounced effect on the transition of the boundary layer from the laminar to the turbulent state. It appears that three important results have been obtained thus far. First, the experimental points for the convex and concave side of the sheet are consistent with each other and, with the experimental accuracy involved, lie on the same curve. Second, this curve can be approximated by a straight line. Third, the order of magnitude of the variation is such that the curvature ordinarily used on the upper surface of an airplane wing might double the critical Reynolds number.

Clauser, Milton↗

Multiscale Morphological Filtering for Analysis of Noisy and Complex Images

Images acquired with passive sensing techniques suffer from illumination variations and poor local contrasts that create major difficulties in interpretation and identification tasks. On the other hand, images acquired with active sensing techniques based on monochromatic illumination are degraded with speckle noise. Mathematical morphology offers elegant techniques to handle a wide range of image degradation problems. Unlike linear filters, morphological filters do not blur the edges and hence maintain higher image resolution. Their rich mathematical framework facilitates the design and analysis of these filters as well as their hardware implementation. Morphological filters are easier to implement and are more cost effective and efficient than several conventional linear filters. Morphological filters to remove speckle noise while maintaining high resolution and preserving thin image regions that are particularly vulnerable to speckle noise were developed and applied to SAR imagery. These filters used combination of linear (one-dimensional) structuring elements in different (typically four) orientations. Although this approach preserves more details than the simple morphological filters using two-dimensional structuring elements, the limited orientations of one-dimensional elements approximate the fine details of the region boundaries. A more robust filter designed recently overcomes the limitation of the fixed orientations. This filter uses a combination of concave and convex structuring elements. Morphological operators are also useful in extracting features from visible and infrared imagery. A multiresolution image pyramid obtained with successive filtering and a subsampling process aids in the removal of the illumination variations and enhances local contrasts. A morphology-based interpolation scheme was also introduced to reduce intensity discontinuities created in any morphological filtering task. The generality of morphological filtering techniques in extracting information from a wide variety of images obtained with active and passive sensing techniques is discussed. Such techniques are particularly useful in obtaining more information from fusion of complex images by different sensors such as SAR, visible, and infrared.

Kher, A.↗

Geomorphology of Ma'adim Vallis, Mars,and Associated Paleolake Basins

Ma'adim Vallis, one of the largest valleys in the Martian highlands, appears to have originated by catastrophic overflow of a large paleola ke located south of the valley heads. Ma'adim Vallis debouched to Gus ev crater, 900 km to the north, the landing site for the Spirit Mars Exploration Rover. Support for the paleolake overflow hypothesis come s from the following characteristics: (I) With a channel width of 3 km at its head, Ma'adim Vallis originates at two (eastern and western) gaps incised into the divide of the approximately 1.1 M km(exp 2) enc losed Eridania head basin, which suggests a lake as the water source. (2) The sinuous course of Ma'adim Vallis is consistent with overland flow controlled by preexisting surface topography, and structural con trol is not evident or required to explain the valley course. (3) The nearly constant approximately 5 km width of the inner channel through crater rim breaches, the anastomosing course of the wide western tri butary, the migration of the inner channel to the outer margins of be nds in the valley's lower reach, a medial sedimentary bar approximate ly 200 m in height, and a step-pool" sequence are consistent with modeled flows of 1-5 x l0 (exp 6) m(exp 3)/s. Peak discharges were likely higher but are poorly constrained by the relict channel geometry. (4 ) Small direct tributary valleys to Ma'adim Vallis have convex-up lon gitudinal profiles, suggesting a hanging relationship to a valley that was incised quickly relative to the timescales of tributary developm ent. (5) The Eridania basin had adequate volume between the initial d ivide and the incised gap elevations to carve Ma'adim Vallis during a single flood. (6) The Eridania basin is composed of many overlapping , highly degraded and deeply buried impact craters. The floor materials of the six largest craters have an unusually high internal relief ( approximately 1 km) and slope (approximately 0.5-1.5 degrees) among d egraded Martian craters, which are usually flat-floored. Long-term, fluvial sediment transport appears to have been inhibited within these craters, and the topography is inconsistent with basaltic infilling. (7) Fluvial valleys do not dissect the slopes of these deeper crater floor depressions, unlike similar slopes that are dissected at higher levels in the watershed. These characteristics (6, 7) suggest that wa ter mantled at least the lower parts of the Eridania basin floor thro ughout the period of relatively intense erosion early in Martian hist ory. The lake level increased and an overflow occurred near the close of the Noachian (age determined using >5 km crater counts). Initially , the Eridania basin debouched northward at two locations into the in termediate basin, a highly degraded impact crater approximately 500 k m in diameter. As this intermediate basin was temporarily filled with water, erosion took place first along the lower (northern) reach of Ma'adim Vallis, debouching to Gusev crater. The western overflow point was later abandoned, and erosion of the intermediate basin interior was concentrated along the eastern pathway. Subsequent air fall depos ition, impact gardening, tectonism, and limited fluvial erosion modified the Eridania basin region, so evidence for a paleolake is restrict ed to larger landforms that could survive post-Noachian degradation p rocesses.

Irwin, Rossman, P., III↗

The shape of asteroid 1917 Cuyo

Lightcurves obtained for 1917 Cuyo at solar phase angles near 54 degrees have an amplitude delta m = 0.44 mag. However, convex-profile inversion of the lightcurves yields an estimate of the asteroid's mean cross section (C, a 2-D average of the 3-D shape) that is only slightly noncircular, with an elongation approximately 1.15. The estimate of C undoubtedly contains systematic errors, the most severe of which could arise from non-equatorial viewing/illumination geometry. However, Cuyo's radar echo shows very little variation in bandwidth vs. rotation phase, supporting the hypothesis that this asteroid's elongation is rather modest.

Ostro, Steven J.↗

Optical Testing of Diamond Machined, Aspheric Mirrors for Groundbased, Near-IR Astronomy

The Infrared Multi-Object Spectrometer (IRMOS) is a facility-class instrument for the Kitt Peak National Observatory 4 and 2.1 meter telescopes. IRMOS is a near-IR (0.8-2.5 micron) spectrometer and operates at approximately 80 K. The 6061-T651 aluminum bench and mirrors constitute an athermal design. The instrument produces simultaneous spectra at low- to mid-resolving power (R=lambda/delta lambda= 300-3000) of approximately 100 objects in its 2.8 x 2.0 arcmin field. We describe ambient and cryogenic optical testing of the IRMOS mirrors across a broad range in spatial frequency (figure error, mid-frequency error, and microroughness). The mirrors include three rotationally symmetric, off-axis conic sections, one off-axis biconic, and several flat fold mirrors. The symmetric mirrors include convex and concave prolate and oblate ellipsoids. They range in aperture from 94x86 mm to 286x269 mm and in f-number from 0.9 to 2.4. The biconic mirror is concave and has a 94x76 mm aperture, R(sub x)=377 mm, k(sub x)=0.0778, R(sub y)=407 mm, and k(sub y)=0.1265 and is decentered by -2 mm in X and 227 mm in Y. All of the mirrors have an aspect ratio of approximately 6:1. The surface error fabrication tolerances are less than 10 nm RMS microroughness, 'best effort' for mid-frequency error, and less than 63.3 nm RMS figure error. Ambient temperature (approximately 293 K) testing is performed for each of the three surface error regimes, and figure testing is also performed at approximately 80 K. Operation of the ADE Phaseshift MicroXAM white light interferometer (micro-roughness) and the Bauer Model 200 profilometer (mid-frequency error) is described. Both the sag and conic values of the aspheric mirrors make these tests challenging. Figure testing is performed using a Zygo GPI interferometer, custom computer generated holograms (CGH), and optomechanical alignment fiducials. Cryogenic CGH null testing is discussed in detail. We discuss complications such as the change in prescription with temperature and thermal gradients. Correction for the effect of the dewar window is also covered. We discuss the error budget for the optical test and alignment procedure. Data reduction is accomplished using commercial optical design and data analysis software packages. Results from CGH testing at cryogenic temperatures are encouraging thus far.

Chambers, V. John↗

Surface reconstruction from scattered data through pruning of unstructured grids

This paper describes an algorithm for reconstructing a surface from a randomly digitized object. Scan data (treated as a cloud of points) is first tesselated out to its convex hull using Delaunay triangulation. The line-of-sight between each surface point and the scanning device is traversed, and any tetrahedra which are pierced by it are removed. The remaining tetrahedra form an approximate solid model of the scanned object. Due to the inherently limited resolution of any scan, this algorithm requires two additional procedures to produce a smooth, polyhedral surface: one process removes long, narrow tetrahedra which span indentations in the surface between digitized points; the other smooths sharp edges. The results for a moderately resolved sample body and a highly resolved aircraft are displayed.

Maksymiuk, C. M.↗

Experimental data and model for the turbulent boundary layer on a convex, curved surface

Experiments were performed to determine how boundary layer turbulence is affected by strong convex curvature. The data gathered on the behavior of the Reynolds stress suggested the formulation of a simple turbulence model. Data were taken on two separate facilities. Both rigs had flow from a flat surface, over a convex surface with 90 deg of turning and then onto a flat recovery surface. The geometry was adjusted so that, for both rigs, the pressure gradient along the test surface was zero. Two experiments were performed at delta/R approximately 0.10, and one at weaker curvature with delta/R approximately 0.05. Results show that after a sudden introduction of curvature the shear stress in the outer part of the boundary layer is sharply diminished and is even slightly negative near the edge. The wall shear also drops off quickly downstream. When the surface suddenly becomes flat again, the wall shear and shear stress profiles recover very slowly towards flat wall conditions. A simple turbulence model, which was based on the theory that the Prandtl mixing length in the outer layer should scale on the velocity gradient layer, was shown to account for the slow recovery.

Gillis, J. C.↗

Late Cenozoic flexural deformation of the middle U.S. Atlantic passive margin

Despite the century-long recognition of regional epeirogeny along the middle Atlantic passive margin, relatively few studies have focused on understanding postrift uplift mechanisms. Here, we demonstrate that epeirogenic uplift of the central Appalachian Piedmont and subsidence of the Salisbury Embayment represent first-order, flexural isostatic processes driven by continental denudation and offshore deposition. Our results show that regional epeirogenic processes, present on all Atlantic-type passive margins, are best resolved by specific stratigraphic and geomorphic relationships, rather than topography. A simple one-dimensional geodynamic model, constrained by well-dated Baltimore Canyon trough, Coastal Plain, and lower Susquehanna River (piedmont) stratigraphy, simulates flexural deforamtion of the U.S. Atlantic margin. The model represents the passive margin lithosphree as a uniformly thick elastic plate, without horizontal compressive stresses, that deforms flexurally under the stress of strike-averaged, vertically applied line loads. Model results illustrate a complex interaction among margin stratigraphy and geomorphology, the isostatic repsonse to denudational and depositional processes, and the modulating influence of exogenic forces such as eustasy. The current elevation, with respect to modern sea level, of fluvial terraces and correlateive Coastal Plain deposits or unconformities is successfully predicted through the synthesis of paleotopography, eustatic change, and margin flexure. Results suggest that the middle U.S. Atlantic margin landward of East Coast Magnetic Anomaly is underlain by lithoshpere with an average elastic thickness of 40 km (flexural rigidity, D = 4 X 10(exp 23) N m), the margin experience an average, long-term denudation rate of approximately 10m/m.y., and the Piedmont has been flexurally upwaped between 35 and 130 meters in the last 15 m.y. Long term isostatic continental uplift resulting rom denudation and basin subsidence resulting rom sediment loading are accomodated primately by a convex-up flexural hinge, physiographically represented by the Fall Zone. Our results elucidate an inherent danger in using topography alone to constrain late-stage passive margin deformation mechanisms. Only through careful synthesis of field stratigraphic and geomorphic elements such as fluvial terraces, Coastal Plain deposits, and offshore stratigraphy can age control be extended from the offshore depositional setting to the erosionally dominated continent. This sudy demonstrates that despite a relatively subdued topography, the middle U.S. Atlantic margin experiences progressive flexural isostatic deformation similar to that proposed for high-relief margins characterized by great escarpments. Thus margin topographic diversity remains a function of other factors, such as lithospheric composition and/or structure, supracrustal stratigraphy and structure, degree of drainage integration, drainage divide migration and climate.

Pazzaglia, Frank J.↗

Efficient approximation concepts using second order information

The application of second derivative information for solving structural optimization problems is considered. In the present method, rather than building approximate nonlinear forms for the objective function and constraints, only linear approximations are used. A separable quadratic approximation of the Lagrangian function is included in the subproblem statement. The method has been successfully used for simple problems that can be solved in closed form, in addition to the sizing optimization of trusses, and it is shown to converge faster than the convex linearization method or the method of moving asymptotes.

Fleury, Claude↗

Theory, computation, and application of exponential splines

A generalization of the semiclassical cubic spline known in the literature as the exponential spline is discussed. In actuality, the exponential spline represents a continuum of interpolants ranging from the cubic spline to the linear spline. A particular member of this family is uniquely specified by the choice of certain tension parameters. The theoretical underpinnings of the exponential spline are outlined. This development roughly parallels the existing theory for cubic splines. The primary extension lies in the ability of the exponential spline to preserve convexity and monotonicity present in the data. Next, the numerical computation of the exponential spline is discussed. A variety of numerical devices are employed to produce a stable and robust algorithm. An algorithm for the selection of tension parameters that will produce a shape preserving approximant is developed. A sequence of selected curve-fitting examples are presented which clearly demonstrate the advantages of exponential splines over cubic splines.

Mccartin, B. J.↗

Design, Fabrication, Optical Testing, and Performance of Diamond Machined Aspheric Mirrors for Ground-Based Near-IR Astronomy

Challenges in fabrication and testing have historically limited the choice of surfaces available for the design of reflective optical instruments. Spherical and conic mirrors are common, but, for future science instruments, more degrees of freedom are necessary to meet challenging performance and packaging requirements. These instruments will be composed of unusual aspheres located far off-axis with large spherical departure, and some designs will require asymmetric surface profiles. In particular, single-surface astigmatism correction in spectrographs necessitates a toroidal surface, which lacks an axis of rotational symmetry. We describe the design, fabrication, optical testing, and performance of three rotationally symmetric, off-axis, aspheric mirrors and one toroidal, off-axis, biconic camera mirror on aluminum substrates for the Infrared Multi-Object Spectrograph (IRMOS) instrument. IRMOS is a facility instrument for the Kitt Peak National Observatory's Mayall Telescope (3.8 m) and an engineering prototype for a possible design of the Next Generation Space Telescope/Multi-Object Spectrograph. The symmetric mirrors range in aperture from 94x86 mm to 286x269 mm and in f-number from 0.9 to 2.4. They are various off-axis, convex and concave, prolate and oblate ellipsoids. The concave biconic mirror has a 94x76 mm aperture, Rx=377 mm, kx=0.0778, Ry=407 mm, and ky=0.1265 and is decentered. by -2 mm in x and 227 mm in y. The mirrors have an aspect ratio of approximately 4:1. The surface error fabrication tolerances are less than 63.3 nm RMS figure error and less than 10 nm RMS microroughness. The mirrors are attached to the instrument bench via a semi-kinematic, integral flexure mount. We describe mirror design, diamond machining, the results of figure testing using computer-generated holograms, and imaging and scattered light modeling and performance.

Ohl, Raymond G.↗

RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

Mizuno, Yosuke↗

Prediction of large negative shaded-side spacecraft potentials

A calculation by Knott, for the floating potential of a spherically symmetric synchronous-altitude satellite in eclipse, was adapted to provide simple calculations of upper bounds on negative potentials which may be achieved by electrically isolated shaded surfaces on spacecraft in sunlight. Large (approximately 60 percent) increases in predicted negative shaded-side potentials are obtained. To investigate effective potential barrier or angular momentum selection effects due to the presence of less negative sunlit-side or adjacent surface potentials, these expressions were replaced by the ion random current, which is a lower bound for convex surfaces when such effects become very severe. Further large increases in predicted negative potentials were obtained, amounting to a doubling in some cases.

Prokopenko, S. M. L.↗

Development of MLPG and LBIE Methods for Nonlinear Problems of Fracture

The truly meshless local Petrov-Galerkin (MLPG) method holds a great promise in solving boundary value problems, using a local symmetric weak form as a natural approach. In the present paper, in the context of MLPG and the meshless interpolation of a moving least squares (MLS) type, a method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The sub-domains for the shape functions in the MLS approximation are defined only from the primary nodes, and the secondary nodes use the same sub-domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and/or moved without an additional mesh. The essential boundary conditions can be imposed exactly, and non-convex boundaries can be treated without special techniques. Several numerical examples are presented to illustrate the performance of the present method.

Atluri, Satya N.↗

Data reduction using cubic rational B-splines

A geometric method is proposed for fitting rational cubic B-spline curves to data that represent smooth curves including intersection or silhouette lines. The algorithm is based on the convex hull and the variation diminishing properties of Bezier/B-spline curves. The algorithm has the following structure: it tries to fit one Bezier segment to the entire data set and if it is impossible it subdivides the data set and reconsiders the subset. After accepting the subset the algorithm tries to find the longest run of points within a tolerance and then approximates this set with a Bezier cubic segment. The algorithm uses this procedure repeatedly to the rest of the data points until all points are fitted. It is concluded that the algorithm delivers fitting curves which approximate the data with high accuracy even in cases with large tolerances.

Chou, Jin J.↗

Emplacement of multiple flow units on very shallow slopes, east Kawelu Planitia flow field, Venus

Lobate flows with variable radar brightness occur at several locations on Venus. Here we present preliminary mapping results for a lava flow field on the eastern side of Kawelu Planitia, centered at approximately 41 deg N, 251 deg E. Numerous individual flows of varying brightness can be traced for up to 400 km at this location. Megallan image F-MIDRP.40N251;1 contains two major flow complexes as well as a distinct minor flow field associated with a volcanic dome. The southern flow field (unit F) changes from radar dark to radar bright in a convex arc from west to east. The flow path is diverted by a ridge system (unit B) and associated N-S faulted zone (unit A) evident in the Magellan altimetry data. The 130-km-long radar-bright flow field located at 41.5 deg N and 249-250 deg E is comprised of a narrow, channeled flow (unit E) overlying an intermediate brightness flow (unit D), all of which appears to emanate from a solitary dome at 249 deg E. The southern lobe of unit E disappears under a broad leveed channel (41.3 deg N, 250.1 deg E) in unit 5, only to reappear south of the unit 5 boundary. The source of this largest flow field is not evident, but could be associated with the region of small domes west of the F-MIDRP image (centered at 41 deg N, 246 deg E). Another possible source is the fractured area south of unit C. The flow complex is clearly a conglomeration of many flow episodes, not necessarily from a single source.

Helgerud, Michael B.↗

Numerical Estimation of the Curvature of Biological Surfaces

Many biological systems may profitably be studied as surface phenomena. A model consisting of isotropic growth of a curved surface from a flat sheet is assumed. With such a model, the Gaussian curvature of the final surface determines whether growth rate of the surface is subharmonic or superharmonic. These properties correspond to notions of convexity and concavity, and thus to local excess growth and local deficiency of growth. In biological models where the major factors controlling surface growth are intrinsic to the surface, researchers thus gained from geometrical study information on the differential growth undergone by the surface. These ideas were applied to an analysis of the folding of the cerebral cortex, a geometrically rather complex surface growth. A numerical surface curvature technique based on an approximation to the Dupin indicatrix of the surface was developed. A metric for comparing curvature estimates is introduced, and considerable numerical testing indicated the reliability of this technique.

Todd, P. H.↗

Measurements of local convective heat transfer coefficients on ice accretion shapes

The thin-skin heat rate technique was used to determine local convective heat transfer coefficients for four representative ice accretion shapes. The shapes represented three stages of glaze ice formation and one rime ice formation; the ice models had varying degrees of surface roughness. In general, convective heat transfer was higher in regions where the model's surfaces were convex and lower in regions where the surfaces were concave. The effect of roughness was different for the glaze and rime ice shapes. On the glaze ice shapes, roughness increased the maximum Nu by 80 percent, but the other Nu values were virtually unchanged. On the rime ice shape, the Nu numbers near the stagnation point were unchanged. The maximum Nu value increased by 45 percent, and the Nu number downstream of the peak increased by approximately 150 percent.

Arimilli, R. V.↗