Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mix Testing”

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 55 records · Page 3

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗

Improving GCM Predictability of Mixed-Phase Clouds and Aerosol Interactions at High Latitudes with ARM Observations

The goal of this project is to improve the predictability of mixed-phase clouds and aerosol interactions in the Community Atmosphere Model version 6 (CAM6) through comparison with the ARM observations. There are three main objectives of the proposed study: (1) Improve the representation of ice microphysical processes in mixed-phase clouds; (2) Test the performance of ice microphysics in CESM-CAM6 with the ARM observations in northern and southern high latitudes; and (3) Examine mixed-phase cloud microphysics-aerosol-turbulence-radiation interactions in CESM-CAM6. In this project, we have (1) Improved the representation of ice microphysical processes in mixed-phase clouds in CESM-CAM6 by implementing the marine organic aerosol (MOA) and treating the ice nucleating particles (INPs) from MOA and its impacts on mixed-phase clouds. We made a first attempt to represent different secondary ice production (SIP) mechanisms in a GCM (CESM2-CAM6). We found that misrepresentation of these ice formation processes in global climate models (GCMs) leads to too weak negative cloud feedback over the Southern Ocean (SO) and too high climate sensitivity in the models. In addition to ice formation processes, we improved the treatment of ice depositional growth through the Wegener–Bergeron–Findeisen (WBF) process by considering the subgrid heterogeneous distributions between liquid droplets and ice crystals in mixed-phase clouds; (2) Tested the performance of our improved representations of ice microphysics with the ARM observations at high latitudes.

54 ENVIRONMENTAL SCIENCES↗

Comprehensive Assessment of the Accuracy of the Ideal Adsorbed Solution Theory for Predicting Binary Adsorption of Gas Mixtures in Porous Materials

Quantifying the adsorption of chemical mixtures in porous adsorbents is critical to developing these materials for useful separation applications. The ideal adsorbed solution theory (IAST) is the most widely applied mixing theory for predicting mixture adsorption using single-component adsorption data, but a perceived lack of experimental data has limited previous efforts to explore the accuracy of IAST in a systematic way. In this paper, we take advantage of a large collection of binary experimental data for gas adsorption that became available recently to tackle this issue. Specifically, we identify more than 400 examples in which binary adsorption data and single-component data are available in the same publication and apply IAST to all these examples. This analysis includes experimental data from 63 gas mixtures of 37 different molecular species and 174 different adsorbents. In addition to being the most systematic evaluation to date of the accuracy of IAST for gas adsorption, these data will be valuable for future efforts to test or develop mixing theories that improve upon IAST.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of oxidation rates and volatile oxidation products for HTGR graphite matrix material exposed to steam atmospheres

High-temperature gas-cooled reactors (HTGRs) in operation use tristructural isotropic (TRISO) particles embedded in graphite and carbonized resin matrix to form the fuel element. This graphite matrix material serves as a supportive structural element, heat transfer medium, and neutron moderator. In HTGR designs, fuel compacts are exposed to helium coolant, which facilitates high outlet temperatures (750°C 2 O, 800°C 2 , and H 2 , are quantified in varied oxidant atmospheres using a coupled thermogravimetric analyzer and mass spectrometer. Furthermore, oxidation rates reported here for varied steam (H 2 O [g]) atmospheres are predominantly linear and comparable with literature values in the range of tested temperatures (800–1200°C). Changes in dominant matrix oxidation products from primarily CO to a mixture of CO, CO 2 , and H 2 were observed at higher temperatures (≥1000°C) and steam atmospheres (≥5 kPa pH 2 O). Kinetic data indicates that there was no shift in oxidation regime with chemical oxidation occurring at all temperatures and H2O (g) atmospheres tested. These data provide insight into the oxidation behavior of graphite matrix material and will inform future testing conditions, notably mixed atmospheric conditions, of HTGR fuel elements.

36 MATERIALS SCIENCE↗

Recovery of High Purity Rare Earth Elements (REE) from Coal Ash via a Novel Electrowinning Process (Final Report)

The goal of this project was to produce a high purity, separated (>90%) rare earth oxide (REO) product from coal-based sources. This product was achieved using the following three steps: 1. Battelle’s ADP: The ADP process involves pretreatment of ash (milling and caustic leaching), leaching of pretreated ash with nitric acid, roasting of loaded acid solution containing REE, and water leaching of residual from roasting step, resulting in a REE pregnant solution. 2. Solvent Extraction (SX) Upgrading: The SX involves the extraction of REE from a pregnant REE solution using an organic extractant (note after extraction, the aqueous phase was a residual solution and the organic phase contained REE), and stripping of REE from the organic phase product obtained from the extraction using an acid solution (the stripped solution was the final purified product containing REE, organic traces, and traces of other metals for further separation in the electrowinning process). 3. Rare Earth Salts’ (RES) Electrochemical Separation and Purification Process: The staged electrochemical process was separated REO products from the mixed REE solution. Testing by RES was also included in the investigation of options to minimize the solvent extraction steps. This report covers laboratory testing, production of a high purity coal-based REO, process design of the overall REE recovery and purification process, technoeconomic assessment (TEA), and a commercialization plan discussing the overall REE recovery and purification process.

01 COAL, LIGNITE, AND PEAT↗

Evaluation of High Level Waste Sludge Processing Behavior

The U.S. Department of Energy’s (DOE) Hanford Site has 177 underground storage tanks that contain wastes from past nuclear fuel reprocessing and waste-management operations. Over 20% of this waste is in the form of an insoluble sludge that will require slurry modification before its transfer to the Waste Treatment and Immobilization Plant (WTP). Specific WTP acceptance criteria for waste feed delivery describe the physical and chemical characteristics of the waste that must be met before the waste is transferred to the WTP. One challenging requirement relates to the undissolved solids (UDS) composition in a waste feed because the waste contains solid particles that settle, and their concentration and relative proportion can change during the transfer of the waste in individual batches. A key uncertainty is the ability to transfer and mix wastes with large variations in UDS concentrations and resulting settling rates. To address this uncertainty, a number of small scale mixing and settling tests have been conducted to determine the mobilization performance of variable chemistry simulants. Comparison of the size and density of the particulate for each simulant to that of southeast area Hanford sludge was made using metrics for particle mobilization, suspension, settling, and pipeline transfer where dependance on particle size and density may be different, including: 1. Settling velocity, 2. Critical shear stress for erosion, 3. Just-suspended impeller speed, and 4. Pipeline critical transport velocity. Existing high-level waste sludge data has shown the effect that increasing Al concentration has on resulting settled solids. This differential settling of particles in the sludge has the possibility of resulting in solids segregation during feed preparation and uneven particle distribution during pipeline transportation or mixer jet pump operations. Understanding the predictive capabilities of HLW solids settling and transport as well as potential remedies for addressing disparate sludge behaviors can help provide technical guidance during HLW flowsheet planning.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Cellulose Nanofibrils Dewatered with Poly(Lactic Acid) for Improved Bio-Polymer Nanocomposite Processing

Cellulose nanofibrils (CNFs) have theoretically ideal properties for bio-based composite applications; however, the incorporation of these materials into polymers is made challenging by the strong binding of water to CNFs and the irreversible agglomeration of CNFs during drying. Previous methods used “contact dewatering”, wherein the addition of wood flour (WF) to CNFs facilitated the mechanical removal of water from the system via cold pressing, which showed potential in producing dried CNF fibrils attached to wood fibers for biocomposite applications. In this work, the viability of contact dewatering with poly(lactic) acid (PLA) powder for PLA/CNF composites was evaluated. The energy efficiency of dewatering, preservation of nanoscale CNF morphology, and mechanical properties were examined by mixing wet CNFs with pulverized PLA at various loading levels, pressing water out of the system, and compression molding and shear mixing composites for testing. The most impressive results from this dewatering method were the preservation of micron-to-nanoscale fibers with high aspect ratios in PLA-CNF composites; increased strength and modulus of 1.7% and 4.2%, respectively, compared to neat PLA; equivalent or better properties than spray-dried nanocellulose at similar loading levels; and an 11-194x reduction in drying energy compared to spray-drying CNFs.

Chemistry↗

Leveraging PHIL for Inverter Functionality Requirement Evaluation to Ensure a Reliable Grid

This presentation showcases NREL's ongoing research on advanced Multi-point Power Hardware-in-the-Loop (PHIL) systems, enabling comprehensive evaluation of interoperability, stability, and wide-area stability in complex power grids. Key features include high-power PHIL capabilities, seamless PHIL Interfaces for effortless Grid-Following (GFL) and Grid-Forming (GFM) mode switching, and advanced multi-domain PHIL/Controller Hardware-in-the-Loop (CHIL) capabilities for evaluating diverse technology mixes, facilitating rigorous testing and validation of emerging power systems for reliable integration, enhanced resilience, and optimal performance.

lab capabilities↗

Comparative study with two reduced test facilities for air-cooled RCCS scaling law

Korea Atomic Energy Research Institute (KAERI) and Argonne National Laboratory (ANL) conducted steady-state tests at constant heat flux conditions to investigate the scale effect on heat removal behavior of the air-cooled Reactor Cavity Cooling System (RCCS). Two differently scaled-down test facilities were used in this study, KAERI’s ¼ scale with 4.0-m cavity height, and ANL’s ½ scale with 6.8-m cavity height. Two scaling laws were proposed to simulate the radiation across the cavity and the buoyancy-driven duct flow in the riser, respectively. The test matrix for KAERI and ANL test facilities focused on the scaling effect of the heated riser length with constant heat flux at the PMR200 RCCS design. The test results showed that mixed convection in the riser duct is an important factor for accurately extrapolating the thermo-fluid behavior in the prototype from the test results in the scale-down facilities. The system analysis code, GAMMA+, with improved heat transfer models predicted fairly well the air-cooled RCCS test data from two facilities. GAMMA+ analysis results showed that the predicted radiation fractions on the heated plate in the scale-down test conditions were larger than those on the reactor vessel in the prototype. The scaling law for air-cooled RCCS was improved by considering the mixed convection in the riser duct and the same radiation fraction on the heated plate. The mixed convection effect was calculated by the height ratio. The same radiation fraction provided a conservative extrapolation from the scale-down test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

Nuclide inventory characterization of EBR-II MOX fuel test pins

Characterization of 3 mixed oxide (MOX) fuel pins from the Experimental Breeder Reactor II (EBR-II) SPA-2/-2B experiment program was performed in support of the Advanced Reactor Experiments for Sodium Fast Reactor Fuels (ARES) joint project between Idaho National Laboratory (INL) and the Japan Atomic Energy Agency (JAEA). The characterization efforts were performed to enable transient testing in the Transient Reactor Test (TREAT) facility at INL and will include post irradiation examination (PIE) and other analyses. MOX fuel pins UW02009, 02011, and 06024 were historically irradiated from 1989 to 1994 for a total of 822.42 effective full power days (EFPDs) and peak burnup of approximately 14.3 at.% (atomic percent of heavy metal atoms, ∼134.4 GWd/t). They had since been placed in storage until a future use was identified. Characterization of the fission products and activation of each fuel pin was necessary to enable removal from storage and subsequent irradiation testing or PIE activities. The initial fuel pin geometry and composition were established using historical documentation. Irradiation history details were utilized to compute irradiated and decayed fuel pin masses, atomic densities, and activities. A pellet-wise axial neutron flux was also computed to assist in evaluation of the axial distribution of {sup 235}U and {sup 239}Pu, per pellet, within each of the test fuel pins. The results computed herein supported experiment design for advanced high-burnup MOX fuel designs under slow transient overpower (MOXTOP) conditions, which was Phase I of the ARES project. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The effect of current and lambda on white-etch-crack failures

White etching cracks (WECs) have been associated with premature failure of wind turbine roller bearings. Various drivers for the generation of WECs have been identified such as loading conditions, slip, steel quality, lubrication, hydrogen embrittlement, corrosion fatigue cracking, and stray electrical currents passing through the surface. Here, in this work, a benchtop test rig utilizing a three-ring-on-roller test configuration was used to investigate the effect of electrical current and operation in different lubricating regimes, defined by lambda (λ), on high-quality bearing steel samples tested in a commercially available power transmission EP gear lubricant. It was observed that there is an inverse correlation between the magnitude of electric current applied to the ring/ roller system and time-to-failure. Higher current magnitudes lead to shorter time-to-failure than lower current magnitudes, with macropitting as the main failure mode. Sub-surface investigation revealed the presence of WECs in all cases. For the same current magnitude, tests conducted in boundary and mixed lubrication regimes showed that time-to-failure increased as lambda increased, and the tests resulted in WEC related macropits, whereas tests conducted in near-hydrodynamic regime resulted in surface damage with no macropit. It was also noted that a shift toward near-hydrodynamic lubrication resulted in a distinct surface distress on the roller surface. Furthermore, there seems to be a transition in the mixed regime during which the surface distress occurred. The damage on the surface of the test samples resembled non-spatially, periodic, groove-like corrugations and, in some cases, crater-like depressions. Sub-surface imaging, performed by sequential sectioning, revealed the presence of WECs in all cases, and broad, branching cracks that were more prevalent under the more severe boundary conditions.

17 WIND ENERGY↗

Probing the role of solids loading and mix procedure on the properties of acoustically mixed materials for additive manufacturing

We report resonant acoustic mixing has been of particular interest for use in additive manufacturing since viscous, solids-loaded materials can be difficult to mix and inhomogeneity has adverse effects on print quality. In this study, we detail a method to iterate through different formulations and mix procedures and assess mixture quality. The approach utilizes a constant pressure-driven flow test to collect statistics on flow rate through a standard geometry. This test is first applied to well mixed formulations containing particulate solids and find that the results are sensitive to viscosity produced by solids content. We then consider the formulations at various stages of mixing and find that the spread in the volume measurements is indicative of the mixing quality; poorly mixed material yields measurements with wide distributions. We believe this testing approach can be useful when screening new formulations, developing mixing processes, or for quality control.

36 MATERIALS SCIENCE↗

Addressing Low-Cost Methane Sensor Calibration Shortcomings with Machine Learning

Quantifying methane emissions is essential for meeting near-term climate goals and is typically carried out using methane concentrations measured downwind of the source. One major source of methane that is important to observe and promptly remediate is fugitive emissions from oil and gas production sites but installing methane sensors at the thousands of sites within a production basin is expensive. In recent years, relatively inexpensive metal oxide sensors have been used to measure methane concentrations at production sites. Current methods used to calibrate metal oxide sensors have been shown to have significant shortcomings, resulting in limited confidence in methane concentrations generated by these sensors. To address this, we investigate using machine learning (ML) to generate a model that converts metal oxide sensor output to methane mixing ratios. To generate test data, two metal oxide sensors, TGS2600 and TGS2611, were collocated with a trace methane analyzer downwind of controlled methane releases. Over the duration of the measurements, the trace gas analyzer’s average methane mixing ratio was 2.40 ppm with a maximum of 147.6 ppm. The average calculated methane mixing ratios for the TGS2600 and TGS2611 using the ML algorithm were 2.42 ppm and 2.40 ppm, with maximum values of 117.5 ppm and 106.3 ppm, respectively. A comparison of histograms generated using the analyzer and metal oxide sensors mixing ratios shows overlap coefficients of 0.95 and 0.94 for the TGS2600 and TGS2611, respectively. Overall, our results showed there was a good agreement between the ML-derived metal oxide sensors’ mixing ratios and those generated using the more accurate trace gas analyzer. This suggests that the response of lower-cost sensors calibrated using ML could be used to generate mixing ratios with precision and accuracy comparable to higher priced trace methane analyzers. This would improve confidence in low-cost sensors’ response, reduce the cost of sensor deployment, and allow for timely and accurate tracking of methane emissions.

03 NATURAL GAS↗

Measuring the weak mixing angle at SBND

The weak mixing angle provides a sensitive test of the Standard Model. We study SBND's sensitivity to the weak mixing angle using neutrino-electron scattering events. We perform a detailed simulation, paying particular attention to background rejection and estimating the detector response. We find that SBND can provide a reasonable constraint on the weak mixing angle, achieving 8% precision for $10^{21}$ protons on target, assuming an overall flux normalization uncertainty of 10%. This result is superior to those of current neutrino experiments and is relatively competitive with other low-energy measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energy Exascale Earth System Model v2.0.1

First patch release of v2.0.0 Changes since v2.0.0 [Important change] Fix ocean threading bug seen in debug cases on Chrysalis with Intel 20.0.4. Was introduced around time of v2.0.0 tag. Does not change v2.0.0 answers on Chrysalis because those didn't use threading or debugging. [EAM] Add semi-lagrangian tracer transport for theta-l (F90 and C++), add new algorithm for finding tropopause, add DSCREAM to allow v2 and SCREAM settings in same code such as adjust_ps [EAM-MMF] 60L default, allow C++ back end of RRTMGP (EAM too). [EAMxx] add nu-top functionality, fix forcing functor, add ttype9 and dcmip2012 tests 2.1, 2.2, and 3 HOMME: remove obsolete remap algs, option to specify dynamics alg indep of tracer, new sponge layer, add imex tests [ELM] Add topography-based subgrid (topounits), add FATES-ELM Nitro., Phos. and CH4 coupling, add land-use ts for NARRM, add lulc for SSP3 RCP7, Fix nutrient fertilization exp test and carbon isotope flux, Fix xactive lnd dry deposition, add lake water storage option, fix plant hydraulics 2d params, fix carbon budget calc, fix soil nutrient conc. bug, fix mosart dam bug, add test for new ELM, MOSART features, fix bug in O3 dry dep stomatal resistances, fix plant hydraulics restart BFB error, update mkmapdata. [MOSART] fix bug for reading the latitude from an unstructured input file, fix oversat in bubble test. [MPAS-ocean] Add CFC11, CFC12 tracers, add 2D spherical transport tests, fix del4 tracer mixing, add MARBL ocean tracer mixing, modify harmonic analysis options, add GPU port of vmix routines, fix calc of ML-averaged BV freq. [MPAS-seaice] Change extents of initial polar disks for oRRS18to6v3 grid, fix ice BGC with MARBL, update spherical test cases, fix DON coupling, Remove Cf from sea ice constants. [MPAS-landice] add CRYO1850-4xCO2 compset [CIME] add GCP, ANL GCE, Spock, Perlmutter, deprecate config_compilers.xml, fix and clean-up cmake macros, fix slurm bindings, refactor CIME internal testing, cleanup SCORPIO perf data, allow position independent compset naming, [also] update v2 benchmarking suite, extend e3sm_prod with throughput and memory checks

E3SM Project, DOE↗

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗