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 271 records · Page 15

Improved boundary layer heat transfer calculations near a stagnation point

A thermal design of a solar receiver has been developed for the solutions of problems involving phase-change thermal energy storage and natural convection loss. Two dimensional axisymmetrical solidification and melting of materials contained between two concentric cylinders of finite length has been studied for thermal energy storage analysis. For calculation of free convection loss inside receiver cavity, two dimensional axisymmetrical, laminar, transient free convection including radiation effects has been studied using integral/finite difference method. Finite difference equations are derived for the above analysis subject to constant or variable material properties, initial conditions, and boundary conditions. The validity of the analyses has been substantiated by comparing results of the present general method with available analytic solutions or numerical results reported in the literature. Both explicit and implicit schemes are tested in phase change analysis with different number of nodes ranging from 4 to 18. The above numerical methods have been applied to the existing solar receiver analyzing computer code as additional subroutines. The results were computed for one of the proposed Brayton cycle solar receiver models running under the actual environmental conditions. Effect of thermal energy storage on the thermal behavior of the receiver has been estimated. Due to the thermal energy storage, about 65% reduction on working gas outlet temperature fluctuation has been obtained; however, maximum temperature of thermal energy storage containment has been increased about 18%. Also, effect of natural convection inside a receiver cavity on the receiver heat transfer has been analyzed. The finding indicated that thermal stratification occurs during the sun time resulting in higher receiver temperatures at the outlet section of the gas tube, and lower temperatures at the inlet section of the gas tube when compared with the results with no natural convection. Due to heat supply from the air during the shade time, minimum temperature has been increased, while maximum temperature has been reduced due to convection loss to air. Consequently, cyclic temperature fluctuation has been reduced 29% for working gas and 16% for thermal energy storage containment. On the other hand, despite the presence of the natural convection the time-averaged temperatures for receiver components were found to be similar for two cases with/without natural convection (maximum difference was 1.8%).

Ahn, Kyung Hwan↗

Investigation into Spectroscopic Techniques for Thermal Barrier Coating Spall Detection

Spectroscopic methods are proposed for detection of thermal barrier coating (TBC) spallation from engine hot zone components. These methods include absorption and emission of airborne marker species originally embedded in the TBC bond coat. In this study, candidate marker materials for this application were evaluated. Thermochemical analysis of candidate marker materials combined with additional constraints such as toxicity and uniqueness to engine environment, provided a short list of four potential species: platinum, copper oxide, zinc oxide. and indium. The melting point of indium was considered to be too low for serious consideration. The other three candidate marker materials, platinum, copper oxide, and zinc oxide were placed in a high temperature furnace and emission and absorption properties were measured over a temperature range from 800-1400 C and a spectral range from 250 to 18000 nm. Platinum did not provide the desired response, likely due to the low vapor Pressure of the metallic species and the low absorption of the oxide species. It was also found, however. that platinum caused a broadening of the carbon dioxide absorption at 4300 nm. The nature of this effect is not known. Absorption and emission caused by sodium and potassium impurities in the platinum were found in the platinum tests. Zinc oxide did not provide the desired response, again, most likely due to the low vapor pressure of the metallic species and the low absorption of the oxide species. Copper oxide generated two strongly temperature dependent absorption peaks at 324.8 and 327.4 nm. The melting point of copper oxide was determined to be too low for serious consideration as marker material.

deGroot, Wim↗

Real-Time Optical Monitoring and Simulations of Gas Phase Kinetics in InN Vapor Phase Epitaxy at High Pressure

Understanding the kinetics of nucleation and coalescence of heteroepitaxial thin films is a crucial step in controlling a chemical vapor deposition process, since it defines the perfection of the heteroepitaxial film both in terms of extended defect formation and chemical integrity of the interface. The initial nucleation process also defines the film quality during the later stages of film growth. The growth of emerging new materials heterostructures such as InN or In-rich Ga(x)In(1-x)N require deposition methods operating at higher vapor densities due to the high thermal decomposition pressure in these materials. High nitrogen pressure has been demonstrated to suppress thermal decomposition of InN, but has not been applied yet in chemical vapor deposition or etching experiments. Because of the difficulty with maintaining stochiometry at elevated temperature, current knowledge regarding thermodynamic data for InN, e.g., its melting point, temperature-dependent heat capacity, heat and entropy of formation are known with far less accuracy than for InP, InAs and InSb. Also, no information exists regarding the partial pressures of nitrogen and phosphorus along the liquidus surfaces of mixed-anion alloys of InN, of which the InN(x)P(1-x) system is the most interesting option. A miscibility gap is expected for InN(x)P(1-x) pseudobinary solidus compositions, but its extent is not established at this point by experimental studies under near equilibrium conditions. The extension of chemical vapor deposition to elevated pressure is also necessary for retaining stoichiometric single phase surface composition for materials that are characterized by large thermal decomposition pressures at optimum processing temperatures.

Dietz, Nikolaus↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

Development and Application of Computational/In Vitro Toxicological Methods for Chemical Hazard Risk Reduction of New Materials for Advanced Weapon Systems

The development of quantitative structure-activity relationship (QSAR) is essential for reducing the chemical hazards of new weapon systems. The current collaboration between HEST (toxicology research and testing), MLPJ (computational chemistry) and PRS (computational chemistry, new propellant synthesis) is focusing R&D efforts on basic research goals that will rapidly transition to useful products for propellant development. Computational methods are being investigated that will assist in forecasting cellular toxicological end-points. Models developed from these chemical structure-toxicity relationships are useful for the prediction of the toxicological endpoints of new related compounds. Research is focusing on the evaluation tools to be used for the discovery of such relationships and the development of models of the mechanisms of action. Combinations of computational chemistry techniques, in vitro toxicity methods, and statistical correlations, will be employed to develop and explore potential predictive relationships; results for series of molecular systems that demonstrate the viability of this approach are reported. A number of hydrazine salts have been synthesized for evaluation. Computational chemistry methods are being used to elucidate the mechanism of action of these salts. Toxicity endpoints such as viability (LDH) and changes in enzyme activity (glutahoione peroxidase and catalase) are being experimentally measured as indicators of cellular damage. Extrapolation from computational/in vitro studies to human toxicity, is the ultimate goal. The product of this program will be a predictive tool to assist in the development of new, less toxic propellants.

Frazier, John M.↗

Analysis of longitudinal thickness-shear stiffness of a monolayer filamentary composite

The longitudinal thickness-shear behavior of a rectangular cross-section prism containing circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Saint-Venant flexure problem modeled by a tip-loaded cantilever composite beam and solution is obtained by the boundary-point least-squares method. It is shown that the longitudinal thickness-shear modulus is much greater than the in-plane shear modulus for a monofilament composite and therefore an independent analysis is necessary to predict the longitudinal thickness-shear behavior of the composite. Numerical results are presented in graphical form for various values of fiber volume fraction and material parameters typical of modern advanced composites.

Chang, S.↗

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

Analysis of Radiation-natural Convection Interactions in 1-g and low-g Environments using the Discrete Exchange Factor Method

A new numerical method is presented for the analysis of combined natural convection and radiation heat transfer with applications in many engineering situations such as materials processing, combustion and fire research. Because of the recent interest in the low gravity environment of space, attention is devoted to both 1-g and low-g applications. The two-dimensional mathematical model is represented by a set of coupled nonlinear integro-partial differential equations. Radiative exchange is formulated using the Discrete Exchange Factor method (DEF). This method considers point to point exchange and provides accurate results over a wide range of radiation parameters. Numerical results show that radiation significantly influences the flow and heat transfer in both low-g and 1-g applications. In the low-g environment, convection is weak, and radiation can easily become the dominant heat transfer mode. It is also shown that volumetric heating by radiation gives rise to an intricate cell pattern in the top heated enclosure.

Kassemi, M.↗

Impact conditions required for formation of melt by jetting in silicates

It is demonstrated that the process of jetting which occurs when particles collide at oblique angles may produce melt at much lower velocities than are required for melt production in head-on collisions. The minimum velocities of impact required for jetting in aluminum, bronzitite, dunite, and quartz are calculated by the method of shock polars. The analysis, which depends on stated assumptions about attained pressures and the occurrence of jetting, uses shock-velocity particle-velocity equations of state with three shock regimes. The treatment indicates that jetting should arise in bronzitite, dunite, and quartz at relative velocities as low as 1-2 km/sec. At such velocities material which passes near the stagnation point in the jet-forming region is subjected to sufficiently high pressures so that it is probably melted.

Kieffer, S. W.↗

NESC GN&C TDT Workshop on 2D Image Motion Optical Transfer Functions, Pointing Performance Analysis, and Requirements

What You Will Learn: The focus is on payload imaging performance due to pointing motion. Some historical background on pointing performance analysis is given. The Optical Transfer Function (OTF) and Modulation Transfer Function (MTF) are defined. The imaging performance due to pointing motion is measured by image motion optical transfer functions (IM OTF). IM OTFs are defined for displacement, smear, and jitter motions, which are all rigorously defined. Deterministic and Statistical IM OTFs are briefly derived and graphically illustrated and compared. The IM OTFs are parameterized by pointing error metrics(PEM), which are means and covariances of displacement, smear, and jitter. Emphasis is on procedures and algorithms to evaluate the image motion optical transfer functions and pointing error metrics. Three procedures are covered, which depend on whether the pointing error data is from time-domain simulation, frequency-domain analysis, or stochastic modeling. A method to evaluate the relative contribution of disturbance sources and to identify the most significant contributors is presented. The presentation includes pertinent discussion of flexible structures and control-structure interaction. No single book can adequately cover this subject, so a book is not required for the course. A list of selected articles, reports, documents, and books is provided for reference and further study. Mathis kept to the minimum necessary to convey principles; lengthy derivations are left to the reference material. Graphics are used to illustrate concepts. As with any such learning endeavor, the knowledge gained will be retained and strengthened through actual practice.c©2019–2021 Mark E. Pittelkau— 5

NASA Engineering and Safety Center (NESC)↗

A quantitative method for photovoltaic encapsulation system optimization

It is pointed out that the design of encapsulation systems for flat plate photovoltaic modules requires the fulfillment of conflicting design requirements. An investigation was conducted with the objective to find an approach which will make it possible to determine a system with optimum characteristics. The results of the thermal, optical, structural, and electrical isolation analyses performed in the investigation indicate the major factors in the design of terrestrial photovoltaic modules. For defect-free materials, minimum encapsulation thicknesses are determined primarily by structural considerations. Cell temperature is not strongly affected by encapsulant thickness or thermal conductivity. The emissivity of module surfaces exerts a significant influence on cell temperature. Encapsulants should be elastomeric, and ribs are required on substrate modules. Aluminum is unsuitable as a substrate material. Antireflection coating is required on cell surfaces.

Garcia, A., III↗

In-Situ Measurement of Hall Thruster Erosion Using a Fiber Optic Regression Probe

One potential life-limiting mechanism in a Hall thruster is the erosion of the ceramic material comprising the discharge channel. This is especially true for missions that require long thrusting periods and can be problematic for lifetime qualification, especially when attempting to qualify a thruster by analysis rather than a test lasting the full duration of the mission. In addition to lifetime, several analytical and numerical models include electrode erosion as a mechanism contributing to enhanced transport properties. However, there is still a great deal of dispute over the importance of erosion to transport in Hall thrusters. The capability to perform an in-situ measurement of discharge channel erosion is useful in addressing both the lifetime and transport concerns. An in-situ measurement would allow for real-time data regarding the erosion rates at different operating points, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over a thruster s operating envelope would also be useful in the modeling of the detailed physics inside the discharge chamber. There are many different sensors and techniques that have been employed to quantify discharge channel erosion in Hall thrusters. Snapshots of the wear pattern can be obtained at regular shutdown intervals using laser profilometry. Many non-intrusive techniques of varying complexity and sensitivity have been employed to detect the time-varying presence of erosion products in the thruster plume. These include the use quartz crystal microbalances, emission spectroscopy, laser induced flourescence, and cavity ring-down spectroscopy. While these techniques can provide a very accurate picture of the level of eroded material in the thruster plume, it is more difficult to use them to determine the location from which the material was eroded. Furthermore, none of the methods cited provide a true in-situ measure of erosion at the channel surface while the thruster is in operation (i.e. none yield a continuous channel erosion measurement). A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor was tested using a linear Hall thruster geometry (see Fig. 1), which served as a means of producing plasma erosion of a ceramic discharge chamber. The mass flow rate, discharge voltage, and applied magnetic field strength could be varied, allowing for erosion measurements over a broad thruster operating envelope. Results are presented demonstrating the ability of the REAST sensor to capture not only the insulator erosion rates but also changes in these rates as a function of the discharge parameters.

Polzin, Kurt↗

Comparison tests and experimental compliance calibration of the proposed standard round compact plane strain fracture toughness specimen

Standard round specimen fracture test results compared satisfactorily with results from standard rectangular compact specimens machined from the same material. The location of the loading pin holes was found to provide adequate strength in the load bearing region for plane strain fracture toughness testing. Excellent agreement was found between the stress intensity coefficient values obtained from compliance measurements and the analytic solution proposed for inclusion in the standard test method. Load displacement measurements were made using long armed displacement gages and hollow loading cylinders. Gage points registered on the loading hole surfaces through small holes in the walls of the loading cylinders.

Fisher, D. M.↗

Freeze-drying for morphological control of high performance semi-interpenetrating polymer networks. III

The feasibility of using a freeze-drying (solvent removal by sublimation) approach for controlling the morphology of a high-performance semi-IPN is assessed. A high-performance thermoplastic polyimide and commercially available 4,4'-bismaleimide diphenylenemethane were dissolved in a solvent, 1,3,5-trioxane. The solvent was removed from the constituents by freeze-drying. For purposes of comparison, the constituents were dissolved in a high-boiling-point solvent, N,N-dimethylformamide. The solvent was removed from the solution by evaporation. The physical and mechanical properties and phase morphology of the neat resins and composites prepared by freeze-drying and traditional solution methods are presented and compared. It is concluded that the TG is higher and that the magnitude of minor constituent separation is less in the freeze-dry processed materials than for the processed solution.

Hsiung, H. J.↗

Kalman Filtering and RTS Smoothing for Arc-Jet Sample Edge Tracking

Accurate and temporally consistent measurements of test-article recession are required to characterize material response during arc-jet ablation experiments. However, image-based boundary measurements are often affected by segmentation noise, brightness variations, and frame-to-frame variability. This work extends arcjetCV [1] with filtering methods for tracking the evolving boundaries of hemispherical and ISO-Q test articles. Two boundary-tracking approaches based on Kalman filtering [2] were implemented. The first applies a point-wise Kalman filter followed by a Rauch–Tung–Striebel smoother [3] to individual boundary-point locations. The second applies the same filtering and smoothing framework to a reduced set of geometry-dependent shape parameters. The point-wise method reduces local frame-to-frame fluctuations while preserving spatial details along the detected boundary. The shape-parameter method provides a compact and geometrically constrained estimate of the sample contour. Figure 1 illustrates the two approaches for a hemispherical test article. Both methods improve temporal consistency and support more robust estimation of surface recession from arc-jet imagery. The two filtering approaches provide complementary representations of boundary evolution and improve the reliability of image-based recession measurements during arc-jet experiments.

recession measurement↗

Method for the Collection, Gravimetric and Chemical Analysis of Nonvolatile Residue (NVR) on Surfaces

Nonvolatile residue (NVR), sometimes referred to as molecular contamination is the term used for the total composition of the inorganic and high boiling point organic components in particulates and molecular films deposited on critical surfaces surrounding space structures, with the particulate and NVR contamination originating primarily from pre-launch operations. The "nonvolatile" suggestion from the terminology NVR implies that the collected residue will not experience much loss under ambient conditions. NVR has been shown to have a dramatic impact on the ability to perform optical measurements from platforms based in space. Such contaminants can be detected early by the controlled application of various detection techniques and contamination analyses. Contamination analyses are the techniques used to determine if materials, components, and subsystems can be expected to meet the performance requirements of a system. Of particular concern is the quantity of NVR contaminants that might be deposited on critical payload surfaces from these sources. Subsequent chemical analysis of the contaminant samples by infrared spectroscopy and gas chromatography mass spectrometry identifies the components, gives semi-quantitative estimates of contaminant thickness, indicates possible sources of the NVR, and provides guidance for effective cleanup procedures. In this report, a method for the collection and determination of the mass of NVR was generated by the authors at NASA Langley Research Center. This report describes the method developed and implemented for collecting NVR contaminants, and procedures for gravimetric and chemical analysis of the residue obtained. The result of this NVR analysis collaboration will help pave the way for Langley's ability to certify flight hardware outgassing requirements in support of flight projects such as Stratospheric Aerosol and Gas Experiment III (SAGE III), Clouds and the Earth's Radiant Energy System (CERES), Materials International Space Station Experiment - X (MISSE-X), and Doppler Aerosol Wind Lidar (DAWN).

Gordon, Keith↗

Acoustic emission monitoring crack propagation in single crystal silicon

The feasibility of acoustic emission (AE) monitoring of cracking and crack propagation in Si semiconductor materials was evaluated experimentally. A double torsion load relaxation method was employed wherein the propagation velocity and the AE levels in precracked (but not notched) boron-doped wafers were recorded simultaneously. A numerical model for the critical stress intensity factor (KIC) was used to relate the crack growth velocity, the instantaneous load and the load relaxation rate. All specimens were monitored with acoustic transducers at six points and examined with SEM after failure. The AE levels reached a peak amplitude of 70 dB at a KIC of 0.997 MNm to the -3/2 for cracking in the 111 plane. No AE was detected before the load reached the KIC, indicating that no subcritical crack growth occurs in Si. The results support the use of AE for monitoring crack propagation in crystal Si.

Chen, C. P.↗