SEARCH · Search NASA
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.
In this talk, we will present interface-aware methods for diffusion equation. The focus is on the necessity of interface-aware, mixed cells, interface reconstruction, discontinuity of material diffusion coefficients, order of accuracy in time, property for large time steps, and simpler interface-aware methods. Numerical examples are provided to demonstrate the features on these aspects of the methods
The design and sizing of adhesively bonded joints has always been a major bottleneck in the design of composite vehicles. Dense finite element (FE) meshes are required to capture the full behavior of a joint numerically, but these dense meshes are impractical in vehicle-scale models where a course mesh is more desirable to make quick assessments and comparisons of different joint geometries. Analytical models are often helpful in sizing, but difficulties arise in coupling these models with full-vehicle FE models. Therefore, a joint FE was created which can be used within structural FE models to make quick assessments of bonded composite joints. The shape functions of the joint FE were found by solving the governing equations for a structural model for a joint. By analytically determining the shape functions of the joint FE, the complex joint behavior can be captured with very few elements. This joint FE was modified and used to consider adhesives with functionally graded material properties to reduce the peel stress concentrations located near adherend discontinuities. Several practical concerns impede the actual use of such adhesives. These include increased manufacturing complications, alterations to the grading due to adhesive flow during manufacturing, and whether changing the loading conditions significantly impact the effectiveness of the grading. An analytical study is conducted to address these three concerns. Furthermore, proof-of-concept testing is conducted to show the potential advantages of functionally graded adhesives. In this study, grading is achieved by strategically placing glass beads within the adhesive layer at different densities along the joint. Furthermore, the capability to model non-linear adhesive constitutive behavior with large rotations was developed, and progressive failure of the adhesive was modeled by re-meshing the joint as the adhesive fails. Results predicted using the joint FE was compared with experimental results for various joint configurations, including double cantilever beam and single lap joints.
The diffraction by a material discontinuity in a thick dielectric/ferrite layer is considered by modeling the layer as a distributed current sheet obeying generalized sheet transition conditions (GSTC's). The sheet currents are then formulated and solved via the standard dual integral equation approach. This yields the diffracted field in terms of unknown constants which underscore the non-uniqueness of the GSTC current sheet representation. The constants are dependent on the geometry and properties of the discontinuity and are determined by enforcing field continuity across the material junction. This requires the field internal to the slab which are determined from the external ones via analytic continuity. Results are given which validate the solution and demonstrate the importance of the constants.
Basaltic rocks are thought to constitute a volumetrically significant rock type on the Moon, Mercury, Mars, and Venus, in addition to the Earth. Spacecraft images of surfaces with known or suspected basaltic composition on these bodies, particularly on Venus, indicate that these rocks have been deformed in the brittle regime to form faults and perhaps dilatant cracks, in addition to folding and more distributed types of deformation. Predictions of brittle fracture or other types of deformation are made by comparing calculated stresses from a tectonic model to some criterion for rock strength. Common strength criteria used in the planetary science literature for near-surface deformation include a Griffith tensile-strength criterion for intact rock, a Mohr envelope for intact basalt, and a brittle strength envelope based on Byerlee's law of rock frictional resistance. However, planetary terrains of basaltic composition consist of much more than just intact basaltic rock. The aggregate basaltic material, termed the 'rock mass,' consists of both the intact rock and the associated fracture, faults, lithologic contacts, and other discontinuous surfaces. A basaltic rock mass is the relevant material for which strength properties must be defined and calculated model stresses must be compared to in order to more accurately predict brittle deformation. For example, the various strengths of a rock mass are less than that of intact material of the same composition. This means that tectonic models which compare stresses to intact failure strengths overestimate the stresses required for fracture and so underestimate the extent and magnitude of brittle deformation predicted in these models. On the other hand, rock mass shear strength can be greater than that predicted from Byerlee's law. The concept of rock mass strength is central to many engineering design studies in which calculated stresses are used to predict brittle fracture, and this experience indicates that brittle strength envelopes which assume properties for intact rock (Griffith parabolas) or sliding along a single, continuous surface (Byerlee's law) inadequately characterize the tensile, compressive, and shear strengths of rock masses. The criterion adopted here to relate stresses to rock mass fracture is based on a Griffith-type curve for tensile normal stress and a concave downward curve for compressive normal stress. It is the only available criterion that explicitly considers the weakening effects of discontinuities within the rock mass on the stress state required for fracture.
Fiber length distribution (FLD), in part, governs mechanical properties in discontinuous fiber composites, yet manual measurement methods limit the high-throughput characterization needed for materials design optimization. This study compares deep learning segmentation approaches for automated FLD measurement in large-field microscopy, evaluating how method choice affects the microstructural descriptors used in structure-property-processing relationships. A critical challenge is that high-resolution microscopy images (10,000×10,000 pixels) must be tiled for deep learning analysis, fragmenting fibers at boundaries. We demonstrate that segmentation method proves crucial for measurement accuracy. For example, instance segmentation with Slicing Aided Hyper Inference (SAHI) preserves individual fiber integrity across tiles while semantic segmentation prioritizes speed. Comparing against manual measurement of extracted carbon fibers, YOLOv11-SAHI matched manual ground truth (238 μm weighted mean) with 40x speedup (4.5 vs 167 minutes per image). U-Net provides rapid quantification although it is at the cost of reduced accuracy due only reliably measuring stand-alone fibers. Our comparative analysis reveals that instance segmentation with SAHI better preserves length measurements while semantic segmentation prioritizes speed, providing empirical guidance for method selection. The characterization provides essential inputs for mechanical property prediction models and inverse design workflows, accelerating composite materials development cycles.
This work focused on preferentially aligning discontinuous carbon fibers in wet-laid or air-laid processes. It is well known that aligned fibers provides higher directional strength and stiffness. Discontinuous fibers further allow higher degree of draw and formability as the gaps in the fibers allow for higher material movement. The current processes are limited in their ability to align carbon fibers during processing. The aligned fibers have several benefits - (a) in applications where chopped fibers can replace continuous fibers for targeted strength and stiffness metrics, but at a substantially reduced cost; (b) they can tolerate deeper draws than continuous fiber composites in thermo-stamping and compression molding processes; (c) they can be tailored for pultrusion and unidirectional applications. Although pultrusion is primarily a process that adopts continuous fibers, stitch bonded entangled discontinuous fibers can provide unique intermediates. This is analogous to natural coir fibers which get aligned and entangled to produce ropes/rods for example, (d) they can be processed in cross-ply and multi-directional formats, like composite laminates. In this work Neenah Paper partnered with IACMI, UT and ORNL to evaluate structure-process-property relationships with Zoltek carbon fiber. A few process parameters such as machine speed, weight basis, fiber length, effect of fiber sizing, direction of mat lay-up etc. were investigated. The produced mats were converted to thermoplastic composite laminates using polyamide 6 (PA6, nylon) resin. The specific objective of this project is to produce a wet-laid nonwoven carbon fiber mat with a high degree of unidirectional fiber alignment, using discontinuous carbon fibers. The report provides details about the processing, characterization, and lower-upper bound properties.
Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.
The improvement in weld quality by the friction stir welding (FSW) process invented by TWI of Cambridge, England, patented in 1991, has prompted investigation of this process for advanced structural materials including Al metal matrix composite (Al-MMC) materials. Such materials can have high specific stiffness and other potential beneficial properties for the extreme environments in space. Developments of discontinuous reinforced Al-MMCs have found potential space applications and the future for such applications is quite promising. The space industry has recognized advantages of the FSW process over conventional welding processes such as the absence of a melt zone, reduced distortion, elimination of the need for shielding gases, and ease of automation. The process has been well proven for aluminum alloys, and work is being carried out for ferrous materials, magnesium alloys and copper alloys. Development work in the FSW welding process for joining of Al-MMCs is relatively recent and some of this and related work can be found in referenced research publications. NASA engineers have undertaken to spear head this research development work for FSW process investigation of Al-MMCs. Some of the reported related work has pointed out the difficulty in fusion welding of particulate reinforced MMCs where liquid Al will react with SiC to precipitate aluminum carbide (Al4C3). Advantages of no such reaction and no need for joint preparation for the FSW process is anticipated in the welding of Al-MMCs. The FSW process has been best described as a combination of extrusion and forging of metals. This is carried out as the pin tool rotates and is slowly plunged into the bond line of the joint as the pin tool's shoulder is in intimate contact with the work piece. The material is friction-stirred into a quality weld.
The sub-surface structures of water-bearing icy worlds (e.g., Europa, Enceladus, and Titan) and the structures of icy regions on Mars contain critical information about their origin and evolution. These bodies have been identified as high-priority targets in the NRC Planetary Science Decadal Survey, which noted a lack of technology readiness as one impediment to related missions. An Elastic Wave Analyzer for Icy Sub-surfaces (EWAIS) instrument is being developed to acquire data that can provide answers to fundamental science questions related to these bodies. The EWAIS instrument generates and receives elastic waves using an array of piezoelectric transducers capable of operating at temperatures as low as 30K. A developed breadboard acquires reflected signals from discontinuities in the traveled materials’ impedance (the product of wave velocity and density) along the wave path. Analysis of the acquired data can be used to determine the elastic properties, presence of cracks, locations of cavities and other discontinuities, as well as the thicknesses of ice and liquid layers. The EWAIS instrument is being developed with a novel dual-frequency transmitter array that will enable reaching distances of kilometers through ice, as well as resolving smaller acoustic reflectors nearer the EWAIS instrument. The selected dual frequencies correspond to two depth ranges that can be analyzed, with resolutions of tens of centimeters in the region of tens to hundreds of meters and tens of meters on the kilometers scale, respectively. The EWAIS will be applicable for operation from any in-situ platform, including surface assets (e.g., lander and rover) and sub-surface assets (e.g., melting probes).
Epitaxial Sc x Al 1- x N thin films of ~100 nm thickness grown on metal polar GaN substrates are found to exhibit significantly enhanced relative dielectric permittivity (ε r ) values relative to AlN. ε r values of ~17–21 for Sc mole fractions of 17%–25% ( x = 0.17–0.25) measured electrically by capacitance–voltage measurements indicate that Sc x Al 1- x N has the largest relative dielectric permittivity of any existing nitride material. Since epitaxial Sc x Al 1- x N layers deposited on GaN also exhibit large polarization discontinuity, the heterojunction can exploit the in situ high-K dielectric property to extend transistor operation for power electronics and high-speed microwave applications.
We have initiated an experimental program to develop an ultrasound system for non-intrusively imaging the thermal field in opaque fluids under an externally imposed temperature gradient. Many industrial processes involve opaque fluids, such as molten metals, semiconductors, and polymers, often in situations in which thermal gradients are important. For example, one may wish to understand semiconductor crystal growth dynamics in a Bridgman apparatus. Destructive testing of the crystal after the process is completed gives only indirect information about the fluid dynamics of the formation process. Knowledge of the coupled thermal and velocity fields during the growth process is then essential. Most techniques for non-intrusive velocity and temperature measurement in fluids are optical in nature, and hence the fluids studied must be transparent. In some cases (for example, LDV (laser Doppler velocimetry) and PIV (particle imaging velocimetry)) the velocities of small neutrally buoyant seed particles suspended in the fluid, are measured. Without particle seeding one can use the variation of the index of refraction of the fluid with temperature to visualize, through interferometric, Schlieren or shadowgraph techniques, the thermal field. The thermal field in turn gives a picture of the pattern existing in the fluid. If the object of study is opaque, non-optical techniques must be used. In this project we focus on the use of ultrasound, which propagates easily through opaque liquids and solids. To date ultrasound measurements have almost exclusively relied on the detection of sound scattered from density discontinuities inside the opaque material of interest. In most cases it has been used to visualize structural properties, but more recently the ultrasound Doppler velocimeter has become available. As in the optical case, it relies on seed particles that scatter Doppler shifted sound back to the detector. Doppler ultrasound techniques are, however, not useful for studying convective fluid flow in crystal growth, because particle seeding is unacceptable and flow velocities are typically too low to be resolved, and may be even lower in microgravity conditions where buoyancy forces are negligible. We will investigate a different use of ultrasound to probe the flows of opaque fluids non-intrusively and without the use of seed particles: our goal is to ultrasonically visualize the thermal field of opaque fluids with relatively high spatial resolution. The proposed technique relies upon the variation of sound speed with temperature of the fluid. A high frequency ultra-sound pulse passing through a fluid-filled chamber will traverse the chamber in a time determined by the relevant chamber dimension and the temperature of the fluid through which the pulse passes. With high time-resolution instrumentation that compares the excitation signal with the received pulse we can detect the influence of the fluid temperature on the pulse travel time. This is effectively an interferometric system, which in its optical form is an extremely sensitive approach to measuring thermal fields in fluids. Moreover, the temperature dependence of sound velocity in liquid metals is comparable to the temperature dependence of the speed of light required for accurate interferometric thermal images in transparent fluids. With an array of transducers scanned electronically a map of the thermal field over the chamber could be produced. An alternative approach would be to use the ultrasound analog of the shadowgraph technique. In the optical version, collimated light passes through the fluid, where it is focused or defocused locally by temperature field induced variations of the index of refraction. The resulting image reveals the thermal field through the spatial pattern of light intensity variations. By analogy, an ultrasound plane wave traversing an opaque fluid sample would be also locally focused or defocused depending on the speed of sound variations, giving rise to spatial variations in sound intensity that will reveal the thermal field pattern. These approaches could be applied to any situation in which temperature differences are expected to occur, and will rapidly provide information about the flow that simply cannot be obtained by any current intrusive or non-intrusive diagnostic technique. As materials processing in microgravity matures it will become increasingly important to have available simple and versatile diagnostic tools, such as we will develop, for studying the flows of opaque fluids under thermal forcing.
In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.
Correlating superconducting qubit performance with advanced materials analysis is a key strategy for improving coherence. Existing diagnostics for key properties, such as dielectric loss, structural discontinuity, and interface heterogeneity, often rely on destructive electron microscopy or low-throughput millikelvin measurements. Here, in this study, we demonstrate noninvasive terahertz (THz) nano-imaging/spectroscopy of encapsulated niobium transmon qubits as a high-throughput proxy for performance evaluation. We identify large variations in sidewall near-field signals, implicating sidewall loss and discontinuity as major coherence limiters, and also use THz hyperspectral line scans to probe dielectric responses and field participation at Al junction interfaces.
In a low- or zero-gravity environment the free surface of a liquid can behave in striking, unexpected ways. For example, in a partially filled container, a free surface that is well behaved under terrestrial conditions can rise to an arbitrarily large height or even fail to exist when gravity is absent. An answer is sought to the central mathematical questions: Under what conditions can free surfaces exist and what are their properties? -- and experimental questions of what means can be devised to observe and to measure the surfaces quantitatively. Current activity focuses on material selection and the design of optical diagnostic methods for in-space experiments. These experiments will test mathematical predictions of discontinuous transition from existence to nonexistence of capillary free surfaces in certain geometric configurations under zero gravity conditions. Liquids and container materials suitable for achieving the required contact-angle range and optical properties are being investigated, as well as the effects of contaminants and fluid motion on contact angle. The optical diagnostic technique investigation concerns laser-induced fluorescence, including emphasis on data on acquisition, sources and detectors, reliability, data management, and computer control.
The basic principles of bifurcation analysis have been established by several investigators, however several issues remain unresolved, specifically how do stress level, grain size distribution, and boundary conditions affect general bifurcation phenomena in pressure sensitive and dilatant materials. General geometrical and kinematics conditions for moving surfaces of discontinuity was derived and applied to problems of instability of solids. In 1962, the theoretical framework of bifurcation by studying the acceleration waves in elasto-plastic (J2) solids were presented. Bifurcation analysis for more specific forms of constitutive behavior was examined by studying localization in pressure-sensitive, dilatant materials, however, analyses were restricted to plane deformation states only. Bifurcation analyses were presented and applied to predict shear band formations in sand under plane strain condition. The properties of discontinuous bifurcation solutions for elastic-plastic solids under axisymmetric and plane strain loading conditions were studied. The study focused on theory, but also references and comparisons to experiments were made. The current paper includes a presentation of a summary of bifurcation analyses for biaxial and triaxial (axisymmetric) loading conditions. The Coulomb model is implemented using incremental piecewise scheme to predict the constitutive relations and shear band inclination angles. Then, a comprehensive evaluation of bifurcation phenomena is presented based on data from triaxial experiments performed under microgravity conditions aboard the Space Shuttle under very low effective confining pressure (0.05 to 1.30 kPa), in which very high peak friction angles (47 to 75 degrees) and dilatancy angles (30 to 31 degrees) were measured. The evaluation will be extended to include biaxial experiments performed on the same material under low (10 kPa) and moderate (100 kPa) confining pressures. A comparison between the behavior under biaxial and triaxial loading conditions will be presented, and related issues concerning influence of confining pressure will be discussed.
Computer aided manufacturing (CAM) techniques for directed energy deposition (DED) affect the material properties of the manufactured component based on the scan strategy used. In this study, we investigate the material characteristics of turning-style toolpath strategies to generate axisymmetric components with additive manufacturing (AM), referred to in this research as additive turning. This novel approach leverages existing CAM technology for turning, where the component rotates around a stationary cutting tool, to generate toolpath trajectories for DED with varying wall-thicknesses and controlled deposition angles. This strategy allows for entire components to be deposited in one continuous deposition, resulting in reduced cycle-time and improved material usage efficiency compared to conventional AM strategies where the beam is switched off at the end of every layer. Results from this study show that the use of additive turning can produce over 99 % dense components with less variation and anisotropy in texture and hardness, as well as a lower variation in elongation to failure when compared to conventional strategies. This research highlights that various CAM strategies could be deployed for AM to improve process efficiency or enable localized control over part performance.
Adhesives with functionally graded material properties are being considered for use in adhesively bonded joints to reduce the peel stress concentrations located near adherend discontinuities. Several practical concerns impede the actual use of such adhesives. These include increased manufacturing complications, alterations to the grading due to adhesive flow during manufacturing, and whether changing the loading conditions significantly impact the effectiveness of the grading. An analytical study is conducted to address these three concerns. An enhanced joint finite element, which uses an analytical formulation to obtain exact shape functions, is used to model the joint. Furthermore, proof of concept testing is conducted to show the potential advantages of functionally graded adhesives. In this study, grading is achieved by strategically placing glass beads within the adhesive layer at different densities along the joint.
A three phase program has been conducted to investigate DuPont's Long Discontinuous Fiber (LDF) composites. Additional tests were conducted to compare LDF composites against toughened thermosets and a baseline thermoset system. Results have shown that the LDF AS4/PEKK offers improved interlaminar (flange bending) strength with little reduction in mechanical properties due to the discontinuous nature of the fibers. In the third phase, a series of AS4/PEKK LDF C-section curved frames (representing a typical rotorcraft light frame) were designed, manufactured and tested. Specimen reconsolidation after 'stretch forming' and frame thickness were found to be key factors in this light frame's performance. A finite element model was constructed to correlate frame test results with expected strain levels determined from material property tests. Adequately reconsolidated frames performed well and failed at strain levels at or above baseline thermoset material test strains. Finally a cost study was conducted which has shown that the use of LDF for this frame would result in a significant cost savings, for moderate to large lot sizes compared with the hand lay-up of a thermoset frame.