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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1962

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↗
DOE OSTI2020

Charge density wave order on a $\textit{π}$-flux square lattice

The effect of electron-phonon coupling on Dirac fermions was recently explored numerically on a honeycomb lattice, leading to precise quantitative values for the finite temperature and quantum critical points. In this paper, we use the unbiased determinant quantum Monte Carlo method to study the Holstein model on a half-filled staggered-flux square lattice and compare the results with the honeycomb lattice geometry, presenting results for a range of phonon frequencies 0.1 ≤ $\omega$ ≤ 2.0. In this work, we find that the interactions give rise to charge density wave order, but only above a finite coupling strength $\gamma_{\text{crit}}$. The transition temperature is evaluated and presented in a $T_{c^{–λ}}$ phase diagram. An accompanying mean-field theory calculation also predicts the existence of a quantum phase transition, but at a substantially smaller coupling strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Comparison of vapour pressure osmometry, freezing point osmometry and direct membrane osmometry for determining the osmotic pressure of concentrated solutions

Osmometry is the analytical method of measuring osmotic pressure ($π$) and related colligative properties. Accurately measuring $π$ is essential for applications that use osmosis, including drug delivery systems and water filtration methods such as forward osmosis. The three most common forms of directly or indirectly measuring $π$ are vapour pressure osmometry (VPO), freezing point osmometry (FPO), and direct membrane osmometry (DMO). In this report, the performance of VPO, FPO, and DMO are compared for various aqueous solutions, including small molecule and polymer solutes, and potentially more difficult samples such as concentrated, carbonated, or viscous solutions. Furthermore, potential problems with each of the three methods are examined. Good practices can avoid inaccuracies in DMO due to external concentration polarization and solute diffusion through the membrane. Similarly, it is important to be aware of and avoid weaknesses in VPO and FPO in order to obtain reliable results. Even with care, however, the three methods still give different values of $π$ for concentrated aqueous solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Structural and Improper Ferroelectric Properties of TbInO 3 Single Crystal Grown by Laser Floating Zone

The honeycomb TbInO 3 has attracted wide research attention due to its fascinating physical properties. However, TbInO 3 single crystal was difficult to grow owing to the high melting point and serious volatilization of indium during the crystal growth. In this study, the volatilization of the indium element was effectively suppressed by controlling the growth atmosphere and pressure. The excess ratio of indium oxide was determined, and pure hexagonal TbInO 3 crystal was obtained by the laser floating zone method. Systematic studies on the crystal structure and optical and ferroelectric properties were carried out. The structure distortion resulted in the improper geometric ferroelectric revealed by single crystal diffraction and Raman spectrum measurements. The topological vortex domains and P-E hysteresis loop demonstrated the presence of ferroelectricity. TbInO 3 crystal has great potential application in vortex memory.

36 MATERIALS SCIENCE↗
NASA NTRS2000

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↗
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗
DOE OSTI2022

Stirred-Reactor Coupon Analysis: An International Round Robin Study

The objective of this task was to determine the precision of the SRCA technique when used to determine the dilute condition corrosion rate. To this end, an interlaboratory round robin study was conducted per the instructions in ASTM Practice E691 to measure the precision with which the SRCA test method can be conducted. Twelve independent labs from eleven different institutions each evaluated four glass compositions in three different conditions. The ASTM procedures recommend at least 6 labs participate in a round robin testing the same 3 materials in the same conditions to determine precision. In this case, 12 labs each performed 12 independent tests. This was only possible thanks to the multi-glass testing capability of the SRCA test. A total of 108 duplicate pairs were used to calculate the repeatability of the tests, with the same glass tested in the same vessel, producing as identical conditions as possible for the replicates. These test results were quite tightly clustered, with a median difference from the average value of the pair of only 2.9%. Based on the calculations outlined in ASTM E177-20 and a measured standard deviation of 4.74%, the intralaboratory repeatability limit (r) was calculated to be within 13.3% of the expected value with 95% confidence level. The reproducibility limit (R) of the test was examined using all 277 data points from the round robin. Because of the differences in dissolution rates due to pH variability and intrinsically for the 12 conditions tested, the reproducibility limits for each condition and overall were calculated from the percent relative residual value for each test. The SRCA test is expected to be reproducible within 64% of the expected value with a 95% confidence level.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1993

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2013

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↗
DOE OSTI2026

LaueMatching: an approach for rapid and robust indexing of Laue diffraction patterns

Traditional Laue diffraction pattern indexing often struggles with noisy data, weak signals, peak overlap and missing reflections, particularly from complex or deformed microstructures. Here, we introduce LaueMatching, a high-throughput indexing algorithm designed to overcome these limitations. LaueMatching utilizes a fundamentally different approach based on direct pattern correlation: experimentally pre-processed images are compared against a comprehensive pre-computed library of simulated diffraction patterns corresponding to a dense grid of possible orientations. This approach bypasses the need for explicit peak identification and fitting, steps that are often a failure point for traditional methods. The algorithm rapidly and robustly indexes multiple crystallographic orientations and crystal systems simultaneously, even from challenging patterns. LaueMatching's effectiveness and accuracy have been rigorously tested and validated on diverse experimental (Ni, Al, EuAl 2 O 4 ) and simulated diffraction patterns, demonstrating high-fidelity orientation refinement. Code to implement this approach on both CPU and GPU resources can be downloaded from https://github.com/AdvancedPhotonSource/LaueMatching.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗
DOE OSTI2023

Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.

36 MATERIALS SCIENCE↗
NASA NTRS1994

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↗
NASA NTRS1994

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↗
NASA NTRSDate not supplied

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↗