Homogenization of radiation transport calculations for multi-physics simulations
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.
In this study, under NASA’s Rapid and Analysis Manufacturing Propulsion Technology (RAMPT) project, effects of different surface post-treatments on surface texture and fatigue behavior of L-PBF Alloy 214 were investigated. Various subtractive SPTs including abrasive flow machining (AFM), machining (M), shot peening (SP), vapormatt (VM), dry electropolishing (DE), chemical milling (CM), chemical-mechanical polishing (CMP), electro chemical (ECP) were applied to study the variations of the surface texture, microstructure and uniaxial fatigue behavior under fully reversed strain-controlled condition at four different strain amplitudes of 0.005, 0.004, 0.003, and 0.0025 mm/mm. The insights gained regarding the relationship between structural and mechanical properties for different SPTs can help establish guidelines for selecting the most effective process to enhance the performance of Laser Powder Bed Fusion (L-PBF) Alloy 214, especially for critical aerospace applications.
Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.
Achieving core-edge compatible divertor detachment is a critical requirement for stable operation in future fusion devices. This study compares nitrogen and neon seeding in KSTAR H-mode plasmas with carbon walls, combining experiments and SOLPS-ITER modelling to evaluate their radiative dissipation and core-edge compatibility. Experimentally, N seeding achieved stronger divertor detachment, with larger reductions in target particle and heat fluxes, a higher divertor radiation fraction, and a lower core radiation fraction than Ne. In contrast, Ne seeding triggered a significant rise in core radiation followed by H–L back transitions, limiting the maximum total radiated power fraction to roughly half that of N. SOLPS-ITER simulations reproduced the experimental trends and revealed that the better core-edge compatibility of N arises from its higher divertor retention in addition to its higher cooling factor. The relative positions of the stagnation points of impurity poloidal velocity and ionization sources did not explain the different divertor compression. Instead, in the present modelling, the higher impurity parallel particle flux, resulting from the higher parallel impurity velocity, explains the stronger nitrogen impurity compression in the divertor region. The parallel temperature distribution with N was more favourable for achieving higher impurity parallel velocity than with Ne, because the impurity velocity is governed by modifications of the main ion flow due to friction and thermal forces, both of which strongly depend on the temperature. Ultimately, this behaviour is attributed to the strongly divertor-localized radiation of N. These results demonstrate that N is more effective than Ne in achieving radiative divertor detachment while maintaining low core contamination in KSTAR, consistent with observations in other present tokamaks.
Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.
Combine the Dronebase aerial imagery with corresponding sites in the NLR Photovoltaic (PV) Fleets database. By combining these two data sources in an aggregated, anonymized fashion, we can perform the following analyses: quantifying power loss due to outages caused by stuck trackers, string outages, and shading/snow, validate site metadata, including tilt and azimuth, and correlate.
Explore the source record for details and available documents.
Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.
Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.
Explore the source record for details and available documents.
In response to the elevated-temperature and weight-reduction demands of modern aerospace applications, a novel oxide-dispersion-strengthened low-density niobium alloy (LDNb-ODS) was fabricated using laser powder bed fusion (L-PBF). To overcome powder procurement barriers, L-PBF feedstock was produced by blending commercial Nb521, Ti64, and Cr powder with Y2O3 nanoparticles via resonant acoustic mixing. Following L-PBF and a 1400°C vacuum heat treatment, the alloy achieved a density of 6.73 g/cc and a fine mean grain size of 4.62 µm stabilized by uniform ~30 nm yttria dispersoids. Microstructural analysis revealed a chemically inhomogeneous build with lack-of-fusion defects and a titanium (Ti) shift from a nominal 31.5 wt% in the starting powder blend to 24.8 wt% in the printed part due to preferential Ti loss during printing. Elevated-temperature tensile testing demonstrated that LDNb-ODS maintained a superior specific yield strength of 60-85 MPa/(g/cc) up to 800°C, outperforming nickel-based alloys Ni625, Ni230, and GRX-810. Between 870°C and 950°C, its specific strength surpassed both Ni718 and Ni625. In rapid stress-rupture testing at 1093°C and 20.7 MPa, uncoated LDNb-ODS survived 21.4 hours (a tenfold increase over legacy C-103) while an R512E silicide coating extended rupture life to 84.8 hours, confirming that oxidation accelerates low-stress failure. These findings demonstrate that additive manufacturing of LDNb-ODS provides a viable, lightweight alternative to nickel-based superalloys for high-temperature (>850°C) aerospace components.
Multipactor discharge is a persistent challenge in high-power microwave (HPM) and accelerator systems, where secondary electron avalanches can cause heating, vacuum degradation, and failure. This work presents the first supervised machine learning (ML) framework for multipactor prediction, trained on high-fidelity 3D Particle-in-Cell (PIC) simulation data in planar geometries. The model maps operational, geometric, and material-dependent secondary electron yield (SEY) parameters to the time-averaged electron growth rate, enabling rapid reconstruction of susceptibility charts. Among the models evaluated, tree-based ensemble methods such as Random Forest and Extra Trees demonstrate superior generalization to unseen materials compared to neural networks such as multilayer perceptron (MLP). Performance metrics, including Intersection over Union (IoU), Structural Similarity Index Measure (SSIM), and Pearson correlation, show close agreement with simulation benchmarks. Principal Component Analysis attributes generalization limits to material feature-space disjointedness.
We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.
Composite materials are increasingly being used in aerospace applications due to their superior strength-to-weight ratio compared to commonly used metals. A current limitation to widespread adoption is the certification of adhesively bonded joints. One approach to improving adhesive bonding in composites is accurately measuring the thickness of adhesive bondlines in composite laminates. Precise bondline thickness control is essential for aerospace applications where adhesive layer thickness directly affects joint fracture properties and structural performance. This study focused on implementing machine learning techniques to determine the ultrasonic time of flight (directly correlated to thickness) in adhesive bondlines throughout autoclave cure cycles. A high-temperature (use up to 180°C) ultrasonic scanning system was deployed in an autoclave to provide time of flight data through composite panels. Three experiments were conducted on the curing of 305 mm × 305 mm unidirectional composite panels. In the first experiment, a piecewise function was fit for the temperature correction factor to account for changing autoclave temperatures. Due to deficiencies in the first calibration experiment, a second experiment was run, and the results were used to train a machine learning model. The revised experiment, in combination with the machine learning model, significantly increased the accuracy of the bondline time of flight predictions (~14% error reduced to <1%). Data was processed using the Regression Learner Application in MATLAB®, with a Support Vector Machine selected for the model. The result was a machine learning algorithm capable of reliably quantifying ultrasonic time of flight through adhesive bondlines. The third experiment provided independent test data for the machine learning model, demonstrating that the model produces accurate predictions from data beyond its training set.
We show that trapped ions in virtual cathode potential wells can raise the transmitted current of emitted electrons into a plasma much closer to the full emission than is predicted by cathode sheath theories without trapped ions. The transmitted current is controlled by the well barrier voltage, which must adjust to balance the creation of low-energy ions within the well, and their loss. Our model considers the case of a plasma-facing cathode where trapped ions are created passively via charge-exchange collisions and lost passively via thermal leakage over the well. We quantify these rates and estimate the current in terms of system parameters for thermionic emission into a plasma with several cathode geometries. A general prediction is that the current as a function of emitted flux does not saturate at the traditional space charge limit (the onset of a well) but can reach far higher values until the trapped ion balance breaks down, causing instability. The maximum stable current depends on parameters but in principle can be arbitrarily high if active techniques are used to manipulate the trapped ion balance. We conclude that major improvements in plasma technologies with hot cathodes might be achieved by optimizing the current enhancement enabled by trapped ions.
Ultra-high-temperature ceramics (UHTCs), including TiC, Ti 11 B₂, and Zr 11 B₂, show great potential for plasma-facing components due to their excellent high-temperature properties prior to irradiation. However, their response to neutron irradiation remains insufficiently understood, limiting robust assessment of their viability for fusion energy applications. Here, this study examines the thermal conductivity, dimensional stability and microstructure of TiC, TiB₂, and ZrB₂ following neutron irradiation at temperatures of 200–1000 °C and fast neutron fluences of 2.0 × 10 25 to 1.1 × 10 26 n/m 2 (E > 0.1 MeV). Lattice swelling measured by synchrotron X-ray diffraction in all three UHTCs was maximized at 200 °C and decreased with increasing irradiation temperature, with no evidence of amorphization observed at 200 °C. Above 600 °C, significant macroscopic volume swelling was observed in irradiated Ti 11 B₂ and Zr 11 B₂, but not in TiC, likely due to cavity formation in the diborides. The post-irradiation thermal conductivity, measured at the irradiation temperature, ranged from 28 to 45 W/m·K, representing a 34–45% reduction relative to the unirradiated material. Notably, neutron-irradiated UHTCs exhibit recoverable thermal conductivity at elevated temperatures, comparable to ferritic–martensitic steels and potentially superior to W when transmutation effects are considered, highlighting promise for shielding or armor plasma-facing components. At 600 °C, both thermal conductivity degradation and lattice swelling saturated at doses exceeding 2–4 dpa.