SEARCH · Engineering Papers
Results for “Density matrix”
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.
A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling
Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.
Microstructural development in DED stainless steels: applying welding models to elucidate the impact of processing and alloy composition
Austenitic stainless steel microstructures produced by directed energy deposition (DED)are analogous to those developed during welding, particularly high energy density welding. To better understand microstructural development during DED, theories of microstructural evolution,which have been established to contextualize weld microstructures, are applied in this study to microstructural development in DED austenitic stainless steels. Phenomenological welding models that describe the development of oxide inclusions, compositional microsegregation, ferrite,matrix austenite grains, and dislocation substructures are utilized to clarify microstructural evolution during deposition of austenitic stainless steels. Two different alloys, 304L and 316L, arecompared to demonstrate the broad applicability of this framework for understanding microstmctural development during the DED process. Despite differences in grain morphology and solidification mode for these two alloys (which can be attributed to compositional differences),similar tensile properties are achieved. It is the fine-scale compositional segregation and dislocation structures that ultimately determine the strength of these materials. The evolution of microsegregation and dislocation structures is shown to be dependent on the rapid solidification and thermomechanical history of the DED processing method and not the composition of the starting material.
Chemical compatibility of hollow ceramic cenospheres as thermal insulation for high-temperature thermal energy storage applications with molten nitrate salt
An effective insulation material that is both thermally and chemically stable in molten salts could transform the design of molten-salt-based thermal energy storage (TES) tanks. Most current molten salt TES tanks hold the metallic tank structure in direct contact with hot salt inventory, a design which leads to thermal expansion of the tank and triggers stresses that can lead to thermomechanical failures. With an internal insulation to lower the temperature at the tank structure, the extent of thermal expansion can be reduced, thereby reducing expansion-induced stresses and allowing for consideration of lower-cost tank structure materials. Conventional insulation materials are either 1) too porous and allow molten salts to permeate into the matrix, which significantly increases the thermal conductivity or 2) too dense and have a thermal conductivity that cannot provide sufficient thermal insulation. This paper presents an alternative thermal insulation concept using cenospheres which have an alumino-silicate structure. The cost analysis suggests that the low-density cenospheres can be one of the cheapest materials to provide cost-effective thermal insulation. The chemical compatibility of cenospheres is investigated in molten 60 wt% NaNO3/40 wt% KNO3 salt which is close to industrial-grade Solar Salt. This paper shows that diffusion of the sodium and potassium cations from the salt into the cenospheres occurs based on weight analysis, energy dispersive spectroscopy (EDS), X-ray diffraction (XRD) and Fourier-transform infrared (FTIR) spectroscopy. The cation diffusion breaks the bridging oxygen bonds and causes volume expansion of the microstructure which is responsible for the failure of the cenosphere particles. The chemical composition of the cenospheres is found to affect their compatibility with molten nitrate salt. A cenosphere product with low iron content showed the best compatibility with an average survival rate of 77.9% +/- 9.8% after 7 days of immersion in the molten nitrate salt. While even the low-iron cenospheres appear to require protection from direct molten salt contact, their slow degradation rate, closed-cell porosity, and low cost hold potential for effective use as internal tank insulation.
Multi-energy reconstructions, central electron temperature measurements, and early detection of the birth and growth of runaway electrons using a versatile soft x-ray pinhole camera at MST
A multi-energy soft x-ray pinhole camera has been designed, built, and deployed at the Madison Symmetric Torus to aid the study of particle and thermal transport, as well as MHD stability physics. This novel imaging diagnostic technique employs a pixelated x-ray detector in which the lower energy threshold for photon detection can be adjusted independently on each pixel. The detector of choice is a PILATUS3 100 K with a 450 μm thick silicon sensor and nearly 100 000 pixels sensitive to photon energies between 1.6 and 30 keV. An ensemble of cubic spline smoothing functions has been applied to the line-integrated data for each time-frame and energy-range, obtaining a reduced standard-deviation when compared to that dominated by photon-noise. The multi-energy local emissivity profiles are obtained from a 1D matrix-based Abel-inversion procedure. Central values of Te can be obtained by modeling the slope of the continuum radiation from ratios of the inverted radial emissivity profiles over multiple energy ranges with no a priori assumptions of plasma profiles, magnetic field reconstruction constraints, high-density limitations, or need of shot-to-shot reproducibility. In tokamak plasmas, a novel application has recently been tested for early detection, 1D imaging, and study of the birth, exponential growth, and saturation of runaway electrons at energies comparable to 100 × Te,0; thus, early results are also presented.
Processing and Mechanical Properties of Various Zirconia/Alumina Composites for Fuel Cell Applications
Various electrolyte materials for solid oxide fuel cells were fabricated by hot pressing 10 mol% yttria-stabilized zirconia (10-YSZ) reinforced with two different forms of alumina, particulates and platelets, each containing 0 to 30 mol% alumina. Flexure strength and fracture toughness of both particulate and platelet composites at ambient temperature increased with increasing alumina content, reaching a maximum at 30 mot% alumina. For a given alumina content, strength of particulate composites was greater than that of platelet composites, whereas, the difference in fracture toughness between the two composite systems was negligible. No virtual difference in elastic modulus and density was observed for a given alumina content between particulate and platelet composites. Thermal cycling up to 10 cycles between 200 to 1000 C did not show any effect on strength degradation of the 30 mol% platelet composites, indicative of negligible influence of CTE mismatches between YSZ matrix and alumina grains.
Properties and microstructure evolution of silicon nitride and zirconium nitride following Ni ion irradiation
We report that silicon nitride and zirconium nitride have been proposed as potential materials for multiple nuclear applications (inert matrix fuels, accident tolerant fuels, space nuclear power, fusion reactor design), yet knowledge on their behavior under irradiation remains limited. Ion irradiations were performed using 15 MeV Ni 5+ ions on Si 3 N 4 and ZrN samples, with midrange doses (around 3 µm) from 1 to 50 dpa and temperatures from 300 to 700°C. Volumetric lattice swelling was determined by grazing incidence X-ray diffraction, defect production and evolution were tracked using Transmission Electron Microscopy, and nanoindentation was performed to quantify the ceramics’ mechanical properties evolution. The results from these irradiation studies on nitride ceramics help fill the current gap present in the literature. Behavior consistent with past work on irradiated Si 3 N 4 was observed with respect to mechanical properties and defect formation up to 15 dpa and 500°C. Failure of the grain boundary sintering aid in Si 3 N 4 was observed above these conditions. Different behavior was observed in both nitrides at 50 dpa and 700°C, where lattice swelling increased past potential saturation values. Unreported cavity formation was witnessed in both materials under all irradiation conditions, with stable number density and slight size increase above 15 dpa. The mechanism for the cavity formation remains to be determined.
Solving the electronic structure problem for over 100000 atoms in real space
Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.
Effect of heat treatment on the tensile behavior of selective laser melted Ti-6Al-4V by in situ X-ray characterization
Selective Laser Melted Ti-6Al-4V (as-SLMed) exhibits decreased yield strength, increased work hardening, and increased ductility after heat treatment at 730 degrees C (HT-730) or 900 degrees C (HT-900) for 2 h. To understand the change of mechanical properties, in situ high energy X-ray diffraction (HEXRD) is used to examine the phase composition, load partitioning, slip system activity, and dislocation density evolution in all three specimens. The as-SLMed specimen is dominated by martensitic α'. After heat treatment, α' partly or fully decomposes into α+β, reducing the yield strength. In HT-730, beta precipitates with confined size show much higher lattice strain than the α'/α matrix during deformation; in HT-900, the lattice strain difference is mostly eliminated. This is a key reason for the increased ductility in HT-900. From the anisotropic lattice strain development, basal slip is identified as the easiest slip system in α'/α. Using an elasto-plastic self-consistent (EPSC) model, the critical resolved shear stress ratio between prismatic slip and basal slip (CRSS prism /CRSS basal ) is estimated to be 1.31 and 1.16 in the as -SLMed and the HT-900 specimens, respectively. Further, the alpha phase in HT-900 is able to activate multiple slip systems and accumulate more dislocations during plastic deformation. This explains why HT-900 has better ductility and higher work hardening rate than the other two specimens.
Fission track astrology of three Apollo 14 gas-rich breccias
The three Apollo 14 breccias 14301, 14313, and 14318 all show fission xenon due to the decay of Pu-244. To investigate possible in situ production of the fission gas, an analysis was made of the U-distribution in these three breccias. The major amount of the U lies in glass clasts and in matrix material and no more than 25% occurs in distinct high-U minerals. The U-distribution of each breccia is discussed in detail. Whitlockite grains in breccias 14301 and 14318 found with the U-mapping were etched and analyzed for fission tracks. The excess track densities are much smaller than indicated by the Xe-excess. Because of a preirradiation history documented by very high track densities in feldspar grains, however, it is impossible to attribute the excess tracks to the decay of Pu-244. A modified track method has been developed for measuring average U-concentrations in samples containing a heterogeneous distribution of U in the form of small high-U minerals. The method is briefly discussed, and results for the rocks 14301, 14313, 14318, 68815, 15595, and the soil 64421 are given.
Nonlinear filtering for spacecraft attitude estimation
Nonlinear filtering techniques are applied to spacecraft attitude estimation using quaternion parameterization for the attitude kinematics. By replacing the angular velocity vector by the gyro output vector, a state dependent noise vector is introduced in the seven-dimensional system equations. The resulting conditional probability density function from the Ito differential rule is governed by the Fokker Planck partial differential equation which is approximated by the second order mean and covariance differential equations. In order to minimize computer loading, the covariance propagation is carried out in six-dimensional state space using a matrix transformation. The star tracker data is used to update the covariance matrix in the seven-dimensional space. The algorithm is simulated for an earth pointing spacecraft mission, using Monte Carlo samples of gyro and star measurements. The performance of the second order filter is compared with the extended Kalman Filter through several simulation runs and drift rates have been identified.
Fluctuations at the blue edge of saturated wind lines in IUE spectra of O-type stars
We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-Lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles and then derived from formal solution integration using this source function. The more appropriate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10 percent or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.
Advanced particulate fibrous composite for thermal control of re-entry vehicles
Combining the selectively reflective properties of particulate reflective coatings with the low density and high temperature capability of ceramic fibrous materials results in a novel composite material with useful optical properties. Numerical predictions of the radiative properties of zirconia and silica composites for a range of mixture fractions and wavelengths show the improvements attainable by dispersing spheres throughout a fiber matrix. The radiation conductance is controlled by the percentage of spheres in the composite. Results show that specific wavelength domains are controlled by either the spheres or the fibers, so the optical properties can be tailored for reentry applications.
On the synthesis of resonance lines in dynamical models of structured hot-star winds
We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles are then derived from formal solution integration using this source function. Two more approximate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10% or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.
Ceramic Aerogel Composite Materials and Characterization
Aerogels a.k.a "Solid Smoke" are gels with the liquid phase replaced by gas, leaving behind a highly porous material with a nanoscale framework. Due to the porous, nanoscale structure, aerogels have the lowest known density and conductivity of solids. Aerogels have the potential for being a breakthrough material because of their extremely light weight and unique properties. In this paper, we address overcoming their most profound weaknesses: mechanical fragility and very high surface activity, which leads to a lowered sintering temperature. A matrix of ceramic aerogel composite materials was produced to investigate their properties and functionality. Mechanical property measurements and Scanning Electron Micrographs are used to identify trends and structure of these ceramic composite materials. Thermal cycling was used to identify the sintering points of the materials.
Exploring the speciation of actinide salts in the presence of contaminants related to molten salt reactors
Molten salt reactors (MSR) are up-and-coming Generation IV nuclear reactors with either coolant and/or radioactive fuel in molten salt form. Due to their improved safety, efficiency, affordability and convenient waste processing system, these reactors are considered a superior alternative to conventional ones. However, prior to commercial implementation, it is necessary to have an extensive understanding of the reactions occurring in the MSR. The safety parameters of the MSR are determined by the rheological properties such as density and viscosity of the molten salt and these properties are governed by the local structure of the molten salts. Due to the high operational temperatures of MSR's, corrosion plays an important role. Because of their influence, the container corrosion products become part of the fuel salt matrix. In this work, we are investigating how these corrosion products, e.g., Mo, Ni, Cr, would impact the speciation of the actinide metal centers. For the initial studies, lanthanides are used as a surrogate for actinides. Lanthanide chloride salts (LnCl3) are mixed with transition metals and excess of alkali or alkali earth metal salts or salt mixtures (LiCl, NaCl, KCl). These are added to an alumina crucible, placed in a quartz tube, and sealed under a vacuum. This is heated up to around 1000 °C and slowly cooled to a temperature above the melting point of the salt and at that temperature, the reaction is taken out of the furnace. The excess salt is removed, and the resulting product is analyzed using X-ray diffraction techniques, scanning electron microscopy, and spectroscopic methods such as Raman, UV-Visible, and FTIR spectroscopy.
Constructing LDPC Codes from Loop-Free Encoding Modules
A method of constructing certain low-density parity-check (LDPC) codes by use of relatively simple loop-free coding modules has been developed. The subclasses of LDPC codes to which the method applies includes accumulate-repeat-accumulate (ARA) codes, accumulate-repeat-check-accumulate codes, and the codes described in Accumulate-Repeat-Accumulate-Accumulate Codes (NPO-41305), NASA Tech Briefs, Vol. 31, No. 9 (September 2007), page 90. All of the affected codes can be characterized as serial/parallel (hybrid) concatenations of such relatively simple modules as accumulators, repetition codes, differentiators, and punctured single-parity check codes. These are error-correcting codes suitable for use in a variety of wireless data-communication systems that include noisy channels. These codes can also be characterized as hybrid turbolike codes that have projected graph or protograph representations (for example see figure); these characteristics make it possible to design high-speed iterative decoders that utilize belief-propagation algorithms. The present method comprises two related submethods for constructing LDPC codes from simple loop-free modules with circulant permutations. The first submethod is an iterative encoding method based on the erasure-decoding algorithm. The computations required by this method are well organized because they involve a parity-check matrix having a block-circulant structure. The second submethod involves the use of block-circulant generator matrices. The encoders of this method are very similar to those of recursive convolutional codes. Some encoders according to this second submethod have been implemented in a small field-programmable gate array that operates at a speed of 100 megasymbols per second. By use of density evolution (a computational- simulation technique for analyzing performances of LDPC codes), it has been shown through some examples that as the block size goes to infinity, low iterative decoding thresholds close to channel capacity limits can be achieved for the codes of the type in question having low maximum variable node degrees. The decoding thresholds in these examples are lower than those of the best-known unstructured irregular LDPC codes constrained to have the same maximum node degrees. Furthermore, the present method enables the construction of codes of any desired rate with thresholds that stay uniformly close to their respective channel capacity thresholds.
Coupled Gel/Seeded Growth of Cs2AgBiBr6 Single Crystals for Radiation Detection
Abstract— The all-inorganic lead free Cs2AgBiBr6 perovskite represents a highly promising environment-friendly radiation detector material. The first attempt to grow Cs2AgBiBr6 single crystals using the unique gel growth method is reported in this study. This technique enables controlled reagent diffusion in a gel medium at room temperature, fostering the growth of visible single crystals. It offers advantages over conventional methods, leading to fewer equilibrium and non-equilibrium defects and simplified execution. Interestingly, the as-grown crystals exhibited variations in size and colors within the gel matrix, indicating the influence of Bi solutes in Bi-poor and Bi-rich Cs2AgBiBr6 crystals. Once the crystals were confirmed by powder xray diffraction, a series of comprehensive material characterizations was performed. The resistivity of Cs2AgBiBr6 single crystals was measured at 1.12x107 Ω•cm. Trap density analysis yielded a value of 2.82 x1010 cm-3. The charge carrier mobility was determined to be 32.9 cm2V-1S-1, suggesting the suitability for radiation detector applications. Additional growth optimization is going on and further characterization will be reported during the presentation. Our work helps accelerate the development of lead-free perovskite materials toward radiation detector applications.