Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Doubling method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Strong-field nonsequential double photoionization using virtual-detector theory with path summation

We present an ab initio study of the nonsequential strong-field ionization dynamics of a model two-electron atom with helium character. Single- and double-ionization events are characterized and displayed using detector signals extracted at different points in the two-electron two-dimensional space. The double photoelectron momentum distribution is calculated via coherent path-summation over virtual-particle trajectories. A comparison is made between the momentum distributions obtained with the virtual-detector method, the Schrödinger wave function, and the time-dependent surface flux method developed by Tao and Scrinzi [New J. Phys. 14, 013021 (2012)]. Insights into different ionization and electron recollision pathways are gained from detailed virtual-particle tracking and energy-time readouts. Furthermore, this study demonstrates the extension of virtual-detector theory to strong-field multielectron quantum dynamics and highlights the importance of the evolving quantum phase in quasiclassical electron propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

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

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2021

Validation of process simulation workflow for thermosetting prepreg platelet molding compounds

Continuous carbon fiber prepreg slit and cut into rectangular platelets has proven to be a useful material for net shape molding of semi-structural and structural components in the aerospace and automotive industries. Furthermore, to assist the designer in use of these prepreg platelet molding compounds, sometimes called carbon fiber sheet molding compounds, simulation tools are required that can predict the as-manufactured fiber orientation state which has a significant impact on the resulting performance. Herein, an analysis workflow for design-enabling predictions is demonstrated for a double dome geometry with two different initial charge configurations. Here, the workflow is validated through comparison with experimental short shots, orientation state, and stiffness trends. Significantly, to complete the validation, a method is proposed for determining the confidence bounds on measured orientation state enhancing the results of optical microscopy which can only produce a small sample of platelet orientations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Strong field ionization and dissociation dynamics of vinyl bromide (C 2 H 3 Br) initiated by few-cycle pulses

Photoelectron-photoion (PEPICO) and photoion-photoion (PIPICO) coincidence measurements, coupled with a few-cycle intense laser field, were employed to investigate strong field ionization/dissociation dynamics of vinyl bromide (C₂H₃Br). The angular streaking technique was used to map the recoil-frame angle-dependent ionization rates to identify ionizing orbitals. Dissociative single ionization was mainly attributed to ionizing the lower lying molecular orbitals, suggesting that direct ionization dominates with few-cycle pulses. Three dissociative double ionization channels, both prompt and delayed fragmentation, were identified. Theoretical analysis using time-dependent configuration interaction with complex absorbing potential (TDCI-CAP) and high-level electronic structure methods was performed to elucidate the underlying ionization and dissociation mechanisms.

Olowolafe, Temitayo A. [Wayne State University, De↗
DOE OSTI2023

Second-Phase Hepatitis C Plasma Viral Kinetics Directly Reflects Reduced Intrahepatic Burden of Hepatitis C Virus

Mathematical models explain how antivirals control viral infections. Hepatitis C virus (HCV) treatment results in at least 2 phases of decline in viremia. Here, the first phase reflects clearance of rapidly produced virions. The second phase is hypothesized to derive from loss of infected cells but has been challenging to prove. Using single-cell methods, we quantified the number of hepatitis C virus (HCV)-infected hepatocytes in liver biopsies taken before and within 7 days of initiating direct-acting antivirals (DAAs) in a double-blinded randomized controlled trial testing 2 (sofosbuvir-velpatasvir) versus 3 (sofosbuvir-velpatasvir-voxilaprevir) DAAs. We employed thousands of intrahepatic measurements in 10 persons with chronic genotype 1a HCV infection: median proportion of infected hepatocytes declined from 11.3% (range, 1.3%–59%) to 0.6% (range, <0.3%–5.8%), a loss of 75%–95% infected hepatocytes. Plasma viremia correlated with numbers of HCV-infected hepatocytes (r = 0.77; P < .0001). Second-phase plasma dynamics and changes in infected hepatocytes were indistinct (P = .16), demonstrating that second-phase viral dynamics derive from loss of infected cells. DAAs led to a decline in intracellular HCV RNA and interferon-stimulated gene expression (P < .05 for both). We proved that second-phase viral dynamics reflect decay of intrahepatic burden of HCV, partly due to clearance of HCV RNA from hepatocytes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Strong-field ionization of water: Nuclear dynamics revealed by varying the pulse duration

Polyatomic molecules in strong laser fields can undergo substantial nuclear motion within tens of femtoseconds. Ion imaging methods based on dissociation or Coulomb explosion therefore have difficulty faithfully recording the geometry dependence of the field ionization that initiates the dissociation process. Here we compare the strong-field double ionization and subsequent dissociation of water (both H 2 O and D 2 O) in 10-fs and 40-fs 800-nm laser pulses. We find that 10-fs pulses turn off before substantial internuclear motion occurs, whereas rapid internuclear motion can take place during the double ionization process for 40-fs pulses. The short-pulse measurements are consistent with a simple tunnel ionization picture whose predictions help interpret the motion observed in the long-pulse measurements.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Riverbed Temperature and 4D ERT Monitoring Reveals Heterogenous Horizontal and Vertical Groundwater-Surface Water Exchange Flows Under Dynamic Stage Conditions

Groundwater surface water exchange plays a critical role in physical, biological, and geochemical function of coastal and riverine systems. Observing exchange flow behavior in heterogeneous systems is a primary challenge, particularly when flows are governed by dynamic river stage or tidal variations. In this paper we demonstrate a novel application of time-lapse 3D electrical resistivity tomography and temperature monitoring where an array of thermistors installed beneath a riverbed double as resistivity electrodes. We use the array to monitor stage driven exchange flows over a 6-day period in a dynamic, stage-driven high order stream. We present a method for addressing the otherwise confounding effects of the moving river-surface boundary on the raw resistivity data, thereby enabling successful tomographic imaging. Temperature time-series at each thermistor location and time-lapse 3D images of changes in bulk electrical conductivity together provide a detailed description of exchange dynamics over a 10-meter by 45-meter section of the riverbed, to a depth of approximately 5 m. Results reveal highly variable flux behavior throughout the monitoring domain including both horizontal and vertical exchange flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development of the Mobile Systems for Conditioning of Disused Sealed Radioactive Sources in Serbia - 20105

Sealed radioactive sources (SRS) are being used worldwide in the field of medicine, agriculture, industry and research. They can be found in mobile as well as stationary devices. SRS contains radioactive material that is (a) permanently sealed in a capsule or (b) closely bounded and in a solid form. The capsule or material of an SRS should be strong enough to maintain leak tightness under the conditions of use and purpose for which the source was designed, also in case of accidents. In this case only emitted radiation is utilized. Firstly, the hazard from external radiation has to be considered, but the possibility of contamination due to fracture of the capsule should not be disregarded. The radioactive sources are composed of the radiating isotope contained in the filling medium, the single or double isotope holder that partially or totally surrounds the filling medium, the outer cover that contains the parts mentioned above and the capsule closed airtightly by welding or using some other method. The capsule must be tested for leakage periodically. If the SRS is no longer needed (e.g. replaced by a different technique) or it becomes useless for the intended application (e.g. the activity becomes too weak, the equipment containing the source works poorly or becomes obsolete, the source is damaged or leaking) it is considered disused. Disused sealed radioactive sources (DSRS) are typically conditioned and disposed if a facility is available. If the disposal option is not available, conditioned DSRS should be stored under proper conditions. In some cases, the radionuclide(s) in DSRS can be recovered/recycled or the DSRS can be repurposed for other applications. Conditioning of DSRS ensures containment of the radioactive material, provides confinement for leaking sources, provides sufficient radiation shielding, reduces storage/disposal volume by allowing consolidation of multiple sources into a single storage/disposal container, facilitates transport operations and contributes to safety and security as well. Typically, conditioning technologies are deployed either as permanently installed stationary systems in centralized or mobile on-site waste processing facilities, or in a mobile configuration. Centralized stationary facilities provide a single processing location for multiple users that requires transport of the waste to the facility. On the other hand, mobile systems may be provided for the selection and application of the optimum technology for a specific waste stream (such as DSRS) by bringing the process to the point where the waste is generated. In addition, mobile systems could offer additional flexibility by sharing equipment among multiple waste generating sites for processing campaigns that vary in duration, from very short periods to several years. The term 'mobile processing system' refers to any radioactive waste processing system or component which is designed to be transportable and which is not considered permanently installed. Two mobile system for conditioning of disused sealed radioactive sources are developed in the Public Company Nuclear Facilities of Serbia. Development of these mobile systems was supported by SRB9005 national project via Technical Cooperation of the IAEA. The first mobile system, built inside the 20 feet ISO container, will be used for conditioning of DSRS category 3 to 5. The second mobile system, built inside the 7 m long vehicle (Iveco Daily Van), will be used for dismantling of ionizing smoke detectors mostly with Am-241 sources. Designs of the mobile systems were defined in cooperation with two companies from Belgium (Belgoprocess and Leniko) and a Croatian company Ekoteh as well as with the support of the IAEA experts. The generic safety assessment and operational procedures for the mobile systems are developed. Based on safety assessment the acceptance criteria and operational limits and conditions are established. Operational procedures include: (a) equipment and material requirements, (b) assembling procedure of the mobile unit, (c) procedure for acceptance of devices for dismantling and conditioning, (d) dismantling procedure for devices to recover the DSRS, (e) characterization of DSRS, (f) encapsulation procedure of DSRS, (g) disassembling procedure of the mobile unit, and (h) keeping records, identification and traceability. In addition, radiation safety, health safety, security and emergency preparedness plans are prepared. The generic safety assessment and operational procedures could be updated with site specific requirements, DSRS inventory, and different needs for future customers. Developed mobile systems could be used in all situations when it is feasible to perform conditioning of DSRS on the spot in the county and worldwide. Development of these mobile units was just the first step to create the Reference Center for Radioactive Waste Treatment and Disused Radioactive Sources Conditioning for Small Facilities which can become a regional training center in the future, and/or as a tool for comprehensive national search and secure programmes. In the next phases development of e-learning platforms and blended learning packages as well as application for the IAEA Qualified Technical Centre (QTC) for the management of DSRS is foreseen. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Solid-state NMR of unlabeled plant cell walls: high-resolution structural analysis without isotopic enrichment

Abstract Background Multidimensional solid-state nuclear magnetic resonance (ssNMR) spectroscopy has emerged as an indispensable technique for resolving polymer structure and intermolecular packing in primary and secondary plant cell walls. Isotope ( 13 C) enrichment provides feasible sensitivity for measuring 2D/3D correlation spectra, but this time-consuming procedure and its associated expenses have restricted the application of ssNMR in lignocellulose analysis. Results Here, we present a method that relies on the sensitivity-enhancing technique Dynamic Nuclear Polarization (DNP) to eliminate the need for 13 C-labeling. With a 26-fold sensitivity enhancement, a series of 2D 13 C– 13 C correlation spectra were successfully collected using the unlabeled stems of wild-type Oryza sativa (rice). The atomic resolution allows us to observe a large number of intramolecular cross peaks for fully revealing the polymorphic structure of cellulose and xylan. NMR relaxation and dipolar order parameters further suggest a sophisticated change of molecular motions in a ctl1 ctl2 double mutant: both cellulose and xylan have become more dynamic on the nanosecond and microsecond timescale, but the motional amplitudes are uniformly small for both polysaccharides. Conclusions By skipping isotopic labeling, the DNP strategy demonstrated here is universally extendable to all lignocellulose materials. This time-efficient method has landed the technical foundation for understanding polysaccharide structure and cell wall assembly in a large variety of plant tissues and species.

09 BIOMASS FUELS↗
DOE OSTI2024

Robust Tensor Hypercontraction of the Particle–Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle–particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O(N 6 ) to O(N 5 ) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. Here, we find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5x.

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

Distributed Frequency Divider for Power System Bus Frequency Online Estimation Considering Virtual Inertia From DFIGs

In this work, a distributed frequency divider is proposed to estimate power system bus frequency with a limited number of PMUs while considering the inertial contributions from double-fed induction generators (DFIGs). The key idea is to reformulate the original frequency divider by modeling the contributions of DFIGs with inertia emulation and external regional power system through Thevenin equivalents. The distributed frequency divider is general and able to estimate local bus frequencies in a distributed manner. Besides, only the knowledge of the rotor speeds of synchronous generator, boundary bus frequencies and terminal bus frequencies of DFIGs as well as the admittance matrix are required to estimate the frequencies at all buses. This is drastically different from existing approaches that require system observability by PMUs to monitor all bus frequencies. Numerical results carried out on the IEEE 39-bus and modified 118-bus power systems with DFIGs demonstrate and effectiveness and robustness of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Kinetic simulations of ignited mode cesium vapor thermionic converters

Cesium vapor thermionic converters are an attractive method of converting high-temperature heat directly to electricity, but theoretical descriptions of the systems have been difficult due to the multi-step ionization of Cs through inelastic electron–neutral collisions. This work presents particle-in-cell simulations of these converters, using a direct simulation Monte Carlo collision model to track 52 excited states of Cs. Here, these simulations show the dominant role of multi-step ionization, which also varies significantly based on both the applied voltage bias and pressure. The electron energy distribution functions are shown to be highly non-Maxwellian in the cases analyzed here. A comparison with previous approaches is presented, and large differences are found in ionization rates due especially to the fact that previous approaches have assumed Maxwellian electron distributions. Finally, an open question regarding the nature of the plasma sheaths in the obstructed regime is discussed. The one-dimensional simulations did not produce stable obstructed regime operation and thereby do not support the double-sheath hypothesis.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

A double peak in the seasonality of California's photosynthesis as observed from space

Solar-induced chlorophyll fluorescence (SIF) has been shown to be a powerful proxy for photosynthesis and gross primary productivity (GPP). The recently launched TROPOspheric Monitoring Instrument (TROPOMI) features the required spectral resolution and signal-to-noise ratio to retrieve SIF from space. Here, we present a downscaling method to obtain 500 m spatial resolution SIF over California. We report daily values based on a 14 d window. TROPOMI SIF data show a strong correspondence with daily GPP estimates at AmeriFlux sites across multiple ecosystems in California. We find a linear relationship between SIF and GPP that is largely invariant across ecosystems with an intercept that is not significantly different from zero. Measurements of SIF from TROPOMI agree with MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation indices – the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near-infrared reflectance of vegetation index (NIR$_v$) – at annual timescales but indicate different temporal dynamics at monthly and daily timescales. TROPOMI SIF data show a double peak in the seasonality of photosynthesis, a feature that is not present in the MODIS vegetation indices. The different seasonality in the vegetation indices may be due to a clear-sky bias in the vegetation indices, whereas previous work has shown SIF to have a low sensitivity to clouds and to detect the downregulation of photosynthesis even when plants appear green. We further decompose the spatiotemporal patterns in the SIF data based on land cover. The double peak in the seasonality of California's photosynthesis is due to two processes that are out of phase: grasses, chaparral, and oak savanna ecosystems show an April maximum, while evergreen forests peak in June. An empirical orthogonal function (EOF) analysis corroborates the phase offset and spatial patterns driving the double peak. The EOF analysis further indicates that two spatiotemporal patterns explain 84 % of the variability in the SIF data. Results shown here are promising for obtaining global GPP at sub-kilometer spatial scales and identifying the processes driving carbon uptake.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗