Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Lamb Wave Stiffness Characterization of Composites Undergoing Thermal-Mechanical Aging

The introduction of new, advanced composite materials into aviation systems requires a thorough understanding of the long term effects of combined thermal and mechanical loading upon those materials. Analytical methods investigating the effects of intense thermal heating combined with mechanical loading have been investigated. The damage mechanisms and fatigue lives were dependent on test parameters as well as stress levels. Castelli, et al. identified matrix dominated failure modes for out-of-phase cycling and fiber dominated damage modes for in-phase cycling. In recent years, ultrasonic methods have been developed that can measure the mechanical stiffness of composites. To help evaluate the effect of aging, a suitably designed Lamb wave measurement system is being used to obtain bending and out-of-plane stiffness coefficients of composite laminates undergoing thermal-mechanical loading. The system works by exciting an antisymmetric Lamb wave and calculating the velocity at each frequency from the known transducer separation and the measured time-of-flight. The same peak in the waveforms received at various distances is used to measure the time difference between the signals. The velocity measurements are accurate and repeatable to within 1% resulting in reconstructed stiffness values repeatable to within 4%. Given the material density and plate thickness, the bending and out-of-plane shear stiffnesses are calculated from a reconstruction of the dispersion curve. A mechanical scanner is used to move the sensors over the surface to map the time-of-flight, velocity, or stiffnesses of the entire specimen. Access to only one side of the material is required and no immersion or couplants are required because the sensors are dry coupled to the surface of the plate. In this study, the elastic stiffnesses D(sub 11), D(sub 22), A(sub 44), and A(sub 55) as well as time-of-flight measurements for composite samples that have undergone combined thermal and mechanical aging for a duration of 10,000 hours are reported.

Seale, Michael D.↗
DOE OSTI2017

Investigation of Microbunching Instabilities in Modern Recirculating Accelerators

Particle accelerators are machines to accelerate and store charged particles, such as electrons or protons, to the energy levels for various scientific applications. A collection of charged particles usually forms a particle beam. There are three basic types of particle accelerators: linear accelerators (linac), storage-ring (or circular) accelerators, and recirculating accelerators. In a linac, particles are accelerated and pass through once along a linear or straight beamline. Storage-ring accelerators propel particles around a circular track and repetitively append the energy to the stored beam. The third type, also the most recent one in chronology, the recirculating accelerator, is designed to accelerate the particle beam in a short section of linac, circulate the beam, and then either continue to accelerate for energy boost or decelerate it for energy recovery. The beam properties of a linac machine are set at best by the initial particle sources. For storage rings, the beam equilibria are instead determined by the overall machine design. The modern recirculating machines share with linacs the advantages to both accelerate and preserve the beam with high beam quality, as well as efficiently reuse the accelerating components. The beamline design in such a machine configuration can however be much more complicated than that of linacs. As modern accelerators push toward the high-brightness or high-intensity frontier by demanding particles in a highly charged bunch (about nano-Coulomb per bunch) to concentrate in an ever-decreasing beam phase space (transverse normalized emittance about 1 μm and relative energy spread of the order of 10^-5 in GeV beam energy), the interaction amongst particles via their self-generated electromagnetic fields can potentially lead to coherent instabilities of the beam and thus pose significant challenges to the machine design and operation. In the past decade and a half, microbunching instability (MBI) has been one of the most challenging issues for such high-brightness or high-intensity beam transport, as it would degrade lasing performance in the fourth-generation light sources, reduce cooling efficiency in electron cooling facilities, and compromise the luminosity of colliding beams in lepton or lepton-hadron colliders. The dissertation work will focus on the MBI in modern recirculating electron accelerators. It has been known that the collective interactions, the coherent synchrotron radiation (CSR) and the longitudinal space charge (LSC) forces, can drive MBI. The CSR effect is a collective phenomenon in which the electrons in a curved motion, e.g. a bending dipole, emit radiation at a scale comparable to the micro-bunched structure of the bunch distribution. The LSC effect stems from non-uniformity of the charge distribution, acts as plasma oscillation, and can eventually accumulate an amount of energy modulation when the beam traverses a long section of a beamline. MBI can be seeded by non-uniformity or shot noise of the beam, which originates from granularity of the elementary charge. Through the aforementioned collective effects, the modulation of the bunch sub-structure can be amplified and, once the beam-wave interaction formed a positive feedback, can result in MBI. The problem of MBI has been intensively studied for linac-based facilities and for storage-ring accelerators. However, systematic studies for recirculation machines are still very limited and form a knowledge gap. Because of the much more complicated machine configuration of the recirculating accelerators than that of linacs, the existing MBI analysis needs to be extended to accommodate the high-brightness particle beam transport in modern recirculating accelerators. This dissertation is focused on theoretical investigation of MBI in such machine configuration in the following seven themes: (1) Development and generalization of MBI theory The theoretical formulation has been extended so as to be applicable to a general linear beamline lattice including horizontal and vertical transport bending elements, and beam acceleration or deceleration. These featured generalizations are required for MBI analysis in recirculation accelerators. (2) Construction of CSR impedance models In addition to the steady-state CSR interaction, it has been found that the exit transient effect (or CSR drift) can even result in more serious MBI in high-brightness recirculation arcs. The onedimensional free-space CSR impedances, especially the exit transients, are derived. The steady-state CSR impedance is also extended to non-ultrarelativistic beam energy for MBI analysis of low-energy merger sections in recirculating accelerators. (3) Numerical implementation of the derived semi-analytical formulation This includes the development of a semi-analytical Vlasov solver for MBI analysis, and also benchmarking of the solver against massive particle tracking simulations. (4) Exploration of multistage amplification behavior of CSR microbunching development The CSR-induced MBI acts as an amplifier, which amplifies the sub-bunch modulation of a beam. The amplification is commonly quantified by the amplification gain. A beam transport system can be considered as a cascaded amplifier. Unlike the two-stage amplification of four-dipole bunch compressor chicanes employed in linacs, the recirculation arcs, which are usually constituted by several tens of bending magnets, show a distinguishing feature of up to six-stage microbunching amplification for our example arc lattices. That is, the maximal CSR amplification gain can be proportional to the peak bunch current up to sixth power. A method to compare lattice performance has been developed in terms of gain coefficients, which nearly depend on the lattice properties only. This method has also proven to be an effective way to quantify the current dependence of the maximal (5) Control of CSR MBI in multibend transport or recirculation arcs The existing mitigation schemes of MBI mostly aim to linac-based accelerators and may not be practical to the recirculating accelerator facilities. Thus a set of conditions for suppression of CSR MBI was proposed and examined for example lattices from low (~100 MeV) to high (~1 GeV) energies. (6) Study of more aspects of microbunched structures in beam phase spaces For a cascaded amplifier in circuit electronics, the total amplification gain can be estimated as the product of individual gains. In a beam transport line of an accelerator, the (scalar) gain multiplication was examined and found to under-estimate the overall microbunching amplification. The concept of gain matrix was developed, which includes the density, energy and transverse-longitudinal modulations in a beam phase space, and used to analyze MBI for a proposed recirculating machine. Throughout the gain matrix approach, it reasonably gives the upper limit of spectral MBI gain curves. This extended analysis can be employed to study multi-pass recirculation. (7) Study of MBI for magnetized beams Driven by a recent energy-recovery-linac based cooler design for electron cooling at Jefferson Lab Electron-Ion Collider Project, the generalized theoretical formulation for MBI to a transversely coupled beam has been developed and applied to this study. A magnetized beam in general features non-zero canonical angular momentum, thus considered to be a transversely coupled beam. A novel idea of utilizing magnetized beam transport was proposed for improvement of cooling efficiency and possible mitigation of collective effects. A concern of MBI regarding this design was studied and excluded. The large transverse beam size associated with the beam magnetization is found to help suppress MBI via the transverse-longitudinal correlation.

43 PARTICLE ACCELERATORS↗
NASA NTRS1996

Hypersonic Airbreathing Vehicles/Technologies

Hypersonic airbreathing horizontal takeoff and landing (HTOL) vehicles are highly integrated systems involving many advanced technologies. The design environment is variable rich, intricately networked, and sensitivity intensive; as such, it represents a tremendous challenge. Creating a viable design requires addressing three main elements: (1) an understanding of the 'figures of merit' and their relationship, (2) the development of sophisticated configuration discipline prediction methods and a synthesis procedure, and (3) the synergistic integration of advanced technologies across the discipline spectrum. This paper will focus on the vision for hypersonic airbreathing vehicles and the advanced technologies that forge the designs. Airbreathing hypersonics encompass endoatmospheric (airplanes...missiles are a part of the matrix but will not be included in this paper since they are an air force focus) and space access vehicles with speed from Mach 4 up to Mach 25 (orbital). These vehicles can be divided into two classes...cruisers and accelerators. The cruiser designs reflect high lift-to-drag whereas the accelerators reflect low drag per unit inlet capture; thus, the cross section of the accelerator attributes a much larger percentage to propulsion. One of the more design influencing items is fuel. The hydrogen fueled vehicles must be very volumetric efficient to contain the low density fuel and thus tend to be a bit bulgy (more conducive to lifting bodies or wing bodies) whereas with hydrocarbon fueled vehicles, the concern is loading because of the high density fuel; thus, they may tend to be more towards waveriders which are not usually very volumetric efficient. Hydrocarbon fuels (endothermic) are limited in their engine cooling capacity to below Mach 8.

Hunt, James L.↗
NASA NTRS1999

Method of Fabricating Chopped-Fiber Composite Piston

A three-dimensional piston molding is fabricated from a mixture of chopped, carbon tow filaments of variable length, which are prepregged with carbonaceous organic resins and/or pitches and molded by conventional molding processes into a near net shape, to form a carbon-fiber reinforced organic-matrix composite part. Continuous reinforcement in the form of carbon-carbon composite tapes or pieces of fabric can be also laid in the mold before or during the charging of the mold with the chopped-fiber mixture, to enhance the strength in the crown and wrist-pin areas. The molded chopped-fiber reinforced organic-matrix composite parts are then pyrolized in an inert atmosphere, to convert the organic matrix materials to carbon. These pyrolized parts are then densified by reimpregnation with resins or pitches, which are subsequently carbonized. Densification is also accomplished by direct infiltration with carbon by vapor deposition processes. Once the desired density has been achieved, the piston molds are machined to final piston dimensions, and piston ring grooves are added. To prevent oxidation and/or to seal the piston surface or near surface, the chopped-fiber piston is coated with ceramic and/or metallic sealants: and/or coated with a catalyst.

Rivers, H. Kevin↗
DOE OSTI2025

Electrolytic gold plating, stripping, and ion transport dynamics through a solid-state iodide perovskite

The pronounced electrochemical reactivity between halide perovskites and metal electrodes can introduce mobile extrinsic metal ions which can cause device instability or enable novel functionalities. Here we systematically investigate the kinetics of gold cation (Au + ) migration in indium tin oxide (ITO)/methylammonium lead triiodide (MAPbI 3 )/Au model devices under long-term potentiostatic biasing. Scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) analyses reveal that Au + ions, electrochemically generated at the Au anode, traverse the perovskite layer with diffusion coefficients on the order of 10 −11 to 10 −10 cm 2 s −1 and are subsequently reduced at the cathode as Au 0 clusters, resembling metal plating behavior in electrolytic cells and solid-state batteries during charging. Furthermore, reversing the applied bias strips the plated Au 0 , revealing reversibility suitable for bipolar resistive switching devices and providing direct evidence of the electrochemical and ionic nature of Au transport within the perovskite matrix. Quantitatively determining diffusion coefficients and ion concentrations provides foundational inputs for future drift-diffusion modelling opportunities and allows us to relate our findings to implications on long term operation of devices like photovoltaic modules. These results clearly demonstrate the solid-state electrochemical nature of perovskite devices, highlight methods to be more quantitative about ion transport properties, provide and emphasize the importance of disentangling electro-, photo-, photoelectrochemical processes for understanding device performance and unlocking new functionalities.

14 SOLAR ENERGY↗
DOE OSTI2024

In-Plane Liftout and Push-to-Pull for In Situ Mechanical Testing of Irradiated Inconel X-750

A streamlined sample preparation method for nanomechanical testing is needed to improve the quality of specimens, reduce the cost, and increase the versatility of specimen fabrication. This work outlines an in-plane liftout focused ion beam (FIB) fabrication procedure to prepare electron-transparent specimens for in situ transmission electron microscopy (TEM) nanomechanical testing. Ion etching and electron backscatter diffraction (EBSD) techniques were used to lift out a [110] oriented grain from a neutron-irradiated bulk X-750 alloy. Careful control of voltages and currents ensured precision. Top surface thinning sweeps prevented resurfacing and redeposition while dog-bone geometries were shaped with a 1:4 gauge width-to-milling pattern diameter ratio. Nanotensile testing in the TEM with a picoindenter allowed for the estimation of an ultimate tensile strength of 2.41 GPa, and inspection revealed a high density of bubbles in the X-750 matrix. The proposed fabrication procedure is significant for preparing samples from radioactive materials, studying complex structures that are orientation-dependent, and analyzing desired planar areas.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Thermal Neutron Scattering Evaluation Methods Development at ORNL [Slides]

The goal of the project was to develop a procedure to generate covariance matrix for S( α,β ) data that incorporated both computational simulations and experimental data. Experimental data and computer simulation fit is achieved using the Unified Monte Carlo (UMC) method. The framework is material & simulation method independent, and this is demonstrated using light water. A methodology was developed for generalized generation of covariances for thermal scattering files. The framework was subsequently tested on light water by generating covariance of pDOS with encouraging results. Preliminary results suggest temperature effects of polyethylene may warrant further investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Isobaric Reconstruction of the Baryonic Acoustic Oscillation

In this Letter, we report a significant recovery of the linear baryonic acoustic oscillation (BAO) signature by applying the isobaric reconstruction algorithm to the nonlinear matter density field. Assuming only the longitudinal component of the displacement being cosmologically relevant, this algorithm iteratively solves the coordinate transform between the Lagrangian and Eulerian frames without requiring any specific knowledge of the dynamics. For dark matter field, it produces the nonlinear displacement potential with very high fidelity. The reconstruction error at the pixel level is within a few percent and is caused only by the emergence of the transverse component after the shell-crossing. As it circumvents the strongest nonlinearity of the density evolution, the reconstructed field is well described by linear theory and immune from the bulk-flow smearing of the BAO signature. Therefore, this algorithm could significantly improve the measurement accuracy of the sound horizon scale s. For a perfect large-scale structure survey at redshift zero without Poisson or instrumental noise, the fractional error Δs/s is reduced by a factor of ∼2.7, very close to the ideal limit with the linear power spectrum and Gaussian covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Time-dependent density functional theory for open systems with a positivity-preserving decomposition scheme for environment spectral functions

Understanding electronic dynamics on material surfaces is fundamentally important for applications including nanoelectronics, inhomogeneous catalysis, and photovoltaics. Practical approaches based on time-dependent density functional theory for open systems have been developed to characterize the dissipative dynamics of electrons in bulk materials. The accuracy and reliability of such approaches depend critically on how the electronic structure and memory effects of surrounding material environment are accounted for. In this work, we develop a novel squared-Lorentzian decomposition scheme, which preserves the positive semi-definiteness of the environment spectral matrix. The resulting electronic dynamics is guaranteed to be both accurate and convergent even in the long-time limit. The long-time stability of electronic dynamics simulation is thus greatly improved within the current decomposition scheme. The validity and usefulness of our new approach are exemplified via two prototypical model systems: quasi-one-dimensional atomic chains and two-dimensional bilayer graphene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Quantitative characterization of spatial distribution of particles in materials: Application to materials processing

Most engineering materials contain second phase particles or fibers which serve to reinforce the matrix phase. The effect of reinforcements on material properties is usually analyzed in terms of the average volume fraction and spacing of reinforcements, quantities which are global microstructural characteristics. However, material properties can also depend on local microstructural characteristics; for example, on how uniformly the reinforcing phase is distributed in the material. The analysis method will then be applied to a materials processing problem to discover how processing parameters can be selected to maximize redistribution of the reinforcing phase during processing. Several mathematical analysis methods could be adapted to the problem of characterizing the distribution of particles in materials. A tessellation-based method was selected. In the first phase of the investigation, a software package was written to automate the analysis. Typical results are shown. The analysis technique allows the degree to which particles are clustered together, the size and spacing of particle clusters, and the particle density in clusters to be found. The analysis methods were applied to computer-generated distributions and to a few real particle-containing materials. Methods for analyzing a nonuniform particle distribution in a material can be applied to two broad classes of materials science problems: understanding how the resulting particle distribution affects properties. The analysis method described is applied to a materials processing problem: how to select extrusion conditions to maximize the redistribution of reinforcing particles that are initially nonuniformly distributed. In addition, the tessellation-based method to analyze star distributions in spiral galaxies was adapted, illustrating the diverse types of problems to which the analysis method can be applied.

Parse, J. B.↗
DOE OSTI2011

Comparison of the Effects of Carbon Ion and Photon Irradiation on the Angiogenic Response in Human Lung Adenocarcinoma Cells

Purpose: Radiotherapy resistance is a commonly encountered problem in cancer treatment. In this regard, stabilization of endothelial cells and release of angiogenic factors by cancer cells contribute to this problem. In this study, we used human lung adenocarcinoma (A549) cells to compare the effects of carbon ion and X-ray irradiation on the cells' angiogenic response. Methods and Materials: A549 cells were irradiated with biologically equivalent doses for cell survival of either carbon ions (linear energy transfer, 170 keV/{mu}m; energy of 9.8 MeV/u on target) or X-rays and injected with basement membrane matrix into BALB/c nu/nu mice to generate a plug, allowing quantification of angiogenesis by blood vessel enumeration. The expression of angiogenic factors (VEGF, PlGF, SDF-1, and SCF) was assessed at the mRNA and secreted protein levels by using real-time reverse transcription-PCR and enzyme-linked immunosorbent assay. Signal transduction mediated by stem cell factor (SCF) was assessed by phosphorylation of its receptor c-Kit. For inhibition of SCF/c-Kit signaling, a specific SCF/c-Kit inhibitor (ISCK03) was used. Results: Irradiation of A549 cells with X-rays (6 Gy) but not carbon ions (2 Gy) resulted in a significant increase in blood vessel density (control, 20.71 {+-} 1.55; X-ray, 36.44 {+-} 3.44; carbon ion, 16.33 {+-} 1.03; number per microscopic field). Concordantly, irradiation with X-rays but not with carbon ions increased the expression of SCF and subsequently caused phosphorylation of c-Kit in endothelial cells. ISCK03 treatment of A549 cells irradiated with X-rays (6 Gy) resulted in a significant decrease in blood vessel density (X-ray, 36.44 {+-} 3.44; X-ray and ISCK03, 4.33 {+-} 0.71; number of microscopic field). These data indicate that irradiation of A549 cells with X-rays but not with carbon ions promotes angiogenesis. Conclusions: The present study provides evidence that SCF is an X-ray-induced mediator of angiogenesis in A549 cells, a phenomenon that could not be observed with carbon ion irradiation. Thus, in this model system evaluating angiogenesis, carbon ion irradiation may have a therapeutic advantage. This observation should be confirmed in orthotopic lung tumor models.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Add-a-Source Analysis Method for Neutron Coincidence Counting of Waste: Laboratory Exercises

In this session you will setup the INCC code to perform calibration and verification measurements on simulated waste drums using the Waste Drum Assay System (WDAS) detector, which uses the add-a-source (AAS) method. The basis of the AAS method is to measure the change in the counting rate from a small 252 Cf source (~3x10 4 n/s) located outside of the drum and to use this information to correct the count rate from material inside the drum for any changes caused by the matrix inside the drum. The count rate of the material inside the drum is corrected to the empty drum case so that only a single mass calibration curve has to be measured (for the empty drum configuration). In the WDAS system, the AAS neutron source has one location, centered beneath the waste drum. The sample matrix has two primary effects on the neutrons: (1) energy reduction by scattering reactions and (2) neutron absorption of the low-energy neutrons. In general, the counter is designed with an optimum moderator (CH 2 ) thickness to be relatively insensitive to the energy reduction; however, as the hydrogen density in the drum increases, the adsorption process significantly reduces the measured neutron signal.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Machine learning assisted prediction of tungsten heavy alloy plasma facing component performance for fusion energy applications

Tungsten and tungsten heavy alloys (WHAs), known for their remarkably high hardness, durability, and corrosion resistance, play a critical role in the thriving development of nuclear fusion reactors in recent years. However, the exploration in tungsten alloys for the nuclear-related applications has been limited by the difficulty of manufacturing and the complexity of experiments to reproduce the environment of nuclear reaction. Therefore, this project aims to utilize nanoscale simulation methods such as density functional theory (DFT) and molecular dynamics (MD) with the help of machine learning techniques to not only understand the mechanisms of tungsten alloys but also allow us to computationally predict their mechanical behaviors under extreme environments. One critical problem of the application of WHAs in nuclear reactors is the surface melting. In the current design of the SPARC reactor, the WHA, W97Ni2.1Fe0.9 or W97NiFe, is chosen to be the first wall components to confine the plasma where the particles are fiercely moving and colliding into each other to create nuclear fusion reaction. This process will generate extremely high heat flux onto these WHA tiles, leaving high surface temperature that could possibly melt the surface of the WHA tiles, As illustrated in Fig. 1(a). a laser experiment previously done illustrates that a rough surface damage would be made after the surface melting where the matrix area mainly composed of nickel and iron as shown in Fig. 1(b), will first melt and then leave vacancies between these tungsten grains. Unfortunately, these kinds of roughness on the first-wall components could be deadly to the plasma inside a Tokmak reactor because the heat that is supposed to dissipate at a designed ratio through the tiles may in turn be excessively absorbed and accumulated on any uneven area of the surface, which will eventually make the whole nuclear reaction fail. In this project, we will introduce a machine learning potential, Allegro, based on DFT calculation and then build a MD model for W-Ni-Fe alloys.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Metal Matrix Composite LOX Turbopump Housing Via Novel Tool-less Net-Shape Pressure Infiltration Casting Technology

Metal matrix composites (MMC) offer relatively higher specific strength, specific stiffness, lower coefficient of thermal expansion (CTE) and lower density as compared with conventional alloys. These unique properties make them very attractive for aerospace turbomachinery applications where there is ever increasing emphasis to reduce weight and cost, and to increase engine performance. Through a joint effort between NASA and Metal Matrix Cast Composites, Inc., a complex liquid oxygen (LOX) compatible turbopump housing is being redesigned and manufactured from hybrid (particulate and fibers) Aluminum MMC. To this end, a revolutionary toolless pressure infiltration casting technology is being perfected. Ceramic preforms for the composite are 3-dimensionally printed using a stereolithography file, acquired from a CAD model. The preforms are then invested into a refractory material and pressure infiltrated with liquid metal. After casting, the refractory material is washed away leaving behind a near net-shape composite part. Benefits of this process include increased composite uniformity, no mold machining, short time from design to part, properties matching traditional methods, ability to make previously impossible to manufacture parts and no size limitations with a newly developed joining technology. The results of materials, manufacturing and design optimizations, preform joining, and sub-element tests will be presented.

Shah, Sandeep↗
DOE OSTI2022

Regularized inversion of aerosol hygroscopic growth factor probability density function: application to humidity-controlled fast integrated mobility spectrometer measurements

Abstract. Aerosol hygroscopic growth plays an important role in atmospheric particle chemistry and the effects of aerosol on radiation and hence climate. The hygroscopic growth is often characterized by a growth factor probability density function (GF-PDF), where the growth factor is defined as the ratio of the particle size at a specified relative humidity to its dry size. Parametric, least-squares methods are the most widely used algorithms for inverting the GF-PDF from measurements of the humidified tandem differential mobility analyzer (HTDMA) and have been recently applied to the GF-PDF inversion from measurements of the humidity-controlled fast integrated mobility spectrometer (HFIMS). However, these least-squares methods suffer from noise amplification due to the lack of regularization in solving the ill-posed problem, resulting in significant fluctuations in the retrieved GF-PDF and even occasional failures of convergence. In this study, we introduce nonparametric, regularized methods to invert the aerosol GF-PDF and apply them to HFIMS measurements. Based on the HFIMS kernel function, the forward convolution is transformed into a matrix-based form, which facilitates the application of the nonparametric inversion methods with regularizations, including Tikhonov regularization and Twomey's iterative regularization. Inversions of the GF-PDF using the nonparameteric methods with regularization are demonstrated using HFIMS measurements simulated from representative GF-PDFs of ambient aerosols. The characteristics of reconstructed GF-PDFs resulting from different inversion methods, including previously developed least-squares methods, are quantitatively compared. The result shows that Twomey's method generally outperforms other inversion methods. The capabilities of Twomey's method in reconstructing the pre-defined GF-PDFs and recovering the mode parameters are validated.

54 ENVIRONMENTAL SCIENCES↗