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.
Anisotropic materials ductile fracture involving crack initation and propagation, using Dugdale mathematical model
There is a growing need for inorganic anisotropic particles in a variety of materials science applications. Structural, optical, and electrical properties can be greatly augmented by the fabrication of composite materials with anisotropic microstructures or with anisotropic particles uniformly dispersed in an isotropic matrix. Examples include structural composites, magnetic and optical recording media, photographic film, certain metal and ceramic alloys, and display technologies including flat panel displays. While considerable progress has been made toward developing an understanding of the synthesis of powders composed of monodispersed, spherical particles, these efforts have not been transferred to the synthesis of anisotropic nanoparticles. The major objective of the program is to develop a fundamental understanding of the growth of anisotropic particles at organic templates, with emphasis on the chemical and structural aspects of layered organic assemblies that contribute to the formation of anisotropic inorganic particles.
A generalized reliability model was developed for use in the design of structural components made from brittle, homogeneous anisotropic materials such as single crystals. The model is based on the Weibull distribution and incorporates a variable strength distribution and any equivalent stress failure criteria. In addition to the reliability model, an energy based failure criterion for elastically anisotropic materials was formulated. The model is different from typical Weibull-based models in that it accounts for strength anisotropy arising from fracture toughness anisotropy and thereby allows for strength and reliability predictions of brittle, anisotropic single crystals subjected to multiaxial stresses. The model is also applicable to elastically isotropic materials exhibiting strength anisotropy due to an anisotropic distribution of flaws. In order to develop and experimentally verify the model, the uniaxial and biaxial strengths of a single crystal nickel aluminide were measured. The uniaxial strengths of the <100> and <110> crystal directions were measured in three and four-point flexure. The biaxial strength was measured by subjecting <100> plates to a uniform pressure in a test apparatus that was developed and experimentally verified. The biaxial strengths of the single crystal plates were estimated by extending and verifying the displacement solution for a circular, anisotropic plate to the case of a variable radius and thickness. The best correlation between the experimental strength data and the model predictions occurred when an anisotropic stress analysis was combined with the normal stress criterion and the strength parameters associated with the <110> crystal direction.
The anisotropic generalization of the parameterized post-Newtonian (PPN) gravitational metric field is made for the case of theories with energy and momentum conservation laws. Such an anisotropic metric field will generally result in two-tensor or bimetric theories of gravity in an anisotropic universe. New anisotropic 3 x 3 spatial PPN matrices are introduced into the general metric expansion. Earth gravimeter measurements strongly restrict some anisotropies, while anisotropic inertial and gravitational mass for celestial bodies result from other combinations of the PPN matrices.
An existing volume grid generation procedure, AFLR3, was successfully modified to generate anisotropic tetrahedral elements using a directional metric transformation defined at source nodes. The procedure can be coupled with a solver and an error estimator as part of an overall anisotropic solution adaptation methodology. It is suitable for use with an error estimator based on an adjoint, optimization, sensitivity derivative, or related approach. This offers many advantages, including more efficient point placement along with robust and efficient error estimation. It also serves as a framework for true grid optimization wherein error estimation and computational resources can be used as cost functions to determine the optimal point distribution. Within AFLR3 the metric transformation is implemented using a set of transformation vectors and associated aspect ratios. The modified overall procedure is presented along with details of the anisotropic transformation implementation. Multiple two-and three-dimensional examples are also presented that demonstrate the capability of the modified AFLR procedure to generate anisotropic elements using a set of source nodes with anisotropic transformation metrics. The example cases presented use moderate levels of anisotropy and result in usable element quality. Future testing with various flow solvers and methods for obtaining transformation metric information is needed to determine practical limits and evaluate the efficacy of the overall approach.
An in-depth tutorial on the constitutive equations for elastic, anisotropic materials is presented. Basic concepts are introduced that are used to characterize materials, and notions about how anisotropic material deform are presented. Hooke s law and the Duhamel-Neuman law for isotropic materials are presented and discussed. Then, the most general form of Hooke s law for elastic anisotropic materials is presented and symmetry requirements are given. A similar presentation is also given for the generalized Duhamel-Neuman law for elastic, anisotropic materials that includes thermal effects. Transformation equations for stress and strains are presented and the most general form of the transformation equations for the constitutive matrices are given. Then, specialized transformation equations are presented for dextral rotations about the coordinate axes. Next, concepts of material symmetry are introduced and criteria for material symmetries are presented. Additionally, engineering constants of fully anisotropic, elastic materials are derived from first principles and the specialized to several cases of practical importance.
Adaptive unstructured mesh techniques have a limited, but growing impact on production analysis workflows to control discretization error for reliable simulation results. Multiple independent implementations of flow solvers, anisotropic metric construction methods, and anisotropic mesh adaptation mechanics have matured. Goal-based metrics target estimated error in output functions, such as lift and drag, through the guidance of an adjoint solution. A unification of goal-based anisotropic metrics is presented for steady viscous flows, which is an active area of research. These goal-based metrics drive robust and efficient anisotropic mesh adaptation for the calculation of output functions. The super-convergent functional output error behavior of stabilized finite-element methods is exploited without a formal proof, and evidence of super-convergence is shown in numerical experiments. Mesh adapted drag and lift outputs for two simple bodies in compressible viscous flow show convergence of error to less than a single drag count. Asymptotic behavior established for relatively coarse meshes shows the efficiency of this goal-based metric when compared to solution interpolation error control and expert-guided meshing. Anisotropic mesh adaptation techniques are applied to a transport aircraft in a high-lift configuration where variation between approaches decreases with mesh refinement, but asymptotic behavior is not observed with available resources.
For electromagnetic wave propagation and scattering in an anisotropic random medium, the Dyson equation for the mean field and the Bethe-Salpeter equation for the correlation or the covariance of the field were derived. With the random permittivity expressed in a general anisotropic form, the bilocal and the nonlinear approximations are employed to solve the Dyson equation, and the ladder approximation to solve the Bethe-Salpeter equation. The mean dyadic Green's function for a two-layer anisotropic random medium with arbitrary three-dimensional correlation functions has been investigated with the zeroth-order solutions to the Dyson equation under the nonlinear approximation. The effective propagation constants are calculated for the four characteristic waves associated with the coherent vector fields, propagating in an anisotropic random-medium layer, which are the ordinary and extraordinary waves with upward- and downward-propagating vectors.
A hybrid stress finite element method is developed for accurate stress and vibration analysis of problems in linear anisotropic elasticity. A modified form of the Hellinger-Reissner principle is formulated for dynamic analysis and an algorithm for the determination of the anisotropic elastic and compliance constants from experimental data is developed. These schemes were implemented in a finite element program for static and dynamic analysis of linear anisotropic two dimensional elasticity problems. Specific numerical examples are considered to verify the accuracy of the hybrid stress approach and compare it with that of the standard displacement method, especially for highly anisotropic materials. It is that the hybrid stress approach gives much better results than the displacement method. Preliminary work on extensions of this method to three dimensional elasticity is discussed, and the stress shape functions necessary for this extension are included.
An efficient preconditioned conjugate gradient (PCG) technique and a computational procedure are presented for the analysis of symmetric anisotropic structures. The technique is based on selecting the preconditioning matrix as the orthotropic part of the global stiffness matrix of the structure, with all the nonorthotropic terms set equal to zero. This particular choice of the preconditioning matrix results in reducing the size of the analysis model of the anisotropic structure to that of the corresponding orthotropic structure. The similarities between the proposed PCG technique and a reduction technique previously presented by the authors are identified and exploited to generate from the PCG technique direct measures for the sensitivity of the different response quantities to the nonorthotropic (anisotropic) material coefficients of the structure. The effectiveness of the PCG technique is demonstrated by means of a numerical example of an anisotropic cylindrical panel.
A study of deformation mechanisms in the single crystal superalloy PWA 1480 was conducted. Monotonic and cyclic tests were conducted from 20 to 1093 C. Both (001) and near-(123) crystals were tested, at strain rates of 0.5 and 50 percent/minute. The deformation behavior could be grouped into two temperature regimes: low temperatures, below 760 C; and high temperatures, above 820 to 950 C depending on the strain rate. At low temperatures, the mechanical behavior was very anisotropic. An orientation dependent CRSS, a tension-compression asymmetry, and anisotropic strain hardening were all observed. The material was deformed by planar octahedral slip. The anisotropic properties were correlated with the ease of cube cross-slip, as well as the number of active slip systems. At high temperatures, the material was isotropic, and deformed by homogeneous gamma by-pass. It was found that the temperature dependence of the formation of superlattice-intrinsic stacking faults was responsible for the local minimum in the CRSS of this alloy at 400 C. It was proposed that the cube cross-slip process must be reversible. This was used to explain the reversible tension-compression asymmetry, and was used to study models of cross-slip. As a result, the cross-slip model proposed by Paidar, Pope and Vitek was found to be consistent with the proposed slip reversibility. The results were related to anisotropic viscoplastic constitutive models. The model proposed by Walter and Jordan was found to be capable of modeling all aspects of the material anisotropy. Temperature and strain rate boundaries for the model were proposed, and guidelines for numerical experiments were proposed.
An ultrasonic wave is applied to an anisotropic sample material in an initial direction and an angle of flux deviation of the ultrasonic wave front is measured from this initial direction. This flux deviation angle is induced by the unknown applied load. The flux shift is determined between this flux deviation angle and a previously determined angle of flux deviation of an ultrasonic wave applied to a similar anisotropic reference material under an initial known load condition. This determined flux shift is then compared to a plurality of flux shifts of a similarly tested, similar anisotropic reference material under a plurality of respective, known load conditions, whereby the load applied to the particular anisotropic sample material is determined. A related method is disclosed for determining the fiber orientation from known loads and a determined flux shift.
To study the effect of cancellations within long-range interactions on local isotropy at the small scales, we calculate explicitly the degree of cancellation in distant interactions in the simulations of Yeung & Brasseur and Yeung, Brasseur & Wang using the single scale disparity parameter 's' developed by Zhou. In the simulations, initially isotropic simulated turbulence was subjected to coherent anisotropic forcing at the large scales and the smallest scales were found to become anisotropic as a consequence of direct large-small scale couplings. We find that the marginally distant interactions in the simulation do not cancel out under summation and that the development of small-scale anisotropy is indeed a direct consequence of the distant triadic group, as argued by Yeung, et. al. A reduction of anisotropy at later times occurs as a result of the isotropizing influences of more local energy-cascading triadic interactions. Nevertheless, the local-to-nonlocal triadic group persists as an isotropizing influence at later times. We find that, whereas long-range interactions, in general, contribute little to net energy transfer into or out of a high wavenumber shell k, the anisotropic transfer of component energy within the shell increases with increasing scale separations. These results are consistent with results by Zhou, and Brasseur & Wei, and suggest that the anisotropizing influences of long range interactions should persist to higher Reynolds numbers. The residual effect of the forced distant group in this low-Reynolds number simulation is found to be forward cascading, on average.
This paper presents a model with ellipsoidal scatterers for applications to polarimetric remote sensing of anisotropic layered media at microwave frequencies. The physical configuration includes an isotropic layer covering an anisotropic layer above a homogeneous half space. The isotropic layer consists of randomly oriented spheroids. The anisotropic layer contains ellipsoidal scatterers with a preferential vertical alignment and random azimuthal orientations. Effective permittivities of the scattering media are calculated with the strong fluctuation theory extended to account for the nonspherical shapes and the scatterer orientation distributions. On the basis of the analytic wave theory, dyadic Green's functions for layered media are used to derive polarimetric backscattering coefficients under the distorted Born approximation. The ellipsoidal shape of the scatterers gives rise to nonzero cross-polarized returns from the untilted anisotropic medium in the first-order approximation. Effects of rough interfaces are estimated by an incoherent addition method. Theoretical results and experimental data are matched at 9 GHz for thick first-year sea ice with a bare surface and with a snow cover at Point Barrow, Alaska. The model is then used to study the sensitivity of polarimetric backscattering coefficients with respect to correlation lengths representing the geometry of brine inclusions. Polarimetric signatures of bare and snow-covered sea ice are also simulated based on the model to investigate effects of different scattering mechanisms.
The objective of this project is to develop and demonstrate a capability for computing stress intensity factors in generally anisotropic materials. These objectives have been met. The primary deliverable of this project is this report and the information it contains. In addition, we have delivered the source code for a subroutine that will compute stress intensity factors for anisotropic materials encoded in both the C and Python programming languages and made available a version of the FRANC3D program that incorporates this subroutine. Single crystal super alloys are commonly used for components in the hot sections of contemporary jet and rocket engines. Because these components have a uniform atomic lattice orientation throughout, they exhibit anisotropic material behavior. This means that stress intensity solutions developed for isotropic materials are not appropriate for the analysis of crack growth in these materials. Until now, a general numerical technique did not exist for computing stress intensity factors of cracks in anisotropic materials and cubic materials in particular. Such a capability was developed during the project and is described and demonstrated herein.
Recently, a metric-based mesh adaptation capability has been developed for the advanced compressible flow solver US3D. Metric-based mesh adaptation is a solution informed mesh adaptation strategy that typically improves the resolution of the mesh along anisotropic flow features, like strong shock waves or shear layers. An anisotropic error indicator is reconstructed based on a provided Computational Fluid Dynamics (CFD) solution which prescribes the locally desired element size and orientation. The benefits of this mesh adaptation strategy is twofold: first, it simplifies the manual meshing labor that typically needs to be carried out by the CFD engineer because the mesh adaptation is informed by the provided solution. Second, the computational cost is kept to a minimum while improving the resolution of the simulation. The aim here is to demonstrate how anisotropic mesh adaptation can be employed to improve the predictions of surface pressure for atmospheric entry vehicles. In this work, we compare CFD simulations that use this newly developed anisotropic mesh adaptation capability for US3D with pressure data that was collected by Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) during the atmospheric entry phase of Mars2020.
A metric-based mesh adaptation capability has been developed for the advanced compressible flow solver US3D. Metric-based mesh adaptation is a solution informed mesh adaptation strategy that improves the resolution of the mesh along anisotropic flow features, like strong shock waves or shear layers. An anisotropic error indicator is reconstructed based on a provided Computational Fluid Dynamics (CFD) solution which prescribes the locally desired element size and orientation. The benefits of this mesh adaptation strategy is twofold: first, it simplifies the manual meshing labor that typically needs to be carried out by the CFD engineer because the mesh adaptation is informed by the provided solution. Second, the computational cost is kept to a minimum while improving the resolution of the simulation. The aim here is to demonstrate how anisotropic mesh adaptation can be employed to improve the predictions of surface pressure for atmospheric entry vehicles. In this work, we compare CFD simulations that use this newly developed anisotropic mesh adaptation capability for US3D with pressure data that was collected by Mars Entry, Descent, and Landing Instrumentation (MEDLI2) suite.
Anisotropic friction, deformation, and fracture studies were conducted with /0001/, /10(-1)0/, and /11(-2)0/ silicon carbide surfaces in sliding contact with diamond. The experiments were conducted with loads of 0.1, 0.2, and 0.3 N at a sliding velocity of 3 mm/min in mineral oil or in dry argon at room temperature. The 1010 direction on the basal /0001/ plane exhibits the lowest coefficient of friction and the greatest resistance to abrasion for silicon carbide. Anisotropic friction and deformation of the /0001/, /10(-1)0/, and /11(-2)0/ silicon carbide surfaces are primarily controlled by the slip system /10(-1)0/ 11(-2)0. The anisotropic fracture during sliding on the basal plane is due to surface cracking along /10(-1)0/ and subsurface cracking along /0001/. The fracture during sliding on the /11(-2)0/ or /10(-1)0/ surfaces is due to surface cracking along /0001/ and /11(-2)0/ or /10(-1)-1 and to subsurface cracking along /10(-1)0/.