Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape 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 235 records · Page 13

Computational aspects of sensitivity calculations in linear transient structural analysis

A study was performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal of the study was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semi-analytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. In several cases this fixed mode approach resulted in very poor approximations of the stress sensitivities. Almost all of the original modes were required for an accurate sensitivity and for small numbers of modes, the accuracy was extremely poor. To overcome this poor accuracy, two semi-analytical techniques were developed. The first technique accounts for the change in eigenvectors through approximate eigenvector derivatives. The second technique applies the mode acceleration method of transient analysis to the sensitivity calculations. Both result in accurate values of the stress sensitivities with a small number of modes and much lower computational costs than if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, William H.↗

Linear-phase approximation in the triangular facet near-field physical optics computer program

Analyses of reflector antenna surfaces use a computer program based on a discrete approximation of the radiation integral. The calculation replaces the actual surface with a triangular facet representation; the physical optics current is assumed to be constant over each facet. Described here is a method of calculation using linear-phase approximation of the surface currents of parabolas, ellipses, and shaped subreflectors and compares results with a previous program that used a constant-phase approximation of the triangular facets. The results show that the linear-phase approximation is a significant improvement over the constant-phase approximation, and enables computation of 100 to 1,000 lambda reflectors within a reasonable time on a Cray computer.

Imbriale, W. A.↗

On the Retrieval of Ice Cloud Particle Shapes from POLDER Measurements

Shapes of ice crystals can significantly affect the radiative transfer in ice clouds. The angular distribution of the polarized reflectance over ice clouds strongly depends on ice crystal shapes. Although the angular-distribution features of the total or polarized reflectance over ice clouds implies a possibility of retrieving ice cloud particle shapes by use of remote sensing data, the accuracy of the retrieval must be evaluated. In this study, a technique, which applies single ice crystal habit and multidirectional polarized radiance to retrieve ice cloud particle shapes, is assessed. Our sensitivity studies show that the retrieved particle shapes &om this algorithm can be considered good approximations to those in actual clouds in calculation of the phase matrix elements. Although a fractal poly-crystal shape or an inhomogeneous hexagonal column may also produce this type of phase functions, more representative single-scattering properties from combinations of natural particle shapes and size distributions may still be necessary in accurate retrieval of other cloud properties such as optical thickness and particle size.

Sun, W.↗

Isotopic Constraints on the Genesis of Carbonates in Martian Meteorite ALH 84001

Oxygen isotopic analyses in approximately 20 micrometer spots in a chemically diverse suite of carbonates from ALH 84001 show highly variable delta(exp 18)O values from +5.4 to +25.3%. The isotopic data are correlated with the major element composition of the carbonate. The earliest forming (Ca-rich) carbonates have the lowest delta(exp 18)O values and the late-forming Mg-rich carbonates have the highest delta(exp 18)O values. Two models that can explain the isotopic variation were investigated. The carbonates could have formed in a water-rich environment at relatively low, but highly variable temperatures. In this open-system case the lower limit to the temperature variation is approximately 125 C, with fluctuations of over 250 C possible within the constraints of the model, depending on fluid composition. Alternatively the data can be explained by a closed-system model in which carbonates precipitated from a limited amount of a CO2-rich fluid. This scenario can reproduce the range of isotopic values observed, even at relatively high temperatures (greater than 500 C). Thus, the oxygen isotopic compositions do not provide unequivocal evidence for formation of the carbonates at low temperature. Neither of these scenarios is consistent with a biological origin of the carbonates and their associated features. Olivine from ALH 84001 occurs as clusters within orthopyroxene adjacent to fractures containing disrupted carbonate globules and feldspathic shock glass. The inclusions are irregular in shape and range in size from approximately 40 micrometers to submicrometer. The olivine exhibits a limited range of chemical composition from Fo(sub 65) to Fo(sub 66). We measured delta(exp 18)O values of the olivine to be +5.1 +/- 1.4%, indistinguishable within uncertainty from the host orthopyroxene. The data suggest that the olivine formed at high temperature (greater than 800 C), and is probably unrelated to carbonate formation. Instead the olivine probably formed by metamorphic reactions involving either dehydration of hydrous silicates or reduction of opx + spinel. If the reactions took place after carbonate formation, they could have caused devolatilization of siderite, producing the magnetite in the ALH carbonate globules. This scenario is also inconsistent with a biological origin for the features in ALH 84001.

Leshin, Laurie A.↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Model

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology from the small MVD conditions. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Meta-aerogels: Auxetic Shape-memory Polyurethane Aerogels

Shape-memory poly(isocyanurate-urethane) (PIR-PUR) aerogels are low-density monolithic nanoporous solids that remember and return to their permanent shape through a heating actuation step. Herein, through structural design at the macro scale, the shape-memory response is augmented with an auxetic effect manifested by a negative Poisson’s ratio of approximately -0.8 at 15% compressive strain. Thus, auxetic shape-memory PIR-PUR monoliths experience volume contraction upon compression at a temperature above the glass transition temperature of the base polymer (Tg ~ 30 °C), and they can be stowed indefinitely in that temporary shape by cooling below Tg. By heating back above Tg, the compressed/shrunk form expands back to their original shape/size. This technology is relevant to a broad range of industries spanning the commercial, aeronautical, and aerospace sectors. The materials are referred to as meta-aerogels and their potential applications include minimally invasive medical devices, soft robotics, and in situations where volume is at a premium, as for example for storage of deployable space structures and planetary habitats during transport to the point of service.

Sadeq Malakooti↗

Approximation theory for boundary layer suction through individual slits

The basic concepts of influencing boundary layers are summarized, especially the prevention of flow detachment and the reduction of frictional resistance. A mathematical analysis of suction through a slit is presented with two parameters, for thickness and for shape of the boundary layer, being introduced to specify the flow's velocity profile behind the slit. An approximation of the shape parameter produces a useful formula, which can be used to determine the most favorable position of the slit. An aerodynamic example is given.

Walz, A.↗

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗

Crack shapes and stress intensity factors for edge-cracked specimens.

A simple stress intensity factor expression is given for a deep edge crack in a plate in tension. The shapes of cracks opened by tension or bending are approximated by conic sections, and the conic section coefficients are related to plate geometry by very simple empirical equations. The magnitude of the crack displacement is a function of applied load, plate geometry, and the elastic constants of the plate material. The shape of a loaded crack in a semiinfinite plate is, approximately, a portion of an ellipse whose semimajor axis is about three times the crack length. As the crack length (relative to the plate width) increases, the crack shape becomes parabolic, then hyperbolic, the acuity of the hyperbola increasing with the relative crack length.

Orange, T. W.↗

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗

X-ray observations of a flare in NGC4151 from OSO-8

The 2-60 keV flux from NGC4151 changed by a factor of two on a timescale of 1.5 days. No fluctuations were detected in excess of a factor of three on timescales less than four hours. During a total observation of approximately 11 days there were no statistically significant changes in spectral shape. The spectrum was fitted by a power law with photon index alpha approximately 1.42 + or - 0.06 and column density N sub H approximately 7.5 + or - 0.5 x 10 to the 22d power at/cu cm. A 2 sigma residual to this fit implies fluorescent Fe line emission with E. W. approximately 240 eV. Both synchrotron self-Compton and thermal Compton models are consistent with the X-ray data.

Mushotzky, R. F.↗

Constraining Cometary Crystal Shapes from IR Spectral Features

A major challenge in deriving the silicate mineralogy of comets is ascertaining how the anisotropic nature of forsterite crystals affects the spectral features' wavelength, relative intensity, and asymmetry. Forsterite features are identified in cometary comae near 10, 11.05-11.2, 16, 19, 23.5, 27.5 and 33 microns [1-10], so accurate models for forsterite's absorption efficiency (Qabs) are a primary requirement to compute IR spectral energy distributions (SEDs, lambdaF lambda vs. lambda) and constrain the silicate mineralogy of comets. Forsterite is an anisotropic crystal, with three crystallographic axes with distinct indices of refraction for the a-, b-, and c-axis. The shape of a forsterite crystal significantly affects its spectral features [13-16]. We need models that account for crystal shape. The IR absorption efficiencies of forsterite are computed using the discrete dipole approximation (DDA) code DDSCAT [11,12]. Starting from a fiducial crystal shape of a cube, we systematically elongate/reduce one of the crystallographic axes. Also, we elongate/reduce one axis while the lengths of the other two axes are slightly asymmetric (0.8:1.2). The most significant grain shape characteristic that affects the crystalline spectral features is the relative lengths of the crystallographic axes. The second significant grain shape characteristic is breaking the symmetry of all three axes [17]. Synthetic spectral energy distributions using seven crystal shape classes [17] are fit to the observed SED of comet C/1995 O1 (Hale-Bopp). The Hale-Bopp crystalline residual better matches equant, b-platelets, c-platelets, and b-columns spectral shape classes, while a-platelets, a-columns and c-columns worsen the spectral fits. Forsterite condensation and partial evaporation experiments demonstrate that environmental temperature and grain shape are connected [18-20]. Thus, grain shape is a potential probe for protoplanetary disk temperatures where the cometary crystalline forsterite formed. The forsterite crystal shapes (equant, b-platelets, c-platelets, b-columns - excluding a- and c-columns) derived from our modeling [17] of comet Hale- Bopp, compared to laboratory synthesis experiments [18], suggests that these crystals are high temperature condensates. By observing and modeling the crystalline features in comet ISON, we may constrain forsterite crystal shape(s) and link to their formation temperature(s) and environment(s).

anisotropic crystal↗

Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy - Effect of crown shape and mutual shadowing

Many natural vegetation covers can be regarded as assemblages of plant crowns, on a background, which interact with light as discrete objects. On this basis, geometric optics furnish an approach to the modeling of vegetation canopies that captures the most important features of such growths' bidirectional measurements. Attention is presently given to models which approximate these phenomena and relate the size, shape, and count density of plant crowns to viewing and illumination positions and to crown-background reflectance contrasts.

Li, Xiaowen↗

Calculation of linearized supersonic flow over slender cones of arbitrary cross section

Supersonic linearized conical-flow theory is used to determine the flow over slender pointed cones having horizontal and vertical planes of symmetry. The geometry of the cone cross sections and surface velocities are expanded in Fourier series. The symmetry condition permits the uncoupling of lifting and nonlifting solutions. The present method reduces to Ward's theory for flow over a cone of elliptic cross section. Results are also presented for other shapes. Results by this method diverge for cross-sectional shapes where the maximum thickness is large compared with the minimum thickness. However, even for these slender-body shapes, lower order solutions are good approximations to the complete solution.

Mascitti, V. R.↗

Low-cost, air mass 2 solar simulator

A low cost air mass 2 solar simulator for testing flat-plate solar collectors was developed and operated at Lewis. Total cost was less than $10,800/sq m ($1000/sq ft). It consists of an array of 143 tungsten-halogen 300-watt lamps, each having integral dichroic coated reflectors. A second array of 143 respective hexagonal-shaped plastic Fresnel lenses are located approximately 1 focal length from the lamps. The simulator will produce a uniform collimated beam covering an area of 1.2 by 1.2m(4 by 4 ft). Design features, construction details, and component costs are given, as well as measured results on a 12-lamp prototype.

Yass, K.↗

Wind modeling and lateral control for automatic landing

For the purposes of aircraft control system design and analysis, the wind can be characterized by a mean component which varies with height and by turbulent components which are described by the von Karman correlation model. The aircraft aero-dynamic forces and moments depend linearly on uniform and gradient gust components obtained by averaging over the aircraft's length and span. The correlations of the averaged components are then approximated by the outputs of linear shaping filters forced by white noise. The resulting model of the crosswind shear and turbulence effects is used in the design of a lateral control system for the automatic landing of a DC-8 aircraft.

Holley, W. E.↗