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 775 records · Page 43

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.↗

Radiative capture rates at deep defects from electronic structure calculations

We present a methodology to calculate radiative carrier capture coefficients at deep defects in semiconductors and insulators from first principles. Electronic structure and lattice relaxations are accurately described with hybrid density functional theory. Calculations of capture coefficients provide an additional validation of the accuracy of these functionals in dealing with localized defect states. We also discuss the validity of the Condon approximation, showing that even in the event of large lattice relaxations the approximation is accurate. We test the method on GaAs:VGa-TeAs and GaN:CN, for which reliable experiments are available, and demonstrate very good agreement with measured capture coefficients.

36 MATERIALS SCIENCE↗

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)↗

Geometry-aware framework for deep energy method: An application to structural mechanics with hyperelastic materials

Here, in this work, we introduce a novel physics-informed framework named the Geometry-Aware Deep Energy Method (GADEM) for solving structural mechanics problems on different geometries. As the weak form of the physical system equation (or the energy-based approach) has demonstrated clear advantages compared to the strong form for solving solid mechanics problems, GADEM employs the weak form and aims to infer the solution on multiple shapes of geometries. Integrating a geometry-aware framework into an energy-based method results in an effective physics-informed deep learning model in terms of accuracy and computational cost. Different ways to represent the geometric information and to encode the geometric latent vectors are investigated in this work. We introduce a loss function of GADEM which is minimized based on the potential energy of all considered geometries. An adaptive learning method is also employed for the sampling of collocation points to enhance the performance of GADEM. We present some applications of GADEM to solve solid mechanics problems, including a loading simulation of a toy tire involving contact mechanics and large deformation hyperelasticity. The numerical results of this work demonstrate the remarkable capability of GADEM to infer the solution on various and new shapes of geometries using only one trained model.

97 MATHEMATICS AND COMPUTING↗

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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.↗

Bayesian estimations of orientation distribution functions from small-angle scattering enable direct prediction of mechanical stress in anisotropic materials

Properties of soft materials are influenced by their anisotropic structuring under nonequilibrium fields. Although anisotropic structure-property relationships have been extensively explored theoretically, comparison to experiments requires determination of the microstructural orientation probability distribution function (OPDF) of microstructural elements. Small angle scattering (SAS) measurements encode information about the OPDF, but tools to navigate this connection are incomplete. Here, we develop and validate an explicit framework to link arbitrary OPDFs to SAS measurements. Specifically, we propose, validate, and apply a method, maximum a posteriori scattering inference (MAPSI), whereby the OPDF may be obtained from SAS measurements using a Bayesian estimation method. Using this method, we obtain estimates of the full 3D OPDF for two model semidilute fd-virus (rodlike) dispersions at concentrations that are approximately equal to and twice the overlap concentration. From the OPDF, we calculate its second and fourth moments and compare these to predictions for a dilute suspension of rigid rods and to a recent theory for semidilute suspensions. Finally, we use both the theoretical and measured moments to calculate the stress, both for dilute and semidilute suspensions. These predictions are not only compared to each other, but also to measured values of the shear stress, and point to new insights into the behavior of suspensions of highly elongated particles in the transition between dilute and semidilute behavior. We also use this new framework to provide perspective on the connection between scalar parameterizations of scattering and the OPDF that have frequently been used in the past. Here, the new tools developed in this work provide an unprecedented path toward experimental validation of dynamical theories of rodlike colloids and polymers, and for measurement of nonequilibrium structures and stresses of other complex fluids and soft materials with SAS.

36 MATERIALS SCIENCE↗