Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Point Method”

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 343 records · Page 19

Minimizing Segregation During the Controlled Directional Solidification of Dendritic Alloys Publication: Metallurgical and Materials Transactions

Gravity-driven thermosolutal convection that arises during controlled directional solidification (DS) of dendritic alloys promotes detrimental macro-segregation (e.g. freckles and steepling) in products such as turbine blades. Considerable time and effort has been spent to experimentally and theoretically investigate this phenomena; although our knowledge has advanced to the point where convection can be modeled and accurately compared to experimental results, little has been done to minimize its onset and deleterious effects. The experimental work demonstrates that segregation can be. minimized and microstructural uniformity promoted when a slow axial rotation is applied to the sample crucible during controlled directional solidification processing. Numerical modeling utilizing continuation and bifurcation methods have been employed to develop accurate physical and mathematical models with the intent of identifying and optimizing processing parameters.

Grugel, R. N.↗

Capacitive skin characterization

NASA is currently involved in research that utilizes a capacitive sensor that is used for proximity detection of objects. This sensor is sensitive to conductive and dielectric materials including metal objects and humans. The range of the sensor has been found to be about twelve inches. It is the goal of this research project to further characterize the sensor so that it can be tailored to specific requirements. The characterization of the sensor should be with respect to shield size, sensor size, object size, and object distance. The method of finite elements to calculate the capacitance of the sensor while varying different parameters was used. Each of the parameters was varied in turn, often by selecting data points from different runs. The plotted results are shown and an apparent functionality developed for each.

Mcconnell, Robert↗

An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses.

Rose, Cheryl A.↗

Compact Electric Bi-Directional Pulley Actuator

While significant work has been dedicated to the study of potential uses of in situ resources once mined, the methods to mine the raw materials has not been fully explored. This paper presents the first attempt at producing a drop in replacement for hydraulic cylinders on heavy construction equipment with an electromechanical actuator utilizing pulleys for the conversion of a single motor’s rotational motion to bidirectional linear motion. By utilizing a twin block and tackle arrangement and paired spools, the system can generate bidirectional motion nearly identical to a hydraulic cylinder. More importantly, the use of rope builds in shock resistance to the system which is often a weak point for screw style actuators. Significant work remains to refine the concept into flight hardware, but the concept shows promise. The produced actuator can be applied or scaled for near direct replacement of hydraulic cylinders on terrestrial mining and construction equipment for faster development of off Earth use.

Kevin Peterson↗

The Dynamics of Solid Propellant Rocket Motors

From the point of view of dynamics, a solid propellant rocket motor is a unique structure in that it is composed of a substantial mass of propellant material case-bonded to a relatively massless, thin-walled cylinder. The mechanical properties of the propellant are such that it contributes little to the stiffness of the composite structure but it does contribute to the dynamic characteristics of the structure because of its mass. Furthermore, due to the viscoelastic character of the propellant, it can be expected to provide considerable damping to the system. At this point in the development of the art of design of solid propellant rocket motors there are no clear-cut or well-founded methods to evaluate quantitatively the contributions of the propellant to the dynamic responses of the composite structure. It is clear that unless such methods are devised, it will be difficult to arrive at accurate missile designs. In this paper a survey of the dynamic problems of solid propellant rocket motors is presented starting from the simplest model thereof and proceeding, step-by-step, to the consideration of more sophisticated and realistic models. The consideration is restricted to infinitesimal deformations of propellant grains with linear mechanical properties. Substantial progress has been achieved towards the solution of many important dynamical problems. We shall attempt to summarize the pertinent developments and indicate along which lines we feel future study should proceed. A few illustrative solutions are included in areas wherein progress has been substantial.

DYNAMIC RESPONSE↗

Ultrasonic measurement of elastic moduli in slender specimens using extensional and torsional wave pulses.

The pulse-echo method using extensional and torsional waves is described relative to determining elastic properties and other characteristics of specimens which are slender acoustically. Several applications are reviewed, including: glass fibers of diameter 16 to 42 microns; glass rods of diameter about 1 mm, tested to 700 C; iron-base alloys of diameter 3.2 mm, tested from nearly -200 up to +200 C; titanium specimens of rectangular cross section, 1.5 by 3 mm; and tantalum-base alloys of diameter 6.4 mm, tested up to 1400 C. Limitations are pointed out regarding dispersive materials and geometry.

Lynnworth, L. C.↗

Moment method with isoparametric elements for three-dimensional anisotropic scatterers

A novel method for computing the frequency-domain electromagnetic fields scattered from, and penetrating into, arbitrarily shaped, three-dimensional, lossy, inhomogeneous anisotropic scatters is presented. The method is based on a general volume integrodifferential formulation of the scattering problem and consists of the numerical solution of the coupled integral equations by the moment method and point matching. The numerical model of the scatterer is obtained by parametric volume elements, and the basis functions used to represent the field within each element are the same used in the finite-element method. Element integration problems due to the singular kernel of the integral equations are treated in some detail. Numerical results for both the isotropic and the anisotropic spherical scatterer are presented, including comparisons with results obtained by different numerical methods for the isotropic cases considered. The capability of the numerical code presented to deal with cases where the material parameters of the scatterer are given by singular matrices is discussed for two particular examples.

Graglia, Roberto D.↗

Preliminary experimental results for a cryogenic brush seal configuration

Preliminary fluid nitrogen flow data are reported for a five-brush, ceramic-coated-rub-runner brush seal system, where the brushes and the rub runner were placed at each end of a centrally pressurized multifunction tester ('back-to-back' set of brushes) and tested at rotor speeds of 0, 10, 18, and 22.5 krpm. After testing, both the brushes and the ceramic-coated rub runner appeared pristine. The coating withstood both the thermomechanical and dynamic loadings with minor wear track scarring. The bristle tips showed some indication of material shearing (smearing) wear. The Ergun porous flow equation was applied to the brush seal data. The Ergun relation, which required heuristic information to characterize the coefficients, fit the gaseous data but was in poor agreement with the fluid results. The brush seal exit conditions were two phase. Two-phase, choked-flow design charts were applied but required one data point at each rotor speed to define the (C(sub f)A x Constant) flow and area coefficients. Reasonable agreement between prediction and data was found, as expected, but such methods are not to be construed as two-phase-flow brush seal analyses.

Carlile, J. A.↗

Evaluation of Double-Vacuum-Bag Process For Composite Fabrication

A non-autoclave vacuum bag process using atmospheric pressure alone that eliminates the need for external pressure normally supplied by an autoclave or a press is an attractive method for composite fabrication. This type of process does not require large capital expenditures for tooling and processing equipment. In the molding cycle (temperature/pressure profile) for a given composite system, the vacuum application point has to be carefully selected to achieve the final consolidated laminate net shape and resin content without excessive resin squeeze-out. The traditional single-vacuum- bag (SVB) process is best suited for molding epoxy matrix based composites because of their superior flow and the absence of reaction by-products or other volatiles. Other classes of materials, such as polyimides and phenolics, generate water during cure. In addition, these materials are commonly synthesized as oligomers using solvents to facilitate processability. Volatiles (solvents and reaction byproducts) management therefore becomes a critical issue. SVB molding, without additional pressure, normally fails to yield void-free quality composites for these classes of resin systems. A double-vacuum- bag (DVB) process for volatile management was envisioned, designed and built at the NASA Langley Research Center. This experimental DVB process affords superior volatiles management compared to the traditional SVB process. Void-free composites are consistently fabricated as measured by C-scan and optical photomicroscopy for high performance polyimide and phenolic resins.

Hou, T. H.↗

Gigaflop (billion floating point operations per second) performance for computational electromagnetics

Accurate and rapid evaluation of radar signature for alternative aircraft/store configurations would be of substantial benefit in the evolution of integrated designs that meet radar cross-section (RCS) requirements across the threat spectrum. Finite-volume time domain methods offer the possibility of modeling the whole aircraft, including penetrable regions and stores, at longer wavelengths on today's gigaflop supercomputers and at typical airborne radar wavelengths on the teraflop computers of tomorrow. A structured-grid finite-volume time domain computational fluid dynamics (CFD)-based RCS code has been developed at the Rockwell Science Center, and this code incorporates modeling techniques for general radar absorbing materials and structures. Using this work as a base, the goal of the CFD-based CEM effort is to define, implement and evaluate various code development issues suitable for rapid prototype signature prediction.

Shankar, V.↗

Study to define points of entry for potential contaminants in limestone aquifers

Visual examinations of both prints and transparencies from ERTS 1 and U-2 aircraft imagery provided a method for discovering possible points of entry of potential contaminants into the limestone aquifer in Madison County, Alabama. Knowledge of the locations at which contaminants could enter the aquifer is an important consideration in water quality management, particularly for regions that depend, at least partially, on ground water for their water supply. ERTS 1 imagery recorded on December 28, 1972 in the Multispectral Scanner-5 (MSS-5) and MSS-7 bands, and a false-color composite of the MSS-4 (green), MSS-5 (red), and MSS-7 (near infrared) bands were the principal materials used, along with thermography recorded by an RS-7 infrared scanner onboard a U-2 aircraft. The results of the study are discussed in detail, providing information on prominent lineations and major fracture trends which are related to aquifer contamination. Maps depicting the observations are also presented.

Doyle, F. L.↗

Control systems using modal domain optical fiber sensors for smart structure applications

Recently, a new class of sensors has emerged for structural control which respond to environmental changes over a significant gauge length; these sensors are called distributed-effect sensors. These sensors can be fabricated with spatially varying sensitivity to the distributed measurand, and can be configured to measure a variety of structural parameters which can not be measured directly using point sensors. Examples of distributed-effect sensors include piezoelectric film, holographic sensors, and modal domain optical fiber sensors. Optical fiber sensors are particularly attractive for smart structure applications because they are flexible, have low mass, and can easily be embedded directly into materials. In this paper we describe the implementation of weighted modal domain optical fiber sensors. The mathematical model of the modal domain optical fiber sensor model is described and used to derive an expression for the sensor sensitivity. The effects of parameter variations on the sensor sensitivity are demonstrated to illustrate methods of spatially varying the sensor sensitivity.

Lindner, Douglas K.↗

Phase-Shifting Liquid Crystal Point-Diffraction Interferometry

Microgravity fluid physics experiments frequently measure concentration and temperature. Interferometers such as the Twyman Green illustrated have performed full-field measurement of these quantities. As with most such devices, this interferometer uses a reference path that is not common with the path through the test section. Recombination of the test and reference wavefronts produces interference fringes. Unfortunately, in order to obtain stable fringes, the alignment of both the test and reference paths must be maintained to within a fraction of the wavelength of the light being used for the measurement. Otherwise, the fringes will shift and may disappear. Because these interferometers are extremely sensitive to bumping, jarring and transmitted vibration, they are typically mounted on optical isolation tables. Schlieren deflectometers or the more recent Shack-Hartmann wavefront sensors also measure concentration and temperature in laboratory fluid flows. Ray optics describe the operation of both devices. In a schlieren system, an expanded, collimated beam passes through a test section where refractive index gradients deflect rays. A lens focuses the beam to a filter placed in the rear focal plane of the decollimating lens. In a quantitative color schlieren system, gradients in the index of refraction appear as colors in the field of view due to the action of the color filter. Since sensitivity is a function of the focal length of the decollimating lens, these systems are rather long and filter fabrication and calibration is rather difficult. A Shack-Hartmann wavefront sensor is an array of small lenslets. Typical diameters are on the order of a few hundred microns. Since these lenslets divide the test section into resolution elements, the spatial resolution can be no smaller than an individual lenslet. Such a device was recently used to perform high-speed tomography of heated air exiting a 1.27 cm diameter nozzle. While these wavefront sensors are very compact, the limited spatial resolution and the methods required for data reduction suggest that a more useful instrument needs to be developed. The category of interferometers known as common path interferometers can eliminate much of the vibration sensitivity associated with traditional interferometry as described above. In these devices, division of the amplitude of the wavefront following the test section produces the reference beam. Examples of these instruments include shearing and point diffraction interferometers. In the latter case, shown schematically, a lens focuses light passing through the test section onto a small diffracting object. Such objects are typically either a circle of material on a high quality glass plate or a small sphere in a glass cell. The size of the focused spot is several times larger than the object so that the light not intercepted by the diffracting object forms the test beam while the diffracted light generates a spherical reference beam. While this configuration is mechanically stable, phase shifting one beam with respect to the other is difficult due to the common path. Phase shifting enables extremely accurate measurements of the phase of the interferogram using only gray scale intensity measurements and is the de facto standard of industry. Mercer and Creath 2 demonstrated phase shifting in a point diffraction interferometer using a spherical spacer in a liquid crystal cell as the diffracting object. By changing the voltage across the cell, they were able to shift the phase of the undiffracted beam relative to the reference beam generated by diffraction from the sphere. While they applied this technology to fluid measurements, the device shifted phase so slowly that it was not useful for studying transient phenomena. We have identified several technical problems that precluded operation of the device at video frame rates and intend to solve them to produce a phase-shifting liquid crystal point-diffraction interferometer operating at video frame rates. The first task is to produce high contrast fringes. Since the diffracted beam is much weaker than the transmitted beam, interferograms have poor contrast unless a dye is added to the liquid crystal to reduce the intensity of the undiffracted light. Dyes previously used were not rigorously characterized and suffered from hysteresis in both the initial alignment state of the device and the electro-optic switching characteristics. Hence, our initial effort will identify and characterize dyes that do not suffer from these difficulties and are readily soluble in the liquid crystal host. Since the ultimate goal of this research is to produce interferometers capable of phase shifting at video frame rates, we will quantify the difference in switching times between ferroelectric and nematic liquid crystals. While we have more experience with nematic crystals, they typically switch more slowly than ferroelectric cells. As part of that effort, we will investigate the difference in the modulation of the interferograms as a function of the type of liquid crystal in the cell. Because the temporal switching response of a liquid crystal cell is directly related its thickness, we intend to explore techniques required to produce cells that are as thin as possible. However, the cells must still produce a total phase shift of two pi radians.

Griffin, DeVon W.↗

Infrared Observations of the Orion Capsule During EFT-1 Hypersonic Reentry

High-resolution infrared observations of the Orion capsule during its atmospheric reentry on December 5, 2015 were made from a US Navy NP-3D. This aircraft, equipped with a long-range optical sensor system, tracked the capsule from Mach 10 to 7 from a distance of approximately 60 nmi. Global surface temperatures of the capsule's thermal heatshield were derived from near infrared intensity measurements. The global surface temperature measurements complemented onboard instrumentation and were invaluable to the interpretation of the in-depth thermocouple measurements which rely on inverse heat transfer methods and material response codes to infer the desired surface temperature from the sub-surface measurements. The full paper will address the motivations behind the NASA Engineering Safety Center sponsored observation and highlight premission planning processes with an emphasis on aircraft placement, optimal instrument configuration and sensor calibrations. Critical aspects of mission operations coordinated from the NASA Johnson Spaceflight Center and integration with the JSC Flight Test Management Office will be discussed. A summary of the imagery that was obtained and processed to global surface temperature will be presented. At the capsule's point of closest approach relative to the imaging system, the spatial resolution was estimated to be approximately 15-inches per pixel and was sufficient to identify localized temperature increases associated with compression pad support hardware on the heatshield. The full paper will discuss the synergy of the quantitative imagery derived temperature maps with in-situ thermocouple measurements. Comparison of limited onboard surface thermocouple data to the image derived surface temperature will be presented. The two complimentary measurements serve as an example of the effective leveraging of resources to advance the understanding of high Mach number environments associated with an ablated heatshield and provide unique data for the validation of design tools and numerical flight simulation techniques. Collaborative opportunities and technology investments in support of planned observations of NASA's next Orion flight test in 2018 will be explored in the full manuscript.

Horvath, Thomas J.↗

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)↗

SKITTER foot design

A mechanical design team was formed to design a foot for the lunar utility vehicle SKITTER. The primary design was constrained to be a ski pole design compatible with the existing femur-tibia design legs. The lunar environment had several important effects on the foot design. Three materials were investigated for the SKITTER foot: aluminum alloys, cold worked stainless steel alloys, and titanium alloys. Thin film coatings were investigated as a method of wear reduction for the foot. The performance of the foot is dependent on the action of the legs. The range of motion for the legs was determined to be vertical to 15 degrees above horizontal. An impact analysis was performed for the foot movement, but the results were determined to be inconclusive due to unknown soil parameters. The initial foot design configuration consisted of an annulus attached to the pointed pole. The annulus was designed to prevent excess sinkage. Later designs call for a conical shaped foot with a disk at the point of the tibia attachment. The conical design was analyzed for strength and deflection by two different approaches. A deformable body analysis was performed for the foot under crane load in crane position, and also under actuator load in the vertical position. In both cases, the deflection of the foot was insignificant and the stresses well below the strength of the titanium alloy.

Choi, Gene↗

A Unique Computational Algorithm to Simulate Probabilistic Multi-Factor Interaction Model Complex Material Point Behavior

The Multi-Factor Interaction Model (MFIM) is used to evaluate the divot weight (foam weight ejected) from the launch external tanks. The multi-factor has sufficient degrees of freedom to evaluate a large number of factors that may contribute to the divot ejection. It also accommodates all interactions by its product form. Each factor has an exponent that satisfies only two points--the initial and final points. The exponent describes a monotonic path from the initial condition to the final. The exponent values are selected so that the described path makes sense in the absence of experimental data. In the present investigation, the data used was obtained by testing simulated specimens in launching conditions. Results show that the MFIM is an effective method of describing the divot weight ejected under the conditions investigated.

Chamis, Christos C.↗

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗