Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transference number”

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 577 records · Page 32

Iron Single‐Atom Catalyzed N‐Alkylation of Amines with Alcohols via Solvent‐Free Borrowing Hydrogen Strategy

Abstract Industrial hydrogenation is a pivotal process in chemical synthesis. However, it has significant drawbacks, including high cost, safety risks associated with the use of molecular hydrogen gas, and substantial energy demands due to the need for elevated temperatures and pressures to achieve satisfactory yields. The borrowing hydrogen synthesis, which enables the transfer of hydrogen between molecules, offers a promising approach for green, one‐pot synthesis of industrially important chemicals and intermediates. Despite its potential, the broad application remains limited due to the reliance on toxic solvents, expensive noble metal catalysts, and the still restricted efficiency and substrate scope. In this study, the first solvent‐free strategy for theN‐alkylation of amines with alcohols is presented, employing an N‐doped graphene‐supported Fe single‐atom catalyst (Fe SA @N‐G; 1.06 wt.%). This approach achieves superior conversion and selectivity (up to 99%) along with record values for turnover number (TON, 1032.7) and turnover frequency (TOF, 413.1 h -1 ) for the coupling reaction of aniline with benzyl alcohol, surpassing all previously reported catalysts. DFT calculations, combined with experimental data, elucidated the reaction mechanism and identified the Fe 1 (III)‐N 4 active site participating in Fe‐H hydride transfer and containing two pyrrolic and two pyridinic nitrogens bound to the Fe center. The developed technology is further supported by the catalyst's excellent scalability, reusability, and performance under continuous‐flow conditions. Additionally, the exceptional efficiency of the single‐atom catalyst is demonstrated across more than 50 substrates, including reactions involving both aliphatic and aromatic amines with aliphatic and aromatic alcohols. The industrial applicability of this technology is validated through the synthesis of pharmaceutically relevant compounds, including stimulant drugs, antihistamines, and pharmaceutical intermediates.

Chemistry↗

A novel transfer learning framework for sorghum biomass prediction using UAV-based remote sensing data and genetic markers

Yield for biofuel crops is measured in terms of biomass, so measurements throughout the growing season are crucial in breeding programs, yet traditionally time- and labor-consuming since they involve destructive sampling. Modern remote sensing platforms, such as unmanned aerial vehicles (UAVs), can carry multiple sensors and collect numerous phenotypic traits with efficient, non-invasive field surveys. However, modeling the complex relationships between the observed phenotypic traits and biomass remains a challenging task, as the ground reference data are very limited for each genotype in the breeding experiment. In this study, a Long Short-Term Memory (LSTM) based Recurrent Neural Network (RNN) model is proposed for sorghum biomass prediction. The architecture is designed to exploit the time series remote sensing and weather data, as well as static genotypic information. As a large number of features have been derived from the remote sensing data, feature importance analysis is conducted to identify and remove redundant features. A strategy to extract representative information from high-dimensional genetic markers is proposed. To enhance generalization and minimize the need for ground reference data, transfer learning strategies are proposed for selecting the most informative training samples from the target domain. Consequently, a pre-trained model can be refined with limited training samples. Field experiments were conducted over a sorghum breeding trial planted in multiple years with more than 600 testcross hybrids. The results show that the proposed LSTM-based RNN model can achieve high accuracies for single year prediction. Further, with the proposed transfer learning strategies, a pre-trained model can be refined with limited training samples from the target domain and predict biomass with an accuracy comparable to that from a trained-from-scratch model for both multiple experiments within a given year and across multiple years.

36 MATERIALS SCIENCE↗

High-Energy X-ray Tomographic Analysis of Precursor Metal Powders (Ti-6Al-4V) Used for Additive Manufacturing

We utilized high-energy x-ray tomography to characterize Ti-6Al-4V metal powders (both as-received and recycled) used in powder-bed additive manufacturing process. The image processing workflow was developed to process and analyze large amount of data objectively by computer program. The distribution of size and shape of the metal particles as well as defect (mainly porosity) inside the particles was analyzed with the statistical representation and resolution in micrometer. The result revealed that circular-shaped porosity with various sizes could be embedded in the powder particles. These porosities could potentially be transferred to the 3D printed part and critically affect the mechanical performance of the component. The present study shows the effectiveness of characterizing metal powders using x-ray imaging techniques where sufficient number of particles can be sampled within tens of minutes with a minimum sample preparation and high accuracy. Clear structural differences in the as-received and recycled powders were delineated that helps in determining the feasibility of using the recycled powders.

36 MATERIALS SCIENCE↗

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗

Comparative Analysis of Direct and Indirect Cooling of Wide-Bandgap Power Modules and Performance Enhancement of Jet Impingement-Based Direct Substrate Cooling

With the development of high-power and high-torque machines, requirements for high-power density electronics are increasing. Thermal management of such systems requires high heat extraction. Conventional air cooling based heat sinks and cold plate based liquid cooling have their own benefits for various applications but has limitations for high power density applications. The current study explores a jet impingement based direct substrate cooling system that was implemented for a SiC based direct bonded Cu substrate for various power losses. Numerical comparison between jet impingement cooling and conventional horizontal/indirect cooling (pin fin heat sink and genetic algorithm-optimized heat sink) showed that the area weighted average of the heat transfer coefficient (HTC) is high for both horizontal cooling designs, and the local HTC is higher for jet impingement. Design iterations were undertaken to resolve the bottleneck of this cooling system. Increasing the number of nozzles helped to cover more area at the direct bonded Cu bottom plate, which drops the chip temperature considerably. With a constant flow rate, increasing the number of nozzles would decrease local jet velocity, which reduces the heat extraction by jet impingement. This issue can be addressed by reducing the diameter of nozzle but doing so results in a high pressure drop where the design constraint is 2 psi. A flared nozzle design is proposed, which has a higher spreading angle of the jet that increases the flow coverage and reduces the pressure drop of the coolant loop.

Barua, Himel↗

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗

OrthoPhyl—streamlining large-scale, orthology-based phylogenomic studies of bacteria at broad evolutionary scales

Abstract There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis.

59 BASIC BIOLOGICAL SCIENCES↗

Optimization-Inspired Pin-Fin Array for Supercritical Carbon Dioxide Recuperator

Additively manufactured heat exchangers are one possible route to cost-effective sCO 2 power cycles. In this paper, experimental results are obtained for two helical pin fin tubes that were designed following parametric optimization of the fin array. Neither the numerical optimization nor the experimental testing have been previously reported in the literature. The two tube designs were, (1) the optimization inspired design, and (2) the optimization-inspired design with a fin diameter increased by a factor of 2. To characterize the print, the designs were scanned using X-ray computed tomography to measure feature sizes and heat transfer area. The optimization was conducted in a commercial, computational fluid dynamics code. The code solved the Reynolds Averaged Navier-Stokes (RANS) and energy equations with turbulence closure provided by the shear stress transport (SST) k-ω model. In the experiments, the Nusselt number augmentation was measured using the Wilson plot technique and the friction factor was determined with mass flow and pressure drop measurements. The experimental testing indicated that the optimization-inspired design had a friction factor that was four times less than the baseline tube design at equal Nusselt number. Additionally, the optimization-inspired design had a 14% improvement in Nusselt number, at equal friction factor, relative to the best performing computational fluid dynamics (CFD) trial points. Tube design (2), with the larger diameter pin fins, had similar performance, within experimental error, as the tube with the smaller diameter pin fins (1). Both tubes achieved overall fin array efficiencies near 1. A performance factor, V/V 0 , equal to the volume of the enhanced heat exchanger divided by the volume of the baseline (no-fins) heat exchanger, is recommended to quantify internal cooling performance. In conclusion, the experimental shell and tube heat exchanger, using the additively manufactured tube, is competitive with printed circuit heat exchangers in its pressure drop class and could be further improved by optimizing a shell-and-tube heat exchanger utilizing this heat transfer enhancement feature.

42 ENGINEERING↗

Dark Energy Survey Year 3 Results: Measuring the Survey Transfer Function with Balrog

Abstract We describe an updated calibration and diagnostic framework, Balrog , used to directly sample the selection and photometric biases of the Dark Energy Survey (DES) Year 3 (Y3) data set. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant “source” galaxies and magnification biases of nearer “lens” galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ∼1–8 mmag, but for a small subset of objects, we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. We discuss approaches to extend the current methodology to capture more aspects of the transfer function and reach full coverage of the survey footprint for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Leveraging History to Predict Infrequent Abnormal Transfers in Distributed Workflows

Scientific computing heavily relies on data shared by the community, especially in distributed data-intensive applications. This research focuses on predicting slow connections that create bottlenecks in distributed workflows. In this study, we analyze network traffic logs collected between January 2021 and August 2022 at the National Energy Research Scientific Computing Center (NERSC). Based on the observed patterns, we define a set of features primarily based on history for identifying low-performing data transfers. Typically, there are far fewer slow connections on well-maintained networks, which creates difficulty in learning to identify these abnormally slow connections from the normal ones. We devise several stratified sampling techniques to address the class-imbalance challenge and study how they affect the machine learning approaches. Our tests show that a relatively simple technique that undersamples the normal cases to balance the number of samples in two classes (normal and slow) is very effective for model training. This model predicts slow connections with an F1 score of 0.926.

97 MATHEMATICS AND COMPUTING↗

Accelerating Scientific Simulations with Bi-Fidelity Weighted Transfer Learning

High-fidelity modeling is an essential design tool for many engineering applications. However, for complex systems, computational cost can be a limiting factor. Analyzing parameter sensitivity, uncertainty quantification, and design optimization require many model evaluations. Surrogate models are often used to develop the relationship between model parameters and quantities of interest. However, in the case of complex systems, surrogate models require several degrees of freedom and, thus, a large number of data points to determine the correct dependencies. For many applications, this may be prohibitively expensive. The reduction of computational requirements can be achieved by leveraging low-fidelity models. Low-fidelity models represent the system at a coarser resolution with the advantage of computational efficiency. Therefore, a bi-fidelity modeling paradigm, which augments the accuracy of a low-fidelity model in a computationally efficient manner by invoking limited runs of a high-fidelity model, can be leveraged to sufficiently balance the accuracy and computational requirements. In this work, a bi-fidelity weighted transfer learning method using neural networks was applied to a computational fluid dynamics heat transfer modeling problem. The transfer learning advantage was investigated as a function of hyperparameters. Our main finding is that the use of a bi-fidelity modeling paradigm achieves accuracy close to that of a high-fidelity Gaussian process model while significantly reducing computational cost. The bi-fidelity model achieves comparable performance with 90 high-fidelity samples-that is, 60% less than the samples needed to achieve similar accuracy without the use of bi-fidelity modeling,

Borowiec, Katarzyna↗

No Galaxy Left Behind: Measuring the Transfer Function of the Dark Energy Survey with Balrog

In this dissertation, we describe a calibration and diagnostic framework called Balrog which was used to directly sample the selection and photometric biases of the Dark Energy Survey (DES) Year 3 (Y3) dataset. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with traditional generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant "source" galaxies and magnification biases of nearer "lens" galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the fiducial Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ~1-8 millimagnitudes, but for a small subset of objects we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. Finally, we discuss approaches to extend the current methodology to capture more aspects of the transfer function for future analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Delocalized polaron and Burstein-Moss shift induced by Li in α–V 2 O 5 : A DFT + DMFT study

For this work, we performed density functional theory (DFT)+U and dynamical mean field theory (DMFT) calculations with a continuous-time quantum Monte Carlo impurity solver to investigate the electronic properties of V 2 O 5 and Li x V 2 O 5 (x = 0.125 and 0.25). Pristine V 2 O 5 is a charge-transfer insulator with strong O p–V d hybridization, and it exhibits a large band gap (E gap ) as well as a nonzero conduction-band (CB) gap.We show that the band gap, the number of d electrons of vanadium, N d , and the CB gap for V 2 O 5 obtained from our DMFT calculations are in excellent agreement with the experimental values. While the DFT +U approach replicates the experimental band gap, it overestimates the value of N d and underestimates the CB gap. In the presence of low Li doping, the electronic properties of V 2 O 5 are mainly driven by a polaronic mechanism, and electron spin resonance and electron nuclear double resonance spectroscopies observed the coexistence of free and bound polarons. Notably, our DMFT results identify both polaron types, with the bound polaron being energetically preferred, while the DFT +U method only predicts the free polaron. Our DMFT analysis also reveals that increased Li doping leads to electron filling in the conduction band, shifting the Fermi level. This result is consistent with the observed Burstein-Moss shift upon enhanced Li doping, and we thus demonstrate that the DFT + DMFT approach can be used for an accurate and realistic description of strongly correlated materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Review and Analysis of Heat Transfer Correlations For Horizontal Pseudocritical CO 2 Heat Exchanger Applications

There is significant interest in the development of supercritical carbon dioxide (sCO 2 ) power cycles because of the potential for smaller and more energy efficient systems than a steam Rankine cycle. Heat exchanger designs typically use empirical correlations, but the applicability of these correlations near the CO 2 critical point is a potential issue. Though numerous correlations have been proposed in the literature, there are some disagreements when it comes to the accuracy. The current work recognizes the role of thermophysical properties, and its impacts on the heat transfer correlations and cycle efficiency. Heat transfer correlations proposed for horizontal flow inside circular pipes were analyzed with the help of numerical simulations. Steady state RANS simulations were performed using SST k-w turbulence model to evaluate the Nusselt number empirical correlations. It was found that the most of correlations (except Yoon) produced a Nusselt number that differed significantly with the one predicted numerically. Some of the correlations were developed for pure forced convection regime and as mentioned in Lin et al. do not account for mixed convection or free convection effects. Based on the limited observation, it appears that Yoon et al. predictions match well with the numerically predicted Nusselt Numbers. However, further analysis is required understand the applicability of various correlations and is contingent on accurate measurements or predictions of wall temperature profiles in the axial and the circumferential directions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Cell-free bioelectrocatalytic platform for carbon dioxide reduction (Final Technical Report)

The University of Minnesota (UMN) EcoSynBio Team aimed to develop a cell-free, enzyme-based platform for the electro- biocatalytic conversion of CO2 into formate as a platform chemical for further upgrading. This type of bio electrocatalytic process delivers a clean product stream without the need for extensive separation from the electrolyte as in electrochemical synthesis and microbial processes. The reduction reaction is catalyzed by metal-dependent formate dehydrogenases (mFDHs) that are capable of efficient electrocatalytic CO2 reduction without the need of costly co-factors. The development of an efficient, scalable electrobiocatalytic process with high total turnover numbers and viable space time yields, however, was not without its challenges. The UM team has developed a protein-based scaffolding system that facilitates enzyme stabilization and attachment to electrodes along with electron transfer. Yet, although FDHs are highly promising enzymes for cell-free, electrobiochemical CO2 reduction, they are also greatly understudied and especially for applications in electrocatalysis. The UM team used the best described mFDH from Clostridium as its benchmark system and spent significant time and effort in attempting to replicate published data and finally, redesigned a recombinant production system for proper metal co-factor incorporation. The UM team has also identified a small set of new enzyme homologs from extreme microorganisms with superior stabilities that have yielded initial structural data for further engineering. In addition, a new bioelectrocatalytic reactor system has been developed that can be 3D printed and used for enzyme attachment to electrodes. In summary the project has generated critical basic information for the further development of this class of enzymes for the electricity driven reduction of CO2 into formate as platform chemical for upgrading into various other chemicals, including fuels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Narrow-Channel, Fluidized Beds for Effective Particle Thermal Energy Transport and Storage

Colorado School of Mines (Mines) led this program in collaboration with Sandia National Laboratories (Sandia) to characterize narrow-channel fluidized beds of aluminosilicate particles – supplied by Carbo Ceramics – as a means for releasing high-temperature thermal energy in particle heat exchangers and for capturing concentrated solar energy in indirect particle receivers. Single-channel, heat transfer experiments at Mines and reduced-order 1-D models and 3-D two-fluid, CFD models explored many aspects of counterflow, bubbling fluidized beds (net downward particle flow and upward gas flow) for enhancing particle-wall heat transfer at elevated temperatures. Results at Mines on single-channel test sections consistently showed that mild bubbling fluidization increases particle-wall heat transfer coefficients (h T,w ) regularly by more than 4.0x over h T,w values without fluidization at similar conditions (mean particle diameter d p , bed depth Δz b , and bed particle temperatures T p ). Insights from lab-scale tests and modeling studies provided Nusselt number correlations for h T,w and informed the design and fabrication (by Vacuum Process Engineering) of a nominal 40-kWth, particle-sCO 2 plate heat exchanger (HX) with 12 parallel narrow-channel, fluidized beds bounded by stainless-steel walls with embedded microchannels for high-pressure sCO 2 flows. Tests of the 40-kW th HX at the particle-sCO 2 HX test stand at Sandia's National Solar Thermal Test Facility (NSTTF) were limited, due to HX design, to particle inlet temperatures T p,in ≤ 520°C with maximum sCO 2 outlet temperatures T sCO2,out ≈ 440°C, which are well below design conditions for a primary HX in a sCO 2 power cycle for a Gen-3 concentrating solar power (CSP) plant. Total heat transfer $\dot{Q}_{HX}$ remains relatively constant with increased fluidization for fixed particle and sCO 2 inlet conditions because higher h T,w due to fluidization is offset by increased axial dispersion, which suppresses temperature differences between the particles and sCO 2 in the counterflow configuration. The axial dispersion reduces the effective overall heat transfer coefficient U based on T p,in to values around 200 W m -2 K -1 .

14 SOLAR ENERGY↗

Cavity-modified molecular dipole switching dynamics

Polaritonic states, which are formed by resonances between a molecular excitation and the photonic mode of a cavity, have a number of useful properties that offer new routes to control molecular photochemistry using electric fields. To provide a theoretical description of how polaritonic states affect the real-time electron dynamics in molecules, a new method is described where the effects of strong light–molecule coupling are implemented using real-time electronic structure theory. The coupling between the molecular electronic states and the cavity is described by the Pauli–Fierz Hamiltonian, and transitions between polaritonic states are induced via an external time-dependent electric field using time-dependent configuration interaction (TDCI) theory, producing quantum electrodynamics TDCI (QED-TDCI). This method is used to study laser-induced ultrafast charge transfer and dipole-switching dynamics of the LiCN molecule inside a cavity. The increase in cavity coupling strength is found to have a significant impact on the energies and transition dipole moments of the molecule–cavity system. The convergence of the polaritonic state energies as a function of the number of included electronic and photonic basis states is discussed.

Chemistry↗

Data-driven acceleration of thermal radiation transfer calculations with the dynamic mode decomposition and a sequential singular value decomposition

In this work, we present a method for accelerating discrete ordinates radiative transfer calculations for radiative transfer. Our method works with nonlinear positivity fixes, in contrast to most acceleration schemes. The method is based on the dynamic mode decomposition (DMD) and using a sequence of rank-one updates to compute the singular value decomposition needed for DMD. Using a sequential method allows us to automatically determine the number of solution vectors to include in the DMD acceleration. We present results for slab geometry discrete ordinates calculations with the standard temperature linearization. Compared with positive source iteration, our results demonstrate that our acceleration method reduces the number of transport sweeps required to solve the problem by a factor of about 3 on a standard diffusive Marshak wave problem, a factor of several thousand on a cooling problem where the effective scattering ratio approaches unity, and a factor of 20 improvement in a realistic, multimaterial radiating shock problem.

97 MATHEMATICS AND COMPUTING↗