SEARCH · Engineering Papers

Results for “ANISOTROPIC PLATE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

DOE OSTI2024

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗

NASA NTRS2003

Analysis Method for Inelastic, Adhesively Bonded Joints with Anisotropic Adherends

A one-dimensional analysis method for evaluating adhesively bonded joints composed of anisotropic adherends and adhesives with nonlinear material behavior is presented in the proposed paper. The strain and resulting stress field in a general, bonded joint overlap are determined by using a variable-step, finite-difference solution algorithm to iteratively solve a system of first-order differential equations. Applied loading is given by a system of combined extensional, bending, and shear forces that are applied to the edge of the joint overlap. Adherends are assumed to behave as linear, cylindrically bent plates using classical laminated plate theory that includes the effects of first-order transverse shear deformation. Using the deformation theory of plasticity and a modified von-Mises yield criterion, inelastic material behavior is modeled in the adhesive layer. Results for the proposed method are verified against previous results from the literature and shown to be in excellent agreement. An additional case that highlights the effects of transverse shear deformation between similar adherends is also presented.

Smeltzer, Stanley S., III↗

DOE OSTI2020

Nonuniform Ionic and Electronic Transport of Ceramic and Polymer/Ceramic Hybrid Electrolyte by Nanometer-Scale Operando Imaging for Solid-State Battery

Replacing the liquid electrolyte in lithium batteries with solid-state ion conductor is promising for next-generation energy storage that is safe and has high energy density. In this work, nanometer-resolution ionic and electronic transport imaging of Li 3 PS 4 (LPS), a solid-state electrolyte (SSE), is reported. This nm resolution is achieved by using a logarithm-scale current amplifier that enhances the current sensitivity to the fA range. Large fluctuations of ion current—one to two orders of magnitude on the LPS and on the LPS region of a polymer/LPS bulk hybrid SSE—that must be mitigated to eliminate Li dendrite formation and growth, are found. This ion current fluctuation is understood in terms of highly anisotropic transport kinetic barriers along the different crystalline axes due to different grain orientations in the polycrystalline and glass ceramic materials. The results on the bulk hybrid SSE show a sharp transition of ionic and electronic transport at the LPS/polymer boundary and decreases in average ionic current with decreasing polyimine particle size and with extensive cycling. The results elucidate the mechanism of polyimine extension into interparticles to prevent Li dendrite growth. This work opens up novel characterization of charge transport, which relates to Li plating and stripping for solid-state-batteries.

Li3PS4↗

NASA NTRS2013

Micromachined Active Magnetic Regenerator for Low-Temperature Magnetic Coolers

A design of an Active Magnetic Regenerative Refrigeration (AMRR) system has been developed for space applications. It uses an innovative 3He cryogenic circulator to provide continuous remote/distributed cooling at temperatures in the range of 2 K with a heat sink at about 15 K. A critical component technology for this cooling system is a highly efficient active magnetic regenerator, which is a regenerative heat exchanger with its matrix material made of magnetic refrigerant gadolinium gallium garnet (GGG). Creare Inc. is developing a microchannel GGG regenerator with an anisotropic structured bed for high system thermal efficiency. The regenerator core consists of a stack of thin, single-crystal GGG disks alternating with thin polymer insulating layers. The insulating layers help minimize the axial conduction heat leak, since GGG has a very high thermal conductivity in the regenerator s operating temperature range. The GGG disks contain micro channels with width near 100 micrometers, which enhance the heat transfer between the circulating flow and the refrigerant bed. The unique flow configuration of the GGG plates ensures a uniform flow distribution across the plates. The main fabrication challenges for the regenerator are the machining of high-aspect-ratio microchannels in fragile, single-crystal GGG disks and fabrication and assembly of the GGG insulation layers. Feasibility demonstrations to date include use of an ultrashort- pulse laser to machine microchannels without producing unacceptable microcracking or deposition of recast material, as shown in the figure, and attachment of a thin insulation layer to a GGG disk without obstructing the flow paths. At the time of this reporting, efforts were focused on improving the laser machining process to increase machining speed and further reduce microcracking.

Chen, Weibo↗

NASA NTRS1991

Measurements of induced radioactivity in some LDEF samples

Twenty-six stainless steel trunnion samples, five aluminum end support retainer plate samples, two aluminum keel plate samples, and two titanium clips were analyzed. The shielded high purity germanium detectors used had efficiencies of 33, 54, and 80 pcts. at 1332 keV. Detector efficiencies as a function of energy and corrections for self absorption in the samples were determined with calibrated sources, unactivated control samples, and calculations. Several measurements were made on most samples. In the trunnion samples, Mn-54 and Co-57 were seen and limits were obtained for other isotopes. The results agree well with 1-D activation calculations for an anisotropic trapped proton model. In the aluminum samples, Na-22 and Be-7 were seen. Other results are presented.

Moss, C. E.↗

NASA NTRS1992

Measurements of induced radioactivity in some LDEF samples

Twenty-six stainless steel trunnion samples, five aluminum end support retainer plate samples, two aluminum keel plate samples, and two titanium clips were analyzed. The shielded high purity germanium detectors used had efficiencies of 33, 54, and 80 percent at 1332 keV. Detector efficiencies as a function of energy and corrections for self-absorption in the samples were determined with calibrated sources and unactivated control samples. Several measurements were made on most samples. In the trunnion samples, Mn-54 and Co-57 were seen and limits were obtained for other isotopes. The results agree well with 1-D activation calculations for an anisotropic trapped proton model. In the aluminum and titanium samples, Na-22 was detected. Other results are presented.

Moss, C. E.↗

NASA NTRS2010

The Fifteen-Year Attitude History of the Wide Field Planetary Camera 2 Radiator and Collection Efficiencies for Micrometeoroids and Orbital Debris

An examination of the Hubble Space Telescope (HST) Wide Field Planetary Camera 2 (WFPC-2) radiator assembly was conducted at NASA Goddard Space Flight Center (GSFC) during the summer of 2009. Immediately apparent was a distinct biasing of the largest 45 impact features towards one side of the radiator, in contrast to an approximately uniform distribution of smaller impacts. Such a distribution may be a consequence of the HST s attitude history and pointing requirements for the cold radiator, or of environmental effects, such as an anisotropic distribution of the responsible population in that size regime. Understanding the size-dependent spatial distribution of impact features is essential to the general analysis of these features. We have obtained from GSFC a 15 minute temporal resolution record of the state vector (Earth Centered Inertial position and velocity) and HST attitude, consisting of the orientation of the velocity and HST-sun vectors in HST body coordinates. This paper reviews the actual state vector and attitude history of the radiator in the context of the randomly tumbling plate assumption and assesses the statistical likelihood (or collection efficiency) of the radiator for the micrometeoroid and orbital debris environments. The NASA Marshall Space Flight Center s Meteoroid Environment Model is used to assess the micrometeoroid component. The NASA Orbital Debris Engineering Model (ORDEM) is used to model the orbital debris component. Modeling results are compared with observations of the impact feature spatial distribution, and the relative contribution of each environmental component are examined in detail.

Anz-Meador, Phillip D.↗

NASA NTRSDate not supplied

Towards an Automated Unstructured Grid Adaptation Workflow with VULCAN

Early work is presented for an unstructured grid adaptation workflow with VULCAN and refine. Anisotropic simplex grids are iteratively adapted to match a Riemannian metric tensor field describing desired mesh spacing. The Riemannian metric tensor field is obtained from Hessians of CFD solution output scalar sensor fields; both Mach number and static temperature sensor fields are explored. In addition, we describe a Newton-method-based solver recently implemented in VULCAN utilizing Jacobian-Free-Newton-Krylov that can be used to increase flow solver automation on early grids in the adadptation process. Hypersonic flow solutions are presented on a high Reynolds number flat plate and wall heat flux is compared against a highly resolved structured solution. Additionally, complex shock boundary-layer interaction is explored in a high Mach number compression corner and complex 3D flow phenomena are evaluated on the Boundary Layer Transition (BOLT) vehicle.

Matthew O'Connell↗

DOE OSTI2021

Ion temperature anisotropy in the tokamak scrape-off layer

Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field B → . The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e. the temperatures are assumed to be isotropic. Additionally, in certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the two-dimensional (2D) transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

Analyses of composite fuselage structure under various loading conditions

With the characteristics of high strength/weight and stiffness/weight ratios in advanced composite materials, composite plate structures have been successfully applied to secondary load-carrying structural components in the aerospace industry for the past two decades. Recently, filament wound composite shells are being considered for design of primary fuselage structures. To implement the design the structural response under various loading conditions needs to be predicted accurately. Therefore, it is essential to establish a rigorous analytical solution in the area of composite laminated shells. A closed form solution is presented that predicts the response of a composite shell subjected to internal pressure, axial tension, bending, and torsion. The material of the shell is assumed to be general cylindrically anisotropic. Based on the theory of cylindrical anisotropic elasticity coupled partial differential governing equations are developed using Lehknitskii's stress function approach. The general expressions for the stresses and displacements in the composite cylinders under these loading conditions is discussed. Three examples (45) off-axis unidirectional, (45/45) unsymmetric, and (45/45)(sub s) symmetric angle-ply fiber-reinforced laminated shells are shown to illustrate the effect of radius-to-thickness ratios, coupling, and stacking sequence.

Yuan, Fuh-Gwo↗

NASA NTRS1993

Three-dimensional analysis of anisotropic spatially reinforced structures

The material-adaptive three-dimensional analysis of inhomogeneous structures based on the meso-volume concept and application of deficient spline functions for displacement approximations is proposed. The general methodology is demonstrated on the example of a brick-type mosaic parallelepiped arbitrarily composed of anisotropic meso-volumes. A partition of each meso-volume into sub-elements, application of deficient spline functions for a local approximation of displacements and, finally, the use of the variational principle allows one to obtain displacements, strains, and stresses at anypoint within the structural part. All of the necessary external and internal boundary conditions (including the conditions of continuity of transverse stresses at interfaces between adjacent meso-volumes) can be satisfied with requisite accuracy by increasing the density of the sub-element mesh. The application of the methodology to textile composite materials is described. Several numerical examples for woven and braided rectangular composite plates and stiffened panels under transverse bending are considered. Some typical effects of stress concentrations due to the material inhomogeneities are demonstrated.

Bogdanovich, Alexander E.↗

NASA NTRS1993

A displacement gradient BEM for accurate stress computation near boundaries in 2-D anisotropic problems

A displacement gradient method for 2D anisotropic elasticity problems is presented, which effectively minimizes the boundary layer effect through a two-step procedure. First, the boundary integral equations are solved for the unknown boundary displacements and tractions. Second, a direct integral equation for displacement gradients is developed in terms of boundary tractions. Three methods based on different evaluation procedures and locations for determining the displacement gradients are proposed. In the first method the displacement gradients are averaged at nodes common to adjacent elements. The second method stores the gradients element-wise. In the third method, the gradients are evaluated at the nodes of discontinuous elements. The three methods are applied to near-isotropic plates with circular and elliptic cutouts. It is concluded that all three methods can yield accurate stress distributions.

Sistla, R.↗

NASA NTRS2011

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Manchiraju, Sivom↗

NASA NTRS2000

Determination of the Fracture Parameters in a Stiffened Composite Panel

A modified J-integral, namely the equivalent domain integral, is derived for a three-dimensional anisotropic cracked solid to evaluate the stress intensity factor along the crack front using the finite element method. Based on the equivalent domain integral method with auxiliary fields, an interaction integral is also derived to extract the second fracture parameter, the T-stress, from the finite element results. The auxiliary fields are the two-dimensional plane strain solutions of monoclinic materials with the plane of symmetry at x(sub 3) = 0 under point loads applied at the crack tip. These solutions are expressed in a compact form based on the Stroh formalism. Both integrals can be implemented into a single numerical procedure to determine the distributions of stress intensity factor and T-stress components, T11, T13, and thus T33, along a three-dimensional crack front. The effects of plate thickness and crack length on the variation of the stress intensity factor and T-stresses through the thickness are investigated in detail for through-thickness center-cracked plates (isotropic and orthotropic) and orthotropic stiffened panels under pure mode-I loading conditions. For all the cases studied, T11 remains negative. For plates with the same dimensions, a larger size of crack yields larger magnitude of the normalized stress intensity factor and normalized T-stresses. The results in orthotropic stiffened panels exhibit an opposite trend in general. As expected, for the thicker panels, the fracture parameters evaluated through the thickness, except the region near the free surfaces, approach two-dimensional plane strain solutions. In summary, the numerical methods presented in this research demonstrate their high computational effectiveness and good numerical accuracy in extracting these fracture parameters from the finite element results in three-dimensional cracked solids.

Lin, Chung-Yi↗

NASA NTRS1997

Effect of Free Stream Turbulence and Other Vortical Disturbances on a Laminar Boundary Layer

This paper is concerned with the effect of free-stream turbulence on the pretransitional flat-plate boundary layer. It is assumed that either the turbulence Reynolds number or the downstream distance (or both) is small enough so that the flow can be linearized. The dominant disturbances in the boundary layer, which are of the Klebanoff type, are governed by the linearized unsteady boundary-region equations, i.e., the Navier Stokes equations with the streamwise derivatives neglected in the viscous and pressure-gradient terms. The turbulence is represented as a superposition of vortical free-stream Fourier modes, and the corresponding individual Fourier component solutions to the boundary-region equations are obtained numerically. The results are then superposed to compute the root mean square of the fluctuating streamwise velocity in the boundary layer produced by the actual free-stream turbulence. The calculated boundary-layer disturbances are in good quantitative agreement with the experimentally observed Klebanoff modes when strong low-frequency anisotropic effects are included in the free-stream turbulence spectrum. We discuss some additional effects that may need to be accounted for in order to obtain a complete description of the Klebanoff modes.

Goldstein, Marvin E.↗

NASA NTRS2004

Material Modeling of Space Shuttle Leading Edge and External Tank Materials For Use in the Columbia Accident Investigation

Upon the commencement of the analytical effort to characterize the impact dynamics and damage of the Space Shuttle Columbia leading edge due to External Tank insulating foam, the necessity of creating analytical descriptions of these materials became evident. To that end, material models were developed of the leading edge thermal protection system, Reinforced Carbon Carbon (RCC), and a low density polyurethane foam, BX-250. Challenges in modeling the RCC include its extreme brittleness, the differing behavior in compression and tension, and the anisotropic fabric layup. These effects were successfully included in LS-DYNA Material Model 58, *MAT_LAMINATED_ COMPOSITE_ FABRIC. The differing compression and tension behavior was modeled using the available damage parameters. Each fabric layer was given an integration point in the shell element, and was allowed to fail independently. Comparisons were made to static test data and coupon ballistic impact tests before being utilized in the full scale analysis. The foam's properties were typical of elastic automotive foams; and LS-DYNA Material Model 83, *MAT_FU_CHANG_FOAM, was successfully used to model its behavior. Material parameters defined included strain rate dependent stress-strain curves for both loading and un-loading, and for both compression and tension. This model was formulated with static test data and strain rate dependent test data, and was compared to ballistic impact tests on load-cell instrumented aluminum plates. These models were subsequently utilized in analysis of the Shuttle leading edge full scale ballistic impact tests, and are currently being used in the Return to Flight Space Shuttle re-certification effort.

Carney, Kelly↗

NASA NTRS2001

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗