Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing methods”

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,117 records · Page 62

Flux-focusing eddy current probe and method for flaw detection

A flux-focusing electromagnetic sensor which uses a ferromagnetic flux-focusing lens simplifies inspections and increases detectability of fatigue cracks and material loss in high conductivity material is presented. The unique feature of the device is the ferrous shield isolating a high-turn pick-up coil from an excitation coil. The use of the magnetic shield is shown to produce a null voltage output across the receiving coil in the presence of an unflawed sample. A redistribution of the current flow in the sample caused by the presence of flaws, however, eliminates the shielding condition and a large output voltage is produced, yielding a clear unambiguous flaw signal. The maximum sensor output is obtained when positioned symmetrically above the crack. Hence, by obtaining the position of the maximum sensor output, it is possible to track the fault and locate the area surrounding its tip. The accuracy of tip location is enhanced by two unique features of the sensor; a very high signal-to-noise ratio of the probe's output which results in an extremely smooth signal peak across the fault, and a rapidly decaying sensor output outside a small area surrounding the crack tip which enables the region for searching to be clearly defined. Under low frequency operation, material thinning due to corrosion damage causes an incomplete shielding of the pick-up coil. The low frequency output voltage of the probe is therefore a direct indicator of the thickness of the test sample.

Simpson, John W.↗

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials↗

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials↗

New multigrid approach for three-dimensional unstructured, adaptive grids

A new multigrid method with adaptive unstructured grids is presented. The three-dimensional Euler equations are solved on tetrahedral grids that are adaptively refined or coarsened locally. The multigrid method is employed to propagate the fine grid corrections more rapidly by redistributing the changes-in-time of the solution from the fine grid to the coarser grids to accelerate convergence. A new approach is employed that uses the parent cells of the fine grid cells in an adapted mesh to generate successively coaser levels of multigrid. This obviates the need for the generation of a sequence of independent, nonoverlapping grids as well as the relatively complicated operations that need to be performed to interpolate the solution and the residuals between the independent grids. The solver is an explicit, vertex-based, finite volume scheme that employs edge-based data structures and operations. Spatial discretization is of central-differencing type combined with a special upwind-like smoothing operators. Application cases include adaptive solutions obtained with multigrid acceleration for supersonic and subsonic flow over a bump in a channel, as well as transonic flow around the ONERA M6 wing. Two levels of multigrid resulted in reduction in the number of iterations by a factor of 5.

Parthasarathy, Vijayan↗

Theoretical effect of modifications to the upper surface of two NACA airfoils using smooth polynomial additional thickness distributions which emphasize leading edge profile and which vary quadratically at the trailing edge

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of the NACA 64-206 and 64 sub 1 - 212 airfoils. The additional thickness distribution had the form of a continuous mathematical function which disappears at both the leading edge and the trailing edge. The function behaves as a polynomial of order epsilon sub 1 at the leading edge, and a polynomial of order epsilon sub 2 at the trailing edge. Epsilon sub 2 is a constant and epsilon sub 1 is varied over a range of practical interest. The magnitude of the additional thickness, y, is a second input parameter, and the effect of varying epsilon sub 1 and y on the aerodynamic performance of the airfoil was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic airfoils, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.

Merz, A. W.↗

A Globally Optimal Particle Tracking Technique for Stereo Imaging Velocimetry Experiments

An important phase of any Stereo Imaging Velocimetry experiment is particle tracking. Particle tracking seeks to identify and characterize the motion of individual particles entrained in a fluid or air experiment. We analyze a cylindrical chamber filled with water and seeded with density-matched particles. In every four-frame sequence, we identify a particle track by assigning a unique track label for each camera image. The conventional approach to particle tracking is to use an exhaustive tree-search method utilizing greedy algorithms to reduce search times. However, these types of algorithms are not optimal due to a cascade effect of incorrect decisions upon adjacent tracks. We examine the use of a guided evolutionary neural net with simulated annealing to arrive at a globally optimal assignment of tracks. The net is guided both by the minimization of the search space through the use of prior limiting assumptions about valid tracks and by a strategy which seeks to avoid high-energy intermediate states which can trap the net in a local minimum. A stochastic search algorithm is used in place of back-propagation of error to further reduce the chance of being trapped in an energy well. Global optimization is achieved by minimizing an objective function, which includes both track smoothness and particle-image utilization parameters. In this paper we describe our model and present our experimental results. We compare our results with a nonoptimizing, predictive tracker and obtain an average increase in valid track yield of 27 percent

McDowell, Mark↗

Surface fitting three-dimensional bodies

The geometry of general three-dimensional bodies is generated from coordinates of points in several cross sections. Since these points may not be smooth, they are divided into segments and general conic sections are curve fit in a least-squares sense to each segment of a cross section. The conic sections are then blended in the longitudinal direction by fitting parametric cubic-spline curves through coordinate points which define the conic sections in the cross-sectional planes. Both the cross-sectional and longitudinal curves may be modified by specifying particular segments as straight lines and slopes at selected points. Slopes may be continuous or discontinuous and finite or infinite. After a satisfactory surface fit has been obtained, cards may be punched with the data necessary to form a geometry subroutine package for use in other computer programs. At any position on the body, coordinates, slopes and second partial derivatives are calculated. The method is applied to a blunted 70 deg delta wing, and it was found to generate the geometry very well.

Dejarnette, F. R.↗

Remarks on the theory of aerodynamic noise

The accuracy of approximations employed in the Lighthill theory in current aerodynamic noise research is critically evaluated. Based on the method of matched asymptotic expansions, the full Navier-Stokes equations are expanded for small Mach numbers. To the first order, the near-field is generally nonisentropic. The pressure field (pseudo-sound) is generated by the incompressible Reynolds stresses in the turbulent flow and the velocity, pressure perturbation, and their derivatives on the boundaries. In the far-field, the first-order pressure (acoustic) field satisfying a linear wave equation is obtained by matching with the pseudo-sound field. A uniformly valid solution for the pressure field in a stationary or uniformly moving medium is obtained. The solution shows that the generation of the first-order aerodynamic noise does not depend on the viscous, thermal, or entropy effects in the adiabatic flow nor on the shear stress on a smooth rigid boundary.

Pan, Y. S.↗

Mach 5 electroformed nickel nozzle refurbishment FNAS investigation of ultra-smooth surfaces

This task is in support of the Quiet Hypersonic Wind Tunnel effort currently in effect at NASA Langley Research Center, VA. A laminar flow wind tunnel nozzle has been previously fabricated by electroforming pure nickel over a machined and polished two piece stainless steel mandrel. The mandrel was removed leaving the replicate nozzle surface. The nozzle was then pressed into a heavy stainless steel jacket for mounting features and rigidity. The original nickel surface was a replication of the polished mandrel but had degraded due to oxidation. The inside surface requirements are very stringent in order to achieve laminar or quiet flow at the desired pressure and temperature for the specific design of Mach 5. The throat area of the axisymmetric device must have a surface finish with no defects greater than 16 microinches. This requires an rms average background of about four microinches or better for inspection purposes. A coating of nickel-phosphorous alloy was applied by catalytic deposition and then the inside of the nozzle was polished retaining dimensional and surface finish tolerances as specified per drawings supplied. Since the unit is not an optical component, conventional optical inspection methods for surface finish and figure are not readily achieved. Measurements were made using surface profilometry.

Rood, Robert↗

Flex Dynamics Avoidance Control of the NEA Scout Solar Sail Spacecraft's Reaction Control System

The Attitude Control System (ACS) is developed for a Near Earth Asteroid (NEA) Scout mission using a solar sail. The NEA-Scout spacecraft is a 6U cubesat with an 86 square-meter solar sail. NEA Scout will launch on Space Launch System (SLS) Exploration Mission 1 (EM-1), currently scheduled to launch in 2018. The spacecraft will rendezvous with a target asteroid after a two year journey, and will conduct science imagery. The solar sail spacecraft ACS consists of three major actuating subsystems: a Reaction Wheel (RW) control system, a Reaction Control System (RCS), and an Adjustable Mass Translator (AMT) system. The three subsystems allow for a wide range of spacecraft attitude control capabilities, needed for the different phases of the NEA-Scout mission. Because the sail is a flexible structure, care must be taken in designing a control system to avoid exciting the structural modes of the sail. This is especially true for the RCS, which uses pulse actuated, cold-gas jets to control the spacecraft's attitude. While the reaction wheels can be commanded smoothly, the RCS jets are simple on-off actuators. Long duration firing of the RCS jets - firings greater than one second - can be thought of as step inputs to the spacecraft's torque. On the other hand, short duration firings - pulses on the order of 0.1 seconds - can be thought of as impulses in the spacecraft's torque. These types of inputs will excite the structural modes of the spacecraft, causing the sail to oscillate. Sail oscillations are undesirable for many reasons. Mainly, these oscillations will feed into the spacecraft attitude sensors and pointing accuracy, and long term oscillations may be undesirable over the lifetime of the solar sail. In order to limit the sail oscillations, an RCS control scheme is being developed to minimize sail excitations. Specifically, an input shaping scheme similar to the method described in Reference 1 will be employed. A detailed description of the RCS control scheme will be provided with particular emphasis on flexible body excitation. The RCS performance will be provided to show that sail and boom excitation is minimized.

Heaton Andrew↗

Terminal area automatic navigation, guidance, and control research using the Microwave Landing System (MLS). Part 3: A comparison of waypoint guidance algorithms for RNAV/MLS transition

The results of an investigation carried out for the Langley Research Center Terminal Configured Vehicle Program are presented. The investigation generated and compared three path update algorithms designed to provide smooth transition for an aircraft guidance system from DME, VORTAC, and barometric navaids to the more precise MLS by modifying the desired 3-D flight path. The first, called the Zero Cross Track, eliminates the discontinuity in cross track and altitude error by designating the first valid MLS aircraft position as the desired first waypoint, while retaining all subsequent waypoints. The discontinuity in track angle is left unaltered. The second, called the Tangent Path also eliminates the discontinuity in cross track and altitude and choose a new desired heading to be tangent to the next oncoming circular arc turn. The third, called the Continued Track eliminates the discontinuity in cross track, altitude and track angle by accepting the current MLS position and track angle as the desired ones and recomputes the location of the next waypoint. A method is presented for providing a waypoint guidance path reconstruction which treats turns of less than, and greater than, 180 degrees in a uniform manner to construct the desired path.

Pines, S.↗

Recursive mass matrix factorization and inversion: An operator approach to open- and closed-chain multibody dynamics

This report advances a linear operator approach for analyzing the dynamics of systems of joint-connected rigid bodies.It is established that the mass matrix M for such a system can be factored as M=(I+H phi L)D(I+H phi L) sup T. This yields an immediate inversion M sup -1=(I-H psi L) sup T D sup -1 (I-H psi L), where H and phi are given by known link geometric parameters, and L, psi and D are obtained recursively by a spatial discrete-step Kalman filter and by the corresponding Riccati equation associated with this filter. The factors (I+H phi L) and (I-H psi L) are lower triangular matrices which are inverses of each other, and D is a diagonal matrix. This factorization and inversion of the mass matrix leads to recursive algortihms for forward dynamics based on spatially recursive filtering and smoothing. The primary motivation for advancing the operator approach is to provide a better means to formulate, analyze and understand spatial recursions in multibody dynamics. This is achieved because the linear operator notation allows manipulation of the equations of motion using a very high-level analytical framework (a spatial operator algebra) that is easy to understand and use. Detailed lower-level recursive algorithms can readily be obtained for inspection from the expressions involving spatial operators. The report consists of two main sections. In Part 1, the problem of serial chain manipulators is analyzed and solved. Extensions to a closed-chain system formed by multiple manipulators moving a common task object are contained in Part 2. To retain ease of exposition in the report, only these two types of multibody systems are considered. However, the same methods can be easily applied to arbitrary multibody systems formed by a collection of joint-connected regid bodies.

Rodriguez, G.↗

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis↗

Dynamic Forms: Application to Aircraft Guidance - Part 2

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of waypoints through which the aircraft trajectory must pass. The waypoints typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory that satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multidimensional, multiaxis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions are required to be smooth. The proposed guidance algorithm is based on the inversion of the pure feedback approximation, followed by correction for the effects of zero dynamics. The paper describes the structure and major modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

Durability and life prediction modeling in polyimide composites

Sudden appearance of cracks on a macroscopically smooth surface of brittle materials due to cooling or drying shrinkage is a phenomenon related to many engineering problems. Although conventional strength theories can be used to predict the necessary condition for crack appearance, they are unable to predict crack spacing and depth. On the other hand, fracture mechanics theory can only study the behavior of existing cracks. The theory of crack initiation can be summarized into three conditions, which is a combination of a strength criterion and laws of energy conservation, the average crack spacing and depth can thus be determined. The problem of crack initiation from the surface of an elastic half plane is solved and compares quite well with available experimental evidence. The theory of crack initiation is also applied to concrete pavements. The influence of cracking is modeled by the additional compliance according to Okamura's method. The theoretical prediction by this structural mechanics type of model correlates very well with the field observation. The model may serve as a theoretical foundation for future pavement joint design. The initiation of interactive cracks of quasi-brittle material is studied based on a theory of cohesive crack model. These cracks may grow simultaneously, or some of them may close during certain stages. The concept of crack unloading of cohesive crack model is proposed. The critical behavior (crack bifurcation, maximum loads) of the cohesive crack model are characterized by rate equations. The post-critical behavior of crack initiation is also studied.

Binienda, Wieslaw K.↗

Near-Surface Flow Fields Deduced Using Correlation Tracking and Time-Distance Analysis

Near-photospheric flow fields on the Sun are deduced using two independent methods applied to the same time series of velocity images observed by SOI-MDI on SOHO. Differences in travel times between f modes entering and leaving each pixel measured using time-distance helioseismology are used to determine sites of supergranular outflows. Alternatively, correlation tracking analysis of mesogranular scales of motion applied to the same time series is used to deduce the near-surface flow field. These two approaches provide the means to assess the patterns and evolution of horizontal flows on supergranular scales even near disk center, which is not feasible with direct line-of-sight Doppler measurements. We find that the locations of the supergranular outflows seen in flow fields generated from correlation tracking coincide well with the locations of the outflows determined from the time-distance analysis, with a mean correlation coefficient after smoothing of bar-r(sub s) = 0.840. Near-surface velocity field measurements can used to study the evolution of the supergranular network, as merging and splitting events are observed to occur in these images. The data consist of one 2048-minute time series of high-resolution (0.6" pixels) line-of-sight velocity images taken by MDI on 1997 January 16-18 at a cadence of one minute.

DeRosa, Marc↗

Considerations of Methods of Improving Helicopter Efficiency

Recent NASA helicopter research indicates that significant improvements in hovering efficiency, up to 7 percent, are available from the use of a special airfoil section formed by combining an NACA 632A015 thickness distribution with an NACA 230 mean line. This airfoil should be considered for flying-crane-type helicopters. Application of standard leading-edge roughness causes a large drop in efficiency; however, the cambered rotor is shown to retain its superiority over a rotor having a symmetrical airfoil when both rotors have leading-edge roughness. A simple analysis of available rotor static-thrust data indicates a greatly reduced effect of compressibility effects on the rotor profile-drag power than predicted from calculations. Preliminary results of an experimental study of helicopter parasite drag indicate the practicability of achieving an equivalent flat-plate parasite-drag area of less than 4 square feet for a rotor-head-pylon-fuselage configuration (landing gear retracted) in the 2,000-pound minimum-flying-weight class. The large drag penalty of a conventional skid-type landing (3.6 square feet) can be reduced by two-thirds by careful design. Clean, fair, and smooth fuselages that tend to have narrow, deep cross sections are shown to have advantages from the standpoint of drag and download. A ferry range of the order of 1,500 miles is indicated to be practicable for the small helicopter considered.

Dingeldein, Richard C.↗

Modified Geometric Truncation of the Scattering Phase Function

Phase function of light scattering on large atmospheric particles has very strong peak in forward direction constituting a challenge for accurate numerical calculations of radiance required in remote sensing problems. Scaling transformation replaces original phase function with a sum of the delta function and a new regular smooth phase function. Geometric truncation is one of the ways to construct such a smooth function. The replacement phase function coincides with the original one outside the forward cone and preserves the asymmetry parameter. It has discontinuity at the cone. Another simple functional form of the replacement phase function within the cone is suggested. It enables continuity and allows for a number of modifications. Three of them are considered in this study: preserving asymmetry parameter, providing continuity of the 1st derivative of the phase function, and preserving mean scattering angle. Yet another problem addressed in this study is objective selection of the width of the forward cone. That angle affects truncation fraction and values of the phase function within the cone. A heuristic approach providing unambiguous criterion of selection of the truncation angle is proposed. The approach has easy numerical implementation. Suggested modifications were tested on cloud phase function using discrete ordinates and Monte Carlo methods. It was shown that the modifications provide better accuracy of the radiance computation compare to the original geometric truncation with discrete ordinates while continuous derivative approach provides significant gain in computer time with Monte Carlo simulations.

Alexander Radkevich↗