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 487 records · Page 27

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗

Growth of Si spherical crystals and the surface oxidation (M-9)

Nearly 90 percent of semiconductor devices are produced with Si single crystals as the starting materials. For instance, the integrated circuits (IC), which are used in almost all electronic equipments such as TV, tape recorders, audio amplifiers, etc., are made after various processings of Si single crystal wafers. In these wafers, the same controlled amounts of impurities are added and the uniformities in their distributions are extremely important. Growth under microgravity makes it possible to eliminate the buoyancy-driven convection in the melt, which is one of the main origins of convections which results in non-uniformity of the impurity. Another source of convection is known as Marangoni convection which is driven on the free surface when a temperature gradient occurs. One of the merits of microgravity experimentation is that the detailed study of this convection becomes possible. Another important advantage of microgravity is that growth of crystals without a crucible is possible. This makes it possible to study melt growth without the strain which is usually introduced on the ground. Nevertheless, we should repeat and analyze many growth experiments in space to get reliable results. However, since in the FMPT, the time for the experiment is limited, we plan to carry out two kinds of very simple and basic experiments as the first step for the semiconductor growth experiment. In the first experiment, we use single crystal Si sphere as the starting material and as shown, this sphere is heated in the furnace at a slightly higher temperature than the melting point. After the melting front moves nearly half way to its center, the temperature is decreased to stop the melting and to start the growth from the seed for which we use the unmelted solid party of the sphere. The sphere is centered by quartz protuberances inside of the quartz crucible. There exists the possibility of temperature fluctuations being introduced when the molten sphere occasionally touches the protuberances. The total time needed for the melting and the growth processes is estimated to be 30 minutes. Infrared emission from the sphere is monitored in order to prevent the accidental loss of the central solid core. The schematical illustration of the second experient is shown. A single crystal, Si rod is used as the starting material. In the first stage, the rod is melted from one end to obtain a liquid sphere. In the second stage, the single crystal is grown by decreasing the temperature from the unmelted part of the rod which is used as the seed. The second experient somewhat resembles the Czochralski method used on the ground; however, in the space experiment, no crucible is employed and the temperature uniformity is much superior. In both experiments, phosphorous is doped to allow observation of the change in the shape of the liquid solid interface during crystal growth and the impurity striations, if any.

Nishinaga, Tatau↗

Worswick Hot Springs: A Radioactive Hydrothermal Field Site.

Introduction: We report on a systematic characterization of the environmental conditions at Worswick Hot Springs, a hydrothermal system in Idaho, USA. Because localized “hot spots” of elevated radiation and biofilms are easily accessible, various biological studies of radiation resistance and biosignature formation are possible, making this fieldsite relevant for analog field studies that consider microbiology, geo-chemistry, and ionizing radiation. In addition to Worswick being a natural radiation biology laboratory that may also be relevant for space biology applications, we assert that these unusual environmental conditions may inform us about locations on Mars that are also enriched in radioactive elements and their poten-tial for hosting biosignatures. Methods: We carried out repeated temperature and radiation measurements at the same hot spring locations to observe the system over time (Figure 2). Radiation: A Bicron Micro Analyst micro-r-meter (Bicron NE, Saint-Gobain Industrial Ceramics, Inc.) capable of sensing x-rays and gamma-rays (0-5000 µR/hr), was used to gather radiation data at twenty-four locations around and above the region of the two main stream channels. Water Temperature: A digital hand-held infrared thermometer (Oakton WD-39642-00 Mini-Temp Tester) was used for all measurements. ICP-MS aqueous geochemistry: Water samples were analyzed at the ISU Center for Archaeology, Materials and Applied Spectroscopy (CAMAS). Results: We have discovered localized areas of elevated radiation that are approximately 4 to 5 times greater than background radiation, and we have ob-served that both radioactivity and temperature of the spring waters vary over time. ICP-MS reveal the presence of thorium and uranium, which are correlated with elevated radioactivity. Several point sources of elevated radioactivity have been identified in both Stream A and Stream B.

Analog↗

Treatment to Control Adhesion of Silicone-Based Elastomers

Seals are used to facilitate the joining of two items, usually temporarily. At some point in the future, it is expected that the items will need to be separated. This innovation enables control of the adhesive properties of silicone-based elastomers. The innovation may also be effective on elastomers other than the silicone-based ones. A technique has been discovered that decreases the level of adhesion of silicone- based elastomers to negligible levels. The new technique causes less damage to the material compared to alternative adhesion mitigation techniques. Silicone-based elastomers are the only class of rubber-like materials that currently meet NASA s needs for various seal applications. However, silicone-based elastomers have natural inherent adhesive properties. This stickiness can be helpful, but it can frequently cause problems as well, such as when trying to get items apart. In the past, seal adhesion was not always adequately addressed, and has caused in-flight failures where seals were actually pulled from their grooves, preventing subsequent spacecraft docking until the seal was physically removed from the flange via an extravehicular activity (EVA). The primary method used in the past to lower elastomer seal adhesion has been the application of some type of lubricant or grease to the surface of the seal. A newer method uses ultraviolet (UV) radiation a mixture of UV wavelengths in the range of near ultraviolet (NUV) and vacuum ultraviolet (VUV) wavelengths.

deGroh, Henry C., III↗

Theoretical investigation of wave-vector-dependent analytical and numerical formulations of the interband impact-ionization transition rate for electrons in bulk silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E. Scholl, and M. Rudan (1993), and a more exact, numerical method of Y. Wang and K. F. Brennan (1994). Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formular, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with the completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to the overlap intergrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such the transition rate may be overestimated at those points for which 'near vertical,' small change in k, transitions occur.

Kolnik, Jan↗

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

The Large UV/Optical/Infrared Surveyor Decadal Mission Concept Thermal System Architecture

The Large Ultraviolet/Optical/Infrared (LUVOIR) Surveyor is one of four large strategic mission concept studies commissioned by NASA for the 2020 Decadal Survey in Astronomy and Astrophysics. Slated for launch to the second Lagrange point (L2) in the mid-to-late 2030s, LUVOIR seeks to directly image habitable exoplanets around sun-like stars, characterize their atmospheric and surface composition, and search for biosignatures, as well as study a large array of astrophysics goals including galaxy formation and evolution. Two observatory architectures are currently being considered which bound the trade-off between cost, risk, and scientific return: a 15-meter diameter segmented aperture primary mirror in a three-mirror anastigmat configuration, and an 8-meter diameter unobscured segmented aperture design. To achieve its science objectives, both architectures require milli-Kelvin level thermal stability over the optics, structural components, and interfaces to attain picometer wavefront RMS stability. A 270 Kelvin operational temperature was chosen to balance the ability to perform science in the near-infrared band and the desire to maintain the structure at a temperature with favorable material properties and lower contamination accumulation. This paper will focus on the system-level thermal designs of both LUVOIR observatory architectures. It will detail the various thermal control methods used in each of the major components - the optical telescope assembly, the spacecraft bus, the sunshade, and the suite of accompanying instruments - as well as provide a comprehensive overview of the analysis and justification for each design decision. It will additionally discuss any critical thermal challenges faced by the engineering team should either architecture be prioritized by the Astro2020 Decadal Survey process to proceed as the next large strategic mission for development.

Yang, Kan↗

The Large UV/Optical/Infrared Surveyor Decadal Mission Concept Thermal System Architecture

The Large Ultraviolet/Optical/Infrared (LUVOIR) Surveyor is one of four large strategic mission concept studies commissioned by NASA for the 2020 Decadal Survey in Astronomy and Astrophysics. Slated for launch to the second Lagrange point (L2) in the mid-to-late 2030s, LUVOIR seeks to directly image habitable exoplanets around sun-like stars, characterize their atmospheric and surface composition, and search for biosignatures, as well as study a large array of astrophysics goals including galaxy formation and evolution. Two observatory architectures are currently being considered which bound the trade-off between cost, risk, and scientific return: a 15-meter diameter segmented aperture primary mirror in a three-mirror anastigmat configuration, and an 8-meter diameter unobscured segmented aperture design. To achieve its science objectives, both architectures require milli-Kelvin level thermal stability over the optics, structural components, and interfaces to attain picometer wavefront RMS stability. A 270 Kelvin operational temperature was chosen to balance the ability to perform science in the near-infrared band and the desire to maintain the structure at a temperature with favorable material properties and lower contamination accumulation. This paper will focus on the system-level thermal designs of both LUVOIR observatory architectures. It will detail the various thermal control methods used in each of the major components - the optical telescope assembly, the spacecraft bus, the sunshade, and the suite of accompanying instruments - as well as provide a comprehensive overview of the analysis and justification for each design decision. It will additionally discuss any critical thermal challenges faced by the engineering team should either architecture be prioritized by the Astro2020 Decadal Survey process to proceed as the next large strategic mission for development.

Yang, Kan↗

Electron Beam Welding of Pure Tungsten Hex Cans for Nuclear Thermal Propulsion Engines

Nuclear thermal propulsion (NTP) is an in-space propulsion method currently being developed at the NASA Marshall Space Flight Center (MSFC). NTP systems are a high specific impulse (750–1,100 s), high thrust (15,000–250,000 lbf ) method of propulsion which have the potential to allow for faster transit times when optimizing for high ΔV. In the nuclear rocket engine, the heat from the nuclear fission reaction is transferred to a low molecular mass propellant (such as hydrogen). Hot propellant is expanded through a nozzle to generate thrust. Development of ceramic metal (cermet) fuel systems for NTP applications is currently ongoing at MSFC. In cermet fuel systems, ceramic fissile fuel particles such as uranium nitride or uranium dioxide are dispersed within a net-shaped, high-density structural matrix. The composite material is cladded by a protective metal structure to make up an NTP fuel element. Cladding materials must be able to withstand the demanding operating conditions required of the engine as well as retain a hermetic seal to allow for retention of fuel element structural integrity, prevent hydrogen attack or migration of the ceramic fuel, and limit release of fission products during operation. For NTP applications, tungsten is a prime material for both the metal matrix and cladding in cermet fuel systems because of its high melting point, high temperature strength, and compatibility with hot hydrogen. If a weld in tungsten with the capability of holding a hermetic seal is achievable, tungsten becomes a strong candidate for NTP applications. This Technical Memorandum focuses on determining the weldability of pure tungsten using electron beam welding (EBW). Tungsten appears well suited for NTP applications, but it has a high ductile to brittle transition temperature (DBTT) dependent upon chemical composition, structure/stress distribution, and mechanical conditions. Therefore, it is highly subject to brittle fracture. Because of its high susceptibility to brittle fracture, it is very difficult to weld. EBW was chosen for joining pure tungsten because of its low heat input compared to gas tungsten arc welding. Reduced heat input can be directly correlated with an increase in ductility of a tungsten weld. EBW is a high energy density welding process in which a stream of electrons penetrates a weld joint in a deep, narrow spike in contrast to a broad gas tungsten arc weld pool. The investigation initially focused on EBW of tungsten plates of both 0.01 in and 0.03 in thickness to determine if EBW could weld pure tungsten without the presence of visual defects—particularly cracking—in the welds. Variation in the weld procedure and post-weld heat treatment (PWHT) was used to improve the surface appearance of flat EBWs on a pure tungsten sheet. The investigation moved on to weld 0.05-in-thick hexagonal tungsten cans with a weld joint thickness of 0.025 in. The goal for welding the pure tungsten hex cans was to avoid any visual surface defects and generate a weld capable of a hermetic seal. This proved difficult. Cold welds commonly exhibited porosity that leaked air. Hot welds exhibited cracks, typically observed immediately after welding. Later welds were preheated to increase ductility and decrease the likelihood of through-thickness cracking. PWHT was used to arrest microcrack growth both in the flat weld samples and hexagonal weld samples.

Courtright, Z. S.↗

MAC/GMC 4.0 User's Manual: Example Problem Manual

This document is the third volume in the three volume set of User's Manuals for the Micromechanics Analysis Code with Generalized Method of Cells Version 4.0 (MAC/GMC 4.0). Volume 1 is the Theory Manual, Volume 2 is the Keywords Manual, and this document is the Example Problems Manual. MAC/GMC 4.0 is a composite material and laminate analysis software program developed at the NASA Glenn Research Center. It is based on the generalized method of cells (GMC) micromechanics theory, which provides access to the local stress and strain fields in the composite material. This access grants GMC the ability to accommodate arbitrary local models for inelastic material behavior and various types of damage and failure analysis. MAC/GMC 4.0 has been built around GMC to provide the theory with a user-friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermo-mechanical loading, generation of output results, and selection of architectures to represent the composite material, have been automated in MAC/GMC 4.0. Finally, classical lamination theory has been implemented within MAC/GMC 4.0 wherein GMC is used to model the composite material response of each ply. Consequently, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. This volume provides in-depth descriptions of 43 example problems, which were specially designed to highlight many of the most important capabilities of the code. The actual input files associated with each example problem are distributed with the MAC/GMC 4.0 software; thus providing the user with a convenient starting point for their own specialized problems of interest.

Bednarcyk, Brett A.↗

Stress Tuning of Laser Crystals

The topic of stress tunable laser crystals is addressed in this study with the purpose of determining the piezo-optic coefficients of a new laser material. This data was collected using a quadruple pass birefringence technique because of its high degree of sensitivity relative to the other methods examined including fringe shift analysis using a Mach-Zender interferometer. A green He-Ne laser was passed through a light chopper and Glan-Thompson prism before entering a crystal of Erbium doped Yttrium Aluminum Garnet (Er:YAG) (used in order to validate the experimental technique). The Er:YAG crystal is mounted in a press mechanism and the laser is quadruple passed through test specimen before being returned through the prism and the orthogonally polarized portion of the beam measured with a optical sensor. At a later stage, the Er:YAG crystal was replaced with a new crystal in order to determine the piezo-optic coefficients of this uncharacterized material. The applied load was monitored with the use of a 50 lb. load cell placed in line with the press. Light transmission readings were taken using a lock-in amplifier while load cell measurements were taken with a voltmeter from a 5 volt, 0.5 amp power supply. Despite the fact that an effective crystal press damping system was developed, size limitations precluded the use of the complete system. For this reason, data points were taken only once per full turn so as to minimize the effect of non uniform load application on the collected data. Good correlation was found in the transmission data between the experimentally determined Er:YAG and the previously known peizo-optic constants of non-doped crystal with which it was compared. The variation which was found between the two could be accounted for by the aforementioned presence of Erbium in the experimental sample (for which exact empirical data was not known). The same test procedure was then carried out on a Yttrium Gallium Aluminum garnet (YGAG) for the purpose of establishing values of its unknown piezo-optic constant tensor using experimentally collected transmission data. Significant variation between the piezo-optic constants of YAG and YGAG crystals was found however, the excellent data correlation of separate experimental runs carried out on the YGAG sample demonstrates the validity of these results. The data collected during the stressing of the YGAG was of high quality, however the amount of data collected was somewhat limited by a fracture of YGAG specimen which undoubted altered the crystalline lattice structure and hence precluded any further testing.

Carty, Atherton A.↗

Removing Contamination from Genesis Sapphire Collectors by Polishing

The Genesis Solar Wind Sample Re-turn Mission used a variety of materials as solar wind (SW) collectors. “Sapphire”, which is single crystal corundum, was ~6.7% of the collector material flown. But, a much higher percentage of Genesis sapphire collectors (SAP) survived the crash of the sample-return capsule as SAP is chemically inert, physically tough, and hard (MOHS 9, by definition). In fact, 39 sapphire fragments of bulk SW collector having areas ≥60 sq mm are available in JSC’s Genesis online catalog, whereas only 1 similarly sized silicon fragment is available (accessed 1/1/2020). In spite of the ease of accessing larger fragments, sapphire has not been popular among Genesis re-searchers because it is an electrical insulator (making SIMS more difficult) and its chemical inertness, high melting point and high hardness precludes the use of many mass-spectroscopic techniques. Recently, as alternative techniques for measuring Genesis samples are developed (e.g., synchrotron TXRF, INAA, laser ablation mass spectroscopy), more researchers are choosing to analyze sapphire. Thus, a number of researchers are actively studying how to best remove contamination from the Genesis sapphire surfaces. This report extends the task in Schmeling et al. that used polishing compounds to remove contamination from Genesis sapphire and focuses on contaminants (welded silicon, aluminum) that current chemical methods find recalcitrant.

M. Schmeling↗

Overall Traveling-Wave-Tube Efficiency Improved By Optimized Multistage Depressed Collector Design

Depressed Collector Design The microwave traveling wave tube (TWT) is used widely for space communications and high-power airborne transmitting sources. One of the most important features in designing a TWT is overall efficiency. Yet, overall TWT efficiency is strongly dependent on the efficiency of the electron beam collector, particularly for high values of collector efficiency. For these reasons, the NASA Glenn Research Center developed an optimization algorithm based on simulated annealing to quickly design highly efficient multistage depressed collectors (MDC's). Simulated annealing is a strategy for solving highly nonlinear combinatorial optimization problems. Its major advantage over other methods is its ability to avoid becoming trapped in local minima. Simulated annealing is based on an analogy to statistical thermodynamics, specifically the physical process of annealing: heating a material to a temperature that permits many atomic rearrangements and then cooling it carefully and slowly, until it freezes into a strong, minimum-energy crystalline structure. This minimum energy crystal corresponds to the optimal solution of a mathematical optimization problem. The TWT used as a baseline for optimization was the 32-GHz, 10-W, helical TWT developed for the Cassini mission to Saturn. The method of collector analysis and design used was a 2-1/2-dimensional computational procedure that employs two types of codes, a large signal analysis code and an electron trajectory code. The large signal analysis code produces the spatial, energetic, and temporal distributions of the spent beam entering the MDC. An electron trajectory code uses the resultant data to perform the actual collector analysis. The MDC was optimized for maximum MDC efficiency and minimum final kinetic energy of all collected electrons (to reduce heat transfer). The preceding figure shows the geometric and electrical configuration of an optimized collector with an efficiency of 93.8 percent. The results show the improvement in collector efficiency from 89.7 to 93.8 percent, resulting in an increase of three overall efficiency points. In addition, the time to design a highly efficient MDC was reduced from a month to a few days. All work was done in-house at Glenn for the High Rate Data Delivery Program. Future plans include optimizing the MDC and TWT interaction circuit in tandem to further improve overall TWT efficiency.

Vaden, Karl R.↗

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗

Probabilistic Multi-Factor Interaction Model for Complex Material Behavior

Complex material behavior is represented by a single equation of product form to account for interaction among the various factors. The factors are selected by the physics of the problem and the environment that the model is to represent. For example, different factors will be required for each to represent temperature, moisture, erosion, corrosion, etc. It is important that the equation represent the physics of the behavior in its entirety accurately. The Multi-Factor Interaction Model (MFIM) is used to evaluate the divot weight (foam weight ejected) from the external launch 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 were 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. The problem lies in how to represent the divot weight with a single equation. A unique solution to this problem is a multi-factor equation of product form. Each factor is of the following form (1 xi/xf)ei, where xi is the initial value, usually at ambient conditions, xf the final value, and ei the exponent that makes the curve represented unimodal that meets the initial and final values. The exponents are either evaluated by test data or by technical judgment. A minor disadvantage may be the selection of exponents in the absence of any empirical data. This form has been used successfully in describing the foam ejected in simulated space environmental conditions. Seven factors were required to represent the ejected foam. The exponents were evaluated by least squares method from experimental data. The equation is used and it can represent multiple factors in other problems as well; for example, evaluation of fatigue life, creep life, fracture toughness, and structural fracture, as well as optimization functions. The software is rather simplistic. Required inputs are initial value, final value, and an exponent for each factor. The number of factors is open-ended. The value is updated as each factor is evaluated. If a factor goes to zero, the previous value is used in the evaluation.

Abumeri, Galib H.↗

4U 1626-67 as Seen by Suzaku Before and After the 2008 Torque Reversal

Aims. The accretion-powered pulsar 4U 1626-67 experienced a new torque reversal at the beginning of 2008, after about 18 years of steadily spinning down. The main goal of the present work is to study this recent torque reversal that occurred in 2008 February. Methods. We present a spectral analysis of this source using two pointed observations performed by Suzaku in 2006 March and in 2010 September. Results. We confirm with Suzaku the presence of a strong emission-line complex centered on 1 keV, with the strongest line being the hydrogen-like Ne Lya at 1.025(3) keV. We were able to resolve this complex with up to seven emission lines. A dramatic increase of the intensity of the Ne Lya line after the 2008 torque reversal occurred, with the equivalent width of this line reaching almost the same value measured by ASCA in 1993. We also report on the detection of a cyclotron line feature centered at approximately 37 keV. In spite of the fact that an increase of the X-ray luminosity (0.5-100keV) of a factor of approximately 2.8 occurred between these two observations, no significant change in the energy of the cyclotron line feature was observed. However, the intensity of the approximately 1 keV line complex increased by an overall factor of approximately 8. Conclusions. Our results favor a scenario in which the neutron star in 4U 1626-67 accretes material from a geometrically thin disk during both the spin-up and spin-down phases.

Camero-Arranz, A.↗